arXiv Curator

The intricate relationship between high temperature superconductivity and antiferromagnetic order in cuprates, and the fundamental origin of electron pairing remain open questions. By utilizing high-resolution laser-based spatially-resolved angle-resolved photoemission spectroscopy, we investigate the seven-layer $Bi_{2}Sr_{2}Ca_{6}Cu_{7}O_{18+δ}$ (Bi2267) and identify a cuprate system that consists of multiple hole Fermi pockets. The observed Fermi pockets exhibit pronounced momentum-, temperature- and Fermi surface-dependent energy gaps. Crucially, high temperature superconductivity with a critical temperature ($T_{\mathrm{c}}$) of $\sim$75 K emerges in a system with multiple Fermi pockets and the presence of strong antiferromagnetic order and correlations. In particular, substantial electron pairing is observed along the Fermi pocket with an energy gap up to $\sim$42 meV in lightly-doped CuO$_{2}$ planes ($p\sim$0.05). These findings challenge the conventional understanding of the roles of the nodal and antinodal electronic states in driving high-temperature superconductivity. They show that superconductivity and antiferromagnetism can coexist in a cuprate with multiple Fermi pockets, offering further insights into the pairing mechanism in cuprate superconductors.

Recent Scanning Tunneling Microscopy (STM) experiments measuring the superconducting gap features in thin films of superconducting bilayer nickelates La2PrNi2O7 at ambient pressure and compressive strain paved the way to study the Cooper-pairing models and the band-selective identification of the gap features in these systems. Here, using the realistic two-orbital bilayer model and the continuum Green's function formalism, we theoretically analyze orbital and band-selective local density of states as well as the corresponding STM spectra. We find that the multiorbital character and the spatial dependence of the Wannier functions leads to the spectra developing characteristic features depending on the position of the scanning tunneling microscope's tip. This allows for a band-resolved analysis of the superconducting coherence peaks and scattering momenta. We identify a clear path for experimental measurements to not only identify the debated incipiency of the gamma-band, but also identification of the coherence peaks' band origins via distance dependent measurements of the local density of states and its corrections through impurity scattering.

We investigate finite-temperature precursors of Mottness in the repulsive Hubbard model above the spin-ordering temperature, using numerically exact determinant quantum Monte Carlo. We show that the finite-temperature crossover is accompanied by a pronounced suppression of charge fluctuations, despite the presence of a gapless single-particle spectra, demonstrating that Mottness first emerges via two-particle response in an anomalous metallic regime, before appearing in single-particle spectral functions. We further show that a gap formation in the density of states occurs through momentum-resolved redistribution of spectral weight across the Brillouin zone, that begins at the onset of anomalous metallic regime, rather than through gap formation in single-particle spectral functions at individual momenta. Upon doping, the anomalous-metallic regime generates transport and spectroscopic signatures characteristic of doped Mott insulators. This shows that the transport and spectroscopic anomalies of the doped Hubbard model do not require the presence of a fully formed Mott insulator at half-filling, but instead originate from the strong charge-response renormalization and spectral-weight redistribution that develop within the precursor anomalous metallic regime at half-filling.

Fermi arcs observed in underdoped cuprates have sparked debate over whether they represent segments of a large Fermi surface or small Fermi pockets. This ambiguity has long hindered their classification as either the conventional Bardeen-Cooper-Schrieffer (BCS) regime or the strongly coupled Bose-Einstein condensation (BEC) crossover limit. Here, using angle-resolved photoemission spectroscopy and quantum oscillations, we demonstrate the coexistence of a small Fermi pocket and a large superconducting gap in the clean inner CuO2 layers of the four-layer cuprate Ba2Ca3Cu4O8(F,O)2. This coexistence constitutes a hallmark of the BCS-BEC crossover and has remained elusive for decades. Despite the presence of antiferromagnetic (AF) order, the superconducting gap in the small pocket is remarkably large, yielding a gap-to-Fermi-energy ratio (Delta_pocket/e_F ~ 0.6) and a critical-to-Fermi-temperature ratio (Tc/TF ~ 1.3) that reach the theoretical upper bound for two-dimensional superconductivity. Unexpectedly, this BCS-BEC crossover emerges not as the carrier density decreases but as it increases, abruptly within a narrow doping range of less than 1%. These results provide a long-sought microscopic foundation for the d-wave pairing mechanism in doped AF-Mott insulators.

Strongly correlated electron systems often display competing or even intertwined ordering tendencies, hinting to extremely close or degenerate many-electron energies. While degeneracies are directly rooted in the underlying symmetries of the problem under investigation, their multifaceted effects on different response functions and their mutual relations often remain elusive. Here we put this subject on a rigorous basis by investigating the degeneracies of charge, spin, and pairing susceptibilities for the unfrustrated, bipartite Hubbard model. Exploiting its pseudospin symmetry, we analytically derive the mutual relations between these response functions for generic spatial modulations, highly relevant, e.g., for the competition of stripe and superconducting orders. By means of two-particle numerical simulations we demonstrate the occurrence of a simultaneous $d$-wave pairing/$d$-density wave (loop current) instability in the vicinity of the metal-insulator transition, driven by short-ranged spin fluctuations for the two-dimensional case. We show how this degeneracy is gradually lifted by geometrical frustration, which favors superconductivity. Our study provides a general tool for revealing symmetry relations in correlated electron systems and establishes a unifying perspective on the nature of their intermingled charge/loop current, pairing, and spin orders.

We investigate the magnetic ground states and transverse spin excitations of bilayer and trilayer nickelates within a multi-orbital itinerant framework. For the bilayer system, although Hartree-Fock calculations slightly favor a double-stripe order, the calculated excitation spectrum of the single-stripe state, characterized by an anisotropic low-energy cone at $Q_{\text{BL}}$ and isotropic high-energy excitations near $Γ$, exhibits good qualitative agreement with recent RIXS and neutron scattering experiments. We further identify mirror-even optical interlayer modes at $Q_{\text{BL}}$ whose energies match the mirror-odd modes at $Γ$. For the trilayer system, both mirror-odd and mirror-even spin-density-wave states can be stabilized near $Q_{\text{TL}}$, with the mirror-odd state lower in energy in the parameter regime studied. The mirror-odd state hosts an additional nearly gapless mode dominated by the middle layer, while the mirror-even state contains only one acoustic branch together with two gapped optical modes. Comparison with available RIXS data favors the mirror-odd spin-density-wave scenario. Our results show that magnetic excitations provide a sensitive probe of the magnetic order and support a common itinerant origin of magnetism in multilayer nickelates.

Focused helium ion beam (He-FIB) irradiation with 30 keV ions is a key tool for nanoscale patterning and defect engineering in high transition temperature (Tc) cuprate superconducting devices, yet its usable fluence window is constrained by the competing requirements of reliable superconductivity suppression and minimal structural degradation. In this work, we provide a comprehensive dataset on the effects of large-area 30 keV He+ ion exposure on the electric transport and superconducting properties of epitaxial YBa2Cu3O7 (YBCO) thin films. X-ray diffraction shows a fluence-driven loss of crystalline order accompanied by an out-of-plane lattice expansion and an orthorhombic-to-tetragonal transition, culminating at predominant amorphization at the highest fluence of 1 x 10^16 cm^-2. Raman spectra exhibit increasing disorder while lacking signatures of oxygen depletion, indicating that irradiation mainly generates oxygen-related Frenkel defects rather than changing the carrier concentration. Consistently, with increasing fluence, the normal-state resistivity ρ_N(T) at temperature T above Tc increases strongly, while d ρ_N/d T remains nearly unchanged at moderate fluence. The suppression of Tc is accurately described by Abrikosov-Gor'kov pair breaking and reaches complete quenching of superconductivity at 4.5 x 10^15 cm^-2. The anisotropic upper critical fields decrease approximately exponentially with increasing fluence, the vortex activation energy is reduced, and the anisotropy drops, in contrast to oxygen-depleted YBCO. Hall-angle analysis confirms a nearly constant carrier density but a systematic increase in defect scattering and reduced mobility, consistent with a crossover toward the dirty limit at high fluence. These results establish quantitative fluence thresholds and a practical operational window for He-FIB nanopatterning of YBCO quantum circuits.

We investigate the influence of nonlocal interactions on superconducting correlations within the extended Hubbard model. In addition to the on-site Coulomb interaction and nearest-neighbor hopping, we include next-nearest-neighbor hopping together with several physically relevant nonlocal terms, namely nearest-neighbor Coulomb interaction, correlated hopping, exchange interaction, and pair hopping. Using Lanczos exact diagonalization on a $4\times4$ cluster, supported by projector quantum Monte Carlo simulations for selected parameter regimes, we analyze pairing correlations in both the s- and d-wave channels. We demonstrate that nonlocal interactions exert a highly nontrivial and symmetry-dependent influence on superconducting correlations. While the on-site repulsion in cooperation with next-nearest-neighbor hopping enhances d-wave pairing tendencies, correlated hopping and nearest-neighbor Coulomb interaction strongly promote s-wave correlations, whereas exchange and pair-hopping interactions can efficiently suppress superconductivity beyond relatively small critical strengths. When all nonlocal interactions are considered simultaneously, the resulting phase diagrams reveal a complex interplay and competition between different pairing symmetries. Our results highlight the crucial role of extended interactions in shaping the pairing landscape of strongly correlated systems and demonstrate that a comprehensive treatment beyond the minimal Hubbard model is essential for a realistic description of unconventional superconductivity.

We investigate the thermodynamic properties of elemental cerium by assessing the crucial role of phonon free energy within the framework of dynamical mean-field theory (DMFT). While conventional density functional theory (DFT) often fails to capture the intricate energy landscape of $f$-electron materials, our approach integrates many-body electronic correlations with lattice dynamics to achieve a more rigorous description of the equation of state. We calculate the total energy as a function of the lattice constant at both the DFT and DFT+DMFT levels, subsequently incorporating the vibrational free energy derived from the phonon density of states. Our findings reveal that electronic renormalisation of the force constants significantly alters the phonon spectra, particularly in the strongly correlated $γ$-phase. By applying these phonon corrections to the energy profiles, we observe a substantial refinement in the predicted equilibrium volumes. Using principal-component-based machine learning, we interpolate phonon dispersions continuously from a finite set of first-principles calculations and compare them to experiment, finding significantly closer agreement compared to conventional DFT and DFT+U calculations that neglect dynamical many-body correlations. This study underlines the necessity of accounting for both electronic and vibrational entropy when evaluating the phase stability and structural transitions of lanthanide systems under varying pressures and temperatures.

Motivated by recent progress in the 1313-La$_3$Ni$_2$O$_7$ nickelate thin films (Nie et al., Nature {\bf 652}, 628 (2026)), we systematically investigate the effects of both compressive and tensile strain in this system. Self-doping effects between the single-layer (SL) and trilayer (TL) blocks are observed in our studies in both cases, but are most pronounced under tensile strain. We find that superconductivity is unlikely to emerge in the 1313 film on the LSAO substrate even with further hole doping. Remarkably, within the random-phase approximation, under {\it tensile} strain, a robust $s^{\pm}$-wave pairing state emerges in the TL subsystem with sign changes between the small electron-like $σ$ pocket at the $Γ$ point and the small hole-like $γ$ pocket at the M point. These pockets are connected by a wave vector close to $(π,π)$. Our calculations also suggest that superconductivity in 1313-LNO requires an optimally sized $γ$ pocket, because an oversized $γ$ pocket suppresses pairing formation. Overall, our results predict strain-dependent electronic reconstruction in 1313-La$_3$Ni$_2$O$_7$ and provide guiding principles for engineering superconductivity under ambient pressure conditions.

