Adiabatic Connection Research Articles

Nonlocal Functionals Inspired by the Strongly Interacting Limit of DFT: Exact Constraints and Implementation.

Capturing strong correlation effects remains a key challenge for the development of improved exchange-correlation (XC) functionals in density functional theory. The recently proposed multiple radii functional (MRF) [J. Phys. Chem. Lett. 2017, 8, 2799; J. Chem. Theory Comput. 2019, 15, 3580] was designed to capture strong correlation effects seamlessly, as its mathematical structure draws from that of the exact XC functional in the limit of infinite correlations. The MRF functional provides a framework for building approximations along the density-fixed adiabatic connection, delivers accurate XC energy densities in the standard DFT gauge (same as that of the exact exchange energy density), and is free of one-electron self-interaction errors. To facilitate the development of XC functionals based on the MRF, we examine the behavior of the MRF functional when applied to uniform and scaled densities and consider how it can be made exact for the uniform electron gas. These theoretical insights are then used to build improved forms for the fluctuation function, an object that defines XC energy densities within the MRF framework. We also show how the MRF fluctuation function for physical correlation can be easily readjusted to accurately capture the XC functional in the limit of infinite correlations, demonstrating the versatility of MRF for building approximations for different correlation regimes. We describe the implementation of MRF using densities expanded on Gaussian basis sets, which improves the efficiency of previous grid-based MRF implementations.

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.

Strong-correlation density functionals made simple

Recent work on incorporating strong-correlation (sc) corrections into the scLH22t local hybrid functional [A. Wodyński and M. Kaupp, J. Chem. Theory Comput. 18, 6111–6123 (2022)] used a hybrid procedure, applying a strong-correlation factor derived from the reverse Becke–Roussel machinery of the KP16/B13 and B13 functionals to the nonlocal correlation term of a local hybrid functional. Here, we show that adiabatic-connection factors for strong-correlation-corrected local hybrids (scLHs) can be constructed in a simplified way based on a comparison of semi-local and exact exchange-energy densities only, without recourse to exchange-hole normalization. The simplified procedure is based on a comparative analysis of Becke’s B05 real-space treatment of nondynamical correlation and that in LHs, and it allows us to use, in principle, any semi-local exchange-energy density in the variable used to construct local adiabatic connections. The derivation of competitive scLHs is demonstrated based on either a modified Becke–Roussel or a simpler Perdew–Burke–Ernzerhof (PBE) energy density, leading to the scLH23t-mBR and scLH23t-tPBE functionals, which both exhibit low fractional spin errors while retaining good performance for weakly correlated situations. We also report preliminary attempts toward more detailed modeling of the local adiabatic connection, allowing a reduction of unphysical local maxima in spin-restricted bond-dissociation energy curves (scLH23t-mBR-P form). The simplified derivations of sc-factors reported here provide a basis for future constructions and straightforward implementation of exchange-correlation functionals that escape the zero-sum game between low self-interaction and static-correlation errors.

Seven useful questions in density functional theory

We explore a variety of unsolved problems in density functional theory, where mathematicians might prove useful. We give the background and context of the different problems, and why progress toward resolving them would help those doing computations using density functional theory. Subjects covered include the magnitude of the kinetic energy in Hartree–Fock calculations, the shape of adiabatic connection curves, using the constrained search with input densities, densities of states, the semiclassical expansion of energies, the tightness of Lieb–Oxford bounds, and how we decide the accuracy of an approximate density.

Correction to "Noncovalent Interactions from Models for the Møller-Plesset Adiabatic Connection".

