arXiv Curator

The discovery of superconductivity with Tc ~ 80 K in the nickelate Ruddlesden-Popper bilayer La3Ni2O7 at high pressure has opened a new platform for unconventional superconductivity, followed by the subsequent observation of superconductivity in trilayer La4Ni3O10, also at high pressure. Remarkably, ambient-pressure superconductivity was also observed recently in La3Ni2O7 ultra-thin films when grown on substrates that provide compressive strain. This discovery significantly extends the type of experimental techniques that can be used in nickelates, previously limited due to the high-pressure constraint. Discussing the similarities and differences among these nickel oxides will provide new insights into understanding the mechanism of high-Tc superconductivity in correlated electron systems. In this paper, we review the experimental and theoretical progress on RuddlesdenPopper nickelates, with emphasis on thin films, and discuss future perspectives and research directions.

Ruddlesden-Popper nickelates offer a new route to high-temperature superconductivity beyond the cuprates and iron-pnictides. However, the electronic reorganization that enables superconductivity in bilayer nickelates remain unresolved, largely due to the difficulty of directly probing the superconducting phase. Here, we overcome this limitation by stabilizing superconducting (La,Pr)$_3$Ni$_2$O$_7$ thin films with a protective capping layer, thereby enabling direct spectroscopic access via X-ray absorption and resonant inelastic X-ray scattering. We resolve the evolution of in-plane and out-of-plane electronic states, spin and orbital excitations, and spin-density-waves across insulating, superconducting, and metallic regimes. Combining experimental results with theoretical analysis, we show that the in-plane $d_{x^2-y^2}$ states form an itinerant backbone, while superconductivity emerges only when coherent $d_{z^2}$-$p_z$-$d_{z^2}$ interlayer hybridization develops, accompanied by suppressed static spin order and strongly damped spin excitations. Oxygen stoichiometry and epitaxial strain both act on this interlayer channel, placing superconductivity within a narrow window of interlayer coherence and correlation strength. These findings identify the microscopic ingredients required for superconductivity in bilayer nickelates and provide a multiorbital picture of its emergence.

The discovery of superconductivity in bulk Ruddlesden-Popper La$_3$Ni$_2$O$_7$(LNO327) under high hydrostatic pressure has redefined the recent experimental consensus that nickelate superconductivity is restricted to systems with a $3d^9$ electronic configuration and square-planar coordination. However, the structural and electronic prerequisites for stabilizing superconductivity, whether under pressure or at ambient conditions in the case of thin films, remain poorly understood, largely due to the metastable nature of the LNO327 phase. Here, we present a detailed structural study of epitaxial La$_3$Ni$_2$O$_{7-δ}$ thin films by using scanning transmission electron microscopy (STEM) combined with electron energy loss spectroscopy (EELS). Grown via pulsed laser deposition onto SrLaAlO$_4$ substrates, those films exhibit distinct superconducting properties as a function of the different post-annealing conditions used. By correlating the rich landscape of stacking polymorphs with transport behavior, this work establishes a framework for understanding the metastable superconducting phase in bilayer nickelate thin films. Our findings underscore the critical role of homogeneity in oxygen stoichiometry, epitaxial strain and structural motif in stabilizing superconductivity, offering a clear pathway for designing ambient-pressure superconducting nickelates.

The discovery of superconductivity in the Ruddlesden-Popper (RP) nickelates has opened a new chapter in the search for high superconducting transition temperatures ($T_\mathrm{c}$) materials. A central and puzzling feature of this family is the wide variation in $T_\mathrm{c}$ despite their common NiO$_2$ building blocks, as highlighted by the recent observation of superconductivity at $\sim$ 30 K in trilayer $\mathrm{La_4Ni_3O_{10}}$, significantly lower than 80 K reported in bilayer $\mathrm{La_3Ni_2O_7}$. Understanding the factors that control $T_\mathrm{c}$ in this family is therefore of paramount importance. Here, we use resonant inelastic x-ray scattering (RIXS) to investigate the electronic and magnetic excitations of $\mathrm{La_4Ni_3O_{10}}$ in direct comparison with its bilayer counterpart. Our results reveal a markedly different landscape. $\mathrm{La_4Ni_3O_{10}}$ exhibits a more itinerant character, evidenced by broader Ni $dd$ orbital excitations and a strong Ni 3$d$ fluorescence continuum, suggesting weaker electronic correlations than in the bilayer. Despite this, well-defined collective spin excitations persist, including dispersive acoustic and optical magnon branches alongside an incommensurate spin density wave. Using linear spin wave theory, we extract the interlayer superexchange interaction ($J_z$) to be $\sim$ 22 meV, much smaller than that in $\mathrm{La_3Ni_2O_7}$. The weaker correlation and reduced interlayer exchange together provide a consistent explanation for the substantially lower $T_\mathrm{c}$ in the trilayer compound. Our findings establish interlayer magnetic coupling and electronic correlation as key parameters governing superconductivity in layered nickelates and offer critical constraints for understanding the pairing mechanism in this emerging family.

Ruddlesden-Popper (RP) nickelates have recently emerged as a new family of high-temperature superconductors. In bilayer RP nickelates, magnetic excitations with large exchange couplings have been observed, supporting a spin-mediated pairing mechanism. Whether comparable spin correlations persist in trilayer nickelates, however, remains unknown. Here, we present a Ni $L$-edge resonant inelastic X-ray scattering (RIXS) study of La$_4$Ni$_3$O$_{10}$ single crystals. While the orbital excitations remain similar to those of La$_3$Ni$_2$O$_{7}$, the collective spin excitations in La$_4$Ni$_3$O$_{10}$ exhibit a comparable bandwidth of about $60$ meV but substantially suppressed spectral weight, implying a weaker electronic correlation in the trilayer compounds. Our results underscore the three-dimensional and multi-orbital electronic character in La$_4$Ni$_3$O$_{10}$, highlighting important differences from the bilayer nickelates. These findings provide crucial insights into the evolution of magnetism across the RP nickelate family and its connection to superconductivity.

Using DFT+DMFT we study the normal-state electronic structure and magnetic correlations of the recently discovered alternating single-layer bilayer Raddlesden-Popper nickelate La$_5$Ni$_3$O$_{11}$ (1212-LNO). Our results exhibit qualitative differences for the structurally distinct single-layer and bilayer Ni ions, implying the importance of confinement and orbital-dependent correlations. The Ni $e_g$ electronic states originating from the bilayer Ni ions form strongly renormalized quasiparticle bands with a large enhancement factor $m^*/m \sim 3.5$ and 4.2 for the Ni $x^2-y^2$ and $3z^2-r^2$ orbitals, respectively. Moreover, the $e_g$ states of the single-layer Ni ions exhibit an orbital-selective Mott insulating state, with a narrow energy gap for the Ni $3z^2-r^2$ states and metallic, strongly incoherent (non-Fermi-liquid) Ni $x^2-y^2$ ones. Our analysis of magnetic correlations suggests the formation of intertwined spin and charge density wave stripes in the bilayer NiO$_6$ slab, in close similarity to the double-layer material. We refine two major instabilities, a leading one is associated with a wave vector ${\bf Q}=(\frac{1}{3},\frac{1}{3})$ (``up-down-0" spin pattern), competing with the bicollinear $(\frac{1}{4},\frac{1}{4})$ (``up-up-down-down") stripe state. The single-layer Ni $3d$ electrons exhibit instability towards the Néel-type magnetic state. Under pressure, 1212-LNO undergoes an orbital-selective Mott insulator-to-metal phase transition, associated with metallization of the single-layer Ni $e_g$ states. As a result, the single-layer Ni $e_g$ bands exhibit non-Fermi-liquid (bad metal) behavior with strongly incoherent spectral weights near $E_F$. We note that correlation effects result in a reconstruction of magnetic correlations as compared to that obtained within DFT. In fact, we observe a crossover from single-layer to double-layer dominated magnetic correlations.

Magnetic correlations of doped Mott insulators hold the key to the unusual characteristics of many quantum materials. Recent experiments with ultracold atoms in optical lattices have provided new information about the magnetic properties of the Fermi-Hubbard model on a square lattice. We demonstrate that recent measurements indicate that a single doping-dependent energy scale determines both static correlations and dynamical response of these systems. To understand these experimental findings, we employ a self-consistent formalism to describe the coupling between antiferromagnetic magnons and doped holes, and we uncover the emergence of a universal magnetic energy scale at finite doping, which we denote by $J^*$. We present the single- and two-magnon spectral properties at finite doping and discuss the appearance of a bimagnon peak in lattice-modulation spectroscopy, at frequencies set by $J^*$. Furthermore, we argue that this same energy scale sets the onset of pseudogap phenomena, leading to the hypothesis $k_BT^* = c J^*$, with $c$ an order one number. We identify another low-energy scale emerging from our analysis of magnetic excitations, and argue that it controls the stability of Néel order at the lowest temperatures, ultimately driving a transition to an incommensurate spin-density-wave at finite doping. We discuss the relation between this low-energy scale and the nature of fermionic quasiparticles. Our analysis suggests that stability of the commensurate antiferromagentic phase at finite doping can be controlled experimentally by introducing additional quasiparticle broadening via disorder or low-frequency noise.

High temperature superconductivity in cuprate superconductors is generally considered to be generated from doping the Mott insulators. The fundamental nature of the doped parent compounds as well as the microscopic origin of electron pairing remain critical issues in understanding the emergence of superconductivity. Here, using high-resolution spatially-resolved laser angle-resolved photoemission spectroscopy, we investigate the intrinsic electronic structures of the CuO$_2$ planes in multilayer cuprates Bi$_2$Sr$_2$Ca$_{n-1}$Cu$_n$O$_{2n+4+δ}$ (n=5$\sim$8). The inner CuO$_2$ planes are well shielded from the disorders and provide a rare and ideal platform to probe the intrinsic electronic phase diagram. We observe well-defined Fermi pockets with hole doping levels as low as 0.007, demonstrating an abrupt transition from the parent Mott insulator to a metallic state upon the introduction of an infinitesimal amount of doping. The innermost CuO$_2$ planes (IP$_0$) display gapless Fermi pockets, while the second innermost planes (IP$_1$) exhibit anisotropic superconducting gaps up to $\sim$33$\,$meV, indicative of robust electron pairing coexisting with strong antiferromagnetic order. Our findings provide a revised framework for understanding the doping-driven transitions and pairing mechanisms in cuprate superconductors.

Finite-Hamiltonian impurity solvers provide direct real-frequency spectra and a natural route to enlarged impurity Hamiltonians, but their applicability is limited by the rapid Hilbert-space growth with the number of bath or other added one-particle orbitals. We introduce an active-learning extension of adaptive-truncation configuration interaction (AL-ATCI) that identifies the determinant manifold relevant to the correlated state. The approximation is systematically controlled by the query size N_query, which also provides an internal convergence parameter when no external benchmark is available. Over the benchmark range studied here, the computational cost grows only weakly with bath size, because enlarging the bath mainly expands the combinatorial determinant space rather than the physically relevant manifold. In dynamical mean-field-theory benchmarks for the one-dimensional Hubbard model, AL-ATCI reproduces exact-diagonalization accuracy and extends cellular calculations to clusters as large as N_c = 10. For a three-orbital rotationally invariant Sr2RuO4 impurity problem, we demonstrate systematic convergence of dynamical quantities and a highly compressed configuration space as N_b is increased from 9 to 18. These results substantially alleviate the bath-discretization bottleneck of exact-diagonalization- and configuration-interaction-based impurity solvers and make large-bath and enlarged-orbital calculations more practical.

