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.

Introduction

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-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.

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 k0 = 2πn0/λ are the wavelength and the wavenumber of the propagating mode, respectively, and n0 is the homogeneous refraction index. The function k2(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]

being k, k 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 k2(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 using the so-called rotating wave approximation[23]; here, we solve it for any set of parameters. We may rewrite Eq. (1) aswithwhere the momentum operators have the usual definition on configuration space, and . By considering the unitary operatorwe make the transformation , such that we obtain a Schrödinger-like equation for aswherewere obtained from the set of transformations

From Eq. (7c), we choose and we haveand then Equations (7) take the formso that (6) is transformed into the equation of two uncoupled quantum harmonic oscillatorsfor which we know the formal solution.

We may see from Eq. (10a), that for some values of the parameters k, k and g of the inhomogeneity of the medium, the effective frequency of the oscillator, , 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,we have the following three cases:

Figure 1

Plot surface of  as a function of k and g, with k = 5. We see that there is a phase transition from positive to negative values of  when g goes to bigger values respect to the former.

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]where q0, q1 are the positions of mass m0 and mass m1, respectively, and κ0, κ1, κ2 are the three spring constants. Introducing the Hamiltonian

Figure 2

Two masses coupled by springs, which is described by the Hamiltonian (15). This system exhibits the same features found in the Hamiltonian (4) for the light propagation in an inhomogeneous medium.

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 andthe two Hamiltonians, (4) and (15), are identified. We can go also in the other sense, from the mechanical system to the light propagation model, making t → −z 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 is positive, so that we keep ourselves in the regime of an harmonic oscillator. Consider a field at z = 0 given bywhere ψ(x) and ψ(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 and therefore write the solution of Eq. (11) asthat produces the propagated fieldthat is clearly an invariant mode with associated phases depending on the angular quantities and , and in the independent energy levels n and n. The way the operator acts on the Hermite-Gauss functions on (18) is detailed in the Appendices.

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, k, g} for which there is a phase transition and therefore the parameter 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 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 follows.

Figure 3

Plot of the square modulus of (18) for different values of n and n. Upper left: n = n = 1, upper right: n = n = 2, lower left: n = n = 3, lower right: n = n = 4. For all the cases k = 1.2, k = 1.5 and g = 0.25. The propagation is invariant along the propagation axis, i.e., it does not depend on the propagation distance z, despite the selection of the energy levels.

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 tothat may be developed to second order in κ2 = k02 − k − k to obtain the approximationwhere

This equation allows the splitting of fields[29] as the square terms resemble nonlinear quantum optical interactions, namely a quantum Kerr medium[30].

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 operatoris factorized.

We define the operators

The exponential of the operator is nothing but a product of squeeze operators[31-35]. The above operators have the commutators

and thus generate a su(1, 1) algebra. We find the factorization

where f1, f2, f3 are given byor, writing it explicitly in terms of the original operators,

Action of the operator R̂ over a function F(x, y)

It is not difficult to show that

where 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 operator over an arbitrary function F(x, y), we makewhere

Note that the operator is a product of squeeze operators[31-35] in x and y. We can prove that

and the action of the squeeze operators on the variables x and y are