Source author record

E. K. U. Gross

E. K. U. Gross appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.chem-ph cond-mat.mtrl-sci cond-mat.other cond-mat.str-el cond-mat.supr-con cond-mat.mes-hall quant-ph physics.comp-ph physics.atom-ph physics.atm-clus cond-mat.stat-mech math-ph math.MP physics.optics

Catalog footprint

What is connected

89works

14topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2025arXiv

Magnetism-Enhanced Strong Electron-Phonon Coupling in Infinite-Layer Nickelate

Intriguing analogies between the nickelates and the cuprates provide a promising avenue for unraveling the microscopic mechanisms underlying high-$T_c$ superconductivity. While electron correlation effects in the nickelates have been extensively studied, the role of electron-phonon coupling (EPC) remains highly controversial. Here, by taking pristine LaNiO$_2$ as an exemplar nickelate, we present an in-depth study of EPC for both the non-magnetic (NM) and the $C$-type antiferromagnetic ($C$-AFM) phase using advanced density functional theory methods without invoking $U$ or other free parameters. The weak EPC strength $λ$ in the NM phase is found to be greatly enhanced ($\sim$4$\times$) due to the presence of magnetism in the $C$-AFM phase. This enhancement arises from strong interactions between the flat bands associated with the Ni-3$d_{z^2}$ orbitals and the low-frequency phonon modes driven by the vibrations of Ni and La atoms. The resulting phonon softening is shown to yield a distinctive kink in the electronic structure around 15 meV, which would provide an experimentally testable signature of our predictions. Our study highlights the critical role of local magnetic moments and interply EPC in the nickelate.

preprint2022arXiv

Energy, Momentum and Angular Momentum Transfer Between Electrons and Nuclei

The recently developed exact factorization approach condenses all electronic effects on the nuclear subsystem into scalar and vector potentials that appear in an effective time dependent Schrödinger equation. Starting from this equation, we derive subsystem Ehrenfest identities characterizing the energy, momentum and angular momentum transfer between electrons and nuclei. An effective electromagnetic force operator induced by the electromagnetic field corresponding to the effective scalar and vector potentials appears in all three identities. The effective magnetic field has two components that can be identified with the Berry curvature calculated with (a) different cartesian coordinates of the same nucleus and (b) arbitrary cartesian coordinates of two different nuclei. (a) has a classical interpretation as the induced magnetic field felt by the nucleus, while (b) has no classical analog. Subsystem Ehrenfest identities are ideally suited for quantifying energy transfer in electron-phonon systems. With two explicit examples we demonstrate the usefulness of the new identities.

preprint2021arXiv

From Kohn-Sham to many-electron energies via step structures in the exchange-correlation potential

Accurately describing excited states within Kohn-Sham (KS) density functional theory (DFT), particularly those which induce ionization and charge transfer, remains a great challenge. Common exchange-correlation (xc) approximations are unreliable for excited states owing, in part, to the absence of a derivative discontinuity in the xc energy ($Δ$), which relates a many-electron energy difference to the corresponding KS energy difference. We demonstrate, analytically and numerically, how the relationship between KS and many-electron energies leads to the step structures observed in the exact xc potential, in four scenarios: electron addition, molecular dissociation, excitation of a finite system, and charge transfer. We further show that steps in the potential can be obtained also with common xc approximations, as simple as the LDA, when addressed from the ensemble perspective. The article therefore highlights how capturing the relationship between KS and many-electron energies with advanced xc approximations is crucial for accurately calculating excitations, as well as the ground-state density and energy of systems which consist of distinct subsystems.

preprint2021arXiv

Geometric energy transfer in two-component systems

Factoring a wave function into marginal and conditional factors partitions the subsystem kinetic energy into two terms. The first depends solely on the marginal wave function, through its gauge-covariant derivative, while the second depends on the quantum metric of the conditional wave function over the manifold of marginal variables. We derive an identity for the rate of change of the second term.

preprint2020arXiv

Ab initio theory of plasmonic superconductivity within the Eliashberg and density-functional formalisms

We extend the two leading methods for the \emph{ab initio} computational descrip tion of phonon-mediated superconductors, namely Eliashberg theory and density fu nctional theory for superconductors (SCDFT), to include plasmonic effects. Furth ermore, we introduce a hybrid formalism in which the Eliashberg approximation fo r the electron-phonon coupling is combined with the SCDFT treatment of the dynam ically screened Coulomb interaction. The methods have been tested on a set of we ll-known conventional superconductors by studying how the plasmon contribution a ffects the phononic mechanism in determining the critical temperature (\tc). Our simulations show that plasmonic SCDFT leads to a good agreement between predict ed and measured \tc's, whereas Eliashberg theory considerably overestimates the plasmon-mediated pairing and, therefore, \tc. The hybrid approach, on the other hand, gives results close to SCDFT and overall in excellent agreement with exper iments.

preprint2020arXiv

Extending solid-state calculations to ultra long-range length scales

We present a method which enables solid-state density functional theory calculations to be applied to systems of almost unlimited size. Computations of physical effects up to the micron length scale but which nevertheless depend on the microscopic details of the electronic structure, are made possible. Our approach is based on a generalization of the Bloch state which involves an additional sum over a finer grid in reciprocal space around each ${\bf k}$-point. We show that this allows for modulations in the density and magnetization of arbitrary length on top of a lattice-periodic solution. Based on this, we derive a set of ultra long-range Kohn-Sham equations. We demonstrate our method with a sample calculation of bulk LiF subjected to an arbitrary external potential containing nearly 3500 atoms. We also confirm the accuracy of the method by comparing the spin density wave state of bcc Cr against a direct supercell calculation starting from a random magnetization density. Furthermore, the spin spiral state of $γ$-Fe is correctly reproduced and the screening by the density of a saw-tooth potential over 20 unit cells of silicon is verified.

preprint2020arXiv

Many-body Green's function theory of electrons and nuclei beyond the Born-Oppenheimer approximation

The method of many-body Green's functions is developed for arbitrary systems of electrons and nuclei starting from the full (beyond Born-Oppenheimer) Hamiltonian of Coulomb interactions and kinetic energies. The theory presented here resolves the problems arising from the translational and rotational invariance of this Hamiltonian that afflict the existing many-body Green's function theories. We derive a coupled set of exact equations for the electronic and nuclear Green's functions and provide a systematic way to approximately compute the properties of arbitrary many-body systems of electrons and nuclei beyond the Born-Oppenheimer approximation. The case of crystalline solids is discussed in detail.

preprint2017arXiv

How is the derivative discontinuity related to steps in the exact Kohn-Sham potential?

The reliability of density-functional calculations hinges on accurately approximating the unknown exchange-correlation (xc) potential. Common (semi-)local xc approximations lack the jump experienced by the exact xc potential as the number of electrons infinitesimally surpasses an integer, and the spatial steps that form in the potential as a result of the change in the decay rate of the density. These features are important for an accurate prediction of the fundamental gap and the distribution of charge in complex systems. Although well-known concepts, the exact relationship between them remained unclear. In this Letter, we establish the common fundamental origin of these two features of the exact xc potential via an analytical derivation. We support our result with an exact numerical solution of the many-electron Schroedinger equation for a single atom and a diatomic molecule in one dimension. Furthermore, we propose a way to extract the fundamental gap from the step structures in the potential.

preprint2016arXiv

A single-electron approach for many-electron dynamics in high-order harmonic generation

We present a novel ab-initio single-electron approach to correlated electron dynamics in strong laser fields. By writing the electronic wavefunction as a product of a marginal one-electron wavefunction and a conditional wavefunction, we show that the exact harmonic spectrum can be obtained from a single-electron Schrödinger equation. To obtain the one-electron potential in practice, we propose an adiabatic approximation, i.e. a potential is generated that depends only on the position of one electron. This potential, together with the laser interaction, is then used to obtain the dynamics of the system. For a model Helium atom in a laser field, we show that by using our approach, the high-order harmonic generation spectrum can be obtained to a good approximation.

preprint2016arXiv

An Exact Factorization Perspective on Quantum Interferences in Nonadiabatic Dynamics

Nonadiabatic quantum interferences emerge whenever nuclear wavefunctions in different electronic states meet and interact in a nonadiabatic region. In this work, we analyze how nonadiabatic quantum interferences translate in the context of the exact factorization of the molecular wavefunction. In particular, we focus our attention on the shape of the time-dependent potential energy surface - the exact surface on which the nuclear dynamics takes place. We use a one-dimensional exactly-solvable model to reproduce different conditions for quantum interferences, whose characteristic features already appear in one-dimension. The time-dependent potential energy surface develops complex features when strong interferences are present, in clear contrast to the observed behavior in simple nonadiabatic crossing cases. Nevertheless, independent classical trajectories propagated on the exact time-dependent potential energy surface reasonably conserve a distribution in configuration space that mimics the one of the exact nuclear probability density.

