No Arabic abstract
Density-functional theory (DFT) has revolutionized computational prediction of atomic-scale properties from first principles in physics, chemistry and materials science. Continuing development of new methods is necessary for accurate predictions of new classes of materials and properties, and for connecting to nano- and mesoscale properties using coarse-grained theories. JDFTx is a fully-featured open-source electronic DFT software designed specifically to facilitate rapid development of new theories, models and algorithms. Using an algebraic formulation as an abstraction layer, compact C++11 code automatically performs well on diverse hardware including GPUs. This code hosts the development of joint density-functional theory (JDFT) that combines electronic DFT with classical DFT and continuum models of liquids for first-principles calculations of solvated and electrochemical systems. In addition, the modular nature of the code makes it easy to extend and interface with, facilitating the development of multi-scale toolkits that connect to ab initio calculations, e.g. photo-excited carrier dynamics combining electron and phonon calculations with electromagnetic simulations.
A curious behavior of electron correlation energy is explored. Namely, the correlation energy is the energy that tends to drive the system toward that of the uniform electron gas. As such, the energy assumes its maximum value when a gradient of density is zero. As the gradient increases, the energy is diminished by a gradient suppressing factor, designed to attenuate the energy from its maximum value similar to the shape of a bell curve. Based on this behavior, we constructed a very simple mathematical formula that predicted the correlation energy of atoms and molecules. Combined with our proposed exchange energy functional, we calculated the correlation energies, the total energies, and the ionization energies of test atoms and molecules; and despite the unique simplicities, the functionals accuracies are in the top tier performance, competitive to the B3LYP, BLYP, PBE, TPSS, and M11. Therefore, we propose that, as guided by the simplicities and supported by the accuracies, the correlation energy is the energy that locally tends to drive the system toward the uniform electron gas.
Octahedral Fe$^{2+}$ molecules are particularly interesting as they often exhibit a spin-crossover transition. In spite of the many efforts aimed at assessing the performances of density functional theory for such systems, an exchange-correlation functional able to account accurately for the energetic of the various possible spin-states has not been identified yet. Here we critically discuss the issues related to the theoretical description of this class of molecules from first principles. In particular we present a comparison between different density functionals for four ions, namely [Fe(H$_2$O)$_6$]$^{2+}$, [Fe(NH$_3$)$_6$]$^{2+}$, [Fe(NCH)$_6$]$^{2+}$ and [Fe(CO)$_6$]$^{2+}$. These are characterized by different ligand-field splittings and ground state spin multiplicities. Since no experimental data are available for the gas phase, the density functional theory results are benchmarked against those obtained with diffusion Monte Carlo, one of the most accurate methods available to compute ground state total energies of quantum systems. On the one hand, we show that most of the functionals considered provide a good description of the geometry and of the shape of the potential energy surfaces. On the other hand, the same functionals fail badly in predicting the energy differences between the various spin states. In the case of [Fe(H$_2$O)$_6$]$^{2+}$, [Fe(NH$_3$)$_6$]$^{2+}$, [Fe(NCH)$_6$]$^{2+}$, this failure is related to the drastic underestimation of the exchange energy. Therefore quite accurate results can be achieved with hybrid functionals including about 50% of Hartree-Fock exchange. In contrast, in the case of [Fe(CO)$_6$]$^{2+}$, the failure is likely to be caused by the multiconfigurational character of the ground state wave-function and no suitable exchange and correlation functional has been identified.
Time-dependent density-functional theory (TDDFT) is a computationally efficient first-principles approach for calculating optical spectra in insulators and semiconductors, including excitonic effects. We show how exciton wave functions can be obtained from TDDFT via the Kohn-Sham transition density matrix, both in the frequency-dependent linear-response regime and in real-time propagation. The method is illustrated using one-dimensional model solids. In particular, we show that our approach provides insight into the formation and dissociation of excitons in real time. This opens the door to time-resolved studies of exciton dynamics in materials by means of real-time TDDFT.
We present accurate optical spectra of semiconductors and insulators within a pure Kohn-Sham time-dependent density-functional approach. In particular, we show that the onset of the absorption is well reproduced when comparing to experiment. No empirical information nor a theory beyond Kohn-Sham density-functional theory, such as $GW$, is invoked to correct the Kohn-Sham gap. Our approach relies on the link between the exchange-correlation kernel of time-dependent density functional theory and the derivative discontinuity of ground-state density-functional theory. We show explicitly how to relate these two quantities. We illustrate the accuracy and simplicity of our approach by applying it to various semiconductors (Si, GaP, GaAs) and wide-gap insulators (C, LiF, Ar).
Exchange interactions are a manifestation of the quantum mechanical nature of the electrons and play a key role in predicting the properties of materials from first principles. In density functional theory (DFT), a widely used approximation to the exchange energy combines fractions of density-based and Hartree-Fock (exact) exchange. This so-called hybrid DFT scheme is accurate in many materials, for reasons that are not fully understood. Here we show that a 1/4 fraction of exact exchange plus a 3/4 fraction of density-based exchange is compatible with a correct quantum mechanical treatment of the exchange energy of an electron pair in the unpolarized electron gas. We also show that the 1/4 exact-exchange fraction mimics a correlation interaction between doubly-excited electronic configurations. The relation between our results and trends observed in hybrid DFT calculations is discussed, along with other implications.