Psi4 developers create psi4-cbo within Psi4 framework; joint nuclear-photonic Hessian analysis
Enabling Polaritonic Electronic Structure Theory for Psi4:The psi4-cbo package
Enabling Polaritonic Electronic Structure Theory for Psi4: The psi4-cbo package
We present psi4-cbo , an open-source extension of the Psi4 program for simulations within the cavity Born-Oppenheimer approximation (CBOA) . By treating the cavity photon mode on the same footing as the nuclear coordinates, psi4-cbo incorporates light-matter coupling into standard electronic structure workflows. These workflows include self-consistent field calculations, response theory, analytic gradients, joint Hessians, and geometry optimization. The package is designed for reproducibility and is documented with working examples for its main workflows. To evaluate the approximations inherent to the CBOA framework, we use an exact quantum-electrodynamical full configuration interaction (QED-FCI) reference. Resolving the correlation energy into electron-electron and electron-photon channels enables us to precisely determine the point at which the CBOA becomes quantitatively sufficient. Additionally, we highlight that, since self-consistent treatments dress the cavity modes, satisfying experimental resonance conditions is method-dependent. As a practical application, we investigate the chromophore acrolein under vibrational strong coupling and detail how to locate resonances on the dressed photon modes rather than the bare molecular bands.
Tightly confined photon modes in optical or plasmonic cavities enable strong light-matter coupling. When applied to molecules, electronic states or molecular vibrations hybridize with the photon modes and form polariton states [ 1 , 2 ] . This phenomenon has been reported to alter reaction landscapes [ 3 ] , modify enzyme activity [ 4 ] , redirect stereoselectivity [ 5 ] , and influence charge and ionic transport [ 6 , 7 ] . There is increasing evidence that these observed effects are collective in nature, meaning that, rather than acting on isolated molecules, the cavity influences how molecules organize. For example, vibrational strong coupling (VSC) can alter the supramolecular assembly of conjugated polymers [ 8 ] , poly(phenylene ethynylene)s [ 9 ] , and porphyrins [ 10 ] , as well as drive the coassembly of DNA origami [ 11 ] . Similarly, electronic strong coupling (ESC) has been shown to modify ground-state intermolecular interactions in self-assembled chlorins [ 12 ] . Because these phenomena involve ensembles that share a single cavity mode, they place strict demands on theoretical descriptions.
Although VSC and ESC offer promising new tools for chemistry, the underlying mechanisms are not yet fully understood, and experimental clarification of the underlying mechanisms remains challenging [ 13 , 14 ] . To explain these experimental observations, ab initio theory may be required that goes beyond minimal or phenomenological models and includes the light-matter coupling at the level of electronic structures. The Pauli-Fierz Hamiltonian and its cavity Born-Oppenheimer approximation (CBOA) or quantum electrodynamics (QED) reductions provide this foundation [ 15 , 16 , 17 , 18 ] . However, applying these methods to correlated wave functions, analytic derivatives, and large molecular ensembles remains challenging [ 19 , 20 , 21 ] . Recent progress includes detailed analyses of the CBOA [ 22 , 23 ] , its extension to perturbative electron-photon correlation [ 24 ] , and studies on how cavity-modified reactivity holds up under fully quantum treatments [ 25 , 26 , 27 ] .
Efficient software packages have established a strong foundation for polaritonic electronic structure theory. In particular, the Octopus code [ 28 , 29 , 30 , 31 ] provides extensive capabilities with a focus on quantum electrodynamics and solid-state density functional theory, while the e T e^{T} program package [ 32 , 33 ] offers comprehensive implementations for fully quantized photon methods [ 34 , 1 ] and CBOA [ 16 , 2 ] for molecular systems. In this paper, we present a complementary tool, psi4-cbo , an open-source package built on Psi4 [ 8 ] and Psi4NumPy [ 38 ] . It is built to be easily extensible and allows for experimentation and analysis of quantities, such as electronic densities and molecular properties that are typically not exposed in other program packages. To achieve rapid prototyping of new theoretical features, psi4-cbo is written in Python. Additionally, the psi4-cbo package specifically features the joint nuclear-photonic Hessian along with its associated analytical framework. After developing and validating the mean-field ansatz, analytic gradients, and vibro-polaritonic spectra in previous studies [ 2 , 39 , 40 , 10 , 42 ] , the methods presented in these studies are available in psi4-cbo . In addition, the program package also introduces post-mean field methods, diagnostic tools, and a comprehensive documentation.
The paper is organized as follows: Section II defines the underlying Pauli-Fierz Hamiltonian and its mean-field approaches cavity Born-Oppenheimer Hartree-Fock (CBO-HF) and quantum-electrodynamical Hartree-Fock (QED-HF) . Section III builds on this foundation by introducing both cavity Born-Oppenheimer configuration interaction (CBO-CI) and quantum-electrodynamical configuration interaction (QED-CI) , discussing their main differences, and building a suite of quantitative correlation diagnostics for mean-field solutions. Section IV outlines the complete technical capabilities of psi4-cbo , and Section IV.1 shows worked examples of the main workflows with the input needed to reproduce the results. The two applications are as follows. Section V discusses the different sources of correlation, systematic diagnostics across different systems, and coupling strengths. Moreover, the parameter regime where the CBOA remains reliable is discussed and where the quantized photon is required. Section VI applies psi4-cbo to acrolein, a prototypical chromophore used to quantitatively validate photon-mode dressing, and compares CBO-HF against cavity Born-Oppenheimer density functional theory (CBO-DFT) for vibro-polaritonic infrared (IR) spectra. Section VII then discusses its limitations, and Section VIII states how the package is obtained, installed, and tested. The paper concludes with a summary of future developments.
II.1 The Pauli-Fierz Hamiltonian and the cavity Born-Oppenheimer approximation
We consider N mod N_{\rm mod} effective cavity modes, each of which is characterized by a frequency, Ο m \omega_{m} , and a coupling vector, π m \bm{\lambda}{m} along the polarization axis, π ^ m \hat{\mathbf{e}}{m} . The coupling vector π m \bm{\lambda}_{m} is defined as
with its magnitude, the coupling strength Ξ» m \lambda_{m} , is set by the effective mode volume V m V_{m} . Atomic units are used throughout, unless otherwise stated. Within the standard Born-Oppenheimer and long-wavelength approximation in the length gauge, the Pauli-Fierz Hamiltonian [ 15 , 16 ] can be written either in terms of Fock states or photon displacement coordinates. In both representations, the nuclei are treated as classical point charges fixed at a fixed geometry π \mathbf{R} . The final Hamiltonian then consists of the electronic Hamiltonian H ^ el β ( π ) \hat{H}_{\rm el}(\mathbf{R}) , the Hamiltonian of the cavity mode, and a light-matter coupling term. The Pauli-Fierz Hamiltonian in the Fock state representation then reads:
where b ^ m β \hat{b}{m}^{\dagger} and b ^ m \hat{b}{m} are the photon creation and annihilation operators of mode m m , and Ο m \omega_{m} is the frequency of the cavity mode. d ^ m \hat{d}_{m} is the dressed dipole operator of mode m m , which is given by the scalar product of the total (electronic and nuclear) molecular dipole operator π ^ \hat{\bm{\mu}} and the coupling vector of Eq. ( 1 ):
The dressed dipole operator is used both in the linear light-matter coupling term and in the dipole self-energy (DSE) term 1 2 β d ^ m 2 \tfrac{1}{2}\hat{d}{m}^{2} . Alternatively, the Pauli-Fierz Hamiltonian can be expressed in terms of photon displacement coordinates q m q{m} ,
where we have used q ^ m = ( b ^ m + b ^ m β )
2 β Ο m \hat{q}{m}=(\hat{b}{m}+\hat{b}_{m}^{\dagger})
\sqrt{2\omega_{m}} and p ^ q m = i β Ο m
2 β ( b ^ m β β b ^ m ) \hat{p}{q{m}}={\rm i}\sqrt{\omega_{m}
2},(\hat{b}{m}^{\dagger}-\hat{b}{m}) to rewrite Eq. ( 2 ). Equations ( 2 ) and ( 4 ) equivalently describe the identical joint electronic-photonic problem.
The CBOA [ 16 ] applies the Born-Oppenheimer approximation to the photon displacement coordinates. It assumes that the photon displacement coordinates change slowly relative to the electrons, analogous to the nuclear coordinates. The CBOA is thus valid in the VSC regime, when the electronic states and the eigenstates of the cavity modes are well separated. The photon displacement coordinates can then be treated as classical parameters on the same footing as the nuclear coordinates π \mathbf{R} . The kinetic photon displacement term can then be separated from the electronic problem similarly to the nuclear kinetic energy. This yields an effective electronic Hamiltonian
All CBOA methods discussed in this paper are based on Eq. ( 5 ). The set of photon displacement coordinates πͺ β‘ { q m } \mathbf{q}\equiv{q_{m}} serves as variational parameters for the energy. They can be optimized together with the molecular orbitals: their initial values only start the iteration, and in every self-consistent field (SCF) iteration, each πͺ \mathbf{q} is updated to the stationary value of Eq. ( 8 ) for the current density, so that the converged density and πͺ opt \mathbf{q}^{\rm opt} are mutually consistent. Otherwise, πͺ \mathbf{q} is treated as a constant. When geometries are relaxed, it is also optimized alongside π \mathbf{R} . In contrast, the family of QED methods [ 34 ] is based on Eq. ( 2 ) and retains a fully quantized photon field in terms of Fock states.
The two Hamiltonians presented in Eq. ( 2 ) and Eq. ( 5 ) support formally different mean-field ansΓ€tze. In CBO-HF [ 2 ] , a single electronic Slater determinant,
is optimized self-consistently together with the photon displacement coordinates πͺ \mathbf{q} , which act as parameters of the energy functional. Consequently, no photonic wavefunction is obtained using CBO-HF . Because QED-HF [ 34 ] is based on Eq. ( 2 ), the electronic Slater determinant forms a tensor product with a coherent state of each cavity mode, that is, a displaced photon vacuum [ 43 , 34 ] :
with displacement amplitudes z m z_{m} optimized jointly with orbitals [ 34 , 20 ] . The coherent-state transformation removes the coupling of the photon field to the mean molecular dipole [ 34 , 44 ] ; therefore, the amplitudes vanish for a molecule without a static dipole component along π m \bm{\lambda}_{m} [ 34 , 2 ] .
Minimizing either energy with respect to its respective classical or quantum photon variable removes the net linear coupling between the photon and the mean dipole (zero-field condition). Both yield the same stationarity condition, up to a normalization factor relating the two:
Thus, each classical coordinate equals, mode by mode, the expectation value of the corresponding coherent state [ 34 , 2 , 1 ] . Evaluating the energy of each theory at this stationary point yields identical results up to a constant shift [ 1 ] :
This offset corresponds to the zero-point energy of the photo modes. For a given set of cavity frequencies { Ο m } {\omega_{m}} , this shift is a constant. Therefore, it does not alter gradients, vibrational frequencies, or densities. For converged E CBO-HF β ( πͺ opt ) E_{\rm\lx@glossaries@gls@link{acronym}{cbohf}{{{}}CBO-HF}}(\mathbf{q}^{\rm opt}) , the corresponding density and the nuclear gradient do not depend on the mode frequencies { Ο m } {\omega_{m}} , as shown below. The frequency dependence is confined to the photon coordinates and to derivative