-
Magnetic Correlations and Pairing Tendencies of the Hybrid Stacking Nickelate Superlattice La$_7$Ni$_5$O$_{17}$ (La$_3$Ni$_2$O$_7$/La$_4$Ni$_3$O$_{10}$) under Pressure
Authors:
Yang Zhang,
Ling-Fang Lin,
Adriana Moreo,
Thomas A. Maier,
Elbio Dagotto
Abstract:
Motivated by the recent rapid progress in high-$T_c$ nickelate superconductors, we comprehensively study the physical properties of the alternating bilayer trilayer stacking nickelate La$_7$Ni$_5$O$_{17}$. The high-symmetry phase of this material, without the tilting of oxygen octahedra, is not stable at ambient conditions but becomes stable under high pressure, where a small hole pocket $γ_0$, co…
▽ More
Motivated by the recent rapid progress in high-$T_c$ nickelate superconductors, we comprehensively study the physical properties of the alternating bilayer trilayer stacking nickelate La$_7$Ni$_5$O$_{17}$. The high-symmetry phase of this material, without the tilting of oxygen octahedra, is not stable at ambient conditions but becomes stable under high pressure, where a small hole pocket $γ_0$, composed of the $d_{3z^2-r^2}$ states in the trilayer sublattice, appears. This pocket was identified in our previous work for trilayer La$_4$Ni$_3$O$_{10}$ as important to develop superconductivity. Moreover, using random-phase approximation calculations, we find a leading $s^\pm$ pairing state for the high-symmetry phase under pressure with similar pairing strength as that obtained previously for the bilayer La$_3$Ni$_2$O$_7$ compound, suggesting a similar or higher superconducting transition temperature $T_c$. In addition, we find that the dominant magnetic fluctuations in the system driving this pairing state have antiferromagnetic structure both in-plane and between the planes of the top and bottom trilayer and bilayer sublattices, while the middle trilayer is magnetically decoupled.
△ Less
Submitted 14 August, 2024;
originally announced August 2024.
-
Magnetic phase diagram of a two-orbital model for bilayer nickelates varying doping
Authors:
Ling-Fang Lin,
Yang Zhang,
Nitin Kaushal,
Gonzalo Alvarez,
Thomas A. Maier,
Adriana Moreo,
Elbio Dagotto
Abstract:
Motivated by the recently discovered high-$T_c$ bilayer nickelate superconductor La$_3$Ni$_2$O$_7$, we comprehensively research a bilayer $2\times2\times2$ cluster for different electronic densities $n$ by using the Lanczos method. We also employ the random-phase approximation to quantify the first magnetic instability with increasing Hubbard coupling strength, also varying $n$. Based on the spin…
▽ More
Motivated by the recently discovered high-$T_c$ bilayer nickelate superconductor La$_3$Ni$_2$O$_7$, we comprehensively research a bilayer $2\times2\times2$ cluster for different electronic densities $n$ by using the Lanczos method. We also employ the random-phase approximation to quantify the first magnetic instability with increasing Hubbard coupling strength, also varying $n$. Based on the spin structure factor $S(q)$, we have obtained a rich magnetic phase diagram in the plane defined by $n$ and $U/W$, at fixed Hund coupling. We have observed numerous states, such as A-AFM, Stripes, G-AFM, and C-AFM. For half-filling $n=2$ (two electrons per Ni site, corresponding to $N$ = 16 electrons), the canonical superexchange interaction leads to a robust G-AFM state $(π,π,π)$ with antiferromagnetic couplings in plane and between layers. By increasing or decreasing electronic densities, ferromagnetic tendencies emerge from the ``half-empty'' and ``half-full'' mechanisms, leading to many other interesting magnetic tendencies. In addition, the spin-spin correlations become weaker both in the hole or electron doping regions compared with half-filling. At $n = 1.5$ (or $N=12$), density corresponding to La$_3$Ni$_2$O$_7$, we obtained the ``Stripe 2'' ground state (antiferromagnetic coupling in one in-plane direction, ferromagnetic coupling in the other, and antiferromagnetic coupling along the $z$-axis) in the $2\times2\times2$ cluster. In addition, we obtained a much stronger AFM coupling along the $z$-axis than the magnetic coupling in the $xy$ plane. The random-phase approximation calculations with varying $n$ give very similar results as Lanczos. Meanwhile, a state with $q/π= (0.6, 0.6, 1)$ close to the E-phase wavevector is found in our RPA calculations by slightly reducing the filling to $n=1.25$, possibly responsible for the E-phase SDW recently observed in experiments.
△ Less
Submitted 10 August, 2024;
originally announced August 2024.
-
SmoQyDEAC.jl: A differential evolution package for the analytic continuation of imaginary time correlation functions
Authors:
James Neuhaus,
Nathan S. Nichols,
Debshikha Banerjee,
Benjamin Cohen-Stead,
Thomas A. Maier,
Adrian Del Maestro,
Steven Johnston
Abstract:
We introduce the SmoQyDEAC.jl package, a Julia implementation of the Differential Evolution Analytic Continuation (DEAC) algorithm [N. S. Nichols et al., Phys. Rev. E 106, 025312 (2022)] for analytically continuing noisy imaginary time correlation functions to the real frequency axis. Our implementation supports fermionic and bosonic correlation functions on either the imaginary time or Matsubara…
▽ More
We introduce the SmoQyDEAC.jl package, a Julia implementation of the Differential Evolution Analytic Continuation (DEAC) algorithm [N. S. Nichols et al., Phys. Rev. E 106, 025312 (2022)] for analytically continuing noisy imaginary time correlation functions to the real frequency axis. Our implementation supports fermionic and bosonic correlation functions on either the imaginary time or Matsubara frequency axes, and treatment of the covariance error in the input data. This paper presents an overview of the DEAC algorithm and the features implemented in the SmoQyDEAC.jl. It also provides detailed benchmarks of the package's output against the popular maximum entropy and stochastic analytic continuation methods. The code for this package can be downloaded from our GitHub repository at https://github.com/SmoQySuite/SmoQyDEAC.jl or installed using the Julia package manager. The online documentation, including examples, can be accessed at https://smoqysuite.github.io/SmoQyDEAC.jl/stable/.
△ Less
Submitted 30 September, 2024; v1 submitted 5 July, 2024;
originally announced July 2024.
-
Fluctuating charge-density-wave correlations in the three-band Hubbard model
Authors:
P. Mai,
B. Cohen-Stead,
T. A. Maier,
S. Johnston
Abstract:
The high-temperature superconducting cuprates host unidirectional spin- and charge-density-wave orders that can intertwine with superconductivity in non-trivial ways. While the charge components of these stripes have now been observed in nearly all cuprate families, their detailed evolution with doping varies across different materials and at high and low temperatures. We address this problem usin…
▽ More
The high-temperature superconducting cuprates host unidirectional spin- and charge-density-wave orders that can intertwine with superconductivity in non-trivial ways. While the charge components of these stripes have now been observed in nearly all cuprate families, their detailed evolution with doping varies across different materials and at high and low temperatures. We address this problem using non-perturbative determinant quantum Monte Carlo calculations for the three-band Hubbard model. Using an efficient implementation, we can resolve the model's fluctuating spin and charge modulations and map their evolution as a function of the charge transfer energy and doping. We find that the incommensurability of the charge modulations is decoupled from the spin modulations and decreases with hole doping, consistent with experimental measurements at high temperatures. These findings support the proposal that the high-temperature charge correlations are distinct from the intertwined stripe order observed at low-temperature and in the single-band Hubbard model.
△ Less
Submitted 21 May, 2024;
originally announced May 2024.
-
Electronic structure, self-doping, and superconducting instability in the alternating single-layer trilayer stacking nickelates La$_3$Ni$_2$O$_7$
Authors:
Yang Zhang,
Ling-Fang Lin,
Adriana Moreo,
Thomas A. Maier,
Elbio Dagotto
Abstract:
Motivated by the recently proposed alternating single-layer trilayer stacking structure for the nickelate La$_3$Ni$_2$O$_7$, we comprehensively study this system using {\it ab initio} and random-phase approximation techniques. Our analysis unveils similarities between this novel La$_3$Ni$_2$O$_7$ structure and other Ruddlesden-Popper nickelate superconductors, such as a similar charge-transfer gap…
▽ More
Motivated by the recently proposed alternating single-layer trilayer stacking structure for the nickelate La$_3$Ni$_2$O$_7$, we comprehensively study this system using {\it ab initio} and random-phase approximation techniques. Our analysis unveils similarities between this novel La$_3$Ni$_2$O$_7$ structure and other Ruddlesden-Popper nickelate superconductors, such as a similar charge-transfer gap value and orbital-selective behavior of the $e_g$ orbitals. However, different from other Ruddlesden-Popper nickelate superconductors, we do not observe any obvious reconstruction of the Fermi surface from ambient conditions (Cmmm phase) to high pressures (P4/mmm phase). Pressure primarily increases the bandwidths of the Ni $e_g$ bands, suggesting an enhancement of the itinerant properties of those $e_g$ states. Furthermore, the $d_{3z^2-r^2}$ orbital also has a layer-selective behavior because the antibonding-bonding-nonbonding splitting can only be obtained in the trilayer. In addition, we observe a "self-doping" effect from the trilayer to the single-layer sublattices and this effect will be enhanced by overall electron doping. Moreover, we find a leading $d_{x^2-y^2}$-wave pairing state that is restricted to the single-layer. Because the effective coupling between the single layers is very weak -- due to the non-superconducting trilayer in between -- this suggests that the superconducting transition temperature $T_c$ in this structure should be much lower than in the bilayer structure.
△ Less
Submitted 14 August, 2024; v1 submitted 25 April, 2024;
originally announced April 2024.
-
Electrolyte effects on the alkaline hydrogen evolution reaction: a mean-field approach
Authors:
Lucas B. T. de Kam,
Thomas L. Maier,
Katharina Krischer
Abstract:
This paper introduces the combination of an advanced double-layer model with electrochemical kinetics to explain electrolyte effects on the alkaline hydrogen evolution reaction. It is known from experimental studies that the alkaline hydrogen evolution current shows a strong dependence on the concentration and identity of cations in the electrolyte, but is independent of pH. To explain these effec…
▽ More
This paper introduces the combination of an advanced double-layer model with electrochemical kinetics to explain electrolyte effects on the alkaline hydrogen evolution reaction. It is known from experimental studies that the alkaline hydrogen evolution current shows a strong dependence on the concentration and identity of cations in the electrolyte, but is independent of pH. To explain these effects, we formulate the faradaic current in terms of the electric potential in the double layer, which is calculated using a mean-field model that takes into account the cation and anion sizes as well as the electric dipole moment of water molecules. We propose that the Volmer step consists of two activated processes: a water reduction sub-step, and a sub-step in which a hydroxide ion is transferred away from the reaction plane through the double layer. Either of these sub-steps may limit the rate. The proposed models for these sub-steps qualitatively explain experimental observations, including cation effects, pH-independence, and the trend reversal between gold and platinum electrodes. We also assess the quantitative accuracy of the water-reduction-limited current model; we suggest that the predicted functional relationship is valid as long as the hydrogen bonding structure of water near the electrode is sufficiently maintained.
△ Less
Submitted 3 June, 2024; v1 submitted 27 March, 2024;
originally announced March 2024.
-
Prediction of $s^\pm$-wave superconductivity enhanced by electronic doping in trilayer nickelates La$_4$Ni$_3$O$_{10}$ under pressure
Authors:
Yang Zhang,
Ling-Fang Lin,
Adriana Moreo,
Thomas A. Maier,
Elbio Dagotto
Abstract:
Motivated by the recently reported signatures of superconductivity in trilayer La$_4$Ni$_3$O$_{10}$ under pressure, we comprehensively study this system using {\it ab initio} and random-phase approximation techniques. Without electronic interactions, the Ni $d_{3z^2-r^2}$ orbitals show a bonding-antibonding and nonbonding splitting behavior via the O $p_z$ orbitals inducing a ``trimer'' lattice in…
▽ More
Motivated by the recently reported signatures of superconductivity in trilayer La$_4$Ni$_3$O$_{10}$ under pressure, we comprehensively study this system using {\it ab initio} and random-phase approximation techniques. Without electronic interactions, the Ni $d_{3z^2-r^2}$ orbitals show a bonding-antibonding and nonbonding splitting behavior via the O $p_z$ orbitals inducing a ``trimer'' lattice in La$_4$Ni$_3$O$_{10}$, analogous to the dimers of La$_3$Ni$_2$O$_{7}$. The Fermi surface consists of three electron sheets with mixed $e_g$ orbitals, and a hole and an electron pocket made up of the $d_{3z^2-r^2}$ orbital, suggesting a Ni two-orbital minimum model. In addition, we find that superconducting pairing is induced in the $s^{\pm}$-wave channel due to partial nesting between the {\bf M}=$(π, π)$ centered pockets and portions of the Fermi surface centered at the {\bf $Γ$}=$(0, 0)$ point. With changing electronic density $n$, the $s^\pm$ instability remains leading and its pairing strength shows a dome-like behavior with a maximum around $n = 4.2$ ($\sim 6.7\%$ electron doping). The superconducting instability disappears at the same electronic density as that in the new 1313 stacking La$_3$Ni$_2$O$_7$, correlated with the vanishing of the hole pocket that arises from the trilayer sublattice, suggesting that the high-$T_c$ superconductivity of La$_3$Ni$_2$O$_7$ may $not$ originate from a trilayer- and single-layer structure. Furthermore, we predict an interesting spin-density-wave state in La$_4$Ni$_3$O$_{10}$ with an in-plane ($π$, $π$) order and antiferromagnetic coupling between the top and bottom Ni layers, while the middle layer has spin zero.
△ Less
Submitted 23 September, 2024; v1 submitted 7 February, 2024;
originally announced February 2024.
-
Topological Phase Transition without Single-Particle-Gap Closing in Strongly Correlated Systems
Authors:
Peizhi Mai,
Jinchao Zhao,
Thomas A. Maier,
Barry Bradlyn,
Philip W. Phillips
Abstract:
We show here that numerous examples abound where changing topology does not necessarily close the bulk insulating charge gap as demanded in the standard non-interacting picture. From extensive determinantal and dynamical cluster quantum Monte Carlo simulations of the half-filled and quarter-filled Kane-Mele-Hubbard model, we show that for sufficiently strong interactions at either half- or quarter…
▽ More
We show here that numerous examples abound where changing topology does not necessarily close the bulk insulating charge gap as demanded in the standard non-interacting picture. From extensive determinantal and dynamical cluster quantum Monte Carlo simulations of the half-filled and quarter-filled Kane-Mele-Hubbard model, we show that for sufficiently strong interactions at either half- or quarter-filling, a transition between topological and trivial insulators occurs without the closing of a charge gap. To shed light on this behavior, we illustrate that an exactly solvable model reveals that while the single-particle gap remains, the many-body gap does in fact close. These two gaps are the same in the non-interacting system but depart from each other as the interaction turns on. We purport that for interacting systems, the proper probe of topological phase transitions is the closing of the many-body rather than the single-particle gap.
△ Less
Submitted 10 January, 2024; v1 submitted 2 January, 2024;
originally announced January 2024.
-
Electronic structure, magnetic correlations, and superconducting pairing in the reduced Ruddlesden-Popper bilayer La$_3$Ni$_2$O$_6$ under pressure: different role of $d_{3z^2-r^2}$ orbital compared with La$_3$Ni$_2$O$_7$
Authors:
Yang Zhang,
Ling-Fang Lin,
Adriana Moreo,
Thomas A. Maier,
Elbio Dagotto
Abstract:
The recent discovery of superconductivity in bilayer La$_3$Ni$_2$O$_7$ (327-LNO) under pressure stimulated much interest in layered nickelates. However, superconductivity was not found in another bilayer nickelate system, La$_3$Ni$_2$O$_6$ (326-LNO), even under pressure. Using density functional theory and the random phase approximation (RPA), we systematically investigate 326-LNO under pressure.…
▽ More
The recent discovery of superconductivity in bilayer La$_3$Ni$_2$O$_7$ (327-LNO) under pressure stimulated much interest in layered nickelates. However, superconductivity was not found in another bilayer nickelate system, La$_3$Ni$_2$O$_6$ (326-LNO), even under pressure. Using density functional theory and the random phase approximation (RPA), we systematically investigate 326-LNO under pressure. The large crystal-field splitting between the $e_g$ orbitals caused by the missing apical oxygen moves the $d_{3z^2-r^2}$ orbital farther away from the Fermi level, implying that the $d_{3z^2-r^2}$ orbital plays a less important role in 326-LNO than in 327-LNO. This also results in a smaller bandwidth for the $d_{x^2-y^2}$ orbital and a reduced energy gap for the bonding-antibonding splitting of the $d_{3z^2-r^2}$ orbital in 326-LNO, as compared to 327-LNO. Moreover, the in-plane hybridization between the $d_{x^2-y^2}$ and $d_{3z^2-r^2}$ orbitals is found to be small in 326-LNO, while it is much stronger in 327-LNO. The weak inter-layer coupling suggests that $s_{\pm}$-wave pairing is unlikely in 326-LNO. The robust in-plane ferromagnetic coupling also suggests that d-wave superconductivity, which is usually caused by antiferromagnetic fluctuations of the $d_{x^2-y^2}$ orbital, is also unlikely in 326-LNO. These conclusions are supported by our many-body RPA calculations of the pairing behavior. Contrasting with the cuprates, for the bilayer cuprate HgBa$_2$CaCu$_2$O$_6$, we find a strong "self-doping effect" of the $d_{x^2-y^2}$ orbital under pressure, with the charge of Cu being reduced by approximately 0.13 electrons from 0 GPa to 25 GPa. In contrast, we do not observe such a change in the electronic density in 326-LNO under pressure, establishing another important difference between the nickelates and the cuprates.
△ Less
Submitted 26 October, 2023;
originally announced October 2023.
-
How Metal/Insulator Interfaces Enable the Enhancement of the Hydrogen Evolution Reaction Kinetics in Two Ways
Authors:
Thomas L. Maier,
Lucas B. T. de Kam,
Matthias Golibrzuch,
Tina Angerer,
Markus Becherer,
Katharina Krischer
Abstract:
Laterally nanostructured surfaces give rise to a new dimension of understanding and improving electrochemical reactions. In this study, we present a peculiar mechanism appearing at a metal/insulator interface, which can significantly enhance the Hydrogen Evolution Reaction (HER) from water reduction by altering the local reaction conditions in two ways: facilitated adsorption of hydrogen on the me…
▽ More
Laterally nanostructured surfaces give rise to a new dimension of understanding and improving electrochemical reactions. In this study, we present a peculiar mechanism appearing at a metal/insulator interface, which can significantly enhance the Hydrogen Evolution Reaction (HER) from water reduction by altering the local reaction conditions in two ways: facilitated adsorption of hydrogen on the metal catalyst surface and improved transfer of ions through the double layer. The mechanism is uncovered using electrodes consisting of well-defined nanometer-sized metal arrays (Au, Cu, Pt) embedded in an insulator layer (silicon nitride), varying various parameters of both the electrode (size of the metal patches, catalyst material) and the electrolyte (cationic species, cation concentration, pH). In addition, simulations of the electrochemical double layer are carried out, which support the elaborated mechanism. Knowledge of this mechanism will enable new design principles for novel composite electrocatalytic systems.
△ Less
Submitted 5 September, 2023;
originally announced September 2023.
-
Trends in electronic structures and $s_{\pm}$-wave pairing for the rare-earth series in bilayer nickelate superconductor $R_ 3$Ni$_2$O$_7$
Authors:
Yang Zhang,
Ling-Fang Lin,
Adriana Moreo,
Thomas A. Maier,
Elbio Dagotto
Abstract:
The recent discovery of pressure-induced superconductivity in the bilayer La$_3$Ni$_2$O$_7$ (LNO) has opened a new platform for the study of unconventional superconductors. In this publication, we investigate theoretically the whole family of bilayer 327-type nickelates $R_3$Ni$_2$O$_7$ ($R$ = Rare-earth elements) under pressure. From La to Lu, the lattice constants and volume decrease, leading to…
▽ More
The recent discovery of pressure-induced superconductivity in the bilayer La$_3$Ni$_2$O$_7$ (LNO) has opened a new platform for the study of unconventional superconductors. In this publication, we investigate theoretically the whole family of bilayer 327-type nickelates $R_3$Ni$_2$O$_7$ ($R$ = Rare-earth elements) under pressure. From La to Lu, the lattice constants and volume decrease, leading to enhanced in-plane and out-of-plane hoppings, resulting in an effectively reduced electronic correlation $U/W$. Furthermore, the Ni's $t_{2g}$ states shift away from the $e_g$ states, while the crystal-field splitting between $d_{3z^2-r^2}$ and $d_{x^2-y^2}$ is almost unchanged. In addition, six candidates were found to become stable in the Fmmm phase, with increasing values of critical pressure as the atomic number increased. Similar to the case of LNO, the $s_{\pm}$-wave pairing tendency dominates in all candidates, due to the nesting between the {\bf M}=$(π,π)$ and the {\bf X}=$(π,0)$ and {\bf Y}=$(0,π)$ points in the Brillouin zone. Then, $T_c$ is expected to decrease as the radius of rare-earth (RE) ions decreases. Our results suggest that LNO is already the "optimal" candidate, with Ce a close competitor, among the whole RE bilayer nickelates, and to increase $T_c$ we suggest to grow on special substrates with larger in-plane lattice spacings.
△ Less
Submitted 26 October, 2023; v1 submitted 14 August, 2023;
originally announced August 2023.
-
Structural phase transition, $s_{\pm}$-wave pairing and magnetic stripe order in the bilayered nickelate superconductor La$_3$Ni$_2$O$_7$ under pressure
Authors:
Yang Zhang,
Ling-Fang Lin,
Adriana Moreo,
Thomas A. Maier,
Elbio Dagotto
Abstract:
Motivated by the recently discovered high-$T_c$ superconductor La$_3$Ni$_2$O$_7$, we comprehensively study this system using density functional theory and random phase approximation calculations. At low pressures, the Amam phase is stable, containing the Y$^{2-}$ mode distortion from the Fmmm phase, while the Fmmm phase is unstable. Because of small differences in enthalpy and a considerable Y…
▽ More
Motivated by the recently discovered high-$T_c$ superconductor La$_3$Ni$_2$O$_7$, we comprehensively study this system using density functional theory and random phase approximation calculations. At low pressures, the Amam phase is stable, containing the Y$^{2-}$ mode distortion from the Fmmm phase, while the Fmmm phase is unstable. Because of small differences in enthalpy and a considerable Y$^{2-}$ mode amplitude, the two phases may coexist in the range between 10.6 and 14 GPa, beyond which the Fmmm phase dominates. In addition, the magnetic stripe-type spin order with wavevector ($π$, 0) was stable at the intermediate region. Pairing is induced in the $s_{\pm}$-wave channel due to partial nesting between the {\bf M}=$(π, π)$ centered pockets and portions of the Fermi surface centered at the {\bf X}=$(π, 0)$ and {\bf Y}=$(0, π)$ points. This resembles results for iron-based superconductors but has a fundamental difference with iron pnictides and selenides. Moreover, our present efforts also suggest that La$_3$Ni$_2$O$_7$ is qualitatively different from infinite-layer nickelates and cuprate superconductors.
△ Less
Submitted 25 March, 2024; v1 submitted 27 July, 2023;
originally announced July 2023.
-
Weak-Coupling Theory of Neutron Scattering as a Probe of Altermagnetism
Authors:
Thomas A. Maier,
Satoshi Okamoto
Abstract:
Inelastic neutron scattering provides a powerful probe of the magnetic excitations of quantum magnets. Altermagnets have recently emerged as a new class of magnets with vanishing net magnetization characteristic of antiferromagnets and with a spin-split electronic structure typical of ferromagnets. Here we introduce a minimal Hubbard model with two-sublattice orthorhombic anisotropy as a framework…
▽ More
Inelastic neutron scattering provides a powerful probe of the magnetic excitations of quantum magnets. Altermagnets have recently emerged as a new class of magnets with vanishing net magnetization characteristic of antiferromagnets and with a spin-split electronic structure typical of ferromagnets. Here we introduce a minimal Hubbard model with two-sublattice orthorhombic anisotropy as a framework to study altermagnetism. Using unrestricted Hartree-Fock calculations, we find an altermagnetic state for this model that evolves from a metallic state to an insulating state with increasing Hubbard-$U$ Coulomb repulsion. We then examine the inelastic neutron scattering response in these states using random-phase approximation calculations of the dynamic spin susceptibility $χ''({\bf q}, ω)$. We find that the magnetic excitation spectrum depends on its chirality for ${\bf q}$ along certain directions in reciprocal space, an observation that may be used in inelastic neutron scattering experiments as a probe of altermagnetism.
△ Less
Submitted 28 September, 2023; v1 submitted 7 July, 2023;
originally announced July 2023.
-
Magnetism and Metallicity in Moiré Transition Metal Dichalcogenides
Authors:
Patrick Tscheppe,
Jiawei Zang,
Marcel Klett,
Seher Karakuzu,
Armelle Celarier,
Zhengqian Cheng,
Chris A. Marianetti,
Thomas A. Maier,
Michel Ferrero,
Andrew J. Millis,
Thomas Schäfer
Abstract:
The ability to control the properties of twisted bilayer transition metal dichalcogenides in situ makes them an ideal platform for investigating the interplay of strong correlations and geometric frustration. Of particular interest are the low energy scales, which make it possible to experimentally access both temperature and magnetic fields that are of the order of the bandwidth or the correlatio…
▽ More
The ability to control the properties of twisted bilayer transition metal dichalcogenides in situ makes them an ideal platform for investigating the interplay of strong correlations and geometric frustration. Of particular interest are the low energy scales, which make it possible to experimentally access both temperature and magnetic fields that are of the order of the bandwidth or the correlation scale. In this manuscript we analyze the moiré Hubbard model, believed to describe the low energy physics of an important subclass of the twisted bilayer compounds. We establish its magnetic and the metal-insulator phase diagram for the full range of magnetic fields up to the fully spin polarized state. We find a rich phase diagram including fully and partially polarized insulating and metallic phases of which we determine the interplay of magnetic order, Zeeman-field, and metallicity, and make connection to recent experiments.
△ Less
Submitted 27 April, 2023; v1 submitted 23 March, 2023;
originally announced March 2023.
-
PlasmoFAB: A Benchmark to Foster Machine Learning for Plasmodium falciparum Protein Antigen Candidate Prediction
Authors:
Jonas Christian Ditz,
Jacqueline Wistuba-Hamprecht,
Timo Maier,
Rolf Fendel,
Nico Pfeifer,
Bernhard Reuter
Abstract:
Motivation: Machine learning methods can be used to support scientific discovery in healthcare-related research fields. However, these methods can only be reliably used if they can be trained on high-quality and curated datasets. Currently, no such dataset for the exploration of Plasmodium falciparum protein antigen candidates exists. The parasite Plasmodium falciparum causes the infectious diseas…
▽ More
Motivation: Machine learning methods can be used to support scientific discovery in healthcare-related research fields. However, these methods can only be reliably used if they can be trained on high-quality and curated datasets. Currently, no such dataset for the exploration of Plasmodium falciparum protein antigen candidates exists. The parasite Plasmodium falciparum causes the infectious disease malaria. Thus, identifying potential antigens is of utmost importance for the development of antimalarial drugs and vaccines. Since exploring antigen candidates experimentally is an expensive and time-consuming process, applying machine learning methods to support this process has the potential to accelerate the development of drugs and vaccines, which are needed for fighting and controlling malaria. Results: We developed PlasmoFAB, a curated benchmark that can be used to train machine learning methods for the exploration of Plasmodium falciparum protein antigen candidates. We combined an extensive literature search with domain expertise to create high-quality labels for Plasmodium falciparum specific proteins that distinguish between antigen candidates and intracellular proteins. Additionally, we used our benchmark to compare different well-known prediction models and available protein localization prediction services on the task of identifying protein antigen candidates. We show that available general-purpose services are unable to provide sufficient performance on identifying protein antigen candidates and are outperformed by our models that were trained on this tailored data. Availability: PlasmoFAB is publicly available on Zenodo with DOI 10.5281/zenodo.7433087. Furthermore, all scripts that were used in the creation of PlasmoFAB and the training and evaluation of machine learning models are open source and publicly available on GitHub here: https://github.com/msmdev/PlasmoFAB.
△ Less
Submitted 3 May, 2023; v1 submitted 16 January, 2023;
originally announced January 2023.
-
Spin singlet topological superconductivity in the attractive Rashba Hubbard model
Authors:
Peter Doak,
Giovanni Balduzzi,
Pontus Laurell,
Elbio Dagotto,
Thomas A. Maier
Abstract:
Fully gapped, spin singlet superconductors with antisymmetric spin-orbit coupling in a Zeeman magnetic field provide a promising route to realize superconducting states with non-Abelian topological order and therefore fault-tolerant quantum computation. Here we use a quantum Monte Carlo dynamical cluster approximation to study the superconducting properties of a doped two-dimensional attractive Hu…
▽ More
Fully gapped, spin singlet superconductors with antisymmetric spin-orbit coupling in a Zeeman magnetic field provide a promising route to realize superconducting states with non-Abelian topological order and therefore fault-tolerant quantum computation. Here we use a quantum Monte Carlo dynamical cluster approximation to study the superconducting properties of a doped two-dimensional attractive Hubbard model with Rashba spin-orbit coupling in a Zeeman magnetic field. We generally find that the Rashba coupling has a beneficial effect towards s-wave superconductivity. In the presence of a finite Zeeman field, when superconductivity is suppressed by Pauli pair-breaking, the Rashba coupling counteracts the spin imbalance created by the Zeeman field by mixing the spins, and thus restores superconductivity at finite temperatures. We show that this favorable effect of the spin-orbit coupling is traced to a spin-flip driven enhancement of the amplitude for the propagation of a pair of electrons in time-reversed states. Moreover, by inspecting the Fermi surface of the interacting model, we show that for sufficiently large Rashba coupling and Zeeman field, the superconducting state is expected to be topologically non-trivial.
△ Less
Submitted 14 December, 2022;
originally announced December 2022.
-
Robust charge-density wave correlations in the electron-doped single-band Hubbard model
Authors:
Peizhi Mai,
Nathan S. Nichols,
Seher Karakuzu,
Feng Bao,
Adrian Del Maestro,
Thomas A. Maier,
Steven Johnston
Abstract:
There is growing evidence that the hole-doped single-band Hubbard and $t$-$J$ models do not have a superconducting ground state reflective of the high-temperature cuprate superconductors but instead have striped spin- and charge-ordered ground states. Nevertheless, it is proposed that these models may still provide an effective low-energy model for electron-doped materials. Here we study the finit…
▽ More
There is growing evidence that the hole-doped single-band Hubbard and $t$-$J$ models do not have a superconducting ground state reflective of the high-temperature cuprate superconductors but instead have striped spin- and charge-ordered ground states. Nevertheless, it is proposed that these models may still provide an effective low-energy model for electron-doped materials. Here we study the finite temperature spin and charge correlations in the electron-doped Hubbard model using quantum Monte Carlo dynamical cluster approximation calculations and contrast their behavior with those found on the hole-doped side of the phase diagram. We find evidence for a charge modulation with both checkerboard and unidirectional components decoupled from any spin-density modulations. These correlations are inconsistent with a weak-coupling description based on Fermi surface nesting, and their doping dependence agrees qualitatively with resonant inelastic x-ray scattering measurements. Our results provide evidence that the single-band Hubbard model describes the electron-doped cuprates.
△ Less
Submitted 26 October, 2022;
originally announced October 2022.
-
Evidence for chiral superconductivity on a silicon surface
Authors:
F. Ming,
X. Wu,
C. Chen,
K. D. Wang,
P. Mai,
T. A. Maier,
J. Strockoz,
J. W. F. Venderbos,
C. Gonzalez,
J. Ortega,
S. Johnston,
H. H. Weitering
Abstract:
Sn adatoms on a Si(111) substrate with 1/3 monolayer coverage form a two-dimensional triangular adatom lattice with one unpaired electron per site and an antiferromagnetic Mott insulating state. The Sn layers can be modulation hole-doped and metallized using heavily-doped $p$-type Si(111) substrates, and become superconducting at low temperatures. While the pairing symmetry of the superconducting…
▽ More
Sn adatoms on a Si(111) substrate with 1/3 monolayer coverage form a two-dimensional triangular adatom lattice with one unpaired electron per site and an antiferromagnetic Mott insulating state. The Sn layers can be modulation hole-doped and metallized using heavily-doped $p$-type Si(111) substrates, and become superconducting at low temperatures. While the pairing symmetry of the superconducting state is currently unknown, the combination of repulsive interactions and frustration inherent to the triangular adatom lattice opens up the possibility for a chiral order parameter. Here, we study the superconducting state of Sn/Si(111) using scanning tunneling microscopy/spectroscopy and quasi-particle interference imaging. We find evidence for a doping-dependent $T_c$ with a fully gapped order parameter, the presence of time-reversal symmetry breaking, and a strong enhancement of the zero-bias conductance near the edges of the superconducting domains. While each individual piece of evidence could have a more mundane interpretation, our combined results suggest the tantalizing possibility that Sn/Si(111) is an unconventional chiral d-wave superconductor.
△ Less
Submitted 12 October, 2022;
originally announced October 2022.
-
Stripe correlations in the two-dimensional Hubbard-Holstein model
Authors:
Seher Karakuzu,
Andy Tanjaroon Ly,
Peizhi Mai,
James Neuhaus,
Thomas A. Maier,
Steven Johnston
Abstract:
Several state-of-the-art numerical methods have observed static or fluctuating spin and charge stripes in doped two-dimensional Hubbard models, suggesting that these orders play a significant role in shaping the cuprate phase diagram. Many experiments, however, also indicate that the cuprates have strong electron-phonon ($e$-ph) coupling, and it is unclear how this interaction influences stripe co…
▽ More
Several state-of-the-art numerical methods have observed static or fluctuating spin and charge stripes in doped two-dimensional Hubbard models, suggesting that these orders play a significant role in shaping the cuprate phase diagram. Many experiments, however, also indicate that the cuprates have strong electron-phonon ($e$-ph) coupling, and it is unclear how this interaction influences stripe correlations. We study static and fluctuating stripe orders in the doped singleband Hubbard-Holstein model using zero temperature variational Monte Carlo and finite temperature determinant quantum Monte Carlo. We find that the lattice couples more strongly with the charge component of the stripes, leading to an enhancement or suppression of stripe correlations, depending on model parameters like the next-nearest-neighbor hopping $t^\prime$ or phonon energy $Ω$. Our results help elucidate how the $e$-ph interaction can tip the delicate balance between stripe and superconducting correlations in the Hubbard-Holstein model with implications for our understanding of the high-$T_\mathrm{c}$ cuprates.
△ Less
Submitted 30 May, 2022;
originally announced May 2022.
-
Leading superconducting instabilities in three-dimensional models for Sr2RuO4
Authors:
Astrid T. Rømer,
T. A. Maier,
Andreas Kreisel,
P. J. Hirschfeld,
Brian M. Andersen
Abstract:
The superconductor Sr2RuO4 has been the subject of enormous interest over more than two decades, but until now the form of its order parameter has not been determined. Since groundbreaking NMR experiments revealed that the pairs are of dominant spin-singlet character, attention has focused on time-reversal symmetry breaking linear combinations of $s$-, $d$- and $g$-wave one-dimensional (1D) irredu…
▽ More
The superconductor Sr2RuO4 has been the subject of enormous interest over more than two decades, but until now the form of its order parameter has not been determined. Since groundbreaking NMR experiments revealed that the pairs are of dominant spin-singlet character, attention has focused on time-reversal symmetry breaking linear combinations of $s$-, $d$- and $g$-wave one-dimensional (1D) irreducible representations. However, a state of the form $d_{xz}+id_{yz}$ has also been proposed. We present a systematic study of the stability of various superconducting candidate states, assuming that pairing is driven by the fluctuation exchange mechanism, including a realistic three-dimensional Fermi surface, full treatment of both local and non-local spin-orbit couplings, and a wide range of interaction parameters $U,J,U',J'$. The leading superconducting instabilities are found to exhibit nodal even-parity $A_{1g} (s')$ or $B_{1g} (d_{x^2-y^2})$ symmetries, similar to the findings in two-dimensional models without longer-range Coulomb interaction which tends to favor $d_{xy}$ over $d_{x^2-y^2}$. Within the so-called Hund's coupling mean-field pairing scenario, the $E_g (d_{xz}/d_{yz})$ solution can be stabilized for large $J$ and specific forms of the spin-orbit coupling, but for all cases studied here the eigenvalues of other superconducting solutions are significantly larger when the full fluctuation exchange vertex is included in the pairing kernel. Additionally, we compute the spin susceptibility in relevant superconducting candidate phases and compare to recent neutron scattering and NMR Knight shift measurements.
△ Less
Submitted 6 April, 2022;
originally announced April 2022.
-
Enhancing $T_{\mathrm{c}}$ in a composite superconductor/metal bilayer system: a dynamical cluster approximation study
Authors:
Philip M. Dee,
Steven Johnston,
Thomas Maier
Abstract:
It has been proposed that the superconducting transition temperature $T_{\mathrm{c}}$ of an unconventional superconductor with a large pairing scale but strong phase fluctuations can be enhanced by coupling it to a metal. However, the general efficacy of this approach across different parameter regimes remains an open question. Using the dynamical cluster approximation, we study this question in a…
▽ More
It has been proposed that the superconducting transition temperature $T_{\mathrm{c}}$ of an unconventional superconductor with a large pairing scale but strong phase fluctuations can be enhanced by coupling it to a metal. However, the general efficacy of this approach across different parameter regimes remains an open question. Using the dynamical cluster approximation, we study this question in a system composed of an attractive Hubbard layer in the intermediate coupling regime, where the magnitude of the attractive Coulomb interaction $|U|$ is slightly larger than the bandwidth $W$, hybridized with a noninteracting metallic layer. We find that while the superconducting transition becomes more mean-field-like with increasing interlayer hopping, the superconducting transition temperature $T_{\mathrm{c}}$ exhibits a nonmonotonic dependence on the strength of the hybridization $t_{\perp}$. This behavior arises from a reduction of the effective pairing interaction in the correlated layer that out-competes the growth in the intrinsic pair-field susceptibility induced by the coupling to the metallic layer. We find that the largest $T_{\mathrm{c}}$ inferred here for the composite system is below the maximum value currently estimated for the isolated negative-$U$ Hubbard model.
△ Less
Submitted 10 March, 2022;
originally announced March 2022.
-
Strongly anisotropic electronic and magnetic structures in oxide dichlorides RuOCl$_2$ and OsOCl$_2$
Authors:
Yang Zhang,
Ling-Fang Lin,
Adriana Moreo,
Thomas A. Maier,
Gonzalo Alvarez,
Elbio Dagotto
Abstract:
Here, using density functional theory and density matrix renormalization group methods, we investigate the electronic and magnetic properties of RuOCl$_2$ and OsOCl$_2$ with $d^4$ electronic configurations. Different from a previous study using VOI$_2$ with $d^1$ configuration, these systems with $4d^4$ or $5d^4$ do not exhibit a ferroelectric instability along the $a$-axis. Due to the fully-occup…
▽ More
Here, using density functional theory and density matrix renormalization group methods, we investigate the electronic and magnetic properties of RuOCl$_2$ and OsOCl$_2$ with $d^4$ electronic configurations. Different from a previous study using VOI$_2$ with $d^1$ configuration, these systems with $4d^4$ or $5d^4$ do not exhibit a ferroelectric instability along the $a$-axis. Due to the fully-occupied $d_{xy}$ orbital in RuOCl$_2$ and OsOCl$_2$, the Peierls instability distortion disappears along the $b$-axis, leading to an undistorted I${\rm mmm}$ phase (No. 71). Furthermore, we observe strongly anisotropic electronic and magnetic structures along the $a$-axis. The large crystal-field splitting energy (between $d_{xz/yz}$ and $d_{xy}$ orbitals) and large hopping between nearest-neighbor Ru and Os atoms suppresses the spin-orbital effect in $M$OCl$_2$ ($M$ = Ru or Os) with electronic density $n = 4$, resulting in a spin-1 system instead of a $J = 0$ singlet ground state. Moreover, we find staggered antiferromagnetic order with $π$ wavevector along the $M$-O chain direction ($a$-axis) while the magnetic coupling along the $b$-axis is weak. Based on Wannier functions from first-principles calculations, we calculated the relevant hopping amplitudes and crystal-field splitting energies of the $t_{2g}$ orbitals for the Os atoms to construct a multi-orbital Hubbard model for the $M$-O chains. Staggered AFM with $\uparrow$-$\downarrow$-$\uparrow$-$\downarrow$ spin structure dominates in our DMRG calculations, in agreement with DFT calculations.
△ Less
Submitted 14 April, 2022; v1 submitted 28 February, 2022;
originally announced March 2022.
-
Fluctuating intertwined stripes in the strange metal regime of the Hubbard model
Authors:
Edwin W. Huang,
Tianyi Liu,
Wen O. Wang,
Hong-Chen Jiang,
Peizhi Mai,
Thomas A. Maier,
Steven Johnston,
Brian Moritz,
Thomas P. Devereaux
Abstract:
Strongly correlated electron systems host a variety of poorly understood correlations in their high temperature normal state. Unlike ordered phases defined by order parameters, these normal state phases are often defined through unconventional properties such as strange metallic transport or spectroscopic pseudogaps. Characterizing the microscopic correlations in the normal state is necessary to e…
▽ More
Strongly correlated electron systems host a variety of poorly understood correlations in their high temperature normal state. Unlike ordered phases defined by order parameters, these normal state phases are often defined through unconventional properties such as strange metallic transport or spectroscopic pseudogaps. Characterizing the microscopic correlations in the normal state is necessary to elucidate mechanisms that lead to these properties and their connection to ground state orders. Here we establish the presence of intertwined charge and spin stripes in the strange metal normal state of the Hubbard model using determinant quantum Monte Carlo calculations. The charge and spin density waves constituting the stripes are fluctuating and short-ranged, yet they obey a mutual commensurability relation and remain microscopically interlocked, as evidenced through measurements of three-point spin-spin-hole correlation functions. Our findings demonstrate the ability of many-body numerical simulations to unravel the microscopic correlations that define quantum states of matter.
△ Less
Submitted 17 February, 2022;
originally announced February 2022.
-
Coupled Hubbard ladders at weak coupling: Pairing and spin excitations
Authors:
Thomas A. Maier,
Elbio Dagotto
Abstract:
The Hubbard model provides a simple framework in which one can study how certain aspects of the electronic structure of strongly interacting systems can be tuned to optimize the superconducting pairing correlations and how these changes affect the mechanisms giving rise to them. Here we use a weak-coupling random phase approximation to study a two-dimensional Hubbard model with a unidirectional mo…
▽ More
The Hubbard model provides a simple framework in which one can study how certain aspects of the electronic structure of strongly interacting systems can be tuned to optimize the superconducting pairing correlations and how these changes affect the mechanisms giving rise to them. Here we use a weak-coupling random phase approximation to study a two-dimensional Hubbard model with a unidirectional modulation of the hopping amplitudes as the system evolves from the uniform square lattice to an array of weakly coupled two-leg ladders. We find that the pairing correlations retain their dominant $d_{x^2-y^2}$-wave like structure and that they are significantly enhanced for a slightly modulated lattice. This enhancement is traced backed to an increase in the strength of the spin-fluctuation pairing interacting due to favorable Fermi surface nesting in the modulated system. We then use a random-phase approximation BCS framework to examine the evolution of the neutron resonance in the superconducting state. We find that it changes only weakly for moderate modulations, but breaks up into two distinct resonances at incommensurate wave-vectors in the limit of weakly coupled ladders.
△ Less
Submitted 13 December, 2021;
originally announced December 2021.
-
Superconductivity in the bilayer Hubbard model: Are two Fermi surfaces better than one?
Authors:
Seher Karakuzu,
Steven Johnston,
Thomas A. Maier
Abstract:
Fully occupied or unoccupied bands in a solid are often considered inert and irrelevant to a material's low-energy properties. But the discovery of enhanced superconductivity in heavily electron-doped FeSe-derived superconductors poses questions about the possible role of incipient bands (those laying close to but not crossing the Fermi level) in pairing. To answer this question, researchers have…
▽ More
Fully occupied or unoccupied bands in a solid are often considered inert and irrelevant to a material's low-energy properties. But the discovery of enhanced superconductivity in heavily electron-doped FeSe-derived superconductors poses questions about the possible role of incipient bands (those laying close to but not crossing the Fermi level) in pairing. To answer this question, researchers have studied pairing correlations in the bilayer Hubbard model, which has an incipient band for large interlayer hopping $t_\perp$, using many-body perturbation theory and variational methods. They have generally found that superconductivity is enhanced as one of the bands approaches the Liftshiz transition and even when it becomes incipient. Here, we address this question using the nonperturbative quantum Monte Carlo (QMC) dynamical cluster approximation (DCA) to study the bilayer Hubbard model's pairing correlations. We find that the model has robust $s_\pm$ pairing correlations in the large $t_\perp$ limit, which can become stronger as one band is made incipient. While this behavior is linked to changes in the effective interaction, we further find that it is counteracted by a suppression of the intrinsic pair-field susceptibility and does not translate to an increased $T_c$. Our results demonstrate that the highest achievable transition temperatures in the bilayer Hubbard model occur when the system has two bands crossing the Fermi level.
△ Less
Submitted 29 July, 2021;
originally announced July 2021.
-
Intertwined spin, charge and pair correlations in the two-dimensional Hubbard model in the thermodynamic limit
Authors:
Peizhi Mai,
Seher Karakuzu,
Giovanni Balduzzi,
Steven Johnston,
Thomas A. Maier
Abstract:
The high-temperature superconducting cuprates are governed by intertwined striped magnetic and charge orders in addition to superconductivity. Remarkably similar behavior has also been seen in numerical calculations for the Hubbard model describing the copper-oxygen layers in these materials. Finite cluster methods typically find that spin and charge stripe order dominates, while embedded quantum…
▽ More
The high-temperature superconducting cuprates are governed by intertwined striped magnetic and charge orders in addition to superconductivity. Remarkably similar behavior has also been seen in numerical calculations for the Hubbard model describing the copper-oxygen layers in these materials. Finite cluster methods typically find that spin and charge stripe order dominates, while embedded quantum cluster methods, which access the thermodynamic limit, often conclude that superconductivity does. Here, we report the observation of fluctuating spin and charge stripes in an embedded cluster calculation for the Hubbard model. This discovery demonstrates that striped states survive in the thermodynamic limit and allows us to study their influence on the model's superconducting properties, where we find evidence for pair-density-wave correlations intertwined with the stripe correlations.
△ Less
Submitted 9 February, 2022; v1 submitted 3 June, 2021;
originally announced June 2021.
-
Dynamical Cluster Approximation Study of Electron Localization in the Extended Hubbard Model
Authors:
Hanna Terletska,
Sergei Iskakov,
Thomas Maier,
Emanuel Gull
Abstract:
We perform a detailed study of the phase transitions and mechanisms of electron localization in the extended Hubbard model using the dynamical cluster approximation on a $2\times 2$ cluster. We explore the interplay of charge order and Mott physics. We find that a nearest-neighbor Coulomb interaction $V$ causes "screening" effects close to the Mott phase transition, pushing the phase boundary to l…
▽ More
We perform a detailed study of the phase transitions and mechanisms of electron localization in the extended Hubbard model using the dynamical cluster approximation on a $2\times 2$ cluster. We explore the interplay of charge order and Mott physics. We find that a nearest-neighbor Coulomb interaction $V$ causes "screening" effects close to the Mott phase transition, pushing the phase boundary to larger values of $U$. We also demonstrate the different effects of $V$ on correlations in metallic and insulating regimes and document the different correlation aspects of charge order and Mott states.
△ Less
Submitted 19 May, 2021;
originally announced May 2021.
-
Pairing correlations in the cuprates: a numerical study of the three-band Hubbard model
Authors:
Peizhi Mai,
Giovanni Balduzzi,
Steven Johnston,
Thomas A. Maier
Abstract:
We study the three-band Hubbard model for the copper oxide plane of the high-temperature superconducting cuprates using determinant quantum Monte Carlo and the dynamical cluster approximation (DCA) and provide a comprehensive view of the pairing correlations in this model using these methods. Specifically, we compute the pair-field susceptibility and study its dependence on temperature, doping, in…
▽ More
We study the three-band Hubbard model for the copper oxide plane of the high-temperature superconducting cuprates using determinant quantum Monte Carlo and the dynamical cluster approximation (DCA) and provide a comprehensive view of the pairing correlations in this model using these methods. Specifically, we compute the pair-field susceptibility and study its dependence on temperature, doping, interaction strength, and charge-transfer energy. Using the DCA, we also solve the Bethe-Salpeter equation for the two-particle Green's function in the particle-particle channel to determine the transition temperature to the superconducting phase on smaller clusters. Our calculations reproduce many aspects of the cuprate phase diagram and indicate that there is an "optimal" value of the charge-transfer energy for the model where $T_c$ is largest. These results have implications for our understanding of superconductivity in both the cuprates and other doped charge-transfer insulators.
△ Less
Submitted 9 April, 2021; v1 submitted 26 January, 2021;
originally announced January 2021.
-
Orbital structure of the effective pairing interaction in the high-temperature superconducting cuprates
Authors:
Peizhi Mai,
Giovanni Balduzzi,
Steven Johnston,
Thomas A. Maier
Abstract:
The nature of the effective interaction responsible for pairing in the high-temperature superconducting cuprates remains unsettled. This question has been studied extensively using the simplified single-band Hubbard model, which does not explicitly consider the orbital degrees of freedom of the relevant CuO$_2$ planes. Here, we use a dynamic cluster quantum Monte Carlo approximation to study the o…
▽ More
The nature of the effective interaction responsible for pairing in the high-temperature superconducting cuprates remains unsettled. This question has been studied extensively using the simplified single-band Hubbard model, which does not explicitly consider the orbital degrees of freedom of the relevant CuO$_2$ planes. Here, we use a dynamic cluster quantum Monte Carlo approximation to study the orbital structure of the pairing interaction in the three-band Hubbard model, which treats the orbital degrees of freedom explicitly. We find that the interaction predominately acts between neighboring copper orbitals, but with significant additional weight appearing on the surrounding bonding molecular oxygen orbitals. By explicitly comparing these results to those from the simpler single-band Hubbard model, our study provides strong support for the single-band framework for describing superconductivity in the cuprates.
△ Less
Submitted 29 September, 2020;
originally announced September 2020.
-
Evolution of spin excitations from bulk to monolayer FeSe
Authors:
J. Pelliciari,
S. Karakuzu,
Q. Song,
R. Arpaia,
A. Nag,
M. Rossi,
J. Li,
T. Yu,
X. Chen,
R. Peng,
M. Garcia-Fernandez,
A. C. Walters,
Q. Wang,
J. Zhao,
G. Ghiringhelli,
D. Feng,
T. A. Maier,
K. -J. Zhou,
S. Johnston,
R. Comin
Abstract:
The discovery of enhanced superconductivity (SC) in FeSe films grown on SrTiO3 (FeSe/STO) has revitalized the field of Fe-based superconductors. In the ultrathin limit, the superconducting transition temperature Tc is increased by almost an order of magnitude, raising new questions on the pairing mechanism. As in other unconventional superconductors, antiferromagnetic spin fluctuations have been p…
▽ More
The discovery of enhanced superconductivity (SC) in FeSe films grown on SrTiO3 (FeSe/STO) has revitalized the field of Fe-based superconductors. In the ultrathin limit, the superconducting transition temperature Tc is increased by almost an order of magnitude, raising new questions on the pairing mechanism. As in other unconventional superconductors, antiferromagnetic spin fluctuations have been proposed as a candidate to mediate SC in this system. Thus, it is essential to study the evolution of the spin dynamics of FeSe in the ultrathin limit to elucidate their relationship with superconductivity. Here, we investigate and compare the spin excitations in bulk and monolayer FeSe grown on STO using high-resolution resonant inelastic x-ray scattering (RIXS) and quantum Monte Carlo (QMC) calculations. Despite the absence of long-range magnetic order, bulk FeSe displays dispersive magnetic excitations reminiscent of other Fe-pnictides. Conversely, the spin excitations in FeSe/STO are gapped, dispersionless, and significantly hardened relative to the bulk counterpart. By comparing our RIXS results with simulations of a bilayer Hubbard model, we connect the evolution of the spin excitations to the Fermiology of the two systems. The present study reveals a remarkable reconfiguration of spin excitations in FeSe/STO, which is essential to understand the role of spin fluctuations in the pairing mechanism.
△ Less
Submitted 29 September, 2020; v1 submitted 21 August, 2020;
originally announced August 2020.
-
The overdoped end of the cuprate phase diagram
Authors:
Thomas A. Maier,
Seher Karakuzu,
Douglas J. Scalapino
Abstract:
Studying the disappearance of superconductivity at the end of the overdoped region of the cuprate phase diagram offers a different approach for investigating the interaction which is responsible for pairing in these materials. In the underdoped region this question is complicated by the presence of charge and stripe ordered phases as well as the pseudogap. In the overdoped region the situation app…
▽ More
Studying the disappearance of superconductivity at the end of the overdoped region of the cuprate phase diagram offers a different approach for investigating the interaction which is responsible for pairing in these materials. In the underdoped region this question is complicated by the presence of charge and stripe ordered phases as well as the pseudogap. In the overdoped region the situation appears simpler with only a normal phase, a superconducting phase and impurity scattering. Here, for the overdoped region, we report the results of a combined dynamic cluster approximation (DCA) and a weak Born impurity scattering calculation for a $t-t'-U$ Hubbard model. We find that a decrease in the $d$-wave pairing strength of the two-particle scattering vertex is closely coupled to changes in the momentum and frequency structure of the magnetic spin fluctuations as the system is overdoped. Treating the impurity scattering within a disordered BCS $d$-wave approximation, we see how the combined effects of the decreasing $d$-wave pairing strength and weak impurity scattering lead to the end of the $T_c$ dome.
△ Less
Submitted 27 April, 2020;
originally announced April 2020.
-
Continuous momentum dependence in the dynamical cluster approximation
Authors:
Urs R. Hähner,
Thomas A. Maier,
Thomas C. Schulthess
Abstract:
The dynamical cluster approximation (DCA) is a quantum cluster extension to the single-site dynamical mean-field theory that incorporates spatially nonlocal dynamic correlations systematically and nonperturbatively. The DCA$^+$ algorithm addresses the cluster shape dependence of the DCA and improves the convergence with cluster size by introducing a lattice self-energy with continuous momentum dep…
▽ More
The dynamical cluster approximation (DCA) is a quantum cluster extension to the single-site dynamical mean-field theory that incorporates spatially nonlocal dynamic correlations systematically and nonperturbatively. The DCA$^+$ algorithm addresses the cluster shape dependence of the DCA and improves the convergence with cluster size by introducing a lattice self-energy with continuous momentum dependence. However, we show that the DCA$^+$ algorithm is plagued by a fundamental problem when its self-consistency equations are formulated using the bare Green's function of the cluster. This problem is most severe in the strongly correlated regime at low doping, where the DCA$^+$ self-energy becomes overly metallic and local, and persists to cluster sizes where the standard DCA has long converged. In view of the failure of the DCA$^+$ algorithm, we propose to complement DCA simulations with a post-interpolation procedure for single-particle and two-particle correlation functions to preserve continuous momentum dependence and the associated benefits in the DCA. We demonstrate the effectiveness of this practical approach with results for the half-filled and hole-doped two-dimensional Hubbard model.
△ Less
Submitted 25 April, 2020; v1 submitted 17 February, 2020;
originally announced February 2020.
-
Effects of momentum-dependent quasiparticle renormalization on the gap structure of iron-based superconductors
Authors:
Shinibali Bhattacharyya,
P. J. Hirschfeld,
Thomas A. Maier,
Douglas J. Scalapino
Abstract:
We discuss the influence of momentum-dependent correlations on the superconducting gap structure in iron-based superconductors. Within the weak coupling approach including self-energy effects at the one-loop spin-fluctuation level, we construct a dimensionless pairing strength functional which includes the effects of quasiparticle renormalization. The stationary solution of this equation determine…
▽ More
We discuss the influence of momentum-dependent correlations on the superconducting gap structure in iron-based superconductors. Within the weak coupling approach including self-energy effects at the one-loop spin-fluctuation level, we construct a dimensionless pairing strength functional which includes the effects of quasiparticle renormalization. The stationary solution of this equation determines the gap function at $T_c$. The resulting equations represent the simplest generalization of spin fluctuation pairing theory to include the effects of an anisotropic quasiparticle weight. We obtain good agreement with experimentally observed anisotropic gap structures in LiFeAs, indicating that the inclusion of quasiparticle renormalization effects in the existing weak-coupling theories can account for the observed anomalies in the gap structure of Fe-based superconductors.
△ Less
Submitted 17 May, 2020; v1 submitted 22 December, 2019;
originally announced December 2019.
-
Concurrent Expandable AMQs on the Basis of Quotient Filters
Authors:
Tobias Maier,
Peter Sanders,
Robert Williger
Abstract:
A quotient filter is a cache efficient AMQ data structure. Depending on the fill degree of the filter most insertions and queries only need to access one or two consecutive cache lines. This makes quotient filters fast compared to the more commonly used Bloom filters that incur multiple cache misses. However, concurrent Bloom filters are easy to implement and can be implemented lock-free while con…
▽ More
A quotient filter is a cache efficient AMQ data structure. Depending on the fill degree of the filter most insertions and queries only need to access one or two consecutive cache lines. This makes quotient filters fast compared to the more commonly used Bloom filters that incur multiple cache misses. However, concurrent Bloom filters are easy to implement and can be implemented lock-free while concurrent quotient filters are not as simple. Usually concurrent quotient filters work by using an external array of locks -- each protecting a region of the table. Accessing this array incurs one additional cache miss per operation. We propose a new locking scheme that has no memory overhead. Using this new locking scheme we achieve 1.8 times higher speedups than with the common external locking scheme.
Another advantage of quotient filters over Bloom filters is that a quotient filter can change its size when it is becoming full. We implement this growing technique for our concurrent quotient filters and adapt it in a way that allows unbounded growing while keeping a bounded false positive rate. We call the resulting data structure a fully expandable quotient filter. Its design is similar to scalable Bloom filters, but we exploit some concepts inherent to quotient filters to improve the space efficiency and the query speed.
We also propose quotient filter variants that are aimed to reduce the number of status bits (2-status-bit variant) or to simplify concurrent implementations (linear probing quotient filter). The linear probing quotient filter even leads to a lock-free concurrent filter implementation. This is especially interesting, since we show that any lock-free implementation of another common quotient filter variant would incur significant overheads in the form of additional data fields or multiple passes over the accessed data.
△ Less
Submitted 19 November, 2019;
originally announced November 2019.
-
Quantum-classical simulation of two-site dynamical mean-field theory on noisy quantum hardware
Authors:
Trevor Keen,
Thomas Maier,
Steven Johnston,
Pavel Lougovski
Abstract:
We report on a quantum-classical simulation of the single-band Hubbard model using two-site dynamical mean-field theory (DMFT). Our approach uses IBM's superconducting qubit chip to compute the zero-temperature impurity Green's function in the time domain and a classical computer to fit the measured Green's functions and extract their frequency domain parameters. We find that the quantum circuit s…
▽ More
We report on a quantum-classical simulation of the single-band Hubbard model using two-site dynamical mean-field theory (DMFT). Our approach uses IBM's superconducting qubit chip to compute the zero-temperature impurity Green's function in the time domain and a classical computer to fit the measured Green's functions and extract their frequency domain parameters. We find that the quantum circuit synthesis (Trotter) and hardware errors lead to incorrect frequency estimates, and subsequently to an inaccurate quasiparticle weight when calculated from the frequency derivative of the self-energy. These errors produce incorrect hybridization parameters that prevent the DMFT algorithm from converging to the correct self-consistent solution. To avoid this pitfall, we compute the quasiparticle weight by integrating the quasiparticle peaks in the spectral function. This method is much less sensitive to Trotter errors and allows the algorithm to converge to self-consistency for a half-filled Mott insulating system after applying quantum error mitigation techniques to the quantum simulation data.
△ Less
Submitted 9 November, 2021; v1 submitted 21 October, 2019;
originally announced October 2019.
-
Pairing in the Two-Dimensional Hubbard Model from Weak to Strong Coupling
Authors:
Astrid T. Rømer,
Thomas A. Maier,
Andreas Kreisel,
Ilya Eremin,
P. J. Hirschfeld,
Brian M. Andersen
Abstract:
The Hubbard model is the simplest model that is believed to exhibit superconductivity arising from purely repulsive interactions, and has been extensively applied to explore a variety of unconventional superconducting systems. Here we study the evolution of the leading superconducting instabilities of the single-orbital Hubbard model on a two-dimensional square lattice as a function of onsite Coul…
▽ More
The Hubbard model is the simplest model that is believed to exhibit superconductivity arising from purely repulsive interactions, and has been extensively applied to explore a variety of unconventional superconducting systems. Here we study the evolution of the leading superconducting instabilities of the single-orbital Hubbard model on a two-dimensional square lattice as a function of onsite Coulomb repulsion $U$ and band filling by calculating the irreducible particle-particle scattering vertex obtained from dynamical cluster approximation (DCA) calculations, and compare the results to both perturbative Kohn-Luttinger (KL) theory as well as the widely used random phase approximation (RPA) spin-fluctuation pairing scheme. Near half-filling we find remarkable agreement of the hierarchy of the leading pairing states between these three methods, implying adiabatic continuity between weak- and strong-coupling pairing solutions of the Hubbard model. The $d_{x^2-y^2}$-wave instability is robust to increasing $U$ near half-filling as expected. Away from half filling, the predictions of KL and RPA at small $U$ for transitions to other pair states agree with DCA at intermediate $U$ as well as recent diagrammatic Monte Carlo calculations. RPA results fail only in the very dilute limit, where it yields a $d_{xy}$ ground state instead of a $p$-wave state established by diagrammatic Monte Carlo and low-order perturbative methods, as well as our DCA calculations. We discuss the origins of this discrepancy, highlighting the crucial role of the vertex corrections neglected in the RPA approach. Overall, comparison of the various methods over the entire phase diagram strongly suggests a smooth crossover of the superconducting interaction generated by local Hubbard interactions between weak and strong coupling.
△ Less
Submitted 2 September, 2019;
originally announced September 2019.
-
Two-particle Correlation Functions in Cluster Perturbation Theory: Hubbard Spin Susceptibilities
Authors:
P. T. Raum,
G. Alvarez,
Thomas Maier,
V. W. Scarola
Abstract:
Cluster Perturbation Theory (CPT) is a computationally economic method commonly used to estimate the momentum and energy resolved single-particle Green's function. It has been used extensively in direct comparisons with experiments that effectively measure the single-particle Green's function, e.g., angle-resolved photoemission spectroscopy. However, many experimental observables are given by two-…
▽ More
Cluster Perturbation Theory (CPT) is a computationally economic method commonly used to estimate the momentum and energy resolved single-particle Green's function. It has been used extensively in direct comparisons with experiments that effectively measure the single-particle Green's function, e.g., angle-resolved photoemission spectroscopy. However, many experimental observables are given by two-particle correlation functions. CPT can be extended to compute two-particle correlation functions by approximately solving the Bethe-Salpeter equation. We implement this method and focus on the transverse spin-susceptibility, measurable via inelastic neutron scattering or with optical probes of atomic gases in optical lattices. We benchmark the method with the one-dimensional Fermi-Hubbard model at half filling by comparing with known results.
△ Less
Submitted 18 February, 2020; v1 submitted 27 August, 2019;
originally announced August 2019.
-
Analytic Continuation of Noisy Data Using Adams Bashforth ResNet
Authors:
Xuping Xie,
Feng Bao,
Thomas Maier,
Clayton Webster
Abstract:
We propose a data-driven learning framework for the analytic continuation problem in numerical quantum many-body physics. Designing an accurate and efficient framework for the analytic continuation of imaginary time using computational data is a grand challenge that has hindered meaningful links with experimental data. The standard Maximum Entropy (MaxEnt)-based method is limited by the quality of…
▽ More
We propose a data-driven learning framework for the analytic continuation problem in numerical quantum many-body physics. Designing an accurate and efficient framework for the analytic continuation of imaginary time using computational data is a grand challenge that has hindered meaningful links with experimental data. The standard Maximum Entropy (MaxEnt)-based method is limited by the quality of the computational data and the availability of prior information. Also, the MaxEnt is not able to solve the inversion problem under high level of noise in the data. Here we introduce a novel learning model for the analytic continuation problem using a Adams-Bashforth residual neural network (AB-ResNet). The advantage of this deep learning network is that it is model independent and, therefore, does not require prior information concerning the quantity of interest given by the spectral function. More importantly, the ResNet-based model achieves higher accuracy than MaxEnt for data with higher level of noise. Finally, numerical examples show that the developed AB-ResNet is able to recover the spectral function with accuracy comparable to MaxEnt where the noise level is relatively small.
△ Less
Submitted 24 May, 2019;
originally announced May 2019.
-
Effective pairing interaction in a system with an incipient band
Authors:
Thomas A. Maier,
Vivek Mishra,
Douglas J. Scalapino
Abstract:
The nature and mechanism of superconductivity in the extremely electron-doped FeSe based superconductors continues to be a matter of debate. In these systems, the hole-like band has moved below the Fermi energy, and various spin-fluctuation theories involving pairing between states near the electron Fermi surface and states of this incipient band have been proposed. Here, using a dynamic cluster q…
▽ More
The nature and mechanism of superconductivity in the extremely electron-doped FeSe based superconductors continues to be a matter of debate. In these systems, the hole-like band has moved below the Fermi energy, and various spin-fluctuation theories involving pairing between states near the electron Fermi surface and states of this incipient band have been proposed. Here, using a dynamic cluster quantum Monte Carlo calculation for a bilayer Hubbard model we show that the pairing in these systems can be understood in terms of an effective retarded attractive interaction between electrons near the electron Fermi surface.
△ Less
Submitted 29 January, 2019;
originally announced January 2019.
-
DCA++: A software framework to solve correlated electron problems with modern quantum cluster methods
Authors:
Urs R. Hähner,
Gonzalo Alvarez,
Thomas A. Maier,
Raffaele Solcà,
Peter Staar,
Michael S. Summers,
Thomas C. Schulthess
Abstract:
We present the first open release of the DCA++ project, a high-performance research software framework to solve quantum many-body problems with cutting edge quantum cluster algorithms. DCA++ implements the dynamical cluster approximation (DCA) and its DCA$^+$ extension with a continuous self-energy. The algorithms capture nonlocal correlations in strongly correlated electron systems, thereby givin…
▽ More
We present the first open release of the DCA++ project, a high-performance research software framework to solve quantum many-body problems with cutting edge quantum cluster algorithms. DCA++ implements the dynamical cluster approximation (DCA) and its DCA$^+$ extension with a continuous self-energy. The algorithms capture nonlocal correlations in strongly correlated electron systems, thereby giving insight into high-$T_c$ superconductivity. The code's scalability allows efficient usage of systems at all scales, from workstations to leadership computers. With regard to the increasing heterogeneity of modern computing machines, DCA++ provides portable performance on conventional and emerging new architectures, such as hybrid CPU-GPU, sustaining multiple petaflops on ORNL's Titan and CSCS' Piz Daint supercomputers. Moreover, we show how sustainable and scalable development of the code base has been achieved by adopting standard techniques of the software industry. These include employing a distributed version control system, applying test-driven development and following continuous integration.
△ Less
Submitted 2 September, 2019; v1 submitted 4 January, 2019;
originally announced January 2019.
-
Pairfield Fluctuations of the 2D Hubbard Model
Authors:
Thomas A. Maier,
Douglas J. Scalapino
Abstract:
At temperatures above the superconducting transition temperature the pairfield susceptibility provides information on the nature of the pairfield fluctuations. Here we study the d-wave pairfield susceptibility of a 2D Hubbard model for dopings which have a pseudogap (PG) and for dopings which do not. One knows that in both cases there will be a region of Kosterlitz-Thouless fluctuations as the tra…
▽ More
At temperatures above the superconducting transition temperature the pairfield susceptibility provides information on the nature of the pairfield fluctuations. Here we study the d-wave pairfield susceptibility of a 2D Hubbard model for dopings which have a pseudogap (PG) and for dopings which do not. One knows that in both cases there will be a region of Kosterlitz-Thouless fluctuations as the transition at $T_{KT}$ is approached. Above this region we find evidence for Emery-Kivelson phase fluctuations for dopings with a PG and Gaussian amplitude fluctuations for dopings without a PG.
△ Less
Submitted 23 October, 2018;
originally announced October 2018.
-
$d$-wave and $s^\pm$ Pairing Strengths in Ba$_2$CuO$_{3+δ}$
Authors:
Thomas A. Maier,
Tom Berlijn,
Doug J. Scalapino
Abstract:
Using a first principles derived two-orbital model ($d_{x^2-y^2}$, $d_{3z^2-r^2}$) and a random phase approximation treatment of the spin-fluctuation pairing vertex, we calculate the $d$-wave and $s^\pm$-wave pairing strengths for a model of the recently discovered Ba$_2$CuO$_{3+δ}$ superconductor. We find that there is significant pairing strength in both of these channels. These results provide…
▽ More
Using a first principles derived two-orbital model ($d_{x^2-y^2}$, $d_{3z^2-r^2}$) and a random phase approximation treatment of the spin-fluctuation pairing vertex, we calculate the $d$-wave and $s^\pm$-wave pairing strengths for a model of the recently discovered Ba$_2$CuO$_{3+δ}$ superconductor. We find that there is significant pairing strength in both of these channels. These results provide an interesting perspective on the relationship between the cuprates and Fe-based superconductors and the high-$T_c$ pairing mechanism.
△ Less
Submitted 11 September, 2018;
originally announced September 2018.
-
Neutron spin resonance as a probe of Fermi surface nesting and superconducting gap symmetry in Ba$_{0.67}$K$_{0.33}$(Fe$_{1-x}$Co$_{x}$)$_{2}$As$_{2}$
Authors:
Rui Zhang,
Weiyi Wang,
Thomas Maier,
Meng Wang,
Matthew B. Stone,
Songxue Chi,
Barry Winn,
Pengcheng Dai
Abstract:
We use inelastic neutron scattering to study energy and wave vector dependence of the superconductivity-induced resonance in hole-doped Ba$_{0.67}$K$_{0.33}$(Fe$_{1-x}$Co$_{x}$)$_{2}$As$_{2}$ ($x=0,0.08$ with $T_c\approx 37, 28$ K, respectively). In previous work on electron-doped Ba(Fe$_{0.963}$Ni$_{0.037}$)$_2$As$_2$ ($T_N=26$ K and $T_c=17$ K), the resonance is found to peak sharply at the anti…
▽ More
We use inelastic neutron scattering to study energy and wave vector dependence of the superconductivity-induced resonance in hole-doped Ba$_{0.67}$K$_{0.33}$(Fe$_{1-x}$Co$_{x}$)$_{2}$As$_{2}$ ($x=0,0.08$ with $T_c\approx 37, 28$ K, respectively). In previous work on electron-doped Ba(Fe$_{0.963}$Ni$_{0.037}$)$_2$As$_2$ ($T_N=26$ K and $T_c=17$ K), the resonance is found to peak sharply at the antiferromagnetic (AF) ordering wave vector ${\bf Q}_{\rm AF}$ along the longitudinal direction, but disperses upwards away from ${\bf Q}_{\rm AF}$ along the transverse direction. For hole doped $x=0, 0.08$ without AF order, we find that the resonance displays ring-like upward dispersion away from ${\bf Q}_{\rm AF}$ along both the longitudinal and transverse directions. By comparing these results with calculations using the random phase approximation, we conclude that the dispersive resonance is a direct signature of isotropic superconducting gaps arising from nested hole-electron Fermi surfaces.
△ Less
Submitted 18 August, 2018;
originally announced August 2018.
-
Stabilization of $s$-wave superconductivity through arsenic $p$-orbital hybridization in electron-doped BaFe$_2$As$_2$
Authors:
David W. Tam,
Tom Berlijn,
Thomas A. Maier
Abstract:
Using random-phase approximation spin-fluctuation theory, we study the influence of the hybridization between iron $d$-orbitals and pnictide $p$-orbitals on the superconducting pairing state in iron-based superconductors. The calculations are performed for a 16-orbital Hubbard-Hund tight-binding model of BaFe$_2$As$_2$ that includes the As-$p$ orbital degrees of freedom in addition to the Fe-$d$ o…
▽ More
Using random-phase approximation spin-fluctuation theory, we study the influence of the hybridization between iron $d$-orbitals and pnictide $p$-orbitals on the superconducting pairing state in iron-based superconductors. The calculations are performed for a 16-orbital Hubbard-Hund tight-binding model of BaFe$_2$As$_2$ that includes the As-$p$ orbital degrees of freedom in addition to the Fe-$d$ orbitals and compared to calculations for a 10-orbital Fe-$d$ only model. In both models we find a leading $s^\pm$ pairing state and a subleading $d_ {x^2-y^2}$-wave state in the parent compound. Upon doping, we find that the $s^\pm$ state remains the leading state in the 16-orbital model up to a doping level of 0.475 electrons per unit cell, at which the hole Fermi surface pockets at the zone center start to disappear. This is in contrast to the 10-orbital model, where the $d$-wave state becomes the leading state at a doping of less than 0.2 electrons. This improved stability of $s^\pm$ pairing is found to arise from a decrease of $d_{xy}$ orbital weight on the electron pockets due to hybridization with the As-$p$ orbitals and the resulting reduction of near $(π,π)$ spin-fluctuation scattering which favors the competing $d$-wave state. These results show that the orbital dependent hybridization of Fermi surface Bloch states with the usually neglected $p$-orbital states is an important ingredient in an improved itinerant pairing theory.
△ Less
Submitted 11 July, 2018; v1 submitted 10 May, 2018;
originally announced May 2018.
-
Phase competition in a one-dimensional three-orbital Hubbard-Holstein model
Authors:
Shaozhi Li,
Yanfei Tang,
Thomas A. Maier,
Steven Johnston
Abstract:
We study the interplay between the electron-phonon (e-ph) and on-site electron-electron (e-e) interactions in a three-orbital Hubbard-Holstein model on an extended one-dimensional lattice using determinant quantum Monte Carlo. For weak e-e and e-ph interactions, we observe a competition between an orbital-selective Mott phase (OSMP) and a (multicomponent) charge-density-wave (CDW) insulating phase…
▽ More
We study the interplay between the electron-phonon (e-ph) and on-site electron-electron (e-e) interactions in a three-orbital Hubbard-Holstein model on an extended one-dimensional lattice using determinant quantum Monte Carlo. For weak e-e and e-ph interactions, we observe a competition between an orbital-selective Mott phase (OSMP) and a (multicomponent) charge-density-wave (CDW) insulating phase, with an intermediate metallic phase located between them. For large e-e and e-ph couplings, the OSMP and CDW phases persist, while the metallic phase develops short-range orbital correlations and becomes insulating when both the e-e and e-ph interactions are large but comparable. Many of our conclusions are in line with those drawn from a prior dynamical mean field theory study of the two-orbital Hubbard-Holstein model [Phys. Rev. B 95, 12112(R) (2017)] in infinite dimension, suggesting that the competition between the e-ph and e-e interactions in multiorbital Hubbard-Holstein models leads to rich physics, regardless of the dimension of the system.
△ Less
Submitted 22 April, 2018;
originally announced April 2018.
-
Doping evolution of charge and spin excitations in two-leg Hubbard ladders: comparing DMRG and RPA+FLEX results
Authors:
A. Nocera,
Y. Wang,
N. D. Patel,
G. Alvarez,
T. A. Maier,
E. Dagotto,
S. Johnston
Abstract:
We study the magnetic and charge dynamical response of a Hubbard model in a two-leg ladder geometry using the density matrix renormalization group (DMRG) method and the random phase approximation within the fluctuation-exchange approximation (RPA+FLEX). Our calculations reveal that RPA+FLEX can capture the main features of the magnetic response from weak up to intermediate Hubbard repulsion for do…
▽ More
We study the magnetic and charge dynamical response of a Hubbard model in a two-leg ladder geometry using the density matrix renormalization group (DMRG) method and the random phase approximation within the fluctuation-exchange approximation (RPA+FLEX). Our calculations reveal that RPA+FLEX can capture the main features of the magnetic response from weak up to intermediate Hubbard repulsion for doped ladders, when compared with the numerically exact DMRG results. However, while at weak Hubbard repulsion both the spin and charge spectra can be understood in terms of weakly-interacting electron-hole excitations across the Fermi surface, at intermediate coupling DMRG shows gapped spin excitations at large momentum transfer that remain gapless within the RPA+FLEX approximation. For the charge response, RPA+FLEX can only reproduce the main features of the DMRG spectra at weak coupling and high doping levels, while it shows an incoherent character away from this limit. Overall, our analysis shows that RPA+FLEX works surprisingly well for spin excitations at weak and intermediate Hubbard $U$ values even in the difficult low-dimensional geometry such as a two-leg ladder. Finally, we discuss the implications of our results for neutron scattering and resonant inelastic x-ray scattering experiments on two-leg ladder cuprate compounds.
△ Less
Submitted 22 March, 2018;
originally announced March 2018.
-
Realization of a hole-doped Mott insulator on a triangular silicon lattice
Authors:
Fangfei Ming,
Steve Johnston,
Daniel Mulugeta,
Tyler S. Smith,
Paolo Vilmercati,
Geunseop Lee,
Thomas A. Maier,
Paul C. Snijders,
Hanno H. Weitering
Abstract:
The physics of doped Mott insulators is at the heart of some of the most exotic physical phenomena in materials research including insulator-metal transitions, colossal magneto-resistance, and high-temperature superconductivity in layered perovskite compounds. Advances in this field would greatly benefit from the availability of new material systems with similar richness of physical phenomena, but…
▽ More
The physics of doped Mott insulators is at the heart of some of the most exotic physical phenomena in materials research including insulator-metal transitions, colossal magneto-resistance, and high-temperature superconductivity in layered perovskite compounds. Advances in this field would greatly benefit from the availability of new material systems with similar richness of physical phenomena, but with fewer chemical and structural complications in comparison to oxides. Using scanning tunneling microscopy and spectroscopy, we show that such a system can be realized on a silicon platform. Adsorption of one-third monolayer of Sn atoms on a Si(111) surface produces a triangular surface lattice with half-filled dangling bond orbitals. Modulation hole-doping of these dangling bonds unveils clear hallmarks of Mott physics, such as spectral weight transfer and the formation of quasi-particle states at the Fermi level, well-defined Fermi contour segments, and a sharp singularity in the density of states. These observations are remarkably similar to those made in complex oxide materials, including high-temperature superconductors, but highly extraordinary within the realm of conventional sp-bonded semiconductor materials. It suggests that exotic quantum matter phases can be realized and engineered on silicon-based materials platforms.
△ Less
Submitted 22 October, 2017;
originally announced October 2017.
-
d-wave superconductivity in the presence of a near neighbor Coulomb repulsion
Authors:
M. Jiang,
U. R. Hahner,
T. C. Schulthess,
T. A. Maier
Abstract:
Dynamic cluster quantum Monte Carlo calculations for a doped two-dimensional extended Hubbard model are used to study the stability and dynamics of $d$-wave pairing when a near neighbor Coulomb repulsion $V$ is present in addition to the on-site Coulomb repulsion $U$. We find that $d$-wave pairing and the superconducting transition temperature $T_c$ are only weakly suppressed as long as $V$ does n…
▽ More
Dynamic cluster quantum Monte Carlo calculations for a doped two-dimensional extended Hubbard model are used to study the stability and dynamics of $d$-wave pairing when a near neighbor Coulomb repulsion $V$ is present in addition to the on-site Coulomb repulsion $U$. We find that $d$-wave pairing and the superconducting transition temperature $T_c$ are only weakly suppressed as long as $V$ does not exceed $U/2$. This stability is traced to the strongly retarded nature of pairing that allows the $d$-wave pairs to minimize the repulsive effect of $V$. When $V$ approaches $U/2$, large momentum charge fluctuations are found to become important and to give rise to a more rapid suppression of $d$-wave pairing and $T_c$ than for smaller $V$.
△ Less
Submitted 19 July, 2017;
originally announced July 2017.
-
Decrease of d-wave pairing strength in spite of the persistence of magnetic excitations in the overdoped Hubbard model
Authors:
Edwin W. Huang,
Douglas J. Scalapino,
Thomas A. Maier,
Brian Moritz,
Thomas P. Devereaux
Abstract:
Evidence for the presence of high energy magnetic excitations in overdoped La$_{2-x}$Sr$_x$CuO$_4$ (LSCO) has raised questions regarding the role of spin-fluctuations in the pairing mechanism. If they remain present in overdoped LSCO, why does $T_c$ decrease in this doping regime? Here, using results for the dynamic spin susceptibility ${\rm Im}χ(q,ω)$ obtained from a determinantal quantum Monte C…
▽ More
Evidence for the presence of high energy magnetic excitations in overdoped La$_{2-x}$Sr$_x$CuO$_4$ (LSCO) has raised questions regarding the role of spin-fluctuations in the pairing mechanism. If they remain present in overdoped LSCO, why does $T_c$ decrease in this doping regime? Here, using results for the dynamic spin susceptibility ${\rm Im}χ(q,ω)$ obtained from a determinantal quantum Monte Carlo (DQMC) calculation for the Hubbard model we address this question. We find that while high energy magnetic excitations persist in the overdoped regime, they lack the momentum to scatter pairs between the anti-nodal regions. It is the decrease in the spectral weight at large momentum transfer, not observed by resonant inelastic X-ray scattering (RIXS), which leads to a reduction in the $d$-wave spin-fluctuation pairing strength.
△ Less
Submitted 10 May, 2017;
originally announced May 2017.
-
Dynamic Space Efficient Hashing
Authors:
Tobias Maier,
Peter Sanders
Abstract:
We consider space efficient hash tables that can grow and shrink dynamically and are always highly space efficient, i.e., their space consumption is always close to the lower bound even while growing and when taking into account storage that is only needed temporarily. None of the traditionally used hash tables have this property. We show how known approaches like linear probing and bucket cuckoo…
▽ More
We consider space efficient hash tables that can grow and shrink dynamically and are always highly space efficient, i.e., their space consumption is always close to the lower bound even while growing and when taking into account storage that is only needed temporarily. None of the traditionally used hash tables have this property. We show how known approaches like linear probing and bucket cuckoo hashing can be adapted to this scenario by subdividing them into many subtables or using virtual memory overcommitting. However, these rather straightforward solutions suffer from slow amortized insertion times due to frequent reallocation in small increments.
Our main result is DySECT ({\bf Dy}namic {\bf S}pace {\bf E}fficient {\bf C}uckoo {\bf T}able) which avoids these problems. DySECT consists of many subtables which grow by doubling their size. The resulting inhomogeneity in subtable sizes is equalized by the flexibility available in bucket cuckoo hashing where each element can go to several buckets each of which containing several cells. Experiments indicate that DySECT works well with load factors up to 98\%. With up to 2.7 times better performance than the next best solution.
△ Less
Submitted 2 May, 2017;
originally announced May 2017.