Arxiv:1904.12388V2 [Gr-Qc] 5 Oct 2019 Gianluca.Gregori@Physics.Ox.Ac.Uk Orsodn Uhr .Gregori G

Home , Quantum hydrodynamics

arXiv:1904.12388v2 [gr-qc] 5 Oct 2019 [email protected] orsodn uhr .Gregori G. author: Corresponding in xeiet aeas eosrtdta the that demonstrated interpreta- also physical have different Experiments a tion. with albeit Schrödinger usual equation, the from developed be can approaches as ( results mechanics same quantum orthodox the yields mechanics quantum Bohmian de the ( or interpretation alter- mechanics Broglie-Bohm an Bohmian however, called is, possibility, There native and particle. position the of the momentum measure simultaneously is to it because impossible is too), Copenhagen This the space (or mechanics. curved quantum orthodox interpretation) in to anymore true extended not be is path this can defined that well a concept along (a moves New- particle ordinary in a trajectory While theory the tonian follows. particle of a terms that path) in (or expressed usually is problem 1987 c,btas ytennlclfebc fteparticle’s ( the motion of own feedback its non-local on function the wave by mechan- also classical but ordinary in ics, the as by potential, only local not of determined action one, classical the from ferent novd(e .. Ref. e.g., remains motions (see reconcil- classical of of unsolved task equations the the with yet them but ing years, 100 nearly been for have known mechanics quantum governing equations The thsbe hw,i h o-eaiitccs,that case, non-relativistic the in shown, been has It yee sn L using 2019 Typeset 8, October version Draft ,weeytepril olw ent ah dif- path, definite a follows particle the whereby ), h aeeetmyas epeeti h aeo ninhomogeneous an of expans case the the Keywords: in in present acceleration be an also as may appear effect same and additional the matter include Univ of equations a for nature modified equations These local Friedmann’s of particles. cla version massive the th modified zero, in apply a We included derive metric. be to background can metric the effects on quantum transformation the conformal that show and proach, euea lentv nepeaino unu ehnc,based mechanics, quantum of interpretation alternative an use We A T oie remn qain i ofra om–D Broglie De – Bohm conformal via equations Friedmann Modified 1. E X omlg:ter omlg:dr atr—csooy aken dark cosmology: — matter dark cosmology: — theory cosmology: INTRODUCTION twocolumn 2 a-lnkIsiu ¨rKrpyi,Psfc 03 0 69 80, 39 10 Postfach für Kernphysik, Max-Planck-Institut 1 eateto hsc,Uiest fOfr,PrsRa,Ox Road, Parks Oxford, of University Physics, of Department ie n Muft and Hiley om1952 Bohm om1952 Bohm tl nAASTeX62 in style om1952 Bohm .Gregori, G. ,ide both indeed ), ( ; 1995 Bohm ). ) The )). 1 tal. et .Reville, B. ABSTRACT iietr,poie gi u aei ae oensure to taken is care invariance. due Lorentz again ad- parti- provided an the including term, for by ditive The account forces that external with acceptable. actions interaction for is cle’s true motion also of is equations same correct the to where fidohrie where otherwise, ified and ( utbemdfiaino h cinfrafe atceof particle free a for through action manner, mass the usual of the modification in achieved suitable particle’s is the This on potential motion. nonlocal including Bohm’s by of path effects classical the the of reinterpretation a mits ubro ieetapoce oitoueteeffects the introduce to approaches different of number case this (in scalar ( world equation a in of metric function the modify to ( presented. been for accounted is phys- non-locality sound full ( a if have interpretation can ical trajectories Broglie-Bohm de 2014 2004 ohv intr (+ signature have to noteaoeato i oie metric: modified a via action above the into sdsusdin discussed As 2 2005 2011 Haas Mahler n .Larder B. and eew dp n ftelte prahs(see approaches latter the of one adopt we Here h ffcso netra edcnb incorporated be can field external an of effects The g ; ; µν o eiw n eeo omls hc per- which formalism a develop and review) a for ) m arl 2007 Carroll oln 1992 Holland u , h pc-iemti w s nt,uls spec- unless units, use (we metric space-time the µ tal. et = ; 2 edleg Germany Heidelberg, 029 smto oteRobertson-Walker the to method is Cross dx o fteUies.W e that see We Universe. the of ion odO13U UK 3PU, OX1 ford 2016 expansion. em htrsl rmtenon- the from result that terms µ S Goldstein nteBhintaetr ap- trajectory Bohmian the on scleuto fmto i a via motion of equation ssical /ds 1 recnitn fsaa,spin- scalar, of consisting erse = n uvdsaetm ( space-time curved and ) tal. et ; .Etnin oqatmfluids quantum to Extensions ). − Zee , hjiadSoa 2006 Shojai and Shojai stedmnines4-momentum dimensionless the is − m c , Z = − ( 2014 2013 ( ergy , p ~ 1950 − ) h erci taken is metric The 1). = ). g ,fil hoy( theory field ), ,freape xlisa exploits example, for ), µν gravity ,teei osnl way single no is there ), u µ u 1 ν ,adi at any fact, in and ), ds u µ u µ htleads that ) aealso have ) g Dürr Dürr µν Carroll → tal. et al. et g ˜ (1) µν . 2 of electromagnetic and gravitational fields into the action. As an example, it is well known that a simple Exploiting our prior knowledge of the Bohm potential, modification of the Minkowski metric tensor: i.e., letting V = Vext +VB, and carrying out the standard substitution = ∂tS it is immediately observed that η η˜ = diag(1 + 2V/m, 1, 1, 1), (2) H − µν → µν − − − (9) reproduces (4) as expected. Thus, if we acknowledge for a scalar potential V , provides a way to reproduce the presence of the Bohm potential in a given system, Newton’s law in the non-relativistic limit. We will em- it is possible to include its effect in macroscopic systems ploy a similar approach below in order to introduce via an appropriate transformation of the metric. Recent quantum effects into more complex metrics. work that exploits this approach (Carroll 2005), tacitly assumes the equations already describes an ensemble or 2. QUANTUM EFFECTS AS A CONFORMAL cloud of particles. We attempt to make this concept TRANSFORMATION OF THE METRIC more precise in the next section. Writing Schrödinger’s equation, for a particle of mass m, in polar form: ψ = ψ exp (iS) with real function | | S(x,t), Bohm (1952) establishes the set of coupled equations 3. THE MANY-BODY VERSION OF THE BOHM ∂ ψ 2 + ∇ ~v ψ 2 = 0 (3) POTENTIAL t| | · | | and While, so far, we have considered the motion of a single test particle, it is clear that we cannot neglect the effects 1 1 2 ψ ∂ S + ( S)2 ∇ | | + V =0, (4) of all the other particles in the system. A full theoret- t 2m ∇ − 2m ψ ext | | ical description of this quantum system must take into where we have set ~p ∇S = m~v = m d~x/dt, and V account effects such as non-locality and correlations on ≡ ext is the external potential acting on the particle. The for- all the scales of interest. A continuum description of mer equation corresponds to conservation of probability such quantum gas is often more convenient, and it usu- density, while the latter is recognised as the Hamilton- ally done in terms of an average density, by replacing Jacobi equation. The analogy with classical mechanics the product of many-body wavefunctions with an ap- is self-evident on defining the so-called Bohm potential propriate density operator into the relevant equations of motion. In previous work (Carroll 2005), the transi- 1 2 ψ tion from a microscopic to a continuum (macroscopic) VB = ∇ | |. (5) −2m ψ description of the Bohm’s potential was implicitly as- | | sumed. Thus, in Bohm’s interpretation, the particle moves as In order to justify a continuum approach in the Bohm guided by its own wave function. In the limit of VB ≪ picture more rigorously, let us consider an isolated quan- Vext the trajectory is indistinguishable from its classical tum system of N particles of mass m, described by a one. many-body wavefunction ψ = ψ(x , x , ..., x ), where As mentioned earlier, an alternative viewpoint is to 1 2 N x is the spatial coordinate of the ith particle. For sim- rewrite (1) in the Minkowski metric as i plicity, in the following, we take Vext = 0. The N-body Bohm potential can be written as (Bohm 1952): = m g˜ uµuν ds, (6) S − µν Z p where following Zee (2013) we set 2 (N) 1 i ψ 2V V = ∇ | |, (10) g˜ = 1+ η . (7) −2m ψ µν m µν i=1,N | | X To order d~x/dt 2, ignoring additive constants, the | | above action gives the classical Lagrangian where is the gradient with respect to the ith ∇i m d~x 2 particle coordinates. It is first assumed that the L V, (8) many-body wavefunction can be written as a prod- ≃ 2 dt − uct of single-particle wavefunctions: ψ(x1, x2, ..., xN )= th with associated Hamiltonian ϕ1(x1)ϕ2(x2)...ϕN (xN ), where ϕi(xi) is the i parti- 2 cle wavefunction. The expectation value of the Bohm m d~x ~p 2 = + V = | | + V. (9) potential is given by H 2 dt 2m 3

