Showing 1–36 of 36 results for author: Heil, C

Probing the superconducting gap structure of ScRuSi via $μ$SR and first-principles calculations

Authors: K. Panda, A. Bhattacharyya, P. N. Ferreira, Rajib Mondal, A. Thamizhavel, D. T. Adroja, C. Heil, L. T. F. Eleno, A. D. Hillier

Abstract: In this study, we present a thorough investigation into the superconducting state of the ruthenium-based ternary equiatomic compound ScRuSi. Our analysis combines experimental techniques, including muon spin rotation/relaxation ($μ$SR) and low-temperature resistivity measurements, with theoretical insights derived from first-principles calculations. The low-temperature resistivity measurements rev… ▽ More In this study, we present a thorough investigation into the superconducting state of the ruthenium-based ternary equiatomic compound ScRuSi. Our analysis combines experimental techniques, including muon spin rotation/relaxation ($μ$SR) and low-temperature resistivity measurements, with theoretical insights derived from first-principles calculations. The low-temperature resistivity measurements reveal a distinct superconducting phase transition in the orthorhombic structure of ScRuSi at a critical temperature ($T_\text{C}$) of $2.5$ K. Further, the TF-$μ$SR analysis yields a gap-to-critical-temperature ratio of $2Δ/k_\mathrm{B}T_\mathrm{C} = 2.71$, a value consistent with results obtained from previous heat capacity measurements. The temperature dependence of the superconducting normalized depolarization rate is fully described by the isotropic $s$-wave gap model. Additionally, zero-field $μ$SR measurements indicate that the relaxation rate remains nearly identical below and above $T_\text{C}$. This observation strongly suggests the preservation of time-reversal symmetry within the superconducting state. By employing the McMillan-Allen-Dynes equation, we calculate a $T_\text{C}$ of $2.11$ K from first-principles calculations within the density functional theory framework. This calculated value aligns closely with the experimentally determined critical temperature. The coupling between the low-frequency phonon modes and the transition metal d-orbital states play an important role in governing the superconducting pairing in ScRuSi. The combination of experimental and theoretical approaches provides a comprehensive microscopic understanding of the superconducting nature of ScRuSi, offering insights into its critical temperature, pairing symmetry, and the underlying electron-phonon coupling mechanism. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: 10 pages, 7 figures

