Light propagation in inhomogeneous media, coupled quantum harmonic oscillators and phase transitions

This contribution has two main purposes. First, using classical optics we show how to model two coupled quantum harmonic oscillators and two interacting quantized fields. Second, we present classical analogs of coupled harmonic oscillators that correspond to anisotropic quadratic graded indexed media in a rotated reference frame, and we use operator techniques, common to quantum mechanics, to solve the propagation of light through a particular type of graded indexed medium. Additionally, we show that the system presents phase transitions.

The existence of analogies between quantum and classical mechanics has been applied for many years, particularly in the generation of mathematical tools to provide solutions of optical problems and vice versa [1][2][3][4][5][6][7][8][9][10][11] . The reason is that the optical paraxial wave equation is mathematically equivalent to the stationary Schrödinger equation and, on the other hand, Helmholtz equation has a treatment similar to that of time-independent Schrödinger equation. Exploring analogies between different physical problems is useful, as it allows exporting insights from one context to the other.
Classical optics systems have been used to model quantum optical phenomena 6 . There has been also a proposal by Man'ko et al. 12 to realize quantum computation by using quantum-like systems. Coherent random walks have been shown to occur in free propagation provided the initial wave function is tailored properly 13 . The modelling of time-dependent harmonic oscillators has been achieved in a graded indexed (GRIN) media 14 allowing to show the splitting of modes in second order solutions of the Helmholtz equation 15 . The propagation of classical fields in a graded indexed medium with linear dependence has shown to produce and control Airy modes 16 .
On the other hand, much attention has been recently devoted to the study of phase transitions 17,18 in the Rabi Model 19,20 . By means of the Holstein-Primakoff transformation 21 , it may be argued that the interaction of two quantized fields is similar to the Rabi interaction, when we approximate the spin-flip operators as creation and annihilation operators and the energy Pauli matrix as a harmonic oscillator, and thus the two interacting fields can show phase transitions.
Being aware of all of these correspondences, in this contribution we apply methods commonly used in quantum optics to solve a problem of classical optics. Although, we could use variable separation to solve the propagation of light in GRIN media, we believe it is of interest to use operator techniques, commonly used in quantum mechanics, as there may be cases in which the refraction index not only depends on x and y but also on the propagation distance z (see for instance Urzúa et al. 22 ), and then the separation of variables is not possible.
In particular, we show in the next section that classical light propagation in an inhomogeneous medium models a coupled quantum harmonic oscillators or the interaction of two quantized fields. In the same section, we also show the conditions for a phase transition of the coupled system, in the sense that the set of parameters involved lead one of the oscillators go through harmonic to free particle and finally to a repulsive oscillator. In the Section Invariant modes, we show that it is possible to obtain invariant modes in the classical propagation of light in the specific inhomogeneous medium considered. In the Section titled Helmholtz equation, we briefly discuss the possibility of splitting a beam by taking into account the fact that, to second order in the Helmholtz equation, a kind of quantum Kerr medium is modelled. Finally, the last section is left for conclusions. www.nature.com/scientificreports www.nature.com/scientificreports/

