Hopfield dielectric

Hopfield dielectric – in quantum mechanics, a model of dielectric consisting of quantum harmonic oscillators interacting with the modes of the quantum electromagnetic field. The collective interaction of the charge polarization modes with the vacuum excitations, photons leads to the perturbation of both the linear dispersion relation of photons and constant dispersion of charge waves  by the avoided crossing between the two dispersion lines of polaritons. Similar to the acoustic and the optical phonons and far from the resonance one branch is photon-like while the other charge is wave-like. Mathematically the Hopfield dielectric for the one mode of excitation is equivalent to the Trojan wave packet in the harmonic approximation. The Hopfield model of the dielectric predicts the existence of eternal trapped frozen photons similar to the Hawking radiation inside the matter with the density proportional to the strength of the matter-field coupling.

Theory
The Hamiltonian of the quantized Lorentz dielectric consisting of $$N$$ harmonic oscillators interacting with the quantum electromagnetic field can be written in the dipole approximation as:


 * $$H=\sum\limits_{A=1}^N{{p_A}^2\over 2m}+{{m{\omega}^2}\over2}{x_A}^2-

e{x_A}\cdot E(r_A) +\sum\limits_{\lambda=1}^2\int d^3ka_{\lambda k}^+ a_{\lambda k}\hbar c k$$ where
 * $$E(r_A)={i\over L^3}\sum\limits_{\lambda=1}^2\int d^3k [{{c k}\over {2\epsilon_0}}]^{1\over 2}

[e_\lambda(k)a_\lambda(k)\exp(ikr_A)-H.C.]$$ is the electric field operator acting at the position $$r_A$$.

Expressing it in terms of the creation and annihilation operators for the harmonic oscillators we get


 * $$H=\sum\limits_{A=1}^N(a_A^+\cdot a_A)\hbar \omega

-{e\over{{\sqrt{2}}\beta}}(a_A+{a_A}^+)\cdot E(r_A)+\sum_\lambda\sum_k a_{\lambda k}^+a_{\lambda k}\hbar c k$$

Assuming oscillators to be on some kind of the regular solid lattice and applying the polaritonic Fourier transform
 * $$B_k^+={1\over{\sqrt{N}}}\sum\limits_{A=1}^N\exp(ikr_A)a_A^+,$$
 * $$B_k={1\over{\sqrt{N}}}\sum\limits_{A=1}^N\exp(-ikr_A)a_A$$

and defining projections of oscillator charge waves onto the electromagnetic field polarization directions
 * $$B_{\lambda k}^+=e_{\lambda}(k)\cdot B_k^+$$
 * $$B_{\lambda k}=e_{\lambda}(k)\cdot B_k,$$

after dropping the longitudinal contributions not interacting with the electromagnetic field one may obtain the Hopfield Hamiltonian
 * $$H=\sum_{\lambda}\sum_k(B_{\lambda k}^+B_{\lambda

k}+{1\over2})\hbar \omega +\hbar cka_{\lambda k}^+a_{\lambda k} +{ie\hbar\over {\sqrt {\epsilon_0 m\omega}}}\sqrt{N\over V} {\sqrt {ck}} [B_{\lambda k}a_{\lambda -k} +B_{\lambda k}^+a_{\lambda k}-B_{\lambda k}^+a_{\lambda-k}^+ -B_{\lambda k}a_{\lambda k}^+]$$ Because the interaction is not mixing polarizations this can be transformed to the normal form with the eigen-frequencies of two polaritonic branches:
 * $$H=\sum_{\lambda}\sum_k \left[ \Omega_{+}(k)C_{\lambda + k}^+C_{\lambda+ k}+\Omega_{-}(k)C_{\lambda - k}^+C_{\lambda - k}\right]+const$$

with the eigenvalue equation
 * $$[C_{\lambda \pm k},H]=\Omega_{\pm} (k) C_{\lambda \pm k}$$
 * $$C_{\lambda \pm k}=c_{1} a_{\lambda k} + c_{2} a_{\lambda -k} + c_{3} a_{\lambda k}^+ + c_{4} a_{\lambda - k}^++c_{5} B_{\lambda k} + c_{6} B_{\lambda -k} + c_{7} B_{\lambda k}^+ + c_{8} B_{\lambda -k}^+$$

where
 * $$\Omega_{-}(k)^2={   {\omega^2+\Omega^2-\sqrt

{ {(\omega^2-\Omega^2)} ^2+4{g}\omega^2\Omega^2 }\over 2   }},$$ with
 * $$\Omega_{+}(k)^2=$$,
 * $$\Omega(k)=ck,$$

(vacuum photon dispersion) and
 * $$g={{Ne^2}\over{Vm\epsilon_0\omega^2}}$$

is the dimensionless coupling constant proportional to the density $$N/V$$ of the dielectric with the Lorentz frequency $$\omega$$ (tight-binding charge wave dispersion). One may notice that unlike in the vacuum of the electromagnetic field without matter the expectation value of the average photon number $$$$ is non zero in the ground state of the polaritonic Hamiltonian $$C_{k \pm}|\mathbf 0>=0$$ similarly to the Hawking radiation in the neighbourhood of the black hole because of the Unruh-Davies effect. One may readily notice that the lower eigenfrequency $$\Omega_{-}$$ becomes imaginary when the coupling constant becomes critical at $$g>1$$ which suggests that Hopfield dielectric will undergo the superradiant phase transition.