Reduced Density-Matrix Approach to Strong Matter-Photon Interaction

Numerically Exact Solution for a Real Polaritonic System under Vibrational Strong Coupling in Thermodynamic Equilibrium: Loss of Light-Matter Entanglement and Enhanced Fluctuations

The first numerically exact simulation of a full ab initio molecular quantum system (HD+) under strong ro-vibrational coupling to a quantized optical cavity mode in thermal equilibrium is presented. Theoretical challenges in describing strongly coupled systems of mixed quantum statistics (bosons and Fermions) are discussed and circumvented by the specific choice of our molecular system. Our numerically exact simulations highlight the absence of zero temperature for the strongly coupled matter and light subsystems, due to cavity-induced noncanonical conditions. Furthermore, we explore the temperature dependency of light-matter quantum entanglement, which emerges for the ground state but is quickly lost already in the deep cryogenic regime. This is in contrast to predictions from the Jaynes-Cummings model, which is the standard starting point to model collective strong-coupling chemistry phenomenologically. Moreover, we find that the fluctuations of matter remain modified by the quantum nature of the thermal and vacuum-field fluctuations for significant temperatures, e.g., at ambient conditions. These observations (loss of entanglement and coupling to quantum fluctuations) have implications for the understanding and control of polaritonic chemistry and materials science, since a semiclassical theoretical description of light-matter interaction becomes reasonable, but the typical (classical) canonical equilibrium assumption for the nuclear subsystem remains violated. This opens the door for quantum fluctuation-induced stochastic resonance phenomena under vibrational strong coupling, which have been suggested as a plausible theoretical mechanism to explain the experimentally observed resonance phenomena in the absence of periodic driving that has not yet been fully understood.

Cavity-Modified Chemiluminescent Reaction of Dioxetane

Chemiluminescence is a thermally activated chemical process that emits a photon of light by forming a fraction of products in the electronic excited state. A well-known example of this spectacular phenomenon is the emission of light in the firefly beetle, where the formation of a four-membered cyclic peroxide compound and subsequent dissociation produce a light-emitting product. The smallest cyclic peroxide, dioxetane, also exhibits chemiluminescence but with a low quantum yield as compared to that of firefly dioxetane. Employing the strong light-matter coupling has recently been found to be an alternative strategy to modify the chemical reactivity. In the presence of an optical cavity, the molecular degrees of freedom greatly mix with the cavity mode to form hybrid cavity-matter states called polaritons. These newly generated hybrid light-matter states manipulate the potential energy surfaces and significantly change the reaction dynamics. Here, we theoretically investigate the effects of a strong light-matter interaction on the chemiluminescent reaction of dioxetane using the extended Jaynes-Cummings model. The cavity couplings corresponding to the electronic and vibrational degrees of freedom have been included in the interaction Hamiltonian. We explore how the cavity alters the ground- and excited-state path energy barriers and reaction rates. Our results demonstrate that the formation of excited-state products in the dioxetane decomposition process can be either accelerated or suppressed, depending on the molecular orientation with respect to the cavity polarization.

Understanding Polaritonic Chemistry from Ab Initio Quantum Electrodynamics

In this review, we present the theoretical foundations and first-principles frameworks to describe quantum matter within quantum electrodynamics (QED) in the low-energy regime, with a focus on polaritonic chemistry. By starting from fundamental physical and mathematical principles, we first review in great detail ab initio nonrelativistic QED. The resulting Pauli-Fierz quantum field theory serves as a cornerstone for the development of (in principle exact but in practice) approximate computational methods such as quantum-electrodynamical density functional theory, QED coupled cluster, or cavity Born-Oppenheimer molecular dynamics. These methods treat light and matter on equal footing and, at the same time, have the same level of accuracy and reliability as established methods of computational chemistry and electronic structure theory. After an overview of the key ideas behind those ab initio QED methods, we highlight their benefits for understanding photon-induced changes of chemical properties and reactions. Based on results obtained by ab initio QED methods, we identify open theoretical questions and how a so far missing detailed understanding of polaritonic chemistry can be established. We finally give an outlook on future directions within polaritonic chemistry and first-principles QED.

