Search

We present the theory and implementation of a fully variational wave function-density functional theory (DFT) hybrid model, which is applicable to many cases of strong correlation. We denote this model as the multiconfigurational self-consistent on-top pair-density functional theory (MC-srPDFT) model. We have previously shown how the multiconfigurational short-range DFT (MC-srDFT) hybrid model can describe many multiconfigurational cases of any spin symmetry and also state-specific calculations on excited states [Hedegård et al., J. Chem. Phys. 148(21), 214103 (2018)]. However, the srDFT part of the MC-srDFT has some deficiencies that it shares with Kohn-Sham DFT; in particular, (1) self-interaction errors (albeit reduced because of the range separation), (2) that different MS states incorrectly become non-degenerate, and (3) that singlet and non-singlet states dissociating to the same open-shell fragments incorrectly lead to different electronic energies at dissociation. The model that we present in this paper corrects these deficiencies by introducing the on-top pair density as an auxiliary variable replacing the spin density. Unlike other models in the literature, our model is fully variational and employs a long-range version of the on-top pair density. The implementation is a second-order optimization algorithm ensuring robust convergence to both ground and excited states. We show how MC-srPDFT solves the mentioned challenges by sample calculations on the ground state singlet curve of H2, N2, and Cr2 and the lowest triplet curves for N2 and Cr2. Furthermore, the rotational barrier for ethene is investigated for the S0 and T1 states. The calculations show correct degeneracy between the singlet and triplet curves at dissociation and the results are invariant to the choice of the MS value for the triplet curves., (© 2025 Author(s). Published under an exclusive license by AIP Publishing.)

In this paper, benchmark results are presented on the calculation of vertical electronic excitation energies using a long-range second-order polarization propagator approximation (SOPPA) description with a short-range density functional theory description based on the Perdew-Burke-Ernzerhof (PBE) functional. The excitation energies are investigated for 132 singlet states and 71 triplet states across 28 medium-sized organic molecules. The results show that overall SOPPA-srPBE always performs better than PBE and that SOPPA-srPBE performs better than SOPPA for singlet states, but slightly worse than SOPPA for triplet states when CC3 results are the reference values., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.

In this paper, we present the theory and implementation of nuclear magnetic resonance shielding constants with gauge-including atomic orbitals for the hybrid multiconfigurational short-range density functional theory model. As a special case, this implementation also includes Hartree-Fock srDFT (HF-srDFT). Choosing a complete-active space (CAS) wave function as the multiconfigurational parameterization of the wave function, we investigate how well CAS-srDFT reproduces experimental trends of nuclear shielding constants compared to DFT and complete active space self-consistent field (CASSCF). Calculations on the nucleobases adenine and thymine show that CAS-srDFT performs on average the best of the tested methods, much better than CASSCF but only marginally better than HF-srDFT. The performance, compared to regular DFT, is similar when functionals containing exact exchange are used. We generally find that the inclusion of exact exchange is important for an accurate description of the shielding constants. In cases where no exact exchange is included, we observe that the HF- and CAS-srDFT often outperform regular DFT. For calculations on transition metal nuclei in organometallic compounds with significant static correlation, the CAS-srDFT method again outperforms CASSCF compared to experimental shielding constants, and the change from HF-srDFT is substantial. In conclusion, the static correlation posed by the metal complexes seems to be captured by CAS-srDFT, which is promising since this type of correlation is not well described by regular DFT.

The polarizable density embedding model is combined with the multiconfigurational self-consistent field and MC-srDFT electronic structure methods to calculate solvatochromic shifts of the n-π* absorption of acrolein and the π-π* absorption of the para -nitrophenolate anion in aqueous solution. Differences between linear-response (LR) and state-specific (SS) solvent shifts are analyzed by assessing the contributions of different terms in the solvent potential. This comparison shows that the differences are not only due to the intrinsically different response of LR and SS excitation energies to the polarizability of the environment but also due to a different response to the static part of the environment potential. These observations show that even in nonpolarizable environments, LR and SS calculations based on SCF (orbital optimization) methods do not necessarily agree on the spectral shift. The difference can be as large as, or even dominate, the difference due to dynamical polarization.

