Abstract:CONSPECTUS: All-atom methods treat solute and solvent at the same level of electronic structure theory and statistical mechanics. All-atom computation of acidity constants (pKa) and redox potentials is still a challenge. In this Account, we review such a method combining density functional theory based molecular dynamics (DFTMD) and free energy perturbation (FEP) methods. The key computational tool is a FEP based method for reversible insertion of a proton or electron in a periodic DFTMD model system. The free energy of insertion (work function) is computed by thermodynamic integration of vertical energy gaps obtained from total energy differences. The problem of the loss of a physical reference for ionization energies under periodic boundary conditions is solved by comparing with the proton work function computed for the same supercell. The scheme acts as a computational hydrogen electrode, and the DFTMD redox energies can be directly compared with experimental redox potentials. Consistent with the closed shell nature of acid dissociation, pKa estimates computed using the proton insertion/removal scheme are found to be significantly more accurate than the redox potential calculations. This enables us to separate the DFT error from other sources of uncertainty such as finite system size and sampling errors. Drawing an analogy with charged defects in solids, we trace the error in redox potentials back to underestimation of the energy gap of the extended states of the solvent. Accordingly the improvement in the redox potential as calculated by hybrid functionals is explained as a consequence of the opening up of the bandgap by the Hartree-Fock exchange component in hybrids. Test calculations for a number of small inorganic and organic molecules show that the hybrid functional implementation of our method can reproduce acidity constants with an uncertainty of 1-2 pKa units (0.1 eV). The error for redox potentials is in the order of 0.2 V.

Energy Gap from Step Structure of the Analytically Inverted Non-Additive Kinetic Potential

Ionization Potentials and Fundamental Gaps in Atomic Systems from the Ensemble-DFT Approach

Developing Orbital-Dependent Corrections for the Non-Additive Kinetic Energy in Subsystem Density Functional Theory

Ultranonlocality and accurate band gaps from a meta-generalized gradient approximation

Origin of static and dynamic steps in exact Kohn-Sham potentials

Importance of the Kinetic Energy Density for Band Gap Calculations in Solids with Density Functional Theory

Understanding Band Gaps of Solids in Generalized Kohn-Sham Theory.

A new Green's function formalism for kinetic energy density functional for atomic and molecular system: Emergence of $N-$dependence using model potentials

Exchange energies with forces in density-functional theory

Alternative Approach to the Optimized Effective Potential Method

Spin gaps and spin-flip energies in density-functional theory

Generalized energy gap law: An open system dynamics approach to non-adiabatic phenomena in molecules

Computation of the Kohn-Sham Orbital Kinetic Energy Density in the Full-Potential Linearized Augmented Plane-Wave Method

Redox potentials and acidity constants from density functional theory based molecular dynamics.

Upper Bound to the Gradient-Based Kinetic Energy Density of Noninteracting Electrons in an External Potential

Fast and stable tight-binding framework for nonlocal kinetic energy density functional reconstruction in orbital-free density functional calculations

Molecular Binding Energies from Partition Density Functional Theory

Insights into Direct Methods for Predictions of Ionization Potential and Electron Affinity in Density Functional Theory.

Training Machine-Learned Density Functionals on Band Gaps

Exchange-correlation potential built on the derivative discontinuity of electron density

Implementation and benchmark of a long-range corrected functional in the density functional based tight-binding method