-
Anharmonic lattice dynamics and superconductivity in strained bulk and surface niobium
Authors:
Mihir Ranjan Sahoo,
Roman Lucrezi,
Pedro Nunes Ferreira,
Chia-Nien Tsai,
Matthew Julian,
Rohit P. Prasankumar,
Mahmoud I. Hussein,
Christoph Heil
Abstract:
Using first-principles calculations, we investigate how homogeneous strain and crystallographic surface orientation modify the vibrational and superconducting properties of niobium. For bulk Nb, tensile strain strongly softens the phonon spectrum and enhances the electron--phonon coupling, increasing the superconducting transition temperature from 9.5 K at equilibrium to 14.5 K at $\sim\!6\%$ latt…
▽ More
Using first-principles calculations, we investigate how homogeneous strain and crystallographic surface orientation modify the vibrational and superconducting properties of niobium. For bulk Nb, tensile strain strongly softens the phonon spectrum and enhances the electron--phonon coupling, increasing the superconducting transition temperature from 9.5 K at equilibrium to 14.5 K at $\sim\!6\%$ lattice expansion. For the low-index Nb(001), Nb(110), and Nb(111) surfaces, harmonic phonon calculations exhibit imaginary modes, showing that anharmonic lattice effects are essential. To treat these effects efficiently, we train Nb-specific machine-learning interatomic potentials on bulk and slab first-principles configurations and use them to accelerate stochastic self-consistent harmonic approximation calculations, thereby obtaining anharmonically renormalized phonon modes that are combined with density-functional perturbation theory electron--phonon matrix elements to construct the Eliashberg spectral function. Among the clean free-standing slabs considered here, Nb(001) exhibits the strongest electron--phonon coupling and the highest calculated transition temperature of 10.0 K, while Nb(110) and Nb(111) show progressively reduced pairing strength. Finally, by analyzing the Eliashberg spectral function and the functional derivative $δT_\text{c}/δα^2F(ω)$, we identify the phonon energy ranges most effective for superconducting pairing. Our results show that strain, surface termination, and anharmonic phonon renormalization provide complementary and interrelated microscopic routes for tuning superconductivity in Nb.
△ Less
Submitted 1 June, 2026;
originally announced June 2026.
-
Superhydrides on the way to ambient pressure: weak localization and persistent X-ray photoconductivity in BaSiH$_{8}$
Authors:
Dmitrii V. Semenok,
Di Zhou,
Sven Luther,
Toni Helm,
Hirokazu Kadobayashi,
Yuki Nakamoto,
Katsuya Shimizu,
Kirill S. Pervakov,
Andrei V. Sadakov,
Oleg A. Sobolevskiy,
Vladimir M. Pudalov,
Simone Di Cataldo,
Roman Lucrezi,
Lilia Boeri,
Michele Galasso,
Frederico G. Alabarse,
Ivan A. Troyan,
Viktor V. Struzhkin
Abstract:
Reducing the stabilization pressure of superhydrides represents one of the most important challenges in hydrogen-saturated compound chemistry. Moving in this direction, we studied the Ba-Si-H system at 0-142 GPa using transport measurements, 1H nuclear magnetic resonance, single-crystal and powder X-ray diffraction in the temperature range of 4-317 K. We synthesized the previously predicted cubic…
▽ More
Reducing the stabilization pressure of superhydrides represents one of the most important challenges in hydrogen-saturated compound chemistry. Moving in this direction, we studied the Ba-Si-H system at 0-142 GPa using transport measurements, 1H nuclear magnetic resonance, single-crystal and powder X-ray diffraction in the temperature range of 4-317 K. We synthesized the previously predicted cubic BaSiH$_{8}$ at pressures of 18-31 GPa. Remarkably, we demonstrate that BaSiH$_8$ remains stable upon decompression to ambient conditions and can be recovered from the diamond anvil cell. Obtained Ba-Si polyhydrides exhibit metallic and superconducting properties ($\textit{T$_c$}$ = 9 K, $\textit{B$_{c2}$}$(0)=13-16 T) at 142 GPa. However, at pressures below 50 GPa, these hydrides behave as degenerate semiconductors (bandgap < 0.4 meV) or poor metals with weak electron localization, negative magnetoresistance, photovoltaic effect, and persistent photoconductivity in the X-ray and visible range. Our work demonstrates the high-pressure synthesis of Ba-Si polyhydrides that remain stable upon decompression to ambient conditions, overcoming a critical bottleneck in superhydride chemistry and establishing a foundation for practical applications in hydrogen storage.
△ Less
Submitted 14 March, 2026;
originally announced March 2026.
-
Inverse Isotope Effect in the Ternary Perovskite Hydride SrPdH/D$_{2.9}$: A Signature of Quantum Zero-Point Fluctuations
Authors:
Wencheng Lu,
Mihir Sahoo,
Roman Lucrezi,
Michael J. Hutcheon,
Shubham Sinha,
Pedro N. Ferreira,
Chris J. Pickard,
Qiang Zhang,
Matthew N. Julian,
Rohit P. Prasankumar,
Christoph Heil,
Timothy A. Strobel
Abstract:
Guided by first-principles calculations, we demonstrate superconductivity in the ternary perovskite hydride SrPdH$_{3-x}$, synthesized at low pressure. Structural characterization via neutron diffraction reveals the near-stoichiometric composition SrPdD$_{2.9(2)}$ with 96\% deuterium site occupancy. Subsequent transport and magnetic susceptibility measurements establish onset superconducting trans…
▽ More
Guided by first-principles calculations, we demonstrate superconductivity in the ternary perovskite hydride SrPdH$_{3-x}$, synthesized at low pressure. Structural characterization via neutron diffraction reveals the near-stoichiometric composition SrPdD$_{2.9(2)}$ with 96\% deuterium site occupancy. Subsequent transport and magnetic susceptibility measurements establish onset superconducting transitions at $T_\text{c} = \SI{2.1}{K} $ (H) and $T_\text{c} = \SI{2.2}{K} $ (D), exhibiting an inverse isotope effect that our first-principles calculations attribute predominantly to quantum zero-point motion. The excellent agreement between theory and experiment with respect to thermodynamic stability and superconducting properties provides important validation for theory-guided superconductor discovery. This work establishes superconductivity in the perovskite hydride structural prototype -- expanding the limited family of experimentally realized ternary hydride superconductors -- and demonstrates the importance of quantum nuclear motion on the accurate theoretical treatment of low-pressure hydride superconductors.
△ Less
Submitted 27 February, 2026;
originally announced February 2026.
-
First-principles evidence for conventional superconductivity in a quasicrystal approximant
Authors:
Pedro N. Ferreira,
Roman Lucrezi,
Sangmin Lee,
Lucy Nathwani,
Matthew Julian,
Rohit P. Prasankumar,
Warren E. Pickett,
Chris J. Pickard,
Philip Kim,
Christoph Heil
Abstract:
Quasicrystals (QCs) host long-range order without translational symmetry, a regime in which the very foundations of BCS theory are not straightforwardly applicable, yet experiments on QCs and their approximant crystals (ACs) point to conventional, $s$-wave, electron-phonon coupled superconductivity. Here we test the predictive power of the electron-phonon framework in a representative decagonal AC…
▽ More
Quasicrystals (QCs) host long-range order without translational symmetry, a regime in which the very foundations of BCS theory are not straightforwardly applicable, yet experiments on QCs and their approximant crystals (ACs) point to conventional, $s$-wave, electron-phonon coupled superconductivity. Here we test the predictive power of the electron-phonon framework in a representative decagonal AC from first principles. Using state-of-the-art \textit{ab initio} methods, we compute the superconducting properties of the recently discovered AC Al$_{13}$Os$_4$ and quantitatively reproduce its bulk $T_\text{c}$. This constitutes, to our knowledge, the first \textit{ab initio} determination of $T_\text{c}$ for an AC and establishes that the electron-phonon framework is predictive in these systems as well. Using the generalized quasichemical approximation for alloy modeling in the decagonal Al--Os family, we predict tunable superconductivity in Al$_{13}$Os$_{4-x}$Re$_x$ and Al$_{13}$Os$_{4-x}$Ir$_x$; in particular, Al$_{13}$Re$_4$ is dynamically stable and estimated to have a $T_\text{c}$ about 30% above Al$_{13}$Os$_4$. Finally, we discuss the role of ACs as high-fidelity proxies for their parent QCs. Although long-range quasiperiodicity may introduce subtle electronic features, our findings indicate that the key ingredients for superconductivity are already encoded in the local structural motifs preserved by the AC. This places the Al--Os and Al--Re families among the most promising candidates for the highest-$T_\text{c}$ quasicrystalline superconductivity.
△ Less
Submitted 1 April, 2026; v1 submitted 12 November, 2025;
originally announced November 2025.
-
Vacancy-free cubic superconducting NbN enabled by quantum anharmonicity
Authors:
Eva Kogler,
Mihir R. Sahoo,
Chia-Nien Tsai,
Fabian Jöbstl,
Roman Lucrezi,
Peter I. C. Cooke,
Birgit Kunert,
Roland Resel,
Chris J. Pickard,
Matthew N. Julian,
Rohit P. Prasankumar,
Mahmoud I. Hussein,
Christoph Heil
Abstract:
Niobium nitride (NbN) is renowned for its exceptional mechanical, electronic, magnetic, and superconducting properties. The ideal 1:1 stoichiometric $δ$-NbN cubic phase, however, is known to be dynamically unstable, and repeated experimental observations have indicated that vacancies are necessary for its stabilization. In this work, we demonstrate that when the structure is fully relaxed and allo…
▽ More
Niobium nitride (NbN) is renowned for its exceptional mechanical, electronic, magnetic, and superconducting properties. The ideal 1:1 stoichiometric $δ$-NbN cubic phase, however, is known to be dynamically unstable, and repeated experimental observations have indicated that vacancies are necessary for its stabilization. In this work, we demonstrate that when the structure is fully relaxed and allowed to distort under quantum anharmonic effects, a previously unreported stable cubic phase with space group $P\bar{4}3m$ emerges - 65 meV/atom lower in free energy than the ideal $δ$ phase. This discovery is enabled by state-of-the-art first-principles calculations accelerated by machine-learned interatomic potentials. To evaluate the vibrational and superconducting properties with quantum anharmonic effects accounted for, we use the stochastic self-consistent harmonic approximation (SSCHA) and molecular dynamics spectral energy density (SED) methods. Electron-phonon coupling calculations based on the SSCHA phonon dispersion yield a superconducting transition temperature of $T_\text{c}$ = 20 K, which aligns closely with experimentally reported values for near-stoichiometric NbN. These findings challenge the long-held assumption that vacancies are essential for stabilizing cubic NbN and point to the potential of synthesizing the ideal 1:1 stoichiometric phase as a route to achieving enhanced superconducting performance in this technologically significant material.
△ Less
Submitted 4 July, 2025;
originally announced July 2025.
-
IsoME: Streamlining High-Precision Eliashberg Calculations
Authors:
Eva Kogler,
Dominik Spath,
Roman Lucrezi,
Hitoshi Mori,
Zien Zhu,
Zhenglu Li,
Elena R. Margine,
Christoph Heil
Abstract:
This paper introduces the Julia package IsoME, an easy-to-use yet accurate and robust computational tool designed to calculate superconducting properties. Multiple levels of approximation are supported, ranging from the basic McMillan-Allen-Dynes formula and its machine learning-enhanced variant to Eliashberg theory including static Coulomb interactions derived from $GW$ calculations, offering a f…
▽ More
This paper introduces the Julia package IsoME, an easy-to-use yet accurate and robust computational tool designed to calculate superconducting properties. Multiple levels of approximation are supported, ranging from the basic McMillan-Allen-Dynes formula and its machine learning-enhanced variant to Eliashberg theory including static Coulomb interactions derived from $GW$ calculations, offering a fully ab initio approach to determine superconducting properties, such as the critical superconducting temperature ($T_\text{c}$) and the superconducting gap function ($Δ$). We validate IsoME by benchmarking it against various materials, demonstrating its versatility and performance across different theoretical levels. The findings indicate that the previously held assumption that Eliashberg theory overestimates $T_\text{c}$ is no longer valid when $μ^*$ is appropriately adjusted to account for the finite Matsubara frequency cutoff. Furthermore, we conclude that the constant density of states (DOS) approximation remains accurate in most cases. By unifying multiple approximation schemes within a single framework, IsoME combines first-principles precision with computational efficiency, enabling seamless integration into high-throughput workflows through its $T_\text{c}$ search mode. This makes IsoME a powerful and reliable tool for advancing superconductivity research.
△ Less
Submitted 18 June, 2025; v1 submitted 5 March, 2025;
originally announced March 2025.
-
Nevanlinna Analytic Continuation for Migdal-Eliashberg Theory
Authors:
D. M. Khodachenko,
R. Lucrezi,
P. N. Ferreira,
M. Aichhorn,
C. Heil
Abstract:
In this work, we present a method to reconstruct real-frequency properties from analytically continued causal Green's functions within the framework of Migdal-Eliashberg (ME) theory for superconductivity. ME theory involves solving a set of coupled equations self-consistently in imaginary frequency space, but to obtain experimentally measurable properties like the spectral function and quasipartic…
▽ More
In this work, we present a method to reconstruct real-frequency properties from analytically continued causal Green's functions within the framework of Migdal-Eliashberg (ME) theory for superconductivity. ME theory involves solving a set of coupled equations self-consistently in imaginary frequency space, but to obtain experimentally measurable properties like the spectral function and quasiparticle density of states, it is necessary to perform an analytic continuation to real frequency space. Traditionally, the ME Green's function is decomposed into three fundamental complex functions, which are analytically continued independently. However, these functions do not possess the causal properties of Green's functions, complicating or even preventing the application of standard methods such as Maximum Entropy. Our approach overcomes these challenges, enabling the use of various analytic continuation techniques that were previously impractical. We demonstrate the effectiveness of this method by combining it with Nevanlinna analytic continuation to achieve accurate real-frequency results for ME theory, which are directly comparable to experimental data, with applications highlighted for the superconductors MgB$_2$ and LaBeH$_8$.
△ Less
Submitted 14 October, 2024; v1 submitted 4 September, 2024;
originally announced September 2024.
-
Ab initio modeling of superconducting alloys
Authors:
Pedro N. Ferreira,
Roman Lucrezi,
Ivan Guilhon,
Marcelo Marques,
Lara K. Teles,
Christoph Heil,
Luiz T. F. Eleno
Abstract:
Designing new, technologically relevant superconductors has long been at the forefront of solid-state physics and chemistry research. However, developing efficient approaches for modeling the thermodynamics of superconducting alloys while accurately evaluating their physical properties has proven to be a very challenging task. To fill this gap, we propose an ab initio thermodynamic statistical met…
▽ More
Designing new, technologically relevant superconductors has long been at the forefront of solid-state physics and chemistry research. However, developing efficient approaches for modeling the thermodynamics of superconducting alloys while accurately evaluating their physical properties has proven to be a very challenging task. To fill this gap, we propose an ab initio thermodynamic statistical method, the Extended Generalized Quasichemical Approximation (EGQCA), to describe off-stoichiometric superconductors. Within EGQCA, one can predict any computationally accessible property of the alloy, such as the critical temperature in superconductors and the electron-phonon coupling parameter, as a function of composition and crystal growth conditions by computing the cluster occurrence probabilities that minimize the overall mixing Gibbs free energy. Importantly, EGQCA incorporates directly chemical ordering, lattice distortions, and vibrational contributions. As a proof of concept, we applied EGQCA to the well-known Al-doped MgB$_2$ and to niobium alloyed with titanium and vanadium, showing a remarkable agreement with the experimental data. Additionally, we model the near-room temperature sodalite-like Y$_{1-x}$Ca$_x$H$_6$ superconducting solid solution, demonstrating that EGQCA particularly possesses a promising potential for designing in silico high-$T_{\text{c}}$ superhydride alloys. Our approach notably enables the high-throughput screening of complex superconducting solid solutions, intrinsically providing valuable insights into the interplay between synthesis, thermodynamics, and physical properties.
△ Less
Submitted 26 June, 2024; v1 submitted 21 June, 2024;
originally announced June 2024.
-
Full-bandwidth anisotropic Migdal-Eliashberg theory and its application to superhydrides
Authors:
Roman Lucrezi,
Pedro P. Ferreira,
Samad Hajinazar,
Hitoshi Mori,
Hari Paudyal,
Elena R. Margine,
Christoph Heil
Abstract:
Migdal-Eliashberg theory is one of the state-of-the-art methods for describing conventional superconductors from first principles. However, widely used implementations assume a constant density of states around the Fermi level, which hinders a proper description of materials with distinct features in its vicinity. Here, we present an implementation of the Migdal-Eliashberg theory within the EPW co…
▽ More
Migdal-Eliashberg theory is one of the state-of-the-art methods for describing conventional superconductors from first principles. However, widely used implementations assume a constant density of states around the Fermi level, which hinders a proper description of materials with distinct features in its vicinity. Here, we present an implementation of the Migdal-Eliashberg theory within the EPW code that considers the full electronic structure and accommodates scattering processes beyond the Fermi surface. To significantly reduce computational costs, we introduce a non-uniform sampling scheme along the imaginary axis. We demonstrate the power of our implementation by applying it to the sodalite-like clathrates YH$_6$ and CaH$_6$, and to the covalently-bonded H$_3$S and D$_3$S. Furthermore, we investigate the effect of maximizing the density of states at the Fermi level in doped H$_3$S and BaSiH$_8$ within the full-bandwidth treatment compared to the constant-density-of-states approximation. Our findings highlight the importance of this advanced treatment in such complex materials.
△ Less
Submitted 16 January, 2024; v1 submitted 29 September, 2023;
originally announced October 2023.
-
Temperature and quantum anharmonic lattice effects on stability and superconductivity in lutetium trihydride
Authors:
Roman Lucrezi,
Pedro P. Ferreira,
Markus Aichhorn,
Christoph Heil
Abstract:
In this work, we resolve conflicting experimental and theoretical findings related to the dynamical stability and superconducting properties of $Fm\overline{3}m$-LuH$_3$, which was recently suggested as the parent phase harboring room-temperature superconductivity at near-ambient pressures. Including temperature and quantum anharmonic lattice effects in our calculations, we demonstrate that the th…
▽ More
In this work, we resolve conflicting experimental and theoretical findings related to the dynamical stability and superconducting properties of $Fm\overline{3}m$-LuH$_3$, which was recently suggested as the parent phase harboring room-temperature superconductivity at near-ambient pressures. Including temperature and quantum anharmonic lattice effects in our calculations, we demonstrate that the theoretically predicted structural instability of the $Fm\overline{3}m$ phase near ambient pressures is suppressed for temperatures above $200\,\text{K}$. We provide a $p\,\unicode{x2013}\,T$ phase diagram for stability up to pressures of $6\,\text{GPa}$, where the required temperature for stability is reduced to $T>80\,\text{K}$. We also determine the superconducting critical temperature $T_\text{c}$ of $Fm\overline{3}m$-LuH$_3$ within the Migdal-Eliashberg formalism, using temperature- and quantum-anharmonically-corrected phonon dispersions, finding that the expected $T_\text{c}$ for electron-phonon mediated superconductivity is in the range of $50$ $\unicode{x2013}$ $60\,\text{K}$, i.e., well below the temperatures required to stabilize the lattice. When considering moderate doping based on rigidly shifting the Fermi level, $T_\text{c}$ decreases for both hole and electron doping. Our results thus provide evidence that any observed room-temperature superconductivity in pure or doped $Fm\overline{3}m$-LuH$_3$, if confirmed, cannot be explained by a conventional electron-phonon mediated pairing mechanism.
△ Less
Submitted 11 January, 2024; v1 submitted 13 April, 2023;
originally announced April 2023.
-
Quantum lattice dynamics and their importance in ternary superhydride clathrates
Authors:
Roman Lucrezi,
Eva Kogler,
Simone Di Cataldo,
Markus Aichhorn,
Lilia Boeri,
Christoph Heil
Abstract:
The quantum nature of the hydrogen lattice in superconducting hydrides can have crucial effects on the material's properties. Taking a detailed look at the dynamic stability of the recently predicted BaSiH$_8$ phase, we find that the inclusion of anharmonic quantum ionic effects leads to an increase in the critical dynamical pressure to $20\,\text{GPa}$ as compared to $5\,\text{GPa}$ within the ha…
▽ More
The quantum nature of the hydrogen lattice in superconducting hydrides can have crucial effects on the material's properties. Taking a detailed look at the dynamic stability of the recently predicted BaSiH$_8$ phase, we find that the inclusion of anharmonic quantum ionic effects leads to an increase in the critical dynamical pressure to $20\,\text{GPa}$ as compared to $5\,\text{GPa}$ within the harmonic approximation. We identify the change in the crystal structure due to quantum ionic effects to be the main driving force for this increase and demonstrate that this can already be understood at the harmonic level by considering zero-point energy corrections to the total electronic energy. In fact, the previously determined critical pressure of kinetic stability $p_\text{kin} = 30\,\text{GPa}$ still poses a stricter bound for the synthesizability of BaSiH$_8$ and similar hydride materials than the dynamical stability and therefore constitutes a more rigorous and accurate estimate for the experimental realizability of these structures.
△ Less
Submitted 16 October, 2023; v1 submitted 19 December, 2022;
originally announced December 2022.
-
In-silico synthesis of lowest-pressure high-$T_\text{c}$ ternary superhydrides
Authors:
Roman Lucrezi,
Simone Di Cataldo,
Wolfgang von der Linden,
Lilia Boeri,
Christoph Heil
Abstract:
We report the theoretical prediction of two high-performing hydride superconductors BaSiH$_8$ and SrSiH$_8$. They are thermodynamically stable above pressures of $130$ and $174\,\text{GPa}$, respectively, and metastable below that. Employing anharmonic phonon calculations, we determine the minimum pressures of dynamical stability to be around $3\,\text{GPa}$ for BaSiH$_8$ and $27\,\text{GPa}$ for…
▽ More
We report the theoretical prediction of two high-performing hydride superconductors BaSiH$_8$ and SrSiH$_8$. They are thermodynamically stable above pressures of $130$ and $174\,\text{GPa}$, respectively, and metastable below that. Employing anharmonic phonon calculations, we determine the minimum pressures of dynamical stability to be around $3\,\text{GPa}$ for BaSiH$_8$ and $27\,\text{GPa}$ for SrSiH$_8$, and using the fully anisotropic Migdal-Eliashberg theory, we predict $T_\text{c}$'s around $71$ and $126\,\text{K}$, respectively. We also introduce a method to estimate the lowest pressure of synthesis, based on the calculation of the enthalpy barriers protecting the BaSiH$_8$ $Fm\bar{3}m$ structure from decomposition at various pressures. This $kinetic$ pressure threshold is sensibly higher than the one based on $dynamic$ stability, but gives a much more rigorous limit for synthesizability.
△ Less
Submitted 30 May, 2022; v1 submitted 3 December, 2021;
originally announced December 2021.
-
Superconductivity and strong anharmonicity in novel Nb-S phases
Authors:
Roman Lucrezi,
Christoph Heil
Abstract:
In this work we explore the phase diagram of the binary Nb-S system from ambient pressures up to 250 GPa using $ab$ $initio$ evolutionary crystal structure prediction. We find several new stable compositions and phases, especially in the high-pressure regime, and investigate their electronic, vibrational, and superconducting properties. Our calculations show that all materials, besides the low-pre…
▽ More
In this work we explore the phase diagram of the binary Nb-S system from ambient pressures up to 250 GPa using $ab$ $initio$ evolutionary crystal structure prediction. We find several new stable compositions and phases, especially in the high-pressure regime, and investigate their electronic, vibrational, and superconducting properties. Our calculations show that all materials, besides the low-pressure phases of pure sulfur, are metals with low electron-phonon coupling strengths and critical superconducting temperatures below 15 K. Furthermore, we investigate the effects of phonon anharmonicity on lattice dynamics, electron-phonon interactions, and superconductivity for the novel high-pressure phase of $\mathrm{Nb}_2\mathrm{S}$, demonstrating that the inclusion of anharmonicity stabilizes the lattice and enhances the electron-phonon interaction.
△ Less
Submitted 12 January, 2021;
originally announced January 2021.
-
Non-equilibrium Green's functions and their relation to the negative differential conductance in the interacting resonant level model
Authors:
Max E. Sorantin,
Wolfgang von der Linden,
Roman Lucrezi,
Enrico Arrigoni
Abstract:
We evaluate the non-equilibrium single particle Green's functions in the steady state of the interacting resonant level model (IRLM) under the effect of an applied bias voltage. Employing the so-called auxiliary master equation approach, we present accurate nonperturbative results for the non-equilibrium spectral and effective distribution functions, as well as for the current-voltage characterist…
▽ More
We evaluate the non-equilibrium single particle Green's functions in the steady state of the interacting resonant level model (IRLM) under the effect of an applied bias voltage. Employing the so-called auxiliary master equation approach, we present accurate nonperturbative results for the non-equilibrium spectral and effective distribution functions, as well as for the current-voltage characteristics. We find a drastic change of these spectral properties between the regimes of low and high bias voltages and discuss the relation of these changes to the negative differential conductance (NDC), a prominent feature in the non-equilibrium IRLM. The anomalous evolution of the distribution function next to the impurity shown by our calculations suggests a mechanism whereby the impurity gets effectively decoupled from the leads at voltages where the NDC sets in, in agreement with previous renormalization group approaches. This scenario is qualitatively confirmed by a Hartree-Fock treatment of the model.
△ Less
Submitted 23 January, 2019; v1 submitted 20 December, 2018;
originally announced December 2018.