We present a general theoretical framework to capture light-matter interactions beyond the electric-dipole approximation (EDA), applicable to extended nano- and microscale materials interacting with spatially structured electric fields without truncation at finite multipolar order. The approach is based on the Power-Zienau-Woolley (PZW) Hamiltonian for light-matter interactions and a representation of the material's Hamiltonian in a basis of maximally localized Wannier functions (MLWFs), obtainable from first-principles calculations. We utilize this approach to clarify the limitations of the ubiquitous dipole approximation. We consider electric fields with both uniform and non-uniform intensities and a range of ratios of system size to the wavelength of light. Through this analysis, we identify the conditions under which the EDA breaks down, leading to significant errors in the light-induced dynamics. Contrary to conventional belief, we find that the EDA is remarkably robust for uniformly illuminated 1-D or 2-D materials when light propagates perpendicular to the material. For 3-D materials or non-perpendicular illumination of lower-dimensional materials, conventional wisdom holds and the EDA begins to break down when the wavelength becomes comparable to the system size. Furthermore, the EDA fails when the material is illuminated partially or non-uniformly. For slowly varying field intensities this failure can be corrected by finite-order multipolar corrections. However, for fields that vary substantially, correcting via multipolar terms becomes computationally impractical. In contrast, our approach captures beyond-dipole light-matter interactions at the computational cost of a standard dipole calculation. This efficiency enables accurate first-principles simulations of spatially structured light-matter dynamics in nanoscale devices, quantum materials, and interfaces.
We present a general theoretical framework to capture light-matter interactions beyond the electric-dipole approximation (EDA), applicable to extended nano- and microscale materials interacting with spatially structured electric fields without truncation at finite multipolar order. The approach is based on the Power-Zienau-Woolley (PZW) Hamiltonian for light-matter interactions and a representation of the material's Hamiltonian in a basis of maximally localized Wannier functions (MLWFs), obtainable from first-principles calculations. We utilize this approach to clarify the limitations of the ubiquitous dipole approximation. We consider electric fields with both uniform and non-uniform intensities and a range of ratios of system size to the wavelength of light. Through this analysis, we identify the conditions under which the EDA breaks down, leading to significant errors in the light-induced dynamics. Contrary to conventional belief, we find that the EDA is remarkably robust for uniformly illuminated 1-D or 2-D materials when light propagates perpendicular to the material. For 3-D materials or non-perpendicular illumination of lower-dimensional materials, conventional wisdom holds and the EDA begins to break down when the wavelength becomes comparable to the system size. Furthermore, the EDA fails when the material is illuminated partially or non-uniformly. For slowly varying field intensities this failure can be corrected by finite-order multipolar corrections. However, for fields that vary substantially, correcting via multipolar terms becomes computationally impractical. In contrast, our approach captures beyond-dipole light-matter interactions at the computational cost of a standard dipole calculation. This efficiency enables accurate first-principles simulations of spatially structured light-matter dynamics in nanoscale devices, quantum materials, and interfaces.
The photonic environment can significantly influence emission properties and interactions among atomic systems. In such scenarios, frequently the electric dipole approximation is assumed that is justified as long as the spatial extent of the atomic system is negligible compared to the spatial variations of the field. While this holds true for many canonical systems, it ceases to be applicable for more contemporary nanophotonic structures. To go beyond the electric dipole approximation, we propose and develop in this article an analytical framework to describe the impact of the photonic environment on emission and interaction properties of atomic systems beyond the electric dipole approximation. Particularly, we retain explicitly magnetic dipolar and electric quadrupolar contributions to the light-matter interactions. We exploit a field quantization scheme based on electromagnetic Green's tensors, suited for dispersive materials. We obtain expressions for spontaneous emission rate, Lamb shift, multipole-multipole shift and superradiance rate, all being modified with dispersive environment. The considered influence could be substantial for suitably tailored nanostructured photonic environments, as demonstrated exemplarily.
The loss of gauge invariance in models of light-matter interaction which arises from material and photonic space truncation can pose significant challenges to conventional quantum optical models when matter and light strongly hybridize. In structured photonic environments, necessary in practice to achieve strong light-matter coupling, a rigorous model of field quantization within the medium is also needed. Here we use the framework of macroscopic QED by quantizing the fields in an arbitrary material system, with a spatially dependent dispersive and absorptive dielectric, starting from a fundamental light-matter action. We truncate the material and mode degrees of freedom while respecting the gauge principle by imposing a partial gauge-fixing constraint during canonical quantization, which admits a large number of gauges including the Coulomb and multipolar gauges commonly used in quantum optics. We also consider gauge conditions with explicit time dependence, enabling us to unambiguously introduce additional phenomenologically time-dependent light-matter interactions in any gauge. Our results allow one to derive rigorous nonrelativistic models of ultrastrong light-matter interactions in structured photonic environments with no gauge ambiguity. Results for two-level systems and the dipole approximation are discussed, as well as how to go beyond the dipole approximation for effective single-particle models. By comparing with the limiting case of an inhomogeneous dielectric, where dispersion and absorption can be neglected and the fields can be expanded in terms of the generalized transverse eigenfunctions of the dielectric, we show how lossy systems can introduce an additional gauge ambiguity, which we resolve and predict to have fundamental implications for open quantum system models. Finally, we show how observables in mode-truncated systems can be calculated without ambiguity by using a simple gauge-invariant model of photodetection.
The photonic environment can significantly influence emission properties and interactions among atomic systems. In such scenarios, frequently the electric dipole approximation is assumed that is justified as long as the spatial extent of the atomic system is negligible compared to the spatial variations of the field. While this holds true for many canonical systems, it ceases to be applicable for more contemporary nanophotonic structures. To go beyond the electric dipole approximation, we propose and develop in this article an analytical framework to describe the impact of the photonic environment on emission and interaction properties of atomic systems beyond the electric dipole approximation. Particularly, we retain explicitly magnetic dipolar and electric quadrupolar contributions to the light-matter interactions. We exploit a field quantization scheme based on electromagnetic Green’s tensors, suited for dispersive materials. We obtain expressions for spontaneous emission rate, Lamb shift, multipole-multipole shift and superradiance rate, all being modified with dispersive environment. The considered influence could be substantial for suitably tailored nanostructured photonic environments, as demonstrated exemplarily.
In this work, we detail different approaches to treat multi-mode photonic environments within non-relativistic quantum electrodynamics in the long-wavelength approximation efficiently. Specifically we show that for equilibrium properties of coupled light-matter systems, we can approximately capture the effects of multi-mode photonic environments on matter systems by either only keeping the polarization part of the electric field in the length-gauge formulation or by a few effective modes. We present a comprehensive set of approximation methods designed to accurately capture equilibrium phenomena in quantum light-matter systems across a range of complex photonic environments, from weak to strong coupling. These methods are applied to atomic and molecular models as well as to a two-dimensional quantum ring, demonstrating the versatility of our approach and laying the groundwork for first-principles simulations of real materials in cavity quantum electrodynamics.
Abstract Most theoretical studies for correlated light–matter systems are performed within the long-wavelength limit, i.e., the electromagnetic field is assumed to be spatially uniform. In this limit the so-called length-gauge transformation for a fully quantized light–matter system gives rise to a dipole self-energy term in the Hamiltonian, i.e., a harmonic potential of the total dipole matter moment. In practice this term is often discarded as it is assumed to be subsumed in the kinetic energy term. In this work we show the necessity of the dipole self-energy term. First and foremost, without it the light–matter system in the long-wavelength limit does not have a ground-state, i.e., the combined light–matter system is unstable. Further, the mixing of matter and photon degrees of freedom due to the length-gauge transformation, which also changes the representation of the translation operator for matter, gives rise to the Maxwell equations in matter and the omittance of the dipole self-energy leads to a violation of these equations. Specifically we show that without the dipole self-energy the so-called ‘depolarization shift’ is not properly described. Finally we show that this term also arises if we perform the semi-classical limit after the length-gauge transformation. In contrast to the standard approach where the semi-classical limit is performed before the length-gauge transformation, the resulting Hamiltonian is bounded from below and thus supports ground-states. This is very important for practical calculations and for density-functional variational implementations of the non-relativistic QED formalism. For example, the existence of a combined light–matter ground-state allows one to calculate the Stark shift non-perturbatively.
In view of the possibilities of laser spectroscopy with very short pulses and relatively high frequencies, we examine aspects of the theory of atom-field interactions that are related to violation of the condition of the long-wavelength approximation (LWA) according to which kr⪡1, where r is ``of the order of atomic dimensions'' and k is the magnitude of the wave vector. On- and off-resonance transitions are considered, with kr being larger than unity due to the large extent of the two wave functions involved in the coupling matrix element. The implementation of the analysis uses bound-free transition matrix elements with the n=50 and 80 hydrogenic functions as initial states and values of k up to 27eV, which is sufficient to produce kr>1, thereby rendering the LWA inoperative. In spite of this, it is shown that, for on-resonance transitions, the results from the use of the well-known multipole operators resulting from the LWA [e.g., the electric dipole approximation (EDA)] agree with those from the exact expressions derived here from the application of the multipolar Hamiltonian. This numerical agreement is proven analytically. As a test of the kind of convergence of the multipole series expansion for small values of k and large r, it is shown that the lowest-order ratio of electric dipole to quadrupole matrix elements decreases rapidly within 1.0 atomic unit above threshold. Finally, it is shown that off-resonance couplings lead to differences between the full-interaction operator and the EDA which cannot be neglected. In the extreme case of intrashell couplings for the n=50 shell, calculation shows that the 50p state is coupled to angular momentum states up to $\ensuremath{\ell}=21$, compared to the electric dipole coupling of $\ensuremath{\Delta}\ensuremath{\ell}=\ifmmode\pm\else\textpm\fi{}1$.