preprint2016arXiv

Exact factorization-based density functional theory of electrons and nuclei

The ground state energy of a system of electrons and nuclei is proven to be a variational functional of the conditional electronic density $n_R(\mathbf{r})$, the nuclear wavefunction $χ(R)$ and an induced vector potential $A_μ(R)$ and quantum geometric tensor $\mathcal{T}_{μν}(R)$ derived from the conditional electronic wavefunction $Φ_R(r)$ over nuclear configuration space, where $r=\mathbf{r}_1,\mathbf{r}_2,\ldots$ are electronic coordinates and $R=\mathbf{R}_1,\mathbf{R}_2,\ldots$ are nuclear coordinates. The ground state $(n_R,χ,A_μ,\mathcal{T}_{μν})$ can be calculated by solving self-consistently (i) conditional Kohn-Sham equations containing an effective potential $v_{\rm s}(\mathbf{r})$ that depends parametrically on $R$, (ii) the Schrödinger equation for $χ(R)$ and (iii) Euler-Lagrange equations that determine $\mathcal{T}_{μν}$. The theory is applied to the $E\otimes e$ Jahn-Teller model.

preprint2016arXiv

Molecular geometric phase from the exact electron-nuclear factorization

The Born-Oppenheimer electronic wavefunction $Φ_R^{BO}(r)$ picks up a topological phase factor $\pm 1$, a special case of Berry phase, when it is transported around a conical intersection of two adiabatic potential energy surfaces in $R$-space. We show that this topological quantity reverts to a geometric quantity $e^{iγ}$ if the geometric phase $γ= \oint \mathrm{Im} \langle Φ_R |\nabla_μ Φ_R\rangle \cdot d\mathbf{R}_μ$ is evaluated with the conditional electronic wavefunction $Φ_R(r)$ from the exact electron-nuclear factorization $Φ_R(r)χ(R)$ instead of the adiabatic function $Φ_R^{BO}(r)$. A model of a pseudorotating molecule, also applicable to dynamical Jahn-Teller ions in bulk crystals, provides the first examples of induced vector potentials and molecular geometric phase from the exact factorization. The induced vector potential gives a contribution to the circulating nuclear current which cannot be removed by a gauge transformation. The exact potential energy surface is calculated and found to contain a term depending on the Fubini-Study metric for the conditional electronic wavefunction.

preprint2016arXiv

Quantum-Classical Non-Adiabatic Dynamics: Coupled- vs. Independent-Trajectory Methods

Trajectory-based mixed quantum-classical approaches to coupled electron-nuclear dynamics suffer from well-studied problems such as the lack of (or incorrect account for) decoherence in the trajectory surface hopping method and the inability of reproducing the spatial splitting of a nuclear wave packet in Ehrenfest-like dynamics. In the context of electronic non-adiabatic processes, these problems can result in wrong predictions for quantum populations and in unphysical outcomes for the nuclear dynamics. In this paper we propose a solution to these issues by approximating the coupled electronic and nuclear equations within the framework of the exact factorization of the electron-nuclear wave function. We present a simple quantum-classical scheme based on coupled classical trajectories, and test it against the full quantum mechanical solution from wave packet dynamics for some model situations which represent particularly challenging problems for the above-mentioned traditional methods.

preprint2016arXiv

Ultrafast Local Magnetization and Demagnetization in Heusler Alloys

With the goal of pushing Spintronic devices towards faster and faster timescales, we demonstrate, using ab-intio time-dependent density functional theory simulations of bulk Heusler compounds subject to ultrashort intense laser pulses, that the local magnetic moment can increase or decrease in a few femtoseconds. This speed is due to the all optical nature of the process, which transfers spin moment from one sublattice to another. This transfer depends on easily tunable laser parameters. By comparing the spin dynamics of a variety of Heusler (or half-Heusler) compounds, we demonstrate that the density of states explains the observed moment transfer; most the physics of inter sublattice moment transfer is due to the flow of spin current which is governed by availability of states above the Fermi level.

preprint2015arXiv

A coupled-trajectory quantum-classical approach to decoherence in non-adiabatic processes

We present a novel quantum-classical approach to non-adiabatic dynamics, deduced from the coupled electronic and nuclear equations in the framework of the exact factorization of the electron-nuclear wave function. The method is based on the quasi-classical interpretation of the nuclear wave function, whose phase is related to the classical momentum and whose density is represented in terms of classical trajectories. In this approximation, electronic decoherence is naturally induced as effect of the coupling to the nuclei and correctly reproduces the expected quantum behaviour. Moreover, the splitting of the nuclear wave packet is captured as consequence of the correct approximation of the time-dependent potential of the theory. This new approach offers a clear improvement over Ehrenfest-like dynamics. The theoretical derivation presented in the Letter is supported by numerical results that are compared to quantum mechanical calculations.

preprint2015arXiv

Ab initio calculation of a Pb single layer on a Si substrate: two-dimensionality and superconductivity

We report on first principles calculations of superconductivity in a single layer of lead on a silicon substrate including a full treatment of phononic and RPA screened coulomb interactions within the parameter free framework of Density Functional Theory for superconductors. A thorough investigation shows that several approximations that are commonly valid in bulk systems fail in this constrained 2D geometry. The calculated critical temperature turns out to be much higher than the experimental value of 1.86K. We argue that the only plausible explanation for the experimental Tc suppression is the onset of fluctuations of the superconducting order parameter.

preprint2015arXiv

Ab Initio Theory of Superconductivity in a Magnetic Field I. : Spin Density Functional Theory For Superconductors and Eliashberg Equations

We present a first-principles approach to describe magnetic and superconducting systems and the phenomena of competition between these electronic effects. We develop a density functional theory: SpinSCDFT, by extending the Hohenberg-Kohn theorem and constructing the non-interacting Kohn- Sham system. An exchange-correlation functional for SpinSCDFT is derived from the Sham Schlüter connection between the SpinSCDFT Kohn-Sham and a self-energy in Eliashberg approximation. The reference Eliashberg equations for superconductors in the presence of magnetism are also derived and discussed.

preprint2015arXiv

Ab-Initio Theory of Superconductivity in a Magnetic Field II. : Numerical solution

We numerically investigate the Spin Density Functional theory for superconductors (SpinSCDFT) and the approximated exchange-correlation functional, derived and presented in the preceding paper I. As a test system we employ a free electron gas featuring an exchange-splitting, a phononic pairing field and a Coulomb repulsion. SpinSCDFT results are compared with Sarma, the Bardeen Cooper and Schrieffer theory and with an Eliashberg type of approach. We find that the spectrum of the superconducting Kohn-Sham SpinSCDFT system is not in agreement with the true quasi particle structure. Therefore, starting from the Dyson equation, we derive a scheme that allows to compute the many body excitations of the superconductor and represents the extension to superconductivity of the G0W0 method in band structure theory. This superconducting G0 W0 method vastly improves the predicted spectra.

preprint2015arXiv

Demonstration of Optimal Control of Laser Induced Spin-Orbit Mediated Ultrafast Demagnetization

Laser induced ultrafast demagnetization is the process whereby the magnetic moment of a ferromagnetic material is seen to drop significantly on a timescale of $10-100$s of femtoseconds due to the application of a strong laser pulse. If this phenomenon can be harnessed for future technology, it offers the possibility for devices operating at speeds several orders of magnitude faster than at present. A key component to successful transfer of such a process to technology is the controllability of the process, i.e. that it can be tuned in order to overcome the practical and physical limitations imposed on the system. In this paper, we demonstrate that the spin-orbit mediated form of ultrafast demagnetization recently investigated [arXiv:1406.6607] by ab-initio time-dependent density functional theory (TDDFT) can be controlled. To do so we use quantum optimal control theory (OCT) to couple our TDDDT simulations to the optimization machinery of OCT. We show that a laser pulse can be found which maximizes the loss of moment within a given time interval while subject to several practical and physical constraints. Furthermore we also include a constraint on the fluence of the laser pulses and find the optimal pulse that combines significant demagnetization with a desire for less powerful pulses. These calculations demonstrate optimal control is possible for spin-orbit mediated ultrafast demagnetization and lay the foundation for future optimizations/simulations which can incorporate even more constraints.

preprint2015arXiv

Doping induced metal-insulator phase transition in NiO

The insulator to metal phase transition in NiO is studied within the framework of reduced density matrix functional theory and density functional theory. We find that the spectral density obtained using reduced density matrix functional theory is in good agreement with experiments both undoped as well as doped NiO. We find that the physical description of the hole-doping induced phase transition \emph{qualitatively differs} depending on whether NiO is calculated within density functional theory or reduced density matrix functional. In the former case the underlying mechanism of the phase transition is identified to be a rigid shift of chemical potential, while in the latter case a redistribution of the spectral weight drives the transition. These latter results are found to be in good agreement with both experiments and previous many-body calculations.