Based on self-consistent field (SCF) atomic mean-field (amf) quantities, we present two simple yet computationally efficient and numerically accurate matrix-algebraic approaches to correct both scalar-relativistic and spin-orbit two-electron picture-change effects (PCEs) arising within an exact two-component (X2C) Hamiltonian framework. Both approaches, dubbed amfX2C and e(xtended)amfX2C, allow us to uniquely tailor PCE corrections to mean-field models, viz. Hartree-Fock or Kohn-Sham DFT, in the latter case also avoiding the need for a point-wise calculation of exchange-correlation PCE corrections. We assess the numerical performance of these PCE correction models on spinor energies of group 18 (closed-shell) and group 16 (open-shell) diatomic molecules, achieving a consistent ≈10 -5 Hartree accuracy compared to reference four-component data. Additional tests include SCF calculations of molecular properties such as absolute contact density and contact density shifts in copernicium fluoride compounds (CnF n , n = 2,4,6), as well as equation-of-motion coupled-cluster calculations of x-ray core-ionization energies of 5d- and 6d-containing molecules, where we observe an excellent agreement with reference data. To conclude, we are confident that our (e)amfX2C PCE correction models constitute a fundamental milestone toward a universal and reliable relativistic two-component quantum-chemical approach, maintaining the accuracy of the parent four-component one at a fraction of its computational cost.

We present an implementation of time-dependent linear-response equations for strongly orthogonal geminal wave function models: the time-dependent generalized valence bond (TD-GVB) perfect-pairing theory and the antisymmetrized product of strongly orthogonal geminals. The geminal wave functions are optimized using a restricted-step second-order algorithm suitable for handling many geminals, and the linear-response equations are solved in an efficient way using a direct iterative approach. The wave function optimization algorithm features an original scheme to create initial orbitals for the geminal functions in a black-box fashion. The implementation is employed to examine the accuracy of the geminal linear response for singlet excitation energies of small and medium-sized molecules. In systems dominated by dynamic correlation, geminal models constitute only a minor improvement with respect to time-dependent Hartree-Fock. Compared to the linear-response complete active space self-consistent field, TD-GVB either misses or gives large errors for states dominated by double excitations.

In this paper, we report reimplementation of the core algorithms of relativistic coupled cluster theory aimed at modern heterogeneous high-performance computational infrastructures. The code is designed for parallel execution on many compute nodes with optional GPU coprocessing, accomplished via the new ExaTENSOR back end. The resulting ExaCorr module is primarily intended for calculations of molecules with one or more heavy elements, as relativistic effects on the electronic structure are included from the outset. In the current work, we thereby focus on exact two-component methods and demonstrate the accuracy and performance of the software. The module can be used as a stand-alone program requiring a set of molecular orbital coefficients as the starting point, but it is also interfaced to the DIRAC program that can be used to generate these. We therefore also briefly discuss an improvement of the parallel computing aspects of the relativistic self-consistent field algorithm of the DIRAC program.

The multi-configurational short-range (sr) density functional theory has been extended to the calculation of indirect spin-spin coupling constants (SSCCs) for nuclear magnetic resonance spectroscopy. The performance of the new method is compared to Kohn-Sham density functional theory and the ab initio complete active space self-consistent field for a selected set of molecules with good reference values. Two density functionals have been considered, the local density approximation srLDA and srPBE from the GGA class of functionals. All srDFT calculations are of Hartree-Fock-type HF-srDFT or complete active space-type CAS-srDFT. In all cases, the calculated SSCC values are of the same quality for srLDA and srPBE functionals, suggesting that one should use the computationally cost-effective srLDA functionals in applications. For all the calculated SSCCs in organic compounds, the best choice is HF-srDFT; the more expensive CAS-srDFT does not provide better values for these single-reference molecules. Fluorine is a challenge; in particular, the FF, FC, and FO couplings have much higher statistical errors than the rest. For SSCCs involving fluorine and a metal atom CAS-srDFT with singlet, generalized Tamm-Dancoff approximation is needed to get good SSCC values although the reference ground state is not a multi-reference case. For VF 6 -1 , all other considered models fail blatantly.

The Dalton Project provides a uniform platform access to the underlying full-fledged quantum chemistry codes Dalton and LSDalton as well as the PyFraME package for automatized fragmentation and parameterization of complex molecular environments. The platform is written in Python and defines a means for library communication and interaction. Intermediate data such as integrals are exposed to the platform and made accessible to the user in the form of NumPy arrays, and the resulting data are extracted, analyzed, and visualized. Complex computational protocols that may, for instance, arise due to a need for environment fragmentation and configuration-space sampling of biochemical systems are readily assisted by the platform. The platform is designed to host additional software libraries and will serve as a hub for future modular software development efforts in the distributed Dalton community.

