No Arabic abstract
It is generally accepted that the effective magnetic field acting on a magnetic moment is given by the gradient of the energy with respect to the magnetization. However, in ab initio spin dynamics within the adiabatic approximation, the effective field is also known to be exactly the negative of the constraining field, which acts as a Lagrange multiplier to stabilize an out-of-equilibrium, non-collinear magnetic configuration. We show that for Hamiltonians without mean-field parameters both of these fields are exactly equivalent, while there can be a finite difference for mean-field Hamiltonians. For density-functional theory (DFT) calculations the constraining field obtained from the auxiliary Kohn-Sham Hamiltonian is not exactly equivalent to the DFT energy gradient. This inequality is highly relevant for both ab initio spin dynamics and the ab initio calculation of exchange constants and effective magnetic Hamiltonians. We argue that the effective magnetic field and exchange constants have the highest accuracy in DFT when calculated from the energy gradient and not from the constraining field.
Spin relaxation and decoherence is at the heart of spintronics and spin-based quantum information science. Currently, theoretical approaches that can accurately predict spin relaxation of general solids including necessary scattering pathways and capable for ns to ms simulation time are urgently needed. We present a first-principles real-time density-matrix approach based on Lindblad dynamics to simulate ultrafast spin dynamics for general solid-state systems. Through the complete first-principles descriptions of pump, probe and scattering processes including electron-phonon, electron-impurity and electron-electron scatterings with self-consistent spin-orbit couplings, our method can directly simulate the ultrafast pump-probe measurements for coupled spin and electron dynamics over ns at any temperature and doping levels. We apply this method to a prototypical system GaAs and obtain excellent agreement with experiments. We found that the relative contributions of different scattering mechanisms and phonon modes differ considerably between spin and carrier relaxation processes. In sharp contrast to previous work based on model Hamiltonians, we point out that the electron-electron scattering is negligible at room temperature but becomes very important at low temperatures for spin relaxation in n-type GaAs. Most importantly, we examine the applicable conditions of the commonly-used Dyakonov-Perel relation, which may break down for individual scattering processes. Our work provides a predictive computational platform for spin relaxation in solids, which has unprecedented potentials for designing new materials ideal for spintronics and quantum information technology.
The anomalous plasmon linewidth dispersion (PLD) measured in K by vom Felde, Sprosser-Prou, and Fink (Phys. Rev. B 40, 10181 (1989)), has been attributed to strong dynamical electron-electron correlations. On the basis of ab initio response calculations, and detailed comparison with experiment, we show that the PLD of K is, in fact, dominated by decay into particle-hole excitations involving empty states of d-symmetry. For Li, we shed new light on the physics of the PLD. Our all-electron results illustrate the importance of ab initio methods for the study of electronic excitations.
The effects of tetragonal strain on electronic and magnetic properties of strontium-doped lanthanum manganite, La_{2/3}Sr_{1/3}MnO_3 (LSMO), are investigated by means of density-functional methods. As far as the structural properties are concerned, the comparison between theory and experiments for LSMO strained on the most commonly used substrates, shows an overall good agreement: the slight overestimate (at most of 1-1.5 %) for the equilibrium out-of-plane lattice constants points to possible defects in real samples. The inclusion of a Hubbard-like contribution on the Mn d states, according to the so-called LSDA+U approach, is rather ineffective from the structural point of view, but much more important from the electronic and magnetic point of view. In particular, full half-metallicity, which is missed within a bare density-functional approach, is recovered within LSDA+U, in agreement with experiments. Moreover, the half-metallic behavior, particularly relevant for spin-injection purposes, is independent on the chosen substrate and is achieved for all the considered in-plane lattice constants. More generally, strain effects are not seen to crucially affect the electronic structure: within the considered tetragonalization range, the minority gap is only slightly (i.e. by about 0.1-0.2 eV) affected by a tensile or compressive strain. Nevertheless, we show that the growth on a smaller in-plane lattice constant can stabilize the out-of-plane vs in-plane e_g orbital and significatively change their relative occupancy. Since e_g orbitals are key quantities for the double-exchange mechanism, strain effects are confirmed to be crucial for the resulting magnetic coupling.
First-principles calculation predict that olivine Li4MnFeCoNiP4O16 has ferrotoroidic characteristic and ferrimagnetic configuration with magnetic moment of 1.56 muB per formula unit. The ferrotoroidicity of this material makes it a potential candidate for magnetoelectric materials . Based on the orbital-resolved density of states for the transtion-metal ions in Li4MnFeCoNiP4O16, the spin configuration for Mn2+,Fe3+,Co2+, and Ni2+ is t2g3eg2, t2g3eg2,t2g1t2g3eg1eg2, and t2g2t2g3eg1eg2, respectively. Density functional theory plus U (DFT+U) shows a indirect band gap of 1.25 eV in this predicted material, which is not simply related to the electronic conductivity in terms of being used as cathode material in rechargeable Li-ion batteries.
We report ab initio calculations of the melting curve of molybdenum for the pressure range 0-400 GPa. The calculations employ density functional theory (DFT) with the Perdew-Burke-Ernzerhof exchange-correlation functional in the projector augmented wave (PAW) implementation. We present tests showing that these techniques accurately reproduce experimental data on low-temperature b.c.c. Mo, and that PAW agrees closely with results from the full-potential linearized augmented plane-wave implementation. The work attempts to overcome the uncertainties inherent in earlier DFT calculations of the melting curve of Mo, by using the ``reference coexistence technique to determine the melting curve. In this technique, an empirical reference model (here, the embedded-atom model) is accurately fitted to DFT molecular dynamics data on the liquid and the high-temperature solid, the melting curve of the reference model is determined by simulations of coexisting solid and liquid, and the ab initio melting curve is obtained by applying free-energy corrections. Our calculated melting curve agrees well with experiment at ambient pressure and is consistent with shock data at high pressure, but does not agree with the high pressure melting curve deduced from static compression experiments. Calculated results for the radial distribution function show that the short-range atomic order of the liquid is very similar to that of the high-T solid, with a slight decrease of coordination number on passing from solid to liquid. The electronic densities of states in the two phases show only small differences. The results do not support a recent theory according to which very low dTm/dP values are expected for b.c.c. transition metals because of electron redistribution between s-p and d states.