## Massive phonons and gravitational dynamics in a photon-fluid model

Francesco Marino

CNR-Istituto Nazionale di Ottica and INFN, Sez. di Firenze, Via Sansone 1, I-50019 Sesto Fiorentino (FI), Italy

(Dated: February 10, 2021)

We theoretically investigate the excitation dynamics in a photon-fluid with both local and nonlocal interactions. We show that the interplay between locality and an infinite-range nonlocality gives rise to a gapped Bogoliubov spectrum of elementary excitations which, at lower momenta, correspond to massive particles (phonons) with a relativistic energy-momentum relation. In this regime and in the presence of an inhomogeneous flow the density fluctuations are governed by the massive KleinGordon equation on the acoustic metric and thus propagate as massive scalar fields on a curved spacetime. We finally demonstrate that in the non-relativistic limit the phonon modes behave as self-gravitating quantum particles with an effective Schr¨dinger-Newton dynamics, although with o a finite-range gravitational interaction and a non-zero cosmological constant. Our photon-fluid represents a viable alternative to BEC models for 'emergent-gravity' scenarios and offers a promising setting for analogue simulations of semiclassical gravity and quantum gravity phenomenology.

## I. INTRODUCTION

Analogue gravity models provide a powerful test-bed for several aspects of classical and quantum field theories in curved spacetime [1-3]. The general idea is that under appropriate conditions the elementary excitations in condensed-matter systems evolve as fields on a curved spacetime induced by the medium. The paradigmatic example is provided by sound waves in an inhomogeneous flowing fluid [4, 5]. In spite of the fact that the background fluid is non-relativistic, the elementary excitations of the flow (phonons) experience a curved spacetime: their evolution is governed by the Klein-Gordon equation for a massless particle in a curved background, whose geometry is specified by a Lorentzian metric tensor (acoustic metric). As a result, the phonon dynamics exhibits an effective Lorentz invariance with the local speed of sound playing the role of the speed of light. The coefficients of the acoustic metric depend on the fluid density, which determines also the sound speed, and the flow velocity. Hence, by tailoring the properties of the flow it is possible to simulate gravitational spacetimes and related phenomena, such as e.g. Hawking radiation, superradiance and cosmological particle production.

light. Alaser beam propagating through a self-defocusing medium can be described in terms of a weakly interacting Bose gas, where the repulsive photon-photon interaction arises from a third-order nonlinearity [28, 29] and the propagation coordinate acts as an effective time variable. Recent experiments in these systems provided evidence of collective many-photon phenomena, such as superfluidity and its breakdown [30, 31] and nonequilibrium precondensation of classical waves [32]. In analogy with BECs, the collective excitations of the mean flow (i.e. small ripples of the transverse optical field) propagate according to the Bogoliubov dispersion relation [33], as recently demonstrated in thermo-optical [34] and Kerr media [35]. As a consequence, for the longer wavelengths a Lorentz invariant, phononic regime takes place where sound-like waves propagate with a constant speed determined by the photon fluid density. The latter is proportional to the optical intensity while the background flow velocity is controlled via the gradient of the phase profile. All these features make these systems particularly suitable for the realization of analogue gravity experiments [3644].

Analog-gravity scenarios have been proposed and realized in a variety of physical systems, including BoseEinstein condensates (BECs) [6, 7], surface waves [8], supefluid 3 He [9] and Fermi liquids [10], dielectrics, [11, 12], moving- and nonlinear-optical media [13-15] and excitonpolariton condensates [16, 17]. Signatures of the Hawking process have been reported in different setups [18-24] and in a recent experiment the observation of superradiance has been achieved [25].

As an alternative to the above systems, photon-fluids have recently attracted considerable attention. Photon fluids belong to the family of the so-called 'quantum fluids of light' [26], together with exciton-polariton and photon BECs [27]. While the last two are drivendissipative systems based on nonlinear optical cavities, photon fluids simply rely on the nonlinear propagation of

All the above systems are generally characterized by a gapless dispersion relation at small momenta [45] typical of massless collective excitations. Therefore most of the theoretical research in this area, and all ongoing experiments, have naturally focused on the simulation of massless fields propagating through a curved spacetime. A notable exception is the work by Visser and Weinfurtner [46] who first proposed a method to give rise to a spectrum of massive relativistic particles in a BEC system. The model describes a two-component BEC with an additional Raman coupling which deformes the spectrum of normal-mode excitations. Interestingly, in appropriate conditions one of the two phonon modes remains massless while the second acquires effective mass. Subsequent investigations by Liberati and co-workers introduced a modified BEC Hamiltonian with an U (1) symmetry-breaking term [47]. This modification provides a mass to the excitations and gives rise to a kind of

analogue gravitational dynamics. Remarkably, the gravitational potential is sourced by (a function of) the density distribution of the excitations which thus play the role of the matter in this system. These are important extensions with respect to usual analogue models, since such massive excitations fields could enable simulations of quantum-gravity phenomenology (e.g. Lorentz-violating dispersion relations [48]) and emergent-gravity scenarios [49].

In the following we consider a photon-fluid model for light propagating in a defocusing medium with both local and nonlocal optical nonlinearities. At difference with the purely local case, the first-order excitations satisfy a massive version of the Bogoliubov dispersion relation in a Bose gas, with the nonlocal term being the massgenerating mechanism. For the longer wavelengths, the spectrum approximates that of a massive particle with a relativistic energy-momentum relation and, in the presence of inhomogeneous flows, the density fluctuations are described by the massive Klein-Gordon equation on the acoustic metric, thus closely mimicking the propagation of massive scalar fields on a curved spacetime. Even more importantly, in the non-relativistic limit the phonons behaves as massive, self-gravitating quantum particles: their dynamics obeys the Schr¨dinger equation in a gravio tational potential whose source depends on the phononic mass density distribution via a modified Poisson equation. Unlike the Newtonian theory, we find that the range of the gravitational interaction is finite and a cosmological constant is also present. In spite of these significant differences with respect to standard gravity, such photon fluid nonetheless remains an interesting workbench for analogue simulations of semiclassical gravity scenarios. Most analogue-gravity models indeed are dealing with massless excitations that in Newtonian theory cannot act as sources of a gravitational field. This system is thus one of the very few in which a form of semiclassical gravitational dynamics can be shown to emerge.

The paper is organized as follows. In Sec. II, we introduce the modified Nonlinear Schr¨dinger Equao tion (NSE) with local and nonlocal nonlinearities and the related photon-fluid model. We then derive the Bogoliubov-de-Gennes equations governing the dynamics of the first-order fluctuations of the optical field and the corresponding dispersion relation for a generic nonlocal function. In Sec. III we focus on a thermo-optical nonlocal nonlinearity, showing that in the defocusing case the photon fluid is stable and allows for the propagation of massive phonons, while it undergoes a Jeans instability and supports tachyonic excitations in the focusing case. The rest of the paper is devoted to analyze the fully stable defocusing regime. In Sec. IV, we address the problem of inhomogenous flows and derive a massive Klein-Gordon equation on the acustic metric that will provide the basis for the subsequent discussion on the emergent gravitational scenario. In Sec. V, we introduce the Newtonian limit of the acoustic metric which allows us to identify the gravitational potential with inhomogeneities in the photon fluid density. We then derive the non-relativistic phonon dynamics from the Klein-Gordon equation for the optical field excitations. Finally, we show that a (modified) Poisson's equation for the potential is encoded in the backreaction equation describing the first corrections to the mean-field dynamics induced by the fluctuations. The conclusions are presented in Sec. VI.

## II. PHOTON FLUID MODEL AND ELEMENTARY EXCITATIONS

The propagation of a monochromatic optical beam oscillating at angular frequency ω in a 2D nonlinear medium can be described within the paraxial approximation in terms of the Nonlinear Schr¨dinger Equation o (NSE) [50]

$$\partial _ { z } E = \frac { i } { 2 k } \nabla ^ { 2 } E - i \frac { k } { n _ { 0 } } E \Delta n ( | E | ^ { 2 }, { \mathbf r }, z ) \quad \ \ ( 1 )$$

where E is the slowly varying envelope of the optical field, z is the propagation coordinate, k = 2 πn /λ 0 is the wavenumber, λ the vacuum wavelength and n 0 is the linear refractive index. The laplacian term ∇ 2 E defined with respect to the transverse coordinates r = ( x, y ) accounts for diffraction and ∆ n is the nonlinear optical response of the medium. For a local (Kerr) defocusing nonlinearity, ∆ n = n 2 | E | 2 with n 2 &gt; 0, Eq. (1) is formally identical to the 2D Gross-Pitaevskii equation for a dilute boson gas with repulsive contact interactions, where the optical field E corresponds to the complex order parameter and the intensity-dependent refractive index ∆ n provides the interaction potential. The dynamics takes place in the transverse plane ( x, y ) of the laser beam so that the propagation coordinate z plays the role of an effective time variable t = ( n /c z 0 ) , where c is speed of light in vacuum. We remark that the analogy between photon-fluids and condensates is here limited to the level of the mean-field evolution equations: as such, the system is purely classical and the optical field would correspond to the ground state wavefunction of a BEC at zero-temperature.

