Search

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 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.

We construct a generalized-gradient approximation for the exchange-energy density of finite two-dimensional systems. Guided by non-empirical principles, we include the proper small-gradient limit and the proper tail for the exchange-hole potential. The observed performance is superior to that of the two-dimensional local-density approximation, which underlines the usefulness of the approach in practical applications.

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

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.

We derive a self-consistent local variant of the Thomas-Fermi approximation for (quasi-)two-dimensional (2D) systems by localizing the Hartree term. The scheme results in an explicit orbital-free representation of the electron density and energy in terms of the external potential, the number of electrons, and the chemical potential determined upon normalization. We test the method over a variety 2D nanostructures by comparing to the Kohn-Sham 2D-LDA calculations up to 600 electrons. Accurate results are obtained in view of the negligible computational cost. We also assess a local upper bound for the Hartree energy.

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

We study the graphite intercalated compound CaC$_6$ by means of Eliashberg theory, focusing on the anisotropy properties. An analysis of the electron-phonon coupling is performed, and we define a minimal 6-band anisotropy structure. Comparing with Superconducting Density Functional Theory (SCDFT) the condition under which Eliashberg theory is able to reproduce the SCDFT gap structure is determined, and we discuss the role of Coulomb interactions. The Engelsberg-Schrieffer polaron structure is computed by solving the Eliashberg equation on the Matsubara axis and analytically continuing it to the full complex plane. This reveals the polaronic quasiparticle bands anisotropic features as well as the interplay with superconductivity.

In modeling low-dimensional electronic nanostructures, the evaluation of the electron-electron interaction is a challenging task. Here we present an accurate and practical density-functional approach to the two-dimensional many-electron problem. In particular, we show that spin-density functionals in the class of meta-generalized-gradient approximations can be greatly simplified by reducing the explicit dependence on the Kohn-Sham orbitals to the dependence on the electron spin density and its spatial derivatives. Tests on various quantum-dot systems show that the overall accuracy is well preserved, if not even improved, by the modifications.

The exchange-correlation potentials stemming from the local-density approximation and several generalized-gradient approximations are known to have incorrect asymptotic decay. This failure is independent of the dimensionality, but so far the problem has been corrected -- within the mentioned approximations -- only in three dimensions. Here we provide a cured exchange-correlation potential in two dimensions, where the applications have a continuously increasing range in, e.g., semiconductor physics. The given potential is a generalized-gradient approximation, which is as easy to apply as the local-density approximation. We demonstrate that the corrected potential agrees very well with the analytic result of a two-electron quantum dot in the asymptotic regime, and yields plausible exchange-correlation potentials for larger two-dimensional systems.

Finite temperature density functional theory provides, in principle, an exact description of the thermodynamical equilibrium of many-electron systems. In practical applications, however, the functionals must be approximated. Efficient and physically meaningful approximations can be developed if relevant properties of the exact functionals are known and taken into consideration as constraints. In this work, derivations of exact properties and scaling relations for the main quantities of finite temperature density functional theory are presented. In particular, a coordinate scaling transformation at finite temperature is introduced and its consequences are elucidated.

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

Accurate treatment of the electronic correlation in inhomogeneous electronic systems, combined with the ability to capture the correlation energy of the homogeneous electron gas, allows to reach high predictive power in the application of density-functional theory. For two-dimensional systems we can achieve this goal by generalizing our previous approximation [Phys. Rev. B 79, 085316 (2009)] to a parameter-free form, which reproduces the correlation energy of the homogeneous gas while preserving the ability to deal with inhomogeneous systems. The resulting functional is shown to be very accurate for finite systems with an arbitrary number of electrons with respect to numerically exact reference data.

We present a practical and accurate density functional for the exchange-correlation energy of electrons in two dimensions. The exchange part is based on a recent two-dimensional generalized-gradient approximation derived by considering the limits of small and large density gradients. The fully local correlation part is constructed following the Colle-Salvetti scheme and a Gaussian approximation for the pair density. The combination of these expressions is shown to provide an efficient density functional to calculate the total energies of two-dimensional electron systems such as semiconductor quantum dots. Excellent performance of the functional with respect to numerically exact reference data for quantum dots is demonstrated., Comment: to appear in Int. J. Quantum Chem. (2010)