Beyond the quasi-two-dimensional (2D) paradigm of cuprates, the role of the third dimension of the Ruddlesden-Popper bilayer nickelates is essential to decoding their superconducting mechanism. Here, using angle-resolved photoemission spectroscopy (ARPES) with varied photon energies, we systematically investigate the electronic band structures in three dimensions for superconducting (La,Pr,Sm)$_3$Ni$2$O$7$/SrLaAlO$4$ thin films (superconducting onset temperature $T_c^{\text{onset}} \sim 48$ K) transferred via a cryogenic ultra-high vacuum suitcase. We reveal an orbital-dependent dimensionality: while the $d{x^2-y^2}$-dominant bands exhibit a quasi-2D character, the $d{z^2}$-dominant band displays a finite $k_z$ dispersion. Finite energy gaps are identified on all observed bands across multiple high-symmetry directions. Systematic temperature-dependent analysis characterizes the superconducting nature of the gap on the $d{z^2}$-derived band, revealing a large gap $Δ\sim 18$ meV and a ratio $2Δ/k_BT_c\sim 8$ exceeding the weak-coupling BCS limit. The suppression of spectral weight near the Fermi level persists above the superconducting transition temperature. Ubiquitous waterfall-like spectral features evidence the presence of electron interactions. These results underscore the role of the $d_{z^2}$ orbital and correlations, placing constraints on theoretical models for nickelate superconductivity.

We present a systematic study of the electronic structure of strained La$_3$Ni$_2$O$_7$ thin films. We show that biaxial compressive strain mainly elongates the outer apical Ni-O bond while leaving the inner apical Ni-O bond nearly unchanged. As a result, the Jahn-Teller splitting $Δ_{JT}$ is strongly enhanced, whereas the interlayer $d_{z^2}$ hopping $t_\perp^z$ changes only weakly. Since superconductivity is widely believed to emerge only below a critical in-plane lattice constant, our results identify the strain-enhanced $Δ_{JT}$ as the relevant microscopic tuning parameter. Consistently, the calculated Fermi surfaces and Hall response for LaAlO$_3$ and SrLaAlO$_4$ substrates agree with ARPES and Hall measurements. Our results identify Jahn-Teller distortion as a key tuning parameter in strained La$_3$Ni$_2$O$_7$ and support its central role in optimizing superconductivity in bilayer nickelates.

Discovery of high-$T_c$ superconductivity (SC) in the bilayer nickelate series La$_3$Ni$_2$O$_7$ have attracted substantial interest, providing a new platform for exploring unconventional SC. Certain experimental evidence has pointed to a correlated electronic nature, which is the driving force responsible for its high critical temperature ($T_c$). This work reviews the SC in La$_3$Ni$_2$O$_7$, with a particular focus on theoretical understanding of its pairing mechanism driven by this strong-coupling, Hund-assisted scenario. The electronic landscape is governed by two $E_g$-orbitals within the bilayer structure of NiO$_2$ planes. The $3d_{z^2}$ orbital is nearly half-filled and exhibits a stronger localized character, while the $3d_{x^2-y^2}$ is approximately quarter-filled and remains highly itinerant. The localized $3d_{z^2}$ orbitals experience robust interlayer hybridization, mediated by the $2p_z$ orbitals of the inner apical oxygen atoms. This hybridization generates a strong interlayer antiferromagnetic (AFM) exchange. In the strong coupling regime, Hund's rule coupling aligns the spins of the two $E_g$ orbitals on the same nickel site. The strong interlayer AFM exchange is effectively transferred to the itinerant $3d_{x^2-y^2}$ orbital, generating an effective coupling $J_{\perp}$ within this orbital. This mechanism is captured by a minimal strong-coupling bilayer $t$-$J$-$J_{\perp}$ model for the $3d_{x^2-y^2}$ band. Driven by $J_{\perp}$, $3d_{x^2-y^2}$ electrons can form interlayer Cooper pairs, leading to an extended $s$-wave pairing SC with high $T_c$. Meanwhile, the strongly localized $3d_{z^2}$ electrons tend to form interlayer rung singlets. Due to a lack of phase coherence, these singlets do not directly participate in the SC condensate, but instead give rise to a pseudogap phase.

The discovery of superconductivity with Tc ~ 80 K in the nickelate Ruddlesden-Popper bilayer La3Ni2O7 at high pressure has opened a new platform for unconventional superconductivity, followed by the subsequent observation of superconductivity in trilayer La4Ni3O10, also at high pressure. Remarkably, ambient-pressure superconductivity was also observed recently in La3Ni2O7 ultra-thin films when grown on substrates that provide compressive strain. This discovery significantly extends the type of experimental techniques that can be used in nickelates, previously limited due to the high-pressure constraint. Discussing the similarities and differences among these nickel oxides will provide new insights into understanding the mechanism of high-Tc superconductivity in correlated electron systems. In this paper, we review the experimental and theoretical progress on RuddlesdenPopper nickelates, with emphasis on thin films, and discuss future perspectives and research directions.

Ruddlesden-Popper nickelates offer a new route to high-temperature superconductivity beyond the cuprates and iron-pnictides. However, the electronic reorganization that enables superconductivity in bilayer nickelates remain unresolved, largely due to the difficulty of directly probing the superconducting phase. Here, we overcome this limitation by stabilizing superconducting (La,Pr)$_3$Ni$_2$O$_7$ thin films with a protective capping layer, thereby enabling direct spectroscopic access via X-ray absorption and resonant inelastic X-ray scattering. We resolve the evolution of in-plane and out-of-plane electronic states, spin and orbital excitations, and spin-density-waves across insulating, superconducting, and metallic regimes. Combining experimental results with theoretical analysis, we show that the in-plane $d_{x^2-y^2}$ states form an itinerant backbone, while superconductivity emerges only when coherent $d_{z^2}$-$p_z$-$d_{z^2}$ interlayer hybridization develops, accompanied by suppressed static spin order and strongly damped spin excitations. Oxygen stoichiometry and epitaxial strain both act on this interlayer channel, placing superconductivity within a narrow window of interlayer coherence and correlation strength. These findings identify the microscopic ingredients required for superconductivity in bilayer nickelates and provide a multiorbital picture of its emergence.

Based on the discovery of high-temperature superconductivity in the bilayer nickelate La$_3$Ni$_2$O$_7$, several Co-based La$_3$Ni$_2$O$_7$-like materials were theoretically predicted as possible high-temperature superconductors by electron doping. Motivated by these findings and the subsequent discovery of superconductivity in the trilayer nickelate La$_4$Ni$_3$O$_{10}$ under high pressure, we propose and investigate a Co-based La$_4$Ni$_3$O$_{10}$-like material. With electron doping to the high-pressure trilayer cobaltate La$_4$Co$_3$O$_{10}$, using density functional theory combined with dynamical mean-field theory (DFT+DMFT), we find that the resulting compound La$_4$Co$_2$NiO$_8$Cl$_2$ exhibits a crystal structure and a strongly correlated electronic structure similar to those of La$_4$Ni$_3$O$_{10}$ under high pressure. This suggests that this new compound may host high-temperature superconductivity.

The interplay of multiorbital physics and nonlocal self-energy effects is studied within an effective three-orbital model for the high-pressure normal state of superconducting bilayer nickelate La$_3$Ni$_2$O$_7$. The model is solved within an advanced many-body framework capturing $k$-dependent correlations beyond dynamical mean-field theory. Different low-energy scenarios subtly depend on the strength of the interorbital interaction, either placing the notorious flat $γ$ quasiparticle band in the occupied part of the spectrum, or letting it cross the Fermi level. In the latter case, intriguing spin-polaron formation due to the scattering of electrons with paramagnon excitations takes place. This leads to bound states appearing as a shadow band with incoherent low-energy spectral weight below the Fermi level. Our results uncover additional competing states that exist in bilayer nickelates and could explain the controversy of recent angle-resolved photoemission experiments.

Beyond the quasi-two-dimensional (2D) paradigm of cuprates, the role of the third dimension of the Ruddlesden-Popper bilayer nickelates is essential to decoding their superconducting mechanism. Here, using angle-resolved photoemission spectroscopy (ARPES) with varied photon energies, we systematically investigate the electronic band structures in three dimensions for superconducting (La,Pr,Sm)$_3$Ni$_2$O$_7$/SrLaAlO$_4$ thin films (superconducting onset temperature $T_c^{\text{onset}} \sim 48$ K) transferred via a cryogenic ultra-high vacuum suitcase. We reveal an orbital-dependent dimensionality: while the $d{x^2-y^2}$-dominant bands exhibit a quasi-2D character, the $d{z^2}$-dominant band displays a finite $k_z$ dispersion. Finite energy gaps are identified on all observed bands across multiple high-symmetry directions. Systematic temperature-dependent analysis characterizes the superconducting nature of the gap on the $d{z^2}$-derived band, revealing a large gap $Δ\sim 18$ meV and a ratio $2Δ/k_BT_c\sim 8$ exceeding the weak-coupling BCS limit. The suppression of spectral weight near the Fermi level persists above the superconducting transition temperature. Ubiquitous waterfall-like spectral features evidence the presence of electron interactions. These results underscore the role of the $d_{z^2}$ orbital and correlations, placing constraints on theoretical models for nickelate superconductivity.

The recent discovery of Ruddlesden-Popper (RP) nickelate thin-film superconductors has opened a new frontier in unconventional superconductivity. Its realization requires both compressive epitaxial strain and highly oxidative growth conditions, yet the microscopic pathway from the parent phase to the superconducting phase remains elusive. Here, X-ray absorption spectroscopy and resonant inelastic X-ray scattering are employed to track this evolution by independently tuning strain and oxygen content in (La,Pr)$_3$Ni$_2$O$_{7-δ}$ thin films. We uncover a remarkable two-step evolution. First, signatures of delocalization emerge in the same way upon two independent tunings: spectral weight transfers from an upper Hubbard-like peak to the hole-like peak associated with the O $2p_z$ state, and in parallel, the initially localized Ni $3d_{z^2}$ orbital becomes more itinerant, accompanied by the broadening and weakening of $dd$ orbital excitations. Second, as itinerancy increases, long-range spin-density-wave (SDW) order is suppressed in both intensity and correlation length, indicating direct competition with superconductivity. Yet, short-range magnons persist: they become damped, but their bandwidth remains unchanged. Our results paint a coherent picture in which both strain and oxygenation drive the RP bilayer nickelates towards superconducting instability, where the O $2p_z$ and Ni $3d_{z^2}$ orbitals become delocalized. Concomitantly, the long-range magnetic order loses coherence and is suppressed. These findings establish an orbital-selective route to RP nickelate superconductivity, in which the delocalization of the interlayer $3d_{z^2}$-$2p_z$-$3d_{z^2}$ molecular orbital and the robust short-range magnons upon the melting of SDW order are prerequisites, providing strong constraints for theory and a roadmap for designing nickelate superconductors.

The recent discovery of Ruddlesden--Popper (RP) nickelate thin-film superconductors has opened a new frontier in unconventional superconductivity. Its realization requires both compressive epitaxial strain and highly oxidative growth conditions, yet the microscopic pathway from the parent phase to the superconducting phase remains elusive. Here, X-ray absorption spectra and resonant inelastic X-ray scattering are employed to track this evolution by independently tuning strain and oxygen content in (La,Pr)$_3$Ni$_2$O$_{7-δ}$ thin films. We uncover a remarkable two-step narrative. First, signatures of delocalization emerge in the same way upon two independent tunings: Spectral weight transfers from a ''Upper Hubbard''-like peak to the hole-like peak associated with O $2p_z$ state, and in parallel, the initially localized Ni $3d_{z^2}$ orbital becomes more itinerant followed by the broadening and weakening of $dd$ orbital excitations. Second, as itinerancy increases, long-range spin-density-wave (SDW) order is suppressed in both intensity and correlation length, indicating direct competition with superconductivity. Yet, short-range magnons persist: they become damped but their bandwidth stays unchanged. Our results paint a coherent picture that both strain and oxygenation drive the RP bilayer nickelates towards the superconducting instability, where the O $2p_z$ and Ni $3d_{z^2}$ orbitals become delocalized. Concomitantly, the long-range magnetic order loses coherence and gets suppressed. These findings establish an orbital-selective route to RP nickelate superconductivity, in which the delocalization of the $2p_z$ and $3d_{z^2}$ orbitals and the robust short-range magnons upon the melting of SDW order are prerequisites, providing strong constraints for theory and the roadmap for designing nickelate superconductors.