ADVERTISEMENT RETURN TO ISSUEPREVAddition/CorrectionNEXTORIGINAL ARTICLEThis notice is a correctionCorrection to “Noncovalent Interactions from Models for the Møller–Plesset Adiabatic Connection”Timothy J. DaasTimothy J. DaasMore by Timothy J. Daashttps://orcid.org/0000-0002-4323-6908, Eduardo FabianoEduardo FabianoMore by Eduardo Fabianohttps://orcid.org/0000-0002-3990-669X, Fabio Della SalaFabio Della SalaMore by Fabio Della Salahttps://orcid.org/0000-0003-0940-8830, Paola Gori-GiorgiPaola Gori-GiorgiMore by Paola Gori-Giorgihttps://orcid.org/0000-0002-5952-1172, and Stefan Vuckovic*Stefan VuckovicMore by Stefan Vuckovichttps://orcid.org/0000-0002-0768-9176Cite this: J. Phys. Chem. Lett. 2023, 14, 6, 1478Publication Date (Web):February 6, 2023Publication History Received30 January 2023Published online6 February 2023Published inissue 16 February 2023https://doi.org/10.1021/acs.jpclett.3c00264Copyright © 2023 The Authors. Published by American Chemical SocietyRIGHTS & PERMISSIONSACS AuthorChoiceCC: Creative CommonsBY: Credit must be given to the creatorNC: Only noncommercial uses of the work are permittedND: No derivatives or adaptations of the work are permittedArticle Views85Altmetric-Citations-LEARN ABOUT THESE METRICSArticle Views are the COUNTER-compliant sum of full text article downloads since November 2008 (both PDF and HTML) across all institutions and individuals. These metrics are regularly updated to reflect usage leading up to the last few days.Citations are the number of other articles citing this article, calculated by Crossref and updated daily. Find more information about Crossref citation counts.The Altmetric Attention Score is a quantitative measure of the attention that a research article has received online. Clicking on the donut icon will load a page at altmetric.com with additional details about the score and the social media presence for the given article. Find more information on the Altmetric Attention Score and how the score is calculated. Share Add toView InAdd Full Text with ReferenceAdd Description ExportRISCitationCitation and abstractCitation and referencesMore Options Share onFacebookTwitterWechatLinked InReddit PDF (920 KB) Get e-Alerts Get e-Alerts

Toward more accurate adiabatic connection approach for multireference wavefunctions.

A multiconfigurational adiabatic connection (AC) formalism is an attractive approach to compute the dynamic correlation within the complete active space self-consistent field and density matrix renormalization group (DMRG) models. Practical realizations of AC have been based on two approximations: (i) fixing one- and two-electron reduced density matrices (1- and 2-RDMs) at the zero-coupling constant limit and (ii) extended random phase approximation (ERPA). This work investigates the effect of removing the "fixed-RDM" approximation in AC. The analysis is carried out for two electronic Hamiltonian partitionings: the group product function- and the Dyall Hamiltonians. Exact reference AC integrands are generated from the DMRG full configuration interaction solver. Two AC models are investigated, employing either exact 1- and 2-RDMs or their second-order expansions in the coupling constant in the ERPA equations. Calculations for model molecules indicate that lifting the fixed-RDM approximation is a viable way toward improving the accuracy of existing AC approximations.

Geometries and vibrational frequencies with Kohn-Sham methods using σ-functionals for the correlation energy.

Recently, Kohn-Sham (KS) methods with new correlation functionals, called σ-functionals, have been introduced. Technically, σ-functionals are closely related to the well-known random phase approximation (RPA); formally, σ-functionals are rooted in perturbation theory along the adiabatic connection. If employed in a post-self-consistent field manner in a Gaussian basis set framework, then, σ-functional methods are computationally very efficient. Moreover, for main group chemistry, σ-functionals are highly accurate and can compete with high-level wave-function methods. For reaction and transition state energies, e.g., chemical accuracy of 1 kcal/mol is reached. Here, we show how to calculate first derivatives of the total energy with respect to nuclear coordinates for methods using σ-functionals and then carry out geometry optimizations for test sets of main group molecules, transition metal compounds, and non-covalently bonded systems. For main group molecules, we additionally calculate vibrational frequencies. σ-Functional methods are found to yield very accurate geometries and vibrational frequencies for main group molecules superior not only to those from conventional KS methods but also to those from RPA methods. For geometries of transition metal compounds, not surprisingly, best geometries are found for RPA methods, while σ-functional methods yield somewhat less good results. This is attributed to the fact that in the optimization of σ-functionals, transition metal compounds could not be represented well due to the lack of reliable reference data. For non-covalently bonded systems, σ-functionals yield geometries of the same quality as the RPA or as conventional KS schemes combined with dispersion corrections.

