Frozen-Density Embedding for Including Environmental Effects in the Dirac-Kohn–Sham Theory: An Implementation Based on Density Fitting and Prototyping Techniques

Frozen density embedding (FDE) represents an embedding scheme in which environmental effects are included from first-principles calculations by considering the surrounding system explicitly by means of its electron density. In the present paper, we extend the full four-component relativistic Dirac–Kohn–Sham (DKS) method, as implemented in the BERTHA code, to include environmental and confinement effects with the FDE scheme (DKS-in-DFT FDE). The implementation, based on the auxiliary density fitting techniques, has been enormously facilitated by BERTHA’s python API (PyBERTHA), which facilitates the interoperability with other FDE implementations available through the PyADF framework. The accuracy and numerical stability of this new implementation, also using different auxiliary fitting basis sets, has been demonstrated on the simple NH3–H2O system, in comparison with a reference nonrelativistic implementation. The computational performance has been evaluated on a series of gold clusters (Aun, with n = 2, 4, 8) embedded into an increasing number of water molecules (5, 10, 20, 40, and 80 water molecules). We found that the procedure scales approximately linearly both with the size of the frozen surrounding environment (consistent with the underpinnings of the FDE approach) and with the size of the active system (in line with the use of density fitting). Finally, we applied the code to a series of heavy (Rn) and super-heavy elements (Cn, Fl, Og) embedded in a C60 cage to explore the confinement effect induced by C60 on their electronic structure. We compare the results from our simulations, with respect to more-approximate models employed in the atomic physics literature. Our results indicate that the specific interactions described by FDE are able to improve upon the cruder approximations currently employed, and, thus, they provide a basis from which to generate more-realistic radial potentials for confined atoms.

[1]  A. Gomes,et al.  Environment Effects on X-Ray Absorption Spectra With Quantum Embedded Real-Time Time-Dependent Density Functional Theory Approaches , 2021, Frontiers in Chemistry.

[2]  Daniel S. Levine,et al.  Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package , 2021, The Journal of chemical physics.

[3]  R. Schlögl,et al.  Core level photoelectron spectroscopy of heterogeneous reactions at liquid-vapor interfaces: Current status, challenges, and prospects. , 2021, The Journal of chemical physics.

[4]  S. Höfener.  The KOALA program: Wavefunction frozen‐density embedding , 2020 .

[5]  C. Hennig,et al.  The missing pieces of the PuO2 nanoparticle puzzle. , 2020, Nanoscale.

[6]  P. Pyykkö,et al.  The periodic table and the physics that drives it , 2020, Nature Reviews Chemistry.

[7]  L. Visscher,et al.  Environmental Effects with Frozen-Density Embedding in Real-Time Time-Dependent Density Functional Theory Using Localized Basis Functions , 2020, Journal of chemical theory and computation.

[8]  K. Ruud,et al.  ReSpect: Relativistic spectroscopy DFT program package. , 2020, The Journal of chemical physics.

[9]  H. Quiney,et al.  BERTHA: Implementation of a four-component Dirac-Kohn-Sham relativistic framework. , 2020, The Journal of chemical physics.

[10]  H. Quiney,et al.  PyBERTHART: A relativistic real-time four-component TDDFT implementation using prototyping techniques based on Python. , 2020, Journal of chemical theory and computation.

[11]  M. Ratner,et al.  Quantum embedding for material chemistry based on domain separation and open subsystems , 2020 .

[12]  Jógvan Magnus Haugaard Olsen,et al.  The DIRAC code for relativistic molecular calculations. , 2020, The Journal of chemical physics.

[13]  Ning Zhang,et al.  BDF: A relativistic electronic structure program package. , 2020, The Journal of chemical physics.

[14]  C. Jacob,et al.  Frozen‐density embedding‐based many‐body expansions , 2020, International Journal of Quantum Chemistry.

[15]  M. Ratner,et al.  Embedding Methods for Quantum Chemistry: Applications from Materials to Life Sciences. , 2020, Journal of the American Chemical Society.