DIRAC is a freely distributed general-purpose program system for one-, two-, and four-component relativistic molecular calculations at the level of Hartree-Fock, Kohn-Sham (including range-separated theory), multiconfigurational self-consistent-field, multireference configuration interaction, electron propagator, and various flavors of coupled cluster theory. At the self-consistent-field level, a highly original scheme, based on quaternion algebra, is implemented for the treatment of both spatial and time reversal symmetry. DIRAC features a very general module for the calculation of molecular properties that to a large extent may be defined by the user and further analyzed through a powerful visualization module. It allows for the inclusion of environmental effects through three different classes of increasingly sophisticated embedding approaches: the implicit solvation polarizable continuum model, the explicit polarizable embedding model, and the frozen density embedding model.

13C solid-state MAS NMR spectra of a series of paramagnetic metal acetylacetonate complexes; [VO(acac)2] (d1, S = ½), [V(acac)3] (d2, S = 1), [Ni(acac)2(H2O)2] (d8, S = 1), and [Cu(acac)2] (d9, S = ½), were assigned using modern NMR shielding calculations. This provided a reliable assignment of the chemical shifts and a qualitative insight into the hyperfine couplings. Our results show a reversal of the isotropic 13C shifts, δiso(13C), for CH3 and CO between the d1 and d2versus the d8 and d9 acetylacetonate complexes. The CH3 shifts change from about -150 ppm (d1,2) to roughly 1000 ppm (d8,9), whereas the CO shifts decrease from 800 ppm to about 150 ppm for d1,2 and d8,9, respectively. This was rationalized by comparison of total spin-density plots and computed contact couplings to those corresponding to singly occupied molecular orbitals (SOMOs). This revealed the interplay between spin delocalization of the SOMOs and spin polarization of the lower-energy MOs, influenced by both the molecular symmetry and the d-electron configuration. A large positive chemical shift results from spin delocalization and spin polarization acting in the same direction, whereas their cancellation corresponds to a small shift. The SOMO(s) for the d8 and d9 complexes are σ-like, implying spin-delocalization on the CH3 and CO groups of the acac ligand, cancelled only for CO by spin polarization. In contrast, the SOMOs of the d1 and d2 systems are π-like and a large CO-shift results from spin polarization, which accounts for the reversed assignment of δiso(13C) for CH3 and CO.

A theoretical framework for understanding molecular structures is crucial for the development of new technologies such as catalysts or solar cells. Apart from electronic excitation energies, however, only spectroscopic properties of molecules consisting of lighter elements can be computationally described at a high level of theory today since heavy elements require a relativistic framework, and thus far, most methods have only been derived in a non-relativistic framework. Important new technologies such as those mentioned above require molecules that contain heavier elements, and hence, there is a great need for the development of relativistic computational methods at a higher level of accuracy. Here, the Second-Order-Polarization-Propagator-Approximation (SOPPA), which has proven to be very successful in the non-relativistic case, is adapted to a relativistic framework. The equations for SOPPA are presented in their most general form, i.e., in a non-canonical spin-orbital basis, which can be reduced to the canonical case, and the expressions needed for a relativistic four-component SOPPA are obtained. The equations are one-index transformed, giving more compact expressions that correspond to those already available for the four-component RPA. The equations are ready for implementation in a four-component quantum chemistry program, which will allow both linear response properties and excitation energies to be calculated relativistically at the SOPPA level.

Linear response theory for the multiconfigurational short-range density functional theory (MC-srDFT) model is extended to triplet response with a singlet reference wave function. The triplet linear response equations for MC-srDFT are derived for a general hybrid srGGA functional and implemented in the Dalton program. Triplet excitation energies are benchmarked against the CC3 model of coupled cluster theory and the complete-active-space second-order perturbation theory using three different short-range functionals (srLDA, srPBE, and srPBE0), both with full linear response and employing the generalized Tamm-Dancoff approximation (gTDA). We find that using gTDA is required for obtaining reliable triplet excitations; for the CAS-srPBE model, the mean absolute deviation decreases from 0.40 eV to 0.26 eV, and for the CAS-srLDA model, it decreases from 0.29 eV to 0.21 eV. As expected, the CAS-srDFT model is found to be superior to the HF-srDFT model when analyzing the calculated triplet excitations for molecules in the benchmark set where increased static correlation is expected.

