PHYSICAL REVIEW LETTERS week ending PRL 100, 191303 (2008) 16 MAY 2008

Time Drift of Cosmological as a Test of the Copernican Principle

Jean-Philippe Uzan* Institut d’Astrophysique de Paris, Universite´ Pierre and Marie Curie-Paris VI, CNRS-UMR 7095, 98 bis, Bd Arago, 75014 Paris, France

Chris Clarkson† and George F. R. Ellis‡ and Gravity Group, Department of Mathematics and Applied Mathematics, University of Cape Town, Rondebosch 7701, South Africa (Received 10 January 2008; published 16 May 2008) We present the time drift of the cosmological in a general spherically symmetric spacetime. We demonstrate that its observation would allow us to test the Copernican principle and so determine if our is radially inhomogeneous, an important issue in our understanding of . In particular, when combined with distance data, this extra observable allows one to fully reconstruct the geometry of a spacetime describing a spherically symmetric underdense region around us, purely from background observations.

DOI: 10.1103/PhysRevLett.100.191303 PACS numbers: 98.80.Es, 98.62.Py

I. INTRODUCTION of the Lemaıˆtre-Tolman-Bondi (LTB) family (that is, spherically symmetric solution of Einstein equations Cosmological data is usually interpreted under the as- sourced by pressureless matter and no cosmological con- sumption that the Universe is spatially homogeneous and stant). Unfortunately, this simple extension of the RW isotropic. This is justified by the Copernican principle, stating that we are not located at a favored position in depends on two free functions (see below for space. Combined with the observed , this leads to details) so that the reconstruction is underdetermined and a Robertson-Walker (RW) geometry [1], at least on the one must fix one function by hand. Thus, one needs at least scale of the observable Universe. one extra independent observation to reconstruct the ge- This implies that the spacetime metric reduces to a ometry of an LTB universe. A limitation to this reconstruc- single function of the cosmic time, the scale factor at. tion arises because most data lie on our past light cone. This function can be Taylor expanded as at This takes us back to the observational cosmology program a H t t 1 q H2t t 2 , where H is the [10] and the question [11] of how to extract as much 0 0 0 2 0 0 0 0 information as possible about our spacetime from cosmo- Hubble parameter and q0 the deceleration parameter. Low redshift observations [2] combined with the assumption of logical data alone. Among many results, it was demon- almost flatness of the spatial sections, justified mainly by strated [12] that the two free functions of a LTB spacetime the cosmic microwave background data [3,4], lead to the can be reconstructed from the angular distance and number counts, even though evolution effects make it impossible to conclusion that q0 < 0: the expansion is accelerating. This conclusion involves no hypothesis about the theory of be conclusive [13]. gravity or the matter content of the Universe [5], as long Recently, two new ideas were proposed. First, it was as the Copernican principle holds. This has stimulated a realized [14] that the distortion of the Planck spectrum of growing interest in possible explanations [5,6], ranging the cosmic microwave background allows one to test the from new matter fields dominating the dynamics at late Copernican principle. Second, a consistency relation be- times to modifications of general relativity. tween distances on the null cone and Hubble rate measure- While many tests of general relativity on astrophysical ments in RW universes was derived [15], based on the fact scales have been designed [7], the verification of the that the curvature is constant; this also serves as an obser- Copernican principle has attracted little attention, despite vational test of the Copernican principle. the fact that relaxing this assumption may be the most In this Letter, we reconsider the time drift of cosmologi- conservative way, from a theoretical perspective, of ex- cal redshift in spacetimes with less symmetries than the plaining the recent dynamics of the Universe without in- RW universe and we demonstrate how, when combined troducing new physical degrees of freedom [8]. with distance data, it can be used to test the Copernican This possibility that we may be living close to the center principle, mainly because observing the thickening of our (because isotropy around us seems well established past light cone brings new information. As pointed out by observationally) of a large underdense region has attracted Sandage and then McVittie [16], one should expect to considerable interest. In particular, the low redshift (back- observe such a time drift in any expanding spacetime. ground) observations such as the magnitude-redshift rela- This may lead to a better understanding of the physical tion can be matched [9] by a nonhomogeneous spacetime origin of the recent acceleration [17,18], or to tests of the

