Phi-Psi model for Electrodynamics in dielectric media: exact quantisation in the Heisenberg representation

We investigate the quantization in the Heisenberg representation of a model which represents a simplification of the Hopfield model for dielectric media, where the electromagnetic field is replaced by a scalar field $\phi$ and the role of the polarization field is played by a further scalar field $\psi$. The model, which is quadratic in the fields, is still characterized by a nontrivial physical content, as the physical particles correspond to the polaritons of the standard Hopfield model of condensed matter physics. Causality is also taken into account and a discussion of the standard interaction representation is also considered.


Introduction
In recent years, the investigations on possible revelations of the Hawking effect in analogues realised in dielectric media, [1][2][3][4][5][6][7][8][9], have raised the necessity of disposing of a model describing the quantum electromagnetic field interacting with a dispersive medium, reproducing the typical phenomenological dispersion relations.
We recall that in the Hopfield model, a purely phenomenological quantisation of the electromagnetic field in the dielectric medium is replaced by a picture where the electromagnetic field interacts with a set of oscillators reproducing sources for dispersive properties of the electromagnetic field in matter [10][11][12][13]. We stress that we do not take into account absorption in our paper, which is reasonable as far as the phenomena one is interested in are not too close to the absorption a e-mail: francesco.belgiorno@polimi.it b e-mail: sergio.cacciatori@uninsubria.it c e-mail: f.dallapiazza@gmail.com d e-mail: m.doronzo@uninsubria.it region. Including absorption would imply a much more tricky approach (cf. e.g. [14][15][16]), which is far beyond the scope of the present paper.
Our model satisfies the requirement to be fully relativistic covariant. This is in agreement with the necessity to proper simulate electrodynamics of moving media [17][18][19][20][21][22][23][24], where phenomenological electrodynamics is adopted. Also in our case, in order to get a complete analysis, one needs to change the inertial frame passing, for example, from the frame where the medium is at rest, to the frame where a given signal is at rest, or to the lab frame if it does not coincide with one of them. With this in mind, a set of models taking into account the dispersion relations have been developed, based on a covariant reformulation of the Hopfield model [25][26][27][28]. In particular, in [25] the Hopfield model has been presented in a simplified version, where, in a two dimensional model, the electromagnetic field has been replaced by a massless scalar field, linearly coupled to a polarisation field, represented by a field of oscillators with characteristic pulsation ω 0 . This had the aim of simplifying several technical complications, but keeping the main relevant characteristics, as Lorentz covariance and the right dispersion relation. In the latter reference, we were interested in the Hawking effect and we did not focus on a systematic quantisation of the model. A perturbative quantisation in a given gauge and in the lab frame of the whole homogeneous and isotropic relativistic Hopfield model has been presented in [27], and a Lorentz and gauge covariant perturbative quantisation has been provided in [28]. However, this model can be quantised exactly, the exact quantisation is involved with some significant physical characteristics of the given system, and in particular in its spectral properties, so quantizing non-perturbatively in the Heisenberg representation of quantum field theory is far more than a simple and straightforward exercise. The exact quantisation of the relativistic Hopfield model will be presented elsewhere [29]: its construction is quite involved and passes through several technical intricacies, which go beyond the ones characterising the present model.
Since the main steps characteristic of the model are present also in the scalar simplified analogue, without the intricacies due to gauge invariance and the presence of unphysical modes, we will present here the exact quantisation of the scalar Hopfield model, which we dub the − model. Since our aim is to illustrate here the strategy of [29] without hiding it behind technical details, we will mainly present the results without proves, illustrating the main steps at an intuitive level, by mean of precise statements, which will be fully proved, in a more general form, in [29] and [30].
In Sect. 2 we present the model and its quantisation, by showing how the Fock representation for the whole interacting model can be realised. In Sect. 3 we show how the problem of causality is related to Lorentz covariance and show how it can be realised in our model, which is only covariant and not invariant under Lorentz transformations. In Sect. 4 we compute the propagator both directly from the Fock representation and with the path integral method, and state their equivalence. In Sect. 5 we discuss the interaction representation and the Fano diagonalisation method. In Sect. 6 we present some further discussion. All statements are simple consequences of the ones proved in [29] and [30].

