Algorithm for Linear Response Functions at Finite Temperatures: Application to ESR spectrum of s=1/2 Antiferromagnet Cu benzoate

We introduce an efficient and numerically stable method for calculating linear response functions $\chi(\vec{q},\omega)$ of quantum systems at finite temperatures. The method is a combination of numerical solution of the time-dependent Schroedinger equation, random vector representation of trace, and Chebyshev polynomial expansion of Boltzmann operator. This method should be very useful for a wide range of strongly correlated quantum systems at finite temperatures. We present an application to the ESR spectrum of s=1/2 antiferromagnet Cu benzoate.Comment: 4 pages, 4 figure

Itinerant ferromagnetism in half-metallic CoS_2

We have investigated electronic and magnetic properties of the pyrite-type CoS_2 using the linearized muffin-tin orbital (LMTO) band method. We have obtained the ferromagnetic ground state with nearly half-metallic nature. The half-metallic stability is studied by using the fixed spin moment method. The non-negligible orbital magnetic moment of Co 3d electrons is obtained as $\mu_L = 0.06 \mu_B$ in the local spin density approximation (LSDA). The calculated ratio of the orbital to spin angular momenta / = 0.15 is consistent with experiment. The effect of the Coulomb correlation between Co 3d electrons is also explored with the LSDA + U method. The Coulomb correlation at Co sites is not so large, $U \lesssim 1$ eV, and so CoS_2 is possibly categorized as an itinerant ferromagnet. It is found that the observed electronic and magnetic behaviors of CoS_2 can be described better by the LSDA than by the LSDA + U.Comment: 4 pages, 3 postscript figure

Electron energy loss and inelastic x-ray scattering cross sections from time-dependent density-functional perturbation theory

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