Self-Consistent Field Methods for Excited States in Strong Magnetic Fields: a Comparison between Energy- and Variance-Based Approaches

Variance-based wave function optimization within the unrestricted doubly occupied configuration interaction framework: A half-projection treatment

The energy-variance-based optimization procedures have proven to be useful tools to describe N-electron spectra. However, the resulting wave functions usually present spin-contaminant contributions. The goal of this work is to reduce the spin contamination of the results arising from the unrestricted doubly occupied configuration interaction method in its energy variance minimization version [Alcoba et al., J. Chem. Phys. 160, 164107 (2024)]. We propose to incorporate the half-projection technique, which allows removing the spin components with even or odd spin quantum number of an approximate N-electron wave function, into the framework of the unrestricted doubly occupied configuration interaction treatment. This implementation can be carried out following several possible ways, whose results are analyzed in detail, in order to show the behavior of each procedure. Numerical determinations performed on selected strongly correlated N-electron systems, in ground and excited states, allow us to assess the most suitable procedure.

Computation of Magnetic Exchange Couplings in Photoexcited Systems Based on KS-DFT

Many efforts have been made in the study of optically excited spin-coupled molecules due to their appealing features for quantum information sciences. However, the characterization of the magnetic exchange couplings occurring from the photoexcitation is challenging experimentally. In this context, theoretical determinations play a critical role and must provide evaluations with rigorous and cost-effective strategies. This work presents a new approach to compute magnetic exchange couplings in photoexcited systems based on the recently generalized decomposition/recomposition method (David et al., Phys. Chem. Chem. Phys. 2024, 6, 8952-8964). This corresponds to the first application of KS-DFT in this context and offers both a completely general method and a powerful rationalization tool. This strategy is applied to mono- and biradical-based molecules recently synthesized by Kirk and co-workers.

Generalized Spin in the Variance-Based Wave Function Optimization Method within the Doubly Occupied Configuration Interaction Framework

In this work, we implement a generalized spin formulation of the doubly occupied configuration interaction methodology using the energy variance of the N-electron Hamiltonian. We perform the optimization of the N-electron wave functions and calculate their corresponding energies, using a unified variational treatment for ground and excited states based on the energy variance, which allows us to describe the entire energy spectra on an equal footing. We analyze the effects produced by the breakdown of the Ŝ2 and Ŝz symmetries in the spectra of model hydrogenic clusters in terms of energies and spin-related quantities, arising from the restricted, unrestricted, and generalized spin methods. The results are compared with other related methods as well as full configuration interaction.

Generalized spin σ-SCF method

We introduce a generalization of the σ-SCF method to approximate noncollinear spin ground and excited single-reference electronic states by minimizing the Hamiltonian variance. The new method is based on the σ-SCF method, originally proposed by Ye et al. [J. Chem. Phys. 147, 214104 (2017)], and provides a prescription to determine ground and excited noncollinear spin states on an equal footing. Our implementation was carried out utilizing an initial simulated annealing stage followed by a mean-field iterative self-consistent approach to simplify the cumbersome search introduced by generalizing the spin degrees of freedom. The simulated annealing stage ensures a broad exploration of the Hilbert space spanned by the generalized spin single-reference states with random complex element-wise rotations of the generalized density matrix elements in the simulated annealing stage. The mean-field iterative self-consistent stage employs an effective Fockian derived from the variance, which is utilized to converge tightly to the solutions. This process helps us to easily find complex spin structures, avoiding manipulating the initial guess. As proof-of-concept tests, we present results for Hn (n = 3-7) planar rings and polyhedral clusters with geometrical spin frustration. We show that most of these systems have noncollinear spin excited states that can be interpreted in terms of geometric spin frustration. These states are not directly targeted by energy minimization methods, which are meant to converge to the ground state. This stresses the capability of the σ-SCF methodology to find approximate noncollinear spin structures as mean-field excited states.

Unitary coupled-cluster for quantum computation of molecular properties in a strong magnetic field