0031-9007=08=100(19)=191303(4) 191303-1 © 2008 The American Physical Society PHYSICAL REVIEW LETTERS week ending PRL 100, 191303 (2008) 16 MAY 2008 variation of fundamental constants [19]. This measure- implies that the luminosity distance is given by DLy 2 ment, while challenging, may be achieved with extremely 1 z DA. The redshift is given by large telescopes (ELT) and, in particular, it is one of the a uak emission Aw0; 0 main science drivers in design of the Cosmic Dynamics 1 z a ; (2) Experiment (CODEX) spectrograph [20]. Our result may uak observer Aw0;y strengthen the scientific case for this project. where the matter velocity and photon wave vector are a 1 a a 1 a We start by introducing observational coordinates, given by u A w and k AB y, respectively. which allow us to derive the general expression for the We deduce that the isotropic expansion rate, defined by time drift in a spherically symmetric, but not necessarily a 3H rau , is given by spatially homogeneous, universe [see Eq. (6)]. We show that observation of both the luminosity distance and the 1 @wBw; y @wCw; y Hw; y 2 : (3) redshift drift allows one to probe the Copernican principle 3A Bw; y Cw; y at low redshifts, when ‘‘dark energy’’ dominates [see the For the central observer, who sees the Universe isotropic, consistency relation (8)]. We demonstrate that this expres- H is simply the Hubble expansion rate. At small redshifts, sion may be used with distance data to fully reconstruct the Hw; y @wBw;0 Oy, so the Hubble constant is geometry of a LTB spacetime [see Eq. (9)]. Bw;0Aw;0 H0 @wBw0; 0=Bw0; 0Aw0; 0. In the particular case of a dust dominated universe, the II. z_ IN A SPHERICALLY SYMMETRIC UNIVERSE acceleration and vorticity vanish and the fluid four-velocity Observational coordinates.—We consider a spherically can be expressed as the gradient of the proper time along symmetric spacetime in observational coordinates the matter worldlines: ua @at. Since we also have fw; y; #; ’g, where w labels the past light cones of events ua A@aw B@ay we deduce that dt Adw Bdy along the wordline C of the observer, assumed to lie at the so that A @wt and B @yt. The surfaces of simultaneity center so that w is constant on each past light cone, with are thus given by Adw Bdy and we have the integra- a a bility condition @ A @ B 0. u @aw>0, u being the four-velocity of the cosmic fluid, y w a uau 1. y is a comoving radial distance coordinate The covariant derivative of ua is therefore of the form specified down the past light cone of an event O on C. raub Hgab uaubab, where the shear ab is a Many choices are possible, such as the affine parameter symmetric, traceless, and satisfies u ab 0. The scalar 2 ab down the null geodesics from O, the area distance, or the shear ab=2 is consequently the only nonvanish- redshift; whatever choice we assume on the past light cone ing kinematical variable and is given by of O, it is specified on other past light cones through being 2 1 @ B @ C a w w comoving with the cosmic fluid, i.e., y u 0.(#, ’) are w; yp ; (4) ;a 3 A B C angular coordinates based at C and propagated parallelly along the past light cone. The metric in these observational where an arbitrary sign has been chosen. The regularity coordinates is conditions imply w; 00, which is expected since the

2 2 2 expansion is observed to be isotropic about the central ds A w; ydw 2Aw; yBw; ydydw worldline. C2w; yd2; (1) Expression of the redshift drift.—From the expression z (2), it is straightforward to deduce that z_ w w0;y is which is clearly spherically symmetric around the world- given by line C defined by y 0. The requirement that the 2-spheres fw; yg behave regularly around when y ! @wA w0; 0 @wA w0;y const 0 C z_w0;y1 z : (5) implies [10] that Aw; y!Aw; 0 Þ 0, Bw; y! Aw0; 0 Aw0;y 2 Bw; 0 Þ 0, and Cw; yBw; 0y Oy . Now, we can choose w such that Aw0; 01. Then, on There remain two coordinate freedoms in possible re- our past light cone, we can choose y such that scalings of w and y. However, once specified on C, w is @w lnBw0;y@w lnAw0;y. Note however that such a determined on the other worldlines by the condition that choice is not possible for all w. It follows that fw constg are past light cones of events on C. This allows 1 us to arbitrarily choose Aw; 0. Also, once specified on z_w0;y1 zH0 Hw0;yp w0;y: (6) one past light cone, y is determined on all the others 3 because it is a coordinate comoving with the fluid. This This is the general expression for the time drift of the allows us to choose Bw0;y for a given value of w w0. redshift as it would be measured by an observer at the On each past light cone, the cross-sectional area of a center of a spherically symmetric universe. source is related to the solid angle d2 under which it is Robertson-Walker case.—Since 0 for a RW space- 2 2 observed by an observer on C at w w0 by C w0;yd . time, Eq. (6) reduces to the standard Sandage-McVittie This implies that C is the angular distance, DA, i.e., formula. Let us consider the RW metric in conformal DAyCw0;y. The distance duality relation [21] then coordinates 191303-2 PHYSICAL REVIEW LETTERS week ending PRL 100, 191303 (2008) 16 MAY 2008