[16]  A. Gomes,et al.  Investigating solvent effects on the magnetic properties of molybdate ions (MoO$_{4}^{2-}$) with relativistic embedding , 2019, 1912.06192.

[17]  V. Vallet,et al.  Investigating solvent effects on the magnetic properties of molybdate ions () with relativistic embedding , 2019 .

[18]  M. Pavanello,et al.  Nonlocal Subsystem Density Functional Theory. , 2019, The journal of physical chemistry letters.

[19]  J. Neugebauer,et al.  Inter-subsystem charge-transfer excitations in exact subsystem time-dependent density-functional theory. , 2019, The Journal of chemical physics.

[20]  Johannes L. Schönberger,et al.  SciPy 1.0: fundamental algorithms for scientific computing in Python , 2019, Nature Methods.

[21]  J. Neugebauer,et al.  Exact subsystem time-dependent density-functional theory. , 2019, The Journal of chemical physics.

[22]  Frederick R Manby,et al.  Projection-Based Wavefunction-in-DFT Embedding. , 2019, Accounts of chemical research.

[23]  F. Della Sala,et al.  Performance of Semilocal Kinetic Energy Functionals for Orbital-Free Density Functional Theory. , 2019, Journal of chemical theory and computation.

[24]  W. Nazarewicz,et al.  Colloquium : Superheavy elements: Oganesson and beyond , 2019, Reviews of Modern Physics.

[25]  F. Réal,et al.  Predictive Simulations of Ionization Energies of Solvated Halide Ions with Relativistic Embedded Equation of Motion Coupled Cluster Theory. , 2018, Physical review letters.

[26]  A. Wasserman,et al.  Constructing a non-additive non-interacting kinetic energy functional approximation for covalent bonds from exact conditions. , 2018, The Journal of chemical physics.

[27]  Johannes Neugebauer,et al.  Serenity: A subsystem quantum chemistry program , 2018, J. Comput. Chem..

[28]  F. Della Sala,et al.  Semilocal Pauli-Gaussian Kinetic Functionals for Orbital-Free Density Functional Theory Calculations of Solids. , 2018, The journal of physical chemistry letters.

[29]  E. Vauthey,et al.  Ultrafast Elementary Photochemical Processes of Organic Molecules in Liquid Solution. , 2017, Chemical reviews.

[30]  R. Distasio,et al.  eQE: An open‐source density functional embedding theory code for the condensed phase , 2017 .

[31]  Qiming Sun,et al.  Quantum Embedding Theories. , 2016, Accounts of chemical research.

[32]  Radovan Bast,et al.  On the calculation of second-order magnetic properties using subsystem approaches in a relativistic framework. , 2016, Physical chemistry chemical physics : PCCP.

[33]  György G. Ferenczy,et al.  Exact density functional and wave function embedding schemes based on orbital localization , 2016 .

[34]  J. Neugebauer,et al.  Part and whole in wavefunction/DFT embedding , 2015, Theoretical Chemistry Accounts.

[35]  S. Rampino,et al.  Gold-superheavy-element interaction in diatomics and cluster adducts: A combined four-component Dirac-Kohn-Sham/charge-displacement study. , 2015, The Journal of chemical physics.

[36]  Tomasz A Wesolowski,et al.  Frozen-Density Embedding Strategy for Multilevel Simulations of Electronic Structure. , 2015, Chemical reviews.

[37]  D. Ceresoli,et al.  Subsystem real-time time dependent density functional theory. , 2015, The Journal of chemical physics.

[38]  Sergio Rampino,et al.  Full Parallel Implementation of an All-Electron Four-Component Dirac-Kohn-Sham Program. , 2014, Journal of chemical theory and computation.

[39]  A. W. Götz,et al.  Calculation of nuclear spin-spin coupling constants using frozen density embedding. , 2014, The Journal of chemical physics.

[40]  Lucas Visscher,et al.  Solvatochromic shifts from coupled-cluster theory embedded in density functional theory. , 2013, The Journal of chemical physics.