In truncated coupled-cluster (CC) theories, non-variational and/or generally complex ground-state energies can occur. This is due to the non-Hermitian nature of the similarity transformed Hamiltonian matrix in combination with CC truncation. For chemical problems that deal with real-valued Hamiltonian matrices, complex CC energies rarely occur. However, for complex-valued Hamiltonian matrices, such as those that arise in the presence of strong magnetic fields, complex CC energies can be regularly observed unless certain symmetry conditions are fulfilled. Therefore, in the presence of magnetic fields, it is desirable to pursue CC methods that are guaranteed to give upper-bound, real-valued energies. In this work, we present the first application of unitary CC to chemical systems in a strong magnetic field. This is achieved utilizing the variational quantum eigensolver algorithm applied to the unitary coupled-cluster singles and doubles (UCCSD) method. We benchmark the method on the H2 molecule in a strong magnetic field and then calculate UCCSD energies for the H4 molecule as a function of both geometry and field angle. We show that while standard CCSD can yield generally complex energies that are not an upper-bound to the true energy, UCCSD always results in variational and real-valued energies. We also show that the imaginary components of the CCSD energy are largest in the strongly correlated region. Last, the UCCSD calculations capture a large percentage of the correlation energy.

TURBOMOLE: Today and Tomorrow

TURBOMOLE is a highly optimized software suite for large-scale quantum-chemical and materials science simulations of molecules, clusters, extended systems, and periodic solids. TURBOMOLE uses Gaussian basis sets and has been designed with robust and fast quantum-chemical applications in mind, ranging from homogeneous and heterogeneous catalysis to inorganic and organic chemistry and various types of spectroscopy, light-matter interactions, and biochemistry. This Perspective briefly surveys TURBOMOLE's functionality and highlights recent developments that have taken place between 2020 and 2023, comprising new electronic structure methods for molecules and solids, previously unavailable molecular properties, embedding, and molecular dynamics approaches. Select features under development are reviewed to illustrate the continuous growth of the program suite, including nuclear electronic orbital methods, Hartree-Fock-based adiabatic connection models, simplified time-dependent density functional theory, relativistic effects and magnetic properties, and multiscale modeling of optical properties.

Semiempirical Methods for Molecular Systems in Strong Magnetic Fields

A general scheme is presented to extend semiempirical methods to include the effects of arbitrary strength magnetic fields, while maintaining computational efficiency. The approach utilizes three main modifications; a London atomic orbital (LAO) basis set is introduced, field-dependent kinetic energy corrections are added to the model Hamiltonian, and spin-Zeeman interaction energy terms are included. The approach is applied to the widely available density-functional tight-binding method GFN1-xTB. Considering the basis set requirements for the kinetic energy corrections in a magnetic field leads to two variants: a single-basis approach GFN1-xTB-M0 and a dual-basis approach GFN1-xTB-M1. The LAO basis in the latter includes the appropriate nodal structure for an accurate representation of the kinetic energy corrections. The variants are assessed by benchmarking magnetizabilities and nuclear magnetic resonance shielding constants calculated using weak magnetic fields. Remarkably, the GFN1-xTB-M1 approach also exhibits excellent performance for strong fields, |B | ≤ 0.2B0 (B0 = 2.3505 × 105 T), recovering exotic features such as the para- to dia-magnetic transition in the BH molecule and the preferred electronic configuration, molecular conformation, and orientation of benzene. At stronger field strengths, |B | > 0.2B0, a degradation in the quality of the results is observed. The utility of GFN1-xTB-M1 is demonstrated by performing conformer searches in a range of field strengths for the cyclooctatetraene molecule, with GFN1-xTB-M1 capturing the transition from tub to planar conformations at high field, consistent with much more computationally demanding current-density functional theory calculations. Magnetically induced currents are also shown to be well described for the benzene and infinitene molecules, the latter demonstrating the flexibility and computational efficiency of the approach. The GFN1-xTB-M1 approach is a useful tool for the study of structure, conformation, and dynamics of large systems in magnetic fields at the semiempirical level as well as for preoptimization of molecular structure in ab initio calculations, enabling more efficient exploration of complex potential energy surfaces and reactivity in the presence of external fields.

Time-dependent nuclear-electronic orbital Hartree-Fock theory in a strong uniform magnetic field

