after several decades during which the study of systems of many interacting particles has focused on matter fluids such as liquid helium , electron gases in solid - state materials , ultracold atom vapors , nuclear - matter fluids , or quark - gluon plasmas in colliders , photon propagation in suitably designed nonlinear optical systems is presently attracting a growing interest as a novel platform to investigate the physics of interacting bose gases , the so - called _ quantum fluids of light _ ( see ref .
@xcite for a review ) .
the interactions between the photons constituting the fluid of light are mediated by the kerr ( @xmath0 ) optical nonlinearity of the underlying medium .
so far , numerous experimental observations have been performed for semiconductor - planar - microcavity geometries , including the demonstration of a bose - einstein - like condensation @xcite , of a superfluid flow @xcite , and of the hydrodynamic nucleation of nonlinear excitations such as solitons @xcite and quantized vortices @xcite in dilute photon gases . in the meanwhile
, active investigations have addressed the possibility of realizing systems characterized by very large optical nonlinearities .
correspondingly , the induced photon - photon interactions are expected to introduce strong quantum correlations within the fluid of light , generating in turn new quantum phases .    in this quest @xcite , researchers have been faced to ( at least ) two main difficulties . on the one hand , obtaining strong enough nonlinearities in scalable systems to study the dynamics of a strongly interacting photon gas in a spatially extended system turns out to be a major experimental challenge . on the other hand , the dynamics of a light field in devices based on cavities
is intrinsically a driven - dissipative one , which introduces severe complications in the theoretical description of such systems and which is typically very detrimental for the study of purely quantum features .
an alternative platform for studying many - body physics in photon fluids is based on light propagating in a bulk nonlinear ( of kerr type ) optical medium . from elementary classical optics
, it is well known that the paraxial propagation of a spectrally narrow beam of light can be described within the so - called paraxial and slowly - varying - envelope approximations ( see , e.g. , refs .
@xcite ) by a nonlinear wave equation formally analogous to the gross - pitaevskii equation for the order parameter of a dilute bose - einstein condensate @xcite . as originally pioneered in refs .
@xcite , this framework naturally translates upon quantization to a many - body quantum nonlinear schrdinger formalism with the roles of the propagation coordinate and time exchanged .
this framework has been used in a number of theoretical works where laser - physics problems have been reformulated in the hydrodynamics language @xcite , including , e.g. , the investigation of superfluid - like behaviors in the flow of a photon fluid @xcite , of nonlinear phenomena with light waves @xcite , and of the so - called acoustic hawking radiation @xcite . from the experimental point of view
, numerous works have been devoted to the study of nonlinear features that may appear in these systems , with a special attention dedicated to their relation to hydrodynamics and superfluidity aspects @xcite .
a major first step in the very quantum direction of realizing a gas of strongly interacting photons in a propagating geometry has been recently reported using an optically dressed atom gas in the rydberg - eit regime as a bulk nonlinear medium : this has allowed for the experimental observation of strongly repulsing photons @xcite and , soon after , of two - photon molecular bound states @xcite .
these remarkable experimental advances call for a theoretical approach that is able to describe in its full generality the many - body dynamics of strongly interacting photons propagating in a cavityless configuration .    building atop the pioneering theoretical works @xcite and refs .
@xcite , we report here a fully general quantum field theory of the propagating photon fluid . in this approach , the roles played by the optical - axis coordinate @xmath1 and the time parameter @xmath2 are exchanged : light propagation in the @xmath1 direction is naturally described in terms of evolution equations while the @xmath2 direction corresponds to a third spatial dimension in addition to the transverse @xmath3 and @xmath4 directions . in a paraxial - propagation configuration , light diffraction provides an effective mass to the photons in the @xmath5 plane and chromatic dispersion leads to a  typically different  effective mass in the @xmath2 direction . as usual
, the kerr nonlinearity of the medium gives rise to photon - photon interactions .
in contrast to microcavity configurations where driving and dissipation play a major role in the photon - fluid dynamics @xcite , in paraxial - propagation geometries , the quantum fluid of light follows a fully conservative hamilton dynamics starting from an initial condition determined by the incident light field and its coherence properties .
the quantum state of the optical field at the end of the evolution is experimentally accessible _ via _ a measurement of the statistical properties of the transmitted light emerging from the nonlinear optical medium .    as a most remarkable example of application of this formalism
, we then study the transmission of a coherent light across a finite slab of weakly nonlinear medium .
in this very simple configuration , the photons experience a pair of sudden jumps of the interaction parameter upon crossing the front and the back faces of the nonlinear medium . as a result of these two _ quantum quenches _ , the fluid of light gets excited and its quantum state after the conservative evolution across the nonlinear material can be reconstructed from the statistical properties of the transmitted light . in the weak - nonlinearity regime , the main excitation process consists in the emission of pairs of correlated counterpropagating bogoliubov phonons , which reflects in peculiar features in the intensity distribution and in the near- and the far - field two - body correlation functions . in its simplicity , this example illustrates the power of the conservative propagation dynamics in view of generating , detecting , and manipulating strongly correlated quantum phases of matter in photon gases , as well as of studying quantum dynamical features of many - body systems @xcite .
the paper is organized as follows .
first of all , in sec .
[ sec : classicalwaveequation ] , we review the classical propagation equation of a paraxial beam of light in a cavityless nonlinear optical medium of kerr type . on this basis , we present in sec .  [
sec : quantumtheory ] a general theory that makes it possible to describe the evolution of the quantum optical field for generic values of the nonlinear interaction parameter . in sec .
[ sec : bogoliubovtheoryofquantumfluctuations ] , we discuss how the relatively small fluctuations superimposing upon a coherent - light field in the weak - nonlinearity regime can be treated within the framework of the bogoliubov theory of dilute bose - einstein condensates . as an application of this formalism ,
we investigate in sec .
[ sec : responsetoquantumquenchesofthekerrnonlinearity ] the two - body quantum correlations resulting from the propagation of a laser beam across a slab of weakly nonlinear medium and interpret their features in terms of a dynamical casimir emission of bogoliubov collective excitations in a temporally modulated quantum fluid of light .
finally , in sec .
[ sec : conclusionsandoutlooks ] , we draw our conclusions and give prospects to the present work .
we consider the propagation of a laser wave in a dispersive and inhomogeneous kerr dielectric for which the ( frequency - dependent ) electric susceptibility reads [ one writes down @xmath6 , with @xmath7 @xmath8 where the homogeneous contribution @xmath9 takes into account the chromatic dispersion of the medium , the linear modulation @xmath10 comes from the existence of spatial inhomogeneities or of an optical confinement , and @xmath11 is the kerr nonlinear shift of the susceptibility , proportional to the local electric intensity @xmath12 , i.e. , to the square modulus of the ( complex representation of the ) electric field . for simplicity s sake , the dielectric is assumed to be devoid of free charge carriers and nonmagnetic .
we finally suppose that the optical field maintains its polarization in the course of its propagation in the medium so that a scalar approach is valid . as explained , e.g. , in ref .
@xcite , this is possible for paraxial beam of light and provided @xmath13 slowly varies on space scales of the order of the optical wavelength .
spin - orbit - coupling effects resulting from significant deviations from the optical axis and significant spatial variations of the electric susceptibility will be subjected to a future publication @xcite .    introducing the envelope @xmath14 of the laser - wave electric field @xmath15 oscillating at the angular frequency @xmath16 as @xmath17 where `` c.c . ''
stands for `` complex conjugate , '' and making use of the standard paraxial and slowly - varying - envelope approximations ( see , e.g. , refs .
@xcite ) , the maxwell equations supplemented by eq .   lead to the following classical wave equation for @xmath14 : @xmath18 in this equation , @xmath19 denotes the nabla operator in the @xmath20 plane and the functions @xmath21 and @xmath22 are defined as @xmath23=-\beta_{0}\,[\delta\chi(\mathbf{x}),\chi^{(3)}(\mathbf{x})]/[2\,(1+\chi_{0})]$ ] , where @xmath24 ; the parameters @xmath25 , @xmath26 , and @xmath27 are respectively the propagation constant @xmath28^{1/2}\,\omega / c$ ] ( @xmath29 is the vacuum speed of light ) of the laser wave in the @xmath30 direction , the group velocity @xmath31^{-1}$ ] of the photons in the medium , and the group - velocity dispersion @xmath32 evaluated at the carrier s angular frequency @xmath16 .    the hydrodynamic interpretation of the propagation equation is mostly well known in the limiting case of a purely monochromatic wave at @xmath16 @xcite and has offered a transparent physical interpretation to a number of nonlinear - optics experiments @xcite . in this time - independent case ,
the first and second derivatives of the envelope @xmath33 with respect to @xmath2 vanish , in such a way that the propagation of the optical field decribed by eq .
recovers the mean - field dynamics of the order parameter of a dilute two - dimensional bose - einstein condensate , @xmath34 playing the role of some external potential and @xmath35 corresponding to the effective two - dimensional boson - boson interaction constant .
on the other hand , only a very few studies so far @xcite have investigated the consequences of this analogy in time - dependent regimes , where , by extension , the propagation of the photon fluid in the positive-@xmath1 direction can be written in terms of a time evolution in a three - dimensional @xmath36 space where the physical time parameter @xmath2 plays the role of a third spatial coordinate in addition to the transverse variables @xmath3 and @xmath4 . in the following , a special attention will be paid to the new hydrodynamic features that originate from this @xmath2 dependence .    in order to facilitate @xmath1 ( @xmath2 ) to be viewed as a true time ( space ) parameter
, we introduce the following coordinates : @xmath37 respectively homogeneous to a time and a length . in these new variables ,
the paraxial - propagation equation takes the form of a time - dependent gross - pitaevskii equation : @xmath38 where the electric - field envelope @xmath33 now has to be considered as a function of @xmath39 and @xmath40 and @xmath41 as functions of @xmath42 .
this @xmath43 dependence of @xmath40 and @xmath41 corresponds in our language to a temporal dependence , which , as we shall see in the following , opens the way to the study of quantum - quench physics in the framework of paraxial optics . on the other hand
, the fact that the medium properties do not depend on the physical time @xmath2 implies that @xmath40 and @xmath41 are independent on the spatial @xmath44 coordinate .    in the gross - pitaevskii - like equation , the rigid - drift ( in the @xmath45 direction )
term @xmath46 originates from the group velocity of the photons in the medium and the kinetic operator @xmath47 involves a contribution in the ( actual - time ) @xmath45 direction in addition to the usual one in the @xmath20 plane .
the anisotropy of the `` mass '' tensor @xmath48 in eq .   is a natural consequence of the different origins of the effective masses @xmath49 , @xmath50 in the ( spatial ) @xmath3 , @xmath4 directions and @xmath51 in the ( temporal ) @xmath44 one : the former are due to diffraction while the latter originates from dispersion ; note that in vacuum , @xmath52 and so @xmath51 is infinite , by definition .    when the medium is characterized by an anomalous group - velocity dispersion , that is , when @xmath53 , one has @xmath54 and so the mass matrix @xmath55 is positive in all the directions . in the following ( see sec .
[ sec : bogoliubovtheoryofquantumfluctuations ] for a detailed discussion ) , we shall see that the dynamical stability of the fluid of light requires such a negative @xmath56 , but also repulsive photon - photon interactions , which sets @xmath57 and therefore @xmath58 , that is , that the kerr nonlinearity is self - defocusing .
in order to describe the quantum features of a light beam propagating in the nonlinear medium considered in sec .  [ sec : classicalwaveequation ] , the classical photon field @xmath59 verifying the paraxial wave equation must be replaced with a quantum field operator satisfying suitable boson commutation relations . to this purpose , in the present work
, we perform a canonical quantization @xcite of the classical field theory from which eq .   may be derived .
the procedure first consists in rewriting the evolution equation in lagrangian and then hamiltonian form ( sec .
[ subsec : hamiltonianformulation ] ) .
section [ subsec : normalizationconstant ] is dedicated to the precise determination of the global multiplicative constant appearing in the lagrangian and the hamiltonian of the paraxial - propagation problem ; while it plays no role at the classical level , it starts having a crucial importance upon quantization .
the quantization is finally accomplished in sec .
[ subsec : canonicalquantization ] by replacing the conjugate fields of the classical hamiltonian theory with quantum field operators obeying equal-@xmath43 bosonic commutation relations , standardly deduced from the canonical poisson - bracket relations .
of course , similar quantized wave equations describing paraxial light propagation in nonlinear optical media have been considered in the past ( see , e.g. , refs .
@xcite ) and applied to concrete problems related to quantum soliton propagation @xcite : as a crucial addition to these works , we do not restrict our attention to one - dimensional fiber geometries for which only the @xmath45 coordinate matters @xcite , but fully take into account the dynamics of the optical field in the transverse @xmath5 plane .
furthermore , an explicit expression for the normalization constant appearing in the boson commutators is provided .
in contrast to works on the quantization of paraxial electromagnetic fields @xcite , our approach is able to naturally describe spatiotemporal light propagation and to include the possible existence of spatial inhomogeneities and/or of an optical confinement as well as of an optical nonlinearity . as we shall see later
, all these features will be crucial for the theory to be applicable to the problem considered in sec .
[ sec : responsetoquantumquenchesofthekerrnonlinearity ] , which requires an accurate description of the peculiar spatial , angular , and spectral correlations displayed by the laser beam emerging from the back face of the nonlinear medium .    at this stage
, the reader who is just interested in knowing the basic elements structuring the quantum theory in order to apply the latter to concrete quantum - optics problems can skip secs .
[ subsec : hamiltonianformulation ] and [ subsec : normalizationconstant ] and go directly to sec .
[ subsec : canonicalquantization ] , where one finds the commutation relations of the quantum field operator associated to the classical photon field @xmath59 [ eqs .  ] , the many - body quantum hamiltonian of the paraxial - propagation problem [ eq .  ] , and
the corresponding heisenberg evolution equation [ eq .  ] .
let us introduce the lagrangian @xmath60=\int d\mathbf{r}\;\mathcal{l}(\mathbf{r},\tau),\ ] ] with the lagrangian density @xmath61.\right.\end{aligned}\ ] ] it is immediate to notice that the euler - lagrange equations of motion @xmath62 and @xmath63 deduced from eqs .   and coincide with the evolution equation and its complex conjugate , respectively .
the global normalization factor @xmath64 in the definition of the lagrangian density @xmath65 will be rigorously determined in sec .
[ subsec : normalizationconstant ] on the basis of microscopic calculations .
it plays obviously no role at the classical level but is on the contrary crucial to derive the exact commutation relations of the quantum field operators .
equation being a first - order differential equation in @xmath43 , the data of the electric - field envelope @xmath59 at a given time @xmath66 is sufficient to determine the subsequent evolution of the system . as a result ,
the lagrangian @xmath67 , functional of @xmath68 , @xmath33 , and their derivatives with respect to @xmath43 , contains an overabundant number of dynamical variables .
thus , before searching for the conjugate momenta and moving on to the hamiltonian formalism , it is convenient to eliminate the redundant dynamical variables from the lagrangian .
following ref .
@xcite , one starts by rewriting the lagrangian density @xmath65 as a function of @xmath69 $ ] , @xmath70 $ ] , and their spatiotemporal derivatives .
in particular , the term involving @xmath71 and @xmath72 in eq .   reads @xmath73 .
thus , by adding @xmath74 to @xmath67 , we get a lagrangian @xmath75 which does not depend on @xmath76 and which is by definition equivalent to @xmath67 .
its density @xmath77 reads @xmath78+\frac{\partial\mathcal{e}_{1}}{\partial\zeta}\,\mathcal{e}_{2}-\mathcal{e}_{1}\,\frac{\partial\mathcal{e}_{2}}{\partial\zeta}\right .
\\ \label{newlagrangiandensity } & \left.-\,u(\mathbf{x}_{\perp},\tau)\,(\mathcal{e}_{1}^{2}+\mathcal{e}_{2}^{2})-\frac{g(\mathbf{x}_{\perp},\tau)}{2}\,(\mathcal{e}_{1}^{2}+\mathcal{e}_{2}^{2})^{2}\bigg\}.\right.\end{aligned}\ ] ] by means of the euler - lagrange equation relating to @xmath79 , that is , @xmath80 , one may express @xmath79 as a function of @xmath81 and its time derivative . inserting this expression of @xmath79 into @xmath75 , we finally obtain a lagrangian @xmath82 which only involves the dynamical variable @xmath83 and its time derivative @xmath84 .
now let us denote by @xmath85 the conjugate momentum of @xmath81 . by definition , @xmath86}{\delta[\partial_{\tau}\mathcal{e}_{1}(\mathbf{r},\tau)]}.\ ] ] along the extremal `` path , ''
@xmath80 , and as a consequence , @xmath87=\delta l'(\tau)/\delta[\partial_{\tau}\mathcal{e}_{1}(\mathbf{r},\tau)]$ ] .
according to eq .
, this yields @xmath88 as conjugated fields , @xmath81 and @xmath85 obey the canonical relations    [ poissonbracketrelations ] @xmath89    where @xmath90 is the equal-@xmath43 poisson bracket .
we are now able to write the hamiltonian of the system , @xmath91 in eq .
, the hamiltonian @xcite density @xmath92 is by definition obtained by legendre transforming the lagrangian density @xmath93 , @xmath94 , which , in the light of eqs .   and and recalling that @xmath95 , leads to @xmath96.\right.\end{aligned}\ ] ]      before moving on to the canonical quantization of the classical field theory ( [ lagrangiandensity ] , [ hamiltoniandensity ] ) ,
let us determine the multiplicative constant @xmath64 in terms of the optical parameters of the electromagnetic wave .
@xmath64 was introduced as a global factor in eq .   to ensure that @xmath65 has the dimension of an energy per unit volume , as it has to be according to eq .  .
while such a normalization constant plays no role at the classical level , since it cancels out in the euler - lagrange equations of motion , it becomes important upon quantization as it determines the actual spacing between the energy levels of the system .
the role of the quantized action in the bohr - sommerfeld quantization rules is the most well - known examples of such a dependence .
a possible strategy to fix @xmath64 is to determine the total energy @xmath97 carried by the laser wave of complex electric field @xmath98 [ @xmath99 ; see eq .  ] in the simple case where @xmath100 and @xmath101 , and compare the latter with the well - known formula of classical electrodynamics for the time - averaged energy of a quasimonochromatic electromagnetic field propagating through a dispersive , homogeneous , linear , nonmagnetic medium , that is ( see , e.g. , ref .
@xcite ) , @xmath102}{d\omega}(\omega_{0})\,|\tilde{e}|^{2}+\frac{|\tilde{b}|^{2}}{\mu_{0}}\bigg\}.\ ] ] in eq .
, @xmath103 $ ] denotes the ( frequency - dependent ) permittivity of the medium , @xmath104 being the one of free space , @xmath105 , where @xmath106 is the slowly varying envelope of the magnetic field ( deduced , e.g. , from maxwell - faraday s law ) , and @xmath107 is the vacuum permeability .    in order to get an expression for the physical energy @xmath97 in our lagrangian formalism
, one has to perform a legendre transformation of the lagrangian with respect to the actual time coordinate , @xmath2 , instead of the `` propagation one , '' @xmath43 . as the jacobian determinant @xmath108 and
since the action @xmath109=\int d\tau~l(\tau)=\int d\mathbf{r}\,d\tau~\mathcal{l}(\mathbf{r},\tau)$ ] of the optical system can be alternatively defined in the @xmath110 coordinates as @xmath109=\int dt~\tilde{l}(t)=\int d\mathbf{x}\,dt~\tilde{\mathcal{l}}(\mathbf{x},t)$ ] , the lagrangian density stays invariant under the coordinate transformations : @xmath111 .
thus , one finds from eq .   that the conjugate momentum @xmath112 [ @xmath113 of @xmath14 [ @xmath114 in the @xmath110 coordinates is given by    [ conjugatemomentumtilde ] @xmath115}{\delta[\partial_{t}\mathcal{e}(\mathbf{x},t ) ] } \\ \label{conjugatemomentumtilde2 } & = \mathcal{n}\,\bigg[\frac{i}{2\,v_{0}}\,\mathcal{e}^{\ast}(\mathbf{x},t)+\frac{d_{0}}{2}\,\frac{\partial\mathcal{e}^{\ast}}{\partial t}(\mathbf{x},t)\bigg],\end{aligned}\ ] ]    as a consequence of which the hamiltonian density @xmath116-\tilde{\mathcal{l}}(\mathbf{x},t)$ ] associated to @xmath117 in the @xmath110 coordinates reads  when @xmath100 and @xmath101as @xmath118.\ ] ] making the substitution @xmath119 into eq .  , one deduces the hamiltonian @xmath120 as a functional of the total electric field @xmath121 : @xmath122.\right.\end{aligned}\ ] ] the temporal evolution of the envelope @xmath14 is given by the hamilton equation of motion @xmath123 , @xmath124 being the equal-@xmath2 poisson bracket  not to be confused with the equal-@xmath43 poisson bracket @xmath90 defined in sec .
[ subsec : hamiltonianformulation ] . inserting @xmath119 into this evolution equation ,
one ends up with the following equation for the total field @xmath121 :    [ hamiltonequationofmotion ] @xmath125    where , making use of @xmath126 , @xmath127 is the hamiltonian encoding the time evolution of the total electric field @xmath121 . by means of eqs .   and
, we can write it in the form @xmath128 where @xmath129.\right.\end{aligned}\ ] ]    within the framework of the paraxial and slowly - varying - envelope approximations , one can check that @xmath130 defined in eq .
gives a negligible contribution to the total energy @xmath97 of the electromagnetic wave . by comparing eqs .   and , one finally obtains the normalization constant @xmath64 in terms of the angular frequency @xmath16 and the propagation constant @xmath131 of the laser in the medium : @xmath132 brief discussions on the constant normalizing the lagrangian and the hamiltonian in one - dimensional fiber geometries were given in ref .
@xcite and , in a bit more detailed way , in ref .
@xcite . a physical interpretation of expression
will be given in the next section .      in order to carry out the canonical quantization of the ( classical ) hamiltonian theory presented in sec .
[ subsec : hamiltonianformulation ] , one replaces the conjugated fields @xmath81 and @xmath85 with quantum field operators @xmath133 and @xmath134 ( by choice , in the heisenberg picture ) and the poisson bracket @xmath90 with the commutator @xmath135/(i\,\hbar)$ ] , where @xmath136 is the reduced planck constant . on doing so ,
become    [ canonicalcommutationrelations ] @xmath137=i\,\hbar\;\delta(\mathbf{r}-\mathbf{r}'),\right .
\\ \label{canonicalcommutationrelations2 } & \left.[\hat{\mathcal{e}}_{1}(\mathbf{r},\tau),\hat{\mathcal{e}}_{1}(\mathbf{r}',\tau)]=0,\right .
\\ \label{canonicalcommutationrelations3 } \text{and}\quad&\left.[\hat{\pi}(\mathbf{r},\tau),\hat{\pi}(\mathbf{r}',\tau)]=0,\right.\end{aligned}\ ] ]    from which and thanks to eq .
one deduces the following equal - time , that is , equal-@xmath138 , commutation relations :    [ commutationrelations ] @xmath139=\frac{\hbar\,v_{0}}{\mathcal{n}}\,\delta(\mathbf{r}-\mathbf{r}')\right .
\\ \label{commutationrelations2 } \text{and}\quad&\left.[\hat{\mathcal{e}}(\mathbf{r},\tau),\hat{\mathcal{e}}(\mathbf{r}',\tau)]=0,\right.\end{aligned}\ ] ]    @xmath140 being the hermitian conjugate of @xmath141 : @xmath142 .
it is worth noting that , instead of reducing the overabundant number of dynamical variables in the lagrangian @xmath67 in order to avoid dealing with fields which are not independent from each other in the transition from the lagrangian formalism to the hamiltonian one and in the canonical quantization _ via _ the poisson bracket
, we could also have implemented the more basic and robust dirac - bergmann quantization procedure @xcite to pass from the classical theory to the quantum one , as recently used in ref .  @xcite . as detailed in footnote @xcite ,
dirac s procedure produces the same equal - time commutation relations for the quantum field operators as the ones deduced from the canonical quantization method adopted in this work , which fully validates our approach .    from the definition @xmath143 ,
one notes that the commutation rules involve quantum field operators at the same `` propagation '' time @xmath138 but at different physical times @xmath2 and @xmath144 .
such a writing is based on the coordinate transformation @xmath145 $ ] , which is of course legitimate only in the paraxial and slowly - varying - envelope approximations and if back - scattered waves are assumed not to exist , i.e. , if light propagation is assumed to occur only in the positive-@xmath1 direction .
such a quantization procedure imposing equal-@xmath1 and different-@xmath2 ( instead of equal-@xmath2 and different-@xmath1 ) canonical commutation relations was pioneered in refs .
@xcite and critically discussed in ref .
it is also important to insist on the fact that the slow spatiotemporal variation condition that is assumed for the classical electric - field envelope @xmath14 directly applies to its quantum counterpart @xmath141 : a discussion of such a condition in a quantum framework is illustrated in refs .
@xcite .    by inserting the explicit expression of the normalization parameter @xmath64 into the commutation rule
, one easily gets that the multiplicative constant in the right - hand side of eq
.   reads @xmath146 $ ] , yielding a simple interpretation of @xmath64 in terms of the energy density of a single - photon wavepacket of total energy @xmath147 .    in the light of eqs .   and , the quantized hamiltonian of the system reads , after considering the normal order , @xmath148.\right.\end{aligned}\ ] ] equation corresponds to the many - body quantum hamiltonian describing the evolution in time @xmath138 of a many - photon laser beam propagating through the bulk inhomogeneous and nonlinear optical medium of electric susceptibility . in the dielectric ,
the position of a point is referenced by the coordinates @xmath20 and @xmath44 .
the two first contributions are the kinetic terms in the transverse @xmath5 plane and in the @xmath45 direction with different effective masses ( as discussed in sec .
[ sec : classicalwaveequation ] ) , the second line describes the rigid global drift along the @xmath45 axis due the group velocity of light in the medium , and the two last terms respectively account for the spatial modulation of the electric susceptibility and for the two - photon interactions mediated by the kerr nonlinearity of the dielectric .    in the theory of ultracold bose fluids ,
contact interactions are usually considered in place of the actual  but much more complicated  two - body interactions .
this approximation is very helpful in simplifying the many - body quantum problem and is well accurate as long as the inter - particle distance is much larger than the range of the boson - boson interactions ( see ref .
@xcite ) . here , the assumed local form of the kerr optical nonlinearity automatically leads to contact - like interactions between the photons of the light beam , that is , no dilutness condition for the photon gas is required to get the four - field interaction term @xmath149 in eq .  .
since no hypothesis is made on the intensity of the photon - photon interaction parameter @xmath41 , the many - body hamiltonian can describe a quantum fluid of weakly interacting photons , that is , in the gross - pitaevskii / bogoliubov regime , as well as a strongly interacting one . in what follows
, we will focus our attention on the weak - nonlinearity regime , in which the quantum fluctuations of the fluid of light can be described within the framework of the well - known bogoliubov theory of dilute bose - einstein condensates ( see , e.g. , refs .
motivated by the intense experimental investigations that are presently in progress @xcite , the strong - interaction regime will be the subject of future works @xcite .    for the sake of completeness
, it is useful to explicitly write the evolution equation of the operator @xmath141 . in the heisenberg picture
, it is obtained from the system s hamiltonian as @xmath150 $ ] .
using eq .
and taking advantage of the commutation relations , this gives @xmath151 equation , which governs the time evolution of @xmath141 in the @xmath152 space , is simply the quantized version of the classical equation .
as originally pointed out in ref .
@xcite for a one - dimensional waveguide geometry , it has the form of a quantum nonlinear schrdinger equation .
in an illuminated dielectric medium devoid of free charges as the one considered in this paper , quantum noise of the electromagnetic field only arises from the quantum uncertainty of the optical field , that is , in more physical terms , from the discreteness of the photon . in the case of a strongly coherent light beam propagating across a weakly nonlinear three - dimensional bulk medium ,
quantum noise is typically small and can be described in terms of weak - amplitude quantum fluctuations oscillating on top of a strongly classical wave .
mutuating well - known results from the theory of weakly interacting ultracold atomic gases @xcite , one may develop a bogoliubov - like theory based on an expansion of the envelope operator @xmath141 of the form @xmath153 in this expression , the classical field @xmath59 , which satisfies the gross - pitaevskii - like equation , corresponds to the coherent component of the electric - field envelope and @xmath154 is a small quantum correction to @xmath59 . as the whole quantum nature of the optical field
is captured in the fluctuation operator @xmath154 , the equal-@xmath43 commutation relations then totally transfer to the latter , giving    [ commutationrelationsfluctuation ] @xmath155=\frac{\hbar\,v_{0}}{\mathcal{n}}\,\delta(\mathbf{r}-\mathbf{r}')\right . \\ \label{commutationrelationsfluctuation2 } \text{and}\quad&\left.[\delta\hat{\mathcal{e}}(\mathbf{r},\tau),\delta\hat{\mathcal{e}}(\mathbf{r}',\tau)]=0.\right.\end{aligned}\ ] ]    linearizing the heisenberg equation of motion with respect to @xmath154 and its hermitian conjugate , one readily gets the so - called bogoliubov - de gennes equation @xmath156 which encodes the time evolution of the quantum fluctuation @xmath154 in the three - dimensional @xmath152 space and is in fact the heart of the bogoliubov approach . with respect to similar ( classical ) equations considered in the literature @xcite
, this equation explicitly includes the @xmath44 dependence of the field and the corresponding effective mass , given by the group - velocity dispersion @xmath157 at @xmath158 .    in sec .
[ subsec : spatiallyhomogeneoussystem ] , we review the bogoliubov theory of linearized fluctuations in position- and time - independent configurations . although well established and known in the context of matter fluids , it is important to quickly review it within the nonlinear - propagating - geometry context because ( i ) of the nontrivial role of the effective mass @xmath159 in the dynamics of the luminous fluid in the temporal @xmath44 direction and ( ii ) of the existence of on - going experiments which aim at probing the phononic part of the bogoliubov excitation spectrum of a propagating fluid of light @xcite .      in the simplest case where the classical background field @xmath59 is at some time @xmath160 homogeneous in all the ( @xmath3 , @xmath4 , and @xmath45 ) directions ,
@xmath161 with @xmath162 , and when the nonlinear optical medium is spatially homogeneous with a constant kerr coefficient @xmath0 ,    [ uniformsituation ] @xmath163    analytical solutions for the classical wave equation and the bogoliubov - de gennes equation are available ( see , e.g. , refs .
@xcite ) .
in such a configuration , the electric - field envelope @xmath59 follows a simple harmonic evolution with a linearly - evolving ( in time @xmath43 ) global phase : @xmath164 , where @xmath165 . in the theory of
dilute bose - einstein condensates , the wavenumber @xmath166 correponds to the chemical potential of the bose gas @xcite .    in the homogeneous situation ,
the elementary excitations of the fluid of light are plane waves of wavevector @xmath167 in the @xmath168 space , the solution of the bogoliubov - de gennes equation obeying the mode expansion @xmath169.\right.\end{aligned}\ ] ] the mode operators @xmath170 satisfy the same-@xmath43 commutation relations    [ commutationbogoliubovoperators ] @xmath171=(2\pi)^{3}\,\frac{\hbar\,v_{0}}{\mathcal{n}}\,\delta(\mathbf{q}-\mathbf{q}')\right .
\\ \label{commutationbogoliubovoperators2 } \text{and}\quad&\left.[\hat{b}(\mathbf{q},\tau),\hat{b}(\mathbf{q}',\tau)]=0,\right.\end{aligned}\ ] ]    and evolve harmonically as @xmath172,\ ] ] with a frequency @xmath173 $ ] determined by the @xmath174-dependent wavenumber @xmath175},\ ] ] where @xmath176 the so - called bogoliubov amplitudes @xmath177 and @xmath178 [ @xmath179 and @xmath43 independent in the situation ] are finally given by @xmath180 by definition , they satisfy the normalization condition @xmath181 .
because of the photon - photon interactions ( @xmath182 ) , @xmath178 is nonzero : this indicates that the ground state of the bogoliubov theory differs from the trivial vacuum without particles in the @xmath183 modes @xcite .
modulo the @xmath184 contribution to the oscillation pulsation @xmath185 , originating from the presence of the drift term @xmath186 in the bogoliubov - de gennes equation , the wavenumber @xmath187 defined in eq
.   corresponds to the well - known bogoliubov dispersion relation for the elementary excitations propagating on top of a uniform dilute bose - einstein condensate at rest , the quadratic function @xmath188 given by eq .
playing the role of the single - particle kinetic energy of matter - wave superfluids .
note that this result stems from the conservative nature of the considered dynamics and is in contrast to the rich variety of dispersions predicted for driven - dissipative fluids of light in microcavity architectures @xcite . the graphical representation of @xmath187 is given in fig .
[ bogoliubovlaw ] for repulsive photon - photon interactions , i.e. , for @xmath189 , in the two cases @xmath190 .    in the anomalous - dispersion case ,
that is , when @xmath53 , the photon effective mass @xmath159 in the @xmath45 direction is  as the photon effective mass @xmath191 in the transverse @xmath5 plane  positive , and the `` kinetic energy '' @xmath188 stays as a consequence positive for any @xmath174 . if one also considers a self - defocusing nonlinearity , that is , if @xmath192 , the bogoliubov dispersion relation @xmath187 never acquires an imaginary part , which means that the photon - photon collision processes mediated by the underlying nonlinear medium do not give rise to unstable behaviors in the photon fluid .
( in units of @xmath166 ) of the elementary excitations propagating on top of a homogeneous quantum fluid of weakly interacting photons ( with @xmath192 ) at rest . in the anomalous - dispersion case [ @xmath53 , panels ( a ) and ( b ) ] , @xmath187 corresponds to the well - known bogoliubov law , linear at small @xmath174 s ( which is indicated by the black straight lines ) and quadratic at large @xmath174 s ( black parabola ) . in the normal - dispersion case [ @xmath193 , panels ( a ) and ( c ) ] ,
the effective mass @xmath51 in the @xmath45 direction is negative and @xmath187 has an imaginary component [ cf .
panel ( c ) ] : in this case , the quantum fluid of light is dynamically unstable . ]    in the case mentioned above where @xmath53 and @xmath192 , at low wavevectors , that is , when @xmath194 , where @xmath195 and @xmath196 are the healing lengths @xcite of the photon fluid in the @xmath20 plane and in the @xmath45 direction , respectively , the bogoliubov dispersion relation @xmath187 is phonon - like in all the directions [ as shown by the black straight lines in panels ( a ) and ( b ) of fig .
[ bogoliubovlaw ] ] , but with different `` sound velocities '' @xcite @xmath197 and @xmath198 : @xmath199 a pump - probe measurement of @xmath200 was recently reported in ref .
another experiment aiming at measuring @xmath201 in a one - dimensional waveguide configuration is presently in progress @xcite . in the opposite limit , that is , when @xmath202 , the dispersion @xmath187 of the elementary excitations takes a single - particle - like shape [ as illustrated by the black parabola in panels ( a ) and ( b ) of fig .
[ bogoliubovlaw ] ] : @xmath203 in the present paraxial - optics context , the hartree interaction term @xmath166 corresponds to the usual modification of the propagation constant due to the kerr nonlinearity of the underlying medium .
the situation is a bit more complicated in the case of a normal dispersion , that is , when @xmath193 : as presented in panel ( c ) of fig .
[ bogoliubovlaw ] , there exists in that case a range of wavevectors @xmath174 in the @xmath45 direction for which the dispersion law @xmath187 is purely imaginary ; this signals that the quantum fluid of light is dynamically unstable against the formation of a train of pulses ( see , e.g. , refs .
@xcite ) .    to complete the picture , we can also mention that for focusing nonlinearities ( @xmath204 ) the low - wavevector modes in the transverse @xmath5 plane are unstable for all @xmath56 : this gives rise to a modulational instability which , in the optical language , goes under the name of filamentation instability ( see , e.g. , refs .
@xcite ) . on the other hand , in this case
, the instability in the temporal direction only occurs for a anomalous group - velocity dispersion ( @xmath53 ) .
as a first example of application of the general quantum formalism developed in secs .  [
sec : quantumtheory ] and [ sec : bogoliubovtheoryofquantumfluctuations ] , we investigate in the present section the propagation of a wide laser beam across a slab of kerr medium immersed in vacuum ( see fig .  [ experimentalsetup ] ) . having in mind the reformulation of the propagation of the optical field along the @xmath1 axis in terms of a temporal evolution in time @xmath138 ,
it is straightforward to see how such a configuration constitutes a very simple realization of a pair of quantum quenches of the nonlinear interaction parameter in the system s hamiltonian .
as the optical nonlinearity is nonzero only inside the kerr material , the first ( second ) quench occurs at the entrance ( exit ) face of the dielectric , where the value of the photon - photon interaction constant suddenly jumps from @xmath205 to @xmath206 ( from @xmath206 to @xmath205 ) . while the present work focuses on the weak - nonlinearity regime , regime within which the quantum fluctuations of the fluid of light are accurately described by the bogoliubov theory of linearized fluctuations reviewed in the previous section , application of the general theory  to the strongly interacting regimes experimentally realized in refs .
@xcite is definitely possible and will be subjected to future works @xcite .
given the extreme simplicity of the proposed setup as compared to corresponding experiments using ultracold atomic vapors @xcite , this study demonstrates the promise of nonlinear optical systems to investigate different features of conservative quantum dynamics , including the response to quenches in the hamiltonian s parameters and the subsequent thermalization dynamics @xcite : in contrast to microcavity systems where cavity losses play a crucial role in the dynamics and quickly wash out quantum correlations @xcite , in the present propagating geometry , losses can be made arbitrarily small simply by choosing a suitable transparent dielectric medium .
building on decades of expertise in quantum - optics experiments , we discuss below how a detailed information on the response of the photon fluid to the quenches of the interaction parameter can be extracted from the statistical properties of the light emerging from the slab of nonlinear material , in particular in the real - space intensity - correlation signal ( sec .
[ subsubsec : correlationsinpositionspace ] ) and in the far - field angular distribution and two - body correlations of the transmitted light ( sec .
[ subsubsec : correlationsinmomentumspace ] ) .
an intuitive interpretation of our predictions in terms of a dynamical casimir emission of bogoliubov waves on top of a temporally modulated quantum gas of weakly interacting photons is provided and explained in detail , and the analogies and differences with standard four - wave - mixing experiments are discussed .    while working with spatially finite , e.g. , gaussian , beams and extending the whole theory to the general case where @xmath207 and @xmath208 depends on @xmath209 [ cf .
eq .  ] is computationally demanding but conceptually straightforward , for the sake of simplicity we will restrict our attention to homogeneous geometries : in an actual experimental implementation , this assumption requires working with continuous - wave fields with a wide top - hat spatial profile propagating in a homogeneous medium , so that a uniform bulk region can be identified in the fluid of light , as done in the recent experiment @xcite .
, constant kerr coefficient @xmath0 , and finite thickness @xmath210 in the @xmath1 direction is immersed in vacuum and illuminated by a wide monochromatic - plane - wave laser beam parallel to the @xmath1 axis . in such a configuration ,
the photon - photon interaction constant undergoes a sudden gate - shaped modulation along the optical axis .
the resulting quantum correlations appearing in the transmitted light reveal the spontaneous emission of correlated counterpropagating excitations in the photon fluid . ]    as sketched in fig .
[ experimentalsetup ] , the slab of nonlinear material consists in a homogeneous dielectric layer of linear electric susceptibility @xmath9 and uniform kerr coefficient @xmath0 .
it is assumed to be infinite in the @xmath3 and @xmath4 directions and to have a uniform thickness @xmath210 in the @xmath1 one .
the @xmath211 ( @xmath212 ) surface defines the front ( back ) interface between the dielectric and vacuum . the monochromatic laser beam of angular frequency @xmath16 which illuminates perpendicularly the medium is supposed to have a wide uniform profile in the transverse @xmath5 plane so that it can be legitimately seen as an infinite plane wave .
the propagation occurs along the increasing-@xmath1 axis and the propagation constant in vacuum is denoted as @xmath213 . at the angular frequency @xmath16 , the material is supposed to have an anomalous group - velocity dispersion ( @xmath53 ) and a negative kerr coefficient [ @xmath58 , i.e. , @xmath192 ; see eq .
] , so that the fluid of light is dynamically stable inside the nonlinear medium ( see sec .
[ subsec : spatiallyhomogeneoussystem ] ) .
back - propagating light waves originating from reflection on the interfaces separating the kerr material from vacuum would spoil the reformulation of the laser propagation in the positive-@xmath1 direction in terms of an effective time evolution . in order to avoid dealing with their existence ,
we assume that the @xmath214 surfaces of the nonlinear - dielectric layer are treated with a perfect antireflection coating .
since its characteristic thickness in the @xmath1 direction ( of the order of a few optical wavelengths ) is typically much shorter than all the other lengths of the problem , its effect on light transmission can be summarized as a simple boundary condition guaranteeing the conservation of the energy - flux density of the electromagnetic wave  the so - called poynting vector  at the @xmath214 interfaces of the dielectic slab .
mathematically speaking , this may be explicitly expressed as @xmath215 where @xmath216 is an arbitrarily - small positive parameter .
equation allows for matching the quantum fluctuations of the optical field in vacuum ( @xmath217 ) to the ones propagating in the bulk nonlinear medium ( @xmath218 ) ; see sec .
[ subsubsec : matching ] .        outside the kerr layer ( when @xmath217 ) , that is , in vacuum , the ( quantized ) electric field @xmath219 of the laser wave admits the usual fourier expansion ( see , e.g. , ref .
@xcite ) @xmath220}\,\hat{\alpha}(\mathbf{k})+\mathrm{h.c.},\ ] ] where `` h.c . '' stands for `` hermitian conjugate . '' in this equation , @xmath221 denotes the photon dispersion law in vacuum , ( we recall that ) @xmath104 is the vacuum permittivity , and the @xmath222 s [ @xmath223 s ] are the photon destruction ( creation ) operators in the state of wavevector @xmath224 , subject to the boson commutation relations    [ photoncommutationrelations ] @xmath225=(2\pi)^{3}\,\delta(\mathbf{k}-\mathbf{k}')\right . \\ \label{photoncommutationrelations2 } \text{and}\quad&\left.[\hat{\alpha}(\mathbf{k}),\hat{\alpha}(\mathbf{k}')]=0.\right.\end{aligned}\ ] ]    in vacuum , the second - quantized version of eq
.   reads @xmath226 , where @xmath213 is the propagation constant of the light wave along the positive-@xmath1 direction in free space .
thus , by identification , one deduces from eq .   that the envelope @xmath227 of the electric field in vacuum is given by @xmath228}{\varepsilon_{0}}}\right .
\\ \label{envelopevacuum } & \left.\times\,e^{i[\mathbf{k}_{\perp}\cdot\mathbf{x}_{\perp}+\delta k_{z}z-\delta\omega(\mathbf{k}_{\perp},\delta k_{z})t]}\,\hat{\alpha}(\mathbf{k}_{\perp},\delta k_{z}),\right.\end{aligned}\ ] ] equation in which we made the variable change @xmath229 and in which we defined @xmath230 and @xmath231 .
the @xmath232 s obey the commutation rules , assuming the exchanges @xmath233 and @xmath234 .
furthermore , @xmath235 , @xmath236 , and @xmath237 in eq .
are related through @xmath238 which is obtained by expanding at leading order the vacuum dispersion law @xmath221 around @xmath239 , that is , within the framework of the paraxial ( @xmath240 ) and slowly - varying - envelope ( @xmath241 and @xmath242 ) approximations considered in this paper . taking advantage of the dispersion law
, one can transform the integration over @xmath243 in eq .   into an integration over @xmath244 , yielding , at leading order in the paraxial and slowly - varying - envelope approximations , @xmath245}\,\hat{\alpha}[\mathbf{k}_{\perp},\delta k_{z}(\mathbf{k}_{\perp},\delta\omega)].\right.\end{aligned}\ ] ] in this way , the expression of the electric - field envelope outside the nonlinear medium
is easily inserted in the @xmath246 mapping language discussed in the previous sections .    defining in the integral    [ timeevolutionreformulation ]
@xmath247,\right.\end{aligned}\ ] ]    where @xmath248 denotes the vacuum value of the normalization constant , one easily shows that eq .
can be put into the form    [ fluctuationvacuumbis ] @xmath249    where , introducing the input / output quantum mode operators @xmath250\,h/2\}\,\hat{a}(\mathbf{q})$ ] , @xmath251\right .
\\ \label{ladderoperatorsinout } & \left.\times\,\hat{a}_{\left(\begin{smallmatrix}\mathrm{in}\\\mathrm{out}\end{smallmatrix}\right)}(\mathbf{q}).\right.\end{aligned}\ ] ] by construction , the ladder operators @xmath252 satisfy the equal-@xmath43 commutation relations [ use eqs .   and ]    [ photoncommutationrelationsbis ] @xmath253=(2\pi)^{3}\,\frac{\hbar\,c}{\mathcal{n}_{\mathrm{vac}}}\,\delta(\mathbf{q}-\mathbf{q}')\right .
\\ \label{photoncommutationrelationsbis2 } \text{and}\quad&\left.[\hat{a}(\mathbf{q},\tau),\hat{a}(\mathbf{q}',\tau)]=0.\right.\end{aligned}\ ] ]    in order to keep some homogeneity in the notations used outside and inside the nonlinear material , the coordinate referencing a point along the @xmath2 axis is conveniently defined in vacuum in the same way as in the kerr layer : @xmath44 [ eq .  ] , where @xmath254 denotes the group velocity in the kerr medium ( i.e. , in the @xmath218 region ) at the laser s angular frequency @xmath16 .
finally , note that the @xmath255 s [ @xmath256 s ] are related to the usual input ( output ) photon operators @xmath257 s [ @xmath258 s ] in the @xmath259 ( @xmath260 ) vacuum through eq .  .
in the specific configuration investigated here , with a coherent monochromatic - plane - wave incident light beam , the optical field in vacuum can be expanded as @xmath261 i.e. , as a classical background field @xmath262 ( with @xmath263 ) plus a small quantum - fluctuation term @xmath154 .
the latter may be decomposed according to the same plane - wave expansion as the one , eq .
, of the total electric - field envelope : the assumption of having a purely coherent incident beam translates into the condition that all the incident modes other than the coherent - pump one are in their vacuum state @xmath264 , by definition such that @xmath265 , which correspondingly reads in terms of the @xmath255 operators as @xmath266 .      inside the nonlinear medium ( when @xmath218 ) , the ( quantized ) electric - field envelope @xmath141 evolves in time @xmath138 in the @xmath143 space according to the heisenberg equation of motion with the prescription .    within the bogoliubov weakly nonlinear regime , the background electric - field envelope @xmath59 is of constant amplitude @xmath33 and of global phase @xmath165 ( see sec .
[ subsec : spatiallyhomogeneoussystem ] ) , where @xmath267 is the time which ( naturally ) initializes the evolution of the optical field in the medium . from the continuity of the poynting vector at the entrance face of the nonlinear material , i.e. , from eq .   for @xmath211 [ @xmath268
, one finds @xmath269 on the other hand , from the continuity of the poynting vector at the exit face of the medium , i.e. , from eq .   for @xmath212 [ @xmath270
, one gets that the phase @xmath271 of the optical field in the @xmath260 vacuum is locked at , using eqs .
,    [ continuityconditionsbackground2 ] @xmath272    finally , the small quantum fluctuation @xmath154 which superimposes upon the mean - field solution @xmath164 is described by eqs .
 , with @xmath33 and @xmath273 satisfying the constraints .
the relation between the fluctuations after ( @xmath260 ) and before ( @xmath259 ) the gate - shaped modulation of the optical nonlinearity along the radiation axis are obtained by matching eqs .   and ( [ fluctuationvacuumbis ] , [ lasercontribution ] ) at the @xmath214 interfaces using the continuity equation .
at the entrance face of the dielectric plate , that is , at @xmath211 [ @xmath268 , one has @xmath274\,e^{-i\beta_{0}h/2},\right.\end{aligned}\ ] ] in such a way that , taking advantage of @xmath181 and of the second of eqs .  , @xmath275.\right.\end{aligned}\ ] ] at the back face of the kerr layer , that is , at @xmath212 [ @xmath270 , the continuity of the poynting vector yields @xmath276\,e^{i\beta_{0}h/2},\right.\end{aligned}\ ] ] from which and using @xmath173 $ ] , eq .  , and eq .   one finally obtains @xmath277 where we have defined    [ newbogoliubovamplitudes ] @xmath278    with    [ newbogoliubovamplitudesbis ] @xmath279.\right.\end{aligned}\ ] ]      from the input - output relation , it is immediate to extract predictions for all the coherence properties of the transmitted light field past the slab of nonlinear medium . as already mentioned in the end of sec .
[ subsubsec : outside ] , our assumption of a perfectly coherent pump incident on the dielectric corresponds to having all the fluctuation modes of the photon field in the vacuum state @xmath264 [ @xmath266 ] . as a consequence , denoting @xmath281 the quantum average in the state @xmath264 , one initially has , just before entering the nonlinear material , @xmath282 and , thanks to the commutation relation , @xmath283      while the average light intensity remains equal to the incident one in vacuum , the effect of the quenches of the photon - photon interaction constant along the optical axis is visible in the intensity - intensity correlation function @xmath284 at a fixed point @xmath260 of the laser - beam axis , i.e. , at a fixed time @xmath285 after the quench sequence . in eq .
, @xmath286 is the ( quantum ) intensity operator at @xmath287 and the colon symbol denotes the normal - ordering operation .
since we are dealing with a strongly coherent light wave , the electric - field envelope @xmath141 at a time @xmath285 after the second quench of the photon - photon interaction constant is given by the expansion , where the weak quantum modulation @xmath154 obeys the plane - wave decomposition . at leading order in the small fluctuations @xmath154 and @xmath288
, the intensity operator @xmath289 may be expressed as    [ intensityoperatorbogoliubov ] @xmath290    is a small quantum correction to the background light intensity @xmath291 outside the nonlinear medium . taking advantage of eqs .   and , the intensity correlator becomes , at leading order in the expansion , @xmath292 physically , the second term in the right - hand side of eq .
is the shot - noise term originating from the discreteness , i.e. , the quantum nature , of the photon in vacuum ; mathematically , it is the consequence of the normal order considered in eq .  .
straightforward manipulations using eqs .
, , , , , , and lead to @xmath293\right .
\\ \label{intensitycorrelations } & \left.+\,\mathcal{v}_{\mathbf{q}}\exp\bigg[i\,\frac{\mathbf{q}_{\perp}^{2}}{2\,k_{0}}\,\big(v_{0}\,\tau-\tfrac{h}{2}\big)\bigg]\bigg|^{2}-1\bigg\}.\right.\end{aligned}\ ] ] the physics of these intensity correlations is most transparent at the exit face of the kerr layer , i.e. , at @xmath294 . in this case
, the correlation function admits a very simple integral formulation ; manipulating the @xmath295 s and the @xmath296 s given by eqs .   and denoting @xmath297 $ ] , one indeed finds    [ intensitycorrelationsexitface ] @xmath298\sin^{2}[k_{\mathrm{b}}(\mathbf{q})\,h]\right .
\\ \notag & \left.=-\frac{2}{\pi^{2}}\,\frac{2\,v_{0}}{c}\,\frac{\hbar\,\omega_{0}}{\varepsilon_{0}^{\phantom{2}}\,\mathcal{e}_{0}^{2}\,\xi_{\perp}^{2}\,\xi_{\zeta}^{\phantom{2}}}\int_{0}^{\infty}d\kappa\;\mathrm{sinc}(\kappa\,\rho)\right . \\
\label{intensitycorrelationsexitface2 } & \left.\phantom{=}\times\frac{\kappa^{2}}{\kappa^{2}+4}\,\sin^{2}\bigg(\frac{\ell}{2}\,\kappa\,\sqrt{\kappa^{2}+4}\bigg),\right.\end{aligned}\ ] ]    where @xmath195 and @xmath196 are the healing lengths in the @xmath20 plane and in the @xmath45 direction in the dielectric , with @xmath299 the internal light intensity , @xmath300 is the thickness of the nonlinear medium in units of @xmath301 , @xmath302 is the so - called sine cardinal function , and    [ relativedistance ] @xmath303    is the dimensionless relative spatiotemporal distance , @xmath304 and @xmath305 being naturally measured in units of @xmath306 and @xmath307 , respectively .    the graph of the two - body correlator @xmath308 [ eqs .  ] as a function of the dimensionless relative distance @xmath309 [ eqs .  ]
is represented in fig .
[ correlations ] .
as one can note from this plot , @xmath308 is characterized by a pronounced dip localized around @xmath310 , that is , at equal @xmath311 in the transverse plane and equal time @xmath2 : @xmath312 .
this antibunching is a direct consequence of the repulsive character ( @xmath192 ) of the photon - photon interactions in the nonlinear medium and has been widely studied in the context of many - body physics ( see , e.g. , refs .
note that the nontrivial spatial structure of @xmath308 shown in fig .
[ correlations ] makes the latter conceptually different from the single - mode antibunching behavior that is obtained , e.g. , by resonantly scattering light off a two - level emitter @xcite or by photon blockade in a strongly nonlinear cavity @xcite .
( red curve ) and @xmath313 ( black curve ) .
the oscillation pattern appearing after @xmath314 ( indicated by the red dotted line ) or @xmath315 ( black dotted line ) is the signature of the occurrence of a dynamical casimir emission of elementary excitations in the quantum fluid of light . ]
the most important feature of @xmath308 is the fringe pattern that is visible at larger @xmath309 s .
the latter is concentrated in the @xmath316 region and can be interpreted from simple physical arguments @xcite .
when the sudden change of the photon - photon interaction strength from zero to @xmath182 occurs at @xmath211 [ @xmath268 , pairs of correlated bogoliubov waves of opposite wavevectors @xmath317 are spontaneously emitted in the @xmath3 , @xmath4 , and @xmath45 directions . inside the plate
[ that is , for @xmath318 , a bogoliubov excitation of wavevector @xmath174 propagates at the group velocity ( @xmath173 $ ] ; see sec .
[ subsec : spatiallyhomogeneoussystem ] )    [ groupvelocity ] @xmath319\right .
\\ \label{groupvelocity3 } \text{and}\quad&\left.\mathbf{v}_{\zeta}(\mathbf{q})=v_{0}\,\bigg[\frac{\partial k_{\mathrm{b}}}{\partial q_{\zeta}}(\mathbf{q})+1\bigg]\,\mathbf{u}_{\zeta},\right.\end{aligned}\ ] ]    @xmath320 ( @xmath321 , @xmath322 ) being a unit vector in the @xmath3 ( @xmath4 , @xmath45 ) direction .
thus , as @xmath323 is a growing function of @xmath174 [ see the definition of @xmath187 in eq .  ]
, the modes forming a pair of entangled excitations of opposite momenta must be separated at a time @xmath324 [ and such that @xmath325 by a distance    [ transverseseparation ] @xmath326    in the @xmath20 plane and by a distance    [ temporalseparation ] @xmath327    along the @xmath44 direction , where @xmath328 are the `` sound speeds '' ( see sec .
[ subsec : spatiallyhomogeneoussystem ] ) in the transverse and temporal directions inside the nonlinear medium . as a result , according to and , the ( dimensionless ) distance @xmath309 separating two correlated counterpropagating modes has to be @xmath329 and so @xmath330 , with @xmath331 at the back face of the nonlinear medium where @xmath294 . at this point , the correlation between counterpropagating bogoliubov excitations generated by the quench of the optical nonlinearity at @xmath211 results in an oscillatory structure in the @xmath308 function that is concentrated in the @xmath330 region , which we easily verify on the examples of fig .
[ correlations ] where @xmath332 ( red curve ) and @xmath313 ( black curve ) .
the amplitude of the fringe pattern diminishes as the separation distance @xmath309 grows : the correlations between bogoliubov excitations of opposite momenta are naturally all the weaker as the excitations are separated in space and time .
this oscillatory behavior of the correlation function in position space can be interpreted as resulting from a dynamical casimir emission @xcite of elementary excitations on top of a quantum fluid of light presenting a quench of the photon - photon interaction constant in time @xmath138 .
this has been theoretically studied in ref .
@xcite in the case of a weakly interacting bose - einstein condensate presenting a quick time modulation of the atom - atom @xmath333-wave scattering length . from the experimental point of view
, the closely related phenomenon of sakharov oscillations @xcite , resulting from the interference of synchronously generated counterpropagating acoustic waves , has been recently observed in the density fluctuations of cesium superfluids after an abrupt quench in time of the atom interaction strength _ via _
feshbach resonances @xcite : this last observation is of particular interest given its connection to the temperature fluctuations in the cosmic microwave background radiation @xcite .
the emission of correlated - phonon pairs in the fluid of light due to the two quantum quenches of the photon - photon interaction constant can also be detected in momentum space by taking a far - field and spectrally resolved picture of the light emerging from the back face of the nonlinear medium . in the optical language
, this process of emission of photon pairs at wavevectors and frequencies different from the incident ones goes under the name of spontaneous four - wave mixing .
a simple measurement that could be done would be the one of the photon momentum distribution @xmath334 at a fixed point @xmath260 of the radiation axis , that is , at a fixed time @xmath285 , as a function of the wavevector @xmath174 . reminding the definition of @xmath174 , i.e. , @xmath335
, this measurement involves ( i ) an angular resolution to isolate the light deflected with a transverse wavevector @xmath235 and ( ii ) a spectral resolution to isolate the angular - frequency component of the transmitted light at @xmath336 . by construction ,
the population operator @xmath337 in the state of wavevector @xmath174 does not depend on time @xmath43 [ see eq .  ] and one simply has @xmath338 for all @xmath285 .
thus , using eqs .
, one obtains    [ momentumdistribution ] @xmath339^{2}\sin^{2}[k_{\mathrm{b}}(\mathbf{q})\,h]\right .
\\ \label{momentumdistribution3 } & \left.=a\,v_{0}\,\delta t\,\frac{\hbar\,\omega_{0}}{\varepsilon_{0}}\,\frac{4}{\kappa^{2}\,(\kappa^{2}+4)}\sin^{2}\bigg(\frac{\ell}{2}\,\kappa\,\sqrt{\kappa^{2}+4}\bigg),\right.\end{aligned}\ ] ]    where @xmath340 and @xmath341 ( @xmath342 ) empirically denote the typical transverse cross section of the beam of light and the typical duration of the measurement @xcite , respectively , @xmath300 is , as before , the normalized slab thickness , and @xmath343 is the dimensionless total excitation wavenumber , @xmath344 and @xmath184 being naturally measured in units of @xmath345 and @xmath346 , respectively .    the plot of the spontaneous - four - wave - mixing intensity distribution @xmath334 is shown in fig .  [ occupationnumber ] for @xmath332 .
the most noticeable feature of the momentum distribution is its oscillation profile originating from @xmath347 : depending on the value of the wavevector @xmath174 , the two - body emission processes at the front and the back interfaces reveal constructive or destructive interferences at well - defined periods .
the latter are determined by the finite thickness @xmath210 of the dielectric layer and the bogoliubov dispersion @xmath187 _ via _ the product @xmath348 : for example , from eq .  , the @xmath174 s for which one has destructive interferences , i.e. , zero minima in the momentum distribution , are such that    [ periodicity ] @xmath349    where @xmath350 is an integer @xmath351 .
the quantum nature of the emitted particles can be assessed _ via _ a measurement of the two - point correlation function in the far - field regime , namely , @xmath352 function which quantifies the correlations between the different angular and spectral components of the four - wave - mixing emission .
one can evaluate by means of wick s theorem [ indeed , since the bogoliubov - de gennes equation is linear , the corresponding hamiltonian  the so - called bogoliubov - de gennes hamiltonian  admits a quadratic dependence on the fluctuations @xmath353 and @xmath354 , and so , on the ladder operators @xmath355 and @xmath356 .
using eq .   which connects the outgoing excitations to the ingoing ones and eqs .   and
, we end up finding , for all @xmath285 , @xmath357.\ ] ]    .
it reveals the occurrence of constructive and destructive interferences at well - defined periods ; the @xmath358 s [ given by eqs .  ]
correspond to the @xmath174 s for which one has destructive - interference effects .
the main plot focuses on the oscillations of the momentum distribution while the inset shows the whole function . ]    as expected on general optical grounds and quantitatively predicted by eq .
, the two - body correlation signals in reciprocal space have to be concentrated along the ( diagonal ) @xmath359 and ( antidiagonal ) @xmath360 lines .
the first one reflects the vacuum fluctuations in the occupation @xmath361 while the second one is a clear signature of a complete correlation between excitations of opposite wavevectors @xmath362 .
the thermal analogue of this phenomenon was recently observed in a quasi - one - dimensional atomic bose - einstein condensate whose speed of sound was modulated _ via _ a suitable time modulation of the transverse confinement @xcite ; even though this resulted in the creation of thermal phonons of equal and opposite velocities , temperature was too high to assess the quantum origin of the observed correlations .
as a final point , it is important to highlight the novelty of our predictions with respect to standard nonlinear - optics four - wave - mixing effects , which are also described by paraxial wave equations @xcite similar @xcite to the classical counterpart of the bogoliubov - de gennes equation . as typical four - wave - mixing experiments benefit from maximizing the intensity of the four - wave - mixing signal , they are typically implemented in a phase - matched regime where both energy and momentum are conserved in the process : in our language , this requires that there exist @xmath244 pairs for which @xmath363 , which implies that the transverse ( @xmath191 ) and temporal [ @xmath159 ] masses have opposite signs . as we have seen in sec .
[ subsec : spatiallyhomogeneoussystem ] , this is often associated to ( absolute or , in some cases , only convective @xcite ) instabilities in the quantum fluid of light .
our angle on this physics is significantly different , as we are interested in the response of the photon fluid to quantum quenches . in an optical language , the energy that is deposed in the fluid by the quenches corresponds to the nonconservation of the normal ( i.e. , @xmath1 ) component of the wavevector at the @xmath214 interfaces . while the four - wave - mixing intensity that results from the quantum quenches is substantially weaker than the one expected for a perfectly phase - matched configuration , it encodes a detailed information on the conservative quantum dynamics of the fluid of light in response to the sudden jumps of the interaction parameter , which is the main object of the present work . as a specific and most illuminating example of such features ,
the evolution of the photon fluid is ruled by the bogoliubov dispersion law @xmath187 defined in eq .  , which involves the optical nonlinearity in a very nonperturbative way .
the measurement of the spatial @xmath5 component of @xmath187 , based on a bulk nonlinear medium , was recently reported in ref .
@xcite ; another experiment in a one - dimensional waveguide geometry is underway to measure the temporal @xmath44 part of the bogoliubov dispersion @xcite .
in this article , we have used a general quantum theory describing paraxial light propagation in a bulk @xmath0 medium in terms of quantum nonlinear schrdinger evolution equations  as pioneered in refs .
@xcite  to theoretically investigate the quantum statistical properties of a laser beam emerging from a finite slab of kerr material . by mapping light propagation onto a temporal evolution and by viewing the physical time as an effective spatial coordinate
, our predictions may be straightforwardly interpreted in the many - body - physics language as the response of a gas of many interacting photons to a pair of quantum quenches of the system s hamiltonian in the nonlinear interaction parameter .
this exemplifies the potential of nonlinear optics in propagating geometry as a novel platform for experimental studies of time - dependent problems in many - body physics and quantum statistical mechanics @xcite .    in the first sections of the present paper
, we have reviewed the theory underlying the space @xmath365 time mapping and we have extended it to a fully three - dimensional configuration .
we have then provided a comprehensive overview of its features and its power from a many - body - physics perspective .
in contrast to quantum fluids of light in planar - microcavity geometries that have been considered and widely studied during the last decade @xcite , the evolution of the quantum photon field in a propagating configuration is fully conservative and not subject to the unavoidable radiative ( or nonradiative ) losses inherent to cavity devices : starting from an initial many - body state that is determined by the statistical properties of the incident beam , quantum coherence can be maintained in the quantum fluid of light over macroscopic evolution times .
this gives access to a _
plethora _ of many - body dynamical quantum phenomena that can be eventually reconstructed from the statistical properties of the transmitted light using standard quantum - optics techniques .    in the second part , as a simplest , yet most remarkable example of application of the quantum formalism , we have studied in a detailed way the response of a spatially homogeneous quantum fluid of light to the pair of quenches of the interaction constant that it feels upon crossing the front and the back faces of a @xmath0 material immersed in vacuum . in the standard optical configuration where the incident light beam is coherent and the kerr medium weakly nonlinear
, we have demonstrated that the quantum quenches of the optical nonlinearity lead to the emission of pairs of correlated counterpropagating bogoliubov waves on top of the laser fluid . in particular ,
we have pointed out how the peculiar quantum features of this process can be experimentally identified in the angular and spectral distribution as well as in the two - body correlation functions of the transmitted light in both the near and the far field .    from a general standpoint ,
our results illustrate the power of well - known nonlinear - optics setups as a workhorse to investigate many - body phenomena presently of great interest in quantum statistical mechanics .
while we have so far focused on the configuration of easiest experimental implementation , this work paves the way to many other interesting directions : as a few examples , one can mention the study of the strong quantum fluctuations of the photon - field phase in reduced dimensions @xcite , of the so - called hawking radiation in analog black - hole configurations @xcite , and , on a longer run , of the strong quantum correlations that naturally appear in the presence of significant nonlinearities at the single - photon level @xcite .
we acknowledge dimitris g. angelakis , daniele faccio , stefano finazzi , pjotrs grisins , christian miniatura , and luciano vanzo for stimulating discussions .
this work was supported by the erc through the qgbe grant , by the eu - fet proactive grant aqus , project no .  640800 , and by the autonomous province of trento , partly through the siquro project ( `` on silicon chip quantum optics for quantum computing and secure communications '' ) .
j. kasprzak , m. richard , s. kundermann , a. baas , p. jeambrun , j. m. j. keeling , f. m. marchetti , m. h. szymaska , r. andr , j. l. staehli , v. savona , p. b. littlewood , b. deveaud , and le si dang , _ bose - einstein condensation of exciton polaritons _ , http://www.nature.com/nature/journal/v443/n7110/full/nature05131.html[nature ( london ) * 443 * , 409 ( 2006 ) ] .
a. amo , j. lefrre , s. pigeon , c. adrados , c. ciuti , i. carusotto , r. houdr , e. giacobino , and a. bramati , _ superfluidity of polaritons in semiconductor microcavities _ , http://www.nature.com/nphys/journal/v5/n11/full/nphys1364.html[nat .
* 5 * , 805 ( 2009 ) ] .
a. amo , s. pigeon , d. sanvitto , v. g. sala , r. hivet , i. carusotto , f. pisanello , g. lemnager , r. houdr , e. giacobino , c. ciuti , and a. bramati , _ polariton superfluids reveal quantum hydrodynamic solitons _ , http://www.sciencemag.org/content/332/6034/1167[science * 332 * , 6034 ( 2011 ) ] .
g. nardin , g. grosso , y. lger , b. pitka , f. morier - genoud , and b. deveaud - pldran , _ hydrodynamic nucleation of quantized vortex pairs in a polariton quantum fluid _ , http://www.nature.com/nphys/journal/v7/n8/full/nphys1959.html[nat . phys .
* 7 * , 635 ( 2011 ) ] .
d. sanvitto , s. pigeon , a. amo , d. ballarini , m. de giorgi , i. carusotto , r. hivet , f. pisanello , v. g. sala , p. s. s. guimaraes , r. houdr , e. giacobino , c. ciuti , a. bramati , and g. gigli , _ all - optical control of the quantum flow of a polariton condensate _ , http://www.nature.com/nphoton/journal/v5/n10/full/nphoton.2011.211.html[nat .
photonics * 5 * , 610 ( 2011 ) ] .
y. lai and h. a. haus , _ quantum theory of solitons in optical fibers .
i. time - dependent hartree approximation _ , http://journals.aps.org/pra/abstract/10.1103/physreva.40.844[phys
a * 40 * , 844 ( 1989 ) ] .
f. p. mattar and j. teichmann , _ fluid formulation of high intensity laser beam propagation using lagrangian coordinates _ , http://www.sciencedirect.com/science/article/pii/0010465581900758[comput .
* 22 * , 1 ( 1981 ) ] .
g. dekel , v. fleurov , a. soffer , and c. stucchio , _ temporal dynamics of tunneling : hydrodynamic approach _ , http://journals.aps.org/pra/abstract/10.1103/physreva.75.043617[phys .
a * 75 * , 043617 ( 2007 ) ] .
e. g. khamis , a. gammal , g. a. el , yu .
g. gladush , and a. m. kamchatnov , _ nonlinear diffraction of light beams propagating in photorefractive media with embedded reflecting wire _ , http://journals.aps.org/pra/abstract/10.1103/physreva.78.013829[phys .
a * 78 * , 013829 ( 2008 ) ] .
g. dekel , o. v. farberovich , a. soffer , and v. fleurov , _ nonlinear dynamic phenomena in macroscopic tunneling _ , http://www.sciencedirect.com/science/article/pii/s0167278908002649[physica d * 238 * , 1475 ( 2009 ) ] .
e. cohen , s. muenzel , o. v. farberovich , j. fleischer , v. fleurov , and a. soffer , _ jetlike tunneling from a trapped vortex _ , http://journals.aps.org/pra/abstract/10.1103/physreva.88.043833[phys . rev . a * 88 * , 043833 ( 2013 ) ] .      v. fleurov and r. schilling , _ regularization of fluctuations near the sonic horizon due to the quantum potential and its influence on hawking radiation
a * 85 * , 045602 ( 2012 ) ] .
s. bar - ad , r. schilling , and v. fleurov , _ nonlocality and fluctuations near the optical analog of a sonic horizon _ , http://journals.aps.org/pra/abstract/10.1103/physreva.87.013802[phys .
a * 87 * , 013802 ( 2013 ) ] .
m. vaupel , k. staliunas , and c. o. weiss , _ hydrodynamic phenomena in laser physics : modes with flow and vortices behind an obstacle in an optical channel _ , http://journals.aps.org/pra/abstract/10.1103/physreva.54.880[phys .
a * 54 * , 880 ( 1996 ) ] .
w. wan , s. muenzel , and j. w. fleischer , _ wave tunneling and hysteresis in nonlinear junctions _ , http://journals.aps.org/prl/abstract/10.1103/physrevlett.104.073903[phys .
lett . * 104 * , 073903 ( 2010 ) ] .
d. vocke , t. roger , f. marino , e. m. wright , i. carusotto , m. clerici , and d. faccio , _ experimental characterisation of nonlocal photon superfluids _ , https://www.osapublishing.org/optica/abstract.cfm?uri=optica-2-5-484[optica * 2 * , 484 ( 2015 ) ] .
t. peyronel , o. firstenberg , qi - yu liang , s. hofferberth , a. v. gorshkov , t. pohl , m. d. lukin , and v. vuleti , _ quantum nonlinear optics with single photons enabled by strongly interacting atoms _ , http://www.nature.com/nature/journal/v488/n7409/full/nature11361.html[nature ( london ) * 488 * , 57 ( 2012 ) ] .
o. firstenberg , t. peyronel , qi - yu liang , a. v. gorshkov , m. d. lukin , and v. vuleti , _ attractive photons in a quantum nonlinear medium _ , http://www.nature.com/nature/journal/v502/n7469/full/nature12512.html[nature ( london ) * 502 * , 71 ( 2013 ) ] .
e. m. wright , _ quantum theory of soliton propagation in an optical fiber using the hartree approximation _ , http://journals.aps.org/pra/abstract/10.1103/physreva.43.3836[phys .
a * 43 * , 3836 ( 1991 ) ] .
b. crosignani , p. di porto , and a. treppiedi , _ finite - dimensional hamiltonian approach to the quantum theory of soliton propagation _ , http://iopscience.iop.org/1355-5111/7/2/001/[quantum semiclassic
* 7 * , 73 ( 1995 ) ] .
a. b. matsko and v. v. kozlov , _ second - quantized models for optical solitons in nonlinear fibers : equal - time versus equal - space commutation relations _ , http://journals.aps.org/pra/abstract/10.1103/physreva.62.033811[phys .
a * 62 * , 033811 ( 2000 ) ] .
a. polkovnikov , k. sengupta , a. silva , and m. vengalattore , _ colloquium : nonequilibrium dynamics of closed interacting quantum systems _ , http://journals.aps.org/rmp/abstract/10.1103/revmodphys.83.863[rev .
* 83 * , 863 ( 2011 ) ] .
m. elazar , v. fleurov , and s. bar - ad , _ all - optical event horizon in an optical analog of a laval nozzle _ , http://journals.aps.org/pra/abstract/10.1103/physreva.86.063821[phys .
a * 86 * , 063821 ( 2012 ) ] .
m. elazar , s. bar - ad , v. fleurov , and r. schilling , _ an all - optical event horizon in an optical analogue of a laval nozzle _ , http://link.springer.com/chapter/10.1007/978-3-319-00266-8_12[lect .
notes phys .
* 870 * , 275 ( 2013 ) ] .                our choice of referring to @xmath366 as a hamiltonian follows from its interpretation as an evolution operator
however , as it generates propagation in the @xmath1 direction , some reader may object that it is related to the @xmath1 component of the linear momentum of the underlying field theory . for the sake of transparency , throughout this paper , we shall stick to the first wording .        in dirac
s approach @xcite , the fields @xmath33 and @xmath68 are canonically conjugated ( in the sense of the same-@xmath43 poisson bracket @xmath90 ) to @xmath367 and @xmath368 , respectively . as a result , the theory is characterized by two primary constraints , @xmath369 and @xmath370 , which may be shown to be of second class .
therefore , @xmath90 becomes irrelevant and must be replaced in the evolution equations by the equal-@xmath43 dirac bracket @xmath371 .
making use of the explicit definition of @xmath371 in terms of @xmath90 @xcite , straightforward calculations yield @xmath372 and @xmath373 , which gives the quantum equal-@xmath43 commutation relations after the replacement of @xmath33 with @xmath374 , @xmath68 with @xmath375 , and @xmath371 with @xmath135/(i\,\hbar)$ ] .
the results deduced from dirac s theory hence exactly coincide with the ones deduced from the canonical quantization method , legitimating our approach .
r. ganapathy , b. a. malomed , and k. porsezian , _ modulational instability and generation of pulse trains in asymmetric dual - core nonlinear optical fibers _ , http://www.sciencedirect.com/science/article/pii/s0375960106001940[phys .
a * 354 * , 366 ( 2006 ) ] .
l. berg , c. goudard , j. schjdt - eriksen , and h. ward , _ filamentation patterns in kerr media vs. beam shape robustness , nonlinear saturation and polarization states _ , http://www.sciencedirect.com/science/article/pii/s0167278902007406[physica d * 176 * , 181 ( 2003 ) ] .
i. carusotto , s. fagnocchi , a. recati , r. balbinot , and a. fabbri , _ numerical observation of hawking radiation from acoustic black holes in atomic bose - einstein condensates _ , http://iopscience.iop.org/1367-2630/10/10/103001[new j. phys . * 10 * , 103001 ( 2008 ) ] .
p. deuar , a. g. sykes , d. m. gangardt , m. j. davis , p. d. drummond , and k. v. kheruntsyan , _ nonlocal pair correlations in the one - dimensional bose gas at finite temperature _ , http://journals.aps.org/pra/abstract/10.1103/physreva.79.043619[phys .
a * 79 * , 043619 ( 2009 ) ] .
t. jacqmin , j. armijo , t. berrada , k. v. kheruntsyan , and i. bouchoule , _ sub - poissonian fluctuations in a 1d bose gas : from the quantum quasicondensate to the strongly interacting regime _ , http://journals.aps.org/prl/abstract/10.1103/physrevlett.106.230405[phys .
lett . * 106 * , 230405 ( 2011 ) ] .
larr , a. recati , i. carusotto , and n. pavloff , _ quantum fluctuations around black hole horizons in bose - einstein condensates _ , http://journals.aps.org/pra/abstract/10.1103/physreva.85.013621[phys
a * 85 * , 013621 ( 2012 ) ] .
p. calabrese and j. cardy , _ time dependence of correlation functions following a quantum quench _ , http://journals.aps.org/prl/abstract/10.1103/physrevlett.96.136801[phys .
. lett . * 96 * , 136801 ( 2006 ) ] .
g. t. moore , _ quantum theory of the electromagnetic field in a variable - length one - dimensional cavity _ , http://scitation.aip.org/content/aip/journal/jmp/11/9/10.1063/1.1665432[j .
phys . * 11 * , 2679 ( 1970 ) ] .
s. a. fulling and p. c. w. davies , _ radiation from a moving mirror in two dimensional space - time : conformal anomaly _ , http://rspa.royalsocietypublishing.org/content/348/1654/393[proc .
london , ser .
a * 348 * , 393 ( 1976 ) ] .
i. carusotto , r. balbinot , a. fabbri , and a. recati , _ density correlations and analog dynamical casimir emission of bogoliubov phonons in modulated atomic bose - einstein condensates _ , http://link.springer.com/article/10.1140%2fepjd%2fe2009-00314-3[eur .
j. d * 56 * , 391 ( 2010 ) ] .
a. d. sakharov , _ the initial stage of an expanding universe and the appearance of a nonuniform distribution of matter _ , zh . eksp .
teor . fiz . *
49 * , 345 ( 1965 ) ; http://www.jetp.ac.ru/cgi-bin/e/index/e/22/1/p241?a=list[sov .
jetp lett .
* 22 * , 241 ( 1966 ) ] .
. hung , v. gurarie , and c. chin , _ from cosmology to cold atoms : observation of sakharov oscillations in a quenched atomic superfluid _ , http://www.sciencemag.org/content/341/6151/1213[science * 341 * , 1213 ( 2013 ) ] .      to treat the infinite factor @xmath376 appearing in eq .
, one replaces the @xmath377 function in momentum space with @xmath378 , where @xmath379 denotes a large cylindrical volume in the @xmath380 space whose cross section @xmath340 in the @xmath5 plane and length @xmath381 in the @xmath2 direction have the natural interpretations given in the main text .
the same reasoning holds for eq .  .
jaskula , g. b. partridge , m. bonneau , r. lopes , j. ruaudel , d. boiron , and c. i. westbrook , _ acoustic analog to the dynamical casimir effect in a bose - einstein condensate _ , http://journals.aps.org/prl/abstract/10.1103/physrevlett.109.220401[phys .
lett . * 109 * , 220401 ( 2012 ) ] .
q. lin , o. j. painter , and g. p. agrawal , _ nonlinear optical phenomena in silicon waveguides : modeling and applications _ , https://www.osapublishing.org/oe/abstract.cfm?uri=oe-15-25-16604[opt .
exp . * 15 * , 16604 ( 2007 ) ] .
x. zhang , s. gao , and s. he ,
_ optimal design of a silicon - on - insulator nanowire waveguide for broadband wavelength conversion _ , http://www.jpier.org/pier/pier.php?paper=08120601[prog .
89 * , 183 ( 2009 ) ] .
s. gao , z. li , e .- k .
tien , s. he , and o. boyraz , _ performance evaluation of nondegenerate wavelength conversion in a silicon nanowire waveguide _ , http://www.osapublishing.org/jlt/abstract.cfm?uri=jlt-28-21-3079[j .
lightwave technol .
* 28 * , 3079 ( 2010 ) ] .
m. blaauboer , a. g. kofman , a. e. kozhekin , g. kurizki , d. lenstra , and a. lodder , _ superluminal optical phase conjugation : pulse reshaping and instability _ , http://journals.aps.org/pra/abstract/10.1103/physreva.57.4905[phys .
a * 57 * , 4905 ( 1998 ) ] .
g. izs , m. santagiustina , m. s. miguel , and p. colet , _ pattern formation in the presence of walk - off for a type ii optical parametric oscillator _ , https://www.osapublishing.org/josab/abstract.cfm?uri=josab-16-9-1592[j .
b * 9 * , 1592 ( 1999 ) ] .
m. blaauboer , d. lenstra , and a. lodder , _ optical josephson effect using phase - conjugation mirrors : an analogy with superconductors _ , http://www.nit.eu / archive?itemid=71&option = com_czasopisma&view = wynikiszukania&jezyk = en&nr_mies=&rok=&autor = blaauboer&tytul= & opis=&slowakl=&limit=50[jtit * 1**2 * , 24 ( 2000 ) ] .