Search

We present an ab initio method for eletcronic structure calculations, which accounts for the interaction of electrons and magnons in ferromagnets. While it is based on a many body perturbation theory we approximate numerically complex quantities with quantities from time dependent density functional theory. This results in a simple and affordable algorithm which allows us to consider more complex materials than those usually studied in this context ($3d$ ferromagnets) while still being able to account for the non-locality of the self energy. Furthermore, our approach allows for a relatively simple way to incorporate self-consistency. Our results are in a good agreement with experimental and theoretical findings for iron and nickel. Especially the experimental exchange splitting of nickel is predicted accurately within our theory. Additionally, we study the halfmetallic ferromagnet NiMnSb concerning its non-qusiparticle states appearing in the bandgap due to spin-flip excitations.

We apply quantum continuum mechanics to the calculation of the excitation spectrum of a coupled electron-hole bilayer. The theory expresses excitation energies in terms of ground-state intra- and inter-layer pair correlation functions, which are available from Quantum Monte Carlo calculations. The final formulas for the collective modes deduced from this approach coincide with the formulas obtained in the "quasi-localized particle approximation" by Kalman et al., and likewise, the theory predicts the existence of gapped excitations in the charged channels, with the gap arising from electron-hole correlation. An immediate consequence of the gap is that the static density-density response function of the charged channel vanishes as $q^2$ for wave vector $q \to 0$, rather than linearly in $q$, as commonly expected. In this sense, the system is {\it incompressible}. This feature, which has no analogue in the classical electron-hole plasma, is consistent with the existence of an excitonic ground state and implies the existence of a discontinuity in the chemical potential of electrons and holes when the numbers of electrons and holes are equal. It should be experimentally observable by monitoring the densities of electrons and holes in response to potentials that attempt to change these densities in opposite directions., Comment: 14 pages, 7 figures

Recent experiments on twisted bilayer graphene (tBG) close to magic angle show that a small relative rotation in a van der Waals heterostructure greatly alters its electronic properties. We consider various scattering mechanisms and show that the carrier transport in tBG is dominated by a combination of charged impurities and acoustic gauge phonons. Charged impurities still dominate at low temperature and densities because of the inability of Dirac fermions to screen long-range Coulomb potentials at charge neutrality; however, the gauge phonons dominate for most of the experimental regime because although they couple to current, they do not induce charge and are therefore unscreened by the large density of states close to magic angle. We show that the resistivity has a strong monotonically decreasing carrier density dependence at low temperature due to charged impurity scattering, and weak density dependence at high temperature due to gauge phonons. Away from charge neutrality, the resistivity increases with temperature, while it does the opposite close to the Dirac point. A non-monotonic temperature dependence observed only at low temperature and carrier density is a signature of our theory that can be tested in experimentally available samples., Comment: 5 pages, 4 figures

Ab initio calculations of the magnon dispersion in ferromagnetic materials typically rely on the adiabatic local density approximation (ALDA) in which the effective exchange-correlation field is everywhere parallel to the magnetization. These calculations, however, tend to overestimate the "magnon stiffness", defined as the curvature of the magnon frequency vs. wave vector relation evaluated at zero wave vector. Here we suggest a simple procedure to improve the magnon dispersion by taking into account gradient corrections to the ALDA at the exchange-only level. We find that this gradient correction always reduces the magnon stiffness. The surprisingly large size of these corrections ($\sim 30\%$) greatly improves the agreement between the calculated and the observed magnon stiffness for cobalt and nickel, which are known to be overestimated within the ALDA., Comment: 6 pages, 4 figures

A semi-relativistic density-functional theory that includes spin-orbit couplings and Zeeman fields on equal footing with the electromagnetic potentials, is an appealing framework to develop a unified first-principles computational approach for non-collinear magnetism, spintronics, orbitronics, and topological states. The basic variables of this theory include the paramagnetic current and the spin-current density, besides the particle and the spin density, and the corresponding exchange-correlation (xc) energy functional is invariant under local U(1)$\times$SU(2) gauge transformations. The xc-energy functional must be approximated to enable practical applications, but, contrary to the case of the standard density functional theory, finding simple approximations suited to deal with realistic atomistic inhomogeneities has been a long-standing challenge. Here, we propose a way out of this impasse by showing that approximate gauge-invariant functionals can be easily generated from existing approximate functionals of ordinary density-functional theory by applying a simple {\it minimal substitution} on the kinetic energy density, which controls the short-range behavior of the exchange hole. Our proposal opens the way to the construction of approximate, yet non-empirical functionals, which do not assume weak inhomogeneity and should therefore have a wide range of applicability in atomic, molecular and condensed matter physics.

