Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userATLAS: A Real-Space Finite-Difference Implementation of Orbital-Free Density Functional Theory
351
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Orbital-free density functional theory (OF-DFT) is a promising method for large-scale quantum mechanics simulation as it provides a good balance of accuracy and computational cost. Its applicability to large-scale simulations has been aided by progress in constructing kinetic energy functionals and local pseudopotentials. However, the widespread adoption of OF-DFT requires further improvement in its efficiency and robustly implemented software. Here we develop a real-space finite-difference method for the numerical solution of OF-DFT in periodic systems. Instead of the traditional self-consistent method, a powerful scheme for energy minimization is introduced to solve the Euler--Lagrange equation. Our approach engages both the real-space finite-difference method and a direct energy-minimization scheme for the OF-DFT calculations. The method is coded into the ATLAS software package and benchmarked using periodic systems of solid Mg, Al, and Al$_{3}$Mg. The test results show that our implementation can achieve high accuracy, efficiency, and numerical stability for large-scale simulations.
rate research
Read More
We present a computational scheme for orbital-free density functional theory (OFDFT) that simultaneously provides access to all-electron values and preserves the OFDFT linear scaling as a function of the system size. Using the projector augmented-wave method (PAW) in combination with real-space methods we overcome some obstacles faced by other available implementation schemes. Specifically, the advantages of using the PAW method are two fold. First, PAW reproduces all-electron values offering freedom in adjusting the convergence parameters and the atomic setups allow tuning the numerical accuracy per element. Second, PAW can provide a solution to some of the convergence problems exhibited in other OFDFT implementations based on Kohn-Sham codes. Using PAW and real-space methods, our orbital-free results agree with the reference all-electron values with a mean absolute error of 10~meV and the number of iterations required by the self-consistent cycle is comparable to the KS method. The comparison of all-electron and pseudopotential bulk modulus and lattice constant reveal an enormous difference, demonstrating that in order to assess the performance of OFDFT functionals it is necessary to use implementations that obtain all-electron values. The proposed combination of methods is the most promising route currently available. We finally show that a parametrized kinetic energy functional can give lattice constants and bulk moduli comparable in accuracy to those obtained by the KS PBE method, exemplified with the case of diamond.
We present and motivate an efficient way to include orbital dependent many--body correlations in trial wave function of real--space Quantum Monte Carlo methods for use in electronic structure calculations. We apply our new orbital--dependent backflow wave function to calculate ground state energies of the first row atoms using variational and diffusion Monte Carlo methods. The systematic overall gain of correlation energy with respect to single determinant Jastrow-Slater wave functions is competitive with the best single determinant trial wave functions currently available. The computational cost per Monte Carlo step is comparable to that of simple backflow calculations.
Classical density-functional theory provides an efficient alternative to molecular dynamics simulations for understanding the equilibrium properties of inhomogeneous fluids. However, application of density-functional theory to multi-site molecular fluids has so far been limited by complications due to the implicit molecular geometry constraints on the site densities, whose resolution typically requires expensive Monte Carlo methods. Here, we present a general scheme of circumventing this so-called inversion problem: compressed representations of the orientation density. This approach allows us to combine the superior iterative convergence properties of multipole representations of the fluid configuration with the improved accuracy of site-density functionals. Next, from a computational perspective, we show how to extend the DFT++ algebraic formulation of electronic density-functional theory to the classical fluid case and present a basis-independent discretization of our formulation for molecular classical density-functional theory. Finally, armed with the above general framework, we construct a simplified free-energy functional for water which captures the radial distributions, cavitation energies, and the linear and non-linear dielectric response of liquid water. The resulting approach will enable efficient and reliable first-principles studies of atomic-scale processes in contact with solution or other liquid environments.
The Breit correction, the finite-light-speed correction for the Coulomb interaction of the electron-electron interaction in $ O left( 1/ c^2 right) $, is introduced to density functional theory (DFT) based on the non-relativistic reduction with the local density approximation. Using this newly developed relativistic DFT, it is found that the possible outer-most electron of lawrencium atom is the $ p $ orbital instead of the $ d $ orbital, which is consistent with the previous calculations based on wave-function theory. A possible explanation of the anomalous behavior of its first ionization energy is also given. This DFT scheme provides a practical calculation method for the study of properties of super-heavy elements.
Real-time time-dependent density functional theory (rt-TDDFT) with hybrid exchange-correlation functional has wide-ranging applications in chemistry and material science simulations. However, it can be thousands of times more expensive than a conventional ground state DFT simulation, hence is limited to small systems. In this paper, we accelerate hybrid functional rt-TDDFT calculations using the parallel transport gauge formalism, and the GPU implementation on Summit. Our implementation can efficiently scale to 786 GPUs for a large system with 1536 silicon atoms, and the wall clock time is only 1.5 hours per femtosecond. This unprecedented speed enables the simulation of large systems with more than 1000 atoms using rt-TDDFT and hybrid functional.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
