Towards an exact orbital-free single-particle kinetic energy density for the inhomogeneous electron liquid in the Be atom

Holas and March (Phys. Rev. A51, 2040 (1995)) wrote the gradient of the one-body potential V(r) in terms of low-order derivatives of the idempotent Dirac density matrix built from a single Slater determinant of Kohn-Sham orbitals. Here, this is first combined with the study of Dawson and March (J. Chem. Phys. 81, 5850 (1984)) to express the single-particle kinetic energy density of the Be atom ground-state in terms of both the electron density n(r) and potential V(r). While this is the more compact formulation, we then, by removing V(r), demonstrate that the ratio t(r)/n(r) depends, though non-locally, only on the single variable n'(r)/n(r), no high-order gradients entering for the spherical Be atom.Comment: Submitted to Journal of Mathematical Chemistr

Accurate interaction energies at DFT level by means of an efficient dispersion correction

This paper presents an approach for obtaining accurate interaction energies at the DFT level for systems where dispersion interactions are important. This approach combines Becke and Johnson's [J. Chem. Phys. 127, 154108 (2007)] method for the evaluation of dispersion energy corrections and a Hirshfeld method for partitioning of molecular polarizability tensors into atomic contributions. Due to the availability of atomic polarizability tensors, the method is extended to incorporate anisotropic contributions, which prove to be important for complexes of lower symmetry. The method is validated for a set of eighteen complexes, for which interaction energies were obtained with the B3LYP, PBE and TPSS functionals combined with the aug-cc-pVTZ basis set and compared with the values obtained at CCSD(T) level extrapolated to a complete basis set limit. It is shown that very good quality interaction energies can be obtained by the proposed method for each of the examined functionals, the overall performance of the TPSS functional being the best, which with a slope of 1.00 in the linear regression equation and a constant term of only 0.1 kcal/mol allows to obtain accurate interaction energies without any need of a damping function for complexes close to their exact equilibrium geometry

Local softness, softness dipole and polarizabilities of functional groups: application to the side chains of the twenty amino acids

The values of molecular polarizabilities and softnesses of the twenty amino acids were computed ab initio (MP2). By using the iterative Hirshfeld scheme to partition the molecular electronic properties, we demonstrate that the values of the softness of the side chain of the twenty amino acid are clustered in groups reflecting their biochemical classification, namely: aliphatic, basic, acidic, sulfur containing, and aromatic amino acids . The present findings are in agreement with previous results using different approximations and partitioning schemes [P. Senet and F. Aparicio, J. Chem. Phys. 126,145105 (2007)]. In addition, we show that the polarizability of the side chain of an amino acid depends mainly on its number of electrons (reflecting its size) and consequently cannot be used to cluster the amino acids in different biochemical groups, in contrast to the local softness. Our results also demonstrate that the global softness is not simply proportional to the global polarizability in disagreement with the intuition that "a softer moiety is also more polarizable". Amino acids with the same softness may have a polarizability differing by a factor as large as 1.7. This discrepancy can be understood from first principles as we show that the molecular polarizability depends on a "softness dipole vector" and not simply on the global softness

eQE: An open‐source density functional embedding theory code for the condensed phase

AbstractIn this work, we present the main features and algorithmic details of a novel implementation of the frozen density embedding (FDE) formulation of subsystem density functional theory (DFT) that is specifically designed to enable ab initio molecular dynamics (AIMD) simulations of large‐scale condensed‐phase systems containing 1000s of atoms. This code (available at http://eqe.rutgers.edu) has been given the moniker of embedded Quantum ESPRESSO (eQE) as it is a generalization of the open‐source Quantum ESPRESSO (QE) suite of programs. The strengths of eQE reside in a hierarchical parallelization scheme that allows for an efficient and fully self‐consistent treatment of the electronic structure (via the addition of an additional DIIS extrapolation layer) while simultaneously exploiting the inherent symmetries and periodicities in the system (via sampling of subsystem‐specific first Brillouin zones and utilization of subsystem‐specific basis sets). While bulk liquids and molecular crystals are two classes of systems that exemplify the utility of the FDE approach (as these systems can be partitioned into weakly interacting subunits), we show that eQE has significantly extended this regime of applicability by outperforming standard semilocal Kohn–Sham DFT (KS‐DFT) for large‐scale heterogeneous catalysts with quite different layer‐specific electronic structure and intrinsic periodicities. eQE features very favorable strong parallel scaling for a model system of bulk liquid water composed of 256 water molecules, which allows for a significant decrease in the overall time to solution when compared to KS‐DFT. We show that eQE achieves speedups greater than one order of magnitude ( ) when performing AIMD simulations of such large‐scale condensed‐phase systems as: (1) molecular liquids via bulk liquid water represented by 1024 independent water molecules (3072 atoms with a 25.3× speedup over KS‐DFT), (2) polypeptide/biomolecule solvation via (gly)6 solvated in (H2O)395 (1230 atoms with a 38.6× speedup over KS‐DFT), and (3) molecular crystals via a 3 × 3 × 3 periodic supercell of pentacene (1940 atoms with a 12.0× speedup over KS‐DFT). These results represent a significant improvement over the current state‐of‐the‐art and now enable subsystem DFT‐based AIMD simulations of realistically sized condensed‐phase systems of interest throughout chemistry, physics, and materials science