preprint2015arXiv

High temperature superconductivity in sulfur and selenium hydrides at high pressure

Due to its low atomic mass hydrogen is the most promising element to search for high-temperature phononic superconductors. However, metallic phases of hydrogen are only expected at extreme pressures (400 GPa or higher). The measurement of a record superconducting critical temperature of 190 K in a hydrogen-sulfur compound at 200 GPa of pressure[1], shows that metallization of hydrogen can be reached at significantly lower pressure by inserting it in the matrix of other elements. In this work we re-investigate the phase diagram and the superconducting properties of the H-S system by means of minima hopping method for structure prediction and Density Functional theory for superconductors. We also show that Se-H has a similar phase diagram as its sulfur counterpart as well as high superconducting critical temperature. We predict SeH3 to exceed 120 K superconductivity at 100 GPa. We show that both SeH3 and SH3, due to the critical temperature and peculiar electronic structure, present rather unusual superconducting properties.

preprint2015arXiv

L1$_0$ stacked binaries as candidates for hard magnets: FePt, MnAl and MnGa

A novel strategy of stacking binary magnets to enhance the magneto crystalline anisotropy is explored. This strategy is used in the search for hard magnets by studying FePt/MnGa and FePt/MnAl stacks. The choice of these binaries is motivated by the fact that they already possess large magneto crystalline anisotropy. Several possible alternative structures for these materials are explored in order to reduce the amount of Pt owing to its high cost.

preprint2015arXiv

Laser-induced electron localization in H$_2^+$: Mixed quantum-classical dynamics based on the exact time-dependent potential energy surface

We study the exact nuclear time-dependent potential energy surface (TDPES) for laser-induced electron localization with a view to eventually developing a mixed quantum-classical dynamics method for strong-field processes. The TDPES is defined within the framework of the exact factorization [A. Abedi, N. T. Maitra, and E. K. U. Gross, Phys. Rev. Lett. 105, 123002 (2010)] and contains the exact effect of the couplings to the electronic subsystem and to any external fields within a scalar potential. We compare its features with those of the quasistatic potential energy surfaces (QSPES) often used to analyse strong-field processes. We show that the gauge-independent component of the TDPES has a mean-field-like character very close to the density-weighted average of the QSPESs. Oscillations in this component are smoothened out by the gauge-dependent component, and both components are needed to yield the correct force on the nuclei. Once the localization begins to set in, the gradient of the exact TDPES tracks one QSPES and then switches to the other, similar to the description provided by surface-hopping between QSPESs. We show that evolving an ensemble of classical nuclear trajectories on the exact TDPES accurately reproduces the exact dynamics. This study suggests that the mixed quantum-classical dynamics scheme based on evolving multiple classical nuclear trajectories on the exact TDPES will be a novel and useful method to simulate strong field processes.

preprint2015arXiv

Nuclear velocity perturbation theory for vibrational circular dichroism: An approach based on the exact factorization of the electron-nuclear wave function

The nuclear velocity perturbation theory (NVPT) for vibrational circular dichroism (VCD) is derived from the exact factorization of the electron-nuclear wave function. This new formalism offers an exact starting point to include correction terms to the Born-Oppenheimer (BO) form of the molecular wave function, similarly to the complete-adiabatic approximation. The corrections depend on a small parameter that, in a classical treatment of the nuclei, is identified as the nuclear velocity. Apart from proposing a rigorous basis for the NVPT, we show that the rotational strengths, related to the intensity of the VCD signal, contain a new contribution beyond-BO that can be evaluated with the NVPT and that only arises when the exact factorization approach is employed. Numerical results are presented for chiral and non-chiral systems to test the validity of the approach.

preprint2015arXiv

Reduced Density Matrix Functional Theory at Finite Temperature: Theoretical Foundations

We present an ab-initio approach for grand canonical ensembles in thermal equilibrium with local or nonlocal external potentials based on the one-reduced density matrix. We show that equilibrium properties of a grand canonical ensemble are determined uniquely by the eq-1RDM and establish a variational principle for the grand potential with respect to its one-reduced density matrix. We further prove the existence of a Kohn-Sham system capable of reproducing the one-reduced density matrix of an interacting system at finite temperature. Utilizing this Kohn-Sham system as an unperturbed system, we deduce a many-body approach to iteratively construct approximations to the correlation contribution of the grand potential.

preprint2015arXiv

Spectrum for non-magnetic Mott insulators from power functional within Reduced Density Matrix Functional Theory

A fully first principles theory capable of treating strongly correlated solids remains the outstanding challenge of modern day materials science. This is exemplified by the transition metal oxides, prototypical Mott insulators, that remain insulating even in the absence of long range magnetic order. Capturing the non-magnetic insulating state of these materials presents a difficult challenge for any modern electronic structure theory. In this paper we demonstrate that reduced density matrix functional theory, in conjunction with the power functional, can successfully treat the non-magnetic insulating state of the transition metal oxides NiO and MnO. We show that the electronic spectrum retains a gap even in the absence of spin order. We further discuss the detailed way in which RDMFT performs for Mott insulators and band insulators, finding that for the latter occupation number minimization alone is required, but for the former full minimization over both occupation numbers and natural orbitals is essential.

preprint2015arXiv

Superconductivity in metastable phases of phosphorus-hydride compounds under high pressure

Hydrogen-rich compounds have been extensively studied both theoretically and experimentally in the quest for novel high-temperature superconductors. Reports on sulfur-hydride attaining metallicity under pressure and exhibiting superconductivity at temperatures as high as 200 K have spurred an intense search for room-temperature superconductors in hydride materials. Recently, compressed phosphine was reported to metallize at pressures above 45 GPa, reaching a superconducting transition temperature (T$_{c}$) of 100 K at 200 GPa. However, neither the exact composition nor the crystal structure of the superconducting phase have been conclusively determined. In this work the phase diagram of PH$_n$ ($n=1,2,3,4,5,6$) was extensively explored by means of {\it ab initio} crystal structure predictions using the Minima Hopping Method (MHM). The results do not support the existence of thermodynamically stable PH$_n$ compounds, which exhibit a tendency for elemental decomposition at high pressure even when vibrational contributions to the free energies are taken into account. Although the lowest energy phases of PH$_{1,2,3}$ display T$_{c}$'s comparable to experiments, it remains uncertain if the measured values of T$_{c}$ can be fully attributed to a phase-pure compound of PH$_n$.

preprint2015arXiv

The exact forces on classical nuclei in non-adiabatic charge transfer

The decomposition of electronic and nuclear motion presented in~[A. Abedi, N. T. Maitra, and E. K. U. Gross, Phys. Rev. Lett. 105, 123002 (2010)] yields a time-dependent potential that drives the nuclear motion and fully accounts for the coupling to the electronic subsystem. Here we show that propagation of an ensemble of independent classical nuclear trajectories on this exact potential yields dynamics that are essentially indistinguishable from the exact quantum dynamics for a model non-adiabatic charge transfer problem. We point out the importance of step and bump features in the exact potential that are critical in obtaining the correct splitting of the quasiclassical nuclear wave packet in space after it passes through an avoided crossing between two Born-Oppenheimer surfaces, and analyze their structure. Lastly, an analysis of the exact potentials in the context of trajectory surface hopping procedure is presented, including preliminary investigations of velocity-adjustment, and the force-induced decoherence effect.

preprint2014arXiv

Ab-initio theory of Iron based superconductors

We report the first-principles study of superconducting critical temperature and superconducting properties of Fe-based superconductors taking into account on the same footing phonon, charge and spin-fluctuation mediated Cooper pairing. We show that in FeSe this leads to a modulated s$\pm$ gap symmetry, and that the antiferromagnetic paramagnons are the leading mechanism for superconductivity in FeSe, overcoming the strong repulsive effect of both phonons and charge pairing.

preprint2014arXiv

Almost Exact Exchange At Almost No Cost

A recently developed semiclassical approximation to exchange in one dimension is shown to be almost exact, with essentially no computational cost. The variational stability of this approximation is tested, and its far greater accuracy relative to local density functional calculations demonstrated. Even a fully orbital-free potential-functional calculation (no orbitals of any kind) yields little error relative to exact exchange, for more than one orbital.

preprint2014arXiv

An efficient algorithm for time propagation within time-dependent density functional theory

An efficient algorithm for time propagation of the time-dependent Kohn-Sham equations is presented. The algorithm is based on dividing the Hamiltonian into small time steps and assuming that it is constant over these steps. This allows for the time-propagating Kohn-Sham wave function to be expanded in the instantaneous eigenstates of the Hamiltonian. The stability and efficiency of the algorithm are tested not just for non-magnetic but also for fully non-collinear magnetic systems. We show that even for delicate properties, like magnetization density, large time-step sizes can be used indicating the stability and efficiency of the algorithm.