We discuss the potential advantages of calculating the effective mass of quasiparticles in the interacting electron liquid from the low-temperature free energy vis-a-vis the conventional approach, in which the effective mass is obtained from approximate calculations of the self-energy, or from a quantum Monte Carlo evaluation of the energy of a variational "quasiparticle wave function". While raw quantum Monte Carlo data are presently too sparse to allow for an accurate determination of the effective mass, the values estimated by this method are numerically close to the ones obtained in previous calculations using diagrammatic many-body theory. In contrast to this, a recently published parametrization of quantum Monte Carlo data for the free energy of the homogeneous electron liquid yields effective masses that considerably deviate from previous calculations and even change sign for low densities, reflecting an unphysical negative entropy. We suggest that this anomaly is related to the treatment of the exchange energy at finite temperature., Comment: 7 pages, 4 figures

We review the progress that has been recently made in the application of time-dependent density functional theory to thermoelectric phenomena. As the field is very young, we emphasize open problems and fundamental issues. We begin by introducing the formal structure of \emph{thermal density functional theory}, a density functional theory with two basic variables -- the density and the energy density -- and two conjugate fields -- the ordinary scalar potential and Luttinger's thermomechanical potential. The static version of this theory is contrasted with the familiar finite-temperature density functional theory, in which only the density is a variable. We then proceed to constructing the full time-dependent non equilibrium theory, including the practically important Kohn-Sham equations that go with it. The theory is shown to recover standard results of the Landauer theory for thermal transport in the steady state, while showing greater flexibility by allowing a description of fast thermal response, temperature oscillations and related phenomena. Several results are presented here for the first time, i.e., the proof of invertibility of the thermal response function in the linear regime, the full expression of the thermal currents in the presence of Luttinger's thermomechanical potential, an explicit prescription for the evaluation of the Kohn-Sham potentials in the adiabatic local density approximation, a detailed discussion of the leading dissipative corrections to the adiabatic local density approximation and the thermal corrections to the resistivity that follow from it., Comment: 34 pages, 6 figures

Quantized spin waves, or magnons, in a magnetic insulator are assumed to interact weakly with the surroundings, and to flow with little dissipation or drag, producing exceptionally long diffusion lengths and relaxation times. In analogy to Coulomb drag in bilayer two dimensional electron gases, in which the contribution of the Coulomb interaction to the electric resistivity is studied by measuring the interlayer resistivity (transresistivity), we predict a nonlocal drag of magnons in a ferromagnetic bilayer structure based on semiclassical Boltzmann equations. Nonlocal magnon drag depends on magnetic dipolar interactions between the layers and manifests in the magnon current transresistivity and the magnon thermal transresistivity, whereby a magnon current in one layer induces a chemical potential gradient and/or a temperature gradient in the other layer. The largest drag effect occurs when the magnon current flows parallel to the magnetization, however for oblique magnon currents a large transverse current of magnons emerges. We examine the effect for practical parameters, and find that the predicted induced temperature gradient is readily observable., Comment: 5 pages, 3 figures, PRL in press

We analyze the short-time behavior of the heat and charge currents through nanoscale conductors exposed to a temperature gradient. To this end, we employ Luttinger's thermomechanical potential to simulate a sudden change of temperature at one end of the conductor. We find that the direction of the charge current through an impurity is initially opposite to the direction of the charge current in the steady-state limit. Furthermore, we investigate the transient propagation of energy and particle density driven by a temperature variation through a conducting nanowire. Interestingly, we find that the velocity of the wavefronts of, both, the particle and the energy wave have the same constant value, insensitive to changes in the average electronic density. In the steady-state regime, we find that, at low temperatures, the local temperature and potential, as measured by a floating probe lead, exhibit characteristic oscillations due to quantum interference, with a periodicity that corresponds to half the Fermi wavelength of the electrons., Comment: 10 pages, 11 figures

We study spin relaxation in dilute magnetic semiconductors near a ferromagnetic transition, where spin fluctuations become strong. An enhancement in the scattering rate of itinerant carriers from the spin fluctuations of localized impurities leads to a change in the dominant spin relaxation mechanism from Dyakonov-Perel to spin flips in scattering. On the ferromagnetic side of the transition, we show that due to the presence of two magnetic components -- the itinerant carriers and the magnetic impurities -- with different gyromagnetic ratios, the relaxation rate of the total magnetization can be quite different from the relaxation rate of the spin. Following a disturbance of the equilibrium magnetization, the spin is initially redistributed between the two components to restore the equilibrium magnetization. It is only on a longer time scale, controlled by the spin-orbit interaction, that the total spin itself relaxes to its equilibrium state., Comment: 12 pages, 5 figures