The discovery of superconductivity in Ruddlesden-Popper bilayer nickelates under both high pressure and ambient conditions has opened a new paradigm for exploring unconventional superconductivity. This review provides a brief survey of theoretical progress on bilayer nickelate superconductors. Drawing from the key experimental observations, we summarize essential physical ingredients including the hybridized Ni-3$d_{x^2-y^2}$ and 3$d_{z^2}$ electronic structure, orbital-dependent electronic correlation, Hund's coupling, and strong interlayer magnetic coupling. The fundamental theoretical models including the bilayer two-orbital Hubbard model and its minimal $t$-$J$ variants are introduced. Starting from the atomic-limit interlayer valence bond picture of the half-filled $d_{z^2}$ orbital, we elaborate on strong correlation interlayer pairing mechanisms based on different limiting considerations. Specific emphasis is placed on the hybridization mechanism, where the $d_{z^2}$ local singlet pairs provide the pairing energy and their hybridization with itinerant $d_{x^2-y^2}$ promotes superconducting phase coherence. We further analyze the pairing symmetry, the dependence of $T_c$ on various internal and external parameters, the nontrivial normal state properties including the Fermi liquid, non-Fermi liquid, weakly insulating and pseudogap behaviors. Effects of pressure tuning, oxygen content, and Kondo scattering induced by oxygen vacancies are also discussed. Finally, weak correlation theories based on spin fluctuations associated with Fermi surface nesting are briefly covered.

The discovery of superconductivity at $15\,\mathrm{K}$ in the infinite-layer nickelate $(\mathrm{Nd},\mathrm{Sr})\mathrm{NiO}_2$, followed by superconductivity at $80\,\mathrm{K}$ in the Ruddlesden--Popper phase $\mathrm{La}_3\mathrm{Ni}_2\mathrm{O}_7$, has ushered in a new era of nickelate research. Despite this progress, large discrepancies between reports exist. Here, we investigate the complete series of bulk-stable $\mathrm{La}_{n+1}\mathrm{Ni}_n\mathrm{O}_{3n+1}$ compounds using a comprehensive set of experimental techniques, including PXRD, single-crystal XRD, electron microscopy, heat capacity, differential scanning calorimetry, magnetic susceptibility, and transport measurements, over a broad temperature range from $2$ to $1000\,\mathrm{K}$. By studying high-quality single crystals, we identify a previously underappreciated high-temperature phase transition in Ruddlesden--Popper nickelates $\mathrm{La}_{n+1}\mathrm{Ni}_n\mathrm{O}_{3n+1}$ distinct from the one going to a tetragonal phase.

The low-pressure phase diagram of La$_3$Ni$_2$O$_7$ provides an important reference for understanding its pressure-induced high-temperature superconductivity. While the spin-density-wave transition at $T_{\text{SDW}}\approx150$ K is increasingly well established, the origin of the second density-wave transition at $T_{\text{DW}}\approx130$ K has remained unresolved. Here, we perform unrestricted Hartree-Fock calculations to investigate the potential origin of the second transition. {Within the orthorhombic phase, the degeneracy between possible ordering wavevectors at $\boldsymbol{Q}_{Y}=(0,π)$ and at $\boldsymbol{Q}_{X}=(π,0)$ is lifted and the electronic system} develops a double-stripe spin-density wave with ordering vector $\boldsymbol{Q}_{Y}=(0,π)$. We identify that the pure double stripe spin state is unstable in La$_3$Ni$_2$O$_7$ towards a commensurate charge-density wave instability, which favors a spin-modulated double stripe order with intertwined charge and spin instabilities and establish the hierarchy of ordered states in La$_3$Ni$_2$O$_7$, providing an important link between its ambient-pressure and superconducting high-pressure phases. We further discuss our results in the context of available experimental literature and propose further experimental tests to elucidate the origin of the SDW/DW states in this system.

We leverage $s$-orbital non-resonant inelastic X-ray scattering ($s$-NIXS) to perform orbital imaging on three bulk rare-earth nickelates spanning a range of formal nickel valence (3$d$ electron filling) from Ni$^{3+}$ (3$d^7$) to Ni$^{1+}$ (3$d^9$). Our results directly reveal the ground states of these compounds all with minimal theoretical input. In particular, we demonstrate the low-spin orbital configuration of trivalent LaNiO$_3$, the $d_{x^2-y^2}$ configuration of monovalent LaNiO$_2$, and resolve the effective $e_g$ crystal field splitting in the distorted octahedral environment of divalent La$_2$NiO$_4$. This work illustrates the potential of $s$-NIXS to study the ground state and excited states of strongly correlated materials without needing complex theoretical analysis of spectroscopic data.

Motivated by the discovery of superconductivity in Sr-doped infinite-layer nickelate films on SrTiO$_3$(001), we explore the broader landscape of $AB$O$_2$ oxides through comprehensive high-throughput first-principles simulations. Specifically, delafossites and their ordered rock-salt (111) variants stand out as intriguing layered oxides that share the infinite-layer $AB$O$_2$ stoichiometry and simultaneously retain a perovskite-like octahedral motif. This positions them as a unique structural bridge between these two phases and as promising candidates for novel correlated electronic states. We compile a phase diagram that compares the relative stability of these four distinct oxides across the periodic table. Surprisingly, we find that the delafossite structure rivals the infinite-layer phase in thermodynamic stability for the nickelates, and even more for the recently suggested palladate and platinate analogs. Comparison of the respective electronic structures reveals that the delafossite compounds, which we find to be characterized by reversed cation order, exhibit a strongly $d_{z^2}$-dominated Fermi surface, in stark contrast to the $d_{x^2-y^2}$ character observed in the infinite-layer phases. Among all candidates, the La-Ni combination stands out as a thermodynamic optimum for stabilizing the infinite-layer motif. Furthermore, we show that hole doping via Ca, Sr, and Ba systematically enhances the stability of the infinite-layer phase in all three transition-metal families. These results reveal fundamental challenges in realizing bulk substrate-free infinite-layer oxides, and simultaneously offer guidance for future experimental synthesis efforts targeting novel superconducting compounds.

We investigate finite-temperature precursors of Mottness in the repulsive Hubbard model above the spin-ordering temperature, using numerically exact determinant quantum Monte Carlo. We show that the finite-temperature crossover is accompanied by a pronounced suppression of charge fluctuations, despite the presence of a gapless single-particle spectra, demonstrating that Mottness first emerges via two-particle response in an anomalous metallic regime, before appearing in single-particle spectral functions. We further show that a gap formation in the density of states occurs through momentum-resolved redistribution of spectral weight across the Brillouin zone, that begins at the onset of anomalous metallic regime, rather than through gap formation in single-particle spectral functions at individual momenta. Upon doping, the anomalous-metallic regime generates transport and spectroscopic signatures characteristic of doped Mott insulators. This shows that the transport and spectroscopic anomalies of the doped Hubbard model do not require the presence of a fully formed Mott insulator at half-filling, but instead originate from the strong charge-response renormalization and spectral-weight redistribution that develop within the precursor anomalous metallic regime at half-filling.

We investigate the interplay between spin-density-wave (SDW) order and superconductivity in the bilayer nickelate La$_3$Ni$_2$O$_7$ using the functional renormalization group~(fRG) applied to multiorbital weak-coupling models of both the ambient- and high-pressure crystal structures. As Hund's coupling increases, the leading instability evolves from superconductivity to an SDW state with ordering vector $\mathbf{Q}_1 \approx (π/2,π/2)$ (equivalently $\mathbf{Q}_Y \approx (0,π)$ in the orthorhombic $Amam$ unit cell), in agreement with experimental observations. Surprisingly, the ambient- and high-pressure structures exhibit nearly identical non-interacting susceptibilities and leading fRG instabilities, indicating that the emergence of superconductivity under pressure cannot be explained solely by changes in the low-energy electronic structure. Instead, our results identify the suppression of orthorhombicity as a key ingredient for superconductivity. As the system approaches the tetragonal limit, symmetry-related SDW fluctuations become nearly degenerate, frustrating long-range magnetic order while enhancing pairing interactions. These findings highlight lattice symmetry as a central tuning parameter of the competing ordered states in bilayer nickelates and suggest that reducing orthorhombicity through uniaxial strain could stabilize bulk superconductivity already at ambient pressure.

Unraveling the mechanism that binds electrons into Cooper pairs in cuprate high-temperature superconductors remains one of the most fundamental challenges in condensed-matter physics. While both magnetic interactions and lattice vibrations are known to govern key electronic properties, their possible cooperation has never been directly observed. We investigate HgBa$_2$Ca$_2$Cu$_3$O$_{8+δ}$ (Hg1223) - the cuprate with the highest $T_{\mathrm{c}}$ at ambient pressure - as a magnifying glass to probe the possible entwining of the excitations at the core of the pairing. Using resonant inelastic X-ray scattering, we find that the charge response is dominated by dynamic charge density fluctuations (CDF) extending up to several hundred meV, where magnetic excitations reside. At the same momentum where CDF are most intense, the paramagnon energy exhibits a pronounced softening, revealing a strong interplay among charge, lattice, and spin excitations. Our results point to a cooperative mechanism in which dynamic charge fluctuations mediate the coupling between lattice, charge and spin degrees of freedom-shedding new light on the fundamental origin of high-$T_{\mathrm{c}}$ superconductivity.

Ruddlesden--Popper nickelates have emerged as a promising platform for high-temperature superconductivity, yet the role of lattice degrees of freedom in their correlated normal state remains largely unexplored. Here, we report the observation of a finite phonon thermal Hall effect in the trilayer nickelate La$_4$Ni$_3$O$_{10}$ at ambient pressure. Remarkably, the thermal Hall response is strongly enhanced below the density-wave transition at $T^*\approx140$ K, exhibiting two distinct plateaus in the thermal Hall resistivity. The characteristic energy scale extracted from the thermal Hall response ($\sim4.1$ meV) closely matches the magnon--phonon crossing span energy ($\sim3.2$ meV), pointing to magnon--phonon hybridization as the primary mechanism enhancing the thermal Hall effect. These results provide new insight into the interplay between lattice and spin excitations in nickelates, with implications for understanding both their superconductivity and the multiple possible origins of insulating thermal Hall signals.

We present a comparative Ni core level photoemission study of La$_3$Ni$_2$O$_7$, Nd$_3$Ni$_2$O$_7$, and LaNiO$_3$ using both the Ni $2p$ and the Ni $1s$. We address the challenges in analyzing the widely investigated Ni $2p$ spectra arising from the substantial overlap in energy of the Ni $2p$ with the La $3d$. We show that on the other hand the deep Ni $1s$ core level does provide a clean view on the intrinsic electronic excitations and we highlight its potential to resolve detailed differences in the electronic structure within the strongly correlated Ruddlesden-Popper series La$_{n+1}$Ni$_n$O$_{3n+1}$.

The topological heavy-fermion model provides a minimal framework for describing the coexistence of localized moments and itinerant Dirac electrons in magic-angle twisted bilayer graphene (MATBG). Several analytical and numerical methods have been applied to this model; however, whether they provide a realistic description of MATBG remains incompletely understood. In this work, we develop an Hubbard operator approach that incorporates non-local correlations beyond the single-site limit. We benchmark the approximate calculations against numerically exact determinant quantum Monte Carlo simulations of a lattice-regularized model. We show that commonly used local approximations, such as Hubbard-I, fail to capture the coupling between localized and itinerant degrees of freedom, leading to incorrect spectral properties in the local-moment regime. In contrast, the Hubbard operator method provides a controlled description of both correlation functions and spectral features over a regime of parameters, in good agreement with exact numerical methods.

We study a two-orbital Hubbard-Kanamori model, which has been originally proposed for iron-based superconductors, using variational Monte Carlo. We span the nonmagnetic sector at both hole-doping and electron-doping, with respect to the half-filled case $n=2$. We report the presence of a superconductive region with a $s^{\pm}$ symmetry only when the half-filled system is in a Mott state, while orbital selectivity is absent. These results are qualitatively different from what was reported in the three-orbital Hubbard-Kanamori model, where a more extended superconductive region was observed with a concomitant development of orbital selectivity, and they are to some extent more reminiscent of the single-band Hubbard model.

