View source for Physics:Superradiant phase transition

{{Short description|Process in quantum optics}}
[[File:Superradiant transition.gif|thumb|Schematic plot of the order parameter of the Dicke transition, which is zero in the normal phase and finite in the superradiant phase. The inset shows the free energy in the normal and superradiant phases]]

In [[Physics:Quantum optics|quantum optics]], a '''superradiant phase transition''' is a [[Physics:Phase transition|phase transition]] that occurs in a collection of fluorescent emitters (such as atoms), between a state containing few electromagnetic excitations (as in the [[Physics:Vacuum state|electromagnetic vacuum]]) and a [[Physics:Superradiance|superradiant state]] with many electromagnetic excitations trapped inside the emitters. The superradiant state is made thermodynamically favorable by having strong, coherent interactions between the emitters.

The superradiant phase transition was originally predicted by the [[Physics:Dicke model|Dicke model]] of [[Physics:Superradiance|superradiance]], which assumes that atoms have only two energetic levels and that these interact with only one mode of the electromagnetic field.<ref>{{cite journal|doi= 10.1016/0003-4916(73)90039-0|title=On the superradiant phase transition for Molecules in a Quantized Radiation Field: Dicke Maser Model|author=Hepp, Klaus|author2=Lieb, Elliott H.|journal=Annals of Physics|volume=76|pages=360–404|bibcode =  1973AnPhy..76..360H|year = 1973|issue=2 }}</ref><ref>{{cite journal|doi= 10.1103/PhysRevA.7.831|title=Phase Transition in the Dicke Model of Superradiance|author=Wang, Y. K.|author2=Hioe, F. T|journal=Physical Review A|volume=7|pages=831–836|bibcode =  1973PhRvA...7..831W|year = 1973|issue=3 }}</ref>
The phase transition occurs when the strength of the interaction between the atoms and the field 
is greater than the energy of the non-interacting part of the system.  (This is similar to the case of [[Physics:Superconductivity|superconductivity]] in [[Physics:Ferromagnetism|ferromagnetism]], which leads
to the dynamic interaction between ferromagnetic atoms and the spontaneous ordering of excitations below the critical temperature.)
The collective [[Physics:Lamb shift|Lamb shift]], relating to the system of atoms interacting with the vacuum fluctuations, becomes comparable to the energies of atoms alone, and the vacuum
fluctuations cause the spontaneous self-excitation of matter.

The transition can be readily understood by the use of the Holstein-Primakoff transformation<ref>{{cite journal|doi= 10.1103/PhysRevA.87.023813|title=Superradiant phase transitions with three-level systems|author=Baksic, Alexandre|author2=Nataf, Pierre|author3=Ciuti, Cristiano|journal=Physical Review A|volume=87|pages=023813–023813–5|year = 2013|issue=2 |arxiv = 1206.3213 |bibcode=2013PhRvA..87b3813B|s2cid=7999910 }}</ref> applied to a [[Physics:Two-state quantum system|two-level atom]]. 
As a result of this transformation, the atoms become Lorentz [[Physics:Quantum harmonic oscillator|harmonic oscillators]] 
with frequencies equal to the difference between the energy levels. The whole system then simplifies to a system of interacting harmonic oscillators of atoms, and the field known as [[Physics:Hopfield dielectric|Hopfield dielectric]] which further predicts in the normal state 
[[Physics:Polaron|polaron]]s for photons or [[Physics:Polariton|polariton]]s.
If the interaction with the field is so strong that the system collapses in the harmonic approximation 
and [[Physics:Dispersion relation|complex polariton frequencies]] (soft modes) appear, then the physical system with nonlinear terms of the higher order
becomes the system with the [[Physics:Spontaneous symmetry breaking|Mexican hat-like potential]], and will undergo 
[[Physics:Quantum paraelectricity|ferroelectric-like]] phase transition.<ref>{{cite journal|doi= 10.1016/0375-9601(76)90411-4|title=Appearance of Collective Polarisation as a Result of Phase Transition in an Ensemble of Two-level Atoms Interacting Through Electromagnetic Field|author=Emaljanov, V. I.|author2=Klimontovicz, Yu. L.|journal=Physics Letters A|volume=59|pages=366–368|bibcode = 1976PhLA...59..366E|issue=5 |year = 1976}}</ref>
In this model, the system is mathematically equivalent for one [[Physics:Normal mode|mode]] of excitation to the [[Physics:Trojan wave packet|Trojan wave packet]], 
when the circularly polarized field intensity corresponds to the electromagnetic coupling constant. Above the critical value, it changes to the unstable motion of the [[Physics:Ionization|ionization]].