V = dx dx ... dx V (N) ψ 2 = dx V (N)(x ) ϕ (x ) 2 = dx V (N)(x) ϕ (x) 2, (11) h i 1 2 N | | i i i | i i | i | i | Z Z Z i=1X,N Z Z i=1X,N

where, is the single-particle Bohm potential. The total particle 2 2 density is n(x) = i ϕi = i ni, where ni is the (N) 1 ϕi(x) | | V (x)= ∇ | | (12) probability distribution for the ith particle. This gives i −2m ϕ (x) P P | i |

2 2 1 √ni 1 √n V = dx ni ∇ dx n(x) ∇ , (13) h i −2m √ni ≈ −2m √n Z i=1X,N Z

where the last step applies when each particle’s wave- candidates. With a view to exploring the Bohm poten- function significantly extends over their inter-particle tial associated with such dark matter candidates, we separation, and so the single-particle’s amplitudes are apply the previous analysis to the Klein-Gordon equa- assumed to be the same throughout space (that is, tion (Nikolić2005, 2007; Shojai and Shojai 2006; Carroll ϕ (x) ϕ (x) ). This is the linearization approxi- 2005) | i | ≡ | 0 | mation of quantum hydrodynamics (Manfredi and Haas ✷ ψ + m2ψ = 0 (15) 2001; Mitcha et al. 2015; Moldabekov et al. 2015), and where ✷ = gµν is the d’Alembertian, and the it strictly applies when the number density of particles ∇µ∇ν ∇µ covariant derivative with respect to the metric g . In varies smoothly in space. However, it can also be shown µν the Minkowski metric ✷ = ∂2 2. (Mitcha et al. 2015; Moldabekov et al. 2015) that even t − ∇ Following de Broglie (1960), writing the Klein-Gordon in presence of strong inter-particle correlations, equation equation in polar form, and separating out the real part (13) continues to apply approximately, to order (1). O gives By comparing equations (11) and (13) we see that the many-body Bohm potential can be thus written as a 2V gµν S S = 1+ B m2 2 , (16) function of a single spatial coordinate as ∇µ ∇ν m ≡M 1 2n1/2 V (x)= ∇ . (14) where VB is now in covariant form −2m n1/2 1 ✷ ψ The many-body quantum potential discussed here and V = | | . (17) B 2m ψ its application to extended systems of electrons and | | ions has recently been tested against more established The question, however, arises as Eq. (16) can not be solutions of the Schödinger’s equation based on den- adopted as the relativistic Hamilton-Jacobi equation sity function theory (DFT) (Hohenberg & Kohn 1964; (Holland 1992; Carroll 2005). The fact that is not Kohn & Sham 1965). In fact, the quantum poten- M positive definite is a serious concern, and furthermore tial approach shows excellent agreement with DFT, (16) leads to tachyonic solutions, which we must exclude with substantial reduction in computational speed on physical grounds – although tachyonic Bohm theo- (Larder et al. 2019). ries have been discussed in the literature (González-D´ıaz The complexity of quantum mechanics is not circum- 2004). The correct approach is that suggested by vented in Bohm’s interpretation; the particles’ trajec- Shojai and Shojai (2006), giving tories are determined by this nonlocal potential, which depends on the particles’ wave function itself. However, 2V 2 = m2 exp B . (18) what Bohm’s description offers is a recipe that allows us M m to obtain quantum corrections to the classical equation. This reduces to the previous case in the weak field limit. 4. THE RELATIVISTIC VERSION OF THE BOHM To see better how the exponential form for the ef- POTENTIAL fective mass, , arises, we follow the approach of M Extensions to cosmological models also requires a rela- Shojai and Shojai (2006). By minimizing the relativis- tivistic and covariant treatment for relevant dark matter tic action, we easily obtain the equation of motion for a 4 particle of variable mass = (x) in the Minkowski approach of Shojai and Shojai (2006) generalises this M M metric, result to ensure it satisfies the correct non-relativistic du µ = (η u u ) ∂ν , (19) limit. Identifying the four momentum Pµ = µS, the µν µ ν µν 2 ∇ M dτ − M relativistic energy equation g PµPν = m is thus re- 2 µ where dτ = dxµdx is the proper time. By taking the covered as the weak quantum potential limit of a more µν 2 non-relativistic limit (i.e., assuming small velocities and general equationg ˜ ˜ µS ˜ ν S = m , where ˜ µ is the ∇ ∇ ∇ t = τ), the space components of the equation of motion covariant derivative with respect to the metricg ˜µν . can be written as Semi-relativistic approaches, where the quantum po-