The pairing mechanism of the bilayer nickelates La$_3$Ni$_2$O$_7$ remains a hotly-debated open question. Existing strong-coupling theories are divided into class favoring intralayer d-wave pairing and that favoring interlayer s-wave pairing, with the latter further divided into $d_{z^2}$ orbital dominated mechanism driven by orbital hybridization and $d_{x^2-y^2}$ orbital dominated mechanism driven by Hund's rule. Recent angle-resolved-photoemission-spectrum (ARPES) and scanning-tunneling-microscope (STM) combinedly reveal a nodeless full pairing gap with low anisotropy, supporting the s-wave pairing. Here we propose that the pairing gap along the Brillouin zone (BZ) diagonal can serve as a useful probe of pairing mechanism. Symmetry analysis suggests that orbital hybridization vanishes along the BZ diagonal, rendering that the pairing gaps on the $γ$- and $α/β$- pockets reflect the $d_{z^2}$- and $d_{x^2-y^2}$- orbital pairing strength respectively. Under the $d_{z^2}$ orbital dominated pairing mechanism driven by orbital hybridization, gap nodes are inevitable on the $α$- and $β$- pockets along the BZ diagonal, which conflicts with the full gap revealed by ARPES and the U-shaped dI/dV curve observed by STM. The Hund's rule driven pairing mechanism instead leads to a full pairing gap, which well fits the ARPES and STM results. Furthermore, through a random-phase-approximation based calculation, we show that the weak-coupling theory, which tends to yield a $d_{z^2}$-orbital dominated pairing, also leads to nodes or near-nodes on the $α$- and $β$- pockets along the BZ diagonal, conflicting with experiments. This analysis clarifies the dominant role of $d_{x^2-y^2}$ orbital in the pairing and establishes the Hund's rule driven pairing mechanism as the most relevant one in La$_3$Ni$_2$O$_7$.

Distilling the intrinsic physics of the superconducting CuO2 plane from the complexities of charge-reservoir layers is a defining challenge in high-temperature superconductivity. While superconducting electron-doped infinite-layer cuprates have been synthesized, controllable and uniform hole doping has long remained elusive despite exploratory attempts, limiting spectroscopic insights. Here, we realize bipolar doping across infinite-layer (Sr,Eu)CuO2 and (Ca,Li)CuO2+δ single-crystalline thin films, mapping the electronic phase diagram. Both electron- and hole-doped films show pronounced electrical resistance anisotropy, indicating the quasi-two-dimensional nature of the CuO2 planes. Angle-resolved photoemission spectroscopy across electron- and hole-doped regimes reveals persistent antiferromagnetic band folding coexisting with superconductivity. Remarkably, at a hole doping ~0.07 determined by Luttinger volume, the antiferromagnetic folding emerges from Fermi arcs within the film's single Fermi surface, with the onset superconducting transition temperature exceeding 60 K. These findings redefine the interplay between magnetic order and superconductivity and establish a definitive platform to investigate the intrinsic mechanism of high-temperature superconducting cuprates.

A breakdown of calculations with exact impurity solvers in the dynamical mean field theory in multiband systems easily occurs due to the expensive numerical cost. To overcome this practical difficulty, three of the present authors developed an inexpensive and reliable impurity solver by combining the iterative perturbation theory (IPT) and parquet equation, and named it IPT+parquet [R. Mizuno, et al., Phys. Rev. B 104, 035160 (2021).]. In this study, we validate IPT+parquet focusing on the orbital fluctuation by comparing the numerically exact impurity solvers. We confirm that IPT+parquet can capture competition between orbital fluctuation channels, which the conventional IPT cannot capture.

Higher-order correlations are essential for understanding exotic phases and uncovering universal aspects of quantum dynamics. While inelastic neutron scattering provides well-established access to two-particle correlations, measuring correlations in solids involving more than two particles remains a major challenge. Focusing on magnetic excitations in the Heisenberg spin-$\frac{1}{2}$ chain, we demonstrate that resonant inelastic X-ray scattering (RIXS) can selectively probe multi-spin dynamical correlations by exciting distinct intermediate states through energy detuning. Through theoretical modeling, we isolate the two- and multi- spin responses, and establish that: (i) The resonant energies of the two- and multi- spin dynamical correlations are separated by $\approx \frac{3}{4}J$; implying that the multi-spin part of the cross-section comes from intermediate states that contain spin flips; (ii) The spectral weight of the two-spin channel exhibits a Lorentzian resonant energy profile consistent with a single dominant electronic configuration, whereas the multi-spin channel response shows a Gaussian-like profile, implying contributions from multiple intermediate spin configurations. These characteristics are reproduced by exact diagonalization of the $t-J$ Hamiltonian, which further reveals that the widths of Lorentzian and Gaussian resonant energy profiles depend primarily on core-hole lifetime ($Γ/2$) and $\frac{2J}Γ$, respectively. Controlled access to multi-spin dynamics, using RIXS energy detuning as a knob, can open new pathways to explore many-body dynamics in quantum materials.

The doping and temperature dependencies of the electronic structure and the dynamic spin susceptibility of the Hubbard model on a square lattice are studied within cluster perturbation theory (CPT). The most general features of both electronic and spin spectra qualitatively agree with the experimental data on cuprates. Generalized mean-field calculations of the electronic structure using static short-range magnetic correlations from CPT are also implemented and compared with the results of CPT, allowing us to discuss the role of short-range antiferromagnetism in the formation of the pseudogap.

Triple-layer cuprate superconductor $\mathrm{HgBa_2Ca_2Cu_3O_8}$ (Hg1223) keeps the record of the highest superconducting (SC) critical temperature $T_{c}\sim 134$K among all the existing materials at ambient pressure. $T_{c}$ further increases under pressure up to $T_{c}\sim 160$K. However, its microscopic mechanism remains to be elucidated. We solve {\it ab initio} Hamiltonians for Hg1223 using a variational solver supplemented by a neural network. The pressure dependence of the $d$-wave SC order parameter and estimated $T_{c}$ show a dome-like structure in essential agreement with the experimental indications. The origin of the strong SC amplitude at ambient pressure is identified as strong local Coulomb repulsion $U$ attributed to poor screening. Further increase in $T_{c}$ under pressure is understood from interplay of three elements, namely increased electron hopping $t$, decreased $U$ and more importantly, strongly reduced offsite Coulomb repulsion $V$ with increasing pressure. Pairing mechanism is identified as the emergent local attraction counterintuitively generated from the originally strong local repulsion $U$. The emergent attraction is interpreted from ``attraction from reduced repulsion'', originating from the release of the fluctuating doubly-occupied sites characterized from the ``false vacuum'' in the Mott insulator to the double-occupation-free $d$-wave SC states upon carrier doping. This instantaneous attraction is in contrast with the conventional BCS SC mediated by bosonic glues. The local attraction is consistent with the electron fractionalization supported in experimental analyses. The coexistence of the SC and antiferromagnetic order is also demonstrated as a characteristic feature of the multi-layer system. The microscopic understanding of Hg1223 offers a new route explicitly using this emergent attraction to design and optimize SC materials.

The three-dimensional repulsive Hubbard model is investigated using the strong-coupling diagram technique. For half-filling, the boundary between paramagnetic and antiferromagnetic states is determined for the range of the Hubbard repulsion $4t\leq U\leq12t$, where $t$ is the hopping integral between neighboring sites. Along this boundary, the density of states is calculated, and it demonstrates the Mott transition at $U\approx9t$. For $U\geq6t$ and half-filling, the density of states has the shape inherent in the strong electron repulsion, while for $U=4t$, its shape points to weak coupling. The dependence of the Néel temperature $T_{\rm N}$ on the electron concentration $\bar{n}$ is investigated for the cases $U=4t$ and $12t$. In the former case, $T_{\rm N}$ decreases monotonously with $\bar{n}$, while in the latter case, there is a plateau in the dependence near $\bar{n}=0.87$. The plateau is connected to a reconstruction of the density of states caused by an effective weakening of electron coupling due to electron depopulation. For $U=12t$ and half-filling, the magnetic critical exponent $γ\approx1.4$, which is close to the value in the Heisenberg model. Some features resembling the first-order phase transition, revealing themselves in a finite crystal, are discussed.

We consider the spectral properties of the two-dimensional Hubbard model, describing the electronic properties of high-$T_c$ compounds, within the SU(2) gauge theory, which assumes the separation of electronic degrees of freedom into those of spinon and chargon subsystems. We use the dynamic mean-field theory (DMFT) approach to describe magnetic long-range order in the chargon subsystem while also treating spinon fluctuations on top of this state. We show that DMFT supplemented by long-wavelength magnetic fluctuations is essential for describing the asymmetry in the damping between the inner and outer regions of the hole pockets and the resulting formation of Fermi arcs in the underdoped regime, especially at low hole doping. The underlying hole pockets in the chargon subsystem can be associated with those observed in quantum oscillation measurements.

Strontium ruthenate (Sr$_2$RuO$_4$, SRO) has been the subject of intense research as a model quasi-two-dimensional metal with strong electronic correlations and potential exotic multi-component superconductor. Yet the nature of the superconducting state and its emergence remain debated, despite highly detailed knowledge of the normal-state electronic properties. Here we use diffuse neutron and x-ray scattering to uncover self-organized structural inhomogeneity on the picosecond timescale in SRO. We show that these structural correlations do not originate from extrinsic disorder but rather involve correlated displacements of oxygen atoms in the quintessential RuO$_2$ planes. Moreover, the observed displacement pattern is consistent with distortions due to orbital order that we obtain in first-principles calculations, which suggests that orbital effects could play a pivotal role in the physics of SRO. The appearance of such dynamic inhomogeneity may be relevant for a wide range of prominent oxides with similar lamellar structures, such as the cuprates and nickelates.

``Infinite layer'' nickelates (ILNs) ${\cal R}$NiO$_2$ (${\cal R}$=rare earth elements), having empty apical O sites, become superconducting upon hole doping. They display a secondary electron Fermi surface (FS), giving hole doping, arising not from atomic orbitals but from a band based on interstitial density. Newly reported La$_3$Ni$_2$O$_5$F, formally Ni$^{1+}$, provides an unexpected example of ILN with essentially ideal two dimensional character. A partially occupied single band $E^*$, based on interstitial density, has distinct properties, as its strongly anisotropic shape extends over the three ``apical'' layers and leads to a cylindrical electron FS giving self-doping. This interstitial density is associated with a {\it network of valence bands}, including a Ni $d_{xz},d_{yz}$ pair that partners with $E^*$ to provide an incipient non-analytic Dirac point, leading to an unusual type of interstitial density--$d$ band coupling. The $E^*$ electron band and the conventional Ni $dpσ$ band will display a dichotomy of hole and electron quasiparticle behavior in normal state transport and far-IR properties, and likely resulting in unconventional superconducting state properties even for nickelates.

Superconducting spin qubits, realized as Yu-Shiba-Rusinov spin-doublet states in quantum-dot-superconductor systems, represent a cornerstone of current research in quantum technologies. We analyze these ground states of quantum impurities in superconducting pseudogap systems, namely systems with a pseudogap tunneling density of states $ρ(ε) \sim |ε|^r$ for energies $|ε|\gg Δ$ ($Δ$ being a $s$-wave pairing potential). For $r=1$, these hosts are realized as Dirac materials (graphene or 3D topological insulator surfaces) in proximity to conventional superconductors, or as $d+i s$ superconductors. Using effective field theory and numerical renormalization group, we map the phase diagram against the pseudogap exponent $r > 0$ and particle-hole symmetry-breaking perturbations. At particle-hole symmetry, increasing $r$ also increases the critical value, $J_c$, of the Kondo coupling that triggers the transition from spin doublet to singlet. Unlike the gapless pseudogap Kondo systems, numerical and analytical evidence suggest that Andreev reflection stabilizes a singlet ground state at $J$ for all $r > 0$. Breaking particle-hole symmetry -- by potential scattering or chemical potential -- eventually restores the transition at lower $J_c$. Our results indicate that coupling to superconducting hosts with large pseudogap exponents enhances the stability of spin qubits at large Kondo coupling.

We present realistic calculations within the density functional plus cluster dynamical mean-field formalism indicating that the physics of the the bilayer Ruddlesden-Popper nickelates is to a significant extent controlled by interlayer "dynamical singlets'' which are formed from the $3z^{2}-r^{2}$ orbitals singly occupied by electrons and are hybridized with itinerant planar $x^{2}-y^{2}$ orbitals. The hybridization is found to respond differently to hydrostatic pressure and to epitaxial strain, capturing the experimentally observed dichotomy between bulk single crystals and epitaxial thin films and reproducing several experimental results including angle-resolved photoemission and transport measurements.