The localized Hartree-Fock potential has proven to be a computationally efficient alternative to the optimized effective potential, preserving the numerical accuracy of the latter and respecting the exact properties of being self-interaction free and having the correct $-1/r$ asymptotics. In this paper we extend the localized Hartree-Fock potential to fractional particle numbers and observe that it yields derivative discontinuities in the energy as required by the exact theory. The discontinuities are numerically close to those of the computationally more demanding Hartree-Fock method. Our potential enjoys a "direct-energy" property, whereby the energy of the system is given by the sum of the single-particle eigenvalues multiplied by the corresponding occupation numbers. The discontinuities $c_\uparrow$ and $c_\downarrow$ of the spin-components of the potential at integer particle numbers $N_\uparrow$ and $N_\downarrow$ satisfy the condition $c_\uparrow N_\uparrow+c_\downarrow N_\downarrow=0$. Thus, joining the family of effective potentials which support a derivative discontinuity, but being considerably easier to implement, the localized Hartree-Fock potential becomes a powerful tool in the broad area of applications in which the fundamental gap is an issue., Comment: 11 pages, 6 figures

Magnetic fluctuations in transition metal oxides are a subject of intensive research because of the key role they are expected to play in the transition from the Mott insulator to the unconventional metallic phase of these materials, and also as drivers of superconductivity. Despite much effort, a clear link between magnetic fluctuations and the insulator-to-metal transition has not yet been established. Here we report the discovery of a compelling link between magnetic fluctuations and the insulator-to-metal transition in Ca(Ir$_{1-x}$Ru$_{x}$)O$_{3}$ perovskites as a function of the doping coefficient x. We show that when the material turns from insulator to metal, at a critical value of x$\sim$ 0.3, magnetic fluctuations change their character from antiferromagnetic, a Mott insulator phase, to ferromagnetic, an itinerant electron state with Hund's orbital coupling. These results are expected to have wide-ranging implications for our understanding of the unconventional properties of strongly correlated electrons systems, Comment: 7 pages, 5 figures

We analyze the effect known as "spin current swapping" (SCS) due to electron-impurity scattering in a uniform spin-polarized two-dimensional electron gas. In this effect a primary spin current $J_i^a$ (lower index for spatial direction, upper index for spin direction) generates a secondary spin current $J_a^i$ if $i \neq a$, or $J_j^j$, with $j\ne i$, if $i= a$. Contrary to naive expectation, the homogeneous spin current associated with the uniform drift of the spin polarization in the electron gas does not generate a swapped spin current by the SCS mechanism. Nevertheless, a swapped spin current will be generated, if a magnetic field is present, by a completely different mechanism, namely, the precession of the spin Hall spin current in the magnetic field. We refer to this second mechanism as Hanle spin Hall effect, and we notice that it can be observed in an experiment in which a homogeneous drift current is passed through a uniformly magnetized electron gas. In contrast to this, we show that an unambiguous observation of SCS requires inhomogeneous spin currents, such as those that are associated with spin diffusion in a metal, and no magnetic field. An experimental setup for the observation of the SCS is therefore proposed., Comment: 8 pages, 5 figures

The broken inversion symmetry at the surface of a metallic film (or, more generally, at the interface between a metallic film and a different metallic or insulating material) greatly amplifies the influence of the spin-orbit interaction on the surface properties. The best known manifestation of this effect is the momentum-dependent splitting of the surface state energies (Rashba effect). Here we show that the same interaction also generates a spin-polarization of the bulk states when an electric current is driven through the bulk of the film. For a semi-infinite jellium model, which is representative of metals with a closed Fermi surface, we prove as a theorem that, regardless of the shape of the confinement potential, the induced surface spin density at each surface is given by ${\bf S} =-\gamma \hbar {\bf \hat z}\times {\bf j}$, where ${\bf j}$ is the particle current density in the bulk, ${\bf \hat z}$ the unit vector normal to the surface, and $\gamma=\frac{\hbar}{4mc^2}$ contains only fundamental constants. For a general metallic solid $\gamma$ becomes a material-specific parameter that controls the strength of the interfacial spin-orbit coupling. Our theorem, combined with an {\it ab initio} calculation of the spin polarization of the current-carrying film, enables a determination of $\gamma$, which should be useful in modeling the spin-dependent scattering of quasiparticles at the interface., Comment: 5 pages, 2 figures

The spin transport and spin-to-charge current conversion properties of bismuth are investigated using permalloy/copper/bismuth (Py/Cu/Bi)lateral spin valve structures. The spin current is strongly absorbed at the surface of Bi, leading to ultra-short spin diffusion lengths. A spin-to-charge current conversion is measured, which is attributed to the inverse Rashba-Edelstein effect at the Cu/Bi interface. The spin-current-induced charge current is found to change direction with increasing temperature. A theoretical analysis relates this behavior to the complex spin structure and dispersion of the surface states at the Fermi energy. The understanding of this phenomenon opens novel possibilities to exploit spin-orbit coupling to create, manipulate, and detect spin currents in 2D systems., Comment: 10 pages, 4 figures for the main text + 6 pages, 6 figures for the supplementary information