The − model and its quantisation
We consider the D + 1 dimensional − model whose classical dynamics is defined by the action where v v v is the spacetime velocity of the rest frame for the ψ field. The conjugate momenta are so that the Hamiltonian is The classical equations of motion in the Fourier space are where ω := k μ v μ . The dispersion relation is given by This defines the support of the solutions in the momentum space, with two positive branches, corresponding to the two solutions having positive ω, which we will indicate with k 0 (a) , a = 1, 2. Everywhere the suffix (a) will mean "evaluated at k 0 = k 0 (a) (k)". We will also use the symbol k for the spatial component of a spacetime vectork k k, and similar for all vectors, whereas k 2 := k k k · k k k.
The classical solutions of the equations of motion are where and is endowed with the conserved scalar product A basis of positive norm plane waves is Notice that which gives and then and all other commutators vanish. This result has been obtained by imposing the requirement that the quantum fields satisfy the equal time canonical commutation relations, which, as a consequence, are satisfied. This way, one can proceed in the usual way in constructing the Fock space, starting from the vacuum state , the unique normalised state that is annihilated by all a (a) (k), we can realise the Fock space as the completion of the set of states generated by all polynomial actions of the creator fields a † (a) (k). This is standard and free of particular difficulties, apart from the fact that the vacuum state is unique only after fixing a choice of v v v, since the theory is not invariant under the whole Poincaré group, but only under the subgroup leaving v v v invariant.
Here, we simply notice that the fact that the CCR are satisfied, together with Lorentz covariance, allows one to prove that the principle of causality is satisfied. Since the explicit presence of v v v breaks the Lorentz invariance, the question of the covariance is a little bit delicate and requires a careful analysis.

Causality and covariance
The algebra of quantum fields is generated by the canonical commutation relations (CCR): where we indicated only the non-zero contributions. Causality conditions are apparently a little bit stronger: for any pair of points x x x, y y y spatially separated, (x x x − y y y) 2 < 0. If the theory is Lorentz invariant, then (20) follow from the CCRs, since we can change frame into the one where x 0 = y 0 = t and then employ the CCRs in order to prove their vanishing. Our model Lagrangian is associated with equations of motion which are covariant with respect to the Lorentz group. This fact does not correspond to a full Lorentz invariance, due to the fact that the Lorentzian metric is not the only absolute object of the theory, but a further absolute object [31] appears: the velocity v v v of the medium. This implies that our theory, and any covariant theory of a dielectric medium (cf. e.g. [32]), is involved with a preferred frame, which corresponds to the rest frame of the medium. The explicit presence of the vector v v v in the Lagrangian implies that, in general, boosts are no more symmetries. This breaks the Poincaré symmetry group down to the subgroup leaving v v v invariant. This behaviour should not be a surprise, as it is common to all the cases where e.g. Klein-Gordon equation is studied in the presence of an external potential. Loss of Lorentz invariance is evident, but, at the same time, the field equations are covariant, and solutions are transformed into solutions of the Klein-Gordon equation by the Lorentz group, provided that external potential is transformed too. See e.g. [33] (p. 516). In this sense, also our theory remains Lorentz covariant, as the equations of motions are, and the above argument proving causality applies again if the covariance is respected at the level of the representation of the quantum theory. Covariance ensures that different inertial observers perceive the same physics, i.e. are involved with the same processes with the same probability. In particular, the number of polaritons in the process remain the same. Unitary maps between Fock spaces of different inertial observer are a natural consequence of Poincaré covariance, and is a consequence of the fact that the Lagrangian is invariant under simultaneous transformations of the fields and of the vector v v v under Poincaré group. So, at the quantum level there exists a Fock space F v v v for any v v v and a set of unitary maps where U is the set of all possible isometric maps among Fock spaces F v v v and F w w w , and U ( ) is defined by and where o is intended to be any one among the operators a (a) (k), a † (a) (k), a = 1, 2.