Light-responsive organic polaritons from first principles

Controlling the optical properties of light-responsive organic molecules is essential for their application in photonics. We demonstrate how light-responsive organic polaritons formed inside an optical cavity can be used to modify these properties based on first principles. Specifically, we study the excited state properties of the trans-azobenzene molecule and the free base tetraphenyl porphyrin (H2TPP) molecule under weak to strong light-matter coupling. Our results show that the cavity can modulate the dispersion and absorption properties of organic molecules. Compared to the case outside the cavity, the anomalous dispersion of the trans-azobenzene molecule inside the cavity is suppressed and this suppression decreases with increasing coupling strength, showing the potential of strong light-matter coupling in manipulating the optical dipole trap of organic molecules. Moreover, by adjusting the cavity parameters to tune the strength of the light-matter coupling, we achieve free switching between symmetric Lorentz and asymmetric Fano line shapes for H2TPP polaritonic excitations. During the switching between these spectral features, we also find that the cavity can be used to control the spontaneous radiation of organic molecules via the Purcell effect. These findings provide a new pathway to manipulate the optical properties of light-responsive organic molecules.

Electron Correlation in the Iron(II) Porphyrin by Natural Orbital Functional Approximations

The relative stability of the singlet, triplet, and quintet spin states of iron(II) porphyrin (FeP) represents a challenging problem for electronic structure methods. While it is currently accepted that the ground state is a triplet, multiconfigurational wave function-based methods predict a quintet, and density functional approximations vary between triplet and quintet states, leading to a prediction that highly depends on the features of the method employed. The recently proposed Global Natural Orbital Functional (GNOF) aims to provide a balanced treatment between static and dynamic correlation, and together with the previous Piris Natural Orbital Functionals (PNOFs), allowed us to explore the importance of each type of correlation in the stability order of the states of FeP with a method that conserves the spin of the system. It is noteworthy that GNOF correlates all electrons in all available orbitals for a given basis set; in the case of the FeP with a double-ζ basis set as used in this work, this means that GNOF can properly correlate 186 electrons in 465 orbitals, significantly increasing the sizes of systems amenable to multiconfigurational treatment. Results show that PNOF5, PNOF7s, and PNOF7 predict the quintet to have a lower energy than the triplet state; however, the addition of dynamic correlation via second-order Møller-Plesset corrections (NOF-MP2) turns the triplet state to be lower than the quintet state, a prediction also reproduced by GNOF that incorporates much more dynamic correlation than its predecessors.

A perspective on ab initio modeling of polaritonic chemistry: The role of non-equilibrium effects and quantum collectivity

This Perspective provides a brief introduction into the theoretical complexity of polaritonic chemistry, which emerges from the hybrid nature of strongly coupled light-matter states. To tackle this complexity, the importance of ab initio methods is highlighted. Based on those, novel ideas and research avenues are developed with respect to quantum collectivity, as well as for resonance phenomena immanent in reaction rates under vibrational strong coupling. Indeed, fundamental theoretical questions arise about the mesoscopic scale of quantum-collectively coupled molecules when considering the depolarization shift in the interpretation of experimental data. Furthermore, to rationalize recent findings based on quantum electrodynamical density-functional theory (QEDFT), a simple, but computationally efficient, Langevin framework is proposed based on well-established methods from molecular dynamics. It suggests the emergence of cavity-induced non-equilibrium nuclear dynamics, where thermal (stochastic) resonance phenomena could emerge in the absence of external periodic driving. Overall, we believe that the latest ab initio results indeed suggest a paradigmatic shift for ground-state chemical reactions under vibrational strong coupling from the collective quantum interpretation toward a more local, (semi)-classically and non-equilibrium dominated perspective. Finally, various extensions toward a refined description of cavity-modified chemistry are introduced in the context of QEDFT, and future directions of the field are sketched.