2 tential is simply added to the equation of motions with- d ~x out a conformal transformation, have also been studied 2 = , (20) M dt −∇M (Das 2014), and lead to solutions that are different from or, alternatively, as standard cosmology (Ali and Das 2015). Following the same arguments as earlier, we can ex- d2~x tend this previous result to the many-body system, such m = m ln M , (21) dt2 −∇ m that the conformal transformation we should apply has the form with m the bare mass, and m is some arbitrary mass scale. Now, if we require that eq. (21) describes the mo- g˜ =eQ g , (24) tion of a quantum particle of mass m, then consistency µν µν with the Hamiltonian (9), implies where 1 2 ψ = m exp ∇ | | . (22) M −2m2 ψ 1 ✷ n1/2 | | Q = . (25) m2 n1/2 The relativistic generalization would then lead to Thus, given Q, quantum effects on classical bodies moving in space-time can be calculated via this prescrip- 1 ✷2 ψ tion. In the next section, this is applied to the Einstein = m exp | | , (23) M 2m2 ψ field equations. | | which is the same as (18) if we choose the mass scale m m. 5. EXTENSION TO COSMOLOGY ≡ As further pointed out by de Broglie (1960) (chap. 10), the inclusion of quantum effects are entirely equiva- We proceed by setting gµν to be the Robertson-Walker lent to the change of the space-time metric g g˜ = metric of an expanding Universe. The line element in µν → µν ( 2/m2)g , that is, a conformal transformation. The such a metric is given by M µν

dr2 ds2 = dt2 R2(t) + r2(dθ2 + sin2 θ dφ2) dt2 R2(t) dΩ2, (26) − 1 κr2 ≡ − − where R(t) is the scale factor, set to unity at present Assuming a matter-like field n 1/R3, the quantum ∼ time, and κ is the curvature parameter. This allows potential has the form some simple scaling relations for the quantum potential 3 1 Q to be obtained. In the homogeneous approximation, Q = RR/¨ R˙ 2 + H2/m2 , (28) −2 2 the density is a function of time only. For the Robertson- Walker metric the d’Alembertian operator is where one can identify q = RR/¨ R˙ 2 as the usual de- d − celeration parameter and H = R/R˙ is the Hubble parameter. We now consider the conformal transformationg ˜µν = Q e gµν on the Robertson-Walker metric. From the above consideration, we see that the conformal factor has a 1 ∂ ∂ ✷ n1/2 = R3 n1/2 . (27) non-trivial dependence on time. In fact, its scaling with R3 ∂t ∂t 5 time has been discussed in the literature (Canuto et al. stant. The second term, instead, is the vacuum energy 1977; Maeder 2017a) – while it is required by all scale- contribution that arise from field theory. Estimating invariant gravity theories, it has been only justified via the actual value for ρ requires further approxima- h vaci arguments needed to match the present value for the cos- tions. In Minkowski space time, such estimate can be mological constant (Canuto and Lee 1977). The modi- obtained by summing up the zero point contributions fied Einstein field equation now reads as (Canuto et al. for each normal mode of a scalar field of mass m; we 1977; Maeder 2017a) have (Weinberg 1989), 1 ˜ g˜ ˜ =8πG ˜ +Λ˜g , (29) Rµν − 2 µν R Tµν µν kc 4πk2dk 1 k4 ˜ 2 2 c where µν is the Ricci tensor with respect to the modi- ρvac = 3 k + m 2 , (30) R h i 0 (2π) 2 ≈ 16π fied metricg ˜µν , ˜ is the Ricci scalar, ˜µν is the energy- Z R T p momentum tensor, and Λ is the effective cosmological constant (Weinberg 1989). Since the left hand side of eq. (29) is manifestly scale-invariant, the same must be with k m some wavelength cut-off which is usu- c ≫ true for the right hand side, in particular the energy- ally taken to be set by the reduced Planck length, −1/2 momentum tensor (Canuto et al. 1977; Maeder 2017a), that is kc = 1/ℓP = (8πG) . Since cosmological ˜ that is, µν = µν . The stress tensor, in general, con- observations tells us that Λ 8πG ρvac (Weinberg T T (M) (Q) ≪ h i tains two terms (Carroll 2007) ˜ = ˜µν + ˜µν . 1989), then it must be that Λ and 8πG ρ can- Tµν T T E h vaci The first one is related to the matter contribution, cels out each other to a very high degree of accu- (M) µν = (p + ρ)u u pg (for a perfect fluid), where racy, and, in fact, Λ = 0 may be a good approxima- T µ ν − µν p is the pressure and ρ is the energy density of the mat- tion (Weinberg 1989; Sahni and Krasiski 2008; Coleman ter. A consequence of the scale invariance of the energy- 1988; Hawking 1984). The latter statement can be jus- momentum tensor is that both the pressure and energy tified from a thermodynamic argument (Volovik 2006; density of matter are not scale invariant, as discussed in Bucher and Spergel 1999). If we consider a vacuum Uni- Maeder (2017a). verse, the total pressure is p = Λ. Interestingly, the − The second term of the energy-momentum tensor, same relation applies for a quantum fluid at zero tem- (Q) ˜µν , instead, arises from the energy density of the perature, i.e., in the absence of any finite temperature T quantum potential. This is a quantum gravity modi- excitations. If there are no external forces acting on the fication, implying that quantum effects back react on fluid, then it must be p =0andsoΛ=0(Volovik 2001; the metric. As these terms appear as corrections of or- Volovik 2006). While Λ = 0 is still hypothetical, it re- der (m2/m2 ) (Pinto-Neto and Struyve 2018), where mains plausible in different approaches of quantum grav- O P m = 2.4 1018 GeV is the reduced Planck mass, we ity (Coleman 1988; Hawking 1984; Barrow and Shaw P × can ignore them for particle masses much smaller than 2011; Ng and van Dam 1990). mP . The solution of Einstein field equations (29) gives The effective cosmological constant contains two modified Friedmann equations (Canuto et al. 1977; terms, Λ = Λ +8πG ρ . The first one is a con- Maeder 2017a): E h vaci