Majorana zero modes in fermionic wires coupled by Aharonov-Bohm cages

We devise a number-conserving scheme for the realization of Majorana zero modes in an interacting fermionic ladder coupled by Aharonov-Bohm cages. The latter provide an efficient mechanism to cancel single-particle hopping by destructive interference. The crucial parity symmetry in each wire is thus encoded in the geometry of the setup, in particular, its translation invariance. A generic nearest-neighbor interaction generates the desired correlated hopping of pairs. We exhibit the presence of an extended topological region in parameter space, first in a simplified effective model via bosonization techniques, and subsequently in a larger parameter regime with matrix-product-states numerical simulations. We demonstrate the adiabatic connection to previous models, including exactly solvable ones, and we briefly comment on possible experimental realizations in synthetic quantum platforms, like cold atomic samples.

Automatic Generation of Local Vibrational Mode Parameters: From Small to Large Molecules and QM/MM Systems.

LModeAGen, a new protocol for the automatic determination of a nonredundant, complete set of local vibrational modes is reported, which is based on chemical graph concepts. Whereas local mode properties can be calculated for a selection of parameters targeting specific local modes of interest, a complete set of nonredundant local mode parameters is requested for the adiabatic connection scheme (ACS), relating each local vibrational mode with a normal mode counterpart, and for the decomposition of normal modes (CNM) in terms of local mode contributions, a unique way to analyze vibrational spectra. So far, nonredundant parameter sets have been generated manually following chemical intuition or from a set of redundant parameters in a trial-and-error fashion, which has hampered the study of larger systems with hundreds of parameters. LModeAGen was successfully applied for a test set of 11 systems, ranging from small molecules to the large QM (>100 atoms) subsystem of carbomonoxy-neuroglobin protein, described with a hybrid QM/MM method. The ωB97X-D/aug-cc-pVDZ, M06L/def2-TZVP, and QM/MM ωB97X-D/6-31G(d,p)/AMBER model chemistries were adopted for the description of the molecules in the test set. Our new protocol is an important step forward for a routine ACS and CNM analysis of the vibrational spectra of complex and large systems with hundreds of atoms, providing new access to important encoded electronic structure information.

Local Hybrid Functional Applicable to Weakly and Strongly Correlated Systems.

The recent idea (Wodyński, A.; Arbuznikov, A. V.; Kaupp M. J. Chem. Phys. 2021, 155, 144101) to augment local hybrid functionals by a strong-correlation (sc) factor obtained from the adiabatic connection in the spirit of the KP16 model has been extended and applied to generate the accurate sc-corrected local hybrid functional scLH22t. By damping small values of the ratio between nondynamical and dynamical correlation entering the correction factor, it has become possible to avoid double counting of nondynamical correlation for weakly correlated situations and thereby preserve the excellent accuracy of the underlying LH20t local hybrid for such cases almost perfectly. On the other hand, scLH22t improves substantially over LH20t in reducing fractional-spin errors (FSEs), in providing improved spin-restricted bond dissociation curves, and in treating some typical systems with multireference character. The obtained FSEs are similar to those of the KP16/B13 model and slightly larger than for B13, but performance for weakly correlated systems is better than for these two related methods, which are also difficult to use self-consistently. The recent DM21 functional based on the training of a deep neural network still performs somewhat better than scLH22t but allows no physical insights into the origins of reduced FSEs. Examination of local mixing functions (LMFs) for the corrected scLH22t and uncorrected LH20t functionals provides further insights: in weakly correlated situations, the LMF remains essentially unchanged. Strong-correlation effects manifest in a reduction of the LMF values in certain regions of space, even to the extent of producing negative LMF values. It is suggested that this is the mechanism by which also DM21, which may be viewed as a range-separated local hybrid, is able to reduce FSEs.