We consider an optical medium with both local and nonlocal third-order nonlinearities ∆ n ( | E | 2 , r , z ) = n 2 | E | 2 +ˆ n nl | E | 2 , where ˆ n nl is the convolution operator

$$\overset { \kappa } { \mathrm e } _ { \mathrm f } \equiv \gamma \left ( R * \right ) = \gamma \theta \int d \mathbf r ^ { \prime } d z ^ { \prime } R ( \mathbf r - \mathbf r ^ { \prime }, z - z ^ { \prime } ) \quad ( 2 )$$

with γ being a coefficient that depends on the specific nonlocal process and R ( r , z ) is the medium response function. In the following we take n 2 &gt; 0, since local repulsive interactions are required to observe a dynamically stable photon fluid on which sound waves can propagate, while θ = 1( -1) corresponds to a defocusing (focusing) nonlocal term, respectively.

The optical nonlocality originates from the fact that the nonlinear change in refractive index at any given position depends both on the local and on surrounding field

intensity through the convolution kernel R . Similar nonlinear responses arise in semiconductor materials with both Kerr and thermo-optical nonlinearities [51], nematic liquid crystals with competing orientational and thermal effects [52] and appear also in BECs with simultaneous local and long-range (e.g. dipolar) interactions [53, 54].

The corresponding hydrodynamic formulation of the NSE is obtained by means of the Madelung transform E = ρ 1 / 2 e iφ ,

$$\partial _ { t } \rho + \nabla \cdot ( \rho \mathbf v ) = 0 \quad ( 3 ) \quad \text{ we } \mathrm d \epsilon \\ \partial _ { t } \psi + \frac { 1 } { 2 } v ^ { 2 } = - \frac { c ^ { 2 } } { n _ { 0 } ^ { 3 } } n _ { 2 } \rho - \frac { c ^ { 2 } } { n _ { 0 } ^ { 3 } } \hat { \mathbf n } _ { n l } \rho + \frac { c ^ { 2 } } { 2 k ^ { 2 } n _ { 0 } ^ { 2 } } \frac { \nabla ^ { 2 } \rho ^ { 1 / 2 } } { \rho ^ { 1 / 2 } } \quad ( 4 ) \quad \text{ chara} \\ \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot$$

where the optical intensity ρ corresponds to the fluid density and v = c kn 0 ∇ φ ≡ ∇ ψ is the flow velocity. On the right-hand side of (4), the first term provides the local repulsive interactions related to the positive bulk pressure P = c 2 n 2 2 n 3 0 ρ 2 . The second one gives rise to a nonlocal interaction potential, while the last term, directly related to diffraction, is the analogue of the Bohm quantum potential whose gradient corresponds to the so-called quantum pressure.

We finally observe how a z -dependent response function would actually lead to a 'non-causal' photon-fluid, as the nonlocal interactions would depend on both directions of z = ( c/n 0 ) . t Such non-causality is actually fictitious since it originates from the mapping of the spatial z -direction into a time-coordinate. However, it suggests that a safe interpretation of the propagation coordinate in terms of a time variable would require a z -independent response kernel.

## A. Bogoliubov-de Gennes equations and excitation spectrum

The first-order complex fluctuations ε ( r , t ) of the optical field can be described in terms of Bogoliubov excitations on top of the photon fluid. Linearizing Eq. (1) around a background solution, E = E 0 (1 + ε + ... ) with E 0 = ρ 1 / 2 0 e iφ 0 , we obtain the nonlocal Bogoliubovde Gennes equations

$$( \partial _ { T } - i \frac { c } { 2 k n _ { 0 } } \partial _ { S } ) \varepsilon = - i \frac { \omega } { n _ { 0 } } ( n _ { 2 } + \hat { n } _ { n l } ) \rho _ { 0 } ( \varepsilon + \varepsilon ^ { * } ) \quad ( 5 ) \quad \text{ble o} ] \\ ( \partial _ { T } + i \frac { c } { 2 k n _ { 0 } } \partial _ { S } ) \varepsilon ^ { * } = i \frac { \omega } { n _ { 0 } } ( n _ { 2 } + \hat { n } _ { n l } ) \rho _ { 0 } ( \varepsilon + \varepsilon ^ { * } ) \quad ( 6 ) \quad \text{grou}$$