preprint2014arXiv

Classical nuclear motion coupled to electronic non-adiabatic transitions

We present a detailed derivation and numerical tests of a new mixed quantum-classical scheme to deal with non-adiabatic processes. The method is presented as the zero-th order approximation to the exact coupled dynamics of electrons and nuclei offered by the factorization of the electron-nuclear wave function [A. Abedi, N. T. Maitra and E. K. U. Gross, Phys. Rev. Lett., 105 (2010)]. Numerical results are presented for a model system for non-adiabatic charge transfer in order to test the performance of the method and to validate the underlying approximations.

preprint2014arXiv

Controlling observables in time-dependent quantum transport

The theory of time-dependent quantum transport addresses the question: How do electrons flow through a junction under the influence of an external perturbation as time goes by? In this paper, we invert this question and search for a time-dependent bias such that the system behaves in a desired way. This can, for example, be an observable that is forced to follow a certain pattern or the minimization of an objective function which depends on the observables. Our system of choice consists of quantum dots coupled to normal or superconducting leads. We present results for junctions with normal leads where the current, the density or a molecular vibration is optimized to follow a given target pattern. For junctions with two superconducting leads, where the Josephson effect triggers the current to oscillate, we show how to suppress the Josephson oscillations by suitably tailoring the bias. In a second example involving superconductivity, we consider a Y shaped junction with two quantum dots coupled to one superconducting and two normal leads. This device is used as a Cooper pair splitter to create entangled electrons on the two quantum dots. We maximize the splitting efficiency with the help of an optimized bias.

preprint2014arXiv

Enhancing and controlling single-atom high-harmonic generation spectra: a time-dependent density-functional scheme

High harmonic generation (HHG) provides a flexible framework for the development of coherent light sources in the extreme-ultraviolet and soft x-ray regimes. However it suffers from low conversion efficiencies as the control of the HHG spectral and temporal characteristics requires manipulating electron trajectories on attosecond time scale. The phase matching mechanism has been employed to selectively enhance specific quantum paths leading to HHG. A few important fundamental questions remain open, among those how much of the enhancement can be achieved by the single-emitter and what is the role of correlations (or the electronic structure) in the selectivity and control of HHG generation. Here we address those questions by examining computationally the possibility of optimizing the HHG spectrum of isolated Hydrogen and Helium atoms by shaping the slowly varying envelope of a 800 nm, 200-cycles long laser pulse. The spectra are computed with a fully quantum mechanical description, by explicitly computing the time-dependent dipole moment of the systems using a first-principles time-dependent density-functional approach (exact for the case of H). The sought optimization corresponds to the selective enhancement of single harmonics, which we find to be significant. This selectivity is entirely due to the single atom response, and not due to any propagation or phase-matching effect. In fact, this single-emitter enhancement adds to the phase-matching techniques to achieving even larger HHG enhancement factors. Moreover, we see that the electronic correlation plays a role in the determining the degree of optimization that can be obtained.

preprint2014arXiv

First-principles study of superconducting Rare-earth doped CaFe2As2

We report a systematic and ab-initio electronic structure calculation of Ca0.75 M0.25 Fe2 As2 with M = Ca, Sr, Eu, La, Ce, Pr, Nd, Pm, Sm, Na, K, Rb. The recently reported experimentally observed structural trends in rare earths-doped CaFe2 As2 compounds are successfully predicted and a complete theoretical description of the pressure induced orthorhombic to collapsed tetragonal transition is given. We demonstrate that the transition pressure is reduced by electron doping and rises linearly with the ionic size of the dopants. We discuss the implications of our description for the realization of a superconducting phase.

preprint2014arXiv

How to interpret the spectral density of the Keldysh nonequilibrium Green's function

This paper is devoted to the study and interpretation of the spectral function $\mathbf{A}(ω, T)$ of the Keldysh nonequilibrium Green's function. The spatial diagonal of the spectral function is often interpreted as a time-dependent local density of states. We show that this object can take negative values implying that a simple probability interpretation as a time-dependent density of states is not possible. The same issue also occurs for the Wigner function $P(x,p)$ where it is solved by taking the uncertainty principle into account. We follow the same path and incorporate the time-energy uncertainty relation to define a convoluted spectral function that allows for a probability interpretation. The usefulness of this quantity as a interpretative tool is demonstrated by visualizing the charge dynamics in a quantum dot coupled to superconducting leads.

preprint2014arXiv

How to represent crystal structures for machine learning: towards fast prediction of electronic properties

High-throughput density-functional calculations of solids are extremely time consuming. As an alternative, we here propose a machine learning approach for the fast prediction of solid-state properties. To achieve this, LSDA calculations are used as training set. We focus on predicting metallic vs. insulating behavior, and on predicting the value of the density of electronic states at the Fermi energy. We find that conventional representations of the input data, such as the Coulomb matrix, are not suitable for the training of learning machines in the case of periodic solids. We propose a novel crystal structure representation for which learning and competitive prediction accuracies become possible within an unrestricted class of spd systems. Due to magnetic phenomena learning on d systems is found more difficult than in pure sp systems.

preprint2014arXiv

Is the molecular Berry phase an artifact of the Born-Oppenheimer approximation?

We demonstrate that the molecular Berry phase and the corresponding non-analyticity in the electronic Born-Oppenheimer wavefunction is, in general, not a true topological feature of the exact solution of the full electron-nuclear Schrödinger equation. For a numerically exactly solvable model we show that a non-analyticity, and the associated geometric phase, only appear in the limit of infinite nuclear mass, while a perfectly smooth behavior is found for any finite nuclear mass.

preprint2014arXiv

Laser induced ultrafast demagnetization: an \emph{ab-initio} perspective

Time-dependent density functional theory is implemented in an all electron solid-state code for the case of fully non-collinear spins. We use this to study laser induced demagnetization in Fe, Co and Ni. It is shown that this demagnetization is a two-step process: excitation of a fraction of electrons followed by spin-flip transitions of the remaining localized electrons. These results successfully explain several experimental features such as the time-lag between the start of the pulse and demagnetization and spin-flip excitations dominating the physics. We further show that it is possible to control the moment loss by tunable laser parameters like frequency, duration and intensity.

preprint2014arXiv

Mixed quantum-classical dynamics from the exact decomposition of electron-nuclear motion

We present a novel mixed quantum-classical approach to the coupled electron-nuclear dynamics based on the exact factorization of the electron-nuclear wave function, recently proposed in [A. Abedi, N. T. Maitra, and E. K. U. Gross, Phys. Rev. Lett. 105, 123002 (2010)]. In this framework, classical nuclear dynamics is derived as the lowest order approximation of the time dependent Schrödinger equation that describes the evolution of the nuclei. The effect of the time dependent scalar and vector potentials, representing the exact electronic back-reaction on the nuclear subsystem, is consistently derived within the classical approximation. We examine with an example the performance of the proposed mixed quantum-classical scheme in comparison with exact calculations.

preprint2014arXiv

Optimizing the Cooper pair splitting efficiency in a Y-shaped junction

This letter is devoted to the optimization of the Cooper pair splitting efficiency in a Y-shaped junction. The latter consists of two quantum dots, one superconducting and two normal leads. We tailor the bias in the two normal leads such that the Cooper pairs leaving the superconductor are split up resulting in entangled electrons, one on each quantum dot. We are able to achieve a splitting efficiency of more than 99% which is significantly better than the efficiencies obtained in experiments so far.

preprint2014arXiv

Superconducting pairing mediated by spin-fluctuations from first principles

We present the derivation of an ab-initio and parameter free effective electron-electron interaction that goes beyond the screened RPA and accounts for superconducting pairing driven by spin-fluctuations. The construction is based on many body perturbation theory and relies on the approximation of the exchange-correlation part of the electronic self-energy within time dependent density functional theory. This effective interaction is included in an exchange correlation kernel for superconducting density functional theory, in order to achieve a completely parameter free superconducting gap equation. First results from applying the new functional to a simplified two-band electron gas model are consistent with experiments.

preprint2014arXiv

The virial theorem and exact properties of density functionals for periodic systems

In the framework of density functional theory, scaling and the virial theorem are essential tools for deriving exact properties of density functionals. Preexisting mathematical difficulties in deriving the virial theorem via scaling for periodic systems are resolved via a particular scaling technique. This methodology is employed to derive universal properties of the exchange-correlation energy functional for periodic systems.

preprint2014arXiv

Thermal Anharmonic Effects in PbTe from First Principles