The Becke-Johnson exchange potential [J. Chem. Phys. 124, 221101 (2006)] has been successfully used in electronic structure calculations within density-functional theory. However, in its original form the potential may dramatically fail in systems with non-Coulombic external potentials, or in the presence of external magnetic or electric fields. Here, we provide a system-independent correction to the Becke-Johnson approximation by (i) enforcing its gauge invariance and (ii) making it exact for any single-electron system. The resulting approximation is then better designed to deal with current-carrying states, and recovers the correct asymptotic behavior for systems with any number of electrons. Tests of the resulting corrected exchange potential show very good results for a Hydrogen chain in an electric field and for a four-electron harmonium in a magnetic field.

We consider density functionals for exchange and correlation energies in two-dimensional systems. The functionals are constructed by making use of exact constraints for the angular averages of the corresponding exchange and correlation holes, respectively, and assuming proportionality between their characteristic sizes. The electron current and spin are explicitly taken into account, so that the resulting functionals are suitable to deal with systems exhibiting orbital currents and/or spin polarization. Our numerical results show that in finite systems the proposed functionals outperform the standard two-dimensional local spin-density approximation, still performing well also in the important limit of the homogeneous two-dimensional electron gas.

We study the properties of the lower bound on the exchange-correlation energy in two dimensions. First we review the derivation of the bound and show how it can be written in a simple density-functional form. This form allows an explicit determination of the prefactor of the bound and testing its tightness. Next we focus on finite two-dimensional systems and examine how their distance from the bound depends on the system geometry. The results for the high-density limit suggest that a finite system that comes as close as possible to the ultimate bound on the exchange-correlation energy has circular geometry and a weak confining potential with a negative curvature.

We extend the Becke-Johnson approximation [J. Chem. Phys. 124, 221101 (2006)] of the exchange potential to two dimensions. We prove and demonstrate that a direct extension of the underlying formalism may lead to divergent behavior of the potential. We derive a cure to the approach by enforcing the gauge invariance and correct asymptotic behavior of the exchange potential. The procedure leads to an approximation which is shown, in various quasi-two-dimensional test systems, to be very accurate in comparison with the exact exchange potential, and thus a considerable improvement over the commonly applied local-density approximation., Comment: submitted to Phys. Rev. B on July 9th, 2009

We derive a non-empirical, orbital-free density functional for the total energy of interacting electrons in two dimensions. The functional consists of a local formula for the interaction energy, where we follow the lines introduced by Parr for three-dimensional systems [R. G. Parr, J. Phys. Chem. 92, 3060 (1988)], and the Thomas-Fermi approximation for the kinetic energy. The freedom from orbitals and from the Hartree integral makes the proposed approximation numerically highly efficient. The total energies obtained for confined two-dimensional systems are in a good agreement with the standard local-density approximation within density-functional theory, and considerably more accurate than the Thomas-Fermi approximation.

Electronic structure calculations are routinely carried out within the framework of density-functional theory, often with great success. For electrons in reduced dimensions, however, there is still a need for better approximations to the exchange-correlation energy functional. Furthermore, the need for properly describing current-carrying states represents an additional challenge for the development of approximate functionals. In order to make progress along these directions, we show that simple and efficient expressions for the exchange energy can be obtained by considering the short-range behavior of the one-body spin-density matrix. Applications to several two-dimensional systems confirm the excellent performance of the derived approximations, and verify the gauge-invariance requirement to be of great importance for dealing with current-carrying states.

Superconductivity in Pb, H under extreme pressure and CaBeSi, in the framework of the density functional theory for superconductors, is discussed. A detailed analysis on how the electron-phonon and electron-electron interactions combine together to determine the superconducting gap and critical temperature of these systems is presented. Pb, H under pressure and CaBeSi are multigap superconductors. We will address the question under which conditions does a system exhibits this phenomenon. The presented results contribute to the understanding of multiband and anisotropic superconductivity, which has received a lot of attention since the discovery of MgB$_2$, and show how it is possible to describe the superconducting properties of real materials on a fully {\it ab-initio} basis.

