Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userHighly accurate local basis sets for large-scale DFT calculations in CONQUEST
222
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Given the widespread use of density functional theory (DFT), there is an increasing need for the ability to model large systems (beyond 1,000 atoms). We present a brief overview of the large-scale DFT code Conquest, which is capable of modelling such large systems, and discuss approaches to the generation of consistent, well-converged pseudo-atomic basis sets which will allow such large scale calculations. We present tests of these basis sets for a variety of materials, comparing to fully converged plane wave results using the same pseudopotentials and grids.
rate research
Read More
We introduce numerical optimization of multi-site support functions in the linear-scaling DFT code CONQUEST. Multi-site support functions, which are linear combinations of pseudo-atomic orbitals on a target atom and those neighbours within a cutoff, have been recently proposed to reduce the number of support functions to the minimal basis while keeping the accuracy of a large basis [J. Chem. Theory Comput., 2014, 10, 4813]. The coefficients were determined by using the local filter diagonalization (LFD) method [Phys. Rev. B, 2009, 80, 205104]. We analyse the effect of numerical optimization of the coefficients produced by the LFD method. Tests on crystalline silicon, a benzene molecule and hydrated DNA systems show that the optimization improves the accuracy of the multi-site support functions with small cutoffs. It is also confirmed that the optimization guarantees the variational energy minimizations with multi-site support functions.
Density Functional Theory (DFT) has become the quasi-standard for ab-initio simulations for a wide range of applications. While the intrinsic cubic scaling of DFT was for a long time limiting the accessible system size to some hundred atoms, the recent progress with respect to linear scaling DFT methods has allowed to tackle problems that are larger by many orders of magnitudes. However, as these linear scaling methods were developed for insulators, they cannot, in general, be straightforwardly applied to metals, as a finite temperature is needed to ensure locality of the density matrix. In this paper we show that, once finite electronic temperature is employed, the linear scaling version of the BigDFT code is able to exploit this locality to provide a computational treatment that scales linearly with respect to the number of atoms of a metallic system. We provide prototype examples based on bulk Tungsten, which plays a key role in finding safe and long-lasting materials for Fusion Reactors. We believe that such an approach might help in opening the path towards novel approaches for investigating the electronic structure of such materials, in particular when large supercells are required.
We survey the underlying theory behind the large-scale and linear scaling DFT code, Conquest, which shows excellent parallel scaling and can be applied to thousands of atoms with exact solutions, and millions of atoms with linear scaling. We give details of the representation of the density matrix and the approach to finding the electronic ground state, and discuss the implementation of molecular dynamics with linear scaling. We give an overview of the performance of the code, focussing in particular on the parallel scaling, and provide examples of recent developments and applications.
Density Functional Theory calculations traditionally suffer from an inherent cubic scaling with respect to the size of the system, making big calculations extremely expensive. This cubic scaling can be avoided by the use of so-called linear scaling algorithms, which have been developed during the last few decades. In this way it becomes possible to perform ab-initio calculations for several tens of thousands of atoms or even more within a reasonable time frame. However, even though the use of linear scaling algorithms is physically well justified, their implementation often introduces some small errors. Consequently most implementations offering such a linear complexity either yield only a limited accuracy or, if one wants to go beyond this restriction, require a tedious fine tuning of many parameters. In our linear scaling approach within the BigDFT package, we were able to overcome this restriction. Using an ansatz based on localized support functions expressed in an underlying Daubechies wavelet basis -- which offers ideal properties for accurate linear scaling calculations -- we obtain an amazingly high accuracy and a universal applicability while still keeping the possibility of simulating large systems with only a moderate demand of computing resources.
One of the goals in the development of large scale electronic structure methods is to perform calculations explicitly for a localised region of a system, while still taking into account the rest of the system outside of this region. An example of this in surface physics would be to embed an adsorbate and a few surface atoms into an extended substrate, hence considerably reducing computational costs. Here we apply the constrained electron density method of embedding a Kohn-Sham system in a substrate system (first described by P. Cortonacite{1} and T.A. Wesolowskicite{2}), within a plane-wave basis and pseudopotential framework. This approach divides the charge density of the system into substrate and embedded charge densities, the sum of which is the charge density of the actual system of interest. Two test cases are considered. First we construct fcc bulk aluminium by embedding one cubic lattice of atoms within another. Second, we examine a model surface/adsorbate system of aluminium on aluminium and compare with full Kohn-Sham results.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