2 2 2 2 2 2 ds a d d fKd ; (7) This solution depends on three arbitrary functions of r only, Er, Mr, and T0r. Their choice determines where is the radial comoving coordinate and f K the model completely. For instance, E; M; T0 sin; ; sinh according to the curvature of the spatial 2 3 K0r ;M0r ; 0 corresponds to a RW universe. One can sections. Now, with w and y , this leads to further use the freedom in the radial coordinate to fix one of the form (1) for the metric (7) with A B aw y and the three functions at will so that one effectively has only C aw yfKy. H is thus constant on each constant two arbitrary independent functions. Assume we fix Mr. time hypersurface and 0 everywhere. Since We want to determine fEr;T0rg to reproduce some @wAjyconst @a, Eq. (5)givesz_ 1 zH0 Hz observables on our past light cone. This can be represented where we have shifted to cosmic time (dt ad), using parametrically as frz;Ez;T0zg. 1 z a0=a. Let us sketch the reconstruction and use r as the inte- We understand from this exercise why the observational gration coordinate, instead of z. Our past light cone is coordinates are adapted to the computation of z_ in an defined as t t^r and we set RrRt^r;r. The arbitrary spacetime. They are just a generalization to less time derivative of R is given by R_ t^r;rR symmetric spacetimes of the conformal coordinates, where p 1 3 0 ^ the expression is easily obtained in the RW case. 2M0r =Rr2Er. Then we get R tr;r ^ 0 We also conclude that since the three functions (A, B, C) R2rfRr3trT0rR1r=2gE =E 0 are expressed in terms of a unique function, all background R1rT0rRr=r. Finally, more algebra leads to _ 0 3 observations depend on some function of H, z_ being no R t^r;rR3rfR1r3M0r t^rT0r= 2 0 3 0 2 exception. R rgE r=2ErM0r T0r=R R1r=r. Thus, 0 0 Consistency relation.—It follows that, in a RW universe, R_ , R , and R_ evaluated on the light cone are just functions we can determine a consistency relation between several of Rr, Er, T0r, and their first derivatives. Now, the observables. From the metric (7) and the relation for null geodesic equation gives that 1 0 z that follows, one deduces that H z D z 1 ^ r dz 1 z 2 2 1=2 dt R2 H D z , where a prime stands for @ and Dz p ; p R3r; K0 0 z dr 1 2Er dr 1 2Er DLz=1 z. This relation is the basis of the test of the RW structure, as recently proposed in Ref. [15], arguing and there that knowledge of Hz at different redshifts, from, e.g., baryon acoustic oscillations or differential age esti- dR R1r