Inspired by the possibility that superconducting properties may be altered by applying a perpendicular electric field in the Ruddlesden-Popper (RP) bilayer nickelate La$_3$Ni$_2$O$_7$, we investigated the imbalanced two-orbital bilayer Hubbard model using dynamical cluster quantum Monte Carlo calculations. Focusing on the pairing symmetries induced by the electric field and their evolution with field strength in the undoped, hole-doped, and electron-doped regimes, we found that the $s^\pm$-wave pairing originating from the $d_{z^2}$ orbital is suppressed; while a pairing symmetry transition from $s^\pm$-wave to $d$-wave pairing occurs, driven by the interlayer $d_{z^2}$ orbital mismatch and the transfer of electrons into the $d_{x^2-y^2}$ orbital under the applied electric field. Intriguingly, the $d$-wave pairing arising from the $d_{x^2-y^2}$ orbital exhibits dome-like behavior with the electric field. Our large-scale many-body calculations align with the previous expectation from weak-coupling methods and provide further insight into the superconducting mechanism in RP nickelates.

Hydrogenation is known to tune superconductivity in a wide range of materials. While its microscopic role has been clarified in phonon-mediated superconductors such as hydrogenated MgB2, LaH10, and H3S, much less is known for hydrogenated cuprates and iron-based superconductors, where even the underlying structural motifs remain elusive. Using hydrogenated FeSe as a prototypical example, we reveal how hydrogen affects superconductivity in the presence of strong electronic correlations: correlation-induced orbital renormalization shifts hydrogen-derived spectral weight from the high-energy region toward the Fermi surface (FS), remarkably enhancing the electron-phonon coupling (EPC). We predict a structurally stable FeSeH phase where, compared to bare FeSe, hydrogen incorporation reshapes the FS topology and increases the number of channels for electron-phonon scattering, while simultaneously introducing high-frequency phonons that strengthen pairing. First-principles EPC calculations combined with dynamical mean field theory (DMFT) yield a superconducting transition temperature (Tc) exceeding 40 K. Fully anisotropic Eliashberg theory reveals a two-gap superconducting state, consistent with the gap structure experimentally observed in doped FeSe. Our findings identify correlation-enhanced EPC as a plausible microscopic mechanism for iron-based superconductivity and offer a new perspective on pairing in strongly correlated systems. In addition, this work establishes hydrogenated FeSe as a promising platform for engineering two-dimensional superconductors and superconducting quantum devices.

The interplay between charge and spin dynamics is at the heart of strongly correlated materials. Inspired by recent quantum simulation experiments, we develop a conserving diagrammatic method to describe the Fermi-Hubbard model for strong repulsion and small hole doping away from the half-filled anti-ferromagnetic ground state. We show that doping leads to four hole pockets in the Brillouin zone formed by magnetic polarons, which become increasingly damped with hole concentration. Likewise, the magnon spectrum of the anti-ferromagnet softens and dampens with doping due to hole-induced magnetic frustration. This gives rise to a suppression of the anti-ferromagnetic correlations in agreement with recent experiments. We then calculate the response of the system to a lattice modulation and recover the qualitative difference between in-phase and out-of-phase modulations seen in experiments, which was interpreted as signs of pseudogap physics. Our results indicate that the complex competition between spin and charge degrees of freedom and the emergence of the pseudogap phase may be usefully analyzed for small dopings, where systematic theories can be developed.

Altermagnets (AMs), characterized by zero net magnetization and momentum-dependent spin splitting, are anticipated to hold significant potential for generating multiple exotic and uncommon superconducting states. However, superconductivity has not yet been realized in AMs to date. Recently, two-dimensional (2D) V$_2$Ch$_2$O (Ch = Se, Te) monolayers, as well as AV$_2$Ch$_2$O (A = K, Rb, Cs) crystals containing [V$_2$Ch$_2$O]$^{δ-}$ building layers, have been predicted and/or demonstrated to be promising altermagnetic materials. Our preliminary attempts to explore superconductivity in these materials by applying pressure or chemical doping were unsuccessful. Here we report the discovery of superconductivity at a relatively high transition temperature of ~ 16.3 K in a newly synthesized layered compound, Na$_{2-x}$V$_2$Se$_2$O, a variant of AV$_2$Ch$_2$O. In this structure, the [V$_2$Ch$_2$O]$^{δ-}$ layers are interspersed with double layers of Na$^+$ instead of a single layer of A$^+$, with sodium sites being only half-filled. This new family of layered vanadium oxychalcogenides, lacking inversion symmetry, represents an intriguing platform for exploring altermagnetic superconductors, and holds the potential to reveal novel phenomena, such as topological states, van Hove singularities, and finite-momentum superconductivity. Furthermore, this material acts as a "bridge" between the cuprate/nickelate and iron-pnictide high temperature superconductors, providing new hope and opportunity to expand the category of layered superconductors with higher critical temperatures (T$_c$) and enhancing our understanding of the underlying mechanisms in these systems.

The Grassmann time-evolving matrix product operator (GTEMPO) method, which represents the Feynman-Vernon influence functional as a temporal matrix product state, has been shown to be a flexible and potentially scalable solution for fermionic quantum impurity problems. In this work, we extend GTEMPO to solve fermionic impurity problems in the Nambu formalism, in which the impurity is coupled to a superconducting bath. A key insight is that by employing the Bogoliubov transformation for the superconducting bath, one could obtain the analytic expression of the Feynman-Vernon influence functional in a similar form to the case of a normal bath, after which the core algorithms of GTEMPO can be straightforwardly adapted. We demonstrate the accuracy of our method by benchmarking it against exact diagonalization in several exactly solvable cases, and against the continuous-time quantum Monte Carlo method using converged dynamical mean field theory (DMFT) iterations on the imaginary contour in the non-integrable case. In all cases, we perform both imaginary- and real-time calculations to illustrate the flexibility of our method. These results illustrate that our method could be potentially useful as an impurity solver in DMFT as well as its non-equilibrium extension for fermionic impurity problems in the Nambu formalism.

We investigate the ground-state properties of correlated three-leg ladders near 1/3 filling. We apply the density-matrix renormalization group method to the three-leg t-J ladder with strong interchain couplings and evaluate its pairing nature. When holes are doped into the spin-gapped state at 1/3 filling, we find that pair correlations develop with power-law decays while spin correlations decay exponentially. On the other hand, doping of electrons into the 1/3-filled state does not give rise to substantial pair correlations. We also discuss the hole-doped state in the three-leg Hubbard model to compare it with the pairing state in the t-J model. Our numerical demonstrations provide insights into the electronic properties of trilayer nickelate superconductors.

The discovery of ambient-pressure superconductivity in compressively strained (La,Pr)$_3$Ni$_2$O$_7$ thin films has intensified efforts to identify the pairing mechanism. However, the symmetry of the superconducting order parameter and the character of the normal state remain unsettled. Here we combine bulk-sensitive terahertz (THz) time-domain spectroscopy with THz third-harmonic generation to present spectroscopic insights into these issues. Linear THz spectroscopy reveals a bulk superconducting response in the (La,Pr)$_3$Ni$_2$O$_7$ films, evidenced by the suppression of low-frequency spectral weight below the onset critical temperature, $T_\mathrm{c}^{\mathrm{onset}}$. A weak coherence peak near $T_\mathrm{c}^{\mathrm{onset}}$, together with substantial residual low-frequency conductivity as $T\to 0$, is consistent with disordered $s_{\pm}$-wave pairing. In the nonlinear regime, the third-harmonic signal rises sharply on cooling through $T_\mathrm{c}^{\mathrm{onset}}$, providing an independent signature of the transition. Strikingly, the nonlinear response persists above $T_\mathrm{c}^{\mathrm{onset}}$, pointing to either disorder-enhanced nonlinearity or a distinct correlated normal state. Motivated by angle-resolved photoemission spectroscopy on similarly grown films that identifies a comparable temperature scale, we associate the anomalous normal-state terahertz nonlinearity with a pseudogap. These results establish (La,Pr)$_3$Ni$_2$O$_7$ as a bulk superconductor with $s_{\pm}$-like pairing that coexists with, and may compete with, a distinct ordered state, providing a platform for exploring unconventional superconductivity beyond cuprates and pnictides.

We present a two-channel extension of the Allen-Dynes framework that unifies phonon-mediated and spin-fluctuation-mediated pairing channels for predicting superconducting critical temperatures. Channel 1 employs the standard Allen-Dynes formula with material-specific electron-phonon coupling; Channel 2 incorporates a spin-fluctuation coupling parameter extracted from inelastic neutron scattering data. Blind predictions for 19 materials spanning conventional superconductors, MgB2, iron pnictides, iron chalcogenides, heavy fermions, cuprates, and hydrides achieve R-squared = 0.96 across five orders of magnitude in Tc (0.4-250 K) without free parameters. We further demonstrate a quantum-metric no-go result: the Peotta-Torma geometric superfluid weight, while essential for flat-band systems, cannot serve as a universal predictor of Tc because it correlates with band-structure topology rather than pairing strength. The framework identifies the spin-fluctuation channel as the dominant contributor to Tc enhancement in unconventional superconductors, providing quantitative design rules for materials with Tc above 100 K.

Ruthenium trihalides Ru$X_3$ ($X$ = Cl, Br, I) provide a tunable platform for Kitaev magnetism in two-dimensional van der Waals materials. Despite their similar crystal structures and zigzag antiferromagnetic order, RuBr$_3$ exhibits a higher Néel temperature ($T_N$) than RuCl$_3$, suggesting their distinct proximity to the Kitaev quantum spin liquid phase. Using Ru $L_3$-edge resonant x-ray scattering, we show that, while the long-range zigzag order in RuBr$_3$ disappears at $T_N$, the zigzag correlations that persist well above $T_N$ show a pronounced spectral weight redistribution along the interlayer direction. These results suggest that the enhanced interlayer magnetic interactions driven by the extended Br 4$p$ orbitals stabilize three-dimensional zigzag correlations in RuBr$_3$.

Strongly-correlated electrons in transition-metal oxides give rise to intriguing emergent phenomena, including high-temperature superconductivity in cuprates. While simplified one-band Hubbard models capture some aspects, explicitly describing the interplay of copper and oxygen orbitals -- as in the three-band Emery model -- is essential to capture the full phenomenology of cuprates. Quantum simulators based on ultracold atoms offer a promising route to study such systems in a controlled setting, but realizing realistic multi-orbital Hubbard models remains challenging. Here we propose an optical superlattice architecture that implements the three-band Emery model with ultracold fermions. By combining lattice beams with controllable interference, we engineer orbital degrees of freedom that reproduce key features of the cuprate band structure, while enabling independent control of orbital-dependent interactions and charge-transfer energy. We show that single-particle quantum walks can benchmark the resulting tight-binding model. Using determinant quantum Monte Carlo, we further investigate thermodynamic properties in the undoped regime and find a finite-temperature metal-insulator crossover accompanied by the onset of antiferromagnetic correlations accessible in current experiments. Finally, we apply a Hamiltonian learning protocol enabling to infer effective single-band Hubbard models from experimental realizations of Emery models. Our results provide a practical pathway to simulate multi-orbital Hubbard physics with quantum gas microscopes.

The orbital degree of freedom in rare-earth nickelates is typically inactive across the temperature-driven metal-insulator transition, where the system develops two inequivalent Ni sites associated with Ni-O bond disproportionation and breathing-mode distortions of NiO$_6$ octahedra. Here, we show that orbital polarization can be induced by optical excitation with linearly polarized light. Using an interacting multiband tight-binding model combined with real-time simulations of coupled electron-ion-spin dynamics, we find that photoinduced $d$-$d$ transitions reduce the local magnetic moments at Ni sites and effectively suppress Hund's coupling $J$ in the excited state. Importantly, these transitions can be made strongly orbital-selective by tuning the light polarization, leading to an imbalance in $e_g$ orbital occupancies. The resulting nonequilibrium state, characterized by reduced effective $J$ and unequal orbital populations, becomes unstable toward Jahn-Teller (JT) distortions, driving structural relaxation along coherently excited JT modes. Our results demonstrate that polarization-controlled optical excitation provides a pathway to access hidden nonthermal phases with emergent orbital order, enabling coherent control of coupled charge, spin, and lattice degrees of freedom on ultrafast timescales.

