skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Extending density functional embedding theory for covalently bonded systems

; ORCiD logo
Publication Date:
Sponsoring Org.:
USDOE Office of Science (SC), Basic Energy Sciences (BES) (SC-22)
OSTI Identifier:
Grant/Contract Number:
Resource Type:
Journal Article: Published Article
Journal Name:
Proceedings of the National Academy of Sciences of the United States of America
Additional Journal Information:
Journal Volume: 114; Journal Issue: 51; Related Information: CHORUS Timestamp: 2017-12-19 14:46:01; Journal ID: ISSN 0027-8424
Proceedings of the National Academy of Sciences
Country of Publication:
United States

Citation Formats

Yu, Kuang, and Carter, Emily A. Extending density functional embedding theory for covalently bonded systems. United States: N. p., 2017. Web. doi:10.1073/pnas.1712611114.
Yu, Kuang, & Carter, Emily A. Extending density functional embedding theory for covalently bonded systems. United States. doi:10.1073/pnas.1712611114.
Yu, Kuang, and Carter, Emily A. 2017. "Extending density functional embedding theory for covalently bonded systems". United States. doi:10.1073/pnas.1712611114.
title = {Extending density functional embedding theory for covalently bonded systems},
author = {Yu, Kuang and Carter, Emily A.},
abstractNote = {},
doi = {10.1073/pnas.1712611114},
journal = {Proceedings of the National Academy of Sciences of the United States of America},
number = 51,
volume = 114,
place = {United States},
year = 2017,
month =

Journal Article:
Free Publicly Available Full Text
This content will become publicly available on June 19, 2018
Publisher's Version of Record

Save / Share:
  • A key element in the density functional embedding theory (DFET) is the embedding potential. We discuss two major issues related to the embedding potential: (1) its non-uniqueness and (2) the numerical difficulty for solving for it, especially for the spin-polarized systems. To resolve the first issue, we extend DFET to finite temperature: all quantities, such as the subsystem densities and the total system’s density, are calculated at a finite temperature. This is a physical extension since materials work at finite temperatures. We show that the embedding potential is strictly unique at T > 0. To resolve the second issue, wemore » introduce an efficient iterative embedding potential solver. We discuss how to relax the magnetic moments in subsystems and how to equilibrate the chemical potentials across subsystems. The solver is robust and efficient for several non-trivial examples, in all of which good quality spin-polarized embedding potentials were obtained. We also demonstrate the solver on an extended periodic system: iron body-centered cubic (110) surface, which is related to the modeling of the heterogeneous catalysis involving iron, such as the Fischer-Tropsch and the Haber processes. This work would make it efficient and accurate to perform embedding simulations of some challenging material problems, such as the heterogeneous catalysis and the defects of complicated spin configurations in electronic materials.« less
  • We evaluate the accuracy of density functional theory quantum calculations of biomolecular subsystems using a simple electrostatic embedding scheme. Our scheme is based on dividing the system of interest into a primary and secondary subsystem. A finite difference discretization of the Kohn-Sham equations is used for the primary subsystem, while its electrostatic environment is modeled with a simple one-electron potential. Force-field atomic partial charges are used to generate smeared Gaussian charge densities and to model the secondary subsystem. We illustrate the utility of this approach with calculations of truncated dipeptide chains. We analyze quantitatively the accuracy of this approach bymore » calculating atomic forces and comparing results with fullQMcalculations. The impact of the choice made in terminating dangling bonds at the frontier of the QM region is also investigated.« less
  • We describe a novel approach for the calculation of local electric dipole moments for periodic systems. Since the position operator is ill-defined in periodic systems, maximally localized Wannier functions based on the Berry-phase approach are usually employed for the evaluation of local contributions to the total electric dipole moment of the system. We propose an alternative approach: within a subsystem-density functional theory based embedding scheme, subset electric dipole moments are derived without any additional localization procedure, both for hybrid and non-hybrid exchange–correlation functionals. This opens the way to a computationally efficient evaluation of local electric dipole moments in (molecular) periodicmore » systems as well as their rigorous splitting into atomic electric dipole moments. As examples, Infrared spectra of liquid ethylene carbonate and dimethyl carbonate are presented, which are commonly employed as solvents in Lithium ion batteries.« less
  • We analyze the sources of error in quantum embedding calculations in which an active subsystem is treated using wavefunction methods, and the remainder using density functional theory. We show that the embedding potential felt by the electrons in the active subsystem makes only a small contribution to the error of the method, whereas the error in the nonadditive exchange-correlation energy dominates. We test an MP2 correction for this term and demonstrate that the corrected embedding scheme accurately reproduces wavefunction calculations for a series of chemical reactions. Our projector-based embedding method uses localized occupied orbitals to partition the system; as withmore » other local correlation methods, abrupt changes in the character of the localized orbitals along a reaction coordinate can lead to discontinuities in the embedded energy, but we show that these discontinuities are small and can be systematically reduced by increasing the size of the active region. Convergence of reaction energies with respect to the size of the active subsystem is shown to be rapid for all cases where the density functional treatment is able to capture the polarization of the environment, even in conjugated systems, and even when the partition cuts across a double bond.« less
  • We present here the coupling of a polarizable embedding (PE) model to the recently developed multiconfiguration short-range density functional theory method (MC-srDFT), which can treat multiconfigurational systems with a simultaneous account for dynamical and static correlation effects. PE-MC-srDFT is designed to combine efficient treatment of complicated electronic structures with inclusion of effects from the surrounding environment. The environmental effects encompass classical electrostatic interactions as well as polarization of both the quantum region and the environment. Using response theory, molecular properties such as excitation energies and oscillator strengths can be obtained. The PE-MC-srDFT method and the additional terms required for linearmore » response have been implemented in a development version of DALTON. To benchmark the PE-MC-srDFT approach against the literature data, we have investigated the low-lying electronic excitations of acetone and uracil, both immersed in water solution. The PE-MC-srDFT results are consistent and accurate, both in terms of the calculated solvent shift and, unlike regular PE-MCSCF, also with respect to the individual absolute excitation energies. To demonstrate the capabilities of PE-MC-srDFT, we also investigated the retinylidene Schiff base chromophore embedded in the channelrhodopsin protein. While using a much more compact reference wave function in terms of active space, our PE-MC-srDFT approach yields excitation energies comparable in quality to CASSCF/CASPT2 benchmarks.« less