In an ultrastrong magnetic field, with field strength B ≈ B₀ = 2.35 × 10⁵ T, molecular structure and dynamics differ strongly from that observed on the Earth. Within the Born-Oppenheimer (BO) approximation, for example, frequent (near) crossings of electronic energy surfaces are induced by the field, suggesting that nonadiabatic phenomena and processes may play a more important role in this mixed-field regime than in the weak-field regime on Earth. To understand the chemistry in the mixed regime, it therefore becomes important to explore non-BO methods. In this work, the nuclear-electronic orbital (NEO) method is employed to study protonic vibrational excitation energies in the presence of a strong magnetic field. The NEO generalized Hartree-Fock theory and time-dependent Hartree-Fock (TDHF) theory are derived and implemented, accounting for all terms that result as a consequence of the nonperturbative treatment of molecular systems in a magnetic field. The NEO results for HCN and FHF^- with clamped heavy nuclei are compared against the quadratic eigenvalue problem. Each molecule has three semi-classical modes owing to the hydrogen-two precession modes that are degenerate in the absence of a field and one stretching mode. The NEO-TDHF model is found to perform well; in particular, it automatically captures the screening effects of the electrons on the nuclei, which are quantified through the difference in energy of the precession modes.

Gold(i)-containing light-emitting molecules with an inverted singlet–triplet gap†

Delayed fluorescence from molecules with an inverted singlet–triplet gap (DFIST) is the consequence of the unusual reverse order of the lowest excited singlet (S₁) and triplet (T₁) states of thermally activated delayed fluorescence (TADF) emitters. Heptazine (1,3,4,6,7,9,9b-heptaazaphenalene) derivatives have an inverted singlet–triplet gap thanks to the combination of multiple resonance (MR) effects and a significant double excitation character. Here, we study computationally the effect of gold(i) metalation and coordination on the optical properties of heptazine (molecule 4) and the phosphine-functionalized 2,5,8-tris(dimethylphosphino)heptazine derivatives (molecules 1–3). Ab initio calculations at the approximate second-order coupled cluster (CC2) and extended multiconfigurational quasi degenerate perturbation theory at the second order (XMC-QDPT2) levels show that molecules 1–4 have an inverted singlet–triplet gap due to the alternating spatial localization of the electron and hole of the exciton in the heptazine core. A non-vanishing one-electron spin–orbit coupling operator matrix element between T₁ and and a fast S₁ ← T₁ intersystem crossing rate constant (k_ISC) calculated at the XMC-QDPT2(12,12) level of theory for molecule 4 suggest that this new family of complexes may be the first organometallic DFIST emitters reported.

Substitution with gold(i)-containing moieties results in non-vanishing oscillator strengths and spin–orbit coupling leading to fast intersystem crossing in light-emitting heptazine derivates with an inverted singlet–triplet gap.

An Embedded Fragment Method for Molecules in Strong Magnetic Fields

An extension of the embedded fragment method for calculations on molecular clusters is presented, which includes strong external magnetic fields. The approach is flexible, allowing for calculations at the Hartree-Fock, current-density-functional theory, Møller-Plesset perturbation theory, and coupled-cluster levels using London atomic orbitals. For systems consisting of discrete molecular subunits, calculations using London atomic orbitals can be performed in a computationally tractable manner for systems beyond the reach of conventional calculations, even those accelerated by resolution-of-the-identity or Cholesky decomposition methods. To assess the applicability of the approach, applications to water clusters are presented, showing how strong magnetic fields enhance binding within the clusters. However, our calculations suggest that, contrary to previous suggestions in the literature, this enhanced binding may not be directly attributable to strengthening of hydrogen bonding. Instead, these results suggest that this arises for larger field strengths as a response of the system to the presence of the external field, which induces a charge density build up between the monomer units. The approach is embarrassingly parallel and its computational tractability is demonstrated for clusters of up to 103 water molecules in triple-ζ basis sets, which would correspond to conventional calculations with more than 12 000 basis functions.

On the origin of the inverted singlet-triplet gap of the 5th generation light-emitting molecules