The superradiant phase transition was the subject of a wide discussion as to whether or not it is only a result of the simplified model of the matter-field interaction; and if it can occur for the real physical parameters of physical systems (a [[Physics:No-go theorem|no-go theorem]]).<ref>{{cite journal|doi= 10.1103/PhysRevLett.35.432|title=Phase Transitions, Two Level Atoms, and the <math>A^2</math> Term|author=Rzążewski, K.|author2=Wódkiewicz, K. T|journal=Physical Review Letters|volume=35|pages=432–434|bibcode = 1975PhRvL..35..432R|issue=7 |year = 1975}}</ref><ref>{{cite journal|doi= 10.1103/PhysRevA.19.301|title=No-go theorem concerning the superradiant phase transition in atomic systems|author=Bialynicki-Birula, Iwo|author2=Rzążewski, Kazimierz|journal=Physical Review A|volume=19|pages=301–303|bibcode = 1979PhRvA..19..301B|issue=1 |year = 1979}}</ref> However, both the original derivation and the later corrections leading to nonexistence of the transition – due to [[Physics:Oscillator strength|Thomas–Reiche–Kuhn sum rule]] canceling for the harmonic oscillator the needed inequality to impossible negativity of the interaction – were based on the assumption that the quantum field operators are commuting numbers, and the atoms do not interact with the static Coulomb forces. This is generally not true like in case of [[Physics:Bohr–van Leeuwen theorem|Bohr–van Leeuwen theorem]] and the classical non-existence of Landau diamagnetism. The negating results were also the consequence of using the simple Quantum Optics models of the electromagnetic field-matter interaction but 
not the more realistic [[Physics:Condensed matter physics|Condenced Matter]] models like for example the superconductivity model of the 
[[Physics:BCS theory|BCS]] but with the [[Physics:Phonon|phonon]]s replaced by [[Physics:Photons|photons]] to first obtain the collective [[Physics:Polariton|polariton]]s. The return of the transition basically occurs because the inter-atom dipole-dipole or generally the electron-electron Coulomb interactions are never negligible in the condensed and even more in the superradiant matter density regime and the Power-Zienau unitary transformation eliminating the quantum vector potential in the minimum-coupling Hamiltonian transforms the Hamiltonian exactly to the form used when it was discovered and without the square of the vector potential which was later claimed to prevent it. Alternatively within the full quantum mechanics including the electromagnetic field  the generalized [[Physics:Bohr–van Leeuwen theorem|Bohr–van Leeuwen theorem]] does not work and
the electromagnetic interactions cannot be eliminated while they only change the <math>\mathbf{p} \cdot \mathbf{A}</math> vector potential coupling to the electric field <math>\mathbf{x} \cdot \mathbf{E}</math> coupling and alter the effective electrostatic interactions. It can be observed in model systems like  [[Physics:Bose–Einstein condensate|Bose–Einstein condensate]]s<ref>{{cite journal|title=Dicke quantum phase transition with a superfluid gas in an optical cavity|author=Baumann, Kristian|author2=Guerlin, Christine |author3= Brennecke, Ferdinand|author4= Esslinger,Tilman  |journal=Nature|volume=464|pages=1301–1306|year = 2010|issue=7293 |doi=10.1038/nature09009 |pmid=20428162 |arxiv=0912.3261|bibcode=2010Natur.464.1301B |s2cid=205220396 }}</ref> 
and artificial atoms.<ref>{{cite journal|doi= 10.1103/PhysRevA.87.013616|title=Finite-temperature Dicke phase transition of a Bose-Einstein condensate in an optical cavity|author=Zhang, Yuanwei|author2=Lian, Jinling|author3=Liang, J.-Q.|author4=Chen, Gang|author5=Zhang, Chuanwei|author6=Suotang, Jia|journal=Physical Review A|volume=87|pages=013616–013616–6|bibcode = 2013PhRvA..87a3616Z|year = 2013|issue=1 |arxiv = 1202.4125 |s2cid=38789923 }}</ref><ref>{{cite journal|title=Superradiant Phase Transitions and the Standard Description of Circuit QED|author=Viehmann, Oliver|author2=von Delft, Jan|author3=Marquard, Florian|journal=Physical Review Letters|volume=107|pages=113602–113602–5|bibcode =2011PhRvL.107k3602V |issue=7 |year = 2011|doi=10.1103/physrevlett.107.113602 |pmid=22026666 |arxiv=1103.4639|s2cid=22747713 }}</ref>