Bounds on the exchange-correlation energy of many-electron systems are derived and tested. By using universal scaling properties of the electron-electron interaction, we obtain the exponent of the bounds in three, two, one, and quasi-one dimensions. From the properties of the electron gas in the dilute regime, the tightest estimate to date is given for the numerical prefactor of the bound, which is crucial in practical applications. Numerical tests on various low-dimensional systems are in line with the bounds obtained, and give evidence of an interesting dimensional crossover between two and one dimensions.

Superconductivity in intercalated graphite CaC6 and H under extreme pressure, in the framework of superconducting density functional theory, is discussed. A detailed analysis on how the electron-phonon and electron-electron interactions combine together to determine the superconducting gap and critical temperature of these systems is presented. In particular, we discuss the effect on the calculated Tc of the anisotropy of the electron-phonon interaction and of the different approximations for screening the Coulomb repulsion. These results contribute to the understanding of multigap and anisotropic superconductivity, which has received a lot of attention since the discovery of MgB2, and show how it is possible to describe the superconducting properties of real materials on a fully ab-initio basis., Comment: 13 pages, 5 figures

We derive a generalized gradient approximation to the exchange energy to be used in density functional theory calculations of two-dimensional systems. This class of approximations has a long and successful history, but it has not yet been fully investigated for electrons in two dimensions. We follow the approach originally proposed by Becke for three-dimensional systems [Int. J. Quantum Chem. 23, 1915 (1983), J. Chem. Phys. 85, 7184 (1986)]. The resulting functional depends on two parameters that are adjusted to a test set of parabolically confined quantum dots. Our exchange functional is then tested on a variety of systems with promising results, reducing the error in the exchange energy by a factor of 4 with respect to the simple local density approximation.

The capability of density-functional theory to deal with the ground-state of strongly correlated low-dimensional systems, such as semiconductor quantum dots, depends on the accuracy of functionals developed for the exchange and correlation energies. Here we extend a successful approximation for the correlation energy of the three dimensional inhomogeneous electron gas, originally introduced by Becke [J. Chem. Phys. {\bf 88}, 1053 (1988)], to the two-dimensional case. The approach aims to non-empirical modeling of the correlation-hole functions satisfying a set of exact properties. Furthermore, the electron current and spin are explicitly taken into account. As a result, good performance is obtained in comparison with numerically exact data for quantum dots with varying external magnetic field, and for the homogeneous two-dimensional electron gas, respectively.

Quantum rings can be characterized by a specific radius and ring width. For this rich class of physical systems, an accurate approximation for the exchange-hole potential and thus for the exchange energy is derived from first principles. Excellent agreement with the exact-exchange results is obtained regardless of the ring parameters, total spin, current, or the external magnetic field. The description can be applied as a density functional outperforming the commonly used local-spin-density approximation, which is here explicitly shown to break down in the quasi-one-dimensional limit. The dimensional crossover, which is of extraordinary importance in low-dimensional systems, is fully captured by our functional.

We derive a local approximation for the correlation energy in two-dimensional electronic systems. In the derivation we follow the scheme originally developed by Colle and Salvetti for three dimensions, and consider a Gaussian approximation for the pair density. Then, we introduce an ad-hoc modification which better accounts for both the long-range correlation, and the kinetic-energy contribution to the correlation energy. The resulting functional is local, and depends parametrically on the number of electrons in the system. We apply this functional to the homogeneous electron gas and to a set of two-dimensional quantum dots covering a wide range of electron densities and thus various amounts of correlation. In all test cases we find an excellent agreement between our results and the exact correlation energies. Our correlation functional has a form that is simple and straightforward to implement, but broadly outperforms the commonly used local-density approximation.