8πG κ R˙ 2 λ˙ R˙ λ˙ 2 Λλ2 ρ = + +2 + , (31) 3 R2 R2 λR λ2 − 3 κ R¨ R˙ 2 λ¨ λ˙ R˙ λ˙ 2 8πGp = +2 + +2 +4 Λλ2, (32) − R2 R R2 λ λR − λ2 −

where λ2 =eQ. It is easy to see that the modiﬁed Fried- mann’s equations, (31) and (32), reduce to the usual ones if λ = 1. This is perhaps more transparent if we R¨ 4πG Λ˜ λ˙ R˙ λ˙ 2 λ¨ = (ρ +3p)+ + , (34) combine the equations to read R − 3 3 − λR λ2 − λ R˙ 2 8πG Λ˜ λ˙ R˙ λ˙ 2 κ = ρ + 2 , (33) with Λ=Λ˜ λ2. R2 3 3 − λR − λ2 − R2 6

We notice that similar equations also appear in quan- with ρc the critical density of matter. Observations tum gravity theories where the infinite degrees of free- suggest Ω 0.3 at the present time, accounting m ∼ dom associated to the solution of the Wheeler-DeWitt for both luminous and dark matter (Riess et al. 1998; equation are collapsed under symmetry constraints, a Perlmutter et al. 1999; Goobar and Leibundgut 2011; so called gravity-matter mini-superspace (Vink 1992; Krauss and Chaboyer 2003). Pinto-Neto and Struyve 2018). The simplest mini- The implication of Eq. (38) in the context of confor- superspace model is described by the metric mal gravity has been discussed in detail in Ref. Maeder (2017a,b). Here we consider for simplicity a Universe 2 2 2 2 2 ds = (t) dt R (t) dΩ , (35) where, in the weak field limit, Λ˜ Λ = 0, and similarly N − ∼ the space curvature term is set to zero (κ = 0). The lat- where (t) is the lapse function (which sets the time N ter approximation is consistent with both observational gauge). Quantum effects are then introduced via the limits (Goobar and Leibundgut 2011) and existing in- Wheeler-DeWitt equation, leading to a modified equa- flationary models (Guth 1981). Proceeding with this tion of motion of the form (Pinto-Neto and Struyve assumption that Ω , Ω 1, Eq. (38) can be conve- 2018) Λ κ ≪ niently re-written in the form R˙ 2 Λ κ = + QWDW, (36) 2 2 2R2 6 − 2R2 ˙ N Q where, for simplicity, we have assumed an empty Uni- Ωm = 1+ , (43) 2H ! verse, and QWDW is the Wheeler-DeWitt quantum potential. Equation (36) can be re-cast in the same form implying that Q˙ < 0, as Ωm < 1. Using our previ- as (33) if we take (t) = 1, and we set N ous definition of Q, from Eq. (28), the above can be expressed as λ˙ R˙ 1 λ˙ 2 Λ(λ2 1) Q + − . (37) WDW ≡−λR − 2 λ2 6 d 1 4 q H2 = m2H 1 Ω1/2 , (44) dt 2 − d 3 − m However, Equation (33) is also reproduced if we take (t) = λ(t), and make a transformation R(t) a non-linear third-order differential equation in R. N → λ(t)R(t) in (35), meaning that quantum effects are Numerical integration of this equation requires three equivalent to a coordinate change of the metric, as noted boundary conditions, at least one of which is uncer- earlier. This also means that, at the level of our approx- tain. Nevertheless, we can still comment on the general imations, the back-reaction of quantum effects onto the behaviour. As the right hand side of this equation is metric is ignored. clearly positive, consistency requires that Equations (33) and (34) are very similar to those of q˙d standard cosmology, with some noticeable differences. < (qd + 1)(2qd 1) (45) Firstly, there is now a modified cosmological constant Λ˜ H − which also depends on the quantum potential. Secondly, indicating a number of different interesting regimes: there is an additional term that produces an accelera- i) q > 1/2: in this case, the Universe is decelerat- tion, as with dark energy, but it is not related to the d ing, and consistent solutions for both positive and properties of the vacuum. Instead, this new term de- negativeq ˙ can be found. Such large deceleration pends on the non-locality of quantum interactions of a d parameters are however disfavoured by current ob- Universe filled with matter. servations. Equation (33) can be rewritten in a more familiar form as ii) q < 1: i.e. the universe is accelerating rapidly. d − Ωm +ΩΛ +Ωκ +ΩQ =1, (38) q˙d again can have both positive and negative val- where, ues but as previously, this regime is also disfavoured by observations. 8πGρ ρ Ωm = 2 = , (39) iii) 1 qd 1/2: this is physically the most inter- 3H ρc − ≤ ≤ esting regime, and can be subdivided further into Λ˜ Ω = , (40) two sub-cases: one accelerating (qd < 0), one de- Λ 3H2 κ celerating (qd > 0). However,q ˙d < 0 for both Ωκ = , (41) implying that a transition from case (iii) to case −R2H2 (i) is not possible. A transition to case (ii) cannot Q˙ Q˙ 2 Ω = , (42) be ruled out. Q −H − 4H2 7