We investigate the harmonic and anharmonic contributions to the phonon spectrum of lead telluride, and perform a complete characterization of how the anharmonic effects dominate the phonons in PbTe as temperature increases. This effect is the strongest factor in the favorable thermoelectric properties of PbTe: an optical-acoustic phonon band crossing reduces the speed of sound and the intrinsic thermal conductivity. We present the detailed temperature dependence of the dispersion relation and compare our calculated neutron scattering cross section with recent experimental measurements. We analyze the thermal resistivity's variation with temperature and clarify misconceptions about existing experimental literature. This quantitative prediction opens the way to phonon phase space engineering, to tailor the lifetimes of crucial heat carrying phonons.

preprint2013arXiv

A discontinuous functional for linear response time-dependent density functional theory: the exact-exchange kernel and approximate forms

We present a detailed study of the exact-exchange (EXX) kernel of time-dependent density functional theory with an emphasis on its discontinuity at integer particle numbers. It was recently found that this exact property leads to sharp peaks and step features in the kernel that diverge in the dissociation limit of diatomic systems [Hellgren and Gross, Phys. Rev. A, 022514 (2012)]. To further analyze the discontinuity of the kernel we here make use of two different approximations to the EXX kernel: the PGG approximation and a common energy denominator approximation (CEDA). It is demonstrated that whereas the PGG approximation neglects the discontinuity the CEDA includes it explicitly. By studying model molecular systems it is shown that the so-called field counter-acting effect in the density functional description of molecular chains can be viewed in terms of the discontinuity of the static kernel. The role of the frequency dependence is also investigated, highlighting its importance for long-range charge transfer excitations as well as inner-shell excitations.

preprint2013arXiv

Bootstrap kernel for organic low dimensional systems; PPV, pentacene and picene

We apply the bootstrap kernel within time dependent density functional theory to study one-dimensional chain of organic polymer poly-phenylene-vinylene and molecular crystals of picene and pentacene. The behaviour of this kernel in the presence and absence of local field effects is studied. The absorption spectra of poly-phenylene-vinylene has a bound excitonic peak which is well reproduced by the bootstrap kernel. Pentacene and picene, electronically similar materials, have remarkably different excitonic physics which is also captured properly by the bootstrap kernel. Inclusion of local-field effects dramatically change the spectra for both picene and pentacene. We highlight the reason behind this change. This also sheds light on the reasons behind the discrepancy in results between two different previous Bethe-Salpeter calculations.

preprint2013arXiv

Exact Electronic Potentials in Coupled Electron-Ion Dynamics

We develop a novel approach to the coupled motion of electrons and ions that focuses on the dynamics of the electronic subsystem. Usually the description of electron dynamics involves an electronic Schrödinger equation where the nuclear degrees of freedom appear as parameters or as classical trajectories. Here we derive the exact Schrödinger equation for the subsystem of electrons, staying within a full quantum treatment of the nuclei. This exact Schrödinger equation features a time-dependent potential energy surface for electrons (e-TDPES). We demonstrate that this exact e-TDPES differs significantly from the electrostatic potential produced by classical or quantum nuclei.

preprint2013arXiv

Mixed quantum-classical dynamics on the exact time-dependent potential energy surface: A fresh look at non-adiabatic processes

The exact nuclear time-dependent potential energy surface arises from the exact decomposition of electronic and nuclear motion, recently presented in [A. Abedi, N. T. Maitra, and E. K. U. Gross, Phys. Rev. Lett. 105, 123002 (2010)]. Such time-dependent potential drives nuclear motion and fully accounts for the coupling to the electronic subsystem. We investigate the features of the potential in the context of electronic non-adiabatic processes and employ it to study the performance of the classical approximation on nuclear dynamics. We observe that the potential, after the nuclear wave-packet splits at an avoided crossing, develops dynamical steps connecting different regions, along the nuclear coordinate, in which it has the same slope as one or the other adiabatic surface. A detailed analysis of these steps is presented for systems with different non-adiabatic coupling strength. The exact factorization of the electron-nuclear wave-function is at the basis of the decomposition. In particular, the nuclear part is the true nuclear wave-function, solution of a time-dependent Schroedinger euqation and leading to the exact many-body density and current density. As a consequence, the Ehrenfest theorem can be extended to the nuclear subsystem and Hamiltonian, as discussed here with an analytical derivation and numerical results.

preprint2013arXiv

Optimal control of strong-field ionization with time-dependent density-functional theory

We show that quantum optimal control theory (OCT) and time-dependent density-functional theory (TDDFT) can be combined to provide realistic femtosecond laser pulses for an enhanced ionization yield in many-electron systems. Using the H$_2$-molecule as a test case, the optimized laser pulse from the numerically exact scheme is compared to pulses obtained from OCT+TDDFT within the TD exact-exchange (TDEXX) and the TD local-density approximation (TDLDA). We find that the TDDFT-pulses produces an ionization yield of up to 50% when applied to the exact system. In comparison, pulses with a single frequency but the same fluence typically reach to yields around 5-15%, unless the frequency is carefully tuned into a Fano-type resonance that leads to $\sim 30%$ yield. On the other hand, optimization within the exact system alone leads to yields higher than 80%, demonstrating that correlation effects beyond the TDEXX and TDLDA can give rise to even more efficient ionization mechanisms.

preprint2013arXiv

Optimal control theory for quantum-classical systems: Ehrenfest Molecular Dynamics based on time-dependent density-functional theory

We derive the fundamental equations of an optimal control theory for systems containing both quantum electrons and classical ions. The system is modeled with Ehrenfest dynamics, a non-adiabatic variant of molecular dynamics. The general formulation, that needs the fully correlated many-electron wave function, can be simplified by making use of time-dependent density-functional theory. In this case, the optimal control equations require some modifications that we will provide. The abstract general formulation is complemented with the simple example of the H$_2^+$ molecule in the presence of a laser field.

preprint2013arXiv

Potential functionals versus density functionals

Potential functional approximations are an intriguing alternative to density functional approximations. The potential functional that is dual to the Lieb density functional is defined and properties given. The relationship between Thomas-Fermi theory as a density functional and as a potential functional is derived. The properties of several recent semiclassical potential functionals are explored, especially in their approach to the large particle number and classical continuum limits. The lack of ambiguity in the energy density of potential functional approximations is demonstrated. The density-density response function of the semiclassical approximation is calculated and shown to violate a key symmetry condition.

preprint2013arXiv

Reply to "Comment on 'Correlated electron-nuclear dynamics: Exact factorization of the molecular wavefunction' [J. Chem. Phys. 137, 22A530 (2012) ]''

This is a reply to the comment by J. L. Alonso et al. (arXiv:1302.6146; to appear in J. Chem. Phys) on our paper A. Abedi, N. T. Maitra, and E. K. U. Gross, J. Chem. Phys. 137, 22A530 (2012) [arXiv: 1208.4388]

preprint2013arXiv

Steps in the exact time-dependent potential energy surface

We study the exact Time-Dependent Potential Energy Surface (TDPES) in the presence of strong non-adiabatic coupling between the electronic and nuclear motion. The concept of the TDPES emerges from the exact factorization of the full electron-nuclear wave-function [A. Abedi, N. T. Maitra, and E. K. U. Gross, Phys. Rev. Lett. \textbf{105}, 123002 (2010)]. Employing a 1D model-system, we show that the TDPES exhibits a dynamical step that bridges between piecewise adiabatic shapes. We analytically investigate the position of the steps and the nature of the switching between the adiabatic pieces of the TDPES.

preprint2013arXiv

Thermal Density Functional Theory in Context

This chapter introduces thermal density functional theory, starting from the ground-state theory and assuming a background in quantum mechanics and statistical mechanics. We review the foundations of density functional theory (DFT) by illustrating some of its key reformulations. The basics of DFT for thermal ensembles are explained in this context, as are tools useful for analysis and development of approximations. We close by discussing some key ideas relating thermal DFT and the ground state. This review emphasizes thermal DFT's strengths as a consistent and general framework.

preprint2013arXiv

Transverse spin gradient functional for non-collinear spin density functional theory

We present a novel functional for spin density functional theory aiming at the description of non-collinear magnetic structures. The construction of the functional employs the spin-spiral-wave state of the uniform electron gas as reference system. We show that the functional depends on transverse gradients of the spin magnetization, i.e., in contrast to the widely used local spin density approximation, the functional is sensitive to local changes of the direction of the spin magnetization. As a consequence the exchange-correlation magnetic field is not parallel to the spin magnetization and a local spin-torque is present in the ground state of the Kohn-Sham system. As a proof-of-principle we apply the functional to a Chromium mono-layer in the non-collinear ${120^{\circ}}$-N{é}el state.

preprint2012arXiv

Ab-initio angle and energy resolved photoelectron spectroscopy with time-dependent density-functional theory