[41]  Lucas Visscher,et al.  Towards systematically improvable models for actinides in condensed phase: the electronic spectrum of uranyl in Cs2UO2Cl4 as a test case. , 2013, Physical chemistry chemical physics : PCCP.

[42]  K. Maher,et al.  Environmental speciation of actinides. , 2013, Inorganic chemistry.

[43]  F. Della Sala,et al.  Semilocal and hybrid density embedding calculations of ground-state charge-transfer complexes. , 2013, The Journal of chemical physics.

[44]  A. Türler,et al.  Advances in the production and chemistry of the heaviest elements. , 2013, Chemical reviews.

[45]  Frederick R. Manby,et al.  A Simple, Exact Density-Functional-Theory Embedding Scheme , 2012, Journal of chemical theory and computation.

[46]  Christoph R. Jacob,et al.  Quantum-chemical embedding methods for treating local electronic excitations in complex chemical systems , 2012 .

[47]  P. Pyykkö.  Relativistic effects in chemistry: more common than you thought. , 2012, Annual review of physical chemistry.

[48]  Lucas Visscher,et al.  Molecular properties via a subsystem density functional theory formulation: a common framework for electronic embedding. , 2012, The Journal of chemical physics.

[49]  Andreas W. Götz,et al.  PyADF — A scripting framework for multiscale quantum chemistry , 2011, J. Comput. Chem..

[50]  A. Wasserman,et al.  Molecular binding energies from partition density functional theory. , 2011, The Journal of chemical physics.

[51]  H. Quiney,et al.  Recent advances and perspectives in four-component Dirac-Kohn-Sham calculations. , 2011, Physical chemistry chemical physics : PCCP.

[52]  M. N. Guimarães,et al.  A study of the electron structure of endohedrally confined atoms using a model potential , 2011 .

[53]  F. Della Sala,et al.  Frozen density embedding with hybrid functionals. , 2010, The Journal of chemical physics.

[54]  Frederick R Manby,et al.  Exact nonadditive kinetic potentials for embedded density functional theory. , 2010, The Journal of chemical physics.

[55]  Kenneth Ruud,et al.  Arbitrary-Order Density Functional Response Theory from Automatic Differentiation. , 2010, Journal of chemical theory and computation.

[56]  Samuel Fux,et al.  Accurate frozen-density embedding potentials as a first step towards a subsystem description of covalent bonds. , 2010, The Journal of chemical physics.

[57]  Sergio Rampino,et al.  Efficient Parallel All-Electron Four-Component Dirac-Kohn-Sham Program Using a Distributed Matrix Approach II. , 2010, Journal of chemical theory and computation.

[58]  Lucas Visscher,et al.  Performance of Kinetic Energy Functionals for Interaction Energies in a Subsystem Formulation of Density Functional Theory. , 2009, Journal of chemical theory and computation.

[59]  K. Dyall.  Relativistic double-zeta, triple-zeta, and quadruple-zeta basis sets for the 4s, 5s, 6s, and 7s elements. , 2009, The journal of physical chemistry. A.

[60]  Johannes Neugebauer,et al.  On the calculation of general response properties in subsystem density functional theory. , 2009, The Journal of chemical physics.

[61]  V. Dolmatov,et al.  Confinement and electron correlation effects in photoionization of atoms in endohedral anions: Ne@Cz−60 , 2009, 0908.2465.

[62]  Kieron Burke,et al.  Partition Density Functional Theory , 2009, 0901.0942.

[63]  Lucas Visscher,et al.  Calculation of local excitations in large systems by embedding wave-function theory in density-functional theory. , 2008, Physical chemistry chemical physics : PCCP.

[64]  H. Quiney,et al.  All-electron four-component Dirac-Kohn-Sham procedure for large molecules and clusters containing heavy elements , 2008 .

[65]  H. Chakraborty,et al.  Photoionization of C60: a model study , 2008 .

[66]  Christoph R. Jacob,et al.  A flexible implementation of frozen‐density embedding for use in multilevel simulations , 2008, J. Comput. Chem..

[67]  Emily A Carter,et al.  Advances in correlated electronic structure methods for solids, surfaces, and nanostructures. , 2008, Annual review of physical chemistry.