Current observations favour the latter scenario, with vations in this interpretation as it gives the present time qd < 0 (Riess et al. 1998; Perlmutter et al. 1999), a very special place in assigning the mass of such a parti- although this conclusion is not universally accepted cle, set only by the current value of the Hubble constant, (Colin et al. 2018). Proceeding on the assumption that an example of fine-tuning. the universe is in fact accelerating, we conclude that the term ΩQ behaves as an acceleration that opposes 6. BEYOND THE ISOTROPIC AND gravity. It is equivalent to what is usually referred to as HOMOGENEOUS UNIVERSE dark energy (Riess et al. 1998; Perlmutter et al. 1999). The model presented so far is still quite ideal. In par- In this context, as pointed out earlier, it represents the ticular, it is assumed that the Universe is homogeneous background energy associated to the non-local quantum and isotropic. However, it is clear from observations nature of matter. that the Universe, mostly at smaller scales, is neither Using dimensional analysis on equation (44), one can homogeneous nor isotropic. Locally, it consists of grav- 0 infer the present value of the dark energy density ΩDE = itationally bound structures, such as cluster of galax- 0 2 2 −1 −1 ΩQ H0 /m (1), where H0 = 70 kms Mpc = ies, and voids. Since the equations of general relativity ∼ − ∼ O 2 10 33 eV (Weinberg 1989). This would imply the are non-linear, it is plausible that local density inhomo- × −33 existence of a field particle of mass m mg 2 10 geneities could produce significant changes to the cos- ≡ ∼ × eV. We note that this result is, in fact, analogous to what mological evolution as described by the standard Fried- obtained in quintessence (scalar-field) models of dark en- mann’s equations. To investigate the effects of an inho- ergy (Copeland et al. 2006), whereby the present value mogeneous distribution of mass, let us consider the case of the cosmological constant is explained by a scalar field of a spherically symmetric dust Universe (as seen from of mass mφ mg H0. The existence of such a field is our location at the center) described by the Lemaˆıtre- ≡ ∼ hypothetical, nevertheless the value given for its mass is Tolman-Bondi metric. To further simplify the analysis, consistent with the more stringent bounds on the gravi- as we did before for the case of the Robertson-Walker ton mass (Goldhaber 2010). There are, however, reser- metric, we take Λ = 0 and κ = 0. The line element is thus (Enqvist 2008):

ds2 = dt2 [∂ F (r, t)]2 dr2 + [F (r, t)]2 (dθ2 + sin2 θ dφ2), (46) − r where F (r, t) is some (generally unknown) function to the Robertson-Walker metric (for Λ = 0 and κ = 0) of the radial coordinate and time. To constrain the if f(r) = 0. model, we now assume that the function F (r, t) can As we have done before for the Robertson-Walker met- be parametrized in the following form: F (r, t) = ric, in the presence of non-local effects induced by the R(t)[r + f(r)], where R(t) is the same scale factor as in quantum nature of matter, we also need to include the the Robertson-Walker metric and f(r) is a function that conformal factor λ2 = eQ(r,t), where we have explicitly describe the departure from homogeneity. This reduces included both a radial and time dependence in the Bohm potential. With this metric, the solution of Einstein’s field equation gives a modified version of Equation (33):

R˙ 2 8πG (∂ Q)R˙ (∂ Q)R 2+4f ′(r) + 2[f ′(r)]2 [r + f(r)]f ′′(r) (∂ Q) (∂ Q)2 + 4(∂2Q) = ρ t t +1 + { − } r + r r . (47) R2 3 − R ˙ 3R2[r + f(r)][1 + f ′(r)]3 12R2[1 + f ′(r)]2 4R