Non-Hermitian Cavity Quantum Electrodynamics - Configuration Interaction Singles Approach for Polaritonic Structure with ab initio Molecular Hamiltonians

We combine ab initio molecular electronic Hamiltonians with a cavity quantum electrodynamics model for dissipative photonic modes and apply mean-field theories to the ground- and excited-states of resulting polaritonic systems. In particular, we develop a non-Hermitian configuration interaction singles theory for mean-field ground- and excited-states of the molecular system strongly interacting with a photonic mode, and apply these methods to elucidating the phenomenology of paradigmatic polaritonic systems. We leverage the Psi4Numpy framework to yield open-source and accessible reference implementations of these methods.

Density Functional Theory Transformed into a One-Electron Reduced-Density-Matrix Functional Theory for the Capture of Static Correlation

Density Functional Theory (DFT), the most widely adopted method in modern computational chemistry, fails to describe accurately the electronic structure of strongly correlated systems. Here we show that DFT can be formally and practically transformed into a one-electron reduced-density-matrix (1-RDM) functional theory, which can address the limitations of DFT while retaining favorable computational scaling compared to wave function based approaches. In addition to relaxing the idempotency restriction on the 1-RDM in the kinetic energy term, we add a quadratic 1-RDM-based term to DFT's density-based exchange-correlation functional. Our approach, which we implement by quadratic semidefinite programming at DFT's computational scaling of O(r³), yields substantial improvements over traditional DFT in the description of static correlation in chemical structures and processes such as singlet biradicals and bond dissociations.

Stochastic Variational Approach to Small Atoms and Molecules Coupled to Quantum Field Modes in Cavity QED

In this work, we present a stochastic variational calculation (SVM) of energies and wave functions of few particle systems coupled to quantum fields in cavity QED. The spatial wave function and the photon spaces are optimized by a random selection process. Using correlated basis functions, the SVM approach solves the problem accurately and opens the way to the same precision that is reached the nonlight coupled quantum systems. Examples for a two-dimensional trion and confined electrons as well as for the He atom and the H_{2} molecule are presented showing that the light-matter coupling drastically changes the electronic states.

Foundation of One-Particle Reduced Density Matrix Functional Theory for Excited States

In Phys. Rev. Lett. 2021, 127, 023001 a reduced density matrix functional theory (RDMFT) was proposed for calculating energies of selected eigenstates of interacting many-Fermion systems. Here, we develop a solid foundation for this so-called w-RDMFT and present the details of various derivations. First, we explain how a generalization of the Ritz variational principle to ensemble states with fixed weights w in combination with the constrained search would lead to a universal functional of the one-particle reduced density matrix. To turn this into a viable functional theory, however, we also need to implement an exact convex relaxation. This general procedure includes Valone's pioneering work on ground state RDMFT as the special case w = (1,0, ···). Then, we work out in a comprehensive manner a methodology for deriving a compact description of the functional's domain. This leads to a hierarchy of generalized exclusion principle constraints which we illustrate in great detail. By anticipating their future pivotal role in functional theories and to keep our work self-contained, several required concepts from convex analysis are introduced and discussed.

Ensemble Reduced Density Matrix Functional Theory for Excited States and Hierarchical Generalization of Pauli's Exclusion Principle

We propose and work out a reduced density matrix functional theory (RDMFT) for calculating energies of eigenstates of interacting many-electron systems beyond the ground state. Various obstacles which historically have doomed such an approach to be unfeasible are overcome. First, we resort to a generalization of the Ritz variational principle to ensemble states with fixed weights. This in combination with the constrained search formalism allows us to establish a universal functional of the one-particle reduced density matrix. Second, we employ tools from convex analysis to circumvent the too involved N-representability constraints. Remarkably, this identifies Valone's pioneering work on RDMFT as a special case of convex relaxation and reveals that crucial information about the excitation structure is contained in the functional's domain. Third, to determine the crucial latter object, a methodology is developed which eventually leads to a generalized exclusion principle. The corresponding linear constraints are calculated for systems of arbitrary size.

Controlling the Photostability of Pyrrole with Optical Nanocavities