[68]  H. Quiney,et al.  Poisson-transformed density fitting in relativistic four-component Dirac-Kohn-Sham theory. , 2008, The Journal of chemical physics.

[69]  M. Kaupp,et al.  A fully relativistic method for calculation of nuclear magnetic shielding tensors with a restricted magnetically balanced basis in the framework of the matrix Dirac-Kohn-Sham equation. , 2008, The Journal of chemical physics.

[70]  Walter Curtis Johnson,et al.  Relativistic Quantum Theory of Atoms and Molecules: Theory and Computation , 2008 .

[71]  H. Takagi,et al.  Superconductivity of doped Ar@C60. , 2007, Chemical communications.

[72]  Kieron Burke,et al.  Partition theory: a very simple illustration. , 2007, The journal of physical chemistry. A.

[73]  Johannes Neugebauer,et al.  Couplings between electronic transitions in a subsystem formulation of time-dependent density functional theory. , 2007, The Journal of chemical physics.

[74]  P. Calaminici,et al.  Density functional theory optimized basis sets for gradient corrected functionals: 3d transition metal systems. , 2007, The Journal of chemical physics.

[75]  D. Suter,et al.  Two-qubit gates between noninteracting qubits in endohedral-fullerene-based quantum computation , 2007 .

[76]  I. P. Grant,et al.  Relativistic quantum theory of atoms and molecules , 2006 .

[77]  Emily A Carter,et al.  Self-consistent embedding theory for locally correlated configuration interaction wave functions in condensed matter. , 2006, The Journal of chemical physics.

[78]  Marcella Iannuzzi,et al.  Density functional embedding for molecular systems , 2006 .

[79]  H. Quiney,et al.  Electron density fitting for the Coulomb problem in relativistic density-functional theory. , 2006, The Journal of chemical physics.

[80]  K. Dyall.  Relativistic Quadruple-Zeta and Revised Triple-Zeta and Double-Zeta Basis Sets for the 4p, 5p, and 6p Elements , 2006 .

[81]  H. Takagi,et al.  Superconductivity of doped Ar@C60 , 2006 .

[82]  Johannes Neugebauer,et al.  An explicit quantum chemical method for modeling large solvation shells applied to aminocoumarin C151. , 2005, The journal of physical chemistry. A.

[83]  J. Tomasi,et al.  Quantum mechanical continuum solvation models. , 2005, Chemical reviews.

[84]  M. Amusia,et al.  Dramatic distortion of the 4d giant resonance by the C60 fullerene shell , 2005 .

[85]  H. Bachau,et al.  Electronic correlation effects in a model of endohedral Mg (Mg@C60) , 2005 .

[86]  K. Dyall.  Relativistic double-zeta, triple-zeta, and quadruple-zeta basis sets for the 5d elements Hf–Hg , 2004 .

[87]  Andreas M Köster,et al.  Calculation of exchange-correlation potentials with auxiliary function densities. , 2004, The Journal of chemical physics.

[88]  H. Quiney,et al.  Relativistic density functional theory using Gaussian basis sets , 2002 .

[89]  Trygve Helgaker,et al.  Four‐component relativistic Kohn–Sham theory , 2002, J. Comput. Chem..

[90]  Kimihiko Hirao,et al.  A new implementation of four-component relativistic density functional method for heavy-atom polyatomic systems , 2001 .

[91]  F. Matthias Bickelhaupt,et al.  Chemistry with ADF , 2001, J. Comput. Chem..

[92]  Guido Rossum,et al.  Python Reference Manual , 2000 .

[93]  F. J. Luque,et al.  Theoretical Methods for the Description of the Solvent Effect in Biomolecular Systems. , 2000, Chemical reviews.

[94]  Eberhard Engel,et al.  Four-component relativistic density functional calculations of heavy diatomic molecules , 2000 .

[95]  A. Baltenkov.  Resonances in photoionization cross sections of inner subshells of atoms inside the fullerene cage , 1999 .