Spin-orbit interactions in two-dimensional electron liquids are responsible for many interesting transport phenomena in which particle currents are converted to spin polarizations and spin currents and viceversa. Prime examples are the spin Hall effect, the Edelstein effect, and their inverses. By similar mechanisms it is also possible to partially convert an optically induced electron-hole density wave to a spin density wave and viceversa. In this paper we present a unified theoretical treatment of these effects based on quantum kinetic equations that include not only the intrinsic spin-orbit coupling from the band structure of the host material, but also the spin-orbit coupling due to an external electric field and a random impurity potential. The drift-diffusion equations we derive in the diffusive regime are applicable to a broad variety of experimental situations, both homogeneous and non-homogeneous, and include on equal footing "skew scattering" and "side-jump" from electron-impurity collisions. As a demonstration of the strength and usefulness of the theory we apply it to the study of several effects of current experimental interest: the inverse Edelstein effect, the spin-current swapping effect, and the partial conversion of an electron-hole density wave to a spin density wave in a two-dimensional electron gas with Rashba and Dresselhaus spin-orbit couplings, subject to an electric field., Comment: 19 pages, 6 figures

We investigate the effects of inhomogeneities on spin entanglement in many-electron systems from an ab-initio approach. The key quantity in our approach is the local spin entanglement length, which is derived from the local concurrence of the electronic system. Although the concurrence for an interacting systems is a highly nonlocal functional of the density, it does have a simple, albeit approximate expression in terms of Kohn-Sham orbitals. We show that the electron localization function -- well known in quantum chemistry as a descriptor of atomic shells and molecular bonds -- can be reinterpreted in terms of the ratio of the local entanglement length of the inhomogeneous system to the entanglement length of a homogenous system at the same density. We find that the spin entanglement is remarkably enhanced in atomic shells and molecular bonds.

Thermoelectric transport in nanoscale conductors is analyzed in terms of the response of the system to a thermo-mechanical field, first introduced by Luttinger, which couples to the electronic energy density. While in this approach the temperature remains spatially uniform, we show that a spatially varying thermo-mechanical field effectively simulates a temperature gradient across the system and allows us to calculate the electric and thermal currents that flow due to the thermo-mechanical field. In particular, we show that, in the long-time limit, the currents thus calculated reduce to those that one obtains from the Landauer-B\"uttiker formula, suitably generalized to allow for different temperatures in the reservoirs, if the thermo-mechanical field is applied to prepare the system, and subsequently turned off at ${t=0}$. Alternatively, we can drive the system out of equilibrium by switching the thermo-mechanical field after the initial preparation. We compare these two scenarios, employing a model noninteracting Hamiltonian, in the linear regime, in which they coincide, and in the nonlinear regime in which they show marked differences. We also show how an operationally defined local effective temperature can be computed within this formalism., Comment: 17 pages, 9 figures

The Wiedemann-Franz law, connecting the electronic thermal conductivity to the electrical conductivity of a disordered metal, is generally found to be well satisfied even when electron-electron (e-e) interactions are strong. In ultra-clean conductors, however, large deviations from the standard form of the law are expected, due to the fact that e-e interactions affect the two conductivities in radically different ways. Thus, the standard Wiedemann-Franz ratio between the thermal and the electric conductivity is reduced by a factor $1+\tau/\tau_{\rm th}^{\rm ee}$, where $1/\tau$ is the momentum relaxation rate, and $1/\tau_{\rm th}^{\rm ee}$ is the relaxation time of the thermal current due to e-e collisions. Here we study the density and temperature dependence of $1/\tau_{\rm th}^{\rm ee}$ in the important case of doped, clean single layers of graphene, which exhibit record-high thermal conductivities. We show that at low temperature $1/\tau_{\rm th}^{\rm ee}$ is $8/5$ of the quasiparticle decay rate. We also show that the many-body renormalization of the thermal Drude weight coincides with that of the Fermi velocity., Comment: 6 pages, 5 appendices (13 pages)

A normal metallic film sandwiched between two insulators may have strong spin-orbit coupling near the metal-insulator interfaces, even if spin-orbit coupling is negligible in the bulk of the film. In this paper we study two technologically important and deeply interconnected effects that arise from interfacial spin-orbit coupling in metallic films. The first is the spin Hall effect, whereby a charge current in the plane of the film is partially converted into an orthogonal spin current in the same plane. The second is the Edelstein effect, in which a charge current produces an in-plane, transverse spin polarization. At variance with strictly two-dimensional Rashba systems, we find that the spin Hall conductivity has a finite value even if spin-orbit interaction with impurities is neglected and "vertex corrections" are properly taken into account. Even more remarkably, such finite value becomes "universal" in a certain configuration. This is a direct consequence of the spatial dependence of spin-orbit coupling on the third dimension, perpendicular to the film plane. The non-vanishing spin Hall conductivity has a profound influence on the Edelstein effect, which we show to consist of two terms, the first with the standard form valid in a strictly two-dimensional Rashba system, and a second arising from the presence of the third dimension. Whereas the standard term is proportional to the momentum relaxation time, the new one scales with the spin relaxation time. Our results, although derived in a specific model, should be valid rather generally, whenever a spatially dependent Rashba spin-orbit coupling is present and the electron motion is not strictly two-dimensional., Comment: 23 pages, 3 figures

