Optimization of constrained density functional theory

Constrained density functional theory (cDFT) is a versatile electronic structure method that enables ground-state calculations to be performed subject to physical constraints. It thereby broadens their applicability and utility. Automated Lagrange multiplier optimisation is necessary for multiple constraints to be applied efficiently in cDFT, for it to be used in tandem with geometry optimization, or with molecular dynamics. In order to facilitate this, we comprehensively develop the connection between cDFT energy derivatives and response functions, providing a rigorous assessment of the uniqueness and character of cDFT stationary points while accounting for electronic interactions and screening. In particular, we provide a new, non-perturbative proof that stable stationary points of linear density constraints occur only at energy maxima with respect to their Lagrange multipliers. We show that multiple solutions, hysteresis, and energy discontinuities may occur in cDFT. Expressions are derived, in terms of convenient by-products of cDFT optimization, for quantities such as the dielectric function and a condition number quantifying ill-definition in multi-constraint cDFT.Comment: 15 pages, 6 figure

Generalized Wannier functions: a comparison of molecular electric dipole polarizabilities

Localized Wannier functions provide an efficient and intuitive means by which to compute dielectric properties from first principles. They are most commonly constructed in a post-processing step, following total-energy minimization. Nonorthogonal generalized Wannier functions (NGWFs) [Skylaris et al., Phys. Rev. B 66, 035119 11 (2002); Skylaris et al., J. Chem. Phys. 122, 084119 (2005)] may also be optimized in situ, in the process of solving for the ground-state density. We explore the relationship between NGWFs and orthonormal, maximally localized Wannier functions (MLWFs) [Marzari and Vanderbilt, Phys. Rev. B 56, 12847 (1997); Souza, Marzari, and Vanderbilt, ibid. 65, 035109 (2001)], demonstrating that NGWFs may be used to compute electric dipole polarizabilities efficiently, with no necessity for post-processing optimization, and with an accuracy comparable to MLWFs.Comment: 5 pages, 1 figure. This version matches that accepted for Physical Review B on 4th May 201

Subspace representations in ab initio methods for strongly correlated systems

We present a generalized definition of subspace occupancy matrices in ab initio methods for strongly correlated materials, such as DFT+U and DFT+DMFT, which is appropriate to the case of nonorthogonal projector functions. By enforcing the tensorial consistency of all matrix operations, we are led to a subspace projection operator for which the occupancy matrix is tensorial and accumulates only contributions which are local to the correlated subspace at hand. For DFT+U in particular, the resulting contributions to the potential and ionic forces are automatically Hermitian, without resort to symmetrization, and localized to their corresponding correlated subspace. The tensorial invariance of the occupancies, energies and ionic forces is preserved. We illustrate the effect of this formalism in a DFT+U study using self-consistently determined projectors.Comment: 15 pages, 8 figures. This version (v2) matches that accepted for Physical Review B on 15th April 201

Ligand Discrimination in Myoglobin from Linear-Scaling DFT+U

Myoglobin modulates the binding of diatomic molecules to its heme group via hydrogen-bonding and steric interactions with neighboring residues, and is an important benchmark for computational studies of biomolecules. We have performed calculations on the heme binding site and a significant proportion of the protein environment (more than 1000 atoms) using linear-scaling density functional theory and the DFT+U method to correct for self-interaction errors associated with localized 3d states. We confirm both the hydrogen-bonding nature of the discrimination effect (3.6 kcal/mol) and assumptions that the relative strain energy stored in the protein is low (less than 1 kcal/mol). Our calculations significantly widen the scope for tackling problems in drug design and enzymology, especially in cases where electron localization, allostery or long-ranged polarization influence ligand binding and reaction.Comment: 15 pages, 3 figures. Supplementary material 8 pages, 3 figures. This version matches that accepted for J. Phys. Chem. Lett. on 10th May 201

Optimization strategies developed on NiO for Heisenberg exchange coupling calculations using projector augmented wave based first-principles DFT+U+J

High-performance batteries, heterogeneous catalysts and next-generation photovoltaics often centrally involve transition metal oxides (TMOs) that undergo charge or spin-state changes. Demand for accurate DFT modeling of TMOs has increased in recent years, driving improved quantification and correction schemes for approximate DFT's characteristic errors, notably those pertaining to self-interaction and static correlation. Of considerable interest, meanwhile, is the use of DFT-accessible quantities to compute parameters of coarse-grained models such as for magnetism. To understand the interference of error corrections and model mappings, we probe the prototypical Mott-Hubbard insulator NiO, calculating its electronic structure in its antiferromagnetic I/II and ferromagnetic states. We examine the pronounced sensitivity of the first principles calculated Hubbard U and Hund's J parameters to choices concerning Projector Augmented Wave (PAW) based population analysis, we reevaluate spin quantification conventions for the Heisenberg model, and we seek to develop best practices for calculating Hubbard parameters specific to energetically meta-stable magnetic orderings of TMOs. Within this framework, we assess several corrective functionals using in situ calculated U and J parameters, e.g., DFT+U and DFT+U+J. We find that while using a straightforward workflow with minimal empiricism, the NiO Heisenberg parameter RMS error with respect to experiment was reduced to 13%, an advance upon the state-of-the-art. Methodologically, we used a linear-response implementation for calculating the Hubbard U available in the open-source plane-wave DFT code Abinit. We have extended its utility to calculate the Hund's exchange coupling J, however our findings are anticipated to be applicable to any DFT+U implementation.Comment: 19 pages, 6 figures (+1 in SI), 7 tables (+1 in SI