Arxiv:1809.08749V3 [Quant-Ph] 5 Dec 2018 a General Guiding Principle in Building the Theory of Fun- Damental Interactions (See, E.G., [1,3])
Total Page:16
File Type:pdf, Size:1020Kb
Resolution of Gauge Ambiguities in Ultrastrong-Coupling Cavity QED Omar Di Stefano,1 Alessio Settineri,2 Vincenzo Macr`ı,1 Luigi Garziano,1 Roberto Stassi,1 Salvatore Savasta,1, 2, ∗ and Franco Nori1, 3 1Theoretical Quantum Physics Laboratory, RIKEN Cluster for Pioneering Research, Wako-shi, Saitama 351-0198, Japan 2Dipartimento di Scienze Matematiche e Informatiche, Scienze Fisiche e Scienze della Terra, Universit`adi Messina, I-98166 Messina, Italy 3Physics Department, The University of Michigan, Ann Arbor, Michigan 48109-1040, USA Gauge invariance is the cornerstone of modern quantum field theory [1{4]. Recently, it has been shown that the quantum Rabi model, describing the dipolar coupling between a two-level atom and a quantized electromagnetic field, violates this principle [5{7]. This widely used model describes a plethora of quantum systems and physical processes under different interaction regimes [8,9]. In the ultrastrong coupling regime, it provides predictions which drastically depend on the chosen gauge. This failure is attributed to the finite-level truncation of the matter system. We show that a careful application of the gauge principle is able to restore gauge invariance even for extreme light- matter interaction regimes. The resulting quantum Rabi Hamiltonian in the Coulomb gauge differs significantly from the standard model and provides the same physical results obtained by using the dipole gauge. It contains field operators to all orders that cannot be neglected when the coupling strength is high. These results shed light on subtleties of gauge invariance in the nonperturbative and extreme interaction regimes, which are now experimentally accessible, and solve all the long-lasting controversies arising from gauge ambiguities in the quantum Rabi and Dicke models [5, 10{18]. The ultrastrong-coupling (USC) between an effective The coupling to the field, obtained performing the re- two-level system (TLS) and the electromagnetic field has placement ∂µ → Dµ in the free Lagrangian, is called the been realized in several solid state systems [8,9]. In \minimal coupling". this regime of quantum light-matter interaction, going It has been shown by several authors [46{50] that ap- beyond weak and strong coupling, the coupling strength proximate models for light-matter interactions derived in becomes comparable to the transition frequencies in the different gauges may lead to different predictions. It has system. Recently, circuit QED experiments involving a also been shown that the convergence of the two-photon single LC-oscillator mode coupled to a flux qubit su- transition rate depends significantly on the choice of the perconducting quantum circuit have realized coupling gauge [51]. strengths larger then the system transition frequencies When the light-matter interaction becomes very [19, 20]. This extreme interaction regime has been de- strong, different gauges can lead to drastically different noted deep strong coupling (DSC). In these interaction predictions, leading to controversies [5, 10{18]. For ex- regimes, several properties of coupled light-matter sys- ample, in the case of several TLSs interacting with a tems change drastically, opening the way to a wealth of single mode of an optical resonator [52], different gauges new intriguing physical effects (see, e.g. [21{35]) which may even lead to very different predictions, such as the may offer opportunities for the development of quantum presence or the absence of a quantum phase transition. technologies [36{44]. One important conclusion that can be drawn from The form of the electron-photon interaction is gauge these controversies is that, once the light-matter coupling dependent (see, e.g., [45]). However, all physical results becomes non-perturbative, the validity of the two-level must be independent of this choice. Gauge invariance is approximation for the atomic dipoles depends explicitly on the choice of gauge. In particular, it has been shown arXiv:1809.08749v3 [quant-ph] 5 Dec 2018 a general guiding principle in building the theory of fun- damental interactions (see, e.g., [1,3]). Let us consider, [6] that in the electric dipole gauge (the multipolar gauge for example, a particle field whose action is invariant un- in the dipole approximation) the two-level approxima- der a global phase change [U(1) invariance]. If this phase tion can be performed as long as the Rabi frequency re- is allowed to depend on the space-time coordinate x, its mains much smaller than the energies of the off-resonant action is not invariant because the coordinate-dependent higher-energy levels. However, it can dramatically fail in phase φ(x) does not commute with the derivatives (corre- the Coulomb gauge, even for systems with an extremely sponding to the four-momentum) in the action. The sym- anharmonic spectrum [6]. This unexpected result is par- metry can be restored replacing these derivatives with ticularly unsatisfactory, since the general procedure to covariant derivatives: Dµ = (∂µ + iqAµ), where q is the derive the multipolar gauge consists of using the min- charge parameter and Aµ is the gauge potential. If the imal coupling replacement first, and then applying the transformation of the particle field ψ → eiφ(x) is associ- Coulomb gauge [45]. ated to a gauge transformation of the four-vector poten- In a recent work [7], considering arbitrary gauges for tial Aµ → Aµ −∂µθ, the local U(1) symmetry is restored. the interaction between a TLS and the electromagnetic 2 field, it has been shown that physical predictions in the tracting great interest, but also provides a general insight USC and DSC regimes strongly depend on the gauge on gauge invariance in extreme interaction regimes. We choice. The gauge which provides the more accurate re- find that the usual strategy, which consists of taking into sults seems to depend on the detuning between the cavity account the nonlocality of the atomic potential, perform- mode resonance and the atomic transition frequency. ing the minimal coupling replacement and developing From all of these previous studies, it is clear that ap- the resulting interaction Hamiltonian up to second order proximations in the description of the matter system, in the vector potential, fails when the coupling strength as e.g., a finite-level truncation, ruin the gauge invari- reaches a significant fraction of the resonance frequencies ance of the theory. In 1970, Ref. [48] pointed out that of the system. We demonstrate that these gauge ambi- gauge ambiguities in the calculation of atomic oscilla- guities can be eliminated for arbitrary coupling strengths tor strengths, originate from the occurrence of nonlocal only by taking into account the approximation-induced potentials determined by the approximation procedures. nonlocality and keeping the resulting interaction Hamil- Since a nonlocal potential in the coordinate representa- tonian to all orders in the vector potential Aˆ . The re- tion is an integral operator, it does not commute with sults presented here solve all the long-lasting controver- the coordinate operator. Indeed, it is easy to show that sies arising from gauge ambiguities in the quantum Rabi it can be expressed as a local momentum-dependent op- and Dicke models. erator V (ˆr, pˆ). This affects the interaction of light with quantum systems described by approximate Hamiltoni- ans. Specifically, in order to introduce the coupling of the RESULTS matter system with the electromagnetic field, the mini- mal replacement rule pˆ → Πˆ = pˆ − Aˆ (r, t) (Aˆ is the Nonlocal potentials vector potential) has to be applied not only to the ki- netic energy terms, but also to the nonlocal potentials in In order to understand why local potentials become the effective Hamiltonian of the particles in the system. nonlocal when the Hilbert space is truncated, let us con- By applying such a procedure, the ambiguities in the sider a one-dimensional potential Vˆ . In the position ba- calculation of approximate matrix elements for electric sis, it can written as dipole transitions can be removed [48]. Moreover, taking Z 0 0 0 into account the nonlocality of the approximate poten- Vˆ = dxdx hx|Vˆ |x i |xihx | (1) tial, two-photon transition rates involving Wannier exci- tons in semiconductors become gauge invariant [49]. Also If the potential is local, its matrix elements can be ˆ 0 0 0 the microscopic quantum theory of excitonic polaritons is written as hx|V |x i ≡ V (x, x ) = W (x)δ(x − x ). If affected by the presence of nonlocal potentials. The ap- |ni is a complete orthonormal basis, the matrix ele- 0 0 plication of the standard minimal coupling replacement ments can be expressed as V (x, x ) = W (x)δ(x − x ) = P ∗ 0 leads to a total Hamiltonian which, besides the usual Aˆ ·pˆ n,n0 Wn,n0 ψn0 (x )ψn(x), where we defined ψn(x) ≡ term, displays an additional diamagnetic term e2Aˆ 2/2m, hx|ni. Notice that the Dirac delta function can be re- where e and m are the electron charge and mass respec- constructed only by keeping all the infinite vectors of the tively. If a limited number of exciton levels are explic- basis. Hence, any truncation of the complete basis can itly included in the model, the long wavelength solution transform a local potential into a nonlocal one. If only of the polariton dispersion cannot be recovered without two states are included, e.g., the two lowest energy levels, introducing ad hoc the Thomas-Reiche-Kuhn sum rule we obtain and truncating the summation consistently [53]. This 0 ∗ 0 ∗ 0 V (x, x ) = W1,0 [ψ (x )ψ1(x) + ψ (x )ψ0(x)] (2) procedure presents some ambiguity and appears to be 0 1 rather artificial. However, if the nonlocality of the ap- where, for simplicity, we assumed parity symmetry and proximate potential is taken into account and the result- real matrix elements. It is evident that adding the ing additional terms are included, up to second order in two terms in Eq.(2), which are products of two smooth the vector potential, the correct dispersion relation is au- wavefunctions it is not possible to reproduce the Dirac- tomatically recovered [54, 55].