We present an electron-pairs-based method employing a generalized valence bond perfect-pairing (GVB-PP) ansatz that provides a uniformly accurate description of systems where various types of electron correlation play a role and the GVB-PP wave function is a suitable reference. In the proposed EERPA-GVB approach, a GVB-PP energy is amended by adding correlation among electron pairs. The latter is achieved by embedding single pairs or couples of pairs in the environment of the other electron fragments and separately accounting for intra- and interfragment correlation effects. For this purpose, we employ truncated extended random phase approximation equations. Application of EERPA-GVB to systems governed by both short-range (energy barriers) and long-range (molecular interactions) correlation effects proves the good accuracy of the method. Moreover, EERPA-GVB is shown to cure a notorious problem of uncorrelated electron-pair models, namely, spatial symmetry breaking in aromatic molecules, using the example of benzene. We have also successfully applied EERPA-GVB to a challenging problem of a phase transition of the boron chain system, where the correlation changes its character along the reaction path. The accuracy and versatility of EERPA-GVB are accompanied by its attractively low computational cost. By truncation of the extended RPA equations and consideration of only at most two-fragment correlation contributions, the cost of computing the EERPA correlation energy is reduced to scale only quadratically with the number of pairs of electrons.

Understanding the bonding, reactivity and electronic structure of actinides is lagging behind that of the rest of the periodic table. This can be partly explained by the challenges that one faces in experimental studies of such radioactive compounds and also by the need to properly account for relativistic effects in theoretical studies. A further challenge is the very complicated electronic structures encountered in actinide chemistry, as vividly illustrated by the naked diuranium molecule U 2 . Here we report a computational study of this emblematic molecule using state-of-the-art relativistic quantum chemical methods. Notably, the variational inclusion of spin-orbit interactions leads not only to a different electronic ground state, but also to a lower bond multiplicity compared with those in previous studies.

Many chemical systems cannot be described by quantum chemistry methods based on a single-reference wave function. Accurate predictions of energetic and spectroscopic properties require a delicate balance between describing the most important configurations (static correlation) and obtaining dynamical correlation efficiently. The former is most naturally done through a multiconfigurational (MC) wave function, whereas the latter can be done by, e.g., perturbation theory. We have employed a different strategy, namely, a hybrid between multiconfigurational wave functions and density-functional theory (DFT) based on range separation. The method is denoted by MC short-range DFT (MC-srDFT) and is more efficient than perturbative approaches as it capitalizes on the efficient treatment of the (short-range) dynamical correlation by DFT approximations. In turn, the method also improves DFT with standard approximations through the ability of multiconfigurational wave functions to recover large parts of the static correlation. Until now, our implementation was restricted to closed-shell systems, and to lift this restriction, we present here the generalization of MC-srDFT to open-shell cases. The additional terms required to treat open-shell systems are derived and implemented in the DALTON program. This new method for open-shell systems is illustrated on dioxygen and [Fe(H 2 O) 6 ] 3+ .

We present a derivation of linear response theory within polarizable embedding starting from a rigorous quantum-mechanical treatment of a composite system. To this aim, two different subsystem decompositions (symmetric and nonsymmetric) of the linear response function are introduced and the pole structures as well as residues of the individual terms are discussed. In addition to providing a thorough justification for the descriptions used in polarizable embedding models, this theoretical analysis clarifies which form of the response function to use and highlights complications in separating out subsystem contributions to molecular properties. The basic features of the presented expressions and various approximate forms are illustrated by their application to a composite model system.

Most chemistry, including chemistry where relativistic effects are important, occurs in an environment, and in many cases, this environment has a significant effect on the chemistry. In nonrelativistic quantum chemistry, a lot of progress has been achieved with respect to including environments such as a solvent or protein in the calculations, and now is the time to extend the possibilities for also doing this in relativistic quantum chemistry. The polarizable embedding (PE) model efficiently incorporates electrostatic effects of the environment by describing it as a collection of localized electric multipoles and polarizabilities obtained through quantum chemical calculations. In this article, we present the theory and implementation of four- and exact two-component Hamiltonians within a PE framework. We denote the methods the PE-4c-DFT and PE-X2C-DFT models. The models include a linear response formalism to calculate time-dependent (TD) properties: PE-TD-4c-DFT and PE-TD-X2C-DFT. With this first implementation, we calculate the PE-TD-4c-PBE0 excitation energies of the TcO 4 - and ReO 4 - ions in an explicit water solvent. This initial investigation focuses on the relative size of relativistic and solvent contributions to the excitation energies. The solvent effect is divided into an indirect solvent effect due to the structural perturbation of the XO 4 - ion and a direct electrostatic effect. The relativistic effects as well as both types of solvent effects are found to contribute to a shift in the excitation energies, but they do so to different extents depending on the ion and the electronic transition in question.