Exact relations between the electron density and external potential for systems of interacting and noninteracting electrons

The differential virial theorem (DVT) is an explicit relation between the electron density ρ( r ), the external potential, kinetic energy density tensor, and (for interacting electrons) the pair function. The time‐dependent generalization of this relation also involves the paramagnetic current density. We present a detailed unified derivation of all known variants of the DVT starting from a modified equation of motion for the current density. To emphasize the practical significance of the theorem for noninteracting electrons, we cast it in a form best suited for recovering the Kohn–Sham effective potential v_s( r ) from a given electron density. The resulting expression contains only ρ( r ), v_s( r ), kinetic energy density, and a new orbital‐dependent ingredient containing only occupied Kohn–Sham orbitals. Other possible applications of the theorem are also briefly discussed. © 2012 Wiley Periodicals, Inc.     

Investigation of dephasing in an open quantum system under chaotic influence via a fractional Kohn–Sham scheme

In this work, the dynamics of dephasing (without relaxation) in the presence of a chaotic oscillator is theoretically investigated. The time‐dependent density functional theory framework was used in tandem with the Lindblad master equation approach for modeling the dissipative dynamics. Using the Kohn–Sham (K–S) scheme under certain approximations, the exact model for the potentials was acquired. In addition, a space‐fractional K–S scheme was developed (using the modified Riemann–Liouville operator) for modeling the dephasing phenomenon. Extensive analyses and comparative studies were then done on the results obtained using the space‐fractional K–S system and the conventional K–S system. © 2014 Wiley Periodicals, Inc.     

Cross approximation in tensor electron density computations

We propose new tensor approximation algorithms for certain discrete functions related with Hartree–Fock/Kohn–Sham equations. Given a canonical tensor representation for the electron density function (for example, produced by quantum chemistry packages such as MOLPRO), we obtain its Tucker approximation with much fewer parameters than the input data and the Tucker approximation for the cubic root of this function, which is part of the Kohn–Sham exchange operator. The key idea is in the fast and accurate prefiltering of possibly large‐scale factors of the canonical tensor input. The new algorithms are based on the incomplete cross approximation method applied to matrices and tensors of order 3 and outperform other tools for the same purpose. First, we show that the cross approximation method is robust and much faster than the singular value decomposition‐based approach. As a consequence, it becomes possible to increase the resolution of grid and the complexity of molecules that can be handled by the Hartree–Fock chemical models. Second, we propose a new fast approximation method for f^1/3(x, y, z), based on the factor prefiltering method for f(x, y, z) and certain mimic approximation hypothesis. Third, we conclude that the Tucker format has advantages in the storage and computation time compared with the ubiquitous canonical format. Copyright © 2009 John Wiley & Sons, Ltd.     

Alternative Representations of the Correlation Energy in Density‐Functional Theory: A Kinetic‐Energy Based Adiabatic Connection

The adiabatic‐connection framework has been widely used to explore the properties of the correlation energy in density‐functional theory. The integrand in this formula may be expressed in terms of the electron–electron interactions directly, involving intrinsically two‐particle expectation values. Alternatively, it may be expressed in terms of the kinetic energy, involving only one‐particle quantities. In this work, we explore this alternative representation for the correlation energy and highlight some of its potential for the construction of new density functional approximations. The kinetic‐energy based integrand is effective in concentrating static correlation effects to the low interaction strength regime and approaches zero asymptotically, offering interesting new possibilities for modeling the correlation energy in density‐functional theory     

Neural network and ReaxFF comparison for Au properties

We have studied how ReaxFF and Behler–Parrinello neural network (BPNN) atomistic potentials should be trained to be accurate and tractable across multiple structural regimes of Au as a representative example of a single‐component material. We trained these potentials using subsets of 9,972 Kohn‐Sham density functional theory calculations and then validated their predictions against the untrained data. Our best ReaxFF potential was trained from 848 data points and could reliably predict surface and bulk data; however, it was substantially less accurate for molecular clusters of 126 atoms or fewer. Training the ReaxFF potential to more data also resulted in overfitting and lower accuracy. In contrast, BPNN could be fit to 9,734 calculations, and this potential performed comparably or better than ReaxFF across all regimes. However, the BPNN potential in this implementation brings significantly higher computational cost. © 2016 Wiley Periodicals, Inc.     

Numerical integration of the Kohn–Sham equations: Integral method for bound states