Journal ref: Physical Review B 109, 224517 (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.

Synthesis of Mg$_2$IrH$_5$: A potential pathway to high-$T_c$ hydride superconductivity at ambient pressure

Authors: Mads F. Hansen, Lewis J. Conway, Kapildeb Dolui, Christoph Heil, Chris J. Pickard, Anna Pakhomova, Mohammed Mezouar, Martin Kunz, Rohit P. Prasankumar, Timothy A. Strobel

Abstract: Following long-standing predictions associated with hydrogen, high-temperature superconductivity has recently been observed in several hydride-based materials. Nevertheless, these high-$T_c$ phases only exist at extremely high pressures, and achieving high transition temperatures at ambient pressure remains a major challenge. Recent predictions of the complex hydride Mg$_{2}$IrH$_{6}$ may help ove… ▽ More Following long-standing predictions associated with hydrogen, high-temperature superconductivity has recently been observed in several hydride-based materials. Nevertheless, these high-$T_c$ phases only exist at extremely high pressures, and achieving high transition temperatures at ambient pressure remains a major challenge. Recent predictions of the complex hydride Mg$_{2}$IrH$_{6}$ may help overcome this challenge with calculations of high-$T_c$ superconductivity (65 K$~<~T_c~<~$ 170 K) in a material that is stable at atmospheric pressure. In this work, the synthesis of Mg$_{2}$IrH$_{6}$ was targeted over a broad range of $P$-$T$ conditions, and the resulting products were characterized using X-ray diffraction (XRD) and vibrational spectroscopy, in concert with first-principles calculations. The results indicate that the charge-balanced complex hydride Mg$_{2}$IrH$_{5}$ is more stable over all conditions tested up to ca 28 GPa. The resulting hydride is isostructural with the predicted superconducting Mg$_{2}$IrH$_{6}$ phase except for a single hydrogen vacancy, which shows a favorable replacement barrier upon insertion of hydrogen into the lattice. Bulk Mg$_{2}$IrH$_{5}$ is readily accessible at mild $P$-$T$ conditions and may thus represent a convenient platform to access superconducting Mg$_{2}$IrH$_{6}$ via non-equilibrium processing methods. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: 18 pages, 17 Figures

Machine learning assisted prediction of organic salt structure properties

Authors: Ethan P. Shapera, Dejan-Krešimir Bučar, Rohit P. Prasankumar, Christoph Heil

Abstract: We demonstrate a machine learning-based approach which predicts the properties of crystal structures following relaxation based on the unrelaxed structure. Use of crystal graph singular values reduces the number of features required to describe a crystal by more than an order of magnitude compared to the full crystal graph representation. We construct machine learning models using the crystal grap… ▽ More We demonstrate a machine learning-based approach which predicts the properties of crystal structures following relaxation based on the unrelaxed structure. Use of crystal graph singular values reduces the number of features required to describe a crystal by more than an order of magnitude compared to the full crystal graph representation. We construct machine learning models using the crystal graph singular value representations in order to predict the volume, enthalpy per atom, and metal versus semiconducting phase of DFT-relaxed organic salt crystals based on randomly generated unrelaxed crystal structures. Initial base models are trained to relate 89,949 randomly generated structures of salts formed by varying ratios of 1,3,5-triazine and HCl with the corresponding volumes, enthalpies per atom, and phase of the DFT-relaxed structures. We further demonstrate that the base model is able to extrapolate to new chemical systems with the inclusion of 2,000 to 10,000 crystal structures from the new system. After training a single model with a large number of data points, extension can be done at significantly lower cost. The constructed machine learning models can be used to rapidly screen large sets of randomly generated organic salt crystal structures and efficiently downselect the structures most likely to be experimentally realizable. The models can be used either as a stand-alone crystal structure predictor or incorporated into more sophisticated workflows as a filtering step. △ Less

Submitted 14 February, 2024; originally announced February 2024.

Comments: 18 pages, 20 figures

Feasible route to high-temperature ambient-pressure hydride superconductivity

Authors: Kapildeb Dolui, Lewis J. Conway, Christoph Heil, Timothy A. Strobel, Rohit Prasankumar, Chris J. Pickard

Abstract: A key challenge in materials discovery is to find high-temperature superconductors. Hydrogen and hydride materials have long been considered promising materials displaying conventional phonon-mediated superconductivity. However, the high pressures required to stabilize these materials have restricted their application. Here, we present results from high-throughput computation, considering a wide r… ▽ More A key challenge in materials discovery is to find high-temperature superconductors. Hydrogen and hydride materials have long been considered promising materials displaying conventional phonon-mediated superconductivity. However, the high pressures required to stabilize these materials have restricted their application. Here, we present results from high-throughput computation, considering a wide range of high-symmetry ternary hydrides from across the periodic table at ambient pressure. This large composition space is then reduced by considering thermodynamic, dynamic, and magnetic stability, before direct estimations of the superconducting critical temperature. This approach has revealed a metastable ambient-pressure hydride superconductor, Mg$_2$IrH$_6$, with a predicted critical temperature of 160 K, comparable to the highest temperature superconducting cuprates. We propose a synthesis route via a structurally related insulator, Mg$_2$IrH$_7$, which is thermodynamically stable above 15 GPa and discuss the potential challenges in doing so. △ Less

Submitted 30 April, 2024; v1 submitted 11 October, 2023; originally announced October 2023.

Comments: 30 pages, 30 figures, including SM

Journal ref: Phys. Rev. Lett. 132, 166001, 2024

Critical Role of Disorder for Superconductivity in the Series of Epitaxial Ti(O,N) Films

Authors: Fengmiao Li, Oliver Dicks, Myung-Geun Han, Solveig Aamlid, Giorgio Levy, Ronny Sutarto, Chong Liu, Hsiang-Hsi Kung, Oleksandr Foyevstov, Simon Godin, Bruce A. Davidson, Andrea Damascelli, Yimei Zhu, Christoph Heil, Ilya Elfimov, George A. Sawatzky, Ke Zou

Abstract: Experimental manipulation of superconductivity is of paramount importance, not only for practical applications but also for identifying the key factors involved in electron pairing. In this work, we have undertaken a meticulous study of the superconductivity in a series of titanium compounds with a rocksalt structure, synthesized as epitaxial films. We find that substituting nitrogen (N) for oxyge… ▽ More Experimental manipulation of superconductivity is of paramount importance, not only for practical applications but also for identifying the key factors involved in electron pairing. In this work, we have undertaken a meticulous study of the superconductivity in a series of titanium compounds with a rocksalt structure, synthesized as epitaxial films. We find that substituting nitrogen (N) for oxygen (O) in titanium monoxide (TiO) with the stoichiometry close to TiO$_{0.6}$N$_{0.4}$ leads to superconductivity with a transition temperature (T$_c$) of ~2.6 K, about five times higher than that of TiO at ~0.5 K and half as high as the T$_c$ of ~6 K in titanium nitride (TiN). However, Eliashberg theoretical calculations predict similar Tc in TiO, Ti oxynitride and TiN. The analysis of electron mean free path suggests the presence of significant disorder in TiO and a remarkable reduction in the impact of disorder in oxynitrides. Density functional theory (DFT) calculations reveal that disorder decreases the coherence of electronic states for non-zero momenta, which would degrade the influence of electron-phonon. Our findings demonstrate the disorder and superconductivity depend strongly on the N/O ratio, highlighting the critical role of disorder for superconductivity in this series of Ti(O,N) materials. △ Less

Submitted 4 October, 2023; originally announced October 2023.

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.

Comments: 16 pages, 4 figures. Revised manuscript with minor changes. Supplemental Information available on the journal website (open access)

Journal ref: Commun Phys 7, 33 (2024)

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.

Comments: 9 pages, 5 figures. Revised manuscript. Supplemental Information available on the journal website (open access)

Journal ref: Nat Commun 15, 441 (2024)

Search for ambient superconductivity in the Lu-N-H system

Authors: Pedro P. Ferreira, Lewis J. Conway, Alessio Cucciari, Simone Di Cataldo, Federico Giannessi, Eva Kogler, Luiz T. F. Eleno, Chris J. Pickard, Christoph Heil, Lilia Boeri

Abstract: Motivated by the recent report of room-temperature superconductivity at near-ambient pressure in N-doped lutetium hydride by Dasenbrock et al. [Nature 615, 244 (2023)], we performed a comprehensive, detailed study of the phase diagram of the Lu-N-H system, looking for superconducting phases. We combined ab initio crystal structure prediction with ephemeral data-derived interatomic potentials to sa… ▽ More Motivated by the recent report of room-temperature superconductivity at near-ambient pressure in N-doped lutetium hydride by Dasenbrock et al. [Nature 615, 244 (2023)], we performed a comprehensive, detailed study of the phase diagram of the Lu-N-H system, looking for superconducting phases. We combined ab initio crystal structure prediction with ephemeral data-derived interatomic potentials to sample over 200,000 different structures. Out of the more than 150 structures predicted to be metastable within $\sim$ 50 meV from the convex hull we identify 52 viable candidates for conventional superconductivity, for which we computed their superconducting properties from Density Functional Perturbation Theory. Although for some of these structures we do predict a finite superconducting $T_{\text{c}}$, none is even remotely compatible with room-temperature superconductivity as reported by Dasenbrock et al. Our work joins the broader community effort that has followed the report of near-ambient superconductivity, confirming beyond reasonable doubt that no conventional mechanism can explain the reported $T_{\text{c}}$ in Lu-N-H. △ Less

Submitted 6 September, 2023; v1 submitted 10 April, 2023; originally announced April 2023.

Journal ref: Nat Commun 14, 5367 (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.

Comments: 10 pages, 5 figures. Revised manuscript with minor changes. Supplemental Material available on the journal website (open access)

Journal ref: Commun Phys 6, 298 (2023)

Computational Reverse-Engineering Analysis for Scattering Experiments for Form Factor and Structure Factor Determination ('P(q) and S(q) CREASE')

Authors: Christian M. Heil, Yingzhen Ma, Bhuvnesh Bharti, Arthi Jayaraman

Abstract: In this paper we present an open-source machine learning (ML) accelerated computational method to analyze small-angle scattering profiles [I(q) vs. q] from concentrated macromolecular solutions to simultaneously obtain the form factor P(q) (e.g., dimensions of a micelle) and structure factor S(q) (e.g., spatial arrangement of the micelles) without relying on analytical models. This method builds o… ▽ More In this paper we present an open-source machine learning (ML) accelerated computational method to analyze small-angle scattering profiles [I(q) vs. q] from concentrated macromolecular solutions to simultaneously obtain the form factor P(q) (e.g., dimensions of a micelle) and structure factor S(q) (e.g., spatial arrangement of the micelles) without relying on analytical models. This method builds on our recent work on Computational Reverse Engineering Analysis for Scattering Experiments (CREASE) that has either been applied to obtain P(q) from dilute macromolecular solutions (where S(q) ~1) or to obtain S(q) from concentrated particle solution when the P(q) is known (e.g., sphere form factor). This paper's newly developed CREASE that calculates P(q) and S(q), termed as 'P(q) and S(q) CREASE' is validated by taking as input I(q) vs. q from in silico structures of known polydisperse core(A)-shell(B) micelles in solutions at varying concentrations and micelle-micelle aggregation. We demonstrate how 'P(q) and S(q) CREASE' performs if given two or three of the relevant scattering profiles - Itotal(q), IA(q), and IB(q) - as inputs; this demonstration is meant to guide experimentalists who may choose to do small-angle X-ray scattering (for total scattering from the micelles) and/or small-angle neutron scattering with appropriate contrast matching to get scattering solely from one or the other component (A or B). After validation of 'P(q) and S(q) CREASE' on in silico structures, we present our results analyzing small-angle neutron scattering profiles from a solution of core-shell type surfactant coated nanoparticles with varying extents of aggregation. △ Less

Submitted 6 December, 2022; originally announced December 2022.

Comments: 41 pages, 7 figures, 1 ToC figure

Superconductivity in Te-deficient ZrTe$_2$

Authors: L. E. Correa, P. P. Ferreira, L. R. de Faria, V. M. Fim, M. S. da Luz, M. S. Torikachvili, C. Heil, L. T. F. Eleno, A. J. S. Machado

Abstract: We present structural, electrical, and thermoelectric potential measurements on high-quality single crystals of ZrTe$_{1.8}$ grown from isothermal chemical vapor transport. These measurements show that the Te-deficient ZrTe$_{1.8}$, which forms the same structure as the non-superconducting ZrTe$_2$, is superconducting below 3.2\,K. The temperature dependence of the upper critical field (H$_{c2}$)… ▽ More We present structural, electrical, and thermoelectric potential measurements on high-quality single crystals of ZrTe$_{1.8}$ grown from isothermal chemical vapor transport. These measurements show that the Te-deficient ZrTe$_{1.8}$, which forms the same structure as the non-superconducting ZrTe$_2$, is superconducting below 3.2\,K. The temperature dependence of the upper critical field (H$_{c2}$) deviates from the behavior expected in conventional single-band superconductors, being best described by an electron-phonon two-gap superconducting model with strong intraband coupling. For the ZrTe$_{1.8}$ single crystals, the Seebeck potential measurements suggest that the charge carriers are predominantly negative, in agreement with the ab initio calculations. Through first-principles calculations within DFT, we show that the slight reduction of Te occupancy in ZrTe$_2$ unexpectedly gives origin to density of states peaks at the Fermi level due to the formation of localized Zr-$d$ bands, possibly promoting electronic instabilities at the Fermi level and an increase at the critical temperature according to the standard BCS theory. These findings highlight that the Te deficiency promotes the electronic conditions for the stability of the superconducting ground state, suggesting that defects can fine-tune the electronic structure to support superconductivity. △ Less

Submitted 10 April, 2023; v1 submitted 5 December, 2022; originally announced December 2022.

Journal ref: J. Phys. Chem. C 2023, 127, 10, 5162-5168

Mechanism of Structural Colors in Binary Mixtures of Nanoparticle-based Supraballs

Authors: Christian M. Heil, Anvay Patil, Bram Vanthournout, Saranshu Singla, Markus Bleuel, Jing-Jin Song, Ziying Hu, Nathan C. Gianneschi, Matthew D. Shawkey, Sunil K. Sinha, Arthi Jayaraman, Ali Dhinojwala

Abstract: Inspired by structural colors in avian species, various synthetic strategies have been developed to produce non-iridescent, saturated colors using nanoparticle assemblies. Mixtures of nanoparticles varying in particle chemistry (or complex refractive indices) and particle size have additional emergent properties that impact the color produced. For such complex multi-component systems, an understan… ▽ More Inspired by structural colors in avian species, various synthetic strategies have been developed to produce non-iridescent, saturated colors using nanoparticle assemblies. Mixtures of nanoparticles varying in particle chemistry (or complex refractive indices) and particle size have additional emergent properties that impact the color produced. For such complex multi-component systems, an understanding of assembled structure along with a robust optical modeling tool can empower scientists to perform intensive structure-color relationship studies and fabricate designer materials with tailored color. Here, we demonstrate how we can reconstruct the assembled structure from small-angle scattering measurements using the computational reverse-engineering analysis for scattering experiments (CREASE) method and then use the reconstructed structure in finite-difference time-domain (FDTD) calculations to predict color. We successfully, quantitatively predict experimentally observed color in mixtures containing strongly absorbing melanin nanoparticles and demonstrate the influence of a single layer of segregated nanoparticles on color produced. The versatile computational approach presented in this work is useful for engineering synthetic materials with desired colors without laborious trial and error experiments. △ Less

Submitted 9 October, 2022; originally announced October 2022.

Comments: 23 Pages, 5 Figures, 1 ToC Figure

Modeling Structural Colors from Disordered One-Component Colloidal Nanoparticle-based Supraballs using Combined Experimental and Simulation Techniques

Authors: Anvay Patil, Christian M. Heil, Bram Vanthournout, Saranshu Singla, Ziying Hu, Jan Ilavsky, Nathan C. Gianneschi, Matthew D. Shawkey, Sunil K. Sinha, Arthi Jayaraman, Ali Dhinojwala

Abstract: Bright, saturated structural colors in birds have inspired synthesis of self-assembled, disordered arrays of assembled nanoparticles with varied particle spacings and refractive indices. However, predicting colors of assembled nanoparticles, and thereby guiding their synthesis, remains challenging due to the effects of multiple scattering and strong absorption. Here, we use a computational approac… ▽ More Bright, saturated structural colors in birds have inspired synthesis of self-assembled, disordered arrays of assembled nanoparticles with varied particle spacings and refractive indices. However, predicting colors of assembled nanoparticles, and thereby guiding their synthesis, remains challenging due to the effects of multiple scattering and strong absorption. Here, we use a computational approach to first reconstruct the nanoparticles' assembled structures from small-angle scattering measurements and then input the reconstructed structures to a finite-difference time-domain method to predict their color and reflectance. This computational approach is successfully validated by comparing its predictions against experimentally measured reflectance and provides a pathway for reverse engineering colloidal assemblies with desired optical and photothermal properties. △ Less

Submitted 7 October, 2022; v1 submitted 25 May, 2022; originally announced May 2022.

Comments: 14 pages, 3 figures, 1 ToC figure

Journal ref: ACS Materials Letters 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.

Comments: 9 pages, 6 figures. Revised manuscript with minor changes. Supplemental Material available on the journal website (open access)

Journal ref: npj Comput Mater 8, 119 (2022)

Structural Color Production in Melanin-based Disordered Colloidal Nanoparticle Assemblies in Spherical Confinement

Authors: Anvay Patil, Christian M. Heil, Bram Vanthournout, Markus Bleuel, Saranshu Singla, Ziying Hu, Nathan C. Gianneschi, Matthew D. Shawkey, Sunil K. Sinha, Arthi Jayaraman, Ali Dhinojwala

Abstract: Melanin is a ubiquitous natural pigment that exhibits broadband absorption and high refractive index. Despite its widespread use in structural color production, how the absorbing material, melanin, affects the generated color is unknown. Using a combined molecular dynamics and finite-difference time-domain computational approach, this paper investigates structural color generation in one-component… ▽ More Melanin is a ubiquitous natural pigment that exhibits broadband absorption and high refractive index. Despite its widespread use in structural color production, how the absorbing material, melanin, affects the generated color is unknown. Using a combined molecular dynamics and finite-difference time-domain computational approach, this paper investigates structural color generation in one-component melanin nanoparticle-based supra-assemblies (called supraballs) as well as binary mixtures of melanin and silica (non-absorbing) nanoparticle-based supraballs. Experimentally produced one-component melanin and one-component silica supraballs, with thoroughly characterized primary particle characteristics using neutron scattering, produce reflectance profiles similar to the computational analogues, confirming that the computational approach correctly simulates both absorption and multiple scattering from the self-assembled nanoparticles. These combined approaches demonstrate that melanin's broadband absorption increases the primary reflectance peak wavelength, increases saturation, and decreases lightness factor. In addition, the dispersity of nanoparticle size more strongly influences the optical properties of supraballs than packing fraction, as evidenced by production of a larger range of colors when size dispersity is varied versus packing fraction. For binary melanin and silica supraballs, the chemistry-based stratification allows for more diverse color generation and finer saturation tuning than does the degree of mixing/demixing between the two chemistries. △ Less

Submitted 4 July, 2022; v1 submitted 20 November, 2021; originally announced November 2021.

Comments: 40 pages, Figure 6

LaBH$_{8}$: the first high-T$_{c}$ low-pressure superhydride

Authors: Simone Di Cataldo, Christoph Heil, Wolfgang von der Linden, Lilia Boeri

Abstract: In the last five years a large number of new high-temperature superconductors have been predicted and experimentally discovered among hydrogen-rich crystals, at pressures which are way too high to meet any practical application. In this work, we report the computational prediction of a hydride superconductor, LaBH$_{8}$, with a T$_{c}$ of 126 K at a pressure of 50 GPa, thermodynamically stable abo… ▽ More In the last five years a large number of new high-temperature superconductors have been predicted and experimentally discovered among hydrogen-rich crystals, at pressures which are way too high to meet any practical application. In this work, we report the computational prediction of a hydride superconductor, LaBH$_{8}$, with a T$_{c}$ of 126 K at a pressure of 50 GPa, thermodynamically stable above 100 GPa, and dynamically stable down to 40 GPa, an unprecedentedly low pressure for high-T$_{c}$ hydrides. LaBH$_{8}$ can be seen as a ternary sodalite-like hydride, in which a metallic hydrogen sublattice is stabilized by the chemical pressure exerted by the guest elements. The combination of two elements with different atomic sizes in LaBH$_{8}$ realizes a more efficient packing of atoms than in binary sodalite hydrides. A suitable choice of elements may be exploited to further reduce the stabilization pressure to ambient conditions. △ Less

Submitted 24 February, 2021; v1 submitted 22 February, 2021; originally announced February 2021.

Comments: 6 pages, 4 figures, 1 table

Journal ref: Phys. Rev. B 104, 020511 (2021)

The 2021 Quantum Materials Roadmap

Authors: Feliciano Giustino, Jin Hong Lee, Felix Trier, Manuel Bibes, Stephen M Winter, Roser Valentí, Young-Woo Son, Louis Taillefer, Christoph Heil, Adriana I. Figueroa, Bernard Plaçais, QuanSheng Wu, Oleg V. Yazyev, Erik P. A. M. Bakkers, Jesper Nygård, Pol Forn-Diaz, Silvano De Franceschi, J. W. McIver, L. E. F. Foa Torres, Tony Low, Anshuman Kumar, Regina Galceran, Sergio O. Valenzuela, Marius V. Costache, Aurélien Manchon , et al. (4 additional authors not shown)

Abstract: In recent years, the notion of Quantum Materials has emerged as a powerful unifying concept across diverse fields of science and engineering, from condensed-matter and cold atom physics to materials science and quantum computing. Beyond traditional quantum materials such as unconventional superconductors, heavy fermions, and multiferroics, the field has significantly expanded to encompass topologi… ▽ More In recent years, the notion of Quantum Materials has emerged as a powerful unifying concept across diverse fields of science and engineering, from condensed-matter and cold atom physics to materials science and quantum computing. Beyond traditional quantum materials such as unconventional superconductors, heavy fermions, and multiferroics, the field has significantly expanded to encompass topological quantum matter, two-dimensional materials and their van der Waals heterostructures, Moire materials, Floquet time crystals, as well as materials and devices for quantum computation with Majorana fermions. In this Roadmap collection we aim to capture a snapshot of the most recent developments in the field, and to identify outstanding challenges and emerging opportunities. The format of the Roadmap, whereby experts in each discipline share their viewpoint and articulate their vision for quantum materials, reflects the dynamic and multifaceted nature of this research area, and is meant to encourage exchanges and discussions across traditional disciplinary boundaries. It is our hope that this collective vision will contribute to sparking new fascinating questions and activities at the intersection of materials science, condensed matter physics, device engineering, and quantum information, and to shaping a clearer landscape of quantum materials science as a new frontier of interdisciplinary scientific inquiry. △ Less

Submitted 4 February, 2021; originally announced February 2021.

Comments: 93 pages, 27 figures

Journal ref: J. Phys. Mater. 3 042006 (2020)

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.

Comments: 9 pages, 6 figures and Supplemental Material

Journal ref: J. Phys.: Condens. Matter (2021)

Magnetic surface reconstruction in the van-der-Waals antiferromagnet Fe$_{1+x}$Te

Authors: C. Trainer, M. Songvilay, N. Qureshi, A. Stunault, C. M. Yim, E. E. Rodriguez, C. Heil, V. Tsurkan, M. A. Green, A. Loidl, P. Wahl, C. Stock

Abstract: Fe$_{1+x}$Te is a two dimensional van der Waals antiferromagnet that becomes superconducting on anion substitution on the Te site. The parent phase of Fe$_{1+x}$Te is sensitive to the amount of interstitial iron situated between the iron-tellurium layers displaying collinear magnetic order coexisting with low temperature metallic resistivity for small concentrations of interstitial iron $x$ and he… ▽ More Fe$_{1+x}$Te is a two dimensional van der Waals antiferromagnet that becomes superconducting on anion substitution on the Te site. The parent phase of Fe$_{1+x}$Te is sensitive to the amount of interstitial iron situated between the iron-tellurium layers displaying collinear magnetic order coexisting with low temperature metallic resistivity for small concentrations of interstitial iron $x$ and helical magnetic order for large values of $x$. While this phase diagram has been established through scattering [see for example E. E. Rodriguez $\textit{et al.}$ Phys. Rev. B ${\bf{84}}$, 064403 (2011) and S. Rössler $\textit{et al.}$ Phys. Rev. B ${\bf{84}}$, 174506 (2011)], recent scanning tunnelling microscopy measurements [C. Trainer $\textit{et al.}$ Sci. Adv. ${\bf{5}}$, eaav3478 (2019)] have observed a different magnetic structure for small interstitial iron concentrations $x$ with a significant canting of the magnetic moments along the crystallographic $c$ axis of $θ$=28 $\pm$ 3$^{\circ}$. In this paper, we revisit the magnetic structure of Fe$_{1.09}$Te using spherical neutron polarimetry and scanning tunnelling microscopy to search for this canting in the bulk phase and compare surface and bulk magnetism. The results show that the bulk magnetic structure of Fe$_{1.09}$Te is consistent with collinear in-plane order ($θ=0$ with an error of $\sim$ 5$^{\circ}$). Comparison with scanning tunnelling microscopy on a series of Fe$_{1+x}$Te samples reveals that the surface exhibits a magnetic surface reconstruction with a canting angle of the spins of $θ=29.8^{\circ}$. We suggest that this is a consequence of structural relaxation of the surface layer resulting in an out-of-plane magnetocrystalline anisotropy. The magnetism in Fe$_{1+x}$Te displays different properties at the surface when the symmetry constraints of the bulk are removed. △ Less

Submitted 3 January, 2021; originally announced January 2021.

Comments: 11 pages, 6 figures, to be Published in Phys. Rev. B

Journal ref: Phys. Rev. B 103, 024406 (2021)

Electronic, vibrational, and electron-phonon coupling properties in SnSe$_2$ and SnS$_2$ under pressure

Authors: Gyanu Prasad Kafle, Christoph Heil, Hari Paudyal, Elena R. Margine

Abstract: The tin-selenide and tin-sulfide classes of materials undergo multiple structural transitions under high pressure leading to periodic lattice distortions, superconductivity, and topologically non-trivial phases, yet a number of controversies exist regarding the structural transformations in these systems. We perform first-principles calculations within the framework of density functional theory an… ▽ More The tin-selenide and tin-sulfide classes of materials undergo multiple structural transitions under high pressure leading to periodic lattice distortions, superconductivity, and topologically non-trivial phases, yet a number of controversies exist regarding the structural transformations in these systems. We perform first-principles calculations within the framework of density functional theory and a careful comparison of our results with available experiments on SnSe$_2$ reveals that the apparent contradictions among high-pressure results can be attributed to differences in experimental conditions. We further demonstrate that under hydrostatic pressure a $\sqrt{3} \times \sqrt{3} \times 1$ superstructure can be stabilized above 20 GPa in SnS$_2$ via a periodic lattice distortion as found recently in the case of SnSe$_2$, and that this pressure-induced phase transition is due to the combined effect of Fermi surface nesting and electron-phonon coupling at a momentum wave vector $\mathbf{q}$ = $(1/3, 1/3, 0)$. In addition, we investigate the contribution of nonadiabatic corrections on the calculated phonon frequencies, and show that the quantitative agreement between theory and experiment for the high-energy $A_{1g}$ phonon mode is improved when these effects are taken into account. Finally, we examine the nature of the superconducting state recently observed in SnSe$_2$ under nonhydrostatic pressure and predict the emergence of superconductivity with a comparable critical temperature in SnS$_2$ under similar experimental conditions. Interestingly, in the periodic lattice distorted phases, the critical temperature is found to be reduced by an order of magnitude due to the restructuring of the Fermi surface. △ Less

Submitted 22 October, 2020; originally announced October 2020.

Comments: accepted manuscript, 15 pages, 11 figures, and supplemental material

Journal ref: Journal of Materials Chemistry C, 2020

Electronic structure and superconductivity of the non-centrosymmetric Sn$_4$As$_3$

Authors: C. A. Marques, M. J. Neat, C. M. Yim, M. D. Watson, L. C. Rhodes, C. Heil, K. S. Pervakov, V. A. Vlasenko, V. M. Pudalov, A. V. Muratov, T. K. Kim, P. Wahl

Abstract: In a superconductor that lacks inversion symmetry, the spatial part of the Cooper pair wave function has a reduced symmetry, allowing for the mixing of spin-singlet and spin-triplet Cooper pairing channels and thus providing a pathway to a non-trivial superconducting state. Materials with a non-centrosymmetric crystal structure and with strong spin-orbit coupling are a platform to realize these po… ▽ More In a superconductor that lacks inversion symmetry, the spatial part of the Cooper pair wave function has a reduced symmetry, allowing for the mixing of spin-singlet and spin-triplet Cooper pairing channels and thus providing a pathway to a non-trivial superconducting state. Materials with a non-centrosymmetric crystal structure and with strong spin-orbit coupling are a platform to realize these possibilities. Here, we report the synthesis and characterisation of high quality crystals of Sn$_4$As$_3$, with non-centrosymmetric unit cell ($R3m$). We have characterised the normal and superconducting state using a range of methods. Angle-resolved photoemission spectroscopy shows a multiband Fermi surface and the presence of two surface states, confirmed by Density-functional theory calculations. Specific heat measurements reveal a superconducting critical temperature of $T_c\sim 1.14$ K and an upper critical magnetic field of $H_c\gtrsim 7$ mT, which are both confirmed by ultra-low temperature scanning tunneling microscopy and spectroscopy. Scanning tunneling spectroscopy shows a fully formed superconducting gap, consistent with conventional $s$-wave superconductivity. △ Less

Submitted 29 June, 2020; v1 submitted 13 December, 2019; originally announced December 2019.

Comments: 11 pages, 7 figures

Journal ref: New J. Phys. 22 (2020) 063049

Superconductivity in sodalite-like yttrium hydride clathrates

Authors: Christoph Heil, Simone di Cataldo, Giovanni B. Bachelet, Lilia Boeri, ITP-CP, Graz University of Technology, Austria, Dipartimento di Fisica, Sapienza Universita' di Roma, Italy

Abstract: We report ab-initio calculations of the superconducting properties of two high-Tc sodalite-like clathrate yttrium hydrides, YH6 and YH10, within the fully anisotropic ME theory, including Coulomb corrections. For both compounds we find almost isotropic superconducting gaps, resulting from a uniform distribution of the electron-phonon coupling over phonon modes and electronic states of mixed Y and… ▽ More We report ab-initio calculations of the superconducting properties of two high-Tc sodalite-like clathrate yttrium hydrides, YH6 and YH10, within the fully anisotropic ME theory, including Coulomb corrections. For both compounds we find almost isotropic superconducting gaps, resulting from a uniform distribution of the electron-phonon coupling over phonon modes and electronic states of mixed Y and H character. The Coulomb screening is rather weak, resulting in a Morel-Anderson pseudopotential mu*= 0:11, at odds with claims of unusually large Tc in lanthanum hydrides. The corresponding critical temperatures at 300 GPa exceed room temperature (Tc = 290 K and 310 K for YH6 and YH10), in agreement with a previous isotropic-gap calculation. The different response of these two compounds to external pressure, along with a comparison to low-Tc superconducting YH3, may inspire strategies to improve the superconducting properties of this class of hydrides. △ Less

Submitted 13 January, 2019; originally announced January 2019.

Journal ref: Phys. Rev. B 99, 220502 (2019)

Quasiparticle $GW$ band structures and Fermi surfaces of bulk and monolayer NbS$_2$

Authors: Christoph Heil, Martin Schlipf, Feliciano Giustino

Abstract: In this work we employ the $GW$ approximation in the framework of the SternheimerGW method to investigate the effects of many-body corrections to the band structures and Fermi surfaces of bulk and monolayer NbS$_2$. For the bulk system, we find that the inclusion of these many-body effects leads to important changes in the band structure, especially in the low-energy regime around the Fermi level,… ▽ More In this work we employ the $GW$ approximation in the framework of the SternheimerGW method to investigate the effects of many-body corrections to the band structures and Fermi surfaces of bulk and monolayer NbS$_2$. For the bulk system, we find that the inclusion of these many-body effects leads to important changes in the band structure, especially in the low-energy regime around the Fermi level, and that our calculations are in good agreement with recent ARPES measurements. In the case of a free-standing monolayer NbS$_2$, we observe a strong increase of the screened Coulomb interaction and the quasiparticle corrections as compared to bulk. In this case we also perform calculations to include the effect of screening by a substrate. We report in detail the results of our convergence tests and computational parameters, to serve as a solid basis for future studies. △ Less

Submitted 10 August, 2018; v1 submitted 24 July, 2018; originally announced July 2018.

Comments: 15 pages, 18 figures

Journal ref: Phys. Rev. B 98, 075120 (2018)

Unusual pressure-induced periodic lattice distortion in SnSe$_2$

Authors: Jianjun Ying, Hari Paudyal, Christoph Heil, Xiao-Jia Chen, Viktor V. Struzhkin, Elena R. Margine

Abstract: We performed high pressure x-ray diffraction (XRD), Raman, and transport measurements combined with first-principles calculations to investigate the behavior of tin diselenide (SnSe$_2$) under compression. The obtained single-crystal XRD data indicate the formation of a $(1/3,1/3,1)$-type superlattice above 17 GPa. According to our density functional theory results, the pressure-induced transition… ▽ More We performed high pressure x-ray diffraction (XRD), Raman, and transport measurements combined with first-principles calculations to investigate the behavior of tin diselenide (SnSe$_2$) under compression. The obtained single-crystal XRD data indicate the formation of a $(1/3,1/3,1)$-type superlattice above 17 GPa. According to our density functional theory results, the pressure-induced transition to the commensurate periodic lattice distortion (PLD) phase is due to the combined effect of strong Fermi surface nesting and electron-phonon coupling at a momentum wave vector $\mathbf{q}=(1/3,1/3,1)$. In contrast, similar PLD transitions associated with charge density wave (CDW) orderings in transition metal dichalcogenides (TMDs) do not involve significant Fermi surface nesting. The discovered pressure-induced PLD is quite remarkable, as pressure usually suppresses CDW phases in related materials. Our findings, therefore, provide new playgrounds to study the intricate mechanisms governing the emergence of PLD in TMD-related materials. △ Less

Submitted 13 July, 2018; originally announced July 2018.

Comments: 6 pages, 4 figures + Supplemental (12 pages, 14 figures)

Journal ref: Phys. Rev. Lett. 121, 027003 (2018)

Absence of superconductivity in iron polyhydrides at high pressures

Authors: Christoph Heil, Giovanni B. Bachelet, Lilia Boeri

Abstract: Recently, C. M. Pépin \textit{et al.} [Science \textbf{357}, 382 (2017)] reported the formation of several new iron polyhydrides FeH$_x$ at pressures in the megabar range, and spotted FeH$_5$, which forms above 130 GPa, as a potential high-\tc \ superconductor, because of an alleged layer of dense metallic hydrogen. Shortly after, two studies by A.~Majumdar \textit{et al.} [Phys. Rev. B \textbf{96… ▽ More Recently, C. M. Pépin \textit{et al.} [Science \textbf{357}, 382 (2017)] reported the formation of several new iron polyhydrides FeH$_x$ at pressures in the megabar range, and spotted FeH$_5$, which forms above 130 GPa, as a potential high-\tc \ superconductor, because of an alleged layer of dense metallic hydrogen. Shortly after, two studies by A.~Majumdar \textit{et al.} [Phys. Rev. B \textbf{96}, 201107 (2017)] and A.~G.~Kvashnin \textit{et al.} [J. Phys. Chem. C \textbf{122}, 4731 (2018)] based on {\em ab initio} Migdal-Eliashberg theory seemed to independently confirm such a conjecture. We conversely find, on the same theoretical-numerical basis, that neither FeH$_5$ nor its precursor, FeH$_3$, shows any conventional superconductivity and explain why this is the case. We also show that superconductivity may be attained by transition-metal polyhydrides in the FeH$_3$ structure type by adding more electrons to partially fill one of the Fe--H hybrid bands (as, e.g., in NiH$_3$). Critical temperatures, however, will remain low because the $d$--metal bonding, and not the metallic hydrogen, dominates the behavior of electrons and phonons involved in the superconducting pairing in these compounds. △ Less

Submitted 14 June, 2018; v1 submitted 10 April, 2018; originally announced April 2018.

Comments: 8 pages, 8 figures

Journal ref: Phys. Rev. B 97, 214510 (2018)

Manipulating surface magnetic order in iron telluride

Authors: Christopher Trainer, Chi Ming Yim, Christoph Heil, Feliciano Giustino, Dorina Croitori, Vladimir Tsurkan, Alois Loidl, Efrain Rodriguez, Chris Stock, Peter Wahl

Abstract: Control of emergent magnetic orders in correlated electron materials promises new opportunities for applications in spintronics. For their technological exploitation, it is important to understand the role of surfaces and interfaces to other materials and their impact on the emergent magnetic orders. Here, we demonstrate for iron telluride, the nonsuperconducting parent compound of the iron chalco… ▽ More Control of emergent magnetic orders in correlated electron materials promises new opportunities for applications in spintronics. For their technological exploitation, it is important to understand the role of surfaces and interfaces to other materials and their impact on the emergent magnetic orders. Here, we demonstrate for iron telluride, the nonsuperconducting parent compound of the iron chalcogenide superconductors, determination and manipulation of the surface magnetic structure by low-temperature spin-polarized scanning tunneling microscopy. Iron telluride exhibits a complex structural and magnetic phase diagram as a function of interstitial iron concentration. Several theories have been put forward to explain the different magnetic orders observed in the phase diagram, which ascribe a dominant role either to interactions mediated by itinerant electrons or to local moment interactions. Through the controlled removal of surface excess iron, we can separate the influence of the excess iron from that of the change in the lattice structure. △ Less

Submitted 4 March, 2019; v1 submitted 16 February, 2018; originally announced February 2018.

Comments: 8 pages, 6 figures, published version

Journal ref: Sci. Adv. 5, eaav3478 (2019)

Origin of superconductivity and latent charge density wave in NbS$_2$

Authors: Christoph Heil, Samuel Poncé, Henry Lambert, Martin Schlipf, Elena R. Margine, Feliciano Giustino

Abstract: We elucidate the origin of the phonon-mediated superconductivity in 2$H$-NbS$_2$ using the ab initio anisotropic Migdal-Eliashberg theory including Coulomb interactions. We demonstrate that superconductivity is associated with Fermi surface hot spots exhibiting an unusually strong electron-phonon interaction. The electron-lattice coupling is dominated by low-energy anharmonic phonons, which place… ▽ More We elucidate the origin of the phonon-mediated superconductivity in 2$H$-NbS$_2$ using the ab initio anisotropic Migdal-Eliashberg theory including Coulomb interactions. We demonstrate that superconductivity is associated with Fermi surface hot spots exhibiting an unusually strong electron-phonon interaction. The electron-lattice coupling is dominated by low-energy anharmonic phonons, which place the system on the verge of a charge density wave instability. We also provide definitive evidence for two-gap superconductivity in 2$H$-NbS$_2$, and show that the low- and high-energy peaks observed in tunneling spectra correspond to the $Γ$- and $K$-centered Fermi surface pockets, respectively. The present findings call for further efforts to determine whether our proposed mechanism underpins superconductivity in the whole family of metallic transition metal dichalcogenides. △ Less

Submitted 4 August, 2017; originally announced August 2017.

Comments: 6 pages, 5 figures and Supplemental Material

Journal ref: Phys. Rev. Lett. 119, 087003 (2017)

Search for high-Tc conventional superconductivity at megabar pressures in the lithium-sulfur system

Authors: Christian Kokail, Christoph Heil, Lilia Boeri

Abstract: Motivated by the recent report of superconductivity above 200 K in ultra-dense hydrogen sulfide, we search for high-\tc\ conventional superconductivity in the phase diagram of the binary Li-S system, using {\em ab-initio} methods for crystal structure prediction and linear response calculations for the electron-phonon coupling. We find that at pressures higher than 20 GPa, several new compositions… ▽ More Motivated by the recent report of superconductivity above 200 K in ultra-dense hydrogen sulfide, we search for high-\tc\ conventional superconductivity in the phase diagram of the binary Li-S system, using {\em ab-initio} methods for crystal structure prediction and linear response calculations for the electron-phonon coupling. We find that at pressures higher than 20 GPa, several new compositions, besides the known Li$_2$S, are stabilized; many exhibit electride-like {\em interstitial} charge localization observed in other alkali metal compounds. Of all predicted phases, only Li$_3$S at P > 640 GPa displays a sizable \tc, in contrast to what is observed in sulfur and phosphorus hydrides, where several stoichiometries lead to high \tc. We attribute this difference to 2$s$-2$p$ hybridization and avoided core overlap, and predict similar behavior for other alkali metal compounds. △ Less

Submitted 24 May, 2016; originally announced May 2016.

Journal ref: Phys. Rev. B 94, 060502 (2016)

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

Authors: José A. Flores-Livas, Maximilian Amsler, Christoph Heil, Antonio Sanna, Lilia Boeri, Gianni Profeta, Chris Wolverton, Stefan Goedecker, E. K. U. Gross

Abstract: 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 w… ▽ More 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$. △ Less

Submitted 7 December, 2015; originally announced December 2015.

Journal ref: Phys. Rev. B 93, 020508 (2016)

Influence of bonding on superconductivity in high-pressure hydrides

Authors: Christoph Heil, Lilia Boeri

Abstract: The recent reports on high-temperature superconductivity above 190 K in hydrogen sulfide at 200 GPa pressure, exceeding all previously discovered superconductors, has greatly invigorated the interest in dense hydrogen-rich solids. In this paper, we investigate a possible way to optimize the critical temperature in these compounds using first-principles linear-response calculations. We construct hy… ▽ More The recent reports on high-temperature superconductivity above 190 K in hydrogen sulfide at 200 GPa pressure, exceeding all previously discovered superconductors, has greatly invigorated the interest in dense hydrogen-rich solids. In this paper, we investigate a possible way to optimize the critical temperature in these compounds using first-principles linear-response calculations. We construct hypothetical alchemical atoms to smoothly interpolate between elements of the chalcogen group and study their bonding and superconducting properties. Our results show that the already remarkable critical temperatures of H$_3$S could be improved even further by increasing the ionic character of the relevant bonds, i.e. replacing sulfur partially with more electronegative elements. △ Less

Submitted 31 August, 2015; v1 submitted 9 July, 2015; originally announced July 2015.

Comments: 5 pages, 4 figures

Journal ref: Phys. Rev. B 92, 060508 (2015)

Effect of the iron valence in the two types of layers in LiFeO$_2$Fe$_2$Se$_2$

Authors: Christoph Heil, Lilia Boeri, Heinrich Sormann, Wolfgang von der Linden, Markus Aichhorn

Abstract: We perform electronic structure calculations for the recently synthesized iron-based superconductor LiFeO$_2$Fe$_2$Se$_2$. In contrast to other iron-based superconductors, this material comprises two different iron atoms in 3$d^5$ and 3$d^6$ configurations. In band theory, both contribute to the low-energy electronic structure. Spin-polarized density functional theory calculations predict an antif… ▽ More We perform electronic structure calculations for the recently synthesized iron-based superconductor LiFeO$_2$Fe$_2$Se$_2$. In contrast to other iron-based superconductors, this material comprises two different iron atoms in 3$d^5$ and 3$d^6$ configurations. In band theory, both contribute to the low-energy electronic structure. Spin-polarized density functional theory calculations predict an antiferromagnetic metallic ground state with different moments on the two Fe sites. However, several other almost degenerate magnetic configurations exist. Due to their different valences, the two iron atoms behave very differently when local quantum correlations are included through the dynamical mean-field theory. The contributions from the half-filled 3$d^5$ atoms in the LiFeO$_2$ layer are suppressed and the 3$d^6$ states from the FeSe layer restore the standard iron-based superconductor fermiology. △ Less

Submitted 16 October, 2014; v1 submitted 5 August, 2014; originally announced August 2014.

Comments: 9 pages, 11 figures

Journal ref: Phys. Rev. B 90, 165122 (2014)

Accurate bare susceptibilities from full-potential $\textit{ab initio}$ calculations

Authors: Christoph Heil, Heinrich Sormann, Lilia Boeri, Markus Aichhorn, Wolfgang von der Linden

Abstract: Electronic susceptibilities are a very popular tool to study electronic and magnetic properties of materials, both in experiment and theory. Unfortunately, the numerical evaluation of even the bare susceptibility, which depends on the computation of matrix elements and sums over energy bands, is very work-intensive and therefore various approximations have been introduced to speed up the calculati… ▽ More Electronic susceptibilities are a very popular tool to study electronic and magnetic properties of materials, both in experiment and theory. Unfortunately, the numerical evaluation of even the bare susceptibility, which depends on the computation of matrix elements and sums over energy bands, is very work-intensive and therefore various approximations have been introduced to speed up the calculations. We present a reliable and efficient implementation to compute static as well as dynamic bare susceptibilities based on full-potential density functional theory (DFT) calculations. Based on the exact results we will assess the accuracy of replacing the matrix elements with a constant and the impact of truncating the sum over the energy bands. Results will be given for representative and topical materials, such as Cr, a classical transition metal, as well as for FeSe and LaFeAsO, examples of iron-based superconductors. △ Less

Submitted 26 September, 2014; v1 submitted 13 July, 2014; originally announced July 2014.

Comments: 10 pages, 10 figures

Journal ref: Phys. Rev. B 90, 115143 (2014)

Indirect evidence for strong coupling superconductivity in FeSe under pressure from first-principle calculations

Authors: Christoph Heil, Markus Aichhorn, Heinrich Sormann, Ewald Schachinger, Wolfgang von der Linden

Abstract: Whether the pairing mechanism for superconductivity in iron-based superconductors is of itinerant or local character is a question still debated on. In order to investigate the influence of Fermi surface nesting, we calculate from first-principles the static susceptibility of FeSe under applied pressures, both hydrostatic and non-hydrostatic. We show that the electronic band structures are highly… ▽ More Whether the pairing mechanism for superconductivity in iron-based superconductors is of itinerant or local character is a question still debated on. In order to investigate the influence of Fermi surface nesting, we calculate from first-principles the static susceptibility of FeSe under applied pressures, both hydrostatic and non-hydrostatic. We show that the electronic band structures are highly sensitive to the way pressure is applied and confront our theoretical results with conclusions drawn from experiments. Given that the critical temperature of FeSe is quite universal as function of pressure, this is clear evidence that in FeSe the evolution of Fermi surface nesting cannot account for the evolution of the critical temperature with pressure. Hence we argue that the pairing in iron-chalcogenide compounds should be of strong-coupling spin fluctuation origin. △ Less

Submitted 9 October, 2012; originally announced October 2012.

Comments: 5 pages, 7 figures

Steady-state spectra, current and stability diagram of a quantum dot: a non-equilibrium Variational Cluster Approach

Authors: Martin Nuss, Christoph Heil, Martin Ganahl, Michael Knap, Hans Gerd Evertz, Enrico Arrigoni, Wolfgang von der Linden

Abstract: We calculate steady-state properties of a strongly correlated quantum dot under voltage bias by means of non-equilibrium Cluster Perturbation Theory and the non-equilibrium Variational Cluster Approach, respectively. Results for the steady-state current are benchmarked against data from accurate Matrix Product State based time evolution. We show that for low to medium interaction strength, non-equ… ▽ More We calculate steady-state properties of a strongly correlated quantum dot under voltage bias by means of non-equilibrium Cluster Perturbation Theory and the non-equilibrium Variational Cluster Approach, respectively. Results for the steady-state current are benchmarked against data from accurate Matrix Product State based time evolution. We show that for low to medium interaction strength, non-equilibrium Cluster Perturbation Theory already yields good results, while for higher interaction strength the self-consistent feedback of the non-equilibrium Variational Cluster Approach significantly enhances the accuracy. We report the current-voltage characteristics for different interaction strengths. Furthermore we investigate the non-equilibrium local density of states of the quantum dot and illustrate that within the variational approach a linear splitting and broadening of the Kondo resonance is predicted which depends on interaction strength. Calculations with applied gate voltage, away from particle hole symmetry, reveal that the maximum current is reached at the crossover from the Kondo regime to the doubly-occupied or empty quantum dot. Obtained stability diagrams compare very well to recent experimental data [Phys. Rev. B, 84, 245316 (2011)]. △ Less

Submitted 16 October, 2012; v1 submitted 24 July, 2012; originally announced July 2012.

Comments: 13 pages, 7 figures

Journal ref: Phys. Rev. B 86, 245119 (2012)

Strong coupling expansion for the Bose-Hubbard and the Jaynes-Cummings lattice model

Authors: Christoph Heil, Wolfgang von der Linden

Abstract: A strong coupling expansion, based on the Kato-Bloch perturbation theory, which has recently been proposed by Eckardt et al. [Phys. Rev. B 79, 195131] and Teichmann et al. [Phys. Rev. B 79, 224515] is implemented in order to study various aspects of the Bose-Hubbard and the Jaynes-Cummings lattice model. The approach, which allows to generate numerically all diagrams up to a desired order in the i… ▽ More A strong coupling expansion, based on the Kato-Bloch perturbation theory, which has recently been proposed by Eckardt et al. [Phys. Rev. B 79, 195131] and Teichmann et al. [Phys. Rev. B 79, 224515] is implemented in order to study various aspects of the Bose-Hubbard and the Jaynes-Cummings lattice model. The approach, which allows to generate numerically all diagrams up to a desired order in the interaction strength is generalized for disordered systems and for the Jaynes-Cummings lattice model. Results for the Bose-Hubbard and the Jaynes-Cummings lattice model will be presented and compared with results from VCA and DMRG. Our focus will be on the Mott insulator to superfluid transition. △ Less

Submitted 26 July, 2012; v1 submitted 12 May, 2011; originally announced May 2011.

Comments: 29 pages, 21 figures

Journal ref: 2012 J. Phys.: Condens. Matter 24 295601