If the conformal factor is independent of the radial co- 2001; Alnes et al. 2006), but in a way which is not trivial ordinate, equation (43) is recovered. However, in an and depends on the exact structure of those perturba- inhomogeneous Universe, the conformal factor will, in tions. general, depend on density and so ﬂuctuations also act toward changing the local expansion rate (Buchert 2000, 8

7. DARK ENERGY IN AN INHOMOGENEOUS axion” is a well-motivated candidate for the dark mat- UNIVERSE ter (Dine and Fischler 1983; Abbott and Sikivie 1983; To simplify equation (47) still further, let us take Preskill et al. 1983). Because of its small coupling, the axion field decouples with the baryonic matter and be- Q(r, t) = Q(t) [1 + ℓq(r)] and assume that fluctuations in the metric are uncorrelated with fluctuations haves as a Bose-Einstein condensate. The property of axionic dark matter (with ma of the quantum potential. Such fluctuations are as- − ∼ 8 10 23 eV giving the best fit to observations) has sociated to random density perturbations. Locally, × however, there are inhomogeneities and anisotropies been investigated in recent structure formation simu- in the distribution of matter, and a transition to statis- lations (Schive et al. 2014). Interestingly, those simu- tical homogeneity is reached only at sufficiently large lations seem to explain the mass distribution around scales. In standard cosmology this scale is taken dwarf galaxies, which instead cannot be reproduced cor- to be 100 Mpc (Hogg et al. 2005), consistent with rectly by structure formation calculations where dark ∼ baryon oscillations from the Sloan Digital Sky Survey matter is made by cold non-interacting classical par- (Eisenstein et al. 2005). On the other hand, there is ticles (Schive et al. 2014). A characteristic of axionic still the possibility that inhomogeneities of the mat- dark matter simulations is a fuzziness of the mass dis- ter distribution persist at distances exceeding 300 tribution. This is set by the axion Compton wavelength ∼ 1/2 2/3 Mpc (Colin et al. 2018), which has then implications on λCa = (2π/maTa) , where Ta = [2π/ζ(3/2)ma] na the determination of the Hubble constant (Colin et al. is the critical temperature of the Bose-Einstein con- 2018; Heß and Kitaura 2016). Thus when averag- densate (Pitaevskii and Stringari 2016), ζ(3/2) = 2.61 ing over a volume larger than this, we must have is the Riemann zeta function and na = ρDM /ma is 2 ′ the present-day number density of axions, with ρDM = ℓq(r) = ∂rℓq(r) = ∂r ℓq(r) = f(r) = f (r) = 0. h i h i h i h i h i 9.6 10−12 eV−4 the co-moving dark matter mass den- Finally, we take the spatial fluctuations to be small in × amplitude, that is, f ′(r) < 1 and ℓ (r) < 1. Thus, sity (Patrignani et al. 2016). | | | q | the only non vanishing terms in equation (47) result in We do not know much about the size of the density perturbations. If we assume they are typically of 2 2 2 2 2 the same order as the fluctuations in temperature seen Q˙ ℓ Q (∂ ℓ ) Q˙ − Ω h qi + h r q i = 1+ . (48) by the cosmic microwave background, then δ 10 5 m 2 2 2 ∼ − 4H 12R H 2H ! (Wright 2004). Setting L λ , the equation for the ∼ Ca dark matter density at the present epoch becomes Hence, as discussed earlier, using dimensional analysis, 2 2 3 2 Q H /m , Q˙ H /m , and ∂ ℓ δ/L, where 2/3 | | ≈ | | ≈ | r q| ≈ δ2H2ρ L is the scale of the fluctuations of the field m, and Ω0 0 DM . (50) DE ≈ 14/3 δ = ℓ2 1/2. Inhomogeneities dominates the total dark 12 ξ(3/2) ma h qi energy density only if L δ/mR. If we take this to be ≪ Requiring Ω0 0.7, this gives m 10−18 eV. Such the case, the dark energy density has a much simplified DE ≈ a ∼ relation (and even smaller) axion masses are, in fact, what is δ2H2 needed for small scale structure formation (Schive et al. Ω . (49) DE ≈ 12L2R2m4 2014) and they may be detected by oscillations in pulsar timings (De Martino et al. 2017). Moreover, small ul- −21 −17 One possibility we explore now is that this hypo- tralight scalars in the range 10 eV . ma . 10 eV thetical field m is due to an ultralight axion. The have been invoked to explain the cosmological ori- et al. axion is a pseudo Nambu-Goldstone boson of the bro- gin of magnetic fields (Choi 2018) (see also et al. ken U(1) Peccei-Quinn symmetry (Peccei and Quinn Ref. Miniati (2018)). 1977; Wilczek 1978), which was introduced to explain the absence of CP violation in strong interac- 8. SUMMARY AND CONCLUSIONS tions. Extensions of this model to axion-like particles In this work we have argued that the interpretation of (Jaeckel and Ringwald 2010) can admit the existence of the modified Friedmann’s equation in conformal grav- ultralight scalars or pseudo-scalars. The axion remains ity including inhomogeneities in the matter is perhaps effectively massless until the Universe cools below some clarified if the conformal factor is determined by the critical temperature, and after that it acquires a mass quantum potential of Bohmian mechanics. In partic- m and starts oscillating with wavelength 1/m . ular, these results would suggest that the acceleration a ∝ a The energy density of these oscillations is of order the term in the expansion of the Universe need not be set critical density of the Universe, hence the “invisible by the vacuum energy but rather by the non-locality of 9 matter. This is by no means a solution of the cosmo- AUTHOR CONTRIBUTIONS logical constant problem, as the solution is reliant on a This project was conceived by G.G. Calculations were term which requires the existence of as yet undiscovered performed by B.R. and G.G., with B.L. contributing to particle fields. the extension to many-body systems. The paper was written by all the authors. ACKNOWLEDGEMENTS The authors would like to thank Prof Subir Sarkar (University of Oxford) for stimulating discussions related to this work. This research was supported in part by the Engineering and Physical Sciences Re- DATA AVAILABILITY search Council, UK (EP/M022331/1, EP/N014472/1, All data that support the findings of this study are EP/N002644/1, EP/P010059/1). available from the authors upon request.

REFERENCES

Hiley, B.J. & Muft, A.A., Zee, A., Einstein Gravity in a Nutshell (Princeton, 2013) The ontological interpretation of quantum field theory applied inDürr, a cosmological D. & Teufel, context S., Bohmian Mechanics (Springer, 2009) (Springer, 1995) pp. 141–156 Manfredi, G. & Haas, F., Self-consistent fluid model for a Bohm, D., A Suggested Interpretation of the Quantum quantum electron gas, Theory of “Hidden” Variables. I, Physical Review 85, Physical Review B 64, 075316 (2001) 166 (1952) Mitcha, D., Graziani, F. & Bonitz, M., Quantum Bohm, D., Hiley B., & Kaloyerou P., An ontological basis Hydrodynamics for Plasmas – a ThomasFermi Theory for the quantum theory, Physics Reports 144, 321 (1987) Perspective, Mahler, D.H., Rozema, L., Fisher, K., Vermeyden, L., Contributions in Plasma Physics 55, 437–443 (2015) Resch, K.J., Wiseman, H.M., & Steinberg, A., Moldabekov, Zh. A., Bonitz, M. & Ramazanov, T. S., Experimental nonlocal and surreal Bohmian trajectories, Theoretical foundations of quantum hydrodynamics for Science Advances 2, e1501466 (2016) plasmas, Phys. Plasmas 25, 031903 (2018) Haas, F., Quantum Plasmas: An Hydrodynamic Approach Hohenberg, P. & Kohn, W., Inhomogeneous electron gas, (Springer, 2011) Physical Review 136, B864 (1964). Cross, J.E., Reville, B., & Gregori, G., Scaling of Kohn, W. & Sham, L. J., Self-consistent equations magneto-quantum-radiative hydrodynamic equations: including exchange and correlation effects, Physical from laser-produced plasmas to astrophysics, Review 140, A1133 (1965). The Astrophysical Journal 795, 59 (2014) Larder, B., Gericke, D., Richardson, S., et al. Fast Dürr, D., Goldstein, S., Tumulka, R., & Zangh`ı, N., Non-Adiabatic Dynamics of Many-Body Quantum Bohmian Mechanics and Quantum Field Theory, Systems, Science Advances (accepted for publication), Physical Review Letters 93, 090402 (2004) arXiv:1811.08161 (2019). Carroll, R., Metric fluctuations, entropy, and the Nikolić, H., Relativistic Quantum Mechanics and the Wheeler-deWitt equation, Bohmian Interpretation, Theoretical and Mathematical Physics 152, 904 (2007) Foundations of Physics Letters 18, 549 (2005) Dürr, D., Goldstein, S., Norsen, T., Struyve, W., & Nikolić, H., Bohmian mechanics in relativistic quantum Zangh`ı, N., Can Bohmian mechanics be made mechanics, quantum field theory and string theory, relativistic? Proc. R. Soc. A 470, 20130699 (2014) Journal of Physics: Conference Series 67, 012035 (2007) Holland, P.R., The Dirac equation in the de Broglie-Bohm González-D´ıaz, P.F., Subquantum dark energy, theory of motion, Foundations of Physics 22, 1287 (1992) Physical Review D 69, 103512 (2004) Shojai, A. & Shojai, F., About Some Problems Raised by de Broglie, L., Nonlinear wave mechanics: A causal the Relativistic Form of De-BroglieBohm Theory of Pilot interpretation (Elsevier, Amsterdam, 1960) Wave, Physica Scripta 64, 413 (2006) Das, S., Quantum Raychaudhuri equation, Carroll, R., Fluctuations, gravity, and the quantum Physical Review D 89, 084068 (2014) potential, arXiv preprint gr-qc/0501045 (2005) Ali, A. & Das, S., Cosmology from quantum potential, Goldstein, H., Classical Mechanics (Addison-Wesley, 1950) Physics Letters B 741, 276 (2015) 10

Canuto, V., Hsieh, S.H., & Adams, P.J., Scale-Covariant Perlmutter, S., Aldering, G., Goldhaber, G., Knop, R.A., Theory of Gravitation and Astrophysical Applications, Nugent, P., Castro, P.G., Deustua S., Fabbro, S., Physical Review Letters 39, 429 (1977) Goobar, A., Groom, D.E., Hook, I.M., Kim, A.G., Kim, Maeder, A., An alternative to the ΛCDM model: the case M.Y., Lee, J.C., Nunes, N.J., Pain, R., Pennypacker, of scale invariance, C.R., Quimby, R., Lidman, C., Ellis, R.S., Irwin, M., The Astrophysical Journal 834, 194 (2017a) McMahon, R.G., RuizLapuente, P., Walton, N., Canuto, V., & Lee, J., The cosmological constant, broken Schaefer, B., Boyle, B.J., Filippenko, A.V., Matheson, T., gauge theories and 3 K black-body radiation, Fruchter, A.S., Panagia, N., Newberg, H.J.M., Couch, Physics Letters B 72, 281 (1977) W.J., & The Supernova Cosmology Project, Measurements of Ω and Λ from 42 High-Redshift Pinto-Neto, N., & Struyve, W., Bohmian quantum gravity Supernovae, The Astrophysical Journal 517, 565 (1999) and cosmology, arXiv:1801.03353 (2018) Goobar, A., & Leibundgut, B., Supernova Cosmology: Weinberg, S., The cosmological constant problem, Legacy and Future, Reviews of Modern Physics 61, 1 (1989) Annual Review of Nuclear and Particle Science 61, 251 (2011) Sahni, V. & Krasiński, A., Republication of: The Krauss, L.M., & Chaboyer, B., Age Estimates of Globular cosmological constant and the theory of elementary Clusters in the Milky Way: Constraints on Cosmology, particles (By Ya.B. Zeldovich), Science 299, 65 (2003) General Relativity and Gravitation 40, 1557 (2008) Maeder, A., Dynamical Effects of the Scale Invariance of Coleman, S., Why there is nothing rather than something: the Empty Space: The Fall of Dark Matter? A theory of the cosmological constant, The Astrophysical Journal 849, 158 (2017b) Nuclear Physics B 310, 643 (1988) Guth, A.H., Inflationary universe: A possible solution to Hawking, S., The cosmological constant is probably zero, the horizon and flatness problems, PhRvD 23, 347 (1981) Physics Letters B 134, 403 (1984) Colin, J., Mohayaee, R., Rameez, M., & Sarkar, S., Volovik, G.E., From Quantum Hydrodynamics to Quantum Apparent cosmic acceleration due to local bulk flow, Gravity, arXiv:gr-qc/0612134 (2006) arXiv:1808.04597 (2018) Bucher, M., & Spergel, D.M., Is the dark matter a solid? Copeland, E.J., Sami, M., & Tsujikawa, S., Dynamics of 60 Phys. Rev. D , 043505 (1999) dark energy, Volovik, G.E., Superfluid analogies of cosmological International Journal of Modern Physics D 15, 1753 (2006) 351 phenomena, Physics Reports , 195 (2001) Goldhaber, A.S., & Nieto, M.M., Photon and graviton Barrow, J.D., & Shaw, D.J., New Solution of the mass limits, Review of Modern Physics 82, 939 (2010) Cosmological Constant Problems, Enqvist, K., LemaitreTolmanBondi model and accelerating 106 Physical Review Letters , 101302 (2011) expansion, Ng, Y.J., & van Dam, H., Possible solution to the General Relativity and Gravitation 40, 451 (2008) cosmological-constant problem, Buchert, T., On Average Properties of Inhomogeneous Phys. Rev. Lett. 65, 1972 (1990) Fluids in General Relativity: Dust Cosmologies, Vink, J.C., Quantum potential interpretation of the wave General Relativity and Gravitation 32, 105 (2000) function of the universe, Phys. Lett. B 369, 707–728 Buchert, T., On Average Properties of Inhomogeneous (1992) Fluids in General Relativity: Perfect Fluid Cosmologies, Riess, A.G., Filippenko, A.V., Challis, P., Clocchiatti, A., General Relativity and Gravitation 33, 1381 (2001) Diercks, A., Garnavich, P.M., Gilliland, R.L., Hogan, Alnes, H., Amarzguioui, M., & Grøn, O., Inhomogeneous C.J., Jha, S., Kirshner, R.P., Leibundgut, B., Phillips, alternative to dark energy? M.M., Reiss, D., Schmidt, B.P., Schommer, R.A., Smith, Phys. Rev. D 73, 083519 (2006) R.C., Spyromilio, J., Stubbs, C., Suntzeff, N.B., & Hogg, D.W., Eisenstein, D.J., Blanton, M.R., Bahcall, Tonry, J., Observational Evidence from Supernovae for N.A., Brinkmann, J., Gunn, J.E., & Schneider, D.P., an Accelerating Universe and a Cosmological Constant, Cosmic Homogeneity Demonstrated with Luminous Red The Astronomical Journal 116, 1009 (1998) Galaxies, The Astrophysical Journal 624, 54 (2005) 11

Eisenstein, D.J., Zehavi, I., Hogg, D.W., Scoccimarro, R., Jaeckel, J., & Ringwald, A., The Low-Energy Frontier of Blanton, M.R., Nichol, R.C., Scranton, R., Seo, H.-J., Particle Physics, 60 Tegmark, M., Zheng, Z., Anderson, S.F., Annis, J., Annu. Rev. Nucl. Part. Sci. , 405 (2010) Dine, M., & Fischler, W., The not-so-harmless axion, Bahcall, N., Brinkmann, J., Burles, S., Castander, F.J., Physics Letters B 120, 137 (1983) Connolly, A., Csabai, I., Doi, M., Fukugita, M., Frieman, Abbott, L., & Sikivie, P., A cosmological bound on the J.A., Glazebrook, K., Gunn, J.E., Hendry, J.S., invisible axion, Physics Letters B 120, 133 (1983) Hennessy, G., Ivezić, Z., Kent, S., Knapp, G.R., Lin, H., Preskill, J., Wise, M.B., & Wilczek, F., Cosmology of the 120 Loh, Y.-S., Lupton, R.H., Margon, B., McKay, T.A., invisible axion, Physics Letters B , 127 (1983) Schive, H.-Y., Chiueh, T., & Broadhurst, T., Cosmic Meiksin, A., Munn, J.A., Pope, A., Richmond, M.W., structure as the quantum interference of a coherent dark Schlegel, D., Schneider, D.P., Shimasaku, K., wave, Nature Physics 10, 496 (2014) Stoughton, C., Strauss, M.A., SubbaRao, M., Szalay, Pitaevskii, L., & Stringari, S., Bose-Einstein Condensation A.S., Szapudi, I., Tucker, D.L., Yanny, B., & York, D.G., and Superfluidity (Oxford University Press, 2016) Detection of the Baryon Acoustic Peak in the Large-Scale Patrignani, C., & Particle Data Group, Review of Particle 40 Correlation Function of SDSS Luminous Red Galaxies, Physics, Chinese Physics C , 100001 (2016) Wright, E.L., Theoretical Overview of Cosmic Microwave ApJ 633, 560 (2005) Background Anisotropy, astro-ph/0305591 Heß, S., & Kitaura, F.-S., Cosmic flows and the expansion De Martino, I., Broadhurst, T., Tye, S.-H.H., Chiueh, T., of the Schive, H.-Y., & Lazkoz R., Recognizing Axionic Dark local Universe from non-linear phasespace reconstructions, Matter by Compton and de Broglie Scale Modulation of 119 Monthly Notices of the Royal Astronomical Society 456, 4247 (2016)Pulsar Timing, Phys. Rev. Lett. , 221103 (2017) Choi, K., Kim, H., & Sekiguchi, T., Late-Time Peccei, R., & Quinn, H., CP conservation in the presence Magnetogenesis Driven by Axionlike Particle Dark Matter of pseudoparticles, and a Dark Photon, Phys. Rev. Lett. 121, 031102 (2018) 38 Physical Review Letters , 1440 (1977) Miniati, F., Gregori, G., Reville, B., & Sarkar, S., Wilczek, F., Problem of Strong P and T Invariance in the Axion-Driven Cosmic Magnetogenesis during the QCD 121 Presence of Instantons, Phys Rev Lett 40, 279 (1978) Crossover, Phys. Rev. Lett. , 021301 (2018)