1. Introduction
Currently, most of the known information about the Universe outside the Earth is obtained observing electromagnetic waves. This is especially true for the deep cosmos, even though gravitational astronomy [
1] and neutrino telescopes [
2,
3] are on the road of opening the doors of the multi-messenger observations [
4]. So far, however, the experiments that have shaped our vision about the Universe and its evolution, such as Planck [
5], Wilkinson Microwave Anisotropy Probe [
6], WiggleZ Dark Energy Survey [
7], Background Imaging of Cosmic Extragalactic Polarization [
8], All Sky Automated Survey for Super Novae [
9], Sloan Digital Sky Survey [
10], Hubble Space Telescope [
11,
12], James Webb Space Telescope [
13] etc., all detect electromagnetic signals emitted somewhere faraway and reaching us after long journey through the spacetime. From them we know that the Universe is currently in an accelerated expansion epoch [
14]. They also show that the flat Friedman – Lemaître – Robertson – Walker (FLRW) metric [
15] is so far the best candidate for large scale metric (and is considered as a standard model for cosmological observations).
In this work, we consider the problem how electromagnetic waves propagate in FLRW metric, which is quite important to observational cosmology. Same problem is addressed in [
16,
17,
18,
19,
20] and many others. The reason for this is that electrodynamics in curved spacetime is quite peculiar. The strong equivalence principle states that the laws of physics are the same in all inertial frames. However, there is a result [
21] that on a world line one can put any metric to Minkowski form (the inertial coordinate system) and to nullify its first derivatives but not the second and higher ones. As a result any equation in which metric’s second and higher derivatives appear is not the same in Minkowski spacetime and in local inertial coordinate system in more general spacetime. The Electrodynamics in curved spacetimes is an example of such theory because of presence of Ricci tensor in its equation of motion.
Here we consider plane and spherical waves described by their electromagnetic potential in generalized Lorentz gauge. We show that in both cases under consideration each transverse component of the potential decouples from other components. The corresponding equations are quite simple and can be solved exactly. Their solutions demonstrate relativistic redshift, amplitude decrease (fading) and dispersion. We also speculate about the connection between electromagnetic wave equation in FLRW metric and Klein–Gordon equation with mass determined by the Ricci tensor. It turns out that for a special case of equation of state the electromagnetic wave describes a tachyon.
2. Theory
Here our approach is to treat the electromagnetic waves as test particles in predefined metric (flat FLRW one in our case). In other words we consider the electromagnetic field as weak, i.e. we do not take into account the gravitation produced by the wave itself. In addition we suppose the existence of a medium in the spacetime which, first, produce the FLRW metric, and second, is inert to the electromagnetic interaction.
The electromagnetic potential
defines the electromagnetic field tensor
which, due to its anti-symmetry and the symmetry of Christoffel symbols, has the same form in any frame
As a result the electromagnetic potential has the usual gauge symmetry
where
is an arbitrary function. In order to fix this arbitrariness we use the generalized Lorentz gauge
The free equation of motion is
which in the gauge (
2) takes the form
Here
is the d’Alembert operator defined for the covariant derivatives
, and
is the Ricci tensor
The following relations between the Riemann tensor
and Ricci tensor are used in the derivation of eq.(
4) :
We want to stress that eq.(
4) is not the wave equation (
) and that it can not be cast into it with the help of change of coordinates.
3. Results
3.1. Plane electromagnetic wave in FLRW metric
The description of a plane wave can be done most easily in the Cartesian form of the FLRW metric. In this case the invariant length is
where
is the scale "parameter" and we work in units where the speed of light
.
For simplicity, we consider a wave propagating in
z direction, i.e. we suppose that
. For this functional choice eq.(
4) simplifies and each of the potential transverse components
and
(i.e.,
and
components of the potential) decouple from all other components. Note that here the transverse potential determines transverse electric and magnetic fields. Therefore, we deal with a transverse electromagnetic wave which remain such forever. We want to stress that because we are free to choose the spacial coordinate system orientation, so our consideration describe in fact the general solution of a plane wave.
The equation, which both
and
satisfy is:
Here we denote with dots the derivatives over time and with primes — derivatives with respect to z.
We look for a solution of eq.(
7) with separated variables
. In this case eq.(
7) leads to the following two equations:
where
is some constant.
At this point we want to make a small comment about eq.(
8). The self-consistency of the assumption for separation of variables requires that
does not participate in it. In other words the equation in question will have the same form for any
, even for
. So, it will be the same as in flat Minkowski spacetime.
In view of the above comment, it is not a surprise that eq.(
8) is well known. Its general real solution for
(which we suppose) is:
where
are arbitrary real constants.
Now we move to eq.(9). In order to find its solution we make two consecutive ansatzes. First, let
which leads to the following differential equation for
:
The second ansatz is
where the argument of the function
is the so called conformal time and the differential equation for
is
Note that eq.(
14) has the same form as eq.(
8) and, accordingly, has the same type solution. Therefore, the general solution of eq.(9) is
3.2. Spherical electromagnetic wave in FLRW metric
We consider the spherical coordinate system as the appropriate one for description of spherical waves (provided it is positioned and oriented accordingly to the source and observer). In spherical coordinates the invariant length is
and it determines the particular form of eq.(
4) in this case.
We consider the electromagnetic wave propagation in
r direction, supposing that
, and in addition supposing separation of variables. Therefore, we are looking for transverse potential
(the time dependence of both transverse components is indeed the same — see below).
Here we see a flaw in our considerations because we postulate the existence of transverse constant vector field on a sphere. However, such field does not exist. Any transverse vector field on a sphere has to have at least one zero (or singularity). Therefore, one of our assumptions, or both of them, are incorrect. Nevertheless, we continue our analysis supposing it is valid at most only approximately in a small vicinity on the sphere for which we suppose that resides at . This should be sufficient for a cosmological observer who, anyway, can gather information only locally.
It turns out that the evolution of each of the transverse components decouples from all other components as in the plane wave case. The corresponding equations of motion are:
Now the primes denote derivatives with respect to
r coordinate and
. The fact that eqs.(
18,) are different is a surprise, taking in to account that orientation of the coordinate system around the axis of observation is a matter of our choice.
Let us consider eq.(19) first. As it has been already mentioned we look for a solution with separated variables, so that
. It turns out that the differential equation for
thus defined coincides with eq.(9) considered in the plane wave case. For the equation satisfied by the function
we get:
where, as in the case of plane wave,
is some positive parameter. Its general solution is
where
and
are the spherical Bessel functions of first and second kind.
Next, we consider eq.(
18) and we again look for a solution in the form
. Once again, the differential equation for
is exactly eq.(9). For the function
we obtain the following equation:
The general solution of eq.(22) is
where
and
are the Bessel functions of first and second kind. Note however, that
dependence reappears in eq.(23) through the quantity
b. This contradicts to our initial ansatz. So, we reconsider our assumption and now we are looking for a solution in the form
Skipping the details it is possible to show that the equation for at is exactly eq.(19). This resolves the problem with rotational symmetry around observation axis. The result also rules out the possible polarization of the light induced by its propagation through the space which is suggested by eqs.(21,23).
3.3. Cosmological redshift, fading and dispersion
In this section we consider some quantities that are of interest for the electromagnetic waves observers.
First, we define the redshift
z. It can be derived in number of different ways in the cosmological context. Usually, one uses the geodesic of the photon to obtain the so called cosmological redshift [
22]. One can also solve the scalar wave equation for light with certain initial data[
23,
24]. One can also use the Einstein-Maxwell’s equations in the space-time defined with its metric, to obtain the solutions and obtain it from there. This is the approach we will use in this paper, so the redshift is defined as :
where
is the wave (angular) frequency at the moment
and
is the observed frequency at the moment
. Note that time dependence for spherical and plane waves is one and the same and is given by eq.(
15). Therefore, in both cases the angular frequency at moment
t is determined by the wave period
, (
) and for
we have the following equation:
As a result, the frequency at the moment
t is
and so, the redshift is
where subscripts stand for plane and spherical wave.
Next, we define the fading
as the amplitude decrease due to the propagation of the wave:
Here is the amplitude of the wave at the moment and is the observed amplitude at the moment , .
Another interesting characteristics of the waves is their dispersion
. Like the redshift there are lot of different ways to define the dispersion. Probably the simplest is an analogue of the group velocity dispersion:
However, in all cases under consideration here the frequency is linear with respect to wave number, but the spatial part of the wave phase is not in the case of spherical wave (see below). So, our definition for dispersion is
where
is the argument of the sine or cosine part of the wave function. The quantity is determined only at the moment
.
Applying the above definitions for both plane and spherical waves, we obtain the following results for the plane wave
The results for spherical wave need some more work. Having in mind that
we can put eq.(21) which defines the spacial depending component of transverse potential in the following form:
We set the origin of coordinate system at the geometric center of spherical wave. We consider the propagation along
axis of a fixed phase of outgoing wave, so that
which can be used to determine
from
,
and
. Therefore, the fading and dispersion are:
3.4. Photon mass
The Ricci tensor in FLRW metric is diagonal. Its space components
(no summation over
i),
(i.e.
,
and
components in Cartesian coordinates and
,
and
components in spherical coordinates), have one and the same value
. Except for the very first moments of the Universe evolution this value is a very slowly changing function of time and therefore, the term
in eq.(
4) behaves like mass term and eq.(
4) itself — like Klein–Gordon equation. The mass term is zero only for
or for
. These are the only two cases for which eq.(
4) coincides with the wave equation of Minkowski Electrodynamics. Note that power law behavior of the scale factor (
) is an exact solution of the Friedman equations for Universe full with pure fluid with constant equation of state
where
p is the fluid pressure and
is its density. In this more general case the photon mass is
which is positive for
. However, in the case
the mass squared is negative and the photon represents a
tachyon.
It is interesting to see what is the behavior of the photon mass for exponentially growing scale factor
. In the case we get
i.e. it is positive constant (exactly).
4. Discussion
The obtained closed expressions for the electromagnetic plane and spherical wave solutions allow us to estimate the observed redshift, fading and dispersion. The obtained redshift (eq.(27)) coincides with the standard one known in the literature [
15,
25]. The amplitude decrease given in eq.(30) is refereed in [
26] as "adiabatic". The extension to a non-flat Universe considered therein have shown that
decrease is possible. On the other hand the amplitude decrease of spherical wave (eq.(34)) demonstrates the dependence of the fading on the wave number. The fading is increasing when the wavelength is increasing. A result, connected to the above one, is given by eq.(35) and predicts nonzero dispersion of spherical wave.
Finally we want to make some comments about possible connection between eq.(
4) and the Klein–Gordon one. We consider the fact that photons can be massive as very interesting and the fact that they can be tachyons as even more interesting. Note that tachyonic models have a long history in the cosmology. Some of them originate as special cases of k-essence theories with Dirac-Born-Infeld (DBI) action [
27]. On the other hand k-essence theories [
28,
29] are used to describe early inflation and dark energy trough a minimally coupled scalar field with non-canonical kinetic term. In the tachyonic models [
30,
31], universe expansion (possibly accelerated) is produced while the tachyon rolls down towards its minimum. Tachyons have also been discussed in terms of the so called tachyonic preheating [
35] which may lead to explosive particle production. Ghost tachyons (i.e. models with negative sign of
in the Lagrangian) have been shown to cross the phantom line of the equation of state
[
36]. How our result can be positioned in these researches? Well, note that our result is valid not only for electromagnetic field but for any
gauge field. For instance it can be the fundamental
field of electroweak interaction which exists before spontaneous symmetry breaking. Suppose that in some early stage of the Universe the perfect fluid which determines the metric in it is with
equation of state (which can be achieved in some
k-essence models [
37,
38,
39]). In this case the fundamental
gauge field becomes tachyonic and in the spirit of the articles cited above these tachyons can drive inflation. It will continue as long as
and in the same time the tachyons will roll to a energy minimum where they are infinitely fast.
Author Contributions
Conceptualization, M.S.; methodology, M.S..; software, D.S and M.S.; validation, D.S and M.S.; formal analysis, D.S and M.S.; investigation, D.S and M.S.; writing—original draft preparation, M.S.; writing—review and editing, D.S and M.S. All authors have read and agreed to the published version of the manuscript.
Funding
This research was funded by Bulgarian National Science Fund grant KP-06-N58/5.
Conflicts of Interest
The authors declare no conflict of interest.
Abbreviations
The following abbreviations are used in this manuscript:
FLRW |
Friedman – Lemaître – Robertson – Walker |
References
- Abbott, B.; et al. Multi-messenger Observations of a Binary Neutron Star Merger. Astrophys. J. Lett. 2017, 848, L12. [Google Scholar] [CrossRef]
- Aartsen, M.; et al. Observation of High-Energy Astrophysical Neutrinos in Three Years of IceCube Data. Phys. Rev. Lett. 2014, 113, 101101. [Google Scholar] [CrossRef]
- Huang, T.; Li, Z. Neutrino observations of LHAASO sources: Present constraints and future prospects. Mon. Not. Roy. Astron. Soc. 2022, 514, 852–862. [Google Scholar] [CrossRef]
- Addazi, A.; et al. Quantum gravity phenomenology at the dawn of the multi-messenger era —A review. Prog. Part. Nucl. Phys. 2022, 125, 103948. [Google Scholar]
- Aghanim, N.; et al. Planck 2018 results. VI. Cosmological parameters. Astron. Astrophys. 2020, 641, A6. [Google Scholar]
- Hinshaw, G.; et al. Nine-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Parameter Results. Astrophys. J. Suppl. 2013, 208, 19. [Google Scholar] [CrossRef]
- Blake, C.; et al. The WiggleZ Dark Energy Survey: Joint measurements of the expansion and growth history at z < 1. Mon. Not. Roy. Astron. Soc. 2012, 425, 405–414. [Google Scholar]
- Ade, P.; et al. Improved Constraints on Primordial Gravitational Waves using Planck, WMAP, and BICEP/Keck Observations through the 2018 Observing Season. Phys. Rev. Lett. 2021, 127, 151301. [Google Scholar] [CrossRef]
- Neumann, K.; et al. The ASAS-SN Bright Supernova Catalog – V. 2018-2020. Mon. Not. Roy. Astron. Soc. 2023, 520, 4356. [Google Scholar] [CrossRef]
- Zhao, C.; et al. The completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: cosmological implications from multitracer BAO analysis with galaxies and voids. Mon. Not. Roy. Astron. Soc. 2022, 511, 5492–5524. [Google Scholar] [CrossRef]
- Riess, A.; et al. Type Ia supernova discoveries at z > 1 from the Hubble Space Telescope: Evidence for past deceleration and constraints on dark energy evolution. Astrophys. J. 2004, 607, 665–687. [Google Scholar] [CrossRef]
- Freedman, W.; et al. Final results from the Hubble Space Telescope key project to measure the Hubble constant. Astrophys. J. 2001, 553, 47–72. [Google Scholar] [CrossRef]
- Yuan, W.; Riess, A.; Casertano, S.; Macri, L. A First Look at Cepheids in a Type Ia Supernova Host with JWST. Astrophys. J. Lett. 2022, 940, L17. [Google Scholar] [CrossRef]
- Riess, A.; et al. Observational evidence from supernovae for an accelerating universe and a cosmological constant. Astron. J. 1998, 116, 1009–1038. [Google Scholar] [CrossRef]
- Peacock, J.A. Cosmological Physics; Cambridge University Press, 1999. [Google Scholar]
- Mashhoon, B. Electromagnetic Waves in an Expanding Universe. Phys. Rev. D 1973, 8, 4297. [Google Scholar] [CrossRef]
- Cohen, J.; Kegeles, L. Electromagnetic fields in curved spaces - a constructive procedure. Phys. Rev. D 1974, 10, 1070–1084. [Google Scholar] [CrossRef]
- Cabral, F.; Lobo, F. Electrodynamics and Spacetime Geometry: Foundations. Found. Phys. 2017, 47, 208–228. [Google Scholar] [CrossRef]
- Cabral, F.; Lobo, F. Electrodynamics and spacetime geometry: Astrophysical applications. Eur. Phys. J. Plus 2017, 132, 304. [Google Scholar] [CrossRef]
- Asenjo, F.; Hojman, S. Birefringent light propagation on anisotropic cosmological backgrounds. Phys. Rev. D 2017, 96, 044021. [Google Scholar] [CrossRef]
- Lichnerowics, A. Elements of Tensor Calculus; John Wiley and Sons, 1962. [Google Scholar]
- Carroll, S.M. The Cosmological constant. Living Rev. Rel. 2001, 4, 1. [Google Scholar] [CrossRef]
- Abbasi, B.; Craig, W. On the initial value problem for the wave equation in friedmann–robertson–walker space–times. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences 2014, 470, 20140361. [Google Scholar] [CrossRef]
- Petersen, O.L. The mode solution of the wave equation in Kasner spacetimes and redshift Math. Phys. Anal. Geom. 2016, 19, 26. [Google Scholar] [CrossRef]
- R. Utiyama, The Theory of Relativity, Atomizdat, 1979.
- Tsagas, C.G. Electromagnetic fields in curved spacetimes. Class. Quant. Grav. 2005, 22, 393–408. [Google Scholar] [CrossRef]
- Calcagni, G.; Liddle, A.R. Tachyon dark energy models: dynamics and constraints. Phys. Rev. D 2006, 74, 043528. [Google Scholar] [CrossRef]
- Armendariz-Picon, C.; Damour, T.; Mukhanov, V. F. k - inflation. Phys. Lett. B 1999, 458, 209–218. [Google Scholar] [CrossRef]
- Bilic, N.; Tupper, G.B.; Viollier, R.D. Cosmological tachyon condensation. Phys. Rev. D 2009, 80, 023515. [Google Scholar] [CrossRef]
- Gibbons, G.W. Cosmological evolution of the rolling tachyon. Phys. Lett. B 2002, 537, 1–4. [Google Scholar] [CrossRef]
- Sen, A. Field theory of tachyon matter. Mod. Phys. Lett. A 2002, 17, 1797–1804. [Google Scholar] [CrossRef]
- Bruneton, J.-P. On causality and superluminal behavior in classical field theories: Applications to k-essence theories and MOND-like theories of gravity. Phys. Rev. D 2007, 75, 085013. [Google Scholar] [CrossRef]
- Kang, J. U.; Vanchurin, V.; Winitzki, S. ; Attractor scenarios and superluminal signals in k-essence cosmology. Phys. Rev. D 2007, 76, 083511. [Google Scholar] [CrossRef]
- Babichev, E.; Mukhanov, V.; Vikman, A. k-essence, superluminal propagation, causality and emergent geometry. JHEP 2008, 2, 101. [Google Scholar] [CrossRef]
- Felder, G. N.; Garcia-Bellido, J.; Greene, P. B.; Kofman, L.; Linde, A. D.; Tkachev, I. Dynamics of symmetry breaking and tachyonic preheating. Phys. Rev. Lett. 2001, 87, 011601. [Google Scholar] [CrossRef] [PubMed]
- Karami, K.; Fahimi, K. Interacting viscous ghost tachyon, K-essence and dilaton scalar field models of dark energy. Class. Quant. Grav. 2013, 30, 065018. [Google Scholar] [CrossRef]
- Garriga, J.; Mukhanov, V. F. Perturbations in k-inflation. Phys. Lett. B 1999, 458, 219–225. [Google Scholar] [CrossRef]
- Armendariz-Picon, C.; Mukhanov, V. F.; Steinhardt, P.J. Essentials of k essence. Phys. Rev. D 2001, 63, 103510. [Google Scholar] [CrossRef]
- Bonvin, C.; Caprini, C.; Durrer, R. A no-go theorem for k-essence dark energy. Phys. Rev. Lett. 2006, 97, 081303. [Google Scholar] [CrossRef]
|
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).