==Theory==

===Criticality of  linearized Jaynes-Cummings model===

A superradiant phase transition is formally predicted by the critical behavior of the resonant Jaynes-Cummings model, describing the interaction of only one atom with one mode of the electromagnetic field.
Starting from the exact Hamiltonian of the Jaynes-Cummings model at resonance

:<math>\hat{H}_{\text{JC}} = \hbar \omega \hat{a}^{\dagger}\hat{a}
+\hbar \omega \frac{\hat{\sigma}_z}{2}
+\frac{\hbar \Omega}{2} \left(\hat{a}\hat{\sigma}_+
+\hat{a}^{\dagger}\hat{\sigma}_-
+\hat{a}\hat{\sigma}_-
+\hat{a}^{\dagger}\hat{\sigma}_+\right),</math>

Applying the Holstein-Primakoff transformation for two spin levels, 
replacing the spin raising and lowering operators by those for the harmonic oscillators
:<math>\hat{\sigma}_- \approx \hat{b}</math>
:<math>\hat{\sigma}_+\approx \hat{b}^{\dagger}</math>
:<math>\hat{\sigma}_z\approx 2 \hat{b}^{\dagger}\hat{b}</math>
one gets the Hamiltonian of two coupled harmonic-oscillators:
:<math>\hat{H}_{\text{JC}} = \hbar \omega \hat{a}^{\dagger}\hat{a}
+\hbar \omega \hat{b}^{\dagger}\hat{b}
+\frac{\hbar \Omega}{2} \left(\hat{a}\hat{b}^{\dagger}
+\hat{a}^{\dagger}\hat{b}
+\hat{a}\hat{b}
+\hat{a}^{\dagger}\hat{b}^{\dagger} \right),</math>
which readily can be diagonalized.
Postulating its normal form 
:<math>\hat{H}_{\text{JC}}=\Omega_+ \hat{A_+}^{\dagger}\hat{A_+}+\Omega_- \hat{A_-}^{\dagger}\hat{A_-}+C</math>
where
:<math>\hat{A_{\pm}}=c_{\pm1} \hat{a} + c_{\pm2} \hat{a}^{\dagger} + c_{\pm3} \hat{b} + c_{\pm4} \hat{b}^{\dagger}</math>
one gets the eigenvalue equation
:<math>[\hat{A_{\pm}},\hat{H}_{\text{JC}}]=\Omega_{\pm}A</math>
with the solutions
:<math>\Omega_{\pm}=\omega \sqrt{1 \pm \frac{\Omega}{\omega}}</math>
The system collapses when one of the frequencies becomes imaginary, i.e. when
:<math>\Omega>\omega</math>
or when the atom-field coupling is stronger than the frequency of the mode and atom oscillators.
While there are physically higher terms in the true system, the system in this regime will therefore undergo the 
phase transition.

===Criticality of Jaynes-Cummings model===

The simplified Hamiltonian of the Jaynes-Cummings model, neglecting the counter-rotating terms, is

:<math>\hat{H}_{\text{JC}} = \hbar \omega \hat{a}^{\dagger}\hat{a}
+\hbar \omega \frac{\hat{\sigma}_z}{2}
+\frac{\hbar \Omega}{2} \left(\hat{a}\hat{\sigma}_+
+\hat{a}^{\dagger}\hat{\sigma}_-\right),</math>

and the energies for the case of zero detuning are

:<math>E_{\pm}(n) = \hbar\omega \left(n+\frac{1}{2}\right) \pm \frac{1}{2} \hbar\Omega(n),</math>

:<math>\Omega(n) = \Omega \sqrt{n+1}</math>

where <math>\Omega</math> is the [[Physics:Rabi frequency|Rabi frequency]].
One can approximately calculate the canonical partition function

: <math> Z = \sum_{ \pm , n} \mathrm{e}^{- \beta E_{\pm}(n)} \approx \sum_{ \pm} \int \mathrm{e}^{- \beta E_{\pm}(n)} dn=\int \mathrm{e}^{\Phi (n)} dn</math>,
where the discrete sum was replaced by the integral.

The normal approach is that the latter integral is calculated by the Gaussian approximation around
the maximum of the exponent:

: <math>\frac{\partial \Phi(n)} {\partial n}=0</math>