Assessment of the Second-Order Statically Screened Exchange Correction to the Random Phase Approximation for Correlation Energies.

With increasing interelectronic distance, the screening of the electron–electron interaction by the presence of other electrons becomes the dominant source of electron correlation. This effect is described by the random phase approximation (RPA) which is therefore a promising method for the calculation of weak interactions. The success of the RPA relies on the cancellation of errors, which can be traced back to the violation of the crossing symmetry of the 4-point vertex, leading to strongly overestimated total correlation energies. By the addition of second-order screened exchange (SOSEX) to the correlation energy, this issue is substantially reduced. In the adiabatic connection (AC) SOSEX formalism, one of the two electron–electron interaction lines in the second-order exchange term is dynamically screened (SOSEX(W, vc)). A related SOSEX expression in which both electron–electron interaction lines are statically screened (SOSEX(W(0), W(0))) is obtained from the G3W2 contribution to the electronic self-energy. In contrast to SOSEX(W, vc), the evaluation of this correlation energy expression does not require an expensive numerical frequency integration and is therefore advantageous from a computational perspective. We compare the accuracy of the statically screened variant to RPA and RPA+SOSEX(W, vc) for a wide range of chemical reactions. While both methods fail for barrier heights, SOSEX(W(0), W(0)) agrees very well with SOSEX(W, vc) for charged excitations and noncovalent interactions where they lead to major improvements over RPA.

Critical topology and pressure-induced superconductivity in the van der Waals compound AuTe2Br

The study on quantum spin Hall effect and topological insulators formed the prologue to the surge of research activities in topological materials in the past decade. Compared to intricately engineered quantum wells, three-dimensional weak topological insulators provide a natural route to the quantum spin Hall effect, due to the adiabatic connection between them and a stack of quantum spin Hall insulators, and the convenience in exfoliation of samples associated with their van der Waals-type structure. Despite these advantages, both theoretical prediction and experimental identification of weak topological insulators remain scarce. Here, based on first-principles calculations, we show that AuTe2Br locates at the boundary between a strong and a weak topological semimetal state. We identify the key structural parameter that dictates the traversal of the topological transition, which can be easily realized in experiments. More interestingly, the critical topology of AuTe2Br persists up to an applied pressure of ~15.4 GPa before a structural phase transition accompanied by a change of electronic topology and the onset of superconductivity. Our results establish AuTe2Br as a new candidate for an effective tuning between weak and strong topological phases in a single material, with the potential to realize various other topological phases of matter.

Scaled σ-functionals for the Kohn-Sham correlation energy with scaling functions from the homogeneous electron gas.

The recently introduced σ-functionals constitute a new type of functionals for the Kohn-Sham (KS) correlation energy. σ-Functionals are based on the adiabatic-connection fluctuation-dissipation theorem, are computationally closely related to the well-known direct random phase approximation (dRPA), and are formally rooted in many-body perturbation theory along the adiabatic connection. In σ-functionals, the function of the eigenvalues σ of the Kohn-Sham response matrix that enters the coupling constant and frequency integration in the dRPA is replaced by another function optimized with the help of reference sets of atomization, reaction, transition state, and non-covalent interaction energies. σ-Functionals are highly accurate and yield chemical accuracy of 1 kcal/mol in reaction or transition state energies, in main group chemistry. A shortcoming of σ-functionals is their inability to accurately describe processes involving a change of the electron number, such as ionizations or electron attachments. This problem is attributed to unphysical self-interactions caused by the neglect of the exchange kernel in the dRPA and σ-functionals. Here, we tackle this problem by introducing a frequency- and σ-dependent scaling of the eigenvalues σ of the KS response function that models the effect of the exchange kernel. The scaling factors are determined with the help of the homogeneous electron gas. The resulting scaled σ-functionals retain the accuracy of their unscaled parent functionals but in addition yield very accurate ionization potentials and electron affinities. Moreover, atomization and total energies are found to be exceptionally accurate. Scaled σ-functionals are computationally highly efficient like their unscaled counterparts.

