Search

JPhys photonics 4, 034003 (2022). doi:10.1088/2515-7647/ac6ea5 special issue: "Focus on Nanophotonics and Biophotonics for Biomedical and Environmental Applications", The dynamics of biologically relevant molecules exposed to ionizing radiation contains many facets and spans several orders of magnitude in time and energy. In the extreme ultraviolet (XUV) spectral range, multi-electronic phenomena and bands of correlated states with inner-valence holes must be accounted for in addition to a plethora of vibrational modes and available dissociation channels. The ability to track changes in charge density and bond length during ultrafast reactions is an important endeavor toward more general abilities to simulate and control photochemical processes, possibly inspired by those that have evolved biologically. By using attosecond XUV pulses extending up to 35 eV and few-femtosecond near-infrared pulses, we have previously time-resolved correlated electronic dynamics and charge migration occurring in the biologically relevant molecule adenine after XUV-induced sudden ionization. Here, using additional experimental data, we comprehensively report on both electronic and vibrational dynamics of this nucleobase in an energy range little explored to date with high temporal resolution. The time-dependent yields of parent and fragment ions in the mass spectra are analyzed to extract exponential time constants and oscillation periods. Together with time-dependent density functional theory and ab-initio Green's function methods, we identify different vibrational and electronic processes. Beyond providing further insights into the XUV-induced dynamics of an important nucleobase, our work demonstrates that yields of specific dissociation outcomes can be influenced by sufficiently well-timed ultrashort pulses, therefore providing a new route for the control of the multi-electronic and dissociative dynamics of a DNA building block., Published by IOP Publishing, Bristol

A new formalism to describe steady-state electronic and thermal transport in the framework of density functional theory is presented. A one-to-one correspondence is proven between the three basic variables of the theory, i.e., the density on as well as the electrical and heat currents through the junction, and the three basic potentials, i.e., the local potential in as well as the DC bias and thermal gradient across the junction. Consequently, the Kohn-Sham system of the theory requires three exchange-correlations potentials. In linear response, the new formalism leads to exact expressions for the many-body transport coefficients (both electrical and thermal conductances and Seebeck coefficient) in terms of both the corresponding Kohn-Sham coefficients and derivatives of the exchange-correlations potentials. The theory is applied to the Single Impurity Anderson Model, and an accurate analytic parametrization for these derivatives in the Coulomb blockade regime is constructed through reverse engineering., Comment: 11 pages, 5 figures

The interaction of electrons with quantized phonons and photons underlies the ultrafast dynamics of systems ranging from molecules to solids, and it gives rise to a plethora of physical phenomena experimentally accessible using time-resolved techniques. Green's function methods offer an invaluable interpretation tool since scattering mechanisms of growing complexity can be selectively incorporated in the theory. Currently, however, real-time Green's function simulations are either prohibitively expensive due to the cubic scaling with the propagation time or do neglect the feedback of electrons on the bosons, thus violating energy conservation. We put forward a computationally efficient Green's function scheme which overcomes both limitations. The numerical effort scales linearly with the propagation time while the simultaneous dressing of electrons and bosons guarantees the fulfillment of all fundamental conservation laws. We present a real-time study of the phonon-driven relaxation dynamics in an optically excited narrow band-gap insulator, highlighting the nonthermal behavior of the phononic degrees of freedom. Our formulation paves the way to first-principles simulations of electron-boson systems with unprecedented long propagation times., 7 pages, with additional supplementary materials

The ultrafast hole dynamics triggered by the photoexcitation of molecular targets is a highly correlated process even for those systems, such as organic molecules, having a weakly correlated ground state. We provide a unifying framework and a numerically efficient matrix formulation of state-of-the-art nonequilibrium Green's function (NEGF) methods such as second-Born as well as $GW$ and $T$-matrix without and with exchange diagrams. Numerical simulations are presented for a paradigmatic, exactly solvable molecular system, and the shortcomings of the established NEGF methods are highlighted. We then develop a NEGF scheme based on the Faddeev treatment of three-particle correlations; the exceptional improvement over established methods is explained and demonstrated. The Faddeev NEGF scheme scales linearly with the maximum propagation time, thereby opening up prospects for femtosecond simulations of large molecules.