We study the two-dimensional Kondo-Heisenberg model on a square lattice, with the conduction electrons away from half-filling, using neural network quantum states. Mapping the ground-state phase diagram as a function of the Kondo and Heisenberg couplings, we identify (i) at weak Kondo coupling, antiferromagnetic Néel order with a Fermi surface whose enclosed area counts only the conduction electrons and is insensitive to the Néel order, and (ii) at strong coupling, a heavy Fermi liquid with a Fermi surface whose enclosed area counts both the conduction electrons and the spins. In the crossover between these regimes, we find $d_{x^2-y^2}$ superconductivity, evidenced by off-diagonal long-range order in the pair-pair correlations and a pairing-amplitude dome that coexists with the underlying magnetic phase. Our results establish Fermi volume change and unconventional superconductivity as intrinsic features of the two-dimensional Kondo-Heisenberg model.

We present a new method for the efficient imaginary time evolution of quantum many-body wavefunctions represented by matrix product states (MPS). We first show that logarithmic time grids are sufficient to resolve long imaginary time dynamics, yielding an exponential reduction in the number of time steps compared with standard approaches. We then show that A-stable implicit time-stepping methods for ordinary differential equations allow stable propagation for any time step size. The resulting scheme requires only matrix-vector products and linear solves, standard operations in the MPS toolbox. We validate our approach with two examples: a Heisenberg spin chain, which we use to demonstrate a speedup of several orders of magnitude over the standard time-dependent variational principle method with uniform time steps, and a single-site Anderson impurity model with a metallic bath, for which propagation to large imaginary times allows one to observe the exponential dependence of the Kondo temperature on the interaction strength.

Neural-network quantum states (NQS) provide a flexible nonlinear representation of quantum many-body wavefunctions, but their efficiency depends sensitively on whether the architecture reflects the sign structure and constrained Hilbert space of the target state. In this work, we propose a solvable-point-guided strategy: design the architecture at an exactly solvable point where the correct local rules can be read off, then refine to the non-exact regime by enlarging only the kernel size and hidden dimension. The strategy is built from four physics-motivated designs: a stride-matched local-rule convolution, geometric pooling, a sign-resolving $\tanh(x^{2k+1})$ activation, and explicit spin-hole sector separation. We test this approach on quasi-one-dimensional valence-bond-solid (VBS) states and their doped soliton variants (sVBS), the exact ground states of a $t$-$J$-like model with a single mobile hole. In finite-size benchmarks, this architecture reaches high fidelity for the exact sVBS state with substantially fewer parameters than generic fully connected, convolutional, and transformer baselines tested under the same setup. For the spin sector, the learned local rule transfers from small to larger systems without retraining. Away from the solvable point, increasing kernel size and hidden dimension systematically improves accuracy, and the model shows approximately $L^2$ parameter scaling in the gapless regime for system size $L$, compared with approximately $L^4$ for matrix-product states in the same regime. Our work establishes a recipe for compact NQS in sign-structured, constrained Hilbert spaces and paves the pathway to physics-informed architectures for the broader $t$-$J$ and Hubbard families.

The emergence of heavy quasiparticles in a Kondo lattice is usually viewed as the formation of a low-energy hybridization gap. Whether this gap represents a rigid electronic structure or a fragile many-body state that can be dynamically reconfigured remains a central question for heavy-fermion systems near magnetic order, quantum criticality, and unconventional superconductivity. Here we use femtosecond pump-probe reflectivity to interrogate this problem in the weakly hybridized Kondo-lattice compound CePt$_2$In$_7$. At low fluence, a slow quasiparticle relaxation channel emerges below $T^* \sim$ 40 K and follows a Rothwarf-Taylor bottleneck response with a low-energy recombination scale 2$Δ\approx$ 7.4 meV. Coherent optical phonons, independently identified by Raman spectroscopy, act as an internal lattice thermometer and rule out large quasi-equilibrium lattice heating as the origin of the nonlinear electronic response. The phonon-free electronic amplitude $A_{\rm elec}$ reveals a fluence-reversal fingerprint: with cooling from the hybridization-crossover regime, the response evolves from weak-linear behavior to Rothwarf-Taylor-like bottleneck suppression and finally to anomalous high-fluence enhancement at the lowest temperatures. This reversal cannot be accounted for by a rigid fixed-gap bottleneck alone and instead identifies an ultrafast optical signature of photoinduced redistribution of a fragile Kondo-hybridized electronic response.

Superconducting domes in strongly correlated electronic systems are often accompanied by charge density waves and peculiar features in the electronic structure. The appearance of a pseudogap, in particular, and its relation to charge density waves remains insufficiently understood. Here, we investigate the electronic Green's function of the underdoped $t$-$t'$-$J$ model of the cuprate superconductors using tensor network algorithms for finite temperature dynamics on cylinders of width $4$. We find the prominent momentum differentiation, the hallmark of the pseudogap, to be strongly dependent on $t'$, which develops into a momentum-dependent opening of a gap upon decreasing temperature, consistent with the formation of a Fermi arc at intermediate temperatures. The nodal gap around ${\bf k}=(π/2,π/2)$ closes and fills with increasing temperature as coherent stripe order melts into a regime of fluctuating charge clusters.

Recently, layered corrected kagome metal CsCr3Sb5 have garnered significant attention attributed to its flat bands near the Fermi level (EF) and altermagnetic ground state [ Yi Liu et al., Nature 632, 1032 (2024)]. However, the van Hove singularities (vHSs) in bulk CsCr3Sb5 are far away from the EF, while an effective modulation of VHS toward the EF is essential for exploring intriguing electron transport properties. In this work, using first-principles calculations, we investigate electronic structures of two-dimensional (2D) Cr3Sb5 and CsCr3Sb5 monolayers which may be mechanically exfoliated from bulk materials. Notably, it is revealed that both flat bands and vHSs simultaneously reside in close proximity to the EF in Cr3Sb5 monolayer, signifying enhanced electronic correlations. Importantly, a tensile strain further shifts the incipient flat bands and vHSs of two monolayers simultaneously to the vicinity of the EF, suggesting strain tunable electronic correlations and concomitant quantum effects. Furthermore, altermagnetic ground state is also revealed due to retained mirror symmetry between two sublattices in these two monolayers. Thus, this work advances understanding and modulations of electronic properties of 2D CsCr3Sb5 monolayers, strengthening their great potential for exploring unconventional quantum phenomena and altermagnetism.