Strong light-matter coupling provides a new strategy to manipulate the non-adiabatic dynamics of molecules by modifying potential energy surfaces. The vacuum field of nanocavities can couple strongly with the molecular degrees of freedom and form hybrid light-matter states, termed as polaritons or dressed states. The photochemistry of molecules possessing intrinsic conical intersections can be significantly altered by introducing cavity couplings to create new conical intersections or avoided crossings. Here, we explore the effects of optical cavities on the photo-induced hydrogen elimination reaction of pyrrole. Wave packet dynamics simulations have been performed on the two-state, two-mode model of pyrrole, combined with the cavity photon mode. Our results show how the optical cavities assist in controlling the photostability of pyrrole and influence the reaction mechanism by providing alternative dissociation pathways. The cavity effects have been found to be intensely dependent on the resonance frequency. We further demonstrate the importance of the vibrational cavity couplings and dipole-self interaction terms in describing the cavity-modified non-adiabatic dynamics.

Polaritonic Chemistry: Collective Strong Coupling Implies Strong Local Modification of Chemical Properties

A fundamental question in the field of polaritonic chemistry is whether collective coupling implies local modifications of chemical properties scaling with the ensemble size. Here we demonstrate from first-principles that an impurity present in a collectively coupled chemical ensemble features such locally scaling modifications. In particular, we find the formation of a novel dark state for a nitrogen dimer chain of variable size, whose local chemical properties are altered considerably at the impurity due to its embedding in the collectively coupled environment. Our simulations unify theoretical predictions from quantum optical models (e.g., collective dark states and bright polaritonic branches) with the single molecule quantum chemical perspective, which relies on the (quantized) redistribution of charges leading to a local hybridization of light and matter. Moreover, our findings suggest that recently developed ab initio methods for strong light-matter coupling are suitable to access these local polaritonic effects and provide a detailed understanding of photon-modified chemistry.

Role of Cavity Losses on Nonadiabatic Couplings and Dynamics in Polaritonic Chemistry

We present a non-Hermitian formulation of the polaritonic structure of azobenzene strongly coupled to a photonic mode that explicitly accounts for the fleeting nature of the photon-molecule interaction. This formalism reveals that the polaritonic nonadiabatic couplings that facilitate cis-trans isomerization can be dramatically modified by photonic dissipation. We perform Fewest-Switches Surface Hopping dynamics on the surfaces that derive from our non-Hermitian formalism and find that the polaritonic isomerization yields are strongly suppressed for moderate dissipation rates and that cavity-free isomerization dynamics are recovered under large dissipation rates. These findings highlight the important role that the finite lifetime of photonic degrees of freedom play in polaritonic chemistry.

Virial Relations for Electrons Coupled to Quantum Field Modes

In this work, we present a set of virial relations for many electron systems coupled to both classical and quantum fields, described by the Pauli-Fierz Hamiltonian in dipole approximation and using length gauge. Currently, there is growing interest in solutions of this Hamiltonian because of its relevance for describing molecular systems strongly coupled to photonic modes in cavities and in the possible modification of chemical properties of such systems compared to the ones in free space. The relevance of such virial relations is demonstrated by showing a connection to mass renormalization and by providing an exact way to obtain total energies from potentials in the framework of quantum electrodynamical density functional theory.

Chemistry in Quantum Cavities: Exact Results, the Impact of Thermal Velocities, and Modified Dissociation

In recent years tremendous progress in the field of light-matter interactions has unveiled that strong coupling to the modes of an optical cavity can alter chemistry even at room temperature. Despite these impressive advances, many fundamental questions of chemistry in cavities remain unanswered. This is also due to a lack of exact results that can be used to validate and benchmark approximate approaches. In this work we provide such reference calculations from exact diagonalization of the Pauli-Fierz Hamiltonian in the long-wavelength limit with an effective cavity mode. This allows us to investigate the reliability of the ubiquitous Jaynes-Cummings model not only for electronic but also for the case of ro-vibrational transitions. We demonstrate how the commonly ignored thermal velocity of charged molecular systems can influence chemical properties while leaving the spectra invariant. Furthermore, we show the emergence of new bound polaritonic states beyond the dissociation energy limit.