Self-Consistent Implementation of Kohn-Sham Adiabatic Connection Models with Improved Treatment of the Strong-Interaction Limit.

Adiabatic connection models (ACMs), which interpolate between the limits of weak and strong interaction, are powerful tools to build accurate exchange–correlation functionals. If the exact weak-interaction expansion from the second-order perturbation theory is included, a self-consistent implementation of these functionals is challenging and still absent in the literature. In this work, we fill this gap by presenting a fully self-consistent-field (SCF) implementation of some popular ACM functionals. While using second-order perturbation theory at weak interactions, we have also introduced new generalized gradient approximations (GGAs), beyond the usual point-charge-plus-continuum model, for the first two leading terms at strong interactions, which are crucial to ensure robustness and reliability. We then assess the SCF–ACM functionals for molecular systems and for prototypical strong-correlation problems. We find that they perform well for both the total energy and the electronic density and that the impact of SCF orbitals is directly connected to the accuracy of the ACM functional form. For the H2 dissociation, the SCF–ACM functionals yield significant improvements with respect to standard functionals also thanks to the use of the new GGAs for the strong-coupling functionals.

Reduced density matrix functional theory from an ab initio seniority-zero wave function: Exact and approximate formulations along adiabatic connection paths

Currently, there is a growing interest in the development of a new hierarchy of methods based on the concept of seniority, which has been introduced quite recently in quantum chemistry. Despite the enormous potential of these methods, the accurate description of both dynamical and static correlation effects within a single and in-principle-exact approach remains a challenge. In this work, we propose an alternative formulation of reduced density-matrix functional theory (RDMFT) where the (one-electron reduced) density matrix is mapped onto an ab initio seniority-zero wave function. In this theory, the exact natural orbitals and their occupancies are determined self-consistently from an effective seniority-zero calculation. The latter involves a universal higher-seniority density matrix functional for which an adiabatic connection (AC) formula is derived and implemented under specific constraints that are related to the density matrix. The pronounced curvature of the (constrained) AC integrand, which is numerically observed in prototypical hydrogen chains and the Helium dimer, indicates that a description of higher-seniority correlations within second-order perturbation theory is inadequate in this context. Applying multiple linear interpolations along the AC or connecting second-order perturbation theory to a full-seniority treatment via Pad\'{e} approximants are better strategies. Such information is expected to serve as a guide in the future design of higher-seniority density-matrix functional approximations.

Improving the Accuracy of Modelling CO2 Electroreduction on Copper Using Many-Body Perturbation Theory.

Copper (Cu) remains the most important metal catalyst for the carbon dioxide reduction reaction (CO2 RR) into C2 products. Due to limited evidence from in situ experiments, mechanistic studies are often performed in the framework of density functional theory (DFT), using functionals at the generalized gradient approximation (GGA) level, which have fundamental difficulties to correctly describe CO adsorption and surface stability. We employ the adiabatic connection fluctuation dissipation theorem within the random phase approximation (RPA), in combination with the linearized Poisson-Boltzmann equation to describe solvation effects, to investigate the mechanism of CO2 RR on the Cu(100) facet. Qualitatively different from the DFT-GGA results, RPA results propose the formation of *OCCHO as the potential determining step towards C2 products. The results suggest that it is important to use more accurate methods like RPA when modeling reactions involving multiple CO-related species like CO2 RR.

Improving the Accuracy of Modelling CO2 Electroreduction on Copper Using Many‐Body Perturbation Theory

AbstractCopper (Cu) remains the most important metal catalyst for the carbon dioxide reduction reaction (CO2RR) into C2 products. Due to limited evidence from in situ experiments, mechanistic studies are often performed in the framework of density functional theory (DFT), using functionals at the generalized gradient approximation (GGA) level, which have fundamental difficulties to correctly describe CO adsorption and surface stability. We employ the adiabatic connection fluctuation dissipation theorem within the random phase approximation (RPA), in combination with the linearized Poisson–Boltzmann equation to describe solvation effects, to investigate the mechanism of CO2RR on the Cu(100) facet. Qualitatively different from the DFT‐GGA results, RPA results propose the formation of *OCCHO as the potential determining step towards C2 products. The results suggest that it is important to use more accurate methods like RPA when modeling reactions involving multiple CO‐related species like CO2RR.

