Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new usereQE 2.0: Subsystem DFT Beyond GGA Functionals
425
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
By adopting a divide-and-conquer strategy, subsystem-DFT (sDFT) can dramatically reduce the computational cost of large-scale electronic structure calculations. The key ingredients of sDFT are the nonadditive kinetic energy and exchange-correlation functionals which dominate its accuracy. Even though, semilocal nonadditive functionals find a broad range of applications, their accuracy is somewhat limited especially for those systems where achieving balance between exchange-correlation interactions on one side and nonadditive kinetic energy on the other is crucial. In eQE 2.0, we improve dramatically the accuracy of sDFT simulations by (1) implementing nonlocal nonadditive kinetic energy functionals based on the LMGP family of functionals; (2) adapting Quantum ESPRESSOs implementation of rVV10 and vdW-DF nonlocal exchange-correlation functionals to be employed in sDFT simulations; (3) implementing deorbitalized meta GGA functionals (e.g., SCAN-L). We carefully assess the performance of the newly implemented tools on the S22-5 test set. eQE 2.0 delivers excellent interaction energies compared to conventional Kohn-Sham DFT and CCSD(T). The improved performance does not come at a loss of computational efficiency. We show that eQE 2.0 with nonlocal nonadditive functionals retains the same linear scaling behavior achieved in eQE 1.0 with semilocal nonadditive functionals.
rate research
Read More
The key feature of nonlocal kinetic energy functionals is their ability to reduce to the Thomas-Fermi functional in the regions of high density and to the von Weizsacker functional in the region of low density/high density gradient. This behavior is crucial when these functionals are employed in subsystem DFT simulations to approximate the nonadditive kinetic energy. We propose a GGA nonadditive kinetic energy functional which mimics the good behavior of nonlocal functionals retaining the computational complexity of typical semilocal functionals. The new functional reproduces Kohn-Sham DFT and benchmark CCSD(T) interaction energies of weakly interacting dimers in the S22-5 and S66 test sets with a mean absolute deviation well below 1 kcal/mol.
Several recent studies have shown that SCAN, a functional belonging to the meta-generalized gradient approximation (MGGA) family, leads to significantly overestimated magnetic moments in itinerant ferromagnetic metals. However, this behavior is not inherent to the MGGA level of approximation since TPSS, for instance, does not lead to such severe overestimations. In order to provide a broader view of the accuracy of MGGA functionals for magnetism, we extend the assessment to more functionals, but also to antiferromagnetic solids. The results show that to describe magnetism there is overall no real advantage in using a MGGA functional compared to GGAs. For both types of approximation, an improvement in ferromagnetic metals is necessarily accompanied by a deterioration (underestimation) in antiferromagnetic insulators, and vice-versa. We also provide some analysis in order to understand in more detail the relation between the mathematical form of the functionals and the results.
We present the One-orbital Ensemble Self-Consistent Field (OE-SCF) method, an {alternative} orbital-free DFT solver that extends the applicability of DFT to system sizes beyond the nanoscale while retaining the accuracy required to be predictive. OE-SCF is an iterative solver where the (typically computationally expensive) Pauli potential is treated as an external potential and updated after each iteration. Because only up to a dozen iterations are needed to reach convergence, OE-SCF dramatically outperforms current orbital-free DFT solvers. Employing merely a single CPU, we carried out the largest ab initio simulation for silicon-based materials to date. OE-SCF is able to converge the energy of bulk-cut Si nanoparticles as a function of their diameter up to 16 nm, for the first time reproducing known empirical results. We model polarization and interface charge transfer when a Si slab is sandwiched between two metal slabs where lattice matching mandates a very large slab size. Additionally, OE-SCF opens the door to adopt even more accurate functionals in orbital-free DFT simulations while still tackling systems sizes beyond the nanoscale.
The non-local van der Waals density functional (vdW-DF) has had tremendous success since its inception in 2004 due to its constraint-based formalism that is rigorously derived from a many-body starting point. However, while vdW-DF can describe binding energies and structures for van der Waals complexes and mixed systems with good accuracy, one long-standing criticism---also since its inception---has been that the $C_6$ coefficients that derive from the vdW-DF framework are largely inaccurate and can be wrong by more than a factor of two. It has long been thought that this failure to describe the $C_6$ coefficients is a conceptual flaw of the underlying plasmon framework used to derive vdW-DF. We prove here that this is not the case and that accurate $C_6$ coefficient can be obtained without sacrificing the accuracy at binding separations from a modified framework that is fully consistent with the constraints and design philosophy of the original vdW-DF formulation. Our design exploits a degree of freedom in the plasmon-dispersion model $omega_{mathbf{q}}$, modifying the strength of the long-range van der Waals interaction and the cross-over from long to short separations, with additional parameters tuned_ to reference systems. Testing the new formulation for a range of different systems, we not only confirm the greatly improved description of $C_6$ coefficients, but we also find excellent performance for molecular dimers and other systems. The importance of this development is not necessarily that particular aspects such as $C_6$ coefficients or binding energies are improved, but rather that our finding opens the door for further conceptual developments of an entirely unexplored direction within the exact same constrained-based non-local framework that made vdW-DF so successful in the first place.
We present a study of the optical response of compact and hollow icosahedral clusters containing up to 868 silver atoms by means of time-dependent density functional theory. We have studied the dependence on size and morphology of both the sharp plasmonic resonance at 3-4 eV (originated mainly from $sp$-electrons), and the less studied broader feature appearing in the 6-7 eV range (interband transitions). An analysis of the effect of structural relaxations, as well as the choice of exchange correlation functional (local density versus generalized gradient approximations) both in the ground state and optical response calculations is also presented. We have further analysed the role of the different atom layers (surface versus inner layers) and the different orbital symmetries on the absorption cross-section for energies up to 8 eV. We have also studied the dependence on the number of atom layers in hollow structures. Shells formed by a single layer of atoms show a pronounced red shift of the main plasmon resonances that, however, rapidly converge to those of the compact structures as the number of layers is increased. The methods used to obtain these results are also carefully discussed. Our methodology is based on the use of localized basis (atomic orbitals, and atom-centered- and dominant- product functions), which bring several computational advantages related to their relatively small size and the sparsity of the resulting matrices. Furthermore, the use of basis sets of atomic orbitals also brings the possibility to extend some of the standard population analysis tools (e.g., Mulliken population analysis) to the realm of optical excitations. Some examples of these analyses are described in the present work.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