A numerical method is presented that solves the multicenter Kohn–Sham equations. The method couples the resolution of the integral form of the equation at a given energy with an iterative search for the eigenvalues. The validity of the method is checked by comparing some test calculations for diatomics with results in the literature from other numerical methods. For these calculations the wave functions are expanded in partial waves either on one center or on two centers with the help of the partitioning of space in fuzzy cells. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 74: 49–54, 1999     

Kohn-Sham Density Matrix and the Kernel Energy MethodSCI北大核心CSCD

The kernel energy method(KEM) has been shown to provide fast and accurate molecular energy calculations for molecules at their equilibrium geometries.KEM breaks a molecule into smaller subsets,called kernels,for the purposes of calculation.The results from the kernels are summed according to an expression characteristic of KEM to obtain the full molecule energy.A generalization of the kernel expansion to density matrices provides the full molecule density matrix and orbitals.In this study,the kernel expansion for the density matrix is examined in the context of density functional theory(DFT) Kohn-Sham(KS) calculations.A kernel expansion for the one-body density matrix analogous to the kernel expansion for energy is defined,and is then converted into a normalizedprojector by using the Clinton algorithm.Such normalized projectors are factorizable into linear combination of atomic orbitals(LCAO) matrices that deliver full-molecule Kohn-Sham molecular orbitals in the atomic orbital basis.Both straightforward KEM energies and energies from a normalized,idempotent density matrix obtained from a density matrix kernel expansion to which the Clinton algorithm has been applied are compared to reference energies obtained from calculations on the full system without any kernel expansion.Calculations were performed both for a simple proof-of-concept system consisting of three atoms in a linear configuration and for a water cluster consisting of twelve water molecules.In the case of the proof-of-concept system,calculations were performed using the STO-3 G and6-31 G(d,p) bases over a range of atomic separations,some very far from equilibrium.The water cluster was calculated in the 6-31 G(d,p) basis at an equilibrium geometry.The normalized projector density energies are more accurate than the straightforward KEM energy results in nearly all cases.In the case of the water cluster,the energy of the normalized projector is approximately four times more accurate than the straightforward KEM energy result.The KS density matrices of this study are applicable to quantum crystallography.     

The Liu‐Parr power series expansion of the Pauli kinetic energy functional with the incorporation of shell‐inducing traits: Atoms

An approximate expression for the Pauli kinetic energy functional T_p is advanced in terms of the Liu‐Parr expansion [S. Liu, R.G. Parr, Phys. Rev. A 1997 , 55, 1792] which involves a power series of the one‐electron density. We use this explicit functional for T_p to compute the value of the noninteracting kinetic energy functional T_s of 34 atoms, from Li to Kr (and their positive and negative monoions). In particular, we examine the effect that a shell‐by‐shell mean‐square optimization of the expansion coefficients has on the kinetic energy values and explore the effect that the size of the expansion, given by the parameter n, has on the accuracy of the approximation. The results yield a mean absolute percent error for 34 neutral atoms of 0.15, 0.08, 0.04, 0.03, and 0.01 for expansions with n = 3, 4, 5, 6, and 7, respectively (where ). We show that these results, which are the most accurate ones obtained to date for the representation of the noninteracting kinetic energy functional, stem from the imposition of shell‐inducing traits. We also compare these Liu‐Parr functionals with the exact but nonexplicit functional generated in the local‐scaling transformation version of DFT.     

Density-based Globally Convergent Trust-region Methods for Self-consistent Field Electronic Structure Calculations

A theory of globally convergent trust-region methods for self-consistent field electronic structure calculations that use the density matrices as variables is developed. The optimization is performed by means of sequential global minimizations of a quadratic model of the true energy. The global minimization of this quadratic model, subject to the idempotency of the density matrix and the rank constraint, coincides with the fixed-point iteration. We prove that the global minimization of this quadratic model subject to the restrictions and smaller trust regions corresponds to the solution of level-shifted equations. The precise implementation of algorithms leading to global convergence is stated and a proof of global convergence is provided. Numerical experiments confirm theoretical predictions and practical convergence is obtained for difficult cases, even if their geometries are highly distorted. The reduction of the trust region is performed by a strategy that uses the structure of the energy function providing the algorithm with a nice practical behavior. This framework may be applied to any problem with idempotency constraints and for which the derivative of the objective function is a symmetric matrix. Therefore, application to calculations based both on Hartree–Fock or Kohn–Sham density functional theory are straightforward.