Accurately describing strongly correlated electrons in systems such as transition metal complexes requires strict adherence to spin symmetry, a feature largely absent in modern neural-network-based variational wavefunctions. This deficiency can lead to severe spin contamination in simulating systems with near-degenerate spin states. To resolve this limitation, we present a spin-adapted neural network backflow (SA-NNBF) ansatz, formulated in second quantization for fermionic lattice models and ab initio quantum chemistry. Our approach constructs a fully antisymmetric wavefunction by combining a neural-network backflow spatial component with a spin eigenfunction expressed in a sum-of-products form. To address the computational complexity of spin adaptation, we introduce a tensor compression algorithm for spin eigenfunctions, and a more compact wavefunction representation based on the particle-hole duality in second quantization. These advancements enable variational Monte Carlo calculations using SA-NNBF for challenging molecular systems with more than one hundred electrons, including the FeMo-cofactor (FeMoco) in nitrogenase. Applications to prototypical strongly correlated molecules demonstrate that SA-NNBF consistently outperforms standard NNBF with a similar number of parameters. Furthermore, it surpasses the accuracy of the state-of-the-art spin-adapted density matrix renormalization group (SA-DMRG) algorithm for FeMoco with a significantly reduced computational resource. Our work establishes a foundational framework for exploring fully symmetry-preserving neural-network quantum states for interacting fermion problems.

Altermagnetism generates exchange-type spin splitting without net magnetization and, in its $\it d$-wave form, resembles the angular symmetry of unconventional $\it d$-wave superconductivity. Whether this correspondence bears directly on superconducting instabilities in real correlated materials remains open. Here we study the quasi-two-dimensional vanadium oxychalcogenide CsV$_2$Se$_2$O (CVSO), a square-net $\it d$-wave altermagnet candidate, through combined experimental and theoretical investigation of its lattice structure, electronic structure and transport properties. At ambient pressure, CVSO is a weakly insulating parent state with a density-wave-like anomaly near 100 K, and its bulk properties are most consistent with a G-type compensated antiferromagnetic background. Under compression, the density-wave-like feature is suppressed, the magnetoresistance evolves from predominantly negative to positive, and a superconducting-like resistive downturn emerges below about 3 K. This low-temperature anomaly is reproducible across samples and pressure media, and is suppressed by magnetic field. Room-temperature X-ray diffraction reveals no symmetry lowering, whereas does show a pronounced compressibility anomaly over the same pressure range. CVSO thus reveals a pressure-tuned phase diagram in which a reconstructed weakly insulating parent state gives way to strange-metal-like transport and superconducting-like behavior, echoing broader phenomenology associated with unconventional superconductors, including cuprates and nickelates.

We develop a theoretical framework for Floquet resonant X-ray scattering, using Floquet theory combined with the ultrashort core-hole lifetime expansion. We obtain a compact expression for the Floquet components of the resonant inelastic X-ray scattering operator, which shows that Floquet X-ray scattering provides direct access to bond and current correlations that do not directly produce charge Bragg peaks in conventional diffraction. Applying this framework to charge-ordered states on the Kagome lattice, we demonstrate that different symmetry-breaking orders exhibit distinct polarization fingerprints in the Floquet Bragg peaks. Moreover, the relative weight of bond and current contributions can be tuned through the drive frequency. These results establish Floquet X-ray scattering as a symmetry-resolved probe of hidden electronic order or fluctuations in quantum materials.