[96]  K. Dyall.  Relativistic and nonrelativistic finite nucleus optimized triple-zeta basis sets for the 4p, 5p and 6p elements , 1998 .

[97]  Emily A. Carter,et al.  Accurate ab initio energetics of extended systems via explicit correlation embedded in a density functional environment , 1998 .

[98]  Lucas Visscher,et al.  Approximate molecular relativistic Dirac-Coulomb calculations using a simple Coulombic correction , 1997 .

[99]  Michael Dolg,et al.  The Beijing four-component density functional program package (BDF) and its application to EuO, EuS, YbO and YbS , 1997 .

[100]  Eric Schwegler,et al.  Fast assembly of the Coulomb matrix: A quantum chemical tree code , 1996 .

[101]  Jan Almlöf,et al.  THE COULOMB OPERATOR IN A GAUSSIAN PRODUCT BASIS , 1995 .

[102]  A. Lembarki,et al.  Obtaining a gradient-corrected kinetic-energy functional from the Perdew-Wang exchange functional. , 1994, Physical review. A, Atomic, molecular, and optical physics.

[103]  Evert Jan Baerends,et al.  Relativistic total energy using regular approximations , 1994 .

[104]  A. Warshel,et al.  Frozen density functional approach for ab initio calculations of solvated molecules , 1993 .

[105]  Cortona,et al.  Direct determination of self-consistent total energies and charge densities of solids: A study of the cohesive properties of the alkali halides. , 1992, Physical review. B, Condensed matter.

[106]  T. Dunning,et al.  Electron affinities of the first‐row atoms revisited. Systematic basis sets and wave functions , 1992 .

[107]  K. Dyall,et al.  Kinetic balance and variational bounds failure in the solution of the Dirac equation in a finite Gaussian basis set , 1990 .

[108]  T. H. Dunning.  Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[109]  A. Becke,et al.  Density-functional exchange-energy approximation with correct asymptotic behavior. , 1988, Physical review. A, General physics.

[110]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[111]  Senatore,et al.  Density dependence of the dielectric constant of rare-gas crystals. , 1986, Physical review. B, Condensed matter.

[112]  S. H. Vosko,et al.  Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis , 1980 .

[113]  M. Levitt,et al.  Theoretical studies of enzymic reactions: dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme. , 1976, Journal of molecular biology.

[114]  J. C. Slater.  A Simplification of the Hartree-Fock Method , 1951 .

[115]  L. H. Thomas.  The calculation of atomic fields , 1927, Mathematical Proceedings of the Cambridge Philosophical Society.

[116]  Loriano Storchi,et al.  BERTHA and PyBERTHA: State of the Art for Full Four-Component Dirac-Kohn-Sham Calculations , 2019, PARCO.

[117]  Lothar Dunsch,et al.  Endohedral fullerenes. , 2013, Chemical reviews.

[118]  Loriano Storchi,et al.  An Efficient Parallel All-Electron Four-Component Dirac-Kohn-Sham Program Using a Distributed Matrix Approach. , 2010, Journal of chemical theory and computation.

[119]  K. Dyall,et al.  Revised relativistic basis sets for the 5d elements Hf–Hg , 2009 .

[120]  T. Wesołowski,et al.  Linearized orbital-free embedding potential in self-consistent calculations† , 2009 .

[121]  P. Calaminici,et al.  Density functional theory optimized basis sets for gradient corrected functionals : 3 d transition metal systems , 2007 .

[122]  Tomasz Adam Wesolowski,et al.  Generalization of the Kohn–Sham equations with constrained electron density formalism and its time‐dependent response theory formulation , 2004 .

[123]  Emily A. Carter,et al.  Periodic density functional embedding theory for complete active space self-consistent field and configuration interaction calculations: Ground and excited states , 2002 .

[124]  I. P. Grant,et al.  Application of relativistic theories and quantum electrodynamics to chemical problems , 2000 .

[125]  L. Dagum,et al.  OpenMP: an industry standard API for shared-memory programming , 1998 .

[126]  Karsten M. Decker,et al.  Programming Environments for Massively Parallel Distributed Systems , 1994, Monte Verità.