The generalized Kadanoff-Baym ansatz (GKBA) offers a computationally inexpensive approach to simulate out-of-equilibrium quantum systems within the framework of nonequilibrium Green's functions. For finite systems the limitation of neglecting initial correlations in the conventional GKBA approach has recently been overcome [Phys. Rev. B 98, 115148 (2018)]. However, in the context of quantum transport the contacted nature of the initial state, i.e., a junction connected to bulk leads, requires a further extension of the GKBA approach. In this work, we lay down a GKBA scheme which includes initial correlations in a partition-free setting. In practice, this means that the equilibration of the initially correlated and contacted molecular junction can be separated from the real-time evolution. The information about the contacted initial state is included in the out-of-equilibrium calculation via explicit evaluation of the memory integral for the embedding self-energy, which can be performed without affecting the computational scaling with the simulation time and system size. We demonstrate the developed method in carbon-based molecular junctions, where we study the role of electron correlations in transient current signatures., Comment: 12 pages, 6 figures

We calculate the {\em exact} time-resolved ARPES spectrum of a two-band model semiconductor driven out of equilibrium by resonant and nonresonant laser pulses, highlighting the effects of phonon-induced decoherence and relaxation. {\em Resonant} excitations initially yield a replica of the {\em valence} band shifted upward by the energy of the exciton peak in photoabsorption. This phase is eventually destroyed by phonon-induced decoherence: the valence-band replica lowers in energy by the Stokes shift, locating at the energy of the exciton peak in photoluminescence, and its width grows due to phonon dressing. {\em Nonresonant} excitations initially yield a map of the conduction band. Then electrons transfer their excess energy to the lattice and bind with the holes left behind to form excitons. In this relaxed regime a replica of the {\em conduction} band appears inside the gap. At fixed momentum the lineshape of the conduction-band replica versus the photoelectron energy is proportional to the exciton wavefunction in "energy space" and it is highly asymmetric. Although the two-band model represents an oversimplified description of real materials the highlighted features are qualitative in nature; hence they provide useful insights into time-resolved ARPES spectra and their physical interpretation., Comment: 14 pages, 3 figures

Density functional theory can be generalized to mixtures of ground and excited states, for the purpose of determining energies of excitations using low-cost density functional approximations. Adapting approximations originally developed for ground states to work in the new setting would fast-forward progress enormously. But, previous attempts have stumbled on daunting fundamental issues. Here we show that these issues can be prevented from the outset, by using a fluctuation dissipation theorem (FDT) to dictate key functionals. We thereby show that existing exchange energy approximations are readily adapted to excited states, when combined with a rigorous exact Hartree term that is different in form from its ground state counterpart, and counterparts based on ensemble ansatze. Applying the FDT to correlation energies also provides insights into ground state-like and ensemble-only correlations. We thus provide a comprehensive and versatile framework for ensemble density functional approximations.

Density functional theory can be generalized to mixtures of ground and excited states, for the purpose of determining energies of excitations using low-cost density functional approximations. Adapting approximations originally developed for ground states to work in the new setting would fast-forward progress enormously. But, previous attempts have stumbled on daunting fundamental issues. Here we show that these issues can be prevented from the outset, by using a fluctuation dissipation theorem (FDT) to dictate key functionals. We thereby show that existing exchange energy approximations are readily adapted to excited states, when combined with a rigorous exact Hartree term that is different in form from its ground state counterpart, and counterparts based on ensemble ansatze. Applying the FDT to correlation energies also provides insights into ground state-like and ensemble-only correlations. We thus provide a comprehensive and versatile framework for ensemble density functional approximations.

The nontrivial topology of {\it p}-wave superfluids make these systems attractive candidate in information technology. In this work we report on the topological state of a {\it p}-wave nonequilibrium excitonic insulator (NEQ-EI) and show how to steer a nontopological band-insulator with bright $p$-excitons toward this state by a suitable laser pulse, thus achieving a dynamical topological phase transition. The underlying mechanism behind the transition is the broken gauge-symmetry of the NEQ-EI which causes self-sustained persistent oscillations of the excitonic condensate and hence a Floquet topological state for high enough exciton densities. We show the formation of Floquet Majorana modes at the boundaries of the open system and discuss topological spectral signatures for ARPES experiments. We emphasize that the topological properties of a {\it p}-wave NEQ-EI arise exclusively from the electron-hole Coulomb interaction as the system is not driven by external fields., 4+2 pages, 5 figures