We study the equilibrium current density profiles of harmonically trapped ultra-cold Fermi gases in quantum Hall-like states that appear when the quasi-two-dimensional trap is set in fast rotation. The density profile of the gas (in the rotating reference frame) consists of incompressible strips of constant quantized density separated by compressible regions in which the density varies. Remarkably, we find that the atomic currents flow in opposite directions in the compressible and incompressible regions -- a prediction that should be amenable to experimental verification., Comment: 5 pages, 2 figures

We provide a precise microscopic definition of the recently observed "Inverse Edelstein Effect" (IEE), in which a non-equilibrium spin accumulation in the plane of a two-dimensional (interfacial) electron gas drives an electric current perpendicular to its own direction. The drift-diffusion equations that govern the effect are presented and applied to the interpretation of the experiments., Comment: 5 pages, 1 figure

We derive the gradient expansion for the exchange energy of a spin-polarized electron gas by perturbing the uniformly spin polarized state and thus inducing a small non-collinearity that is slowly varying in space. We show that the exchange-energy contribution due to the induced longitudinal gradient of the spin polarization to the exchange energy differs from the contribution due to the transverse gradient. The difference is present at any non-vanishing spin polarization and becomes larger with increasing spin polarization. We argue that improved generalized gradient approximations of Spin-Density-Functional Theory must account for the difference between the longitudinal and transverse spin stiffness., Comment: 9 pages, 6 figures

The charge density relaxation propagator of a two dimensional electron system, which is the slope of the imaginary part of the polarization function, exhibits singularities for bosonic momenta having the order of the spin-orbit momentum and depending on the momentum orientation. We have provided an intuitive understanding for this non-analytic behavior in terms of the inter chirality subband electronic transitions, induced by the combined action of Bychkov-Rashba (BR) and Dresselhaus (D) spin-orbit coupling. It is shown that the regular behavior of the relaxation propagator is recovered in the presence of only one BR or D spin-orbit field or for spin-orbit interaction with equal BR and D coupling strengths. This creates a new possibility to influence carrier relaxation properties by means of an applied electric field., Comment: 4 figures

We introduce a non-equilibrium density-functional theory of local temperature and associated local energy density that is suited for the study of thermoelectric phenomena. The theory rests on a local temperature field coupled to the energy-density operator. We identify the excess-energy density, in addition to the particle density, as the basic variable, which is reproduced by an effective noninteracting Kohn-Sham system. A novel Kohn-Sham equation emerges featuring a time-dependent and spatially varying mass which represents local temperature variations. The adiabatic contribution to the Kohn-Sham potentials is related to the entropy viewed as a functional of the particle and energy density. Dissipation can be taken into account by employing linear response theory and the thermoelectric transport coefficients of the electron gas., Comment: Paper: 5 pages, 1 figure, Supplemental Material: 3 pages

We present a unified Boltzmann-transport theory for the drag resistivity in two-component systems close to a second-order phase transition. We find general expressions for the drag resistivity in two and three spatial dimensions, for arbitrary population and mass imbalance, for particle- and hole-like bands, and show how to incorporate, at the Gaussian level, the effect of fluctuations close to a phase transition. We find that the proximity to the phase transition enhances the drag resistivity upon approaching the critical temperature from above, and we qualitatively derive the temperature dependence of this enhancement for various cases. In addition, we present numerical results for two concrete experimental systems: i) three-dimensional cold atomic Fermi gases close to a Stoner transition and ii) two-dimensional spatially-separated electron and hole systems in semiconductor double quantum wells., Comment: 15 pages, 5 figures

We show that an electric field parallel to the wavefronts of an electron-hole grating in a GaAs quantum well generates, via the electronic spin Hall effect, a spin grating of the same wave vector and with an amplitude that can exceed 1% of the amplitude of the initial density grating. We refer to this phenomenon as "collective spin Hall effect". A detailed study of the coupled-spin charge dynamics for quantum wells grown in different directions reveals rich features in the time evolution of the induced spin density, including the possibility of generating a helical spin grating., Comment: 8 pages, 3 figures

An asymmetric triangular potential well provides the simplest model for the confinement of mobile electrons at the interface between two insulating oxides, such as LaAlO_3 and SrTiO_3 (LAO/STO). These electrons have been recently shown to exhibit a large spin-orbit coupling of the Rashba type, i.e., linear in the in-plane momentum. In this paper we study the intrinsic spin Hall effect due to Rashba coupling in an asymmetric triangular potential well. Besides splitting each subband into two branches of opposite helicity, the spin-orbit interaction causes the wave function in the direction perpendicular to the plane of the quantum well (the z direction) to depend on the in-plane wave vector k. In contrast to the extreme asymmetric case, i.e., the wedge-shaped quantum well, for which the intrinsic spin Hall effect vanishes due to vertex corrections, we find that the asymmetric well supports a non-vanishing intrinsic spin Hall conductivity, which increases in magnitude as the well becomes more symmetric. The spin Hall conductivity is found to be proportional to the square of the spin-orbit coupling constant and, in the limit of low carrier density, depends only on the effective mass renormalization associated with the k-dependence of the wave functions in the z direction. Its origin lies in the non-vanishing matrix elements of the spin current between subbands corresponding to different states of quantized motion perpendicular to the plane of the well., Comment: 22 pages, 5 figures