An attempt is made to give a brief but coherent account of the situation of the theoretical ideas in addressing the mechanism of superconductivity in cuprate high-Tc superconductors. Specifically, the idea of superconductivity from repulsive interactions is discussed as it is gaining ground since the `consensus' paper was written in 2015\cite{kei}. The challenges it faces is also discussed. Three main schools of thought are presented, and an experimental result of 2022 pertaining to Anderson's super-exchange mechanism is also discussed. An updated list of Anderson's ``dogmas" is also presented, as after year 2000, many other universally applicable experimental facts has been discovered. The ``dogmas" are universal facts which are distilled from a variety of complex experimental results, and highlights the key findings that seems to be central to the mechanism of superconductivity in cuprates. These are discussed as a commemoration of 40 years of high-Tc cuprate research.

MnF$_2$ is widely regarded as a candidate altermagnet, but the magnitude and implications of its altermagnetic band splitting remain debated. Using electronic-structure calculations, we construct minimal models that capture the magnetic and electronic properties of MnF$_2$. These models show that the parameters governing the chiral magnon splitting and the spin splitting of the electronic bands are relatively small. Moreover, the electronic system lies in the strong-coupling regime, where most magnetic properties are controlled by the ratio $t/U$ between the characteristic hopping amplitude $t$ and the large on-site Coulomb repulsion $U$. Consequently, all exchange interactions scale as $1/U$, so a small altermagnetic hopping $δt$ produces only a proportionally small exchange term. Upon doping, the altermagnetic contribution to the anomalous Hall effect is likewise suppressed, being smaller than the conventional (non-altermagnetic) contribution by a factor of order $δt/U$. In contrast, the behavior of the conductivity tensor $\hatσ(ω)$ at $\hbar ω\sim U$ differs qualitatively, because $δt$ enters the energies of interband optical transitions \emph{directly} rather than through the reduced ratio $δt/U$. This contribution strongly reshapes $\hatσ(ω)$, leading to a dramatic enhancement of the magneto-optical response.

We study the circularly driven Falicov-Kimball model on a honeycomb lattice within real space Floquet dynamical mean field theory (DMFT). The noninteracting version of this model has been realized experimentally. The noninteracting system hosts an effective Haldane phase at large driving frequencies, while at intermediate frequencies it hosts an anomalous topological phase. We study the effect of two-body interactions $U$ on the stability of these phases. We find that charge pumping does not remain quantized upon increasing $U$, despite the presence of edge modes in the spectrum. This can be attributed to the broadening of the edge modes due to interaction. We also calculate the rate of energy dissipation into the bath and find remarkably different behaviour in the two regimes.

We investigate topological electronic responses in a kagome altermagnetic metal hosting a compensated coplanar $120^\circ$ magnetic texture. Using a minimal tight-binding model incorporating nearest-neighbor hopping, noncollinear exchange coupling, intrinsic spin--orbit coupling, and an emergent orbital chiral flux, we demonstrate that frustrated kagome altermagnets provide a natural platform for realizing momentum-dependent spin splitting and Berry-curvature engineering without net magnetization. The noncollinear exchange field alone generates pronounced altermagnetic spin splitting and spin-polarized Fermi surfaces despite the absence of relativistic effects. However, for a strictly coplanar magnetic state, the system preserves a hidden antiunitary symmetry $\mathcal{T}C_{2z}$, which enforces identically vanishing Berry curvature even in the presence of sizeable spin--orbit coupling. We show that finite Berry curvature emerges only after introducing an orbital chiral flux term that breaks the hidden symmetry and generates effective momentum-space gauge fields analogous to a Haldane-type orbital flux. Remarkably, this mechanism produces local Berry curvature hot spots even in the complete absence of spin--orbit coupling and scalar spin chirality, establishing a purely orbital route toward topological altermagnetism. By systematically analyzing the anomalous Hall conductivity as a function of exchange coupling, spin--orbit interaction, and chiral flux, we identify a hierarchy of competing energy scales governing the transition from a symmetry-protected altermagnetic metal to a topological altermagnetic phase with strong Hall response. Our results demonstrate that frustrated kagome altermagnets constitute a versatile platform for engineering topological transport, Berry curvature, and spin-selective electronic structure in compensated magnetic systems.

We present an \textit{ab initio} framework for predicting resonant inelastic X-ray scattering (RIXS) in optically pumped materials. Our methodology is based on the Kramers-Heisenberg formula for the double-differential cross section formulated using the results of the Bethe-Salpeter equation (BSE) from many-body perturbation theory. To extend this approach to the time domain, we incorporate non-equilibrium charge-carrier distributions obtained from real-time, time-dependent density-functional theory (RT-TDDFT). Generalizing the RIXS implementation with respect to arbitrary polarizations, allows us to consider different orientations of incoming and outgoing light. We demonstrate our method's capabilities by studying RIXS at the K-edge of graphite for various non-equilibrium charge-carrier distributions, representing different delay times after optical pumping. Our results reveal angular dependencies in $π$- and $σ$-orbital-derived spectral regions, in good agreement with experiment.

Trotter product formulas are a fundamental class of methods for Hamiltonian simulation, particularly attractive due to their low qubit requirements. However, they are often overlooked for use with fault-tolerant quantum algorithms, because of their perceived higher gate counts and the difficulty of estimating Trotter error. Here, we introduce Symmetry-Protected Randomized near-Integrable Trotter (SPRINT) formulas, a framework for building optimized product formulas for electronic structure Hamiltonians widely used in quantum chemistry. SPRINT integrates a generalization of classical near-integrability, randomization, symmetry protection, use of QROM, and other techniques into a thoroughly optimized methodology for Hamiltonian simulation. When applied to concrete simulation tasks, we find SPRINT leads to substantial reduction in gate count compared to previous approaches. Alongside SPRINT, we introduce and analyze a Generalized Rank Decomposition (GRADE) of electronic Hamiltonians that generalizes previous factorization methods. We apply these techniques to the task of simulating the X-ray absorption spectrum of Li$_4$Mn$_2$O, a candidate battery cathode material, leveraging recent advances in tight Trotter error estimation to carefully identify the best version of SPRINT for this problem. Using a Trotter error estimation tool developed in the PennyLane software platform, we show that SPRINT reduces the Toffoli gate cost by a factor of $4.5$ relative to the previous state of the art for this problem, with a gate cost only $\times 2.5$ higher than qubitization, while requiring a dramatic $\times 5.5$ fewer logical qubits. These results establish well-designed Trotter product formulas as an attractive Hamiltonian simulation method for industrially relevant problems in chemistry and materials science.

We investigate the ground-state phase diagram of repulsively interacting bosons on a two-leg ladder threaded by a uniform artificial magnetic flux, using the cluster Gutzwiller mean-field method. In the strong-rung-coupling regime, self-consistent calculations are performed on a $2\times4$ cluster. By analyzing the superfluid order parameter, leg-resolved currents, chiral current, the current ratio on adjacent legs, and the density imbalance between the two legs, we distinguish Mott-insulating from superfluid regimes and characterize the observed states as Meissner-like, vortex-like (superfluid or Mott insulating), or biased-ladder. In regions overlapping with previous DMRG studies, our results qualitatively agree with the established phase structure, demonstrating that the cluster Gutzwiller approach balances computational efficiency and physical accuracy. We then construct the first grand-canonical $t$--$μ$ phase diagrams for this system, revealing how the magnetic flux modifies the shape, tilt, and extent of the Mott lobes. We further explore previously inaccessible regimes, including higher fillings $ρ\gtrsim1$ and the intermediate interaction window $U/t\in[7.69,9.09]$. Special attention is paid to $\varphi=π$, where the effective triangular-ladder mapping becomes singular. Owing to the equivalence of $\varphi=π$ and $-π$ modulo $2π$, a combined symmetry forbids net chiral currents, leading to a nonchiral Mott-insulating state, in contrast to the chiral-superfluid tendency expected away from $\varphi=π$. Our results offer a computationally efficient route for mapping the global phase structure of bosonic flux ladders and provide guidance for future ultracold-atom experiments in artificial gauge fields.

A new family of high-temperature superconductors was recently discovered in the $n=2,3$ Ruddlesden-Popper nickelates, where superconductivity emerges concomitant with suppression of parent density waves and structural octahedral rotations under hydrostatic pressure. Intriguingly, compressive strain mimics the structural effects of pressure in the $n=2$ phase, yielding ambient-pressure superconductivity. However, analogous strain-stabilized superconductivity has not been realized in the $n=3$. Here, we use atomically-precise synthesis, transport, picoscale electron microscopy, and synchrotron X-ray diffraction to probe $n=3$ La$_4$Ni$_3$O$_{10}$ thin films. Although compressive strain suppresses density wave order, we do not observe superconductivity even under the largest strain state. Importantly, we identify a structural distortion unique to strained $n=3$ thin films that may inhibit superconductivity: persistent, layer-inequivalent octahedral rotations around the $c$-axis. Our results highlight key differences between the $n=3$ and $n=2$ systems, suggesting that ambient-pressure superconductivity in the $n=3$ may require new methods beyond epitaxial strain engineering.

The discovery of superconductivity in hole-doped Ni$^{1+}$ systems with "infinite layer" NiO$_2$ square-lattices analogous to the Cu$^{2+}$ CaCuO$_2$ cuprate has renewed conflicting pictures of the Cu$^{2+}$$-$Ni$^{1+}$ similarity or distinction. Recent synthesis of formal Ni$^{1+}$ La$_3$Ni$_{2}$O$_{5}$F with two infinite NiO$_{2}$ layers per cell provides a novel member of this class. First principles density functional theory studies reveal an interstitial density derived single band $E^*$ in three layers unrelated to any atom, which provides self-doping to a Ni$^{1.09+}$ ion.The blocking La(O/F)La provides isolation of the NiO$_2$ bilayer and an interstitial $E^*$ density to strictly two-dimensional electronic and magnetic systems. Calculations of magnetic tendencies reveals behavior unlike previous nickelates, including vanishing susceptibility up to a large magnetic field. Two dimensional fluctuations and self-doping away from half-filling can account for the lack of observation of a magnetic transition.

Superconductivity in cuprates emerges out of a complex normal state that hosts density waves, pseudogap physics, and strange metal properties. Here, we access this normal state by synchronizing free-electron laser x-rays with high-magnetic-field pulses up to 44 T. We observe a linear increase in charge order amplitude and correlation length that persists far above the vortex melting transition. This behavior is incompatible with standard phase competition between charge order and superconductivity. By means of conventional hard x-ray diffraction and magnetostriction, we show that applied fields also enhance monoclinic lattice distortions. However, this magnetoelastic response is weaker and an epiphenomenon of the stripe order enhancement. Combined with recent observations of field-linear spin freezing, our results point to a direct coupling between magnetic field and the spin component of stripe order in the high-field normal state -- a mechanism independent of superconductivity suppression that has so far remained hidden from scattering probes.

Quantum impurity models with strong electron correlations, such as the paradigmatic Anderson Impurity Model (AIM), are central to our understanding of a range of physical phenomena including local moment formation, Coulomb blockade and Kondo screening. They describe magnetic atoms and molecules on surfaces, quantum dot circuits, and correlated materials through dynamical mean field theory. The physics of such systems in strongly non-equilibrium conditions is particularly complex and challenging to capture, whereas Gaussian models of free fermions can be easily solved. Here we show that the time-evolving dynamics of the AIM after a quench can be described by a completely non-interacting version of the model, at the expense of coupling to additional static auxiliary degrees of freedom. Starting from the full solution of the quenched AIM using ED and DMRG, we study the properties of this mapping using numerical optimization, and uncover intriguing structure in the auxiliary system. The method allows us to understand interacting non-equilibrium dynamics through the simpler lens of an effective non-interacting system of larger dimension.

We investigate a frustrated four-leg spin-$1/2$ ladder using density matrix renormalization group calculations. The uniform system displays three regimes: short-range antiferromagnetic legs, short-range ferromagnetic legs, and an effective spin-2 Heisenberg chain, separated by a crossover and a first-order transition. The spin-2 regime is confirmed through its finite string order parameter, edge-localized excitations, and excellent agreement with a projected $S_r=2$ effective Hamiltonian. Recasting the model as two frustrated two-leg ladders coupled by rung and diagonal interactions, we track how the trivial and Haldane phases of an isolated ladder evolve as interladder couplings are introduced. The resulting phase diagrams reveal crossover and first-order lines whose locations are captured by the spin-2 projection and show how singlet- and triplet-dominated regimes reorganize when two ladders merge into a four-leg structure, clarifying the emergence of effective spin-1 versus spin-2 behavior.

Recent theoretical progress has established conditions under which machine learning models can efficiently predict ground-state properties of gapped local Hamiltonians when trained on quantum-generated data. Previous experimental demonstrations in this paradigm, however, have largely been limited to small systems or highly structured states, due to the difficulty of preparing many-body ground states on quantum processors. In this work, we demonstrate learning from experimental quantum data generated from approximate ground states of the two-dimensional Heisenberg XXZ model with system sizes up to 115 qubits. We construct a dataset of single-site expectation values, two-point correlations, and 12-body loop correlations across the antiferromagnetic phase. We then train neural networks on this data and show that they can accurately predict spatially resolved observables for previously unseen Hamiltonian parameters, both within the training distribution and in an out-of-distribution regime approaching the phase boundary. Our results demonstrate the practical realization of learning from quantum data for an interacting two-dimensional many-body system at scale, motivating a path toward regimes where quantum processors could provide training data beyond the reach of classical approximation methods.

A new moiré material platform was recently proposed based on twisting two-dimensional atomic monolayers whose low-energy states lie at the three M-points of the Brillouin Zone. Continuum and ab initio modeling suggest that electrons in the conduction bands of these materials realize three-valley Hubbard models with valley-selective, quasi-one-dimensional hopping. Remarkably, the onsite Hubbard repulsion is almost $U(6)$-symmetric without fine-tuning. Here, we show that this class of systems naturally admits sign-free determinantal Quantum Monte Carlo simulations at a filling of three electrons per moiré unit cell. We use these to explore the phase diagram for interactions of various strengths and $U(6)$-breaking anisotropies. We show that for near-isotropic interactions as relevant to, e.g., AA-stacked twisted SnSe$_2$, the system exhibits an extended intermediate-coupling regime in which local-moment formation and itinerancy compete, and the crossover to a putative low-temperature ordered state can be understood in terms of fluctuating $U(6)$ local moments. We argue that many of these features persist beyond the idealized sign-problem-free limit.

Motivated by the recent discovery of superconductivity at 16.3 K in layered oxychalcogenide Na$_{2-x}$V$_2$Se$_2$O we investigate pairing instabilities in the broader family of layered materials composed of V$_2$O planes, believed to exhibit altermagnetic order in their monolayer form. Even though the bulk family members KV$_2$Se$_2$O and Rb$_{1-δ}$V$_{2}$Te$_{2}$O are likely conventional antiferromagnets that show only surface altermagnetism, our analysis predicts exotic equal-spin triplet superconductivity as the dominant pairing instability in these materials. This is a consequence of their unique magnetic and sublattice structure that renders electron bands incompatible with conventional spin-singlet pairing. The predicted triplet superconducting phases are topologically non-trivial and capable of supporting spin-polarized persistent currents, properties potentially useful in technological applications.

The interplay between quantum criticality and Fermi surface reconstruction is central to elucidating the phase diagram of high-temperature cuprate superconductors. While studies on electron-doped T'-structure cuprates suggest an antiferromagnetic origin of this reconstruction, quantitative consensus has been hindered by apical oxygen instabilities and uncontrolled oxygen vacancies. Here, we overcome these limitations by utilizing ozone-assisted molecular beam epitaxy to synthesize high-quality, oxygen-stoichiometric thin films of infinite-layer cuprate Sr1-xNdxCuO2 across its entire superconducting dome. Hall transport measurements reveal a sharp carrier-type transition signaling a Fermi surface reconstruction at a critical doping xc ~ 0.155. We show that a spin-density-wave tight-binding model quantitatively reproduces the transport evolution, supporting an antiferromagnetic origin of this quantum phase transition. Furthermore, upon suppressing superconductivity with magnetic fields, the normal-state resistivity exhibits a pristine strange metal behavior that persists down to 2 K in the vicinity of xc. Our findings establish an intrinsic, universal antiferromagnetic quantum criticality in electron-doped cuprates, positioning the structurally simplest infinite-layer cuprates as a clean benchmark platform for theories of unconventional superconductivity.

We investigate the ground-state phase diagram of a frustrated spin-1/2 Heisenberg ladder in the transverse magnetic field with an anisotropic inter-rung exchange coupling $α$. This bond-anisotropic parameter continuously interpolates among the Ising-like limit $α=0$, the isotropic point $α=1$, and the XY-dominated regime $α\gg 1$. Using density-matrix renormalization group calculations within a matrix-product-state framework, we analyze bipartite entanglement, magnetization, and spin correlations to characterize the emergent quantum phases. We identify six distinct ground-state phases, including rung-singlet, Haldane-like, Tomonaga-Luttinger liquid, canted Ising-ordered, XY-polarized, and ordered ferromagnetic states. While magnetization and local correlations provide a first insight into phase classifications, the finite-size scaling of the entanglement measures offers a more sensitive and unified diagnostics of both gapped and gapless regimes. In the gapless regime, we extract the central charge $c \simeq 1$, confirming Tomonaga-Luttinger liquid behavior, while at the transition between the canted Ising and ferromagnetic phases, we observe $c \simeq 1/2$, consistent with the Ising universality. Finally, we find that both the Haldane-like behavior and the extended critical Tomonaga-Luttinger liquid regime are strongly confined to the vicinity of the isotropic point $α= 1$.

The Kitaev honeycomb model has attracted significant interest due to its quantum spin liquid ground state, fractionalized Majorana excitations, and topological properties. Motivated by these features, we introduce a quasi-one-dimensional Kitaev-like spin chain derived from a truncated honeycomb geometry. The resulting Majorana band structure contains both dispersive and flat bands, with a gap closing at a critical anisotropy that separates trivial and chiral topological phases. Under open boundary conditions, the topological regime hosts zero-energy edge modes protected by chiral symmetry, placing the system in the BDI symmetry class with a quantized winding number. In the excited spectrum, localized plaquette modes emerge near zero energy and can be tuned by introducing domains of negative plaquettes. The dynamical spin correlation function reveals broad continua associated with fractionalized excitations, together with characteristic low-energy spectral weight from edge Majorana modes. These features distinguish the present system from conventional Heisenberg spin chains and provide experimentally relevant fingerprints of chiral topological order. Our model thus establishes a conceptual bridge between the two-dimensional Kitaev spin liquid and Kitaev's one-dimensional quantum wire.

Deconfined quantum critical points (DQCPs) describe continuous transitions between ordered phases beyond the Landau paradigm. A simple example is the Néel antiferromagnet (AFM) to valence bond solid (VBS) transition in a 1D antiferromagnetic $J_1-J_2$ model. In analogy to the spin-Peierls instability of critical spin chains, DQCPs are predicted to be unstable towards lattice distortions below a critical phonon frequency. In this work, we use tensor network simulations to investigate this instability in the antiferromagnetic $J_1-J_2$ model coupled to lattice vibrations. We confirm the stability of DQCP for large phonon frequencies and demonstrate that the transition turns strongly first-order below a critical frequency. The instability is caused by a reduction of the Luttinger parameter due to spin-phonon interactions and we identify the effective theory of the behavior as the double sine-Gordon model. The same effective theory is known to describe the classical Ashkin-Teller model, which enables us to show that the critical endpoint is in the four-state Potts universality class. Furthermore, we provide quantitative numerical scaling results for the phonon spectral function, offering an experimental signature to probe DQCP-phonon coupling in low-dimensional materials.

How do we detect topological phases from a single ground state? Entanglement entropy and spectrum have long been the standard tools -- but the reduced density matrix (RDM) itself contains far more information. We show that the RDM of a 2D topologically ordered system, expressed at the entanglement cut, realizes a 1D mixed-state phase. For the $\mathbb{Z}_2$ toric code phase, it is a 1D $\mathbb{Z}_2$ strong-to-weak spontaneous symmetry breaking (SW-SSB) phase, where deconfinement of anyons manifests as the short-range correlation of both $\mathbb{Z}_2$ charge and $\mathbb{Z}_2$ domain-wall in the RDM. The bulk $e$-$m$ duality translates into a Kramers--Wannier self-duality of the SW-SSB phase. Extending the framework to gapped $\mathbb{Z}_2$ spin liquids, the global spin-rotation symmetry manifests as an additional weak symmetry for the 1D RDM. Spin-$\frac{1}{2}$ spinons result in a cusp on the disorder parameter of spin-rotation at $θ=π$, providing a direct, ground-state signature of symmetry fractionalization. We verify this prediction analytically using the matrix product density operator formalism and numerically for the kagome-lattice resonating valence bond state. The proposed observable requires only a single ground-state wavefunction, making it amenable to quantum simulation platforms.

Spin-wave excitations provide a central probe of magnetic order and electronic correlations in strongly correlated materials. In this work, we develop an adiabatic theory of spin dynamics by combining the Niu-Kleinman formalism with Kotliar-Ruckenstein slave-boson theory (NK+KRSB). For each frozen spin configuration, the constrained slave-boson saddle point is solved self-consistently, allowing the Berry-curvature matrix and energy Hessian entering the linearized adiabatic equations of motion to be extracted directly. Applied to the half-filled single-orbital Hubbard model, the resulting spin-wave dispersion shows substantially improved agreement with determinant quantum Monte Carlo benchmarks compared with the random phase approximation and closely approaches results from the time-dependent Gutzwiller approximation. We further extend the method to a two-orbital model of $\mathrm{La}_2\mathrm{NiO}_4$, demonstrating its applicability to realistic multi-orbital correlated systems. Because the approach only requires saddle-point solutions near the magnetic ground state, it remains computationally efficient while incorporating strong-correlation effects beyond conventional weak-coupling descriptions, providing a practical framework for studying low-energy spin excitations in correlated quantum materials.

Altermagnets are collinear magnetically ordered states that exhibit momentum-dependent spin splitting in the absence of net magnetization and spin-orbit coupling (SOC). Related spin-splitting patterns, however, can also emerge in noncollinear magnetic systems with large SOC. Here we show, via a microscopic model, that the all-in-all-out (AIAO) state in pyrochlore iridates constitutes a noncollinear counterpart of a $d$-wave altermagnet stabilized by strong SOC. Starting from a microscopic $j_{\mathrm{eff}} = 1/2$ tight-binding model on the pyrochlore lattice, we demonstrate that electronic interactions favor the AIAO phase and analyze its symmetry properties. We show that the AIAO order parameter transforms as an $A_{2g}^{-}$ octupolar magnetic moment, breaking time-reversal symmetry while preserving inversion and zero net magnetization. Using group-theory analysis and mean-field calculations, we demonstrate that this symmetry enforces both a spin-polarized momentum-dependent lifting of band degeneracies that is similar to that of a collinear $d$-wave cubic altermagnet, but also a band splitting at zero-momentum. We show that the latter feature is captured by a low-energy model similar to the Luttinger-Kohn model for cubic semiconductors. Our results identify pyrochlore iridates as a platform for noncollinear counterparts of altermagnetism and provide a general symmetry framework for spin-split phenomena in spin-orbit coupled materials.

Altermagnetism defines a third fundamental class of collinear magnetic order, featuring compensated magnetic moments with antiparallel spin alignment, yet lifted Kramers degeneracy without the need for relativistic spin-orbit coupling. Its ability to host spin-polarized electronic bands and unconventional chiral magnons makes it a promising platform for functional materials. However, experimental verification has proven challenging; while circular dichroism in resonant inelastic X-ray scattering (RIXS) has been suggested as a signature of chiral magnons, it remains controversial whether this effect is an intrinsic property of altermagnetism or an artifact of experimental geometry. In this work, we resolve this debate and provide unambiguous experimental evidence of $d$-wave altermagnetism in the strongly correlated Lieb-lattice magnet La$_2$O$_3$Mn$_2$Se$_2$. The RIXS spectra exhibit a $d$-wave-symmetry circular dichroism in the magnetic excitations that vanishes in the paramagnetic phase. Through RIXS-operator symmetry analysis and exact-diagonalization calculations, we prove that the observed dichroism is a direct consequence of altermagnetic symmetry constraints, independent of magnon branch splitting. Our results provide definitive evidence for the experimental realization of $d$-wave altermagnetism in La$_2$O$_3$Mn$_2$Se$_2$ and establish circularly polarized RIXS as a highly symmetry-sensitive spectroscopic framework for detecting magnetic phases that evade conventional probes.

The interaction-driven crossover from a Slater insulator to a Mott insulator in the Néel-ordered ground state of the Hubbard model is a fundamental paradigm of strongly correlated electrons, yet its quantitative characterization has remained elusive. Here we establish a clear and experimentally accessible thermal criterion for this crossover via the sign change of the temperature derivative of double occupancy, $(\partial D/\partial T)_U$, near zero temperature. In the Slater regime, $(\partial D/\partial T)_U>0$ reflects the major role of charge fluctuations; in the Mott regime, the anomalous $(\partial D/\partial T)_U<0$, a manifestation of the Pomeranchuk effect, signals the dominance of low-energy spin superexchange physics. Using exact diagonalization and {\it numerically exact} quantum Monte Carlo simulations, we demonstrate that this criterion determines the crossover boundary at $U_{\rm cross}/t=4.0(2)$ for the half-filled two-dimensional Hubbard model. Furthermore, we obtain a consistent boundary independently from the maximum in the antiferromagnetic correlation length, which also arises from the superexchange physics. These two thermal signatures are theoretically unified through the local minimum of thermal entropy versus interaction $U$ at low temperatures. Our results offer a direct, measurable, and physically intuitive framework to identify the Slater-Mott crossover in optical lattice experiments.

Scaling laws, the power-law relations between loss, architecture size, and compute observed in modern neural networks, offer a quantitative way to characterize the complexity of a learning problem, with the exponent governing the decay of the loss reflecting how rapidly additional resources translate into improved accuracy, and thus how hard the target is to learn. Whether an analogous framework can characterize the complexity of physical problems remains open. We address this question for Neural-Network Quantum States, a leading variational approach for strongly correlated quantum many-body systems. Using transformer wave functions to approximate ground states of the $J_1$-$J_2$ Heisenberg model on triangular and square lattices with up to $20\times 20$ sites, we find that the $V$-score, a measure of accuracy of a variational state, decays as a power law in training compute. Under an appropriate rescaling of compute, results for different system sizes collapse onto a single curve, analogous to scaling collapse in critical phenomena. The resulting power law is, to a good approximation, independent of the number of sites, showing that the transformer Ansatz is size-consistent for the systems considered. The exponent decreases systematically with frustration, identifying it as a quantitative measure of representational difficulty of the ground state and establishing scaling laws as a general framework for benchmarking variational ansätze.

Rare-earth quantum materials provide a promising platform for emergent phenomena ranging from quantum spin liquids with long-range entanglement to topological magnetic textures. However, the strong spin-orbit coupling that stabilizes their low-energy pseudospin degrees of freedom also tends to generate strongly anisotropic exchange interactions, complicating the realization of quasi-isotropic Heisenberg magnetism. Here we investigate the microscopic origin of superexchange in $\mathrm{Ce}^{3+}$- and $\mathrm{Yb}^{3+}$-based insulators with edge-sharing octahedral geometry. Using degenerate perturbation theory for a multiorbital Hubbard model, we show that isotropic exchange originates predominantly from virtual hopping within the ground-state Kramers doublet, whereas anisotropic interactions arise primarily from processes involving excited multiplets. This leads to a simple orbital design principle: quasi-isotropic exchange is promoted when the ground-state doublet has a strong maximal-angular-momentum character with respect to the quantization axis perpendicular to the superexchange plane spanned by rare-earth and ligand ions. We demonstrate this mechanism for both ideal and distorted geometries and show that it is broadly consistent with experimentally studied Yb-based insulators. Our results establish a practical framework for engineering quasi-isotropic interactions in rare-earth quantum materials.

We present a universal framework to harness fermionic ultracold atom platforms for quantum simulation, showing how variational algorithms on existing hardware can simulate many-body systems well beyond the hardware's native Hamiltonian. Our analysis provides evidence that one can quantum simulate the ground-state properties of a broad class of gapless target Hamiltonians of local observables in a quantum evolution time that grows polynomially with the inverse relative error, $T\sim O(\mathrm{poly}(1/ε))$ up to logarithmic corrections, offering an exponential speedup over na{ï}ve classical algorithms such as exact diagonalization. We provide numerical evidence and theoretical argument that this holds for energy density, density-density, and spin-spin correlations in three qualitatively distinct models -- the repulsive Hubbard model; a Hubbard model augmented with nearest-neighbor attractive interactions, which introduces the phenomenon of pairing; and the Hofstadter-Hubbard model, which introduces a gauge field and fractional quantum Hall physics. This work demonstrates quantum simulation using current fermionic platforms far beyond the models natively implemented in the hardware.

The numerical renormalization group (NRG) approach was developed to identify and quantify different equilibrium regimes of quantum impurity systems (QISs) with unprecedented accuracy by a tailored finite size representation. Out of equilibrium, the steady-state density operator is not of the Boltzmannian form but one that is determined by the imposed boundary conditions. We extend the NRG to the nonequilibrium setting by augmenting each Wilson site with a reservoir, whose coupling functions are calculated via a continuous fraction expansion in order to recover the continuum limit exactly. The nonequilibrium parameters such as a finite bias as well as a finite temperature gradient enters through the Bloch-Redfield tensor (BRT), whose zero eigenvector gives the steady-state density operator. We used the resulting open chain full density matrix (OC-FDM) approach with an effective single lead description to investigate the charge and spin transport through a quantum dot (QD) under finite bias and temperature gradient. The influence of lead asymmetry and an external magnetic field on the transport properties are also studied. For completeness, we have also investigated local properties of the QD, such as charge fluctuations and find excellent agreement with real-time quantum Monte Carlo (RT-QMC) data. The OC-FDM approach was able to explore Kondo energy scales as low as $T_K/D \approx 10^{-8}$ in the non-equilibrium regime, as well as show convergence with the established equilibrium benchmarks, such as a quantitative agreement with full density matrix numerical renormalization group (FDM-NRG) and Fermi-liquid scaling at small bias and temperature. Owing to the effective single lead description, a single OC-FDM data point takes orders of magnitude less time on a standard laptop, compared to other state-of-the-art numerical methods.

By using unbiased determinant quantum Monte Carlo simulations, we investigate the interplay between superconductivity and charge density wave order in the Holstein model on a honeycomb lattice with next-nearest-neighbor hopping \(t^{\prime}\). We find that a finite negative \(t^{\prime}\) enhances \(s\)-wave superconducting pairing susceptibility near the van Hove fillings in the weak electron-phonon coupling regime, while it suppresses superconductivity and promotes charge density wave order at intermediate electron-phonon coupling strengths. The effect of \(t^{\prime}\) on a charge density wave is filling-dependent: It suppresses the charge density wave at half filling but enhances it near the van Hove singularities. A spectral analysis reveals the opening of a gap at low temperatures, highlighting the competitive relationship between superconducting and charge density wave orders mediated by electron-phonon coupling and tuned by \(t^{\prime}\).

Linear-in-temperature resistivity is a hallmark for strange metallic transport, and appears universally in many strongly correlated electron systems. However, the focus on the longitudinal channel often overshadows the profound microscopic insights contained within the transverse response. Here, we utilize numerically exact determinantal quantum Monte Carlo simulations of the doped Hubbard model in a magnetic field to calculate longitudinal and transverse transport. We demonstrate that while the resistivity is robustly $T$-linear across parameter sets, the Hall response is highly sensitive to particle-hole asymmetry, Fermi surface topology, and many-body correlation effects. Specifically, the combination of these effects determine a crossover scale in which the system becomes quantum-coherent, and is reflected in the Hall conductivity. Our results demonstrate that while the $T$-linearity in resistivity appears universal, the Hall response reveals a crossover from semi-classical to quantum-coherent transport otherwise masked in the longitudinal channel.

High harmonic generation (HHG) is a sensitive probe for investigating electronic structures and dynamics of materials and a source for attosecond pulses. In particular, HHG with terahertz (THz) light can enable probing of nonlinear responses in correlated materials arising from low-energy many-body interactions. However, THz HHG studies have so far largely focused on topological materials and superconductors, leaving out other potential material systems which could also become efficient THz HHG sources. Here, we report THz third harmonic generation (THG) in rare-earth nickelates -- a prototype material for exploring the Mott insulator-metal transition and related technological applications. We find that the THG amplitude is highly sensitive to the strengths of electronic and magnetic phases of nickelates. In films with sharp phase-transitions, the local maximum and minimum in the temperature-dependent THG amplitude coincide with insulator-metal and magnetic transition temperatures, respectively. While in films with weaker transitions, these features shift toward lower temperatures or even monotonous THG enhancement is observed down to low temperatures. We developed a generalized theory for THz harmonic generation in negative charge-transfer insulators and outlined strategies to enhance the THz nonlinearities further. Our study broadens the scope of THz HHG studies and related applications to strongly correlated materials.

The emergence of non-relativistic spin splitting (NRSS) in altermagnetic systems has introduced a new paradigm in antiferromagnets with vanishing net magnetization. Although sliding-induced valley-polarized phases have recently been demonstrated in two-dimensional altermagnets, the observed valley-polarized state represents only a partial manifestation of altermagnetism, and a comprehensive classification based on spinsplitting characteristics remains lacking. Here, using first-principles calculations, general stacking theory, and spin-Laue symmetry analysis, we propose a coupled quasialtermagnetic state representing a distinct subclass of altermagnetism, in which reversible type-IV NRSS is controlled through interlayer sliding. Accordingly, the sliding-induced phases are classified into two categories: altermagnetic and quasi-altermagnetic states. We establish a direct correspondence between reciprocal-space spin splitting and real-space switching between the two quasi-altermagnetic states. Importantly, the spin-polarized bands in these states remain spin split at Γ point even in the absence of spin-orbit coupling (SOC), distinguishing them within the proposed classification framework. To demonstrate the interplay between altermagnetic and quasi-altermagnetic states, we investigate the two-dimensional Lieb-lattice material Mn2WS4 and its Janus derivative Mn2WS2Se2, analysing how changes in the local environment influence the different magnetic phases. Importantly, the underlying mechanism is broadly applicable to a wide class of twodimensional square-lattice systems. We further investigate the effects of SOC, focusing on spin texture and transport signatures in coupled quasi-altermagnetic states.

Symmetric mass generation is the interaction-induced opening of a fermion gap without spontaneous symmetry breaking. The anomaly-free 3-4-5-0 model of Wang and Wen provides a minimal one-dimensional setting for this phenomenon, but a direct lattice realization faces two obstacles: fermion doubling for local chiral discretizations and perturbative irrelevance of the six-fermion gapping interaction. We address both obstacles. First, we formulate the model on a strictly one-dimensional tangent-fermion lattice, where a nonlocal hopping produces a single chiral branch without a mirror partner while retaining an efficient tensor-network representation. Second, we add a Hubbard-type density-density interaction (Luttinger parameter $K$) that reduces the scaling dimension of the 3-4-5-0 interaction from $5$ to $5K$, making it relevant for $K<2/5$. Density-matrix renormalization group calculations show the opening of an excitation gap in this regime without the appearance of a degenerate ground state, the hallmark of symmetric mass generation.

Quantum impurity models provide a paradigmatic setting for studying Kondo screening, boundary criticality, and impurity entropies. While these phenomena are well understood in unitary systems, their fate in non-Hermitian many-body settings remains largely unexplored. We study a $\mathscr{PT}$-symmetric quantum impurity model consisting of a unitary $SU(2)_1$ Wess--Zumino--Witten bulk coupled to two impurity spins through complex-conjugate boundary Kondo interactions. Using an integrable lattice realization with $\mathscr{PT}$-symmetric boundary impurities, solved by the Bethe Ansatz and benchmarked against finite-temperature matrix-product-state calculations, we determine the impurity contribution to the free energy and entropy. In the Kondo-screened regime, where the spectrum remains entirely real and the impurities are screened by many-body Kondo clouds, we find that the impurity entropy decreases monotonically from $\ln 4$ in the ultraviolet to $0$ in the infrared. This monotonic flow persists despite the nonunitary nature of the boundary interaction, which places the system beyond the standard assumptions of the $g$-theorem.

From the condensed matter physics perspective, the most natural single orbital tight-binding Hamiltonians, and hence the most widely studied, contain two fermionic species, corresponding to spin up and spin down electrons. In cold atom systems, however, SU(N) symmetry, in which $N > 2$ fermionic species reside within a single band, also occurs. In order to understand such experiments, the SU(N) Hubbard model has been increasingly studied. Here we present determinant Quantum Monte Carlo simulations of the SU(N) {\it Holstein} Hamiltonian, in which $N$ fermionic species couple to a single local phonon mode. We show that at half filling it has an insulating charge density wave phase (CDW) at low temperatures, in which empty sites alternate with sites with $N$ particles. We determine the $N=3$ CDW phase diagram in the temperature, $T$, versus electron-phonon coupling, $α$, plane at fixed phonon frequency $ω_0$ and half-filling $ρ=1.5$. The critical temperature $T_c$ for $N=3$ can be as high as twice the maximum attainable for $N=2$. We also obtain the $N$ dependence of $T_c$ for a representative, fixed, $ω_0$ and $α$.

Altermagnets are a newly identified class of magnetic materials that combine vanishing net magnetization within the unit cell with spin-split electronic states. Their theoretical description relies on symmetry properties of the bulk band structure. Surfaces and thin films, however, inherently break these symmetries. Here, we investigate the consequences of such symmetry reduction for the electronic structure of bulk altermagnets near the surface and of thin films. When the surface coincides with a symmetry plane of the bulk altermagnetic order, the resulting two-dimensional Brillouin zone exhibits spin-degenerate bands, corresponding to conventional antiferromagnetic behavior. In all other cases, the symmetry of the altermagnetic order is reduced, leading to modified spin splitting. Remarkably, we discover a thin-film geometry of a $g$-wave altermagnet with a particular surface orientation that enables a $d$-wave spin splitting, which is commonly accompanied by the spin-splitter effect, suggesting the functionalization of non-$d$-wave altermagnets by surfaces. Our findings demonstrate that symmetry breaking at surfaces and in thin films fundamentally reshapes altermagnetic spin textures, providing a tunable platform for controlling spin-dependent electronic phenomena.

Altermagnetism features strong momentum-dependent spin splitting despite zero net magnetization, offering a transformative platform for next-generation spintronics. However, the nonvolatile and deterministic switching between its two equivalent spin-splitting states remains a fundamental bottleneck. Here, we propose a universal layer-engineering paradigm to achieve symmetry-enforced ferroelectric switching of two-dimensional altermagnetism. By sandwiching a conventional antiferromagnetic monolayer between two identical ferroelectric layers, the out-of-plane polarization cleanly breaks the spatial symmetry to induce robust altermagnetic splitting. Crucially, the global combined parity-time symmetry dictates that reversing the ferroelectric polarization exactly inverts the altermagnetic spin-splitting pattern. We rigorously validate this mechanism in the In2Se3/MnPTe3/In2Se3 trilayer using first-principles calculations. As a direct consequence, the ferroelectrically driven spin-splitting reversal deterministically flips the anomalous Hall effect signal, providing an unambiguous transport fingerprint to electrically distinguish the two altermagnetic states. Unconstrained by the stringent symmetry requirements of intrinsic single-phase materials, our findings establish a versatile physical framework for electrically addressable altermagnetic spintronics.

Quantum geometry shapes the effective mass of Bloch particles through the geometric properties of single-particle states. Here we show that this principle extends to paired states. Starting from a generic multiband Hubbard model, we derive an exact effective-mass theorem for two-body bound states and its many-body counterpart for Cooper pairs near the critical temperature within Gaussian fluctuation theory. In both cases, the inverse effective mass separates into a ``conventional'' band-structure contribution and a new geometric contribution, pair quantum geometry, governed by quantum metrics on the pairing manifold, which becomes nontrivial when pairing is non-uniform across sublattices. In the many-body setting, analytic continuation renders the fluctuation kernel non-Hermitian, producing a biorthogonal pair geometry and a generally complex Cooper-pair effective mass whose imaginary part reflects Landau damping. Exact calculations on one-, two-, and three-dimensional lattice models show that pair quantum geometry can make quantitatively significant contributions to the effective mass. These results establish pair quantum geometry as a fundamental ingredient of superconductivity beyond conventional band geometry.

In the Hertz theory, a quantum critical metal is described by the coupling of a Fermi surface to fluctuations of a Landau-damped bosonic field $φ$, which may represent either an order parameter or a Higgs field for a transition without symmetry breaking. Scattering from $φ$ produces non-Fermi-liquid behavior in the normal state, while the same fluctuations mediate enhanced Cooper pairing. By the Harris criterion, symmetry-preserving disorder couples most strongly to the coefficient of the $φ^2$ term, locally tuning the system toward or away from criticality. This random mass (``Harris disorder'') leads to a finite density of localized low-energy $φ$ modes even when the fermionic states remain extended and continue to provide Landau damping. We study the onset of pairing mediated by these localized overdamped bosonic modes. Starting from a real-space linearized Usadel equation for the two-particle propagator, we show that the localized bosonic wave functions generate both a spatially random pairing vertex and an effective random potential for Cooper pairs. Numerical solution, a self-consistent Born approximation, and Lifshitz-tail analysis reveal two regimes of the pairing instability. At high temperatures, pairing nucleates compact superconducting puddles on the most localized bosonic modes. At lower temperatures, an extended pairing eigenstate appears, but its transition scale and spatial structure remain strongly affected by mesoscopic correlation effects and enhanced probability of returns to favorable regions of the localized bosonic glue. The resulting distribution of local pairing scales has a power-law tail, in contrast to the stretched-exponential tails of disordered BCS superconductors. This mechanism provides a route to broad gap inhomogeneity and superconducting puddles in quantum-critical metals, and offers an interpretation of STM measurements of the cuprates.

Bipartite Lieb-Mattis ferrimagnetism and the symmetry-protected Haldane phase are paradigmatic mechanisms in quasi-one-dimensional quantum magnets. Both emerge, in distinct regimes, in a mixed-spin $1$-$\tfrac{1}{2}$-$\tfrac{1}{2}$ Heisenberg trimer chain with antiferromagnetic backbone exchange $J$ and a side spin-$\tfrac{1}{2}$ coupled to each backbone spin by an exchange $J_t$ of either sign. Using the density matrix renormalization group, we compute magnetization curves and the entanglement spectrum and entropy. For $J_t>0$ a robust ferrimagnetic plateau forms at magnetization per unit cell $m=1$, whose multiplet entropy reflects how the conserved magnetization splits between the halves. For $J_t<0$ an $m=0$ plateau opens and grows with $|J_t|$, while the $m=1$ plateau closes. As $J_t\to-\infty$ the chain maps onto a spin-$1$ Heisenberg chain with coupling $J/2$: the $m=0$ width $Δh\simeq 0.196$ matches half the Haldane gap. Exponentially localized spin-$\tfrac{1}{2}$ edge states and the even-fold degeneracy of the entanglement spectrum confirm the Haldane character of the $m=0$ phase.

Generically a quantum system initialized in an unentangled state will, under unitary dynamics, rapidly become entangled, a process closely related to information transport and to thermalization. Disorder can suppress the growth of entanglement and result in memory of initial conditions. In non-interacting systems this arises from localization of single-particle states, the occupancy of which is fixed by the initial condition. In interacting systems similar localized conserved quantities persist, but with the added feature that they are coupled, resulting in entanglement growth which is distinct from both non-interacting localized systems and from generic ergodic systems. The Fermi-Hubbard model has two degrees of freedom per site -- charge and spin -- and disorder may be present in both of these. We study the growth of entanglement in two scenarios -- disorder in charge equal and unequal to that in spin, and determine the distinct contributions of charge and spin degrees of freedom by expanding the Hamiltonian in terms of a set of optimally localized conserved quantities with separate charge and spin character. We find that coupling between charge and spin is significantly weaker than charge-charge and spin-spin coupling. While this decoupling is present in all our results, it is only apparent when the strength of the disorder in the two sectors is different such that there is a separation between the characteristic timescales of the contributions to entanglement made by charge and by spin.

We show that the non-Hermitian skin effect (NHSE) can enhance pairing correlations in moiré Hubbard systems through a channel-selective mechanism: skin-induced localization amplifies the boundary density of states, strengthening local pairing tendencies within an intermediate ``golden window'' of non-reciprocity $γ\in[0.5,1.2]\,t$. Using exact diagonalization of the non-Hermitian Hubbard model on triangular lattices with open boundaries, we map the $(U,γ)$ phase diagram. The double occupancy $D(γ)$ exhibits non-monotonic behavior -- rising by up to 21\% then declining -- reflecting a competition between NHSE-enhanced boundary pairing and over-localization. A decomposition of the pairing susceptibility $χ_{\mathrm{SC}}$ on the $3\times3$ cluster reveals that the NHSE acts \emph{channel-selectively}: it enhances on-site pairing ($+21\%$) while simultaneously suppressing competing antiferromagnetic correlations (22\% reduction), so that the total pairing susceptibility, dominated by the on-site channel, grows by $+98\%$ on that cluster. These trends are corroborated by an independent non-Hermitian DMRG calculation and establish an enhancement of finite-cluster pairing correlations rather than trivial density redistribution. We do not claim long-range superconducting order. A BCS scaling estimate converts the same pairing-response signal into a dome-shaped $T_c(γ)$ fingerprint, suggesting an experimentally distinguishable response in coherent-drive versus reservoir-dominated moiré devices.

Superconductivity mediated by Cooper quartets of charge 4e is a phenomenon of key importance for both the understanding of certain exotic superconductors and the engineering of quantum memories protected at the hardware level by topological order. Josephson junction arrays in the shape of the dice lattice constitute one of the main candidates for the realization of this phase of matter. Here, we analyze numerical signatures of the emergence of this exotic phase when superconducting dice arrays are frustrated by inserting one third of a flux quantum per rhombic plaquette. We adopt simulations of relaxation dynamics to study the critical current of such devices and analyse the Fourier decomposition of their bulk supercurrents. A further characterization of these systems at finite temperature is obtained through Monte Carlo techniques and two-dimensional infinite tensor networks. Our results indicate that the observed peaks of the critical current and temperature at frustration 1/3 correspond to a superconductor-insulator phase transition compatible with the deconfinement of half-vortices, while the low-temperature correlation functions of the model confirm the onset of a 4e superconducting phase mediated by Cooper quartets. We finally address the effects of Josephson energy and flux disorder typical of experimental arrays, and comment on the role of charging energies in the corresponding two-dimensional quantum model.

The nonequilibrium electron-phonon interplay is central to thermalization of solids, yet the microscopic picture of transient states and relaxation pathways remains incomplete. Previous nonequilibrium Green's function (NEGF) studies were restricted to local phonons and local self-energy approximations, leaving momentum-dependent dynamics largely unexplored. In this work, we demonstrate the power of the recently developed quantics-tensor-train (QTT) NEGF framework through large-scale lattice simulations with arbitrary phonon dispersions. QTTs provide a memory-efficient representation of two-time Green's functions, enabling momentum-resolved simulations with full electron-phonon feedback on lattices up to 256x256 sites. Comparing optical and acoustic phonon models, we reveal a hierarchy of relaxation bottlenecks that extends the well-known phonon-window bottleneck effect. For optical phonons, we confirm the main phonon-energy window and uncover a reduced window separating momentum-space regions of excess and deficit electronic population. We also identify a separate bottleneck in phonon thermalization, rooted in the momentum-dependent coupling to the particle-hole continuum. For acoustic phonons, the phonon-energy window acquires pronounced momentum dependence dictated by simultaneous energy-momentum conservation. The reduced window becomes asymmetric; directional scattering between Brillouin-zone regions creates a persistent bottleneck for low-momentum phonon modes. The high momentum and frequency resolution of our spectra further reveals a direct correspondence between phonon relaxation and charge response. Our results establish QTT-NEGF simulations as a scalable and controlled framework for quantitative nonequilibrium electron-phonon dynamics, overcoming previous lattice-size and propagation-time limitations and providing accurate reference data for time-resolved spectroscopies.

Altermagnetism and topological insulators represent two of the most transformative frontiers in modern condensed matter physics, spintronics, and quantum information science. Bringing these two paradigms together opens a largely unexplored route toward fundamentally new quantum phenomena. Here, we predict a topological altermagnetic phase in bilayer tetragonal Fe-based superconductors and reveal it as a highly tunable platform for valley-polarized anomalous Hall physics. Based on first-principles calculations, we show that the characteristic spin-splitting and valley polarization can be effectively tuned via applied strain. Moreover, the resulting valley-polarized anomalous Hall conductivity can be manipulated by shifting the Fermi level. These findings reveal a powerful route for controlling altermagnetism in topological materials and identify a realistic material platform for its experimental realization and technological exploitation.

In correlated insulators, the interplay among coexisting charge, spin, orbital, and lattice degrees of freedom gives rise to rich quantum phenomena, while unraveling the interplay is not straightforward. In the family of cubic $5d^2$ double perovskites, the ground spin-orbit coupled electronic states of $5d$ metal sites are degenerate and couple to the Jahn-Teller active vibrations, whereas no experimental evidence of the symmetry-lowering in the low-temperature ordered phases has been reported. To quantitatively unravel the nature of $5d^2$ centers, we apply equation-of-motion coupled cluster (EOM-CC) theory to analyze the vibronic and magnetic properties of $5d^2$ W sites of Cs$_2$WCl$_6$. We derive the electronic and vibronic model Hamiltonians, calculate the W $L_3$ edge resonant inelastic x-ray scattering (RIXS) spectra, and determine the effective magnetic moment. The simulated RIXS spectra show that vibronic coupling makes several peaks asymmetric. The effective magnetic moments exhibit a temperature dependence similar to that observed experimentally, confirming the validity of the calculated distribution of low-energy levels. Our calculations indicate that the Jahn-Teller effect in Cs$_2$WCl$_6$ is in a weak regime, and noticeable deformation would not occur, whereas the dynamic Jahn-Teller effect modulates the shapes of the RIXS spectra and affects the magnetic moment. This work demonstrates the usefulness of the EOM-CC method for predicting physical phenomena on metal sites in correlated insulating materials.

Quantum hardware platforms are getting increasingly sophisticated in their ability to simulate condensed matter, including but not limited to strongly-correlated, topological, and non-equilibrium phenomena. This review surveys recent progress in quantum-hardware-based simulations of condensed matter, primarily emphasizing gate-based digital quantum computer simulation, with analog experiments discussed as complementary benchmarks. We first review major hardware platforms, including superconducting qubits, trapped-ions, ultracold atoms, Rydberg arrays, photonic systems, and moire quantum materials. We then introduce the basic ingredients of digital quantum simulation. Building on this foundation, we discuss representative applications to condensed-matter physics, spanning ground-state problems, strongly correlated matter, topological phases, non-equilibrium dynamics, open-system physics, and high-energy-physics-inspired simulations. Finally, we summarize key methodological tools used in state-of-the-art quantum-simulation workflows. We emphasize that present noisy quantum simulations serve not only as near-term demonstrations, but also as prototypes for the encodings, diagnostic protocols and error-control strategies required for future fault-tolerant quantum simulation.

Light-induced superconducting-like states have been reported in several classes of correlated materials. From a theoretical point of view, the induction of $η$-pairing is a promising route to nonthermal superconductivity. Numerical studies of photo-doped Mott systems revealed $η$-pairing states with very high effective critical temperatures. These investigations were however restricted to particle-hole symmetric states in large-gap Mott insulators, while the experiments were performed on strongly correlated metallic systems. It is thus relevant to explore if $η$-pairing also exists in non-particle-hole symmetric setups and in photo-excited metallic states. Here we use steady-state nonequilibrium dynamical mean field theory combined with a strong-coupling impurity solver up to third order to investigate this issue. We find that in the strongly correlated regime with large Mott gap, and for low effective doublon and holon temperatures, $η$-pairing is robust against changes in the total filling and an imbalance in the doublon and holon density. An asymmetry in the effective doublon and holon temperatures can however strongly suppress the order parameter. In photo-doped metallic systems with a three-peak structure in the local density of states, $η$-pairing can be realized in set-ups with positive doublon and holon temperatures and a population inversion in the low-energy quasi-particle band.

Non-Hermitian systems host a wide range of unconventional topological phenomena while large-scale simulations in finite three dimensional systems remain challenging because of the rapidly growing number of sites. In particular, higher-order topological corner modes are often studied only in small lattices, where strong finite-size effects can mask their intrinsic behavior. Here, we develop a tensor-network framework that combines quantics tensor cross interpolation with the kernel polynomial method, enabling compact representations of large non-Hermitian tight-binding Hamiltonians and direct calculations of real-space spectral functions for systems exceeding one billion lattice sites. Using this approach, we investigate three-dimensional non-Hermitian higher-order topological insulators with with structured real-space geometries. The unprecedented system size enables direct access to the macroscopic regime and allows corner-mode spectral responses to be resolved in genuinely three-dimensional systems. By tuning the loss strength, we identify distinct in-gap corner modes across weak- and strong-loss regimes. Our results establish tensor-network algorithms as a powerful strategy to perform real-space spectral calculations in exceptionally large non-Hermitian systems.