The journal of physical chemistry letters 11(3), 891 - 899 (2020). doi:10.1021/acs.jpclett.9b03517, Extreme-ultraviolet-induced charge migration in biorelevant molecules is a fundamental step in the complex path leading to photodamage. In this work we propose a simple interpretation of the charge migration recently observed in an attosecond pump–probe experiment on the amino acid tryptophan. We find that the decay of the prominent low-frequency spectral structure with increasing pump–probe delay is due to a quantum beating between two geometrically distinct, almost degenerate charge oscillations. Quantum beating is ubiquitous in these systems, and at least on the few-to-tens of femtosecond time scales, it may dominate over decoherence the line intensities of time-resolved spectra. We also address the experimentally observed phase shift in the charge oscillations of two different amino acids, tryptophan and phenylalanine. Our results indicate that a beyond mean-field treatment of the electron dynamics is necessary to reproduce the correct behavior., Published by ACS, Washington, DC

The nonequilibrium (NEQ) excitonic-insulator (EI) phase is a dynamical state of matter for excited insulators or semiconductors, and it is characterized by self-sustained coherent oscillations of the excitonic condensate. In this Rapid Communication we highlight the main qualitative features of a NEQ-EI time- and angle-resolved photoemission spectroscopy (tr-ARPES) spectrum. We show that monochromatic probes yield a steady-state spectrum with excitonic structures originating from the dressing of conduction states with the coherent condensate. These structures are replicas of the NEQ valence bands but shifted upward in energy by the condensate frequency. Reducing the probe duration below the condensate period, the band structure gradually fades away and the tr-ARPES signal becomes proportional to the excitonic wave function. In addition, the signal amplitude becomes periodic in the impinging time of the probe, with the same period of the oscillating condensate.

We present a computationally efficient method to obtain the spectral function of bulk systems in the framework of steady-state density functional theory (i-DFT) using an idealized Scanning Tunneling Microscope (STM) setup. We calculate the current through the STM tip and then extract the spectral function from the finite-bias differential conductance. The fictitious non-interacting system of i-DFT features an exchange-correlation (xc) contribution to the bias which guarantees the same current as in the true interacting system. Exact properties of the xc bias are established using Fermi-liquid theory and subsequently implemented to construct approximations for the Hubbard model. We show for two different lattice structures that the metal-insulator transition is captured by i-DFT., Comment: 5 pages, 3 figures + supplemental material (5 pages, 3 figures); replaced with revised manuscript

Density functional theory can be generalized to mixtures of ground and excited states, for the purpose of determining energies of excitations using low-cost density functional approximations. Adapting approximations originally developed for ground states to work in the new setting would fast-forward progress enormously. But, previous attempts have stumbled on daunting fundamental issues. Here we show that these issues can be prevented from the outset, by using a fluctuation dissipation theorem (FDT) to dictate key functionals. We thereby show that existing exchange energy approximations are readily adapted to excited states, when combined with a rigorous exact Hartree term that is different in form from its ground state counterpart, and counterparts based on ensemble Ansatzë. Applying the FDT to correlation energies also provides insights into ground statelike and ensemble-only correlations. We thus provide a comprehensive and versatile framework for ensemble density functional approximations.

We study the screened dynamics of the nonequilibrium excitonic consensate forming in a bulk WSe$_{2}$ when illuminated by coherent light resonant with the lowest-energy exciton. Intervalley scattering causes electron migration from the optically populated K valley to the conduction band minimum at $\Sigma$. Due to the electron-hole unbalance at the K point a plasma of quasi-free holes develops, which efficiently screens the interaction of the remaining excitons. We show that this plasma screening causes an ultrafast melting of the nonequilibrium consensate and that during melting coherent excitons and quasi-free electron-hole pairs coexist. The time-resolved spectral function does exhibit a conduction and excitonic sidebands of opposite convexity and relative spectral weight that changes in time. Both the dependence of the time-dependent conduction density on the laser intensity and the time-resolved spectral function agree with recent experiments., Comment: 10 pages, 9 figures

Time-resolved soft-X-ray photoemission spectroscopy is used to simultaneously measure the ultrafast dynamics of core-level spectral functions and excited states upon excitation of excitons in WSe$_2$. We present a many-body approximation for the Green's function, which excellently describes the transient core-hole spectral function. The relative dynamics of excited-state signal and core levels reveals a delayed core-hole renormalization due to screening by excited quasi-free carriers, revealing an excitonic Mott transition. These findings establish time-resolved core-level photoelectron spectroscopy as a sensitive probe of subtle electronic many-body interactions and an ultrafast electronic phase transition., 6 pages, 4 figures