We present a time-dependent density-functional method able to describe the photoelectron spectrum of atoms and molecules when excited by laser pulses. This computationally feasible scheme is based on a geometrical partitioning that efficiently gives access to photoelectron spectroscopy in time-dependent density-functional calculations. By using a geometrical approach, we provide a simple description of momentum-resolved photoe- mission including multi-photon effects. The approach is validated by comparison with results in the literature and exact calculations. Furthermore, we present numerical photoelectron angular distributions for randomly oriented nitrogen molecules in a short near infrared intense laser pulse and helium-(I) angular spectra for aligned carbon monoxide and benzene.

preprint2012arXiv

Correlated electron-nuclear dynamics: Exact factorization of the molecular wavefunction

It was recently shown [Phys. Rev. Lett. 105, 123002 (2010)] that the complete wavefunction for a system of electrons and nuclei evolving in a time-dependent external potential can be exactly factorized into an electronic wavefunction and a nuclear wavefunction. The concepts of an exact time-dependent potential energy surface (TDPES) and exact time-dependent vector potential emerge naturally from the formalism. Here we present a detailed description of the formalism, including a full derivation of the equations that the electronic and nuclear wavefunctions satisfy. We demonstrate the relationship of this exact factorization to the traditional Born-Oppenheimer expansion. A one-dimensional model of the H$_2^+$ molecule in a laser field shows the usefulness of the exact TDPES in interpreting coupled electron-nuclear dynamics: we show how features of its structure indicate the mechanism of dissociation. We compare the exact TDPES with potential energy surfaces from the time-dependent Hartree-approach, and also compare traditional Ehrenfest dynamics with Ehrenfest dynamics on the exact TDPES.

preprint2012arXiv

Correlation potentials for molecular bond dissociation within the self-consistent random phase approximation

Self-consistent correlation potentials for H$_2$ and LiH for various inter-atomic separations are obtained within the random phase approximation (RPA) of density functional theory. The RPA correlation potential shows a peak at the bond midpoint, which is an exact feature of the true correlation potential, but lacks another exact feature: the step important to preserve integer charge on the atomic fragments in the dissociation limit. An analysis of the RPA energy functional in terms of fractional charge is given which confirms these observations. We find that the RPA misses the derivative discontinuity at odd integer particle numbers but explicitly eliminates the fractional spin error in the exact-exchange functional. The latter finding explains the accurate total energy in the dissociation limit.

preprint2012arXiv

Discontinuities of the exchange-correlation kernel and charge-transfer excitations in time-dependent density functional theory

We identify the key property that the exchange-correlation (XC) kernel of time-dependent density functional theory must have in order to describe long-range charge-transfer excitations. We show that the discontinuity of the XC potential as a function of particle number induces a space -and frequency-dependent discontinuity of the XC kernel which diverges as $r\to\infty$. In a combined donor-acceptor system, the same discontinuity compensates for the vanishing overlap between the acceptor and donor orbitals, thereby yielding a finite correction to the Kohn-Sham eigenvalue differences. This mechanism is illustrated to first order in the Coulomb interaction.

preprint2012arXiv

Effect of discontinuities in Kohn-Sham-based chemical reactivity theory

We provide a new derivation of a formula for the Fukui function of density-functional chemical reactivity theory which incorporates the discontinuities in the Kohn-Sham reference system. Orbital relaxations are described in terms of the exchange-correlation (XC) kernel, i.e., the derivative of the XC potential with respect to the density and it is shown that in order to correctly measure the reactivity toward a nucleophilic reagent a discontinuity of the XC kernel has to be taken into account. The importance of this finding is illustrated in model molecular systems.

preprint2012arXiv

Enhanced excitonic effects in the energy loss spectra of LiF and Ar at large momentum transfer

It is demonstrated that the bootstrap kernel [\onlinecite{sharma11}] for finite values of ${\bf q}$ crucially depends upon the matrix character of the kernel and gives results of the same good quality as in the ${\bf q} \rightarrow 0$ limit. The bootstrap kernel is further used to study the electron loss as well as absorption spectra for Si, LiF and Ar for various values of ${\bf q}$. The results show that the excitonic effects in LiF and Ar are enhanced for values of ${\bf q}$ away from the $Γ$-point. The reason for this enhancement is the interaction between the exciton and high energy inter-band electron-hole transitions. This fact is validated by calculating the absorption spectra under the influence of an external electric field. The electron energy loss spectra is shown to change dramatically as a function of ${\bf q}$.

preprint2012arXiv

Ionization potentials and electron affinities from reduced density matrix functional theory

In the recent work of S. Sharma \emph{et al.}, (arxiv.org: arxiv:0912.1118), a single-electron spectrum associated with the natural orbitals was defined as the derivative of the total energy with respect to the occupation numbers at half filling for the orbital of interest. This idea reproduces the bands of various periodic systems using the appropriate functional quite accurately. In the present work we apply this approximation to the calculation of the ionization potentials and electron affinities of molecular systems using various functionals within the reduced density-matrix functional theory. We demonstrate that this approximation is very successful in general and in particular for certain functionals it performs better than the direct determination of the ionization potentials and electron affinities through the calculation of positive and negative ions respectively. The reason for this is identified to be the inaccuracy that arises from different handling of the open- and closed-shell systems.

preprint2012arXiv

Minimization procedure in reduced density matrix functional theory by means of an effective noninteracting system

In this work, we propose a self-consistent minimization procedure for functionals in reduced density matrix functional theory. We introduce an effective noninteracting system at finite temperature which is capable of reproducing the groundstate one-reduced density matrix of an interacting system at zero temperature. By introducing the concept of a temperature tensor the minimization with respect to the occupation numbers is shown to be greatly improved.

preprint2012arXiv

Multiplicity of solutions to GW-type approximations

We show that the equations underlying the $GW$ approximation have a large number of solutions. This raises the question: which is the physical solution? We provide two theorems which explain why the methods currently in use do, in fact, find the correct solution. These theorems are general enough to cover a large class of similar algorithms. An efficient algorithm for including self-consistent vertex corrections well beyond $GW$ is also described and further used in numerical validation of the two theorems.

preprint2012arXiv

Phononic Self energy effects and superconductivity in CaC$_6$

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.

preprint2012arXiv

Reduced density matrix functional theory at finite temperature. II. Application to the electron gas: Exchange only

Using the newly introduced theory of finite-temperature reduced density matrix functional theory, we apply the first-order approximation to the homogeneous electron gas. We consider both collinear spin states as well as symmetry broken states describing planar spin spirals and investigate the magnetic phase diagram as well as the temperature-dependence of the single particle spectra.

preprint2012arXiv

Reduced density matrix functional theory at finite temperature. III. Application to the electron gas: Correlation effects

Based on our derivation of finite temperature reduced density matrix functional theory and the discussion of the performance of its first-order functional this work presents several different correlation-energy functionals and applies them to the homogeneous electron gas. The zero temperature limits of the correlation-energy and the momentum distributions are investigated and the magnetic phase diagrams in collinear spin configuration are discussed.

preprint2012arXiv

Spectral density and metal-insulator phase transition in Mott insulators within RDMFT

We present a method for calculating the spectrum of periodic solids within reduced density matrix functional theory. This method is validated by a detailed comparison of the angular momentum projected spectral density with that of well established many-body techniques, in all cases finding an excellent agreement. The physics behind the pressure induced insulator-metal phase transition in MnO is investigated. The driving mechanism of this transition is identified as increased crystal field splitting with pressure, resulting in a charge redistribution between the Mn $e_g$ and $t_2g$ symmetry projected states.

preprint2012arXiv

Spectrum of extended systems from Reduced Density Matrix Functional Theory

We present a method for calculating the spectrum of extended solids within reduced density matrix functional theory. An application of this method to the strongly correlated transition metal oxide series demonstrates that (i) an insulating state is found in the absence of magnetic order and, in addition, (ii) the interplay between the change transfer and Mott-Hubbard correlation is correctly described. In this respect we find that while NiO has a strong charge transfer character to the electronic gap, with substantial hybridization between $t_{2g}$ and oxygen-$p$ states in the lower Hubbard band, for MnO this is almost entirely absent. As a validation of our method we also calculate the spectra for a variety of weakly correlated materials, finding good agreement with experiment and other techniques.

preprint2011arXiv

Bootstrap approximation for the exchange-correlation kernel of time-dependent density functional theory

A new parameter-free approximation for the exchange-correlation kernel $f_{\rm xc}$ of time-dependent density functional theory is proposed. This kernel is expressed as an algorithm in which the exact Dyson equation for the response as well as a further approximate condition are solved together self-consistently leading to a simple parameter-free kernel. We apply this to the calculation of optical spectra for various small bandgap (Ge, Si, GaAs, AlN, TiO$_2$, SiC), large bandgap (C, LiF, Ar, Ne) and magnetic (NiO) insulators. The calculated spectra are in very good agreement with experiment for this diverse set of materials, highlighting the universal applicability of the new kernel.