We make use of the time-dependent density functional theory to derive a new formally exact expression for the dc resistivity of metals with impurities. This expression takes fully into account the dynamics of electron-electron interactions. Correction to the conventional $T$-matrix (phase-shifts) theory is treated within hydrodynamics of inhomogeneous viscous electron liquid. As a first application, we present calculations of the residual resistivity of aluminum as a function of the atomic number of the impurities. We show that the inclusion of many-body corrections considerably improves the agreement between theory and experiment., Comment: 4 pages, 1 figure

Phase-resolved transient grating spectroscopy in semiconductor quantum wells has been shown to be a powerful technique for measuring the electron-hole drag resistivity $\rho_{eh}$, which depends on the Coulomb interaction between the carriers. In this paper we develop the interacting drift-diffusion theory, from which $\rho_{eh}$ can be determined, given the measured mobility of an electron-hole grating. From this theory we predict a cross-over from a high-excitation-density regime, in which the mobility has the "normal" positive value, to a low-density regime, in which Coulomb-drag dominates and the mobility becomes negative. At the crossover point, the mobility of the grating vanishes., Comment: 5 pages, 5 figures

This comment criticizes the above paper by Xiao-Yin Pan and Viraht Sahni. It is shown that their formulation of "Physical Current Density Functional Theory" is, at best, a garbled reformulation of the Vignale-Rasolt current-density functional theory, and, at worst, a potential source of mistakes insofar as it complicates the formulation of the variational principle and prevents the constrained search construction of the universal functional., Comment: 5 pages

We consider spin transport in a two-component ultracold Fermi gas with attractive interspecies interactions close to the BCS pairing transition. In particular, we consider the spin-transport relaxation rate and the spin-diffusion constant. Upon approaching the transition, the scattering amplitude is enhanced by pairing fluctuations. However, as the system approaches the transition, the spectral weight for excitations close to the Fermi level is decreased by the formation of a pseudogap. To study the consequence of these two competing effects, we determine the spin-transport relaxation rate and the spin-diffusion constant using both a Boltzmann approach and a diagrammatic approach. The former ignores pseudogap physics and finite lifetime effects. In the latter, we incorporate the full pseudogap physics and lifetime effects, but we ignore vertex corrections, so that we effectively calculate single-particle relaxation rates instead of transport relaxation rates. We find that there is qualitative agreement between these two approaches although the results for the transport coefficients differ quantitatively., Comment: 9 pages, 10 figures

Recent advances in the creation and modulation of graphene-like systems are introducing a science of "designer Dirac materials". In its original definition, artificial graphene is a man-made nanostructure that consists of identical potential wells (quantum dots) arranged in a adjustable honeycomb lattice in the two-dimensional electron gas. As our ability to control the quality of artificial graphene samples improves, so grows the need for an accurate theory of its electronic properties, including the effects of electron-electron interactions. Here we determine those effects on the band structure and on the emergence of Dirac points.

Spin-orbit interaction is usefully classified as extrinsic or intrinsic depending on its origin: the potential due to random impurities (extrinsic), or the crystalline potential associated with the band or device structure (intrinsic). In this paper we will show how by using a SU(2) formulation the two sources of spin-orbit interaction may be described in an elegant and unified way. As a result we obtain a simple description of the interplay of the two types of spin-orbit interaction, and a physically transparent explanation of the vanishing of the d.c. spin Hall conductivity in a Rashba two-dimensional electron gas when spin relaxation is neglected, and its reinstatement when spin relaxation is allowed. Furthermore, we obtain an explicit formula for the transverse spin polarization created by an electric current, which generalizes the standard formula obtained by Edelstein and Aronov and Lyanda-Geller by including extrinsic spin-orbit interaction and spin relaxation., Comment: 9 pages

The non-equilibrium spin dynamics of a one-dimensional system of repulsively interacting fermions is studied by means of density-matrix renormalization-group simulations. We focus on the short-time decay of the oscillation amplitudes of the centers of mass of spin-up and spin-down fermions. Due to many-body effects, the decay is found to evolve from quadratic to linear in time, and eventually back to quadratic as the strength of the interaction increases. The characteristic rate of the decay increases linearly with the strength of repulsion in the weak-coupling regime, while it is inversely proportional to it in the strong-coupling regime. Our predictions can be tested in experiments on tunable ultra-cold few-fermion systems in one-dimensional traps., Comment: 10 pages, 12 figures, 4 appendices, final version accepted in PRL; main text and supplementary material merged in a single PRB style document