We construct an effective single-particle Hamiltonian $K_{\mathrm{eff}}$ from Monte Carlo--averaged matrix logarithms of the imaginary-time propagator in determinant quantum Monte Carlo (DQMC). The logarithm maps the multiplicative sign problem into an additive framework where the central limit theorem guarantees convergence, rendering $K_{\mathrm{eff}}$ sign-problem-free: both sign sectors yield identical dispersions to $<1\%$. $K_{\mathrm{eff}}$ captures the exact correlated single-particle spectrum, incorporating all self-energy effects non-perturbatively. Applied to the Hubbard model ($t'/t = -0.30$, $U/t = 4$), $K_{\mathrm{eff}}$ reveals a $d$-wave pseudogap with strong nodal-antinodal dichotomy below a computational phase transition at $T^*$. Three sign-free observables provide evidence consistent with spin-fluctuation pairing: (i) the gap ratio $R_g > 1$ confirms $d$-wave symmetry -- a temperature-independent property of the correlated band structure that provides the medium for pairing; (ii) the superfluid stiffness $ρ_s$ forms a dome across doping at $L = 8$, $10$, and $12$, exceeding the Berezinskii-Kosterlitz-Thouless threshold by $5$-$7\times$ at the dome peak; (iii) $S(π,π)$ is approximately flat across doping, establishing that the dome originates from Fermi-surface geometry responding to uniform spin-fluctuation glue. The pseudogap grows monotonically toward half-filling while $ρ_s$ forms a dome, mirroring cuprate phenomenology where $T_c$ is limited by the superfluid density (Uemura relation). Vertex corrections remain to be quantified.

Odd-parity altermagnetism(ALM) extends compensated collinear magnetism beyond the even-parity spin splitting of conventional altermagnets, but its role in correlated topological phases remains largely unexplored. Using the cluster slave-spin method, we show that the odd-parity ALM appearing in the ALM Chern-insulating phase of Haldane-Hubbard model significantly reconstructs the local topology in the conventional Chern-insulating phase, while the total Chern number remains unchanged compared to the Chern-insulating phase. The Berry curvature becomes spin and valley selective; zigzag ribbons develop chiral-symmetry-breaking edge states; while armchair ribbons remain inversion symmetric. The optical response mirrors this separation between the local reconstruction and the global topology: low-energy spectra are governed by quasiparticles near the gap, whereas the low-frequency Hall conductivity stays quantized, $σ_{\rm T\uparrow}(Ω\to 0)=σ_{\rm T\downarrow}(Ω\to 0)=e^2/h$. These results establish the Haldane-Hubbard model as a minimal correlated platform for odd-parity altermagnetic topology.

Motivated by the continued interest in Fermi-surface reconstruction without symmetry breaking, we present an analytically tractable microscopic model of a fractionalized Fermi liquid (FL$^*$) on a square lattice and discuss its potential relevance to the cuprates. As in ancilla-qubit constructions, the model is related to Kondo lattice systems, but in this case, the conduction electrons interact with a $\mathbb{Z}_2$ spin liquid of the Yao--Lee type, with a Majorana Fermi surface. The associated $\mathbb Z_2$ gauge theory is static so that the model can be analytically solved to leading-logarithic accuracy. There are two phases: one in which the fractionalized fermions of the spin liquid hybridize with conduction electrons to form a common Fermi surface violating the naive Luttinger count, and one in which they remain decoupled. We discuss the salient features of the small Fermi-surface phase, including analytically derived momentum dependent coherence factors responsible for the appearance of Fermi arcs à la Yang-Rice-Zhang. We further discuss the impact of quantum and thermal fluctuations, including a strong diamagnetic response and a logarithmically divergent Sommerfeld coefficient at the onset of the pseudogap.

Maximally-localized Wannier functions are quantum wavefunctions resembling atomic orbitals that are used to describe electrons in condensed matter. Since their introduction in 1997, these functions have become ubiquitous in ab initio materials simulations, including applications in linear-scaling methods, strongly-correlated electron systems, quantum transport, electron-phonon interactions, and topological materials. Despite their widespread adoption in a vast software ecosystem, Wannier functions have not yet attained their fullest potential in the presence of entangled bands, as their optimization remains challenging and labor-intensive. Here, we introduce a universal meta-optimization method that leverages workflow abstraction and machine learning techniques like differential evolution and Bayesian optimization to generate globally optimized Wannier functions without human intervention. We demonstrate this approach through three applications: (i) autonomous interpolation of entangled band structures with millielectronvolt accuracy starting from coarse Brillouin zone grids, (ii) thousand-fold acceleration of fully ab initio Boltzmann transport calculations via the use of minimal coarse Brillouin zone grids, and (iii) ultra-fast high-throughput calculations of high-precision Wannier functions for large materials libraries. This work brings calculations that previously required supercomputers within the reach of personal computers.

Recent years have seen a growing interest in the use of belief propagation - an algorithm originally introduced for performing statistical inference on graphical models - for approximate, but highly efficient, tensor network contraction. Here, we detail how to apply generalized belief propagation (GBP) - where messages are passed within a hierarchy of overlapping regions of the tensor network - to approximately contract tensor networks and obtain accurate results. The original belief propagation algorithm is a corner case of this approach, corresponding to a particularly simple choice of regions of the tensor network. We implement the GBP algorithm for a number of different region choices on a range of two- and three-dimensional, infinite and finite tensor networks, solving the corresponding fixed point equations both numerically and, in certain tractable cases, analytically. Our examples include calculating the partition function of the fully frustrated Ising model, computing the ground state degeneracy of three-dimensional ice models, measuring observables on the deformed AKLT quantum state and evaluating the norm of randomly generated tensor network states.

By employing a combined method of density functional theory and dynamical mean field theory (DFT+DMFT) we investigate the effect of electronic correlations on the magnetic and superconducting properties of the iron-based parent compound LaFeAsO. We find that the static non-local susceptibility $χ({\bf q})$ exhibits a peak at the in-plane wave vector ${\mathbf Q}=(π,π)$, which is strongly enhanced upon inclusion of the vertex corrections in the ladder approximation, leading to magnetic instability. Considering the eigenfunctions of the Bethe-Salpeter equation with the vertex, obtained within the second order perturbation theory, as well as the ladder approach containing dynamic interaction vertices, in agreement with earlier weak-coupling-based studies of LaFeAsO, we obtain a close competition between $d$-wave and $s_{\pm}$ order parameters, dominating in the second-order and ladder approach, respectively. We argue that the dominating $s_{\pm}$ instability in the ladder DFT+DMFT approach is related to the reduced degree of magnetic frustration by itinerant degrees of freedom due to only partially formed local magnetic moments. Our study shows that dynamic correlation effects do not change the type of the leading superconducting instability in LaFeAsO.

Topochemical fluorination offers a low--temperature route for modifying the anion chemistry and electronic ground states of layered transition-metal oxides, providing access to metastable phases and functionalities that are not able to be achieved through conventional solid--state synthesis. Despite extensive work on polycrystalline samples and thin films, topochemical fluorination of bulk single crystals has not been studied, limiting insights into intrinsic structure property relationships. Here, we investigate the topochemical fluorination of optical float zone grown (OFZ) La$_2$NiO$_{4+δ}$ single crystals using polymer-based PTFE, PVDF and inorganic CuF$_{2}$ fluorination agents and compare it to our topochemical pathways of reduction of LaNiO$_{3-x}$. By systematically investigating direct and indirect contact reaction pathways, we can understand fluorination mechanisms, quantify the degree of fluorine incorporation, and evaluate the resulting structural and magnetic modifications in a detail that was not possible in powder and thin films. Powder and single--crystal X-ray diffraction reveal that fluorination proceeds without destroying the Ruddlesden--Popper framework, while inducing lattice parameter changes consistent with anion intercalation in the bulk and ion exchange on the surface. This even induces a clear superstructure, which was not reported before and extends the understanding of anion insertion reactions beyond what is known on stage ordering in nickelates. Energy-dispersive X--ray spectroscopy confirms strong fluorine incorporation on the surface and reduced homogeneity in the bulk. Magnetic susceptibility measurements demonstrate a change in antiferromagnetic ordering upon fluorination.

Transition metal dichalcogenides exhibit a wide range of semiconducting, metallic, correlated, and topological electronic states that arise from strong coupling between lattice structure, dimensionality, and electronic degrees of freedom. In these materials, crystal growth is not merely a preparative step but a thermodynamic boundary condition that establishes chemical potentials, defect populations, polytype stability, and access to metastable phases. As a result, synthesis determines the structural and defect landscape from which collective electronic behavior emerges. In this Review, we develop a unified thermodynamic and kinetic framework that connects growth conditions to phase stability, defect energetics, and microstructure. We examine how chemical potential constraints define stability windows, how supersaturation and mass-transport regimes govern nucleation and morphology, and how nonequilibrium pathways enable kinetic trapping and polymorph selection. Bulk and thin-film synthesis approaches, including chemical vapor transport, flux growth, physical vapor transport, solvent-assisted crystallization, chemical vapor deposition, and molecular beam epitaxy, are placed within a common thermodynamic and kinetic map to clarify how distinct growth regimes produce characteristic disorder profiles and structural phases. By explicitly linking synthesis variables to charge-density-wave order, superconductivity, band topology, and correlation effects, this Review establishes crystal growth as a central parameter in determining the effective electronic Hamiltonian realized in experiment. This perspective provides a physically grounded framework for improving reproducibility and guiding deterministic control of emergent quantum phases in layered materials.

Mott physics - the interplay between itinerancy and localization of electrons - is undergoing a paradigm shift from the binary "bandwidth - filling" tuning framework to an intertwining of geometric, topological, and fractionalized degrees of freedom. Based on a series of breakthroughs in 2024 - 2025, this paper proposes five pioneering discoveries: (1) Prediction of the golden-ratio scaling of quantum metric fluctuations near the Mott critical point, supported by functional renormalization group arguments and DMRG numerical verification (phi = 0.618 +/- 0.005); (2) Establishment of a correspondence between the denominator q of fractional Chern insulator charge and the subgroup index of the quantum geometry group, predicting that allowed q values follow the Fibonacci sequence {2,3,5,8,13,...} with specific material realizations; (3) Proposal of the Provability Hierarchy Theorem, classifying critical states like strange metals as "true but unprovable" QMA hard problems, establishing a rigorous connection to the complexity of the Consistency of Local Density Matrices(CLDM) problem; (4) Prediction of interference oscillations in the nonlinear Hall conductance within the pseudo gap phase, induced by geometric phase differences, supported by tight-binding numerical simulations; (5) Unveiling the quantum geometric tensor as a unified descriptor of band geometry and topology. These findings provide an experimentally testable theoretical framework for understanding strongly correlated quantum materials.

We study a Su-Schrieffer-Heeger chain coupled to a single mode photonic cavity. Considering an off-resonant regime we use the high-frequency expansion in order to obtain an effective fermionic Hamiltonian with cavity-mediated interactions. We characterize the effects of the cavity on topology in a finite size chain by studying three different markers adapted for interacting systems: correlation functions between edges in a chain with open boundary conditions, and a winding number based on the single-particle Green's function and bulk electric polarization via the many-body formula by Resta for a chain with periodic boundary conditions. There is excellent agreement between the winding number and polarization approaches to compute the phase diagram, with the presence of the edge states being confirmed through the calculations of the two-point correlation function. Our approach provides an alternative perspective on cavity-modified topological phases through a study of an effective interacting electronic Hamiltonian and complements methods that treat the full light-matter Hamiltonian directly.

We use segment DMRG on the infinite cylinder to compute energies of charged excitations in gate-screened fractional quantum Hall states. For the $ν=1/3$ Laughlin, $ν=2/5$ Jain, and $ν=5/2$ anti-Pfaffian states, we find screening can bind like-charged anyons into molecules, with a strong dependence on filling-factor, gate distance, and fusion channel. In the Laughlin state, stable $\pm 2e/3$ molecules and larger clusters appear over a broad gate-distance window. The Jain state is molecular throughout the range we consider. In the anti-Pfaffian, binding is strongest on the hole side, where the charge-$e/2$ molecule is fused into the $ψ$ channel over a broad window of gate distances. In all three cases, screening suppresses long-range repulsion and exposes an intermediate-range attraction encoded in the oscillatory density tail of the fundamental anyon. We discuss consequences for addition spectra, interferometry, Wigner crystallization, anyon superconductivity, and entropy measurements.

We investigate superconductivity in a two-band $t$-$J$ model consisting of an itinerant orbital (orbital-0) and a quasi-localized orbital (orbital-1) using variational Monte Carlo. A robust orbital-selective $d$-wave superconducting state is found to emerge exclusively from the itinerant orbital. An analysis of the superexchange energy hierarchy shows that the quasi-localized orbital-1 competes with superconductivity by favoring local inter-orbital bound states, which act as energy defects and disrupt phase coherence. Consistently, the superconducting order parameter is monotonically suppressed as the occupancy of orbital-1 increases. Motivated by superconductivity in nickelate La$_3$Ni$_2$O$_7$, these results highlight the essential role of multi-orbital physics beyond the single-band $t$-$J$ framework and point to a concrete route to enhance $T_c$: suppressing the involvement of localized $d_{z^2}$-derived orbitals.

We study the spin excitation spectra in magnetically ordered phases proximate to the Kitaev quantum spin liquid (KQSL). Although the low-energy universal features should be governed by the magnetic orders, the $\textit{non-universal}$ high-energy features of the KQSL and adjacent phases can be remarkably similar. Therefore, we study the extended Kitaev model within a Stoner-like theory using Majorana partons, and compute the inelastic neutron scattering (INS) intensities in the random phase approximation. First, we benchmark against the antiferromagnetic (AFM) Heisenberg model and recover the AFM order with linear Goldstone modes. We then explore the phase diagram which agrees qualitatively with previous numerical results. In particular, the Majorana parton theory accurately captures Order-by-Disorder effects in the Kitaev-Heisenberg limit. We also find large INS intensities near the associated high-symmetry Brillouin zone (BZ) points of the magnetic orders. At intermediate and high energies, broad multi-spinon continua emerge across the BZ, providing a distinct mechanism for magnon decay and spectral broadening beyond the conventional multi-magnon decay scenario. Finally, we study the model Hamiltonian of candidate Kitaev material $α$-RuCl$_3$. The zigzag ground state agrees qualitatively with experiments, its stability under external magnetic field also exhibits strong anisotropy in the field directions, and broad scattering continua are recovered similar to those observed experimentally.

We investigate the competition between attractive spin-spin interactions and spin-separating external forces in the ground state of a one-dimensional Fermi-Hubbard model. We consider a lattice with open boundary conditions, subject to a linear external potential whose gradient is opposite for the two spin components, so that each spin species sees a potential minimum at a different end of the lattice. Using density-matrix renormalization group (DMRG) simulations, we map the ground-state density distributions and the number of doubly occupied sites as a function of the potential gradient $β$ and interaction strength. We identify three distinct regimes separated by critical threshold gradients: (i) a small-$β$ regime where fermion pairing remains robust against the external potential; (ii) an intermediate-$β$ phase-separated regime characterized by a staircase-like decrease in the doublon number, corresponding to the successive, one-by-one breaking of bound pairs; and (iii) a large-$β$ regime where the two spin components are completely spatially separated. We complement the numerical results with a phenomenological model and a local-density approximation analysis, from which we derive closed-form analytical estimates for these critical threshold values. We also verify that the staircase structure persists under additional harmonic confinement. Our results are directly testable in cold-atom experiments, and demonstrate that a spin-dependent linear potential enables precise, integer-level control of the number of bound fermion pairs.

The two-dimensional Hubbard model at finite doping hosts competing or intertwined orders, resulting in conflicting conclusions from different computational approaches regarding its ground state. We show that a key source of such discrepancies is the bias encoded in the variational ansatz. We consider three different Transformer backflow fermionic wave functions based on a Slater determinant, its particle-hole counterpart, and a Pfaffian, initialized without any mean-field pretraining. We show that, despite achieving nearly degenerate, state-of-the-art variational energies, each ansatz converges to a state with qualitatively different spin, charge, and pairing correlations. Upon improving accuracy via symmetry restoration and variance reduction, however, all three converge to the same physical picture: coexisting superconducting and stripe orders. These results demonstrate that variational energy alone is insufficient to identify the ground state in the presence of competing phases, and highlight the importance of tracking how correlation functions evolve as the wave function is systematically improved before drawing physical conclusions.

Lattice compression has emerged as a fundamental tuning parameter for nickelate superconductivity. Pressure acts as a trigger to induce superconductivity in bulk Ruddlesden-Popper nickelates. For infinite-layer nickelate thin films, compressive epitaxial strain and rare-earth ion chemical pressure have been used to substantially enhance the superconducting transition temperature ($T_c$). Efforts to go further have been constrained by the limits of epitaxial stability or the challenges of measuring thin films in high-pressure environments. Here, we overcome this limitation by developing a technique to incorporate freestanding infinite-layer $\mathrm{Nd_{0.85}Sr_{0.15}NiO_2}$ membranes into a diamond anvil cell. Using this platform, we observe a strong increase in $T_c$ up to our highest measurement pressure of $\sim$90 GPa, where a superconducting downturn can be observed near liquid nitrogen temperatures. Strikingly, we find a simple linear enhancement of $T_c$ at a rate of 0.65 K GPa$^{-1}$, with no signs of saturation. This suggests that the pairing strength in infinite-layer nickelates can be raised to a surprisingly high scale, using an approach that can be broadly applied to many two-dimensional materials.

Neural Quantum States (NQS) are now among the most accurate methods for studying strongly correlated many-fermion systems, outperforming existing many-body approaches for large systems. However, NQS calculations remain extremely resource-intensive. Here, we introduce a new Pareto frontier of efficiency and accuracy for NQS in simulating strongly correlated lattice fermions, defined by two complementary backflow-related architectures: the Sparse Convolutional Ansatz for Lattice Electrons (SCALE) (state-of-the-art efficiency) and the Accurate Convolutional ansatz for lattice Electrons (ACE) (state-of-the-art accuracy), benchmarked on the iconic Hubbard and $t-J$ models for large lattices. SCALE uses a tailored convolutional design enabling efficient local updates via low-rank determinant updates, reducing computational scaling from $O(N^4)$ to $O(N^3)$ in backflow methods and yielding a >40$\times$ practical speed-up in tests while maintaining high variational accuracy. As an application, we study the previously inaccessible 1/8-doped pure Hubbard model up to $32 \times 32$, finding no significant energy difference between horizontal and vertical filled stripe states - contrasting with half-filled stripe states when next-nearest-neighbor hoppings are included. ACE employs a deep convolutional stack to maximize expressive power, achieving unprecedented accuracy on large systems. Extensive benchmarks on Hubbard and $t-J$ models show SCALE delivers variational energies competitive with leading methods at a fraction of the cost, while ACE sets a new accuracy benchmark, surpassing recent results with only 1/6 the runtime for $16 \times 4$ systems. These new NQS approaches provide scalable, affordable, and accurate tools for exploring strongly correlated fermionic physics, such as the microscopic mechanism of unconventional superconductivity.

The pair density wave (PDW) is a novel superconducting state with non-zero center-of-mass momentum Cooper pairing in the absence of external magnetic fields. Its realization in microscopic models as the ground state is very rare and extremely challenging, because a genuine PDW state is free of a uniform component or modulations by a pre-existing spin/charge density wave order at the same wavevector. Here, we report the discovery of a genuine primary PDW phase in a two-orbital Hubbard model on the kagome lattice by state-of-art functional renormalization group studies. It emerges out of competing orders over a wide physical parameter range suitable for realistic material realizations. The key ingredients in favor of the PDW order are the strongly sublattice and orbital polarized Bloch states on multiple Fermi pockets. They force the zero-momentum Cooper pairing to involve the same sublattices and be suppressed by onsite Coulomb repulsion, while pairing between different sublattices to be dominated by different Fermi pockets with nonzero total momentum. The degenerate PDW states at three momenta ${\bf M}_{1,2,3}$ on the Brillouin zone boundary exhibit novel intertwined order and can linearly combine into topologically nontrivial chiral PDW states. We propose that the model can be realized in multiorbital kagome materials such as CsCr$_3$Sb$_5$ as well as cold atom systems.

Resonant inelastic X-ray scattering measurements reveal robust magnetic excitations in the perovskite iridates SrIr$_{1-x}$Sn$_x$O$_3$. We analyzed the dispersions of the magnetic excitation with spin-dilution ratio $x$ = 0, 0.03, 0.06, 0.1, and 0.2, crossing from semi-metal to spin-diluted while antiferromagnetically ordered insulators. The extracted effective magnetic exchange interactions decrease continuously upon increasing spin dilution, and their evolution follows a simple spin-dilution scaling law. These results not only verify the strong electron-correlation nature of the metallic parent SrIrO$_3$, but also reveal the entanglement of the charge and spin dynamics in this system.

Electron-doped cuprates require not only electron doping by chemical substitution but also post-growth reduction annealing for realizing superconductivity. However, electron concentration can also be varied by reduction annealing, making it challenging to disentangle the respective influences of electron concentration and oxygen non-stoichiometry. Here, by combining alkali-metal dosing and angle-resolved photoemission spectroscopy, we monitored changes in the electronic structure of an electron-doped cuprate while supplying additional electrons to its surface without modifying oxygen content. Whereas a Fermi surface reconstruction due to long-range antiferromagnetic order was suppressed by alkali-metal deposition, the pseudogap -- which is associated with short-range spin/charge correlations and can be suppressed by efficient reduction annealing -- was found to persist. The results highlight significant contribution of impurity oxygen atoms to pseudogap formation in electron-doped cuprates.

Neural-network variational Monte Carlo (NNVMC) has emerged as a powerful tool for solving quantum many-body problems, yet systematic pathways for improving its accuracy remain largely heuristic. Here, we introduce a physically motivated basis transformation for NNVMC that enhances variational expressivity without increasing the complexity of the neural-network ansatz itself. By formulating the many-body wave function in a Gaussian basis, we introduce a single learnable locality parameter, $α$, that reshapes the target ground state into a more learnable representation. This approach introduces minimal computational overhead and can be readily combined with existing neural-network architectures. Using the three-dimensional homogeneous electron gas as a benchmark, we show that the optimized basis transformation consistently lowers the variational energy for both FermiNet and message-passing neural-network architectures. Notably, for the latter, it enables a more precise determination of the Fermi liquid to Wigner crystal phase transition. More broadly, our results highlight basis transformation as a new route to improving NNVMC in continuous space, showing that accuracy can be enhanced not only by refining the ansatz but also by making the target ground state easier to represent.

There are few exact results for the Hubbard model on bipartite lattices of spatial dimension $d>1$. Nevertheless, the Hubbard model with transfer integral $t$ and onsite repulsion $U$ on bipartite lattices with $N_a$ sites, such as the square, honeycomb, cubic, body-centered cubic, face-centered cubic, and diamond lattices, provides the simplest toy model for describing electronic correlations in many condensed-matter systems and is therefore a quantum problem of considerable physical interest. Seven exact theorems that provide new physical insight into the model are established. Overall, the exact framework based on physical spins and physical $η$-spins for the Hubbard model on bipartite lattices of spatial dimension $d>1$ introduced in this paper offers a robust foundation for future studies of the model, as well as of the condensed-matter materials, such as cuprate superconductors, graphene and graphene-derived systems, and other quantum systems, that it describes.

The Hubbard model and its extensions are canonical theoretical frameworks for understanding correlated electronic states, including those in high-$T_c$ cuprates. Here, we use state-of-the-art thermal tensor network method to map out the temperature-doping phase diagram of the $t$-$t^\prime$ Hubbard model. On the electron-doped side, we find a $d$-wave superconducting (dSC) regime, supporting the scenario of high-$T_c$ superconductivity. In contrast, on the hole-doped side, no robust dSC phase is detected. Instead, a finite-temperature regime dominated by strong pair-density-wave (PDW) fluctuations emerges, which may eventually give way to charge density wave order upon further cooling. The PDW state exhibits inter-arc pairing with net momentum near $(0, π)$, distinct from the zero-momentum pairing in conventional dSC. Furthermore, these fluctuating PDW states occupy the lower portion of the pseudogap regime on the hole-doped side. We provide a comprehensive finite-temperature perspective consistent with previous ground-state studies, shedding new light on pairing instabilities and exotic electronic states in high-$T_c$ superconductors.

Based on the process tensor framework, we extend the time-evolving matrix product operator (TEMPO) method to solve bosonic quantum impurity problems (QIPs) with off-diagonal system-bath coupling. Our method is a most generic extension of TEMPO, which applies for any QIPs as long as the bath is noninteracting and the system is linearly coupled to the bath. It naturally contains all the current developments of TEMPO in more restricted settings. As an application, we study the real-time dynamics of a spin that is coupled to a sub-ohmic bath via the Jaynes-Cummings-type system-bath coupling, and compare it against that of the standard spin-boson model. Our results show that the commonly used secular approximation could easily fail in presence of a structural bath. Our method provides a unified framework to understand different variants of TEMPO and directly suggests a fermionic generalization which has not been explored so far, it could also be straightforwardly used as an impurity solver in the bosonic dynamical mean field theory.

Understanding the transport properties of cuprate superconductors is one of the central challenges in the physics of strongly correlated electrons. The most common approach is to define and solve a low-energy lattice model, but it is still unclear what the minimal model is to capture all relevant mechanisms and provide quantitative predictions. The main uncertainty concerns the choice of the orbital degrees of freedom to be included in the model, as well as the definition of the effective coupling. In this paper, we study the two most commonly considered models, namely the single-orbital Hubbard model and the three-orbital Emery model. We investigate and compare their spectral and transport properties, and find that the two models present a similar, but not the same, physical picture. We identify several strong quantitative differences which might allow one to discriminate between the two models by comparing theory with experiments. We compare our results for several physical quantities with 7 different experiments on 3 different La$_2$CuO$_4$-based cuprates, and in general find excellent agreement. The dc resistivity and the effective mass results suggest that the coupling constant in the effective Hubbard model is larger than expected. We find several more properties that are sensitive to the precise value of the coupling constant, including the critical doping for the Lifshitz transition, and the local spectral weight in the vicinity of the Fermi level; the latter provides a promising way to estimate the effective coupling constant in future photoemission experiments.

The entangled charge and spin dynamics in strongly electron correlated system has been a fruitful playground for exploring new physical phenomena. Here with resonant inelastic X-ray scattering we studied the spin dynamics of SrIrO$_3$, a half-filled paramagnetic semimetal hosting highly itinerant Dirac Fermions due to its topological band structure. Our results show that its magnetic excitations share much similarity to the ordered compounds upon Sn substitution in exchange strength and AFM instability, while the system maintains spin non-ordered. Further, the non-ordered pristine SrIrO$_3$ hosts even longer lifetime magnetic excitations near the AFM zone center comparing to the Sn substituted ordered compounds, contrary to general expectation. These observations indicate an interesting connection between band topology and electron correlation in SrIrO$_3$.

Altermagnetic (AM) order, a recently discovered magnetic state, has attracted intense research interest for its potential applications in spintronic and quantum technologies. Here, we theoretically investigate the AM state in the Lieb lattice, a prototypical two-dimensional lattice, using the Hubbard model. We show that AM order emerges with only moderate electronic correlations. Strikingly, spin-orbit coupling drives the system into a topological phase exhibiting a new quantum spin Hall effect (QSHE) with spin-biased topological edge states in one-dimensional nanoribbons. These edge states possess different localizations and velocities, and hence may produce spin and charge currents, fundamentally distinct from that in conventional topological insulators with spin degeneracy. This novel spin-biased QSHE in the AM Lieb lattice unveils exciting opportunities for both fundamental studies and innovative device concepts, motivating immediate experimental exploration.

"Hidden" phases, generated using ultrafast laser pulses (few hundred femtoseconds), with properties distinct from thermodynamic equilibrium, are appealing for technologies because they can be long-lived, with lifetimes of hours or weeks, and reversible with temperature sweeping or extra pulses. In this regard, La$_{2/3}$Ca$_{1/3}$MnO$_3$ (LCMO) stands out due to its tunability through epitaxial strain, which can drive the bulk ferromagnetic metal (FMM) into an antiferromagnetic insulator (AFI), and its susceptibility to photo-induced transitions. Indeed, AFI LCMO displays a long-lived photo-induced transition into a putative 'hidden' phase whose exact nature and excitations are still largely unknown. Here, we combine ultrafast photo-excitation in the near infrared with in situ transport, x-ray absorption (XAS), and Resonant Inelastic X-ray Scattering (RIXS) to investigate the excitations (polarons, phonons, and orbital) of the photo-excited phase of LCMO and contrast them with the thermodynamic phases achieved through strain and temperature. In the thermodynamic regime, we establish the correlation between polarons and transport, placing them in the 'strong coupling' regime of the Holstein model. Upon photo-excitation of LCMO-AFI, we uncover a long-lived phase characterized by the softening of the polaron excitations, the partial suppression of the Jahn-Teller distortion, and nearly unchanged phonons, showing the emergence of a photo-excited state absent in the equilibrium phase diagram. Finally, by varying temperature, epitaxial strain, and photo-excitation fluence, we construct a polaron phase diagram and identify the key spectroscopic signatures of each phase. Our laser-RIXS approach establishes a versatile platform for exploring photo-induced 'hidden' phases in quantum materials in non-stroboscopic conditions.

Layered two-dimensional electron systems exhibit both optical and acousticlike plasmons around the Brillouin-zone center. In the layered cuprate La$_{2-x}$Sr$_x$CuO$_4$, resonant inelastic x-ray scattering (RIXS) has detected corresponding acousticlike plasmons in a low-energy regime comparable to that of other collective excitations associated with distinct regions of the cuprate phase diagram. This overlap in energy scale raises the question of whether the acousticlike plasmons are significantly influenced by phase-specific electronic phenomena, including the pseudogap, charge and spin order, superconductivity, and strange-metal behavior. Here we show that a single parameter set of the layered $t$-$J$-$V$ model, which incorporates strong correlations and the long-range Coulomb interaction $V$, consistently describes the acousticlike plasmon dispersion across all currently available RIXS data from the underdoped to the heavily overdoped regime. This transferability of a single parameter set exceeds that of earlier theoretical descriptions and supports a picture in which strong correlations persist into the heavily overdoped regime, while the collective plasmon mode exhibits only limited sensitivity to the phase-specific electronic phenomena that distinguish different regions of the phase diagram.

We investigate Mn substitution at the Fe site in PrFe1-xMnxAsO0.7F0.3 (0 to 0.1) using structural, Raman, density functional theory (DFT), transport, and magnetic measurements. X-ray diffraction and Raman analyses confirm preferential Mn incorporation into the FeAs planes, accompanied by lattice expansion and suppression of Fe-related vibrational modes. Electrical transport reveals a systematic decrease of the superconducting transition temperature from 48 K (x = 0) to complete suppression at x = 0.1, together with low-temperature resistivity upturns evolving toward insulating-like behavior. Magnetization and magnetotransport measurements show degradation of superconducting coherence, critical current density, upper critical field, and vortex activation energy with increasing Mn content. The results demonstrate that Mn acts as an efficient magnetic impurity, strongly perturbing the electronic and magnetic environment of the FeAs layers. Comparative analysis indicates relatively enhanced robustness of superconductivity in the Pr-based system, highlighting the role of rare-earth-dependent electronic correlations in impurity effects.

The three-dimensional Hofstadter model exhibits a critical rational flux at which Weyl points emerge in the single-particle spectrum. We study the superconducting regime of the model in the presence of a Hubbard attractive interaction by tuning the magnetic flux across its critical value. We determine the phase diagram in the plane of the coprime pairs parametrizing the magnetic flux. We show that the system exhibits two distinct regimes separated by a critical flux $Φ_c$: for $Φ>Φ_c$, a semimetal-to-superconductor quantum phase transition occurs at a finite interaction strength ($U_c\neq0$), while for $Φ<Φ_c$ superconductivity arises for arbitrarily weak attraction, with a BCS-like exponential scaling of the gap due to the finiteness of the density of states. Close to the transition, we study the scaling behavior and identify the critical exponents. Our results highlight the interplay between magnetic band topology and attractive pairing in three-dimensional Hofstadter systems.

Motivated by recent experiments on ultracold fermionic spin-1/2 $^6$Li atoms in a Lieb lattice at various Hubbard repulsion $U$ and filling fractions $n$ (Lebrat et al., arXiv:2404.17555), we conduct a density matrix renormalization group (DMRG) analysis of the Hubbard model on a two-leg Lieb ladder. At half-filling, we find a ferrimagnetic Mott insulating ground state, consistent with Lieb's theorem. Away from half-filling, a state with finite total spin $\vec{S}^2 \neq 0$ and vanishing charge gap persists down to filling $n_c \approx 2/3$. For lower, incommensurate fillings, the system is described by a Luttinger liquid with one charge and one spin mode. Intriguingly, in a narrow window near $n_c = 2/3$, close to the onset of ferromagnetic order, we identify a superconducting Luther-Emery phase with dominant $s_{xy}$-wave pairing.

The realization of bulk metallic altermagnetism has remained elusive despite the growing number of candidate materials. Here, we present evidence that moderate cobalt substitution ($\sim$15%) drives the correlated narrow-gap semiconductor FeSb$_2$ into a metallic altermagnetic state persisting up to room temperature. The infrared optical conductivity reveals low-energy interband transitions near 0.1 eV that emerge upon doping and grow with Co concentration. Density functional theory calculations show that these transitions originate exclusively from altermagnetic spin ordering, with spin split bands ($\sim$0.2 meV) of non-relativistic origin, together with spin-orbit coupling induced band splitting of the order of $\sim$5 meV near the Fermi level. Co substitution further leads to Fano lineshapes and mode mixing in the infrared-active phonons, reflecting enhanced electron-phonon coupling and local inversion symmetry breaking, while leaving the altermagnetic spin symmetry intact. Our results establish carrier-tuned FeSb$_2$ as a platform for exploring metallic $d$-wave altermagnetism and its coupling to lattice degrees of freedom.

The relationship between the pseudogap and superconductivity remains a central puzzle in the physics of cuprates. Hydrostatic pressure provides a clean tuning parameter free from chemical disorder, yet probing the microscopic energy scales of these phases under compression has remained experimentally challenging. Here, we utilize ultrafast optical spectroscopy to construct the high-pressure phase diagram of the underdoped cuprate Bi$_2$Sr$_2$CaCu$_2$O$_{8+δ}$ up to 37 GPa. Our results reveal a striking dichotomy within the pseudogap state: while the onset temperature $T^*$ rises monotonically with pressure, the energy gap $Δ_{\mathrm{PG}}$ is continuously suppressed. In contrast, the critical temperature $T_{\mathrm{c}}$ and the superconducting gap $Δ_{\mathrm{SC}}$ trace a correlated dome-like trajectory, demonstrating that superconductivity evolves independently from the pseudogap. Furthermore, an abrupt collapse of the gap ratio $2Δ_{\mathrm{SC}}/k_{\mathrm{B}}T_{\mathrm{c}}$ near 8 GPa marks a pressure-driven dimensional crossover, quenching two-dimensional phase fluctuations to stabilize global three-dimensional coherence. Upon reaching 37 GPa, the superconducting condensate is completely quenched into an insulating-like state. By resolving the extended phase evolution, our findings disentangle the pseudogap and superconducting orders, establishing a rigorous experimental basis for the pairing mechanism of high-temperature superconductivity.

The condensation of electron quartets, known as charge-4e superconductivity (SC), represents a novel quantum state of matter beyond the standard paradigm of Cooper pairing. However, concrete microscopic models realizing this phase in two dimensions remain a central challenge. Here, we introduce a non-engineered and sign-problem-free model, unambiguously demonstrating the emergence of a robust and high-temperature charge-4e SC phase using unbiased quantum Monte Carlo simulations. At zero temperature, the phase diagram reveals that charge-4e SC is the primary ground state in the strong-coupling regime. At finite temperature in the absence of charge-2e SC, we identify charge-4e SC through a Berezinskii-Kosterlitz-Thouless transition, marked by a universal jump in the superfluid stiffness consistent with a condensate of charge 4e. Remarkably, the transition temperature Tc increases nearly linearly with interaction strength, providing a robust mechanism for high-Tc quartet superconductivity. Furthermore, spectral analysis reveals a prominent pseudogap above Tc arising from strong phase fluctuations. Our results establish a canonical and numerically exact model system for charge-4e superconductivity, offering crucial guidance for its realization in experimental platforms such as moiré materials and ultracold atomic systems.

The orbital degree of freedom in rare-earth nickelates is typically inactive across the temperature-driven metal-insulator transition, where the system develops two inequivalent Ni sites associated with Ni-O bond disproportionation and breathing-mode distortions of NiO$_6$ octahedra. Here, we show that orbital polarization can be induced by optical excitation with linearly polarized light. Using an interacting multiband tight-binding model combined with real-time simulations of coupled electron-ion-spin dynamics, we find that photoinduced $d$-$d$ transitions reduce the local magnetic moments at Ni sites and effectively suppress Hund's coupling $J$ in the excited state. Importantly, these transitions can be made strongly orbital-selective by tuning the light polarization, leading to an imbalance in $e_g$ orbital occupancies. The resulting nonequilibrium state, characterized by reduced effective $J$ and unequal orbital populations, becomes unstable toward Jahn-Teller (JT) distortions, driving structural relaxation along coherently excited JT modes. Our results demonstrate that polarization-controlled optical excitation provides a pathway to access hidden nonthermal phases with emergent orbital order, enabling coherent control of coupled charge, spin, and lattice degrees of freedom on ultrafast timescales.

Unlike conventional charge-2e superconductors, a charge-4e superconductor exhibits long-range coherence of electron quartets rather than Cooper pairs. Clear zero-temperature realizations of charge-4e superconductivity remain rare. Here, we investigate the zero-temperature phase diagram of the attractive SU(4) Hubbard model with numerically exact, large-scale quantum Monte Carlo (QMC) simulations overcoming major technical hurdles. We identify both charge-2e and charge-4e superconducting phases. Upon increasing interaction, charge-2e correlations are suppressed and eventually vanish, while the charge-4e correlations remain robust and converge with system size, signaling the onset of a quartet-condensed phase. Interestingly, across the charge-2e--charge-4e transition, single electrons remain gapped, while charge-2e correlations exhibit a scaling behavior inconsistent with a conventional Landau description. These features are naturally captured by a fractionalized framework in which the physical charge-2e order parameter is a composite field coupled to an emergent non-Abelian gauge structure. We formulate an Sp(4) gauge-Higgs theory that realizes deconfined quantum pseudocriticality between the Higgs (charge-2e) phase and the confined (charge-4e) phase. The Sp(4) gauge-Higgs theory yields pseudocriticality through a fixed-point collision, and its one-loop collision-point exponents quantitatively track the QMC results. Our results establish charge-4e superconductivity as a bona fide zero-temperature phase, provide a simple model for future studies in a numerically exact framework, and reveal an unconventional route to superconducting criticality.

The hexagonal 4d ruthenates Ba3MRu2O9 host structural dimers and exhibit a delicate balance of competing interactions. Hund's coupling, trigonal crystal-field splitting, and hopping for $a_{1g}$ and $e_g^π$ orbitals all fall within a narrow energy window. This yields a series of possible ground states, ranging from the localized Mott limit with (anti-) ferromagnetic exchange coupling via orbital-selective behavior to the cluster Mott limit with quasimolecular orbitals that are delocalized over the two dimer sites. Using resonant inelastic x-ray scattering, we show that Ba3CeRu2O9 with four holes per dimer resides in the intricate crossover regime between the localized Mott case and the quasimolecular limit. The spin-orbit entangled singlet ground state predominantly shows a Mott-like charge distribution with two holes per Ru site. At the same time, spin and orbital occupation contradict an exchange-based Mott scenario but agree with a cluster Mott approach. A quasimolecular trial wave function describes more than 70% of the ground state. In this crossover regime, small changes of, e.g., the crystal field may strongly affect the character of electronic states.

Auxiliary Field Quantum Monte Carlo (AFQMC) has emerged as a powerful framework for treating strongly correlated electronic systems, offering a favorable balance between computational cost and accuracy. In this paper, we present a novel AFQMC method that uses the isometric tensor hypercontraction (ITHC) technique to diagonalize the two-body Coulomb interaction of molecular electronic Hamiltonians by introducing additional fictitious fermionic modes. Our method shows reduced theoretical complexity and better practical performance for both propagation and local energy evaluation compared to the standard AFQMC method. We demonstrate the efficacy of this approach by computing the ground-state energies of a linear $\ce{H10}$-chain and the benzene molecule. Our results show that the extended-basis AFQMC recovers many-body correlations with a precision comparable to that of high-level wavefunction methods such as Coupled Clusters (CC) or Density Matrix Renormalization Group (DMRG), while offering significantly improved scaling.

Tensor-network simulations of quantum many-body dynamics are fundamentally limited by entanglement build-up, which leads to exponentially growing computational costs. Furthermore, these classical simulation algorithms are inherently sequential as typically a tensor network representation of the quantum state is updated incrementally at each time step. We build on recently introduced randomized quantum algorithms for time evolution (TE-PAI), and adapt them to the classical simulation context with the purpose of enabling massive parallelisation. Our MPS TE-PAI approach achieves exact time evolution on average (unbiased estimator) and proceeds by representing an ensemble of randomized shallow Trotter-variant circuits as tensor networks. As each circuit instance yields a deterministic quantum state (or observable expecation value), the only source of randomness is the sampling of circuit variants; the absence of shot noise therefore yields a reduced estimator variance relative to quantum hardware implementations of TE-PAI. We simulate representative disordered one-dimensional spin-ring Hamiltonians, and numerically observe reductions in the per-sample gate-count by a factor of up to $10^3$ relative to Trotterized MPS evolution, yielding orders of magnitude reduction in the time-to-solution under realistic levels of parallelisation. Finally, we numerically observe that MPS TE-PAI is substantially more robust against severe bond-dimension truncation than product formulas, potentially making it useful for the simulation of strongly correlated systems where truncation is necessary in practice. We also demonstrate that the approach can be used naturally in combination with existing time evolution algorithms, effectively extending their time depth via parallelisation.

We present a comprehensive experimental investigation of the temperature evolution of magnetic states in triangular-lattice delafossite YbCuSe$_2$. Magnetization measurements on high-quality single crystals reveal easy-plane anisotropy. Specific heat, magnetization, and muon spin relaxation ($μ$SR) establish the absence of magnetic order or spin freezing down to 0.03 K ($\leq J_{\mathrm{avg}}/250$), demonstrating a dynamically fluctuating quantum spin liquid (QSL) ground state. Thermodynamic measurements uncover multiple characteristic energy scales at $T_H \approx 4.5$ K, $T_L \approx 1.8$ K, and $T^* \approx 0.7$ K. Below $T^*$, $μ$SR detects a dynamical phase separation in which the majority of the spins are forming a QSL state whereas the remaining spins form a sporadic, disorder-induced state decoupled from the dominant QSL component. Remarkably, the unconventional temperature dependence of the $μ$SR relaxation rate indicates roton-like excitations emerging between $T_H$ and $T_L$, a feature not previously observed in any QSL system, preceding the stabilization of the low-temperature QSL at 0.3 K. These findings identify YbCuSe$_2$ as a unique QSL platform, providing valuable insights for further experimental and theoretical exploration.

The interplay between unconventional magnetism and band topology in two-dimensional materials has emerged as an important theme in condensed matter physics. Here, we present first-principles calculations that reveal the coexistence of d-wave altermagnetism and topological behavior in Janus FeSeX (X = S, Te) monolayers. The chemical asymmetry of the Janus structure breaks both out-of-plane mirror and inversion symmetries, leading to anisotropic exchange interactions and momentum-dependent spin splittings even in the absence of spin-orbit coupling, the defining signature of altermagnetism. Phonon dispersion analyses confirm the dynamical stability of both compounds, while strain-dependent calculations demonstrate that the magnitude of the altermagnetic exchange splitting ($Δ_s$) can be efficiently tuned by biaxial strain. When spin-orbit coupling is included, a finite topological band gap emerges at the Fermi level, accompanied by quantized spin Hall conductivity plateaus and nontrivial topological invariants (spin Chern number = 1, Z2=1). These findings establish FeSeS and FeSeTe as promising two-dimensional platforms for realizing topological altermagnetism and spin--orbit--driven charge--spin conversion, thus opening new avenues for low-dissipation spintronic devices.

The identification and manipulation of charge-neutral fractionalized quasi-particles, in particular chiral edge modes (CEM), is a long-standing quest in physics. Remarkably, the microscopically mediated interaction between light and charge-neutral excitations in Mott-Hubbard insulators can take an identical form to the Raman coupling between light and particles with electric charge. However, since CEMs are Raman-inactive due to conservation of lattice momentum, Raman probes have been deemed unsuitable for their identification. Here, using the Kitaev quantum spin liquid (KSL) as an illustrative example, we demonstrate that the long-range correlated disorder inherent to a closed edge can lead to a Raman circular dichroism (RCD) signal that avoids suppression by linear and angular momentum selection rules, and exhibits a dependence on experimentally tunable length and energy scales that are characteristic of CEM. Having calculated the low-frequency RCD response of generic KSL, we argue that the interaction of the chiral matter fermion with the $\Ztwo$ boundary charge leaves a unique fingerprint of the KSL via the anisotropic Zeeman field dependence.

The emergence of the antiferromagnetic (AFM) Chern insulator (AFCI) phase in the Kane-Mele-Hubbard (KMH) model with a finite sublattice potential is investigated. The AFCI, characterized by AFM correlations coexisting with quantized Hall conductance, has long raised the question of whether it can exist in the KMH model that respects time-reversal symmetry (TRS). Using exact diagonalization, we analyze the excitation gap, anisotropic AFM correlations along the $z$ axis and in the $xy$ plane, and the fidelity susceptibility under twisted boundary conditions, all of which provide consistent evidence for the AFCI phase. In particular, our numerical evaluation on the (spin) Chern number reveals a breakdown of adiabatic continuity in the twist-angle space, indicating an instability toward TRS breaking driven by Hubbard-induced AFM perturbations. A modified computational scheme is further proposed, which yields a robust quantized Chern number $C=1$ within this phase.

Interior-gap superconductivity has long been discussed as an exotic paired state in the presence of Fermi-surface mismatch, but its realization in canonical strongly correlated models has remained elusive. Here we present evidence that the superconducting phase of one-dimensional Kondo-Heisenberg models realizes an interior-gap pair-density-wave (PDW) state generated by strong correlations. Combining infinite density-matrix-renormalization-group (iDMRG) and finite DMRG calculations for $S=1/2$ and $S=3/2$ chains, we show that the PDW correlation is the dominant bulk superconducting correlation in the spin-gapped regime and that the momentum distribution function $n(k)$ exhibits a reconstructed structure characteristic of interior-gap physics. In particular, while the feature in $n(k)$ for the $S=1/2$ chain is only hump-like, the corresponding structure in the $S=3/2$ chain develops into a clear dip, strongly supporting the interpretation in terms of an interior-gap-like dip structure. Unlike conventional interior-gap scenarios based on a mismatch between preexisting Fermi surfaces, the present system starts from a single bare conduction-electron Fermi surface, and the additional low-energy single-particle structure emerges dynamically together with the dominant PDW correlation through the Kondo coupling. Finite DMRG data further demonstrate that boundary effects can substantially modify real-space correlations in this gapless one-dimensional system, making a direct thermodynamic-limit calculation essential for identifying the intrinsic bulk momentum structure and the dominant correlation channel.

Locally constrained gauge theories underpin our understanding of fundamental interactions in particle physics and the emergent behaviour of quantum materials. In strongly correlated systems, they can give rise to quantum spin liquids that lack conventional order and are defined by coherent superpositions of an extensive number of many-body configurations. Realising and probing such exotic states experimentally is an outstanding challenge both in solid-state and synthetic quantum systems, not least due to the difficulty of detecting the fragile coherences between many-body states. Here, we report a large-scale (>3,000 sites) realisation of a two-dimensional U(1) lattice gauge theory with ultracold atoms in a square optical superlattice and demonstrate non-equilibrium preparation of extended regions of U(1) quantum spin liquids. We demonstrate Gauss's law validity in a quench experiment, enabled by a new microscopy technique for detecting doubly occupied sites. We observe characteristic real-space correlations and momentum-space pinch points, hallmarks of the emergent U(1) gauge structure. Using round-trip interferometric protocols, we directly observe large-scale coherence between many-body configurations, providing strong evidence for quantum spin liquid regions extending over ~100 lattice sites. Our results establish non-equilibrium quantum simulation protocols as a powerful route for accessing and probing exotic, highly-entangled states beyond those hosted by the engineered Hamiltonian in thermal equilibrium.

The superconducting transition temperature $T_c$ of the two-dimensional attractive Hubbard model is computed in the vicinity of both ordinary (logarithmic) and higher-order (power-law) Van Hove singularities using determinant quantum Monte Carlo simulations. For interaction strengths $|U| \lesssim W/3$, where $W$ is the electronic bandwidth, $T_c$ is enhanced in the neighborhood of the Van Hove point, albeit more weakly than expected from weak-coupling BCS theory. Enhancing the Van Hove singularity from logarithmic to power-law yields only a minor additional enhancement of $T_c$. For $|U| \gtrsim W/3$, the maximum $T_c$ shifts away from the Van Hove point and instead occurs at a density unrelated to any features in the non-interacting density of states, consistent with a strong-coupling interpretation. We find that the maximal $T_c$ in the model is achieved at intermediate $U$ and at a density away from the Van Hove point.

The $\mathbb{Z}_N$ cluster-state wavefunction, a paradigmatic example of symmetry-protected topological (SPT) order with $\mathbb{Z}_N \times \mathbb{Z}_N$ symmetry, is expressed in various equivalent ways. We identify the projector-based scheme called the $P$-representation as the efficient way to express cluster and dipolar cluster state's wavefunctions. Employing the restricted Boltzmann machine scheme to re-write the interaction matrix in the $P$-representation in terms of neural weight matrices allows us to develop the neural quantum state (NQS) and the matrix product state (MPS) representations of the same state. The NQS and MPS representations differ only in the way the weight matrices are split and grouped together in a matrix product. For both $\mathbb{Z}_N$ cluster and dipolar cluster states, we derive in closed form the weight function $W(s,h)$ that couples physical spins $s$ to hidden variables $h$, generalizing the previous construction for $Z_2$ cluster states to $\mathbb{Z}_N$. For the dipolar cluster state protected by two charge and two dipole symmetries, the procedure we have developed leads to the tensor product state (TPS) representation of the wavefunction where each local tensor carries three virtual indices connecting a given site to two nearest neighbors and one further neighbor. We benchmark the resulting TPS construction against conventional MPS representation using density-matrix renormalization group simulations and argue that the TPS could offer a more efficient representation for some modulated SPT states. As a by-product of the investigation, we generalize the previous $Z_2$ matrix product operator construction of the Kramers-Wannier (KW) operator to $\mathbb{Z}_N$ and interprets it as the dipolar generalization of the discrete Fourier transform on $\mathbb{Z}_N$ variables. The new interpretation naturally explains why the KW map is non-invertible.

The $\mathbb{Z}_N$ cluster-state wavefunction, a paradigmatic example of symmetry-protected topological (SPT) order with $\mathbb{Z}_N \times \mathbb{Z}_N$ symmetry, is expressed in various equivalent ways. We identify the projector-based scheme called the $P$-representation as the efficient way to express cluster and dipolar cluster state's wavefunctions. Employing the restricted Boltzmann machine scheme to re-write the interaction matrix in the $P$-representation in terms of neural weight matrices allows us to develop the neural quantum state (NQS) and the matrix product state (MPS) representations of the same state. The NQS and MPS representations differ only in the way the weight matrices are split and grouped together in a matrix product. For both $\mathbb{Z}_N$ cluster and dipolar cluster states, we derive in closed form the weight function $W(s,h)$ that couples physical spins $s$ to hidden variables $h$, generalizing the previous construction for $Z_2$ cluster states to $\mathbb{Z}_N$. For the dipolar cluster state protected by two charge and two dipole symmetries, the procedure we have developed leads to the tensor product state (TPS) representation of the wavefunction where each local tensor carries three virtual indices connecting a given site to two nearest neighbors and one further neighbor. We benchmark the resulting TPS construction against conventional MPS representation using density-matrix renormalization group simulations and argue that the TPS could offer a more efficient representation for some modulated SPT states. As a by-product of the investigation, we generalize the previous $Z_2$ matrix product operator construction of the Kramers-Wannier (KW) operator to $\mathbb{Z}_N$ and interprets it as the dipolar generalization of the discrete Fourier transform on $\mathbb{Z}_N$ variables. The new interpretation naturally explains why the KW map is non-invertible.

In recent decades, the study of quantum magnets, which feature unconventional behaviour such as exotic quantum phase transitions and quantum spin liquids, and unconventional magnetic states of matter, has made remarkable progress. However, each of the three foundational pillars -- numerical simulations, analytical methods, and, to a lesser extent, materials synthesis and experiments -- often tends to view itself as the primary driver of the field. Even through the need for collaboration among theory, numerics and experiment to understand the complex phases of quantum magnets is well established, yet, in our view there remains a persistent perception from experts in one area that the other two serve merely as supporting tool, primarily useful for validating the dominant ideas of one specialty, and less relevant to shaping the underlying scientific narrative. In this article, we advocate for a different, more integrated approach to overcome the challenges faced by quantum magnetism researchers. We argue that this alternative mindset has already started to advance the understanding of several important quantum magnetic models and their materials realizations.

Collective excitations in fermionic systems play a crucial role in determining their physical properties. An important challenge is to develop efficient theoretical approaches for describing these excitations and their coupling to fermionic degrees of freedom. In this work, we revisit the problem of quantifying the contributions of individual bosonic modes of collective fluctuations to observable properties of correlated fermion systems within the framework of the Fluctuating Local Field (FLF) method. Whereas the auxiliary field in this method was previously considered only classically, we formulate its systematic extension termed Quantum FLF (Q-FLF) that incorporates selected bosonic Matsubara modes, thus tailoring it to description of quantum collective fluctuations. As a testbed, we apply the approach to a half-filled one-dimensional Hubbard chain and compute the Green's function, the total energy, and the antiferromagnetic susceptibility. Our results demonstrate that the proposed scheme enables an efficient and selective characterization of the contributions of individual bosonic modes. In particular, low Matsubara frequencies are found to have a quantitative impact on integrated observables such as total energy and antiferromagnetic susceptibility. At the same time, an accurate description of single-particle properties requires inclusion of higher-frequency bosonic modes.

Geometric frustration gives rise to vast manifolds of degenerate ground states and competing orders in spin and charge systems. Typically, classical ground states are governed by a local ``zero-sum constraint" that relieves frustrated antiferromagnetic interactions or Coulomb repulsion. To date, the paradigm of geometric frustration has yielded a rich landscape of emergent phases, from spin ices and quantum spin liquids to charge glasses. However, an analogous phase rooted in chemical bonding has yet to be firmly demonstrated. Here we report the discovery of bond-density-wave orders induced by geometric frustration in the kagome metal CeRu$_3$Si$_2$ above room temperature. Through synchrotron X-ray diffraction, real-space transmission electron microscopy, and model calculations, we observe two distinct long-period superlattices with harmonic and anharmonic structural modulations. Crucially, interlayer bonds between kagome planes modulate in a sublattice-selective manner to fulfill the zero-sum constraint on the kagome lattice. We demonstrate the potential of kagome metals to host complex bond-ordered states constrained by geometric frustration and establish chemical bonding as a distinct pathway to frustration physics in quantum materials even above room temperature.

Understanding charge transport in strongly correlated systems remains a central challenge in condensed matter physics, particularly in light of the ubiquitous linear-in-$T$ resistivity observed in strange metals across many platforms from bulk cuprates to twisted bilayer graphene. Here, we investigate charge transport in the one-dimensional Hubbard model in the infinite-interaction limit. Focusing on the dilute limit with a fixed number of doped holes, we first construct the exact \emph{and explicit - i.e. beyond Bethe ansatz} energy spectrum and then derive a closed-form analytical expression for the charge Drude weight at arbitrary temperatures. We further analyze the low-temperature scaling and identify a linear-in-$T$ correction to the Drude weight. Upon regularizing the singular Drude contribution to the DC conductivity, we find that this behavior corresponds to an effective linear-in-$T$ resistivity, which may provide analytical insight into the emergence of strange-metal transport in two-dimensional strongly correlated systems.