preprint2011arXiv

Comparative study of many-body perturbation theory and time-dependent density functional theory in the out-of-equilibrium Anderson model

We study time-dependent electron transport through an Anderson model. The electronic interactions on the impurity site are included via the self-energy approximations at Hartree-Fock (HF), second Born (2B), GW, and T-Matrix level as well as within a time-dependent density functional (TDDFT) scheme based on the adiabatic Bethe-Ansatz local density approximation (ABALDA) for the exchange correlation potential. The Anderson model is driven out of equilibrium by applying a bias to the leads and its nonequilibrium dynamics is determined by real-time propagation. The time-dependent currents and densities are compared to benchmark results obtained with the time-dependent density matrix renormalization group (tDMRG) method. Many-body perturbation theory beyond HF gives results in close agreement with tDMRG especially within the 2B approximation. We find that the TDDFT approach with the ABALDA approximation produces accurate results for the densities on the impurity site but overestimates the currents. This problem is found to have its origin in an overestimation of the lead densities which indicates that the exchange correlation potential must attain nonzero values in the leads.

preprint2011arXiv

Correlation effects in bistability at the nanoscale: steady state and beyond

The possibility of finding multistability in the density and current of an interacting nanoscale junction coupled to semi-infinite leads is studied at various levels of approximation. The system is driven out of equilibrium by an external bias and the non-equilibrium properties are determined by real-time propagation using both time-dependent density functional theory (TDDFT) and many-body perturbation theory (MBPT). In TDDFT the exchange-correlation effects are described within a recently proposed adiabatic local density approximation (ALDA). In MBPT the electron-electron interaction is incorporated in a many-body self-energy which is then approximated at the Hartree-Fock (HF), second-Born (2B) and GW level. Assuming the existence of a steady-state and solving directly the steady-state equations we find multiple solutions in the HF approximation and within the ALDA. In these cases we investigate if and how these solutions can be reached through time evolution and how to reversibly switch between them. We further show that for the same cases the inclusion of dynamical correlation effects suppresses bistability.

preprint2011arXiv

Electronic structure via potential functional approximations

The universal functional of Hohenberg-Kohn is given as a coupling-constant integral over the density as a functional of the potential. Conditions are derived under which potential-functional approximations are variational. Construction via this method and imposition of these conditions are shown to greatly improve the accuracy of the non-interacting kinetic energy needed for orbital-free Kohn-Sham calculations.

preprint2011arXiv

Exact conditions and scaling relations in finite temperature density functional theory

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.

preprint2011arXiv

Optimization Schemes for Selective Molecular Cleavage with Tailored Ultrashort Laser Pulses

We present some approaches to the computation of ultra-fast laser pulses capable of selectively breaking molecular bonds. The calculations are based on a mixed quantum-classical description: The electrons are treated quantum mechanically (making use of time-dependent density-functional theory), whereas the nuclei are treated classically. The temporal shape of the pulses is tailored to maximise a control target functional which is designed to produce the desired molecular cleavage. The precise definition of this functional is a crucial ingredient: we explore expressions based on the forces, on the momenta and on the velocities of the nuclei. The algorithm used to find the optimum pulse is also relevant; we test both direct gradient-free algorithms, as well as schemes based on formal optimal control theory. The tests are performed both on one dimensional models of atomic chains, and on first-principles descriptions of molecules.

preprint2011arXiv

Real-time switching between multiple steady-states in quantum transport

We study transport through an interacting model system consisting of a central correlated site coupled to finite bandwidth tight-binding leads, which are considered as effectively noninteracting. Its nonequilibrium properties are determined by real-time propagation of the Kadanoff-Baym equations after applying a bias voltage to the system. The electronic interactions on the central site are incorporated by means of self-energy approximations at Hartree-Fock, second Born and GW level. We investigate the conditions under which multiple steady-state solutions occur within different self-energy approximations, and analyze in detail the nature of these states from an analysis of their spectral functions. At the Hartree-Fock level at least two stable steady-state solutions with different densities and currents can be found. By applying a gate voltage-pulse at a given time we are able to switch between these solutions. With the same parameters we find only one steady-state solution when the self-consistent second Born and GW approximations are considered. We therefore conclude that treatment of many-body interactions beyond mean-field can destroy bistability and lead to qualitatively different results as compared those at mean-field level.

preprint2011arXiv

Time-dependent density-functional and reduced density-matrix methods for few electrons: Exact versus adiabatic approximations

To address the impact of electron correlations in the linear and non-linear response regimes of interacting many-electron systems exposed to time-dependent external fields, we study one-dimensional (1D) systems where the interacting problem is solved exactly by exploiting the mapping of the 1D $N$-electron problem onto an $N$-dimensional single electron problem. We analyze the performance of the recently derived 1D local density approximation as well as the exact-exchange orbital functional for those systems. We show that the interaction with an external resonant laser field shows Rabi oscillations which are detuned due to the lack of memory in adiabatic approximations. To investigate situations where static correlations play a role, we consider the time-evolution of the natural occupation numbers associated to the reduced one-body density matrix. Those studies shed light on the non-locality and time-dependence of the exchange and correlation functionals in time-dependent density and density-matrix functional theories.

preprint2010arXiv

Exact factorization of the time-dependent electron-nuclear wavefunction

We present an exact decomposition of the complete wavefunction for a system of nuclei and electrons evolving in a time-dependent external potential. We derive formally exact equations for the nuclear and electronic wavefunctions that lead to rigorous definitions of a time-dependent potential energy surface (TDPES) and a time-dependent geometric phase. For the $H_2^+$ molecular ion exposed to a laser field, the TDPES proves to be a useful interpretive tool to identify different mechanisms of dissociation.

preprint2010arXiv

Static and Dynamical Susceptibility of LaO1-xFxFeAs

The mechanism of superconductivity and magnetism and their possible interplay have recently been under debate in pnictides. A likely pairing mechanism includes an important role of spin fluctuations and can be expressed in terms of the magnetic susceptibility chi. The latter is therefore a key quantity in the determination of both the magnetic properties of the system in the normal state, and of the contribution of spin fluctuations to the pairing potential. A basic ingredient to obtain chi is the independent-electron susceptibility chi0. Using LaO1-xFxFeAs as a prototype material, in this report we present a detailed ab-initio study of chi0(q,omega), as a function of doping and of the internal atomic positions. The resulting static chi0(q,0) is consistent with both the observed M-point related magnetic stripe phase in the parent compound, and with the existence of incommensurate magnetic structures predicted by ab-initio calculations upon doping.

preprint2010arXiv

Theoretical Investigation of Optical Conductivity in Ba [Fe(1-x)Co(x)]2 As2

We report on theoretical calculations of the optical conductivity of Ba [Fe(1-x)Co(x)]2 As2, as obtained from density functional theory within the full potential LAPW method. A thorough comparison with experiment shows that we are able to reproduce most of the observed experimental features, in particular a magnetic peak located at about 0.2 eV which we ascribe to antiferromagnetic ordered magnetic stripes. We also predict a large in-plane anisotropy of this feature, which agrees very well with measurements on detwinned crystals. The effect of Co doping as well as the dependence of plasma frequency on the magnetic order is also investigated.

preprint2009arXiv

Femtosecond laser pulse shaping for enhanced ionization

We demonstrate how the shape of femtosecond laser pulses can be tailored in order to obtain maximal ionization of atoms or molecules. For that purpose, we have overlayed a direct-optimization scheme on top of a fully unconstrained computation of the three-dimensional time-dependent Schrodinger equation. The procedure looks for pulses that maintain the same total length and integrated intensity or fluence as a given pulse that serves as an initial guess. It allows, however, for changes in frequencies -- within a certain, predefined range -- and overall shape, leading to enhanced ionization. We illustrate the scheme by calculating ionization yields for the H2+ molecule when irradiated with short (~5 fs), high-intensity laser pulses.

preprint2009arXiv

Non-collinear spin-spiral phase for the uniform electron gas within Reduced-Density-Matrix-Functional Theory

The non-collinear spin-spiral density wave of the uniform electron gas is studied in the framework of Reduced-Density-Matrix-Functional Theory. For the Hartree-Fock approximation, which can be obtained as a limiting case of Reduced-Density-Matrix-Functional Theory, Overhauser showed a long time ago that the paramagnetic state of the electron gas is unstable with respect to the formation of charge or spin density waves. Here we not only present a detailed numerical investigation of the spin-spiral density wave in the Hartree-Fock approximation but also investigate the effects of correlations on the spin-spiral density wave instability by means of a recently proposed density-matrix functional.

preprint2008arXiv

Density-Matrix-Power Functional: Performance for Finite Systems and the Homogeneous Electron Gas