Light-Matter Hybrid-Orbital-Based First-Principles Methods: The Influence of Polariton Statistics

A detailed understanding of strong matter–photon interactions requires first-principle methods that can solve the fundamental Pauli–Fierz Hamiltonian of nonrelativistic quantum electrodynamics efficiently. A possible way to extend well-established electronic-structure methods to this situation is to embed the Pauli–Fierz Hamiltonian in a higher-dimensional light–matter hybrid auxiliary configuration space. In this work we show the importance of the resulting hybrid Fermi–Bose statistics of the polaritons, which are the new fundamental particles of the “photon-dressed” Pauli–Fierz Hamiltonian for systems in cavities. We show that violations of these statistics can lead to unphysical results. We present an efficient way to ensure the correct statistics by enforcing representability conditions on the dressed one-body reduced density matrix. We further present a general prescription how to extend a given first-principles approach to polaritons and as an example introduce polaritonic Hartree–Fock theory. While being a single-reference method in polariton space, polaritonic Hartree–Fock is a multireference method in the electronic space, i.e., it describes electronic correlations. We also discuss possible applications to polaritonic QEDFT. We apply this theory to a lattice model and find that, the more delocalized the bound-state wave function of the particles is, the stronger it reacts to photons. The main reason is that within a small energy range, many states with different electronic configurations are available as opposed to a strongly bound (and hence energetically separated) ground-state wave function. This indicates that under certain conditions coupling to the quantum vacuum of a cavity can indeed modify ground state properties.

Relevance of the Quadratic Diamagnetic and Self-Polarization Terms in Cavity Quantum Electrodynamics

Experiments at the interface of quantum optics and chemistry have revealed that strong coupling between light and matter can substantially modify the chemical and physical properties of molecules and solids. While the theoretical description of such situations is usually based on nonrelativistic quantum electrodynamics, which contains quadratic light-matter coupling terms, it is commonplace to disregard these terms and restrict the treatment to purely bilinear couplings. In this work, we clarify the physical origin and the substantial impact of the most common quadratic terms, the diamagnetic and self-polarization terms, and highlight why neglecting them can lead to rather unphysical results. Specifically, we demonstrate their relevance by showing that neglecting these terms leads to the loss of gauge invariance, basis set dependence, disintegration (loss of bound states) of any system in the basis set limit, unphysical radiation of the ground state, and an artificial dependence on the static dipole. Besides providing important guidance for modeling of strongly coupled light-matter systems, the presented results also indicate conditions under which those effects might become accessible.

Octopus, a computational framework for exploring light-driven phenomena and quantum dynamics in extended and finite systems

Over the last few years, extraordinary advances in experimental and theoretical tools have allowed us to monitor and control matter at short time and atomic scales with a high degree of precision. An appealing and challenging route toward engineering materials with tailored properties is to find ways to design or selectively manipulate materials, especially at the quantum level. To this end, having a state-of-the-art ab initio computer simulation tool that enables a reliable and accurate simulation of light-induced changes in the physical and chemical properties of complex systems is of utmost importance. The first principles real-space-based Octopus project was born with that idea in mind, i.e., to provide a unique framework that allows us to describe non-equilibrium phenomena in molecular complexes, low dimensional materials, and extended systems by accounting for electronic, ionic, and photon quantum mechanical effects within a generalized time-dependent density functional theory. This article aims to present the new features that have been implemented over the last few years, including technical developments related to performance and massive parallelism. We also describe the major theoretical developments to address ultrafast light-driven processes, such as the new theoretical framework of quantum electrodynamics density-functional formalism for the description of novel light-matter hybrid states. Those advances, and others being released soon as part of the Octopus package, will allow the scientific community to simulate and characterize spatial and time-resolved spectroscopies, ultrafast phenomena in molecules and materials, and new emergent states of matter (quantum electrodynamical-materials).