The propagator
Since the theory is Gaussian, it is completely determined by the two point functions. It is given by the matrix We will also write Since G I J v v v (x x x, y y y) − is easily obtained from G I J v v v (x x x, y y y) + , we will write down only the latter: The propagator can be determined also by means of the path integral formulation. After introducing the currents J φ and J ψ , we can define the functional generating the propagators: From this we can formally compute the propagator, which as a result is found to be (32) where M M M is defined in (5), so that Naturally, this is not the complete story, since this expression requires a prescription avoiding the poles defined by the dispersion relation. Such a prescription must respect causality. As a result this can be accomplished by means of a iε Feynman prescription. Indeed, we have the following.

Theorem 1 The propagator is
where M M M −1 iε (k k k) is obtained from (33) by taking the complex shifts k 2 → k 2 + iε, This proposition is a particular case of a more general one proved in [30].

The interaction representation and Fano diagonalisation
In this section we take into account the more standard Interaction representation, and perform the so-called Fano diagonalisation [11] of the full Hamiltonian operator in order to find its eigenmodes. As the Hamiltonian is quadratic in the fields and their conjugate momenta, we are able to obtain an exact result which leads again to polaritons as physical states of the system. It is interesting to stress that this approach, which is pursued both in the original paper by Hopfield [10] and in standard textbooks (see e.g. [12,13]), is in principle apt to perturbation theory and leads to the same result thanks to the diagonalisation process.
In line of principle, the interaction representation is constructed by assuming that g is small and allows a perturbation theory in powers of g. For simplicity, we consider only the case where v v v = (c, 0), and we put c = 1. The Hamiltonian is then characterised by three contributions: two free-field From (49) one obtains the equation which amounts to which again gives the same eigenmodes as in the previous sections. In particular, two positive branches E ± can be obtained, with associated eigenvectors. One can obtain α a , with a = ±, which correspond to the a (a) given in the previous analysis.

Final comments
In the appendix it is shown that D R (a) , a = 1, 2, are always positive, apart in k = 0 for the lower branch. In this case, D R vanishes linearly in |k|, and the integrals defining the fields (in a distributional sense) are well posed for D > 1. Thus, the fields and propagators are well defined as tempered distributions.
On the opposite, in a two dimensional spacetime, the integrals diverge unless the oscillator modes vanish quickly enough at the origin. This divergence at k = 0 cannot be interpreted as an infrared divergence, since it does not occur only in the propagator, but also in the definition of the fields, so that it needs to be eliminated by a suitable choice of the space of test functions. Assuming that these must be chosen inside the set of smooth rapidly decreasing functions S(R 2 ), we must consider functions whose Fourier transform vanishes in k = 0 when k 0 is evaluated on k 0 (2) (k). Since also k 0 vanishes, by employing smoothness, it is sufficient to consider functions in S(R 2 ) whose Fourier transform vanish in k k k = 0 0 0. These are the rapidly decreasing smooth functions having null mean. So, the test functions must be chosen in Thus, it is exactly the same as for any massless free field. After having defined the fields, we have been able to treat them exactly, being the action quadratic. However, despite the theory is Gaussian and, then, essentially a free theory, in a sense the interaction manifests itself through a highly non-trivial dispersion relation.
The second difficulty, related to causality, is covariance, which is realised in a non-trivial way. We have just sketched how this can be done: in practice the quantum algebra is represented on a bundle of Fock spaces over the homogeneous where P is the Poincaré group and G its subgroup that leaves v v v invariants (the little group of v v v). Finally, we have computed the two point function, which characterises the whole theory, being Gaussian. We have done it both starting from the canonical representation and with the path integral method. Here, the iε Feynman-Stückelberg prescription has been introduced and stated to be equivalent to the causal propagator computed in the oscillator representation. The proof can be found in [30].
In the relativistic Hopfield model, the true target of all our efforts, all these difficulties are present and amplified by the presence of a higher number of field components, including non-physical ones, a larger number of spectral branches, the