Local-Scaling Many-Electron Density Functional Theory

Chih-Chuen Lin, Phani Motamarri, Vikram Gavini

Abstract We present a tensor-structured algorithm for efficient large-scale density functional theory (DFT) calculations by constructing a Tucker tensor basis that is adapted to the Kohn–Sham Hamiltonian and localized in real-space. The proposed approach uses an additive separable approximation to the Kohn–Sham Hamiltonian and an L1 localization technique to generate the 1-D localized functions that constitute the Tucker tensor basis. Numerical results show that the resulting Tucker tensor basis exhibits exponential convergence in the ground-state energy with increasing Tucker rank. Further, the proposed tensor-structured algorithm demonstrated sub-quadratic scaling with system-size for both systems with and without a gap, and involving many thousands of atoms. This reduced-order scaling has also resulted in the proposed approach outperforming plane-wave DFT implementation for systems beyond 2000 electrons.

Eugene S. Kryachko, Eduardo V. Ludeña

Swapan K. Ghosh

Keiji Kosaka

The formally exact expression of the exchange-correlation potential in the local-scaling density functional theory (LSDFT) was derived recently by the author, and it is equivalent to that of Holas and March. The basic equation of the LS-DFT was also derived, and it is equivalent to the differential virial theorem (DVT). The LS-DFT has been extended to the local-scaling current density functional theory (LS-CDFT). This has been compared with the formulation of Holas and March. In this note, the LS-CDFT is extended to the timedependent system; i.e., the local-scaling time-dependent current density functional theory (LS-TDCDFT) is formulated. As a straightforward extension of ref. 3, the density and the paramagnetic current density jp are adopted as the basic variables, which is the same procedure used by Colwell et al. It is explicitly shown that the time-dependent exchange-correlation scalar potential vxc and the vector one Axc are same functionals of firstand second-order density matrices as those of the time-independent case, while the total effective potentials are not so. The equation of motion of jp provides the DVT for the time-dependent system in the magnetic field. The time-dependent term of the DVT is equal to that of the nonmagnetic time-dependent system. Throughout the entire formulation, the notation of the tensor product is suppressed. Here, an N-electron system in an external scalar potential vextðr; tÞ, an external vector potential Aextðr; tÞ, and an interaction potential uðr; r0Þ 1⁄4 =jr r0j is considered, where is the coupling constant treated as a variable parameter. In a time-independent theory, the quantity to be minimized is the energy. But in the time-dependent case, the action S plays the role:

Chi Yung Yam, Satoshi Yokojima, GuanHua Chen

A linear-scaling first-principles quantum mechanical method is developed to evaluate the optical responses of large molecular systems. Instead of a many-body wave function, the equation of motion is solved for the reduced single-electron density matrix in the time domain. The locality of the reduced single-electron density matrix is utilized to ensure that computational time scales linearly with system size. The two-electron Coulomb integrals are evaluated with the fast multipole method, and the calculation of exchange-correlation quadratures utilizes the locality of an exchange-correlation functional and the integral prescreening technique. As an illustration, the resulting time-dependent density-functional theory is used to calculate the absorption spectra of polyacetylene oligomers and linear alkanes. The linear-scaling of computational time versus the system size is clearly demonstrated.

Viraht Sahni