Wannier-function-based constrained DFT with nonorthogonality-correcting Pulay forces in application to the reorganization effects in graphene-adsorbed pentacene
arXiv:1802.01669 · doi:10.1103/PhysRevB.97.205120
Abstract
Pulay terms arise in the Hellman-Feynman forces in electronic structure calculations when one employs a basis set made of localized orbitals that move with their host atoms. If the total energy of the system depends on a subspace population defined in terms of the localized orbitals across multiple atoms, then unconventional Pulay terms will emerge due to the variation of the orbital nonorthogonality with ionic translation. Here, we derive the required exact expressions for such terms, which cannot be eliminated by orbital orthonormalization. We have implemented these corrected ionic forces within the linear-scaling density functional theory (DFT) package ONETEP, and have used constrained DFT to calculate the reorganization energy of a pentacene molecule adsorbed on a graphene flake. The calculations are performed by including ensemble DFT, corrections for periodic boundary conditions, and empirical Van der Waals interactions. For this system we find that tensorially invariant population analysis yields an adsorbate subspace population that is very close to integer-valued when based upon nonorthogonal Wannier functions, and also but less precisely when using pseudoatomic functions. Thus, orbitals can provide a very effective population analysis for constrained DFT. Our calculations show that the reorganization energy of the adsorbed pentacene is typically lower than that of pentacene in the gas phase. We attribute this effect to steric hindrance.
arXiv comment12 pages, 6 figures, and 2 tables. As accepted for publication in Physical Review B on 25th April 2018
References in corpus (23)
- Maximally-localized generalized Wannier functions for composite energy bands
- Global Optimization by Basin-Hopping and the Lowest Energy Structures of Lennard-Jones Clusters Containing up to 110 Atoms
- Electronic Structure Calculations with Dynamical Mean-Field Theory: A Spectral Density Functional Approach
- Hubbard-corrected DFT energy functionals: the LDA+U description of correlated systems
- Dynamical Mean-Field Theory within the Full-Potential Methods: Electronic structure of Ce-115 materials
- First-principles Calculations of the Electronic Structure and Spectra of Strongly Correlated Systems: Dynamical Mean-field Theory
- Ensemble density-functional theory for ab-initio molecular dynamics of metals and finite-temperature insulators
- Dynamical mean-field theory using Wannier functions: a flexible route to electronic structure calculations of strongly correlated materials
- Accurate ionic forces and geometry optimisation in linear scaling density-functional theory with local orbitals
- Electrostatic Interactions in Finite Systems treated with Periodic Boundary Conditions: Application to Linear-Scaling Density Functional Theory
- Subspace representations in ab initio methods for strongly correlated systems
- Constrained-DFT method for accurate energy level alignment of metal/molecule interfaces
- Projector self-consistent DFT+U using non-orthogonal generalized Wannier functions
- Towards a full ab initio theory of strong electronic correlations in nanoscale devices
- Electronic correlations and crystal structure distortions in BaBiO3
- Theory of projections with non-orthogonal basis sets: Partitioning techniques and effective Hamiltonians
- Quantum mechanics in an evolving Hilbert space
- Linear-scaling DFT+U with full local orbital optimization
- Optimization of constrained density functional theory
- Supercell convergence of charge-transfer energies in pentacene molecular crystals from constrained DFT
- Hellmann-Feynman Forces within the DFT+U in Wannier functions basis
- Charge transfer energies of benzene physisorbed on a graphene sheet from constrained density functional theory
- Generalized Wannier functions: a comparison of molecular electric dipole polarizabilities