We calculate the optical spectra of silicon and germanium in the adiabatic time-dependent density functional formalism, making use of kinetic energy density-dependent (meta-GGA) exchange-correlation functionals. We find excellent agreement between theory and experiment. The success of the theory on this notoriously difficult problem is traced to the fact that the exchange-correlation kernel of meta-GGA supports a singularity of the form $\alpha/q^2$ (where $q$ is the wave-vector and $\alpha$ is a constant), whereas previously employed approximations (e.g. local density and generalized gradient approximations) do not. Thus, the use of the adiabatic meta-GGA opens a new path for handling the extreme non-locality of the time-dependent exchange-correlation potential in solid-state systems., Comment: 4 pages, 2 figures

We extend a recent formulation of quantum continuum mechanics [J. Tao et. al, Phys. Rev. Lett. {\bf 103}, 086401 (2009)] to many-body systems subjected to a magnetic field. To accomplish this, we propose a modified Lagrangian approach, in which motion of infinitesimal volume elements of the system is referred to the "quantum convective motion" that the magnetic field produces already in the ground-state of the system. In the linear approximation, this approach results in a redefinition of the elastic displacement field $\uv$, such that the particle current $\jv$ contains both an electric displacement and a magnetization contribution: $\jv=\jv_0+n_0\partial_t \uv+\nabla \times (\jv_0\times \uv)$, where $n_0$ and $\jv_0$ are the particle density and the current density of the ground-state and $\partial_t$ is the partial derivative with respect to time. In terms of this displacement, we formulate an "elastic approximation" analogous to the one proposed in the absence of magnetic field. The resulting equation of motion for $\uv$ is expressed in terms of ground-state properties -- the one-particle density matrix and the two-particle pair correlation function -- and in this form it neatly generalizes the equation obtained for vanishing magnetic field., Comment: 13 pages, revised version accepted to PRB

The onset of exciton condensation in a topological insulator thin film was recently predicted. We calculate the critical temperature for this transition, taking into account screening effects. Furthermore, we show that the proximity to this transition can be probed by measuring the Coulomb drag resistivity between the surfaces of the thin film as a function of temperature. This resistivity shows an upturn upon approaching the exciton-condensed state., Comment: 4 pages, 3 figures

Electrons in artificial lattices enable explorations of the impact of repulsive Coulomb interactions in a tunable system. We have trapped two-dimensional electrons belonging to a gallium arsenide quantum well in a nanofabricated lattice with honeycomb geometry. We probe the excitation spectrum in a magnetic field identifying novel collective modes that emerge from the Coulomb interaction in the artificial lattice as predicted by the Mott-Hubbard model. These observations allow us to determine the Hubbard gap and suggest the existence of a novel Coulomb-driven ground state. This approach offers new venues for the study of quantum phenomena in a controllable solid-state system., Comment: 10 pages, 10 figures

The presence of low-symmetry impurities or defect complexes in the zinc-blende direct-gap semiconductors (e.g. interstitials, DX-centers) results in a novel spin-orbit term in the effective Hamiltonian for the conduction band. The new extrinsic spin-orbit interaction is proportional to the matrix element of the defect potential between the conduction and the valence bands. Because this interaction arises already in the first order of the expansion of the effective Hamiltonian in powers of Uext/Eg << 1 (where Uext is the pseudopotential of an interstitial atom, and Eg is the band gap), its contribution to the spin relaxation rate may exceed that of the previously studied extrinsic contributions, even for moderate concentrations of impurities., Comment: extended version, 5+ pages

Quasi-two-dimensional superconductors with tunable spin-orbit coupling are very interesting systems with properties that are also potentially useful for applications. In this manuscript we demonstrate that these systems exhibit undamped collective spin oscillations that can be excited by the application of a supercurrent. We propose to use these collective excitations to realize persistent spin oscillators operating in the frequency range of 10 GHz - 1 THz., Comment: 10 pages, 4 figures, published version

Spin waves in insulating magnets are ideal carriers for spin currents with low energy dissipation. An electric field can modify the dispersion of spin waves, by directly affecting, via spin-orbit coupling, the electrons that mediate the interaction between magnetic ions. Our microscopic calculations based on the super-exchange model indicate that this effect of the electric field is sufficiently large to be used to effectively control spin currents. We apply these findings to the design of a spin-wave interferometric device, which acts as a logic inverter and can be used as a building block for room-temperature, low-dissipation logic circuits., Comment: 4 pages, 3 figures, added the LL equation and the discussion on spin-wave-induced electric field, accepted by PRL

Electrons in double-layer semiconductor heterostructures experience a special type of spin-orbit interaction which arises in each layer from the perpendicular component of the Coulomb electric field created by electron density fluctuations in the other layer. We show that this interaction, acting in combination with the usual spin-orbit interaction, can generate a spin current in one layer when a charge current is driven in the other. This effect is symmetry-wise distinct from the spin Hall drag. The spin current is not, in general, perpendicular to the drive current., Comment: 4 pages, 2 figures

We demonstrate that the plasmon frequency and Drude weight of the electron liquid in a doped graphene sheet are strongly renormalized by electron-electron interactions even in the long-wavelength limit. This effect is not captured by the Random Phase Approximation (RPA), commonly used to describe electron fluids and is due to coupling between the center of mass motion and the pseudospin degree of freedom of the graphene's massless Dirac fermions. Making use of diagrammatic perturbation theory to first order in the electron-electron interaction, we show that this coupling enhances both the plasmon frequency and the Drude weight relative to the RPA value. We also show that interactions are responsible for a significant enhancement of the optical conductivity at frequencies just above the absorption threshold. Our predictions can be checked by far-infrared spectroscopy or inelastic light scattering., Comment: 15 pages, 8 figures, submitted

We calculate the spin-drag relaxation rate for a two-component ultracold atomic Fermi gas with positive scattering length between the two spin components. In one dimension we find that it vanishes linearly with temperature. In three dimensions the spin-drag relaxation rate vanishes quadratically with temperature for sufficiently weak interactions. This quadratic temperature dependence is present, up to logarithmic corrections, in the two-dimensional case as well. For stronger interaction the system exhibits a Stoner ferromagnetic phase transition in two and three dimensions. We show that the spin-drag relaxation rate is enhanced by spin fluctuations as the temperature approaches the critical temperature of this transition from above., Comment: Submitted to New Journal of Physics Focus Issue "Strongly Correlated Quantum Fluids: From Ultracold Quantum Gases to QCD Plasmas"

We present an analytical study of the D'yakonov-Perel' spin relaxation time for degenerate electrons in a photo-excited electron-hole liquid in intrinsic semiconductors exhibiting a spin-split band structure. The D'yakonov-Perel' spin relaxation of electrons in these materials is controlled by electron-hole scattering, with small corrections from electron-electron scattering and virtually none from electron-impurity scattering. We derive simple expressions (one-dimensional and two-dimensional integrals respectively) for the effective electron-hole and electron-electron scattering rates which enter the spin relaxation time calculation. The electron-hole scattering rate is found to be comparable to the scattering rates from impurities in the electron liquid - a common model for n-type doped semiconductors. As the density of electron-hole pairs decreases (within the degenerate regime), a strong enhancement of the scattering rates and a corresponding slowing down of spin relaxation is predicted due to exchange and correlation effects in the electron-hole liquid. In the opposite limit of high density, the original D'yakonov-Perel' model fails due to decreasing scattering rates and is eventually superseded by free precession of individual quasiparticle spins., Comment: 16 pages, 5 figures

Spin-orbit-coupled two-dimensional electron gases (2DEGs) are a textbook example of helical Fermi liquids, i.e. quantum liquids in which spin (or pseudospin) and momentum degrees-of-freedom at the Fermi surface have a well-defined correlation. Here we study the long-wavelength plasmon dispersion and the Drude weight of archetypical spin-orbit-coupled 2DEGs. We first show that these measurable quantities are sensitive to electron-electron interactions due to broken Galileian invariance and then discuss in detail why the popular random phase approximation is not capable of describing the collective dynamics of these systems even at very long wavelengths. This work is focussed on presenting approximate microscopic calculations of these quantities based on the minimal theoretical scheme that captures the basic physics correctly, i.e. the time-dependent Hartree-Fock approximation. We find that interactions enhance the "plasmon mass" and suppress the Drude weight. Our findings can be tested by inelastic light scattering, electron energy loss, and far-infrared optical-absorption measurements., Comment: 18 pages, 11 figures, submitted

The orbital magnetic susceptibility (OMS) of a gas of noninteracting massless Dirac fermions is zero when the Fermi energy is away from the Dirac point. Making use of diagrammatic perturbation theory, we calculate exactly the OMS of massless Dirac fermions to first order in the Coulomb interaction demonstrating that it is finite and positive. Doped graphene sheets are thus unique systems in which the OMS is completely controlled by many-body effects., Comment: 4 pages, 2 figures, submitted

We express the high-frequency (anti-adiabatic) limit of the exchange-correlation kernels of an inhomogeneous electron gas in terms of the following equilibrium properties: the ground-state density, the kinetic stress tensor, the pair-correlation function, and the ground-state exchange-correlation potential. Of these quantities, the first three are amenable to exact evaluation by Quantum Monte Carlo methods, while the last can be obtained from the inversion of the Kohn-Sham equation for the ground-state orbitals. The exact scalar kernel, in this limit, is found to be of very long range in space, at variance with the kernel that is used in the standard local density approximation. The anti-adiabatic xc kernels should be useful in calculations of excitation energies by time-dependent DFT in atoms, molecules, and solids, and provides a solid basis for interpolation between the low- and high-frequency limits of the xc kernels., Comment: 9 pages, 3 figures, to be submitted to PRB