Excitation energies of the lowest singlet and triplet state of molecules whose first excited singlet state lies energetically below the first triplet state have been studied computationally at (time-dependent) density functional theory, coupled-cluster, and second-order multiconfiguration perturbation theory levels. The calculations at the ab initio levels show that the singlet-triplet gap is inverted as compared to the one expected from Hund's rule, whereas all density functionals yield the triplet state as the lowest excited state. Double excitations responsible for the inverted singlet-triplet gap have been identified. Employing the spin-flip and ΔSCF methods, singlet-triplet inversion was obtained at the density functional theory level for some of the studied molecules.

Molecular Dynamics of Linear Molecules in Strong Magnetic Fields

Molecular rotations and vibrations have been extensively studied by chemists for decades, both experimentally using spectroscopic methods and theoretically with the help of quantum chemistry. However, the theoretical investigation of molecular rotations and vibrations in strong magnetic fields requires computationally more demanding tools. As such, proper calculations of rotational and vibrational spectra were not feasible up until very recently. In this work, we present rotational and vibrational spectra for two small linear molecules, H2 and LiH, in strong magnetic fields. By treating the nuclei as classical particles, trajectories for rotations and vibrations are simulated from ab initio molecular dynamics. Born-Oppenheimer potential energy surfaces are calculated at the Hartree-Fock and MP2 levels of theory, using London atomic orbitals to ensure gauge origin invariance. For the calculation of nuclear trajectories, a highly efficient Tajima propagator is introduced, incorporating the Berry curvature tensor accounting for the screening of nuclear charges.

Revealing the exotic structure of molecules in strong magnetic fields

A novel implementation for the calculation of molecular gradients under strong magnetic fields is employed at the current-density functional theory level to optimize the geometries of molecular structures, which change significantly under these conditions. An analog of the ab initio random structure search is utilized to determine the ground-state equilibrium geometries for He n and CH n systems at high magnetic field strengths, revealing the most stable structures to be those in high-spin states with a planar geometry aligned perpendicular to the field. The electron and current densities for these systems have also been investigated to develop an explanation of chemical bonding in the strong field regime, providing an insight into the exotic chemistry present in these extreme environments.

Linear Response of Current-Dependent Density Functional Approximations in Magnetic Fields

This Letter outlines the steps and derivations that are necessary to apply density functional approximations that depend on the current and kinetic energy density rigorously within the framework of linear-response methods, including adiabatic time-dependent current density functional theory. This includes systems with a non-zero current density in the ground state. The necessary exchange-correlation kernel for these density functional approximations is derived, and the matrix elements are given explicitly. Due to the gauge variance of the kinetic energy density in an external magnetic field, having access to the proper current-dependent exchange-correlation kernel is necessary to recover gauge invariance for excited states. As a proof of principle application, the excited states of two small molecules in strong external magnetic fields are calculated using linear-response time-dependent current density functional theory. Finally, the implications of the derived current density-dependent exchange-correlation kernel for systems with strong spin-orbit coupling are discussed.

Energy Landscape of State-Specific Electronic Structure Theory

State-specific approximations can provide a more accurate representation of challenging electronic excitations by enabling relaxation of the electron density. While state-specific wave functions are known to be local minima or saddle points of the approximate energy, the global structure of the exact electronic energy remains largely unexplored. In this contribution, a geometric perspective on the exact electronic energy landscape is introduced. On the exact energy landscape, ground and excited states form stationary points constrained to the surface of a hypersphere, and the corresponding Hessian index increases at each excitation level. The connectivity between exact stationary points is investigated, and the square-magnitude of the exact energy gradient is shown to be directly proportional to the Hamiltonian variance. The minimal basis Hartree-Fock and excited-state mean-field representations of singlet H2 (STO-3G) are then used to explore how the exact energy landscape controls the existence and properties of state-specific approximations. In particular, approximate excited states correspond to constrained stationary points on the exact energy landscape, and their Hessian index also increases for higher energies. Finally, the properties of the exact energy are used to derive the structure of the variance optimization landscape and elucidate the challenges faced by variance optimization algorithms, including the presence of unphysical saddle points or maxima of the variance.