p mates of passively evolving galaxies, could then be used to 1 R2r: dr 1 2Er check this yields the same value of K0. Here, we argue that it can be implemented using z_z as an observational These are three first-order differential equations relating input. Defining five functions Rr, t^r, zr Er, and T0r. To recon- 2 2 struct the free functions we thus need two observational CopDLz; z_z;z1 K0H0Dz relations. RzDAz is the obvious choice. Then, from 2 0 2 H01 zz_z D z ; (8) z_z, we have the new relation 2 2 we must have CopDLz; z_z;z0 whatever the matter 1 R3 R1 1 R3 R1 content of the Universe and the field equations, since it z_ 1 zH0 2 : 9 R2 R 3 R2 R derives from a purely kinematical relation that does not (9) rely on the dynamics (i.e., the Friedmann equations). This is a consistency relation between independent observables It was shown [9] that the observed DLz can be reproduced that holds in any Robertson-Walker spacetime. from the function T0r assuming E 0, or from the Spherically symmetric spacetimes.—Writing the LTB function Er assuming T0 0. Here, we have shown metric in observational coordinates requires the solution that fEr;T0rg can be completely reconstructed from of the null geodesic equation, which is in general possible the data without assumptions. only numerically. Consider an LTB spacetime with metric 2 2 2 2 2 2 ds dt S r; tdr R r; td ; III. DISCUSSION p where Sr; tR0= 1 2Er and R_ 2 2Mr=Rr; t In this Letter, we have shown that observation of both 2Er, using a dot and prime to refer to derivatives with the luminosity distance and time drift of the redshift as a respect to t and r. The Einstein equations can be solved function of z allows one to construct a test of the Mr 0 Copernican principle. We have derived the general expres- parametrically as fRr; ;t; rg f ;T0r Er sion for z_z in a spherically symmetric spacetime [see Mr g, where is defined by sinh Er3=2 Eqs. (5) and (6)]. This extends the standard computation ; 3=6; sin, and Er2E; 2; 2E according which was restricted to RW spacetimes and was extended to whether E is positive, null, or negative. to almost RW spacetimes only recently [22]. 191303-3 PHYSICAL REVIEW LETTERS week ending PRL 100, 191303 (2008) 16 MAY 2008