: <math>\Phi(n)=-\beta \hbar\omega \left(n+\frac{1}{2}\right)+\log 2 \cosh \frac{\hbar \Omega(n) \beta} {2} </math>

This leads to the critical equation

: <math>\tanh \frac{\hbar \Omega(n) \beta} {2}= 4 \frac{\omega}{\Omega}\sqrt{n+1}</math>

This has the solution only if

: <math>\Omega>4 \omega</math>

which means that the normal, and the superradiant phase, exist only if the field-atom coupling is significantly stronger 
than the energy difference between the atom levels.
When the condition is fulfilled, the equation gives the solution for the order parameter <math>n</math> depending on the inverse of the 
temperature <math>1/\beta</math>, which means non-vanishing ordered field mode.
Similar considerations can be done in true thermodynamic limit of the infinite number of atoms.

===Instability of the classical electrostatic model===

The better insight on the nature of the superradiant phase transition as well on the physical value of the critical parameter which must be exceeded 
in order for the transition to occur may be obtained by studying the classical stability of the system of the charged classical harmonic oscillators in the 3D space interacting only with the electrostatic repulsive forces for example between electrons in the locally harmonic oscillator potential. Despite of the original model of the superradiance the quantum electromagnetic field is totally neglected here. The oscillators may be assumed to be placed for example on the [[Cubic crystal system|cubic lattice]] with the lattice constant <math>a</math> in the analogy to the crystal system of the condensed matter.
The worse scenario of the defect of the absence of the two out-of-the-plane motion-stabilizing electrons from the 6-th nearest neighbors of a chosen electron is assumed while the four nearest electrons are first assumed to be rigid in space and producing the anti-harmonic potential in the direction perpendicular to the plane of the all five electrons. The condition of the [[Lyapunov stability|instability of motion]] of the chosen electron is that the net potential being the superposition of the
harmonic oscillator potential and the quadratically expanded Coulomb potential from the four electrons is negative i.e.
: <math>\frac{m \omega^2}{2}- \frac{1}{2} \times 4 \times \frac {e^2}{4 \pi \epsilon_0} \frac {1}{a^3}<0</math>
or
: <math>\frac{e^2}{\pi \epsilon_0 m \omega^2}\frac {1}{a^3}>1</math>
Making it artificially quantum by multiplying the numerator and the denominator of the fraction by the <math>\hbar</math> one obtains 
the condition 
: <math>\frac{2}{\pi} \frac{|D_{12}|^2}{E_{12} \epsilon_0} \left( \frac{N}{V} \right)>1</math>
where 
: <math>|D_{12}|^2 = \frac{e^2 \hbar} {2 m \omega}</math>
is the square of the [[Physics:Sum rule in quantum mechanics|dipole transition strength]] between the ground state and the first excited state of the [[Physics:Quantum harmonic oscillator|quantum harmonic oscillator]], 
: <math>E_{12}=\hbar \omega</math>
is the energy gap between consecutive levels and it is also noticed that 
: <math>\frac {1}{a^3}=\frac{N}{V}</math>
is the spatial density of the oscillators.
The condition is almost identical to this obtained in the original discovery of the superradiant phase transition when replacing the harmonic oscillators with two level atoms with the same distance between the energy levels, dipole transition strength, and the density which means that it occurs in the regime 
when the Coulomb interactions between electrons dominate over locally harmonic oscillatory influence of the atoms. It that sense the free 
[[Physics:Fermi gas|electron gas]] with <math>\omega=0</math>
is also purely superradiant.

The critical inequality rewritten yet differently
: <math>\omega<\sqrt{\frac {e^2}{m \pi \epsilon_0} \frac {N}{V}}\approx\sqrt{\frac {e^2}{m \epsilon_0} \frac {N}{V}}</math>
expresses the fact that superradiant phase transition occurs when the frequency of the binding atomic oscillators is lower than 
so called electron gas [[Physics:Plasma oscillation|plasma frequency]].

==References==
{{reflist}}

==External links==
* [https://www.youtube.com/watch?v=1DxI-0GJ6fw&list=PL9KfozrSC_kEhUoOSy378v0S_S_uqxP6a&index=8 Old Warsaw School of "No-Go" of the Superradiant Phase Transition], the talk by the 2022 [[Physics:Wigner Medal|Wigner Medal]] recipient Iwo Bialynicki-Birula's former PhD student K. Rzążewski

[[Category:Quantum mechanics]]
[[Category:Phase transitions]]

{{Sourceattribution|Superradiant phase transition}}