An exchange correlation energy functional involving fractional power of the one-body reduced density matrix [Phys. Rev. B {\bf 78}, 201103 (2008)] is applied to finite systems and to the homogeneous electron gas. The performance of the functional is assessed for the correlation and atomization energies of the molecules contained in the G2 set and for the correlation energy of the homogeneous electron gas. High accuracy is found for these two very different types of systems.

preprint2008arXiv

Discontinuity of the chemical potential in RDMFT for open-shell systems

We employ reduced density-matrix functional theory in the calculation of the fundamental gap of open-shell systems. The formula for the calculation of the fundamental gap is derived with special attention to the spin of the neutral and the charged systems. We discuss the effects of different functionals as well as the changes due to different basis sets. Also, we investigate the importance of varying the natural orbitals for the calculation of the fundamental gap

preprint2008arXiv

Electronic, dynamical and superconducting properties of CaBeSi

We report first-principles calculations on the normal and superconducting state of CaBe(x)Si(2-x) (x=1), in the framework of density functional theory for superconductors (SCDFT). CaBeSi is isostructural and isoelectronic to MgB2 and this makes possible a direct comparison of the electronic and vibrational properties and the electron-phonon interaction of the two materials. Despite the many similarities with MgB2 (e.g. sigma bands at the Fermi level and a larger Fermi surface nesting), according to our calculations CaBeSi has a very low critical temperature (Tc ~ 0.4 K, consistent with the experiment). CaBeSi exhibits a complex gap structure, with three gaps at Fermi level: besides the two sigma and pi gaps, present also in MgB2, the appearance of a third gap is related to the anisotropy of the Coulomb repulsion, acting in different way on the bonding and antibonding electronic pi states.

preprint2008arXiv

Exchange-correlation orbital functionals in current-density-functional theory: Application to a quantum dot in magnetic fields

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.

preprint2008arXiv

Time-dependent natural orbitals and occupation numbers

We report equations of motion for the occupation numbers of natural spin orbitals and show that adiabatic extensions of common functionals employed in ground-state reduced-density-matrix-functional theory have the shortcoming of leading always to occupation numbers which are independent of time. We illustrate the exact time-dependence of the natural spin orbitals and occupation numbers for the case of electron-ion scattering and for atoms in strong laser fields. In the latter case, we observe strong variations of the occupation numbers in time.

preprint2007arXiv

Performance of one-body reduced density matrix functionals for the homogeneous electron gas

The subject of this study is the exchange-correlation-energy functional of reduced density matrix functional theory. Approximations of this functional are tested by applying them to the homogeneous electron gas. We find that two approximations recently proposed by Gritsenko, Pernal, and Baerends, J. Chem. Phys., {\bf 122}, 204102 (2005), yield considerably better correlation energies and momentum distributions than previously known functionals. We introduce modifications to these functionals which, by construction, reproduce the exact correlation energy of the homogeneous electron gas.

preprint2006arXiv

Optimized Effective Potential Method in Current-Spin Density Functional Theory

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.

E. K. U. Gross

What is connected

Connect this record

See the researcher in context

Building this map preview

89 published item(s)

Magnetism-Enhanced Strong Electron-Phonon Coupling in Infinite-Layer Nickelate

Energy, Momentum and Angular Momentum Transfer Between Electrons and Nuclei

From Kohn-Sham to many-electron energies via step structures in the exchange-correlation potential

Geometric energy transfer in two-component systems

Ab initio theory of plasmonic superconductivity within the Eliashberg and density-functional formalisms

Extending solid-state calculations to ultra long-range length scales

Many-body Green's function theory of electrons and nuclei beyond the Born-Oppenheimer approximation

How is the derivative discontinuity related to steps in the exact Kohn-Sham potential?

A single-electron approach for many-electron dynamics in high-order harmonic generation

An Exact Factorization Perspective on Quantum Interferences in Nonadiabatic Dynamics

Exact factorization-based density functional theory of electrons and nuclei

Molecular geometric phase from the exact electron-nuclear factorization

Quantum-Classical Non-Adiabatic Dynamics: Coupled- vs. Independent-Trajectory Methods

Ultrafast Local Magnetization and Demagnetization in Heusler Alloys

A coupled-trajectory quantum-classical approach to decoherence in non-adiabatic processes

Ab initio calculation of a Pb single layer on a Si substrate: two-dimensionality and superconductivity

Ab Initio Theory of Superconductivity in a Magnetic Field I. : Spin Density Functional Theory For Superconductors and Eliashberg Equations

Ab-Initio Theory of Superconductivity in a Magnetic Field II. : Numerical solution

Demonstration of Optimal Control of Laser Induced Spin-Orbit Mediated Ultrafast Demagnetization

Doping induced metal-insulator phase transition in NiO

High temperature superconductivity in sulfur and selenium hydrides at high pressure

L1$_0$ stacked binaries as candidates for hard magnets: FePt, MnAl and MnGa

Laser-induced electron localization in H$_2^+$: Mixed quantum-classical dynamics based on the exact time-dependent potential energy surface

Nuclear velocity perturbation theory for vibrational circular dichroism: An approach based on the exact factorization of the electron-nuclear wave function

Reduced Density Matrix Functional Theory at Finite Temperature: Theoretical Foundations

Spectrum for non-magnetic Mott insulators from power functional within Reduced Density Matrix Functional Theory

Superconductivity in metastable phases of phosphorus-hydride compounds under high pressure

The exact forces on classical nuclei in non-adiabatic charge transfer

Ab-initio theory of Iron based superconductors

Almost Exact Exchange At Almost No Cost

An efficient algorithm for time propagation within time-dependent density functional theory

Classical nuclear motion coupled to electronic non-adiabatic transitions

Controlling observables in time-dependent quantum transport

Enhancing and controlling single-atom high-harmonic generation spectra: a time-dependent density-functional scheme

First-principles study of superconducting Rare-earth doped CaFe2As2

How to interpret the spectral density of the Keldysh nonequilibrium Green's function

How to represent crystal structures for machine learning: towards fast prediction of electronic properties

Is the molecular Berry phase an artifact of the Born-Oppenheimer approximation?

Laser induced ultrafast demagnetization: an \emph{ab-initio} perspective

Mixed quantum-classical dynamics from the exact decomposition of electron-nuclear motion

Optimizing the Cooper pair splitting efficiency in a Y-shaped junction

Superconducting pairing mediated by spin-fluctuations from first principles

The virial theorem and exact properties of density functionals for periodic systems

Thermal Anharmonic Effects in PbTe from First Principles

A discontinuous functional for linear response time-dependent density functional theory: the exact-exchange kernel and approximate forms

Bootstrap kernel for organic low dimensional systems; PPV, pentacene and picene

Exact Electronic Potentials in Coupled Electron-Ion Dynamics

Mixed quantum-classical dynamics on the exact time-dependent potential energy surface: A fresh look at non-adiabatic processes

Optimal control of strong-field ionization with time-dependent density-functional theory

Optimal control theory for quantum-classical systems: Ehrenfest Molecular Dynamics based on time-dependent density-functional theory

Potential functionals versus density functionals

Reply to "Comment on 'Correlated electron-nuclear dynamics: Exact factorization of the molecular wavefunction' [J. Chem. Phys. 137, 22A530 (2012) ]''

Steps in the exact time-dependent potential energy surface

Thermal Density Functional Theory in Context

Transverse spin gradient functional for non-collinear spin density functional theory

Ab-initio angle and energy resolved photoelectron spectroscopy with time-dependent density-functional theory

Correlated electron-nuclear dynamics: Exact factorization of the molecular wavefunction

Correlation potentials for molecular bond dissociation within the self-consistent random phase approximation

Discontinuities of the exchange-correlation kernel and charge-transfer excitations in time-dependent density functional theory

Effect of discontinuities in Kohn-Sham-based chemical reactivity theory

Enhanced excitonic effects in the energy loss spectra of LiF and Ar at large momentum transfer

Ionization potentials and electron affinities from reduced density matrix functional theory

Minimization procedure in reduced density matrix functional theory by means of an effective noninteracting system

Multiplicity of solutions to GW-type approximations

Phononic Self energy effects and superconductivity in CaC$_6$

Reduced density matrix functional theory at finite temperature. II. Application to the electron gas: Exchange only

Reduced density matrix functional theory at finite temperature. III. Application to the electron gas: Correlation effects

Spectral density and metal-insulator phase transition in Mott insulators within RDMFT

Spectrum of extended systems from Reduced Density Matrix Functional Theory

Bootstrap approximation for the exchange-correlation kernel of time-dependent density functional theory

Comparative study of many-body perturbation theory and time-dependent density functional theory in the out-of-equilibrium Anderson model

Correlation effects in bistability at the nanoscale: steady state and beyond

Electronic structure via potential functional approximations

Exact conditions and scaling relations in finite temperature density functional theory