No Arabic abstract
The Liouville-Lanczos approach to linear-response time-dependent density-functional theory is generalized so as to encompass electron energy-loss and inelastic X-ray scattering spectroscopies in periodic solids. The computation of virtual orbitals and the manipulation of large matrices are avoided by adopting a representation of response orbitals borrowed from (time-independent) density-functional perturbation theory and a suitable Lanczos recursion scheme. The latter allows the bulk of the numerical work to be performed at any given transferred momentum only once, for a whole extended frequency range. The numerical complexity of the method is thus greatly reduced, making the computation of the loss function over a wide frequency range at any given transferred momentum only slightly more expensive than a single standard ground-state calculation, and opening the way to computations for systems of unprecedented size and complexity. Our method is validated on the paradigmatic examples of bulk silicon and aluminum, for which both experimental and theoretical results already exist in the literature.
We present a time-dependent density functional theory (TDDFT) based approach to compute the light-matter couplings between two different manifolds of excited states relative to a common ground state. These quantities are the necessary ingredients to solve the Kramers--Heisenberg equation for resonant inelastic X-ray scattering (RIXS) and several other types of two-photon spectroscopies. The procedure is based on the pseudo-wavefunction approach, where TDDFT eigenstates are treated as a configuration interaction wavefunction with single excitations, and on the restricted energy window approach, where a manifold of excited states can be rigorously defined based on the energies of the occupied molecular orbitals involved in the excitation process. We illustrate the applicability of the method by calculating the 2p4d RIXS maps of three representative Ruthenium complexes and comparing them to experimental results. The method is able to accurately capture all the experimental features in all three complexes, with relative energies correct to within 0.6 eV at the cost of two independent TDDFT calculations.
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.
Using a super-operator formulation of linearized time-dependent density-functional theory, the dynamical polarizability of a system of interacting electrons is given a matrix continued-fraction representation whose coefficients can be obtained from the non-symmetric block-Lanczos method. The resulting algorithm allows for the calculation of the {em full spectrum} of a system with a computational workload which is only a few times larger than that needed for {em static} polarizabilities within time-independent density-functional perturbation theory. The method is demonstrated with the calculation of the spectrum of benzene, and prospects for its application to the large-scale calculation of optical spectra are discussed.
We present a benchmark of the density functional linear response calculation of NMR shieldings within the Gauge-Including Projector-Augmented-Wave method against all-electron Augmented-Plane-Wave$+$local-orbital and uncontracted Gaussian basis set results for NMR shieldings in molecular and solid state systems. In general, excellent agreement between the aforementioned methods is obtained. Scalar relativistic effects are shown to be quite large for nuclei in molecules in the deshielded limit. The small component makes up a substantial part of the relativistic corrections.
We examine the experimental and theoretical electron-energy loss spectra in 2$H$-Cu$_{0.2}$NbS$_2$ and find that the 1 eV plasmon in this material does not exhibit the regular positive quadratic plasmon dispersion that would be expected for a normal broad-parabolic-band system. Instead we find a nearly non-dispersing plasmon in the momentum-transfer range $q<0.35$ AA$^{-1}$. We argue that for a stoichiometric pure 2$H$-NbS$_2$ the dispersion relation is expected to have a negative slope as is the case for other transition-metal dichalcogenides. The presence of Cu impurities, required to stabilize the crystal growth, tends to shift the negative plasmon dispersion into a positive one, but the doping level in the current system is small enough to result in a nearly-non-dispersing plasmon. We conclude that a negative-slope plasmon dispersion is not connected with the existence of a charge-density-wave order in transition metal dichalcogenides.