Diagrammatic perturbation theory is a powerful tool for the investigation of interacting many-body systems, the self-energy operator $\mathrm{\ensuremath{\Sigma}}$ encoding all the variety of scattering processes. In the simplest scenario of correlated electrons described by the $GW$ approximation for the electron self-energy, a particle transfers a part of its energy to neutral excitations. Higher-order (in screened Coulomb interaction $W$) self-energy diagrams lead to improved electron spectral functions (SFs) by taking more complicated scattering channels into account and by adding corrections to lower order self-energy terms. However, they also may lead to unphysical negative spectral functions. The resolution of this difficulty has been demonstrated in our previous works. The main idea is to represent the self-energy operator in a Fermi golden rule form which leads to a manifestly positive definite SF and allows for a very efficient numerical algorithm. So far, the method has only been applied to the three-dimensional electron gas, which is a paradigmatic system, but a rather simple one. Here we systematically extend the method to two dimensions including realistic systems such as monolayer and bilayer graphene. We focus on one of the most important vertex function effects involving the exchange of two particles in the final state. We demonstrate that it should be evaluated with the proper screening and discuss its influence on the quasiparticle properties.

We consider a ground-state wide-gap band insulator turning into a nonequilibrium excitonic insulator (NEQ-EI) upon visiting properly selected and physically relevant highly excited states. The NEQ-EI phase, characterized by self-sustained oscillations of the complex order parameter, neatly follows from a Nonequilibrium Green's Function treatment on the Konstantinov-Perel' contour. We present the first {\em ab initio} band structure of LiF, a ground-state bulk insulator, in different NEQ-EI states and show that these states can be generated by currently available pump pulses. We highlight two general features of time-resolved ARPES spectra: (1) during the pump-driving the excitonic spectral structure undergoes a convex-to-concave shape transition and {\em concomitantly} the state of the system goes through a BEC-BCS crossover; (2) attosecond pulses shone after the pump-driving at different times $t_{\rm delay}$ generate a photocurrent which {\em oscillates} in $t_{\rm delay}$ with a pump-tunable frequency -- we show that this phenomenon is similar to the AC response of an exotic Josephson junction., Comment: 13 pages, 9 figures