paraxial Wave equation for inhomogeneous Media
There exist classical optical systems related to waveguiding of optical modes that may be described by the paraxial wave equation where λ and k 0 = 2πn 0 /λ are the wavelength and the wavenumber of the propagating mode, respectively, and n 0 is the homogeneous refraction index. The function k 2 (x, y) describes the inhomogeneity of the medium responsible for the waveguiding of the optical field E. The inhomogeneity may be physical, in the sense that it is produced in the medium when fabricated; an example of such a medium is a graded index fiber. In particular, we will consider 23 x y 2 0 2 2 2 being k x , k y and g parameters that characterize the inhomogeneous medium. In ref. 23 , it was proposed a way, that could be easy to realize experimentally, of generating the inhomogeneity function k 2 (x, y) by the co-propagation of two modes in a Kerr nonlinear medium.
Equation (1), which is analogous in structure to the Schrödinger equation 24 , together with the inhomogeneity (2), has been solved for the special case of  g k k , x y using the so-called rotating wave approximation 23 ; here, we solve it for any set of parameters. We may rewrite Eq. (1) as we make the transformation  = θ † E R , such that we obtain a Schrödinger-like equation for  as x x y y 0 x y were obtained from the set of transformations www.nature.com/scientificreports www.nature.com/scientificreports/ and then Equations (7) take the form so that (6) is transformed into the equation of two uncoupled quantum harmonic oscillators for which we know the formal solution.
We may see from Eq. (10a), that for some values of the parameters k x , k y and g of the inhomogeneity of the medium, the effective frequency of the oscillator,  k x , goes from positive to negative (see Fig. 1), therefore undergoing a phase transition 17,18 . It basically shows that the oscillator in the variable x goes from a harmonic oscillator to a repulsive harmonic oscillator, passing through a free particle 25 . In other words, if we define a critical value for g, If in Eq. (6), we write x and y in terms of annihilation and creation operators, we see that the so-called counter rotating terms have not been neglected (i.e., the rotating wave approximation was not performed). The presence of those terms is responsible for the phase transition in the propagation of light in the GRIN (Gradient-index) medium we are considering 23 .
A mechanical system equivalent to the light propagation analyzed above is given by two masses connected between them by springs and connected to physical or mathematical constrains (like a pair of walls) also by springs, as it is depicted in Fig. (2). The equations that rule this system are 26 www.nature.com/scientificreports www.nature.com/scientificreports/  in the Hamilton equations, performing some derivatives and some algebra, we get the system (14). So, (15) is indeed a Hamiltonian of the system illustrated in Fig. (2). We can make now the connection between our model of light propagation with this mechanical system. If in Eq. (3) we make z → −t and

invariant Modes
It is easy to show that there exist invariant modes (without dependence on the propagation) for the inhomogeneity we are studying, provided that the parameter  k x is positive, so that we keep ourselves in the regime of an harmonic oscillator. Consider a field at z = 0 given by n n x y where ψ nx (x) and ψ ny (y) are Hermite-Gauss functions; i.e., eigenfunctions of the harmonic oscillators given in (11). This allows us to write the initial transformed field (16) as x y  and therefore write the solution of Eq. (11) as x y x x y y n n 0 0 0 x y  that produces the propagated field Some examples of the behavior of the propagated field E(x, y, z) in (18) are shown in Fig. 3. It is convenient to stress that although the propagation is invariant, there are some values of the coefficients {k x , k y , g} for which there is a phase transition and therefore the parameter  k x goes from positive (an usual harmonic oscillator) to negative (a repulsive harmonic oscillator). In the negative region, see Fig. 1, the Hermite-Gauss functions are not anymore eigenfunctions of the repulsive harmonic oscillator, therefore an invariant mode is not produced. As it was pointed out by Yuce 27 and Muñoz 28 , the eigenfunctions of the repulsive harmonic oscillator are a combination of the eigenfunctions of the free particle due to the existence of a equivalence between the repulsive potential and www.nature.com/scientificreports www.nature.com/scientificreports/ the free evolution through a proper canonical transform. Since these eigenfunctions are a combination of plane waves and waves confined in a box or propagating to a potential barrier, the non invariance argument above when g > g c follows.

Helmholtz equation
We consider now the complete Helmholtz equation For the case we considered above, in the section titled Paraxial wave equation for inhomogeneous media, and doing the same transformation (5), we arrive to x x y y 0 0 0 that may be developed to second order in κ 2 = k 0 2 − k x − k y to obtain the approximation  www.nature.com/scientificreports www.nature.com/scientificreports/

conclusions
We have shown that we can model the interaction of two masses through springs, i.e. two coupled quantum harmonic oscillators and the interaction of two quantized fields in a graded index media. The solutions we have given to second order approximation of the Helmholtz equation show how a beam splitter may be achieved in such media. We have also studied that phase transitions occur in the propagation of light in this medium, by showing that one of the harmonic oscillators is inverted for certain values of the parameters involved. factorization of the operator R̂θ. In this appendix the unitary operator y x is factorized. We define the operators The exponential of the operator K 0 is nothing but a product of squeeze operators [31][32][33][34][35] . The above operators have the commutators 0 and thus generate a su(1, 1) algebra. We find the factorization where θ R is given in (32) and F(x, y) is an arbitrary, but well behaved, function of x and y. To study the action of the θ R operator over an arbitrary function F(x, y), we make Note that the operator Ŝ xy is a product of squeeze operators [31][32][33][34][35] in x and y. We can prove that