Cosmology in the laboratory: an analogy between hyperbolic metamaterials and the Milne universe
David Figueiredo Departamento de F´ısica, CCEN, Universidade Federal da Para´ıba, Caixa Postal 5008, 58051-900, Jo˜aoPessoa, Para´ıba, Brazil
Fernando Moraes∗ Departamento de F´ısica, Universidade Federal Rural de Pernambuco, 52171-900, Recife, Pernambuco, Brazil
S´ebastienFumeron and Bertrand Berche Statistical Physics Group, Laboratoire de Physique et Chimie Th´eoriques, Universit´ede Lorraine, B.P. 70239, 54506 Vandœuvre les Nancy, France (Dated: November 21, 2017) This article shows that the compactified Milne universe geometry, a toy model for the big crunch/big bang transition, can be realized in hyperbolic metamaterials, a new class of nanoengi- neered systems which have recently found its way as an experimental playground for cosmological ideas. On one side, Klein-Gordon particles, as well as tachyons, are used as probes of the Milne geometry. On the other side, the propagation of light in two versions of a liquid crystal-based meta- material provides the analogy. It is shown that ray and wave optics in the metamaterial mimic, respectively, the classical trajectories and wave function propagation, of the Milne probes, leading to the exciting perspective of realizing experimental tests of particle tunneling through the cosmic singularity, for instance.
I. INTRODUCTION in the entire four-dimensional spacetime [2]. The ques- tion we want to address in this work is: can we simulate Initial conditions are always a trouble in cosmology but the passage through the singularity in the laboratory? can be circumvented by cyclic universe models, like an Recent advances in the field of metamaterials suggest endless repetition of big crunches followed by big bangs, this possibility. Arising from a pioneering idea by Vese- for instance. This is in fact an old idea that can be traced lago [5] and developed by Pendry [6], metamaterials are back to ancient mythologies. Even without referring to artificial media structured at subwavelength scales, such initial conditions some issues remain in these models, no- that their permittivity and permeability values can be tably the passage through the singularity, the transition taylored quasiarbitrarily (for instance, they may exhibit from big crunch to big bang. Recently, safe transition negative refractive index). Because analogy is a power- has been proposed [1, 2], where the singularity is noth- ful tool that a physicist possesses to arrive at an under- ing more than the temporary collapse of a fifth dimen- standing of the properties of nature, along with the fact sion. The three space dimensions remain large and time that visualizations of celestial object features in the lab- keeps flowing smoothly. A toy model for the geometry oratory have been a charming subject for human beings of this transition is the compactified 2D Milne universe over centuries, analogue gravity became an active field in arXiv:1706.05470v2 [gr-qc] 19 Nov 2017 [3] which is essentially a double cone in 3D Minkowski physics with the help of metamaterials. Therefore, sev- spacetime. Because there is still a singularity in one spa- eral works based on different results of general relativity tial dimension, a physically correct model should be able were done, including topics like analog spacetimes [7, 8], to describe the propagation of a particle through it. An time travel [9], cosmic strings [10], celestial mechanics example of such a model [4] revealed that the compacti- [11], black holes [12, 13] and wormholes [14], just to cite fied Milne universe seems to model the cosmic singularity a few. However, most of the approaches aforementioned in a satisfactory way. It is important to stress that the rely on transformation optics [15], where the permittiv- collapse here is far less severe than in ordinary 4D general ity ij and permeability µij tensors must be numerically relativity, because it happens just in one spatial dimen- equal (impedance matched). It turns out that sometimes sion, whereas in the ordinary case the collapse happens it is difficult to obtain such a constraint experimentally. Nevertheless, a particularly promising class of such media is that of hyperbolic metamaterials, for which one of the ∗ Also at Departamento de F´ısica, Universidade Federal da eigenvalues of either the permittivity or the permeability Para´ıba, Caixa Postal 5008, 58051-900, Jo˜aoPessoa, Para´ıba, tensors do not share the same sign with the two others Brazil (thus, they do not need to be impedance matched). Hy- 2 perbolic metamaterials are now extensively studied, both for practical purpose (enhanced spontaneous emission, hyperlenses [16–18]) but also for modeling cosmological phenomena (metric signature transition [19], modeling of time [20, 21], and even inflation [22]). Therefore, since inflationary models of the universe can be mimicked through hyperbolic metamaterials, it can be interesting to see if the cyclic/ekpyrotic ones also can. As far as light propagation is concerned, a (2+1) Minkowski spacetime can be simulated with such materials. As will be seen below, the quantum dynamics of Klein-Gordon particles through the big crunch/big bang transition may be experimentally verified with light propagating in a suitable hyperbolic metamaterial. FIG. 1. Plane sector which can be bent around a cone. The deficit angle δ = 2π(1 − sin θ) reduces the range of the polar angle ϑ to [0, 2π sin θ). II. THE MC UNIVERSE
In this section we will summarize the definition and one. Making these changes in Eq. (6), we obtain the 2D features of the Milne universe and the compactified Milne induced line element for the double cone universe, MC . However, first it is useful to analyze the geometry of a rectangular cone in R3 to get some insight. ds2 = dl2 + l2 sin2 θdφ2. (8) Therefore, let θ be the angle between the generatrix and the axis. In Cartesian coordinates (x, y, z) the cone sur- Concerning the cosmological model, consider the usual face will have the following form Robertson-Walker line element with negative spatial cur- vature [24] x2 + y2 = z2 tan2 θ, (1) ds2 = −dt2 + a(t)2 dχ2 + sinh2 χ dΩ2 , (9) with parametric equations where a(t) is the expansion factor of the universe and x = r sin θ cos (ϑ/ sin θ) , (2) dΩ2 = dθ2 + sin2 θdφ2 the solid angle. For a linear ex- y = r sin θ sin (ϑ/ sin θ) , (3) pansion factor a(t) = t, we obtain the Milne universe metric z = r cos θ, (4) ds2 = −dt2 + t2 dχ2 + sinh2 χ dΩ2 , (10) and position vector which was proposed by Edward Arthur Milne in 1933 r = [r sin θ cos (ϑ/ sin θ) , r sin θ sin (ϑ/ sin θ) , r cos θ] , and represents a homogeneous, isotropic and expanding (5) model for the universe [25] which grows faster that simple due to the fact that it is locally isometric to a piece of the cold matter dominated or radiation dominated universes. plane [23]. In the parametrization above 0 ≤ ϑ < 2π sin θ We are interested in the hypersurfaces where dΩ = 0. is the polar angle and 2π(1 − sin θ) is the deficit angle Then, Eq. (10) reduces to (see Fig. 1). Rescaling the polar angle as φ = ϑ/ sin θ one recovers the parametrization in terms of spherical ds2 = −dt2 + t2dχ2. (11) coordinates (r, θ, φ). Thus, the line element at the cone surface is, in terms of φ, Next, we make a coordinate transformation to new variables (T,X) given by ds2 = dr2 + r2 sin2 θdφ2. (6) T = t cosh χ, (12) Then, because θ = constant we have an induced 2D met- X = t sinh χ, (13) ric on the cone with the following metric tensor which leads to the following form for the line element in 1 0 g = , (7) Eq. (11), namely ij 0 r2 sin2 θ ds2 = −dT 2 + dX2, (14) i 1 2 where gij = ei · ej (ei = ∂r/∂x , x = r, x = φ). To extend the above case for the double cone, we mod- which is the usual Minkowskian metric in two dimen- ify the radial coordinate r making r → l, so that l ∈ R. sions. However, one must stress that the Milne universe Thus, l > 0 for the upper cone and l < 0 for the lower only covers half of the Minkowski spacetime. To see why, 3 consider the lines of constant values χ = χ0 in Eqs. (12) Taking into account all previous parametrizations, the and (13), parametric equations (17)–(19) become
X = T tanh χ0. (15) x = t sinh θ cos φ, (21) y = t sinh θ sin φ, (22) As a result, we have straight lines in a Minkowskian di- z = t cosh θ, (23) agram. Taking the limits χ0 → ±∞ in the equation above, one gets the light rays X = ±T , which form the √ where sinh θ = κ, cosh θ = 1 + κ2 and φ ∈ S1. As for boundaries of the past and future light cones from the the metric in M , origin (0, 0). Thus, one is confined in this region where C X2 −T 2 < 0, since the slope of the line given by Eq. (15) ds2 = −dt2 + κ2t2dφ2. (24) goes from 0 (for χ0 = 0) to ±1 (for χ0 → ±∞). Another important point to stress is that the three- The presence of κ2 in the above metric indicates a conic dimensional space for the comoving Milne observers has singularity of the curvature at the origin. As we will see infinite extension. The reason is due to the fact that the below, this has important implications to the geodesics lines of constant t = t0 are hyperbolas, namely and wave functions of particles approaching the singu- larity, which acts as a filter for classical particles and a 2 2 2 T − X = t0, (16) phase eraser for quantum ones. Due to the fact that the ordinary cone surface is embedded in three-dimensional and each one of the hyperbolas has infinite length. This Euclidean space and the MC universe has a cone surface fact is expected since the Milne universe has a negative embedded in three-dimensional Minkowski space, they spatial curvature [25]. share some similarities. For instance, the Milne metric Therefore, in order to compactify the Milne universe, tensor we follow the usual approach [1–3, 26] and let the vari- −1 0 able χ acquire some period Π. The meaning of this is gij = 2 2 . (25) as follows, in the Minkowski diagram (T,X) the lines 0 κ t X = 0 and X = T tanh Π should be identified as one, for is similar to the two-dimensional induced metric tensor instance. Therefore, because one is constrained between given by Eq. (7). Also, as a last comparison, let us con- these two lines, the Milne universe now has a finite length sider the Laplace-Beltrami operator, which will be used and is called compactified Milne universe, MC (see Fig. later in this paper. Thus, 2). Following Refs. [4, 27], one can visualize the MC i 1 p ij ∆ = ∇i∇ = ∂i |g|g ∂j . (26) universe through a mapping into a three-dimensional p|g| Minkowski space with ds2 = −dz2 + dx2 + dy2, with ij where g is the inverse metric tensor and g = det(gij). x = tκ cos (χ/κ) , (17) From Eqs. (8) and (24), we will have y = tκ sin (χ/κ) , (18) 1 ∂ ∂ 1 ∂2 p 2 ∆Cone = l + , (27) z = t 1 + κ , (19) l ∂l ∂l l2sin2θ ∂φ2 1 ∂ ∂ 1 ∂2 where t ∈ 1 and 0 < κ ∈ 1 is a constant parameter R R ∆MC = − t + 2 2 2 , (28) for compactifications (in Refs. [27, 28] it is shown that t ∂t ∂t κ t ∂φ −1 κ is related to the rapidity tanh v = 2πκ of a finite which also look similar due to the resemblances between Lorentz boost). Thus, without loss of generality, we take the two metrics, aside from the timelike behavior in t and the period of χ as 2πκ. Moreover, as in the previous case the hyperbolic cone angle factor κ2 = sinh2 θ. of the ordinary cone, we rescale χ by φ = χ/κ (and hence giving a period of 2π for φ). Solving Eqs. (17)–(19) one gets the following constraint equation III. CLASSICAL PARTICLE IN MC κ2 x2 + y2 = z2. (20) In this section we will obtain the path followed by free 1 + κ2 classical particles (timelike geodesics) in MC . Thus, we start with the relativistic action (in natural units c = 1) This equation is similar to Eq. (1) because the space is given by [29] Euclidean for the planes z = constant and also due to the periodicity of φ. As x, y, z can assume both positive and Z q i j negative values, Eq. (20) represents a double cone with S = −m −gijx˙ x˙ dλ, (29) a vertex at (0, 0, 0) in the 3D Minkowski space (see Fig. 3). However, due to the timelike aspect of z, the cone where the dot notation stands for a derivative with re- angle θ is a hyperbolic angle with tanh2 θ = κ2/1 + κ2. spect to an affine parameter λ along the curve (which can 4
which for the Milne metric (24) gives the following equa- tions X = ± T X = T tanh 2 φ¨ + t˙φ˙ = 0, (34) t t¨+ κ2t φ˙2 = 0. (35) From the equations above, one can see that the nonva- T φ φ nishing Christoffel symbols are Γtφ = Γφt = 1/t and t 2 Γφφ = κ t. Furthermore, since φ is a cyclic coordinate in the Lagrangian, ∂Lkin/∂φ = 0 and the angular momen- ˙ 2 2 ˙ tum pφ = ∂Lkin/∂φ = κ t φ is conserved. This fact is expressed in Eq. (34). 2 X Substituting pφ = κ J = constant in Eq. (35) and after some manipulations, the second geodesic equation becomes FIG. 2. Minkowski diagram showing the M universe (gray C κ2J 2 region) delimited by the identification between the T -axis and t˙2 − = 2E, (36) the blue line X = T tanh Π. The orange lines are showing t2 the light rays X = ±T while the hyperbolas are surfaces for where E is a constant of integration. The left-hand side of constant t in Eq.(16). Eq. (36) is nothing else but minus the square modulus of the velocity ui =x ˙ i. As a result, recalling that in special 8 i relativity one has u ui = −1 and λ is an affine parameter, i 4 it is useful to choose u ui = −2E = −1. Then, Eq. (36) becomes z 0 κ2J 2 t˙2 − = 1. (37) -4 t2 -8 Solving Eq. (37) for t˙, we get -8 8 r -4 4 t2 + κ2J 2 0 ˙ 0 t = ± 2 , (38) y 4 -4 t 8 -8 x which leads to a simple integration of the form FIG. 3. Double cone surface corresponding to Eqs. (1) and Z tdt Z (20) for 3D Euclidean space and 3D Minkowski space, respec- √ = ± dλ, (39) 2 2 2 tively. The upper half part corresponds to l, t > 0 while the t + κ J lower one is for values of l, t < 0. and therefore to the parametric equation t = t(λ) given by i i p i j be the proper time) and L(x , x˙ , λ) = −m −gijx˙ x˙ is p 2 2 2 the Lagrangian. The variation of the action, δS, will be t(λ) = ± ∆λ − κ J , (40)
Z q where ∆λ = λ − λ0, with λ0 being a constant of integra- i j δS = −m δ −gijx˙ x˙ dλ, (30) tion. In order to obtain the parametric equation φ = φ(λ), thus, the least action principle δS = 0 demands that we substitute Eq. (40) in the relation t2φ˙ = J and solve Z q for φ˙. Namely, i j δ −gijx˙ x˙ dλ = 0, (31) J φ˙ = , (41) from which the geodesic equations follow ∆λ2 − κ2J 2 i i j k x¨ + Γjkx˙ x˙ = 0, (32) leading to the following integration i Z Z dλ where Γjk are the Christoffel symbols of the second kind. dφ = J . (42) An alternative, and more direct route to get the ∆λ2 − κ2J 2 geodesic equations and Christoffel coefficients is through As a result, the parametric equation φ = φ(λ) is the Euler-Lagrange equations for the kinetic Lagrangian 1 1 −1 ∆λ L = g x˙ ix˙ j, (33) φ(λ) = φ0 − coth , (43) kin 2 ij κ κJ 5 where φ0 is a constant. Combining Eqs. (40) and (43), the same reasoning applies for particles coming from the one gets the equation of the trajectory orange curve. For the purpose of completeness, we perform the same t t = ± 0 , (44) calculations for spacelike geodesics, which can be world- |sinh κ∆φ| lines of tachyons [32]. The procedure is the same, the major change being the choice of the constant of inte- where t0 = κ |J| and ∆φ = φ−φ0. The above equation is a Poinsot spiral [30], with the + (−) sign corresponding gration in Eq. (36). Therefore, because the momen- to the upper (lower) cone. Furthermore, for J 6= 0 one tum and velocity must be spacelike, it is useful to choose i can see that depending on the choice of the signal (or u ui = −2E = 1. Thus, Eq. (36) becomes sheet of the cone), the particle always remains in the up- κ2J 2 per or lower region with no link between those regions of t˙2 − = −1, (46) spacetime (see Fig. 4), but for J = 0, which corresponds t2 to φ = φ0 and t = ± ∆λ, the geodesics are straight lines. the rest of the procedure being rather straightforward On the other hand, from Eqs. (12) and (13) one can and leading to the following parametric equations see that all the trajectories given by Eq. (44) are straight lines in Minkowski spacetime, p t(λ) = ± −∆λ2 + κ2J 2, (47) t0 1 ∆λ X = ± + T tanh κφ0, (45) φ(λ) = φ + tanh−1 , (48) cosh κφ0 0 κ κJ where for J = 0 one recovers Eq. (15). As a result, the where ∆λ = λ − λ , with λ, λ and φ denoting an affine particle indeed can travel from one cone to the other, but 0 0 0 parameter and constants of integration, respectively. As such trajectories are very unstable since small perturba- for the trajectory t = t(φ), tions in the value of J = 0 cause large deviations on the trajectories. A similar result was found for a classical t0 particle in a double cone in Ref. [31], where a classical t = ± , (49) cosh κ∆φ nonrelativistic particle constrained to a double cone only crosses the vertex through straight lines. where t0 = κ |J| and ∆φ = φ−φ0. Equation (49) is also a However, we are dealing with a toy model for a cyclic Poinsot spiral, with the + (−) sign corresponding to the universe with contraction and expansion phases joined upper (lower) cone. Furthermore, as in the previous case by a cosmic singularity. Thus, as pointed out in Refs. the trajectories are straight lines in Minkowski spacetime, [4, 27], due to the fact that timelike geodesics coincide with the trajectories of test particles, which do not dis- t0 X = ∓ + T coth κφ0, (50) tort the spacetime around them, there is no obstacle for sinh κφ0 such geodesics to reach (leave) the singularity. Further- more, if one postulates that such a particle arriving at the where the − (+) corresponds to the upper (lower) sheet. singularity coming from the lower cone is “annihilated” From Eq. (49) the variable t has t0 as its limiting value. at the singularity, while another one is “created” at the To clarify this result, we remark that it was pointed out upper cone, and considering that the Cauchy problem is by Feinberg [32] that tachyons lose energy√ as they speed i not well defined at t = 0, there are several types of prop- up. Thus, from Eq. (49) the velocity dx dxi/dt (with agation depending on the way a particle travels towards i = 1, 2, 3), the singularity (see Sec. III of Ref. [27]). Although all of them must be consistent both with the constraint given dφ t0 κt = (51) p 2 2 by Eq. (37) and with the fact that the angular momen- dt t0 − t tum J is constant. It is out of the scope of this work to discuss the details and properties of such propagations. ranges from (1, ∞) as t ranges between 0 and t0, respec- Next we present the timelike geodesics in a polar plot, tively (note that for a lightlike interval, κt dφ/dt = 1). which provides a clearer way to visualize the trajectories Therefore, as the tachyon comes accelerating from the (see Fig. 5). As already shown in Fig. 4, the trajecto- singularity, we can see from Eq. (46) that the time com- ries in the lower cone are going towards the cosmic sin- ponent of the momentum p0 = t˙ (which is associated gularity, which means that particles in the red (orange) with the energy) decreases, reaching its minimum value curve have positive (negative) angular momentum J and at t0. From that point on it starts to increase as the therefore are spinning in the counterclockwise (clockwise) tachyon decelerates towards the singularity (see Fig. 6). direction. For the upper cone, particles traveling along This leads to an interpretation of tachyons being created the green (blue) curve have positive (negative) values of and annihilated at the same sheet of the cone, which can J and are spinning counterclockwise (clockwise). As a be the upper sheet or the lower one. Every annihilation result, in a transition from the red to the green curve on a sheet creates a tachyon on the other sheet, in an the angular momentum is conserved, while for a transi- endless cycle. This can be seen in a polar plot as closed tion to the blue one, J would not be conserved. Clearly, curves in spacetime (see Fig. 7). 6
3 t + t + t 2 t 1 y 0 0 arrow of time t / 1