As a by-product, this also allowed us to derive the [4] J.-P. Uzan, U. Kirchner, and G. F. R. Ellis, Mon. Not. R. consistency relation (8) between the observed DLz and Astron. Soc. 344, L65 (2003). z_z, thereby extending the result of Ref. [15] to an alter- [5] J.-P. Uzan, Gen. Relativ. Gravit. 39, 307 (2007). nate observable. That is, while Hz characterizes the local [6] B. Ratra and P.J. E. Peebles, Rev. Mod. Phys. 75, 559 isotropic expansion rate, z_z gives access to the expansion (2003). [7] J.-P. Uzan and F. Bernardeau, Phys. Rev. D 64, 083004 between us and a source. In RW models these are trivially (2001); M. White and C. Kochanek, Astrophys. J. 560, 539 related but in general the shear enters these observables (2001); J.-P. Uzan, Int. J. Theor. Phys. 42, 1163 (2003); C. differently, thereby presenting a test of the Copernican Sealfon, L. Verde, and R. Jimenez, Phys. Rev. D 71, principle using only background observations. We have 083004 (2005); F. Bernardeau, arXiv:astro-ph/0509224; shown that we can extract the shear as a function of z J. Martin et al., Phys. Rev. Lett. 96, 061303 (2006); R. and demonstrated that it allows one to close the recon- Caldwell, A. Cooray, and A. Melchiorri, Phys. Rev. D 76, struction problem for a LTB spacetime. 023507 (2007); B. Jain and P. Zhang, arXiv:0709.2375. DLz can be measured from the observation of type Ia [8] G. F. R. Ellis, Nature (London) 452, 158 (2008). supernovae, particularly with actual projects such as [9] M.-N. Ce´le´rier, Astron. Astrophys. 362, 840 (2000); R. K. JDEM, up to redshifts of order 1.5. z_z has a typical Barrett and C. A. Clarkson, Classical Quantum Gravity 17, amplitude of order z 5 1010 on a time scale of 5047 (2000); K. Tomita, Mon. Not. R. Astron. Soc. 326, 287 (2001); H. Iguchi, T. Nakamura, and K.-I. Nakao, t 10 yr, for a source at redshift z 4. This measure- Prog. Theor. Phys. 108, 809 (2002); H. Alnes, M. ment is challenging, and impossible with present-day fa- Amarzguioui, and Ø. Grøn, Phys. Rev. D 75, 023506 cilities. However, it was recently revisited [23] in the (2007); R. Vanderveld, E. Flanagan, and I. Wasserman, context of ELT, arguing they could measure velocity shifts Phys. Rev. D 74, 023506 (2006); J. W. Moffat, J. Cosmol. of order v 1–10 cm=s over a 10 yr period from the Astropart. Phys. 05 (2006) 001; C. H. Chuang, J. A. Gu, observation of the Lyman- forest. It is one of the science and W. Y. Hwang, arXiv:astro-ph/0512651; K. Bolejko, drivers in design of the CODEX spectrograph [20] for the arXiv:astro-ph/0512103; D. Garfinkle, Classical Quantum future European ELT. The study of the precision to which Gravity 23, 4811 (2006); R. Mansouri, arXiv:astro-ph/ we can check the Copernican principle with these two data 0606703; D. J. Chung and A. E. Romano, Phys. Rev. D sets is beyond the scope of this Letter. Indeed, many 74, 103507 (2006); K. Enqvist and T. Mattsson, J. Cosmol. effects, such as of the sources, local gravi- Astropart. Phys. 02 (2007) 019. [10] G. F. R. Ellis et al., Phys. Rep. 124, 315 (1985). tational potential, or acceleration of the may con- [11] J. Kristian and R. Sachs, Astrophys. J. 143, 379 (1966). tribute to the time drift of the redshift. It was shown [22], [12] W. Stoeger, G. F. R. Ellis, and S. Nel, Classical Quantum however, that these contributions can be brought to a Gravity 9, 509 (1992); R. Maartens et al., Classical 0.1% level so that the cosmological redshift is actually Quantum Gravity 13, 253 (1996); M. E. Arau´jo and measured. W. R. Stoeger, Phys. Rev. D 60, 104020 (1999). Future high precision data may thus allow a test of the [13] N. Mustapha, C. Hellaby, and G. F. R. Ellis, Mon. Not. R. Copernican principle, even though observations are local- Astron. Soc. 292, 817 (1999). ized on our past light cone. While important in its own right [14] J. Goodman, Phys. Rev. D 52, 1821 (1995); R. Caldwell for understanding the foundations of our cosmological and A. Stebbins, arXiv:0711.3459. model, it is also critical for our understanding of the [15] C. Clarkson, B. Bassett, and T. Hui-Ching Lu, acceleration of the Universe—it will permit us to be con- arXiv:0712.3457. [16] A. Sandage, Astrophys. J. 136, 319 (1962); G. McVittie, fident that any such acceleration is not simply a misinter- Astrophys. J. 136, 334 (1962). pretation of the data because of incorrectly assuming the [17] K. Lake, arXiv:astro-ph/0703810; L. Amendola, C. geometry of our Universe at low redshift. Quercellini, and A. Balbi, arXiv:0708.1132. We thank B. Bassett, F. Bernardeau, and Y. Mellier for [18] P.-S. Corasaniti, D. Huterer, and A. Melchiorri, Phys. Rev. discussions. C. A. C. is funded by the NRF (South Africa). D 75, 062001 (2007); A. Balbi and C. Quercellini, arXiv:0704.235; H. Zhang et al., arXiv:0705.4409. [19] J.-P. Uzan, Rev. Mod. Phys. 75, 403 (2003); J.-P. Uzan, arXiv:astro-ph/0409424; P. Molaro et al., Proc. IAU *[email protected] Symposium 198, 1 (2005). †[email protected] [20] L. Pasquini et al., The Messenger 122, 10 (2005). ‡[email protected] [21] I. M. Etherington, Philos. Mag. 15, 761 (1933); B. A. [1] G. F. R. Ellis, Q. J. R. Astron. Soc. 16, 254 (1975). Bassett and M. Kunz, Phys. Rev. D 69, 101305 (2004); [2] A. Riess et al., Astrophys. J. 607, 665 (2004); P. Astier J.-P. Uzan, N. Aghanim, and Y. Mellier, Phys. Rev. D 70, et al., Astron. Astrophys. 447, 31 (2006); D. Eisenstein 083533 (2004). et al., Astrophys. J. 633, 560 (2005); L. Fu et al., [22] J.-P. Uzan, F. Bernardeau, and Y. Mellier, Phys. Rev. D 77, arXiv:0712.0884. 021301(R) (2008). [3] D. Spergel et al., Astrophys. J. Suppl. Ser. 148, 175 [23] A. Loeb, Astrophys. J. 499, L111 (1998). (2003).

191303-4