Recent Advances in Cartesian-Grid DFT in Atoms and Molecules.

In the past several decades, density functional theory (DFT) has evolved as a leading player across a dazzling variety of fields, from organic chemistry to condensed matter physics. The simple conceptual framework and computational elegance are the underlying driver for this. This article reviews some of the recent developments that have taken place in our laboratory in the past 5 years. Efforts are made to validate a viable alternative for DFT calculations for small to medium systems through a Cartesian coordinate grid- (CCG-) based pseudopotential Kohn–Sham (KS) DFT framework using LCAO-MO ansatz. In order to legitimize its suitability and efficacy, at first, electric response properties, such as dipole moment ( μ ), static dipole polarizability ( α ), and first hyperpolarizability ( β ), are calculated. Next, we present a purely numerical approach in CCG for proficient computation of exact exchange density contribution in certain types of orbital-dependent density functionals. A Fourier convolution theorem combined with a range-separated Coulomb interaction kernel is invoked. This takes motivation from a semi-numerical algorithm, where the rate-deciding factor is the evaluation of electrostatic potential. Its success further leads to a systematic self-consistent approach from first principles, which is desirable in the development of optimally tuned range-separated hybrid and hyper functionals. Next, we discuss a simple, alternative time-independent DFT procedure, for computation of single-particle excitation energies, by means of “adiabatic connection theorem” and virial theorem. Optical gaps in organic chromophores, dyes, linear/non-linear PAHs, and charge transfer complexes are faithfully reproduced. In short, CCG-DFT is shown to be a successful route for various practical applications in electronic systems.

Dispersion Interactions in Exciton-Localized States. Theory and Applications to π-π* and n-π* Excited States.

We address the problem of intermolecular interaction energy calculations in molecular complexes with localized excitons. Our focus is on the correct representation of the dispersion energy. We derive an extended Casimir-Polder formula for direct computation of this contribution through second order in the intermolecular interaction operator V̂. An alternative formula, accurate to infinite order in V̂, is derived within the framework of the adiabatic connection (AC) theory. We also propose a new parametrization of the VV10 nonlocal correlation density functional, so that it corrects the CASSCF energy for the dispersion contribution and can be applied to excited-state complexes. A numerical investigation is carried out for benzene, pyridine, and peptide complexes with the local exciton corresponding to the lowest π–π* or n– π* states. The extended Casimir-Polder formula is implemented in the framework of multiconfigurational symmetry-adapted perturbation theory, SAPT(MC). A SAPT(MC) analysis shows that the creation of a localized exciton affects mostly the electrostatic component of the interaction energy of investigated complexes. Nevertheless, the changes in Pauli repulsion and dispersion energies cannot be neglected. We verify the performance of several perturbation- and AC-based methods. Best results are obtained with a range-separated variant of an approximate AC approach employing extended random phase approximation and CASSCF wave functions.

Efficient Adiabatic Connection Approach for Strongly Correlated Systems: Application to Singlet-Triplet Gaps of Biradicals.

Strong electron correlation can be captured with multireference wave function methods, but an accurate description of the electronic structure requires accounting for the dynamic correlation, which they miss. In this work, a new approach for the correlation energy based on the adiabatic connection (AC) is proposed. The ACn method accounts for terms up to order n in the coupling constant, and it is size-consistent and free from instabilities. It employs the multireference random phase approximation and the Cholesky decomposition technique, leading to a computational cost growing with the fifth power of the system size. Because of the dependence on only one- and two-electron reduced density matrices, ACn is more efficient than existing ab initio multireference dynamic correlation methods. ACn affords excellent results for singlet–triplet gaps of challenging organic biradicals. The development presented in this work opens new perspectives for accurate calculations of systems with dozens of strongly correlated electrons.