No Arabic abstract
The computation of excited electronic states with commonly employed (approximate) methods is challenging, typically yielding states of lower quality than the corresponding ground state for a higher computational cost. In this work, we present a mean field method that extends the previously proposed eXcited Constrained DFT (XCDFT) from single Slater determinants to ensemble 1-RDMs for computing low-lying excited states. The method still retains an associated computational complexity comparable to a semilocal DFT calculation while at the same time is capable of approaching states with multireference character. We benchmark the quality of this method on well-established test sets, finding good descriptions of the electronic structure of multireference states and maintaining an overall accuracy for the predicted excitation energies comparable to semilocal TDDFT.
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.
Lanthanide-based single-ion magnetic molecules can have large magnetic hyperfine interactions as well as large magnetic anisotropy. Recent experimental studies reported tunability of these properties by changes of chemical environments or by application of external stimuli for device applications. In order to provide insight onto the origin and mechanism of such tunability, here we investigate the magnetic hyperfine and nuclear quadrupole interactions for $^{159}$Tb nucleus in TbPc$_2$ (Pc=phthalocyanine) single-molecule magnets using multireference ab-initio methods including spin-orbit interaction. Since the electronic ground and first-excited (quasi)doublets are well separated in energy, the microscopic Hamiltonian can be mapped onto an effective Hamiltonian with an electronic pseudo-spin $S=1/2$. From the ab-initio-calculated parameters, we find that the magnetic hyperfine coupling is dominated by the interaction of the Tb nuclear spin with electronic orbital angular momentum. The asymmetric $4f$-like electronic charge distribution leads to a strong nuclear quadrupole interaction with significant non-axial terms for the molecule with low symmetry. The ab-initio calculated electronic-nuclear spectrum including the magnetic hyperfine and quadrupole interactions is in excellent agreement with experiment. We further find that the non-axial quadrupole interactions significantly influence the avoided level crossings in magnetization dynamics and that the molecular distortions affect mostly the Fermi contact terms as well as the non-axial quadrupole interactions.
Kohn-Sham density functional theory (DFT) has become established as an indispensable tool for investigating aqueous systems of all kinds, including those important in chemistry, surface science, biology and the earth sciences. Nevertheless, many widely used approximations for the exchange-correlation (XC) functional describe the properties of pure water systems with an accuracy that is not fully satisfactory. The explicit inclusion of dispersion interactions generally improves the description, but there remain large disagreements between the predictions of different dispersion-inclusive methods. We present here a review of DFT work on water clusters, ice structures and liquid water, with the aim of elucidating how the strengths and weaknesses of different XC approximations manifest themselves across this variety of water systems. Our review highlights the crucial role of dispersion in describing the delicate balance between compact and extended structures of many different water systems, including the liquid. By referring to a wide range of published work, we argue that the correct description of exchange-overlap interactions is also extremely important, so that the choice of semi-local or hybrid functional employed in dispersion-inclusive methods is crucial. The origins and consequences of beyond-2-body errors of approximate XC functionals are noted, and we also discuss the substantial differences between different representations of dispersion. We propose a simple numerical scoring system that rates the performance of different XC functionals in describing water systems, and we suggest possible future developments.
Graphite is the most widely used and among the most widely-studied anode materials for lithium-ion batteries. With increasing demands on lithium batteries to operate at lower temperatures and higher currents, it is crucial to understand lithium intercalation in graphite due to issues associated with lithium plating. Lithium intercalation into graphite has been extensively studied theoretically using density functional theory (DFT) calculations, complemented by experimental studies through X-ray diffraction, spectroscopy, optical imaging and other techniques. In this work, we present a first principles based model using DFT calculations, employing the BEEF-vdW as the exchange correlation functional, and Ising model to determine the phase transformations and subsequently, the thermodynamic intercalation potential diagram. We explore a configurational phase space of about 1 billion structures by accurately determining the important interactions for the Ising model. The BEEF-vdW exchange correlation functional employed accurately captures a range of interactions including vdW, covalent and ionic interactions. We incorporate phonon contributions at finite temperatures and configurational entropy to get high accuracy in free energy and potentials. We utilize the built-in error estimation capabilities of the BEEF-vdW exchange correlation functional and to develop a methodological framework for determining the uncertainty associated with DFT calculated phase diagrams and intercalation potentials. The framework also determines the confidence of each predicted stable phase. The confidence value of a phase can help us to identify regions of solid solutions and phase transformations accurately.
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.