$$\top \, \varepsilon \overline { 2 k n _ { 0 } } \, \varrho _ { S J } \varsigma \, - \, \varepsilon \, \overline { n _ { 0 } } \, \left ( \nu \, \varepsilon \,.$$

in which we have defined the usual comoving derivative ∂ T = ∂ t + v 0 · ∇ , with v 0 = c kn 0 ∇ φ 0 , and the spatial differential operator ∂ S = 1 ρ 0 ∇· ( ρ 0 ∇ ). In the spatially homogeneous case where both the background density ρ 0 and velocity v 0 do not depend on the transverse coordinates, the plane-wave solutions of Eqs. (5)-(6) satisfy the dispersion relation

$$\Omega ^ { 2 } = c _ { s } ^ { 2 } K ^ { 2 } ( 1 + \theta \frac { \gamma } { n _ { 2 } } \tilde { R } ( K, n _ { 0 } \Omega / c ) + \frac { \xi ^ { 2 } } { 4 \pi ^ { 2 } } K ^ { 2 } ). \quad ( 7 ) \quad \underset { \text{due} } { \text{duc} }.$$

where K is the wavenumber of the mode, Ω = Ω ′ -K v · 0 its angular frequency in the locally-comoving background frame and ˜ is the three-dimensional Fourier transform R of the response function R ( r , z ). We remark that here the angular frequency Ω ′ actually corresponds to the longitudinal wavenumber K z expressed in temporal units via Ω ′ = K c/n z 0 while K =( K K x , y ) is the transverse wavevector.

In analogy to purely local BECs and photon-fluids [36], we defined in Eq. (7) the sound speed c 2 s ≡ dP ρ ( 0 ) dρ = c 2 n 2 n 3 0 ρ 0 and the healing length, ξ = λ/ 2 √ n n ρ 0 2 0 as the characteristic length separating the linear (phononic) and quadratic (single-particle) regime of the dispersion relation for γ = 0.

The length ξ determines the critical wavenumber K c = 2 π/ξ associated to the breakdown of Lorentz invariance, generally expected to occur in quantum gravity phenomenology at the Planck scale. Low energy modes with K /lessmuch K c propagate indeed at the invariant universal speed c s , while at higher momenta K /greatermuch K c the terms arising from the quantum pressure become dominant and the group velocity of excitations increases with K .

Within the paraxial approximation and in the presence of nonlocal processes with negligible longitudinal dependence, the main contribution to the non-locality comes from the K -dependence of ˜ and we can thus safely asR sume ˜ ( R K,n 0 Ω /c ) /similarequal ˜ ( R K, 0). This is indeed the case of thermo-optical nonlinearities dominated by the transverse diffusion of heat [34, 55] that we will discuss in the next sections. From now on we shall ignore the zdependence of R .

We finally observe that on the basis of Eq. (7) wave propagation for focusing nonlinearities θ = -1 is allowed only for wavenumbers K such that Ω 2 &gt; 0, i.e. ˜ ( R K ) &lt; n 2 γ (1 + K /K 2 2 c ). Negative values of Ω 2 corresponds to exponentially-growing modes characteristic of linearly-unstable flows. In the defocusing case θ = 1, the system is neutrally stable to perturbations of all wavenumbers, hence supporting travelling waves. While the plane wave solution is always modulationally stable, instabilities and wave-breaking phenomena [56] are expected in the presence of inhomogeneous beams and/or discontinuous response kernels. In spite of this fact, stable operation in nonlocal photon fluids has been experimentally demonstrated even in the presence of background inhomogeneities [30, 42].

## III. THERMO-OPTICAL NONLOCALITY AND MASSIVE EXCITATIONS

The functional form of ˜ ( R K ) depends on the specific nonlocal process under consideration. A case of particular interest is provided by light propagation in thermo-optical media, where the change of refractive index ˆ n nl ρ = n th arises from the temperature increase due to the residual laser absorption. The heat diffuses

through the material and eventually across the boundaries of the medium. As a result, the shape of the response function will strongly depend also on the transverse boundary conditions [57]. This might open interesting perspectives in experiments since it could be possible to tailor the nonlocal response of the medium, e.g. by acting on the geometry of sample, in order to modify the dispersion (7) and in turn the physical properties of the collective excitations [30].

In the limit of an infinite medium in the two transverse dimensions n th is coupled to the optical intensity through the stationary heat equation [58, 59]

$$- \nabla ^ { 2 } n _ { t h } = \frac { \alpha | \beta | } { \kappa } \rho \quad \quad \ \ ( 8 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \quad - \nabla ^ { 2 } n _ { t h } = \frac { \alpha | \beta | } { \kappa } \rho \quad \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where κ is the thermal conductivity of the material, α its linear absorption coefficient and β = | ∂n th /∂T | is the change in the refractive index with respect to the temperature. The heat equation (8) dictates that the corresponding range of the nonlocal interactions between photons is infinite (infinite-range nonlocality) [59]. The nature of the nonlinearity (focusing or defocusing, leading to attractive or repulsive interactions) depends on the sign of β . Here we take the absolute value | β | since the sign is already considered in (2) by the coefficient θ .

Fourier transforming the expression n th = γθR ( r ) ∗ ρ and Eq. (8) one can readily verify that ˜ n th = γθR K ρ ˜ ( )˜ = α β | | κK 2 ˜ and thus ρ ˜ ( R K ) ∝ 1 /K 2 . This implies that the convolution integral ˆ n nl ρ is, up to a constant, the solution of Eq. (8) with R ( r ) being the Green's function of the 2D Laplacian operator. Hence, the following relation holds: γθ ∇ 2 R ( r ) = -( α β /κ δ | | ) ( r -r ′ ).

## A. Defocusing nonlocal nonlinearity: massive phonons

Using the above ˜ ( R K ) = α β | | κγK 2 in Eq. (7) and considering a defocusing nonlinearity θ = 1 we find

$$\Omega ^ { 2 } = \Omega _ { 0 } ^ { 2 } + c _ { s } ^ { 2 } K ^ { 2 } \left ( 1 + \frac { \xi ^ { 2 } } { 4 \pi ^ { 2 } } K ^ { 2 } \right ) \quad \ \ ( 9 ) \quad \ \ \ \ \ t h e _ { \text{res} } ^ { \text{w.w.} }$$

where Ω 0 = c √ α β | | κn 3 0 ρ 0 has indeed the dimensions of a frequency. Hence, we can thus identify /planckover2pi1 Ω 0 with the rest energy of a particle and write /planckover2pi1 Ω 0 = mc 2 s , where m is the rest mass and c s plays the role of the light speed. Defining the excitation momentum p = /planckover2pi1 K and the critical momentum p c = /planckover2pi1 K c = h/ξ , Eq. (9) can be rewritten in the form

$$\mathcal { E } ^ { 2 } = m ^ { 2 } c _ { s } ^ { 4 } + c _ { s } ^ { 2 } p ^ { 2 } \left ( 1 + \frac { p ^ { 2 } } { p _ { c } ^ { 2 } } \right ) \. \quad \quad ( 1 0 ) \quad \text{ in } v$$

The properties of the above dispersion strongly depend on the thermo-optical coefficients of the material used to produce the photon fluid. Here we are interested in investigating the regime in which p c /greatermuch mc s . In this case,

Eq. (10) is a generalization of the Bogoliubov dispersion relation describing massive collective excitations with high-energy, Lorentz-violating corrections. Similar modified dispersion laws with extra momentum-dependent terms appear in several phenomenological approaches to quantum gravity, where p c is typically associated to the Planck momentum [60].

The dispersion curve (10) interpolates between three different regimes depending on the fluctuations momentum.

When p /greatermuch p c , the quartic term dominates and Eq. (10) approximates the free-particle behaviour E ≈ c p /p s 2 c . Using the above definitions of p c and c s and the photon momentum p γ = /planckover2pi1 n k 0 we get E ≈ cp / 2 (2 p γ ) (or equivalently E ≈ p / 2 2 m γ introducing an effective photon mass m γ = p γ /c ). Therefore in analogy to BEC analogue models, the excitation energy tends to the energy of the individual particles forming the background fluid, i.e. in our case the photons.

In the intermediate regime, mc s /lessorsimilar p /lessmuch p c , we obtain the 'relativistic' dispersion relation for a massive particle E ≈ √ p c 2 2 s + m c 2 4 s , with the speed of sound playing the role of the speed of light. These are collective excitations exactly like usual phonons in local quantum fluids, but possessing a finite rest mass.

At lower momenta, p /lessmuch mc s , the phonon modes enter the non-relativistic regime: the energy-momentum relation reduces to E ≈ p / 2 2 m + mc 2 s , where the first term is the kinetic energy of a particle of mass m and the second its constant rest mass energy.

We notice that the rest frequency Ω 0 depends only on the strength of the thermo-optical (nonlocal) nonlinearity. The latter is thus responsible for the generation of the gap in the dispersion relation and hence for the onset of the excitation mass. On the other hand, the 'invariant' limit speed c s is determined solely on the local defocusing effect.

The rest frequency Ω 0 can also be expressed in terms of the sound speed as Ω 0 = c s √ α β | | κn 2 = c /λ s C . The characterisitic length-scale λ C , given by the square root of the ratio between the local and nonlocal coefficients, corresponds to the acoustic analogue of the reduced Compton wavelength of the particle λ C = /planckover2pi1 / mc ( s ). The inverse of this length defines the above non-relativistic limit through K /lessmuch λ -1 C and, as we shall see in Sect. V, it provides also some of the fundamental characteristic scales of the emergent gravitational force.

We conclude the section briefly discussing a more realistic model of the thermo-optical nonlinearity [30, 56, 61] given by

$$- \nabla ^ { 2 } n _ { t h } + \frac { n _ { t h } } { \sigma ^ { 2 } } = \frac { \alpha | \beta | } { \kappa } \rho \, \quad \quad ( 1 1 )$$

in which the effects of the distant boundaries has been included in the distributed loss term n th /σ 2 , where σ is the length-scale of the nonlocal interaction.

Eq. (11) allows us to continuosly describe the transition from an infinite-range to finite-range thermo-optical

nonlocality and has provided a theoretical framework for the phenomenological Lorentzian response adopted in previous experiments (see e.g. Refs. [34, 40]). The Fourier-transformed response associated to (11) γR ˜ = α β | | κ σ 2 1+ σ K 2 2 has indeed a Lorentzian shape, where 2 /σ is its full-width at half-maximum. The dispersion (10) is thus modified as

$$\mathcal { E } ^ { 2 } = m ^ { 2 } c _ { s } ^ { 4 } \frac { p ^ { 2 } } { p _ { n l } ^ { 2 } + p ^ { 2 } } + c _ { s } ^ { 2 } p ^ { 2 } \left ( 1 + \frac { p ^ { 2 } } { p _ { c } ^ { 2 } } \right ) \quad ( 1 2 ) \quad \text{grow} \\ \quad \text{of tl}$$

where we introduced the nonlocal momentum p nl = /planckover2pi1 /σ . The above response kernel reduces to the ideal form of the infinite space model γR K ˜ ( ) = ( α β | | κ ) /K 2 in the limit of σK /greatermuch 1. Such a regime can be reasonably reproduced by means of suitable background optical beams comprising wavevectors only of K /greatermuch 1 /σ . This procedure has been implemented in a lead-doped glass experiment [55]. Before being launched into the nonlinear medium, the laser beam has been passed through a phase mask generating a ring-shaped beam with zero intensity at K = 0 and large-enough transverse wavevectors. Using this technique, the authors demostrated a nonlocal thermo-optical nonlinearity with σK ≈ 20. In this case p /greatermuch p nl , and Eq. (12) well approximates the massive Bogoliubov dispersion (10). However, for finite p nl the two relations will eventually differ at arbitrarily low momenta, as the gap in (12) arises only in the singular limit of an infinite-range nonlocality, p nl = 0 ( σ →∞ ).

## B. Focusing nonlocal nonlinearity: Jeans instability

For θ = -1 the hydrodynamic equations Eq. (3)-(4) together with Eq. (11) describe a (2 + 1)-dimensional quantum fluid with local repulsive and finite-range attractive interactions. In the ideal case of an infinite medium the model reproduces the nonlinear evolution of a self-gravitating BEC [62], where the nonlocal change of refractive index n th , solution of Eq. (8), mimics a Newtonian potential generated by the fluid mass density. In the absence of local interactions, Eq. (1)-(8) are indeed formally equivalent to the Schr¨dinger-Newton (SN) o equation in two spatial dimensions [55, 63], originally proposed by Diosi [64] and Penrose [65] as a model for quantum wave function collapse (see also [66] for further discussion).

Concerning the dynamics of elementary excitations, in the more general case the dispersion relation reads

$$\mathcal { E } ^ { 2 } = - m ^ { 2 } c _ { s } ^ { 4 } \frac { p ^ { 2 } } { p _ { n l } ^ { 2 } + p ^ { 2 } } + c _ { s } ^ { 2 } p ^ { 2 } \left ( 1 + \frac { p ^ { 2 } } { p _ { c } ^ { 2 } } \right ) \. \quad ( 1 3 ) \quad \text{with}$$

The negative sign in front of the rest-energy term originated from to the attractive photon interactions gives rise to two fundamentally different behaviours at high and low momenta. The critical wavenumber K = K J separating these two regimes is implicitly defined by the condition at which the wave frequency (energy) vanishes,

$$\frac { \sigma ^ { 2 } } { \lambda _ { C } ^ { 2 } } \frac { 1 } { 1 + \sigma ^ { 2 } K ^ { 2 } } - \left ( 1 + \frac { K ^ { 2 } } { K _ { c } ^ { 2 } } \right ) = 0 \. \quad ( 1 4 )$$

For high wavenumbers K &gt; K J , the local repulsive interactions and the quantum pressure are sufficiently strong to countebalance the nonlocal attractive forces and the waves are freely oscillating. In the opposite case, we have growing excitation modes revealing the linear instability of the system (see also Refs. [67, 68] for a discussion in self-gravitating BECs).

In the hydrodynamic (Thomas-Fermi) approximation ( K /lessmuch K c ) and for ( σK ) /greatermuch 1, Eq. (14) simply reduces to the ordinary Jeans instability condition, where K J corresponds to the Compton wavenumber of the particle K J = λ -1 C . In astrophysics such instability is thought to be responsible for the collapse of interstellar gas clouds eventually leading to star formation.

In the stable regime, K &gt; λ -1 C , and for K /lessmuch K c Eq. (13) yields the tachyonic dispersion relation E ≈ √ p c 2 2 s -m c 2 4 s , with real energy and momentum and imaginary rest mass. Excitations of any wavenumuber in fact propagate at supersonic group velocities, with the invariant c s being now a lower limit for propagation speeds. From now on we will focus on the fully stable case of a photon-fluid with defocusing non-local nonlinearity.

## IV. INHOMOGENOUS FLOWS AND MASSIVE KLEIN-GORDON EQUATION

For purely local interactions γ = 0, a formal equivalence can be established between phonons propagating on top of the photon fluid and the evolution of scalar fields in curved spacetime [36]. The equation of motion is typically derived by linearizing Eqs. (3)-(4) around a background state since phonons, i.e. the acoustic elementary excitations, are defined as the first-order fluctuations of the quantities describing the mean fluid flow: ρ = ρ 0 + /epsilon1ρ 1 + O /epsilon1 ( 2 ), ψ = ψ 0 + /epsilon1ψ 1 + O /epsilon1 ( 2 ). When the terms arising from quantum pressure are negligible, the phonon dynamics is fully described by a single secondorder equation for the linearized velocity-potential which has the form of the Klein-Gordon equation for a massless scalar field

$$\Box \psi _ { 1 } & \equiv \frac { 1 } { \sqrt { - g } } \partial _ { \mu } \left ( \sqrt { - g } \, g ^ { \mu \nu } \, \partial _ { \nu } \, \psi _ { 1 } \right ) \quad \quad ( 1 5 ) \\ \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \ dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots\dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \delta \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \Dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots$$

propagating in a (2+1)-dimensional curved spacetime whose geometry is described by the acoustic metric g µν , with inverse g µν and determinant g ,

$$g _ { \mu \nu } & = \left ( \frac { \rho _ { 0 } } { c _ { s } } \right ) ^ { 2 } \left ( \begin{array} { c c } - ( c _ { s } ^ { 2 } - v _ { 0 } ^ { 2 } ) & - \mathbf v _ { 0 } ^ { T } \\ - \mathbf v _ { 0 } & \mathbf I \end{array} \right ) \quad \text{(1)}$$

where I is the two-dimensional identity matrix.

The above scenario is deeply modified in the presence of both local and nonlocal nonlinearities. In this context,

it is convenient to derive the acoustic metric directly from the nonlocal Bogoliubov-de Gennes equations (5)-(6).

To this end, we apply the operator ( ∂ T + i c 2 kn 0 ∂ S )( 1 ρ 0 ) to Eq. (5) and we obtain

$$( \partial _ { T } + i \frac { c } { 2 k n _ { 0 } } \partial _ { S } ) \frac { 1 } { \rho _ { 0 } } ( \partial _ { T } - i \frac { c } { 2 k n _ { 0 } } \partial _ { S } ) \varepsilon = \frac { c _ { s } ^ { 2 } } { \rho _ { 0 } } \partial _ { S } \varepsilon \\ - i \gamma ( \frac { \omega } { n _ { 0 } } \partial _ { T } + i \frac { c ^ { 2 } } { 2 n _ { 0 } ^ { 3 } } \partial _ { S } ) \frac { 1 } { \rho _ { 0 } } R * [ \rho _ { 0 } ( \varepsilon + \varepsilon ^ { * } ) ] \quad ( 1 7 ) \\ \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \end{array}$$

As in the local case, we remind that the gravitational analogy holds in the phononic regime in which the dispersion relation takes the relativistic form with a limit propagation speed, i.e. for wavenumbers K /lessmuch K c . The corresponding equation for the excitation field can thus be obtained by ignoring the higher-order spatial derivatives in Eq. (17), which indeed are responsible for the Lorentzbreaking, K 4 -terms in the dispersion relation. A close inspection of Eq. (17) suggests that such approximation corresponds to take the diffractionless limit k →∞ [44] or, equivalently, neglect the quantum pressure terms arising from the linearized hydrodynamic equations [36]. In this limit Eqs. (3)-(4) indeed reduce to the Navier-Stokes equations for a barotropic, irrotational and inviscid fluid, in which the Lorentz simmetry associated to phonon dynamics is not explicitly broken [5].

Under this approximation and using the fact that the background density ρ 0 satisfies the continuity equation (3) with v = v 0 , Eq. (17) can be re-written as

$$\Box \varepsilon \, = \, - i \gamma \frac { \omega } { n _ { 0 } } ( \partial _ { T } + \nabla \cdot \mathbf v _ { 0 } ) R * [ \rho _ { 0 } ( \varepsilon + \varepsilon ^ { * } ) ] \quad \begin{pmatrix} \ln \\ \text{citati} \\ a \text{ rel} \\ \text{ence} \\ \text{Kleir} \\ \text{densi} \end{pmatrix}$$

where

$$\Box \equiv ( \partial _ { T } + \nabla \cdot \mathbf v _ { 0 } ) \partial _ { T } - \nabla \cdot ( c _ { s } ^ { 2 } \nabla \ ) \quad \ ( 1 9 ) \quad \text{from} \\ \text{now}$$

is precisely the d'Alambertian operator associated with the acoustic metric g µν .

In the purely local case γ = 0, we recover the usual Klein-Gordon equation for a massless particle on curved spacetime, here described by the complex field ε . For a spatially homogeneous background, with constant density ρ 0 and constant flow velocity v 0 = c kn 0 ∇ φ 0 , the operator ∂ T = ∂ t + v 0 · ∇ commutes with the convolution operation and using Eqs. (5)-(6) we find that the wave equation (18) becomes independent of ε ∗ ,

$$( \partial _ { T T } ^ { 2 } - c _ { s } ^ { 2 } \nabla ^ { 2 } ) \varepsilon = c _ { s } ^ { 2 } \frac { \gamma } { n _ { 2 } } R ( r ) * \nabla ^ { 2 } \varepsilon \,. \quad \ ( 2 0 ) \quad \stackrel { \lambda } { \text{scril} } \\ \text{som} \\ \text{char}$$

It is immediate to verify that the Fourier transform of Eq. (20) leads to the dispersion law (7).

The complex fluctuations ε can be easily linked to the real density and phase perturbations though the relations ρ 1 = ρ 0 ( ε + ε ∗ ) and φ 1 = ( i/ 2)( ε ∗ -ε ). By means of these expressions and using the relation between the optical phase and velocity-potential of the flow, ψ 1 = ( c/kn 0 ) φ 1 , one can split Eq. (18) into the following system of waveequations

$$\Box \psi _ { 1 } = - \gamma \frac { c ^ { 2 } } { n _ { 0 } ^ { 3 } } ( \partial _ { T } + \nabla \cdot \mathbf v _ { 0 } ) R * \rho _ { 1 } \quad ( 2 1 )$$

$$\Box \left ( \frac { \rho _ { 1 } } { \rho _ { 0 } } \right ) = \gamma \frac { c ^ { 2 } } { n _ { 0 } ^ { 3 } } \nabla \cdot ( \nabla - \nabla \ln \rho _ { 0 } ) R \ast \rho _ { 1 } \quad \ ( 2 2 )$$

For local fluids γ = 0 Eq. (21) reduces to the massless Klein-Gordon equation (15) for the velocity-potential perturbations ψ 1 and an equation of the same form is satisfied also by the relative density fluctuations ρ /ρ 1 0 .

In the ideal case of infinite-range thermo-optical nonlinearity, σ → ∞ , the response function satisfies γ ∇ 2 R ( r ) = -( α β /κ δ | | ) ( r -r ′ ). Using this result and considering a nearly homogeneous background density [69], Eq. (22) takes the form of the massive Klein-Gordon equation in curved spacetime

$$\Box \rho _ { 1 } + \Omega _ { 0 } ^ { 2 } \rho _ { 1 } = 0 \,. \quad \quad \quad ( 2 3 )$$

In the more realistic case of finite-range thermo-optical nonlinearities, Eq. (23) remains basically valid for perturbations with wavenumbers 1 /σ /lessmuch K /lessmuch K c .

## V. EMERGENT GRAVITATIONAL DYNAMICS

In the previous sections we have seen that phonon excitations in our system behave as massive particles with a relativistic energy-momentum relation. In the presence of inhomogeneous flows we also derived a massive Klein-Gordon equation on the acoustic metric for the density fluctuations that thus reproduce the evolution of massive scalar fields on a curved spacetime. The spacetime curvature which mimics the gravitational field arises from the inhomogeneity of the background, whose dynamics however is governed by a non-relativistic nonlinear equations (cf Eq. (1) or, equivalently, Eqs. (3)-(4)). As such the analogy works only at the kinematical level: the fluctuations propagate in a given background solution associated to a specific spacetime configuration. All effects due to gravitational back-reaction are neglected, i.e. the spacetime geometry is not modified by the perturbations propagating on it. While under certain conditions it is possible to extend the analogy and include in a geometric framework even the evolution of the background [70-72], there is no possibility in general to describe the dynamics of the acoustic metric in terms of something similar to Einstein's equations. The situation changes if one considers relativistic models and interesting progresses in this direction have been made, e.g. in the framework of relativistic BECs [73].

Nevertheless, as mentioned before a kind of gravitational dynamics may emerge even in non-relativistic BECs upon suitable modifications of the standard equations to break the U (1) symmetry associated with the

conservation of particle number [47]. In such a modified model, the massive excitations feel a Newtonian gravitational potential whose source is related to the excitation density.

In the following we show that a similar scenario arises also in our nonlocal photon-fluid.

## A. The Newtonian limit

In Sect. IV we treated the general case of an inhomogeneous background, i.e. of an arbitrary curved spacetime simulating a generic gravitational field. Since the background dynamics is non-relativistic we expect to find at most a kind of Newtonian gravity, as previously shown in other non-relativistic frameworks [47]. We thus focus on a nearly-homogeneous background corresponding to a weak gravitational field because, in analogy to the weak-field approaximation of General Relativity (GR), is in this limit that a Newtonian-like gravity is expected to emerge.

In GR the weakness of gravitational field allows for the decomposition of the metric into a flat Minkowski spacetime, η µν , plus a small perturbation g µν = η µν + h µν , | h µν | /lessmuch 1 where h µν represents the weak deviations from flatness g µν i.e. the gravitational field. In this regime the Newtonian potential is related to the metric though the equation g 00 = η 00 + h 00 /similarequal -(1 + 2Φ N /c 2 ), which follows from a non-relativistic limit of the geodesics equation [74].

In analogy to the above, we consider a photon fluid with zero-flow and a spatially-localized inhomogeneity in the density, i.e. E 0 = ρ 1 / 2 ∞ (1 + u r ( )) with u /lessmuch 1 and u → 0 at infinity. We thus assume that only a small region of the fluid deviates from the constant asymptotic value of the density ρ ∞ . This implies a rescaling of the speed of sound c 2 s = c 2 ∞ (1 + 2 u r ( )) and thus of the 00component of the acoustic metric (16). On the basis of the hydrodynamic equations (3)-(4), a density inhomogeneity would also imply an inhomogeneous flow. However as demonstrated in Ref. [47], the velocity perturbations do not contribute at first order to the anologue gravitational potential. In other words, at the first order all the information about the gravitational potential is encoded in the density perturbation. This result is general in acoustic models and does not depend on the specific fluid under consideration. Therefore for the sake of simplicity and without loss of generality, here we assume deviations in the density only.

## B. Non-relativistic phonon dynamics

In order to show the emergence of a gravitational potential term we should derive the equation of motion for excitations in the non-relativistic regime of the Bogoliubov spectrum, i.e. when p /lessmuch mc ∞ , and for an homogeneous background except for the small densityinhomogeneity u r ( ).

To this end, we start directly from the nonlocal waveequation (18) for the complex excitation field ε . Setting ρ 0 = ρ ∞ (1 + 2 u r ( )) and v 0 = 0 we get

$$\Box \varepsilon = c _ { \infty } ^ { 2 } \frac { \gamma } { n _ { 2 } } \nabla ^ { 2 } R ( r ) * \left [ ( 1 + 2 u ( r ) ) \varepsilon \right ] \quad ( 2 4 )$$

In deriving (24) we disregarded terms containing the spatial derivatives of u r ( ) and the products u r ( ) ∇ 2 ε . The former are negligible in the asymptotic region and the latter are suppressed both by the smallness of u and by the fact that we are interested into the non-relativistic regime p /lessmuch mc ∞ . Restricting ourselves to the case of an infinite-range thermo-optical nonlinearity, Eq. (24) further simplifies and reads

$$( \partial _ { t t } ^ { 2 } - c _ { \infty } ^ { 2 } \nabla ^ { 2 } ) \varepsilon + \Omega _ { \infty } ^ { 2 } ( 1 + 2 u ( r ) ) \varepsilon = 0 \quad ( 2 5 )$$

where we have defined the asymptotic rest frequency Ω ∞ = c ∞ √ α β /κn | | 2 .

The non-relativistic limit p /lessmuch mc ∞ (or equivalently, c ∞ → ∞ ) means that the kinetic energy of the particle should be small with respect to its mass energy /planckover2pi1 Ω ∞ = mc 2 ∞ . Making the ansatz ε = ϕexp ( -i Ω ∞ t ) to factor out the rest frequency (i.e. the contribution to the total energy due to the rest energy of the particle) we can approximate [75]

$$\partial _ { t t } ^ { 2 } \varepsilon \simeq ( - 2 i \Omega _ { \infty } \partial _ { t } \varphi - \Omega _ { \infty } ^ { 2 } \varphi ) e ^ { - i \Omega _ { \infty } t } \.$$

Substituting the above expression into Eq. (25) we get the Schr¨dinger equation for a particle of mass o m

$$i \hbar { \partial } _ { t } \varphi = - \frac { \hbar { ^ { 2 } } } { 2 m } \nabla ^ { 2 } \varphi + m c _ { \infty } ^ { 2 } u ( r ) \varphi \text{ \quad \ \ } ( 2 6 )$$

subject to an external potential proportional to u r ( ). The latter can be formally identified as a gravitational potential defining Φ G = c 2 ∞ u r ( ). We finally remark that for finite-range thermo-optical nonlinearities, Eq. (26) would remain approximately valid in the momentum range p nl /lessmuch p /lessmuch mc ∞ .

## C. Modified Poisson equation

In the previous section we have found that the low-energy evolution of massive phonons obeys the Schr¨dinger equation for a nonrelativistic quantum partio cle in an external potential Φ G . Our interpretation of Φ G as a gravitational potential is based on the way it enters in the Schr¨dinger equation and because it is related to o the 00-component of the acoustic metric, similarly to the Newtonian potential in the weak-field approximation of GR. In this framework we should find that in the appropriate limits Φ G also obeys a kind of Poisson's equation [76].

Following the same argument of Ref. [47], since the Newtonian potential is the manifestation of small deviations of the order parameter E from perfect homogeneity, the corresponding Laplace's equation should be encoded in the nonlinear evolution equation (1). On the other hand, we expect the source term to be directly related to the phononic fluctuations. Indeed in Newtonian gravity the only source for the gravitational field is a mass-density distribution and in our system this can originate only from the massive elementary excitations. As a result, the non-relativistic massive phonons should experience a kind of gravitational potential generated by themselves, i.e. they should feel their own gravity.

This self-interaction can thus be derived from the nonlinear equation Eq. (1), but adding the corrections to the mean-field dynamics induced by the fluctuations, in order to see how the phonons backreact over the background fluid.

Backreaction effects can be calculated expanding Eq. (1) up to second order, E = E 0 + η 1 + η 2 = E B + η 1 : here η 1 and η 2 are linear and quadratic quantities in the fluctuation amplitude, respectively, and we introduced new variable E B including the modifications to the zeroorder dynamics, that is to say the backreaction [77].

Substituting the above ansatz in Eq. (1) with ∆ n = n 2 | E | 2 + γR ∗ | E | 2 we obtain

$$\partial _ { t } ( E _ { B } + \eta _ { 1 } ) = \frac { i c } { 2 k n _ { 0 } } \nabla ^ { 2 } ( E _ { B } + \eta _ { 1 } ) - \frac { i \omega } { n _ { 0 } } ( E _ { B } + \eta _ { 1 } ) \Delta n _ { B } \begin{array} { c } K _ { c } = 2 \\ \lambda / 2 \sqrt { n _ { 0 } } \\ ( 2 7 ) & \text{Poisson} \\ \dots \end{array}$$

where ∆ n B = ( n 2 + γR ∗ )( | E B | 2 +2Re( E η ∗ B 1 )+ | η 1 | 2 ) and we used the effective time coordinate t = ( n /c z 0 ) .

Since E B consists only of zero-order and second-order terms in the fluctuation amplitude, all linear quantities in η 1 must vanish. The zeroing of the linear fluctuation terms in (27) leads to the Bogoliubov-de Gennes equations (5)-(6) upon substituting η 1 = E ε 0 and E B = E 0 . What remains is a nonlinear evolution equation for E B in which the fluctuations appear quadratically

$$\partial _ { t } E _ { B } \ & = \ \frac { i c } { 2 k n _ { 0 } } \nabla ^ { 2 } E _ { B } - \frac { i \omega } { n _ { 0 } } E _ { B } [ ( n _ { 2 } + \gamma R * ) | E _ { B } | ^ { 2 } + 2 n _ { 2 } | \eta _ { 1 } | ^ { 2 } ] \quad \underset { n _ { 0 } } { \text{ use } } \quad \text{$\theta_{c}$} \\ & - \ \frac { i \omega } { n _ { 0 } } [ n _ { 2 } E _ { B } ^ { * } \eta _ { 1 } ^ { 2 } + \gamma ( E _ { B } R * | \eta _ { 1 } | ^ { 2 } + 2 \eta _ { 1 } R * \text{Re} ( E _ { B } \eta _ { 1 } ^ { * } ) ) ] \ \text{sidered w} \\ & ( 2 8 ) \ \searrow \ \text{$K$} _ { c },$$

(28)

Here, η 1 is in general time-dependent and for the purpose of simulating the effects of quantum and/or thermal fluctuations on the mean-field dynamics in analogy to real quantum gases, it could be taken as a stochastic variable. However, for our purposes we are now interested in calculating the backreaction effects on stationary background solutions, i.e. stationary spacetime geometries. The corresponding mean-field equation is then obtained by time-averaging Eq. (28) and replacing the fluctuating quantities with their mean values. For γ = 0, the averaged (28) closely resembles the modified GrossPitaevskii equation with beyond-mean-field corrections due to fluctuations [78-80]. In this context the quadratic terms n ( r ) = 〈| η 1 | 2 〉 , and m ( r ) = 〈 η 2 1 〉 play the role of the density of non-condensed particles (i.e. the 'outof-condensate' photons) and of the anomalous density in the Bogoliubov-Popov-Beliaev approximation. The quantities γ η R 〈 1 ∗ Re( E η ∗ B 1 ) 〉 and γR ∗ 〈| η 1 | 2 〉 provide further corrections due to nonlocality.

Setting 〈 E B 〉 = ρ 1 / 2 ∞ (1 + u r ( )) with u /lessmuch 1 in the averaged Eq. (28) yields

$$\frac { c } { 2 k n _ { 0 } } \nabla ^ { 2 } u - \frac { \omega } { n _ { 0 } } \rho _ { \infty } ( n _ { 2 } + \gamma R * ) ( 1 + 2 u ) = \\ \frac { \omega } { n _ { 0 } } [ n _ { 2 } ( 2 \mathfrak { n } ( \mathbf r ) + \mathfrak { m } ( \mathbf r ) ) + \gamma ( R * \mathfrak { n } ( \mathbf r ) + 2 \mathfrak { g } ( \mathbf r ) ) ] \quad ( 2 9 )$$

where g ( r ) = 〈 η 1 [ R ∗ Re( η ∗ 1 )] 〉 .

In the ideal case of our interest of infinite-range thermo-optical nonlinearity, and applying the operator c/ (2 kn 0 ) ∇ 2 to both sides of Eq. (29), we get

$$& \overset { \nu } { 1 } \colon \\ & \underset { \i e } { \i d } \\ & \circ - \\ & \overset { \circ } { - } \\ & = \begin{array} { c } \overset { 1 } { 1 } \colon \\ \i e \\ \i d \end{array} \\ \frac { 1 } { K _ { c } ^ { 2 } } \nabla ^ { 4 } \Phi _ { G } - \nabla ^ { 2 } \Phi _ { G } + \frac { 1 } { \lambda _ { C } ^ { 2 } } \Phi _ { G } + \frac { c _ { \infty } ^ { 2 } } { 2 \lambda _ { C } ^ { 2 } } = \\ & \overset { c _ { \infty } ^ { 2 } } { 2 \lambda _ { C } ^ { 2 } \rho _ { \infty } } [ \lambda _ { C } ^ { 2 } \nabla ^ { 2 } ( 2 \mathfrak { n } ( \mathbf r ) ) - \mathfrak { n } ( \mathbf r ) + 2 \nabla ^ { 2 } \mathfrak { g } ( \mathbf r ) ] \quad ( 3 0 ) \\ & \quad \cdot \\ & \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \end{array} \quad.$$

where here in the definition of the critical wavevector K c = 2 π/ξ ∞ we used the asymptotic healing length ξ ∞ = λ/ 2 √ n n ρ 0 2 ∞ .

Eq. (30) can be interpreted as a modified fourth-order Poisson equation for the potential Φ G , provided that we are able to identify the right-hand side with a genuine source term. To this end, we remark that n , m and g are quadratic terms in η 1 , i.e. they have the dimension of a photon-fluid mass-density and are related to the phononic excitations through the relation η 1 = E ε 0 . Therefore, they can be safely interpreted as a source of the gravitational field.

In contrast to Newtonian gravity, Eq. (30) contains the additional terms ∇ 4 Φ G /K 2 c and Φ G /λ 2 C . Owing to the K 2 c coefficient, the first term suggests modifications of the Poisson equation which would become increasingly important as the analogue of the Planck scale ξ ∞ is approached. Since in the nonrelativistic limit here considered we are dealing with long-wavelength modes with K /lessmuch K c , the variations of Φ G over spatial scales of order ξ ∞ can be neglected.

The second term Φ G /λ 2 C denotes a finite range for the gravitational interaction, with a characteristic length scale given by λ C . Such a finite interaction scale for gravity would translate into a massive graviton with a mass that, in our model, corresponds to that of the massive phonons. A further comparison with the Newtonian limit of Einstein equations allows us to identify the quantity 1 / 2 λ 2 C , which indeed has the dimension of the square of an inverse length, with a cosmological constant Λ. Finally, in the right hand side of the Eq. (30) it is natural to define the analogous of the universal gravitational con-

$$\cdots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots\dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \Delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta
\text{start $G=\frac{c_{\infty}^{2} } { 2 \tilde { \lambda } _ { C } ^ { 2 } \rho _ { \infty } }.$$

In light of the above considerations, Eq. (30) takes the more meaningful form

$$\nabla ^ { 2 } \Phi _ { G } - \frac { 1 } { \lambda _ { C } ^ { 2 } } \Phi _ { G } = c _ { \infty } ^ { 2 } \Lambda + G _ { Q m a t t e r } ( { \mathbf r } ) \quad ( 3 1 ) \quad \text{such} \quad \text{dov}$$

where we introduced the mass-density distribution

$$\varrho _ { m a t t e r } ( \mathbf r ) = - \lambda _ { C } ^ { 2 } \nabla ^ { 2 } [ 2 \mathbf n ( \mathbf r ) + \mathbf m ( \mathbf r ) ] + \mathbf n ( \mathbf r ) - 2 \nabla ^ { 2 } \mathbf g ( \mathbf r ) \ ( 3 2 ) \quad \text{Her} \\ \text{...} \ \. \. \. \. \. \. \. \. \. \. \ \text{with $1$}$$

The mass-density (32) is a complicated function of quadratic fluctuation-terms, which deserves further analysis. We first note that in the non-relativistic limit, λ 2 C K 2 /lessmuch 1, the first laplacian-term is less relevant with respect to the others and can thus be neglected. Moreover, applying the product rule to ∇ 2 g ( r ) we find

$$\nabla ^ { 2 } \mathfrak { g } ( \mathbf r ) \ & = \ \langle [ R * \mathrm R e ( \eta _ { 1 } ^ { * } ) ] \nabla ^ { 2 } \eta _ { 1 } \rangle + 2 \langle \nabla \eta _ { 1 } \cdot \nabla [ R * \mathrm R e ( \eta _ { 1 } ^ { * } ) ] \rangle \quad \text{$\ u y m a u $ } \quad \text{$\ u t i n $ } \quad \text{$\ t i n $ } \quad \text{$\ t i n $ } \quad \text{$\ t i n $ } \quad \text{$\ t i n $ } \quad \text{$\ t i n $ } \quad \text{$\ t i n $ } \quad \text{$\ t i n $ } \quad \text{$\ t i n $ } \quad \text{$\ t i n $ } \quad \text{$\ t i n $ } \quad \text{$\ t i n $ } \quad \text{$\ t i n $ } \quad \text{$\ t i n $ } \quad \text{$\ t i n $ } \quad \text{$\ t i n $ } \quad \text{$\ t i n $ } \quad \text{$\ t i n $ } \quad$$

where we used the previously defined quantities n ( r ) = 〈| η 1 | 2 〉 , and m ( r ) = 〈 η 2 1 〉 and the relation ∇ 2 R ( r ) = -δ ( r -r ′ ), valid in the limit of infinite-range nonlocality. While the first two terms in (33)are difficult to handle in general, we observe that their contribution becomes less important at very low wavevectors, in which case the mass-density distribution approximates

$$\varrho _ { m a t t e r } ( \mathbf r ) \approx 2 \mathbf n ( \mathbf r ) + \mathbf m ( \mathbf r ) \quad \ \ ( 3 4 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

The above expression coincides with the density distribution obtained in Ref. [47] (see Eq. 43) where a Newtonian-like gravity has been shown to emerge in a BEC model modified with a U (1) symmetry-breaking term. The distribution (34) has an immediate physical interpretation: in analogy to BECs, the two terms n ( r ) and m ( r ) indeed correspond to the so-called normal and anomalous density encoding the effects of 'noncondensed particles' on the mean-field dynamics [78-80].

We finally remark that, similarly to the rest energy, also the Newton constant can be expressed solely in terms of the nonlocal coefficients, G = ( c 2 / 2 n 3 0 )( α β /κ | | ). The nonlocal nonlinearity is thus responsible for the emergent gravitational interaction.

To conclude, the Eqs. (26) and (31), mutually coupled via the relations η 1 = E ε 0 and (34), actually gives rise to an effective Schr¨dinger-Newton dynamics describing o the evolution of a quantum mass density experiencing its own gravitational field [66], although here the source is a more complicated function of the mass-density distribution and the gravitational interaction is characterized by a short interaction-range and a nonzero cosmological constant.

## VI. CONCLUSIONS AND FUTURE PERSPECTIVES

Quantum fluids of light such as exciton-polaritons BECs and more recently photon fluids have offered al- ternative platforms for fundamental studies of quantum many-body physics. Recent experiments in these systems provided evidence of collective many-photon phenomena, such as the emergence of a phonon-regime in the Bogoliubov dispersion [34, 35], superfluidity and nucleation of quantized vortices in the flow past a physical obstacle [30, 31] and classical wave condensation [32].

Here, we have theoretically investigated a photon-fluid with both local and nonlocal interactions from the analogue gravity perspective. We have found that collective excitations in this system display a gapped Bogoliubov spectrum which at low energies correspond to that of massive phonons with a relativistic energy-momentum relation. In the presence of an inhomogeneous flow the dynamics of the density fluctuations is equivalent to that of a massive scalar field propagating in a curved spacetime whose geometry is specified by the acoustic metric. This generalizes previous studies in local fluids to the case of massive phonons and provides a quite natural setting for analogue simulations of quantum gravity phenomenology.

The massive nature of the elementary excitations allows us to study their nonrelativistic dynamics in nearlyhomogeneous background that, as explained, corresponds to the case of a weak gravitational field. In this limits we find that the phonon-modes behave as a self-gravitating quantum system. The evolution equations are indeed the Schr¨dinger equation for a massive quantum particle, ino cluding a term that represents the interaction of the particle with its own gravitational field, and a kind of Poisson equation with a source depending on the phononic mass-density distribution. In analogy to the Newtonian limit of GR, the potential in the Poisson equation is related to the background geometry (namely to the 00component of the metric) experienced by the particles propagating on it. Since most of analogue models are dealing with massless excitations that in the framework of Newtonian gravity cannot act as sources of a gravitational field, our system is one of the very few in which a form of semiclassical gravitational dynamics can emerge.

One of the next stages of this investigation will focus on the design of realistic experimental schemes for the implementation of such photon-fluids. Apart from the analogue-gravity side, we expect these experiments to be interesting also from the perspective of the quantum fluids of light, as the interplay between local and nonlocal nonlinearities with different -possibly tunablekernels could unveil new collective many-photon phenomena and hydrodynamic phase transitions.

## Acknowledgements

I wish to thank Daniele Faccio and Antonello Ortolan for stimulating discussions and useful comments on the manuscript.

- [1] C. Barcel`, S. Liberati, M. Visser, Living Rev. Rel. o 14 , 3 (2011).
- [2] C. Barcel`, Nat. Phys. o 15 , 210, (2018).
- [3] D. Faccio, F. Belgiorno, S. Cacciatori, V. Gorini, S. Liberati and U. Moschella, Analogue Gravity Phenomenology: Horizons and from Theory to Experiment , Springer, New York, (2013).
- [4] W. G. Unruh, Phys. Rev. Lett. 46 , 1351 (1981)
- [5] M. Visser, Phys. Rev. Lett. 80 , 3436 (1998); Class. Quantum Grav. 15 , 1767 (1998).
- [6] L. J. Garay, J. R. Anglin, J. I. Cirac, and P. Zoller, Phys. Rev. Lett. 85 , 4643 (2000); Phys. Rev. A 63 , 023611 (2001).
- [7] O. Lahav, A. Itah, A. Blumkin, C. Gordon, S. Rinott, A. Zayats and J. Steinhauer, Phys. Rev. Lett. 105 , 240401 (2010).
- [8] R. Sch¨tzhold u and W. G. Unruh, Phys. Rev. D, 66 044019 (2002).
- [9] T. A. Jacobson and G. E. Volovik, Phys. Rev. D 58 , 064021 (1998).
- [10] S. Giovanazzi, Phys. Rev. Lett. 94 , 061302 (2005).
- [11] B. Reznik, Phys. Rev. D, 62 , 044044 (2000).
- [12] R. Sch¨tzhold, G. Plunien, and G. Soff, Phys. Rev. Lett. u 88 , 061101 (2002).
- [13] U. Leonhardt and P. Piwnicki, Phys. Rev. Lett. 84 , 822 (2000).
- [14] T. G. Philbin, C. Kuklewicz, S. Robertson, S. Hill, F. K¨ onig, and U. Leonhardt, Science 319 , 5868 (2008).
- [15] M. Petev, N. Westerberg, D. Moss, E. Rubino, C. Rimoldi, S. L. Cacciatori, F. Belgiorno, and D. Faccio, Phys. Rev. Lett. 111 , 043902 (2013).
- [16] D. Gerace and I. Carusotto, Phys. Rev. B 86 , 144505 (2012).
- [17] H. S. Nguyen, D. Gerace, I. Carusotto, D. Sanvitto, E. Galopin, A. Lemaitre, I. Sagnes, J. Bloch and A. Amo, Phys. Rev. Lett. 114 , 036402 (2015).
- [18] J. Steinhauer, Nat. Phys. 10 , 864 (2014)
- [19] J. Steinhauer, Nat. Phys. 12 , 959 (2016)
- [20] G. Rousseaux, C. Mathis, P. Maissa, T. G. Philbin, and U. Leonhardt, New J. Phys. 10 , 053015 (2008).
- [21] F. Belgiorno, S. L. Cacciatori, M. Clerici, V. Gorini, G. Ortenzi, L. Rizzi, E. Rubino, V. G. Sala and D. Faccio, Phys. Rev. Lett. 105 , 203901 (2010)
- [22] E. Rubino, F. Belgiorno, S. Cacciatori, M. Clerici, V. Gorini, G. Ortenzi, L. Rizzi, V. Sala, M. Kolesik and D. Faccio, New J. Phys. 13 , 085005 (2011).
- [23] S. Weinfurtner, E. W. Tedford, M. C. J. Penrice, W. G. Unruh, G. A. Lawrence, Phys. Rev. Lett. 106 , 021302 (2011).
- [24] J. Drori, Y. Rosenberg, D. Bermudez, Y. Silberberg and Ulf Leonhardt Phys. Rev. Lett. 122 , 010404 (2019).
- [25] T. Torres, S. Patrick, A. Coutant, M. Richartz, E. W. Tedford and S. Weinfurtner, Nat. Phys. 13 , 833 (2017).
- [26] I. Carusotto and C. Ciuti, Rev. Mod. Phys. 85 , 299 (2013).
- [27] J. Klaers, J. Schmitt, F. Vewinger, and M. Weitz, Nature 468 , 545 (2010).
- [28] T. Frisch, Y. Pomeau, S. Rica, Phys. Rev. Lett. 69 , 1644 (1992).
- [29] Y. Pomeau, S. Rica, Phys. Rev. Lett. 71 , 247 (1993).
- [30] D. Vocke, K. Wilson, F. Marino, I. Carusotto, E. M.
- Wright, T. Roger, B. P. Anderson, P. Ohberg, and D. ¨ Faccio, Phys. Rev. A 94 , 013849 (2016).
- [31] C. Michel, O. Boughdad, M. Albert, P-E. Larr´ and M. e Bellec, Nature Comm. 9 , 2108 (2018).
- [32] N. Santi´, ˘ c A. Fusaro, S. Salem, J. Garnier, A. Picozzi and R. Kaiser Phys. Rev. Lett. 120 , 055301 (2018).
- [33] R. Y. Chiao, J. Boyce, Phys. Rev. A 60 , 4114 (1999).
- [34] D. Vocke, T. Roger, F. Marino, E. M. Wright, I. Carusotto, M. Clerici, and D. Faccio, Optica 2 , 484 (2015).
- [35] Q. Fontaine, T. Bienaim´ e, S. Pigeon, E. Giacobino, A. Bramati and Q. Glorieux, Phys. Rev. Lett. 121 , 183604 (2018)
- [36] F. Marino, Phys. Rev. A 78 , 063804 (2008).
- [37] F. Marino, M. Ciszak, and A. Ortolan, Phys. Rev. A 80 , 065802 (2009).
- [38] I. Fouxon, O. V. Farberovich, S. Bar-Ad, and V. Fleurov, Europhys. Lett. 92 , 14002 (2010).
- [39] M. Elazar, V. Fleurov, and S. Bar-Ad, Phys. Rev. A 86 , 063821 (2012)
- [40] S. Bar-Ad, R. Schilling, and V. Fleurov, Phys. Rev. A 87 , 013802 (2013).
- [41] M. C. Braidotti and C. Conti, arXiv:1708.02623 (2017).
- [42] D. Vocke, C. Maitland, A. Prain, K. E. Wilson, F. Biancalana, E. M. Wright, F. Marino and D. Faccio, Optica 5 , 1099 (2018).
- [43] M. Ornigotti, S. Bar-Ad, A. Szameit and Victor Fleurov Phys. Rev. A 97 , 013823 (2018).
- [44] A. Prain, C. Maitland, D. Faccio and F. Marino, Phys. Rev. D 100 , 024037 (2019).
- [45] The confinement of photons into a cavity such as in polariton BECs gives rise to gapped dispersion relations although with imaginary terms due to the intrinsically dissipative nature of the cavity configuration (see e.g.[16, 26]. Here we are referring to (formally) conservative analog gravity models.
- [46] M. Visser and S. Weinfurtner, Phys. Rev. D 72 , 044020 (2005).
- [47] F. Girelli, S. Liberati and L. Sindoni, Phys. Rev. D 78 , 084013 (2008)
- [48] S. Liberati, M. Visser and S. Weinfurtner, Phys. Rev. Lett. 96 , 151301 (2006)
- [49] S. Finazzi, S. Liberati and L. Sindoni, Phys. Rev. Lett. 108 , 071101 (2012).
- [50] R. W. Boyd, Nonlinear Optics (Academic, Orlando, 2008).
- [51] Y. V. Kartashov, V. A. Vysloukh, L. Torner, Opt. Lett. 34 , 1543 (2009).
- [52] M. Warenghem, J. F. Blach, J.F.Henninot, J. Opt. Soc. Am. B 25 , 1882 (2008).
- [53] A. Griesmaier, J. Stuhler, T. Koch, M. Fattori, T. Pfau, and S. Giovanazzi Phys. Rev. Lett. 97 , 250402 (2006)
- [54] I. Tikhonenkov, B. A. Malomed, A. Vardi, Phys. Rev. A 78 , 043614 (2008).
- [55] T. Roger, C. Maitland, K. Wilson, N. Westerberg, D. Vocke, E. M. Wright, and D. Faccio, Nat. Commun. 7 , 13492 (2016).
- [56] N. Ghofraniha, C. Conti, G. Ruocco, and S. Trillo, Phys. Rev. Lett. 99 , 043903 (2007).
- [57] A. Minovich, D. N. Neshev, A. Dreischuh, W. Krolikowski, and Y. S. Kivshar, Optics Lett. 32 , 1599 (2007).
- [58] A. M. Deykoon and G. A. Swartzlander, Jr., J. Opt. Soc.

- Am. B 18 , 804 (2001).
- [59] C. Rotschild, O. Cohen, O. Manela, M. Segev, and T. Carmon, Phys. Rev. Lett. 95 , 213904 (2005).
- [60] S. Weinfurtner, S. Liberati, M. Visser J.Phys. A 39 , 6807 (2006).
- [61] C. Conti, M. Peccianti and G. Assanto, Phys. Rev. Lett. 91 , 073901 (2003)
- [62] K.R.W. Jones and D. Bernstein, Class. Quantum Grav. 18 , 1513 (2001).
- [63] R. Bekenstein, R. Schley, M. Mutzafi, C. Rotschild and M. Segev, Nat. Phys. 11 , 872 (2015).
- [64] L. Diosi, Phys. Lett. A 105 , 199 (1984).
- [65] R. Penrose, Gen. Relat. Gravit. 28 , 581 (1996).
- [66] M. Bahrami, A. Grossardt, S. Donadi, A. Bassi, New J. Phys. 16 , 115007 (2014).
- [67] P.H. Chavanis, Eur. Phys. J. Plus, 132 , 248 (2017).
- [68] S. Ghosh and N. Chakrabarti Phys. Rev E 84 , 046601 (2011).
- [69] On the basis of the fluid equations (3)-(4), a nearlyconstant density would also imply a nearly-homogeneous flow. Nevertheless, there are specific situations in which the simplifying assumption to consider a nearly-constant density while letting the flow vary provides an useful model of the system (e.g. in the case of vortex flows when the locations of the event horizons and the ergosphere occur in regions of nearly constant density, far away from the vortex core).
- [70] M. Novello and E. Goulart, Class. Quant. Grav. 28 , 145022 (2011).
- [71] E. Goulart, M. Novello, F. T. Falciano, and J. D. Toniato,
- Class. Quantum Grav. 28 , 245008 (2011)
- [72] F. Marino, C. Maitland, D. Vocke, A. Ortolan and D. Faccio, Sci. Reports 6 , 23282 (2016).
- [73] A. Belenchia, S. Liberati and A. Mohd, Phys. Rev. D 90 , 104015 (2014).
- [74] R. M. Wald, General Relativity University of Chicago Press, Chicago IL (1984).
- [75] W. Greiner, Relativistic Quantum Mechanics. Wave Equations , 3rd Edition, Springer-Verlag (2000).
- [76] We remark that the weak-field limit of GR in (2 + 1)dimensions (i.e. the dimensionality of our system) does not correspond to 2D Newtonian gravity. In two space dimensions all vacuum solutions of Einstein's equations are flat and even in the presence of matter, the spacetime curvature remains null everywhere except in the interior of the matter distribution. Similarly to other analogue gravity models, such a restiction does not apply to our system whose spacetime metrics are not solutions of the Einstein field equations. As a consequence, the Newtonian-like dynamics here discussed should thus be hypothetically considered as the limit of a theory that strongly differs from Einstein's gravity.
- [77] R. Balbinot, A. Fabbri, S. Fagnocchi, R. Parentani, Riv. Nuovo Cim. 28 , 1 (2005).
- [78] S. Giorgini, L.P. Pitaevskii, S. Stringari, J. Low Temp. Phys. 109 , 309 (1997).
- [79] S. Giorgini, Phys. Rev. A 61 , 063615 (2000).
- [80] L. Salasnich, Appl. Sci. 8 , 1998 (2018).
Francesco Marino
Massive phonons and gravitational dynamics in a photon-fluid model
We theoretically investigate the excitation dynamics in a photon-fluid with
both local and nonlocal interactions. We show that the interplay between
locality and an infinite-range nonlocality gives rise to a gapped Bogoliubov
spectrum of elementary excitations which, at lower momenta, correspond to
massive particles (phonons) with a relativistic energy-momentum relation. In
this regime and in the presence of an inhomogeneous flow the density
fluctuations are governed by the massive Klein-Gordon equation on the acoustic
metric and thus propagate as massive scalar fields on a curved spacetime. We
finally demonstrate that in the non-relativistic limit the phonon modes behave
as self-gravitating quantum particles with an effective Schr\"{o}dinger-Newton
dynamics, although with a finite-range gravitational interaction and a non-zero
cosmological constant. Our photon-fluid represents a viable alternative to BEC
models for "emergent-gravity" scenarios and offers a promising setting for
analogue simulations of quantum gravity phenomenology and semiclassical
gravity.