Multi-terminal transport setups allow to realize more complex measurements and functionalities (e.g., transistors) of nanoscale systems than the simple two-terminal arrangement. Here the steady-state density functional formalism (i-DFT) for the description of transport through nanoscale junctions with an arbitrary number of leads is developed. In a three-terminal setup and in the ideal STM limit where one of the electrodes (the `STM tip') is effectively decoupled from the junction, the formalism allows to extract its non-equilibrium spectral function (at arbitrary temperature) while a bias is applied between the other two electrodes. Multi-terminal i-DFT is shown to be capable of describing the splitting of the Kondo resonance in an Anderson impurity in the presence of an applied bias voltage, as predicted by numerically exact many-body approaches., 7 pages, 4 figures

We study time-resolved charge transport in a superconducting nanowire using time-dependent Landauer-B{\"u}ttiker theory. We find that the steady-state Majorana zero-bias conductance peak emerges transiently accompanied by characteristic oscillations after a bias-voltage quench. These oscillations are absent for a trivial impurity state that otherwise shows a very similar steady-state signal as the Majorana zero mode. In addition, we find that Andreev bound states or quasi-Majorana states in the topologically trivial bulk phase can give rise to a zero-bias conductance peak, also retaining the transient properties of the Majorana zero mode. Our results imply that (1) time-resolved transport may be used as a probe to distinguish between the Majorana zero mode and interfacial impurity states; and (2) the quasi-Majorana states mimic the transient signatures of the topological Majorana zero modes, indicating of the possibility of utilizing also the quasi-Majorana states in topological quantum computing as anyons with non-Abelian braiding statistics., Comment: 7 pages and 4 figures in the main text, 9 pages and 1 figure in the supplemental material

The Auger decay is a relevant recombination channel during the first few femtoseconds of molecular targets impinged by attosecond XUV or soft X-ray pulses. Including this mechanism in time--dependent simulations of charge--migration processes is a difficult task, and Auger scatterings are often ignored altogether. In this work we present an advance of the current state-of-the-art by putting forward a real--time approach based on nonequilibrium Green's functions suitable for first-principles calculations of molecules with tens of active electrons. To demonstrate the accuracy of the method we report comparisons against accurate grid simulations of one-dimensional systems. We also predict a highly asymmetric profile of the Auger wavepacket, with a long tail exhibiting ripples temporally spaced by the inverse of the Auger energy., Comment: 11 pages, 7 figures

We have recently put forward a steady-state density functional theory (i-DFT) to calculate the transport coefficients of quantum junctions. Within i-DFT it is possible to obtain the steady density on and the steady current through an interacting junction using a fictitious noninteracting junction subject to an effective gate and bias potential. In this work we extend i-DFT to the time domain for the single-impurity Anderson model. By a reverse engineering procedure we extract the exchange-correlation (xc) potential and xc bias at temperatures above the Kondo temperature $T_{\rm K}$. The derivation is based on a generalization of a recent paper by Dittmann et al. [arXiv:1706.04547]. Interestingly the time-dependent (TD) i-DFT potentials depend on the system's history only through the first time-derivative of the density. We perform numerical simulations of the early transient current and investigate the role of the history dependence. We also empirically extend the history-dependent TD i-DFT potentials to temperatures below $T_{\rm K}$. For this purpose we use a recently proposed parametrization of the i-DFT potentials which yields highly accurate results in the steady state., 7 pages, 4 figures

We explore the finite bias DC differential conductance of a correlated quantum dot under the influence of an AC field, from the low-temperature Kondo to the finite temperature Coulomb blockade regime. Real-time simulations are performed using a time-dependent generalization of the steady-state density functional theory (i-DFT) [Nano Lett. {\bf 15}, 8020 (2015)]. The numerical simplicity of i-DFT allows for unprecedented long time evolutions. Accurate values of average current and density are obtained by integrating over several periods of the AC field. We find that (i) the zero-temperature Kondo plateau is suppressed, (ii) the photon-assisted conductance peaks are shifted due to correlations and (iii) the Coulomb blockade is lifted with a concomitant smoothening of the sharp diamond edges., 5 pages, 4 figures

The early stage density oscillations of the electronic charge in molecules irradiated by an attosecond XUV pulse takes place on femto- or subfemtosecond timescales. This ultrafast charge migration process is a central topic in attoscience as it dictates the relaxation pathways of the molecular structure. A predictive quantum theory of ultrafast charge migration should incorporate the atomistic details of the molecule, electronic correlations and the multitude of ionization channels activated by the broad-bandwidth XUV pulse. In this work we propose a first-principles Non Equilibrium Green's Function method fulfilling all three requirements, and apply it to a recent experiment on the photoexcited phenylalanine aminoacid. Our results show that dynamical correlations are necessary for a quantitative overall agreement with the experimental data. In particular, we are able to capture the transient oscillations at frequencies 0.15PHz and 0.30PHz in the hole density of the amine group, as well as their suppression and the concomitant development of a new oscillation at frequency 0.25PHz after about 14 femtoseconds., Comment: 7 pages, 5 figures

The use of ultra-short laser pulses to pump and probe materials activates a wealth of processes which involve the coherent and non coherent dynamics of interacting electrons out of equilibrium. Non equilibrium (NEQ) many body perturbation theory (MBPT) offers an equation of motion for the density- matrix of the system which well describes both coherent and non coherent processes. In the non correlated case there is a clear relation between these two regimes and the matrix elements of the density-matrix. The same is not true for the correlated case, where the potential binding of electrons and holes in excitonic states need to be considered. In the present work we discuss how NEQ-MBPT can be used to describe the dynamics of both coherent and non-coherent excitons in the low density regime. The approach presented is well suited for an ab initio implementation.

We describe an algorithm to efficiently compute the second-Born self-energy of many-body perurbation theory. The core idea consists in dissecting the set of all four-index Coulomb integrals into properly chosen subsets, thus avoiding to loop over those indices for which the Coulomb integrals are zero or negligible. The scaling properties of the algorithm with the number of basis functions is discussed. The computational gain is demonstrated in the case of one-particle Kohn-Sham basis for organic molecules., Comment: 6 pages, contribution to the proceedings of the workshop "Progress in Nonequilibrium Green's Function VII"

Within the non-equilibrium Green's function (NEGF) formalism, the Generalized Kadanoff-Baym Ansatz (GKBA) has stood out as a computationally cheap method to investigate the dynamics of interacting quantum systems driven out of equilibrium. Current implementations of the NEGF--GKBA, however, suffer from a drawback: real-time simulations require {\em noncorrelated} states as initial states. Consequently, initial correlations must be built up through an adiabatic switching of the interaction before turning on any external field, a procedure that can be numerically highly expensive. In this work, we extend the NEGF--GKBA to allow for {\em correlated} states as initial states. Our scheme makes it possible to efficiently separate the calculation of the initial state from the real-time simulation, thus paving the way for enlarging the class of systems and external drivings accessible by the already successful NEGF--GKBA. We demonstrate the accuracy of the method and its improved performance in a model donor-acceptor dyad driven out of equilibrium by an external laser pulse.

We use the nonequilibrium Green function (NEGF) method to perform real-time simulations of the ultrafast electron dynamics of photoexcited donor-C

We present recent advances in Density Functional Theory (DFT) for applications to the field of quantum transport, with particular emphasis on transport through strongly correlated systems. We review the foundations of the popular Landauer-B\"uttiker(LB)+DFT approach. This formalism, when using approximations to the exchange-correlation (xc) potential with steps at integer occupation, correctly captures the Kondo plateau in the zero bias conductance at zero temperature but completely fails to capture the transition to the Coulomb blockade (CB) regime as temperature increases. To overcome the limitations of LB+DFT the quantum transport problem is treated from a time-dependent (TD) perspective using TDDFT, an exact framework to deal with nonequilibrium situations. The steady-state limit of TDDFT shows that in addition to an xc potential in the junction, there also exists an xc correction to the applied bias. Open shell molecules in the CB regime provide the most striking examples of the importance of the xc bias correction. Using the Anderson model as guidance we estimate these corrections in the limit of zero bias. For the general case we put forward a steady-state DFT which is based on the one-to-one correspondence between the pair of basic variables steady density on and steady current across the junction and the pair local potential on and bias across the junction. Like TDDFT, this framework also leads to both an xc potential in the junction and an xc correction to the bias. Unlike in TDDFT, these potentials are independent of history. We highlight the universal features of both xc potential and xc bias corrections for junctions in the CB regime and provide an accurate parametrization for the Anderson model at arbitrary temperatures and interaction strengths thus providing a unified DFT description for both Kondo and CB regimes and the transition between them., Comment: 29 pages, 22 Figures

At zero temperature, the Landauer formalism combined with static density functional theory is able to correctly reproduce the Kondo plateau in the conductance of the Anderson impurity model provided that an exchange-correlation potential is used which correctly exhibits steps at integer occupation. Here we extend this recent finding to multi-level quantum dots described by the constant-interaction model. We derive the exact exchange-correlation potential in this model for the isolated dot and deduce an accurate approximation for the case when the dot is weakly coupled to two leads. We show that at zero temperature and for non-degenerate levels in the dot we correctly obtain the conductance plateau for any odd number of electrons on the dot. We also analyze the case when some of the levels of the dot are degenerate and again obtain good qualitative agreement with results obtained with alternative methods. As in the case of a single level, for temperatures larger than the Kondo temperature, the Kohn-Sham conductance fails to reproduce the typical Coulomb blockade peaks. This is attributed to {\em dynamical} exchange-correlation corrections to the conductance originating from time-dependent density functional theory.

We present a systematic study of vertex corrections in the homogeneous electron gas at metallic densities. The vertex diagrams are built using a recently proposed positive-definite diagrammatic expansion for the spectral function. The vertex function not only provides corrections to the well known plasmon and particle-hole scatterings, but also gives rise to new physical processes such as generation of two plasmon excitations or the decay of the one-particle state into a two-particles-one-hole state. By an efficient Monte Carlo momentum integration we are able to show that the additional scattering channels are responsible for the bandwidth reduction observed in photoemission experiments on bulk sodium, appearance of the secondary plasmon satellite below the Fermi level, and a substantial redistribution of spectral weights. The feasibility of the approach for first-principles band-structure calculations is also discussed.

The single-impurity Anderson model is studied within the i-DFT framework, a recently proposed extension of density functional theory (DFT) for the description of electron transport in the steady state. i-DFT is designed to give both the steady current and density at the impurity, and it requires the knowledge of the exchange-correlation (xc) bias and on-site potential (gate). In this work we construct an approximation for both quantities which is accurate in a wide range of temperatures, gates and biases, thus providing a simple and unifying framework to calculate the differential conductance at negligible computational cost in different regimes. Our results mark a substantial advance for DFT and may inform the construction of functionals applicable to other correlated systems., 5 pages, 4 figures