Implicit and explicit density functionals for the exchange energy in finite two-dimensional systems are developed following the approach of Becke and Roussel [Phys. Rev. A 39, 3761 (1989)]. Excellent agreement for the exchange-hole potentials and exchange energies is found when compared with the exact-exchange reference data for the two-dimensional uniform electron gas and few-electron quantum dots, respectively. Thereby, this work significantly improves the availability of approximate density functionals for dealing with electrons in quasi-two-dimensional structures, which have various applications in semiconductor nanotechnology., Comment: 5 pages, 3 figures

Popular density functionals for the exchange-correlation energy typically fail to reproduce the degeneracy of different ground states of open-shell atoms. As a remedy, functionals which explicitly depend on the current density have been suggested. We present an analysis of this problem by investigating functionals that explicitly depend on the Kohn-Sham orbitals. Going beyond the exact-exchange approximation by adding correlation in the form of the Colle-Salvetti functional we show how current-dependent terms enter the Colle-Salvetti expression and their relevance is evaluated. A very good description of the degeneracy of ground-states for atoms of the first and second row of the periodic table is obtained.

The relative merits of current-spin-density- and spin-density-functional theory are investigated for solids treated within the exact-exchange-only approximation. Spin-orbit splittings and orbital magnetic moments are determined at zero external magnetic field. We find that for magnetic (Fe, Co and Ni) and non-magnetic (Si and Ge) solids, the exact-exchange current-spin-density functional approach does not significantly improve the accuracy of the corresponding spin-density functional results., Comment: 4 pages, 1 fig, 2 tables

Current-spin density functional theory (CSDFT) provides a framework to describe interacting many-electron systems in a magnetic field which couples to both spin- and orbital-degrees of freedom. Unlike in usual (spin-) density functional theory, approximations to the exchange-correlation energy based on the model of the uniform electron gas face problems in practical applications. In this work, explicitly orbital-dependent functionals are used and a generalization of the Optimized Effective Potential (OEP) method to the CSDFT framework is presented. A simplifying approximation to the resulting integral equations for the exchange-correlation potentials is suggested. A detailed analysis of these equations is carried out for the case of open-shell atoms and numerical results are given using the exact-exchange energy functional. For zero external magnetic field, a small systematic lowering of the total energy for current-carrying states is observed due to the inclusion of the current in the Kohn-Sham scheme. For states without current, CSDFT results coincide with those of spin density functional theory., Comment: 11 pages, 3 figures

The description of interacting many-electron systems in external magnetic fields is considered in the framework of the optimized effective potential method extended to current-spin-density functional theory. As a case study, a two-dimensional quantum dot in external magnetic fields is investigated. Excellent agreement with quantum Monte Carlo results is obtained when self-interaction corrected correlation energies from the standard local spin-density approximation are added to exact-exchange results. Full self-consistency within the complete current-spin-density-functional framework is found to be of minor importance., Comment: 5 pages, 2 figures, submitted to PRB

The problem of degenerate ground states of open-shell atoms is investigated in spin-restricted and unrestricted density functional theory using the exact exchange energy functional. For the spin-unrestricted case, spurious energy splittings of the order of 2 to 3 kcal/mol are found for atoms of the second and third period which is larger than the splittings obtained from recently proposed approximate exchange functionals depending explicitly on the current density. In remarkable contrast, for spin-restricted calculations the degeneracy of different atomic ground states is recovered to within less than 0.6 kcal/mol., Comment: 5 pages

A description of non-collinear magnetism in the framework of spin-density functional theory is presented for the exact exchange energy functional which depends explicitly on two-component spinor orbitals. The equations for the effective Kohn-Sham scalar potential and magnetic field are derived within the optimized effective potential (OEP) framework. With the example of a magnetically frustrated Cr monolayer it is shown that the resulting magnetization density exhibits much more non-collinear structure than standard calculations. Furthermore, a time-dependent generalization of the non-collinear OEP method is well suited for an ab-initio description of spin dynamics. We also show that the magnetic moments of solids Fe, Co and Ni are well reproduced., Comment: 4 pages 2 fig accepted in Phys. Rev. Lett. (2007)