arXiv Curator

The discovery of high-temperature superconductivity in Ruddlesden-Popper nickelates has drawn great attention. However, unlike cuprates and iron-based superconductors, Ruddlesden-Popper nickelates exhibit superconductivity either under high pressure in bulk samples or under compressive strain in thin films. Genuine bulk superconductivity under ambient pressure has remained elusive in these materials, precluding key measurements such as specific heat and superfluid density. In this work, we combine density-functional-theory, dynamical-mean-field-theory, and random-phase-approximation to solve the superconducting gap equation for bulk hole-doped bilayer nickelate La$_{3-x}$Sr$_x$Ni$_2$O$_7$ at ambient pressure. We find that hole doping induces a Ni-$d_{3z^2-r^2}$-derived $γ$ pocket on the Fermi surface, and serves as a tuning parameter for both its size and \textit{shape}. As $x$ approaches 0.4, the $γ$ pocket evolves from circular to diamond-shaped and expands to span half of the Brillouin zone, resulting in nearly perfect Fermi surface nesting with the optimal nesting vector $\textbf{Q} = (π, π)$. This, in turn, strongly enhances antiferromagnetic spin fluctuations and substantially increases the leading superconducting eigenvalue to a level at which superconductivity becomes experimentally observable. Our work provides both a robust mechanism and an experimentally feasible route to inducing the long-sought bulk superconductivity in La$_3$Ni$_2$O$_7$ without pressure or strain.

The discovery of high-$T_c$ superconductivity in Ruddlesden-Popper nickelates has sparked substantial effort towards understanding unconventional electronic states beyond a traditional cuprate-like d^9 configurational ground state. An understanding of the interplay between magnetic ground states and multi-orbital physics is key for establishing a microscopic mechanism for superconductivity. In the bilayer nickelates, spin density wave (SDW) order is a prominent feature in the non-superconducting regime. However, its relation to superconducting pairing remains an open question. Here, we use resonant x-ray scattering to examine the existence of SDW order in superconducting bilayer nickelate thin films La$_2$PrNi$_2$O$_7$ (LPNO). Comparing superconducting and oxygen-deficient LPNO thin films, we find that superconductivity occurs in SDW-free, oxygen-stoichiometric regions, whereas oxygen-deficiency promotes SDW order, indicating phase segregation of SDW and superconductivity. Furthermore, Ni-$L_3$ and O-$K$ edge spectroscopy reveals distinct electronic structures - particularly along the $c$-axis - between the two domains. Our results identify oxygen stoichiometry as a key parameter controlling interlayer coupling and thus the electronic structure of bilayer nickelates. In concert with theory, we propose that a ligand hole primarily resides at the inter-bilayer apical oxygen, forming a robust interlayer five-spin polaron state, which serves as the ground state for superconducting bilayer nickelates.

In cuprates, electron doping yields a much lower superconducting $T_c$ than hole doping. For recently discovered nickelate superconductors, the analogous doping strategies become more challenging. Consequently, while hole-doped Ruddlesden-Popper (RP) nickelates have been extensively studied, electron-doped RP nickelates remain rarely explored both experimentally and theoretically. Here we fill this gap by systematically investigating the two-orbital bilayer model for three representative systems: bulk La$_3$Ni$_2$O$_7$ at ambient pressure and 15\,GPa, and a heterostructure La$_3$Ni$_2$O$_7$:La$_3$Al$_2$O$_7$ that provides a feasible experimental route to electron doping. Using first-principle calculations and large-scale dynamical cluster quantum Monte Carlo simulations, we find that electron doping generically enhances $s^\pm$-wave pairing superconductivity (SC) in all three cases, with the heterostructure showing the highest $T_c$ in the underdoped regime. Furthermore, our results suggest an inter-orbital cooperative mechanism that the pairing on the $d_{x^2-y^2}$ orbital, induced by that on the $d_{z^2}$ orbital, plays vital role in the SC. This work provides the first theoretical prediction of enhanced SC in electron-doped RP nickelates and calls for future experimental verification.

The observation of superconductivity in undoped infinite-layer nickelates $R$NiO$_2$ ($R$ = rare earth) challenges our current understanding and calls for a re-examination of the underlying electronic structure of this family of materials. In this context, it is particularly important to extend the investigation of $R$NiO$_2$ compounds from the intensively studied hole-doped regime to the almost unexplored electron-doped one. Here, we use a combination of density-functional theory and dynamical mean-field theory to study the evolution of the electronic structure of infinite-layer nickelates in these two doping regimes. We find a striking asymmetry in the self-doping of the Ni-$d_{x^2-y^2}$ band due to the $R(5d)$ states: while this effect is strongly suppressed upon hole doping, electron doping instead leads to an increase in the size of the $R(5d)$ electron pockets, but without effectively hole-doping the Ni-$d_{x^2-y^2}$ band. This asymmetry has an important impact on the magnetic response as antiferromagnetism is rapidly suppressed upon hole doping, whereas it remains the ground state upon electron doping. Despite these differences, electronic correlations on both sides of the phase diagram are dominated by the Ni $d_{x^2-y^2}$ orbital, suggesting that a single-band description may be appropriate for infinite-layer nickelates in both the electron- and hole-doped regimes.

The synthesis of high-quality Ruddlesden-Popper (RP) nickelates remains challenging due to variations in oxygen content and the prevalence of intergrown RP phases. Precisely controlling the stoichiometry and characterizing the resulting physical properties are essential for understanding the mechanism of high-$T_c$ superconductivity in these materials. In this work, we synthesize a series of La$_3$Ni$_2$O$_{7+δ}$ samples with systematically controlled oxygen content and perform comprehensive structural and compositional analyses. Precise oxygen tuning enables us to tailor the microstructure, yielding a pure bilayer phase, a mixture of bilayer and hybrid single-layer-bilayer phases, and a predominantly bilayer phase containing trilayer intergrowths. High-pressure transport measurements reveal distinct superconducting transitions with contrasting $T_c$ values, corresponding to the bilayer phase, the hybrid phase, and trilayer inclusions. Notably, we find that oxygen content not only governs the phase purity$-$i.e., the presence of intergrowth phases$-$but also directly modulates the upper critical field ($H_{c2}$) of the bilayer superconductivity. By establishing a phase diagram of $T_c$ and $H_{c2}$ as functions of oxygen content in La$_3$Ni$_2$O$_{7+δ}$, this work advances synthetic control and provides new insights into the superconducting mechanism of RP nickelates.

The extent to which electronic dimensionality influences magnetism and superconductivity in Ruddlesden-Popper (RP) nickelates remains unsettled. Here we report high-precision crystallographic-axis-resolved dc transport measurements on high-quality single crystals of bilayer and trilayer RP nickelates. Using a six-terminal geometry, we self-consistently determine the intrinsic in-plane ($ρ_\parallel$) and out-of-plane ($ρ_\bot$) resistivities on the same crystal, while minimizing uncertainties associated with current redistribution in highly anisotropic conductors. We uncover strong intrinsic electronic anisotropy in both bilayer and trilayer nickelates, in contrast to the weak anisotropy inferred from conventional four-probe measurements. Moreover, $ρ_\bot$ exhibits a nonmonotonic temperature dependence, revealing a universal coherent-to-incoherent crossover in interlayer transport. Across the RP nickelate series, the maximum superconducting transition temperature ($T_c$) observed under pressure is inversely correlated with the ambient-pressure resistivity anisotropy, suggesting that stronger interlayer electronic coherence is favorable for superconductivity. In addition, $ρ_\bot$ serves as an exceptionally sensitive and selective probe of magnetic and density-wave orders, exhibiting pronounced anomalies, whereas only weak signatures are observed in $ρ_\parallel$. Our results highlight interlayer coherence as a key organizing parameter that both tracks the relevant magnetic correlations and is closely tied to superconductivity, providing stringent constraints on microscopic theories of high-$T_c$ superconductivity in nickelates.

For treating correlated electronic systems on quantum computers, we propose a quantum-classical hybrid scheme for dynamical mean-field theory (DMFT). In the quantum part of the scheme, we use modified quantum phase estimation (QPE) circuits suitable for the one-particle Green's function (GF) at a finite temperature so that we can extract spectral amplitudes and the excitation energies without knowing the excitation channel invoked at each measurement. In the classical part of the scheme, we adopt an approach that estimates reasonably the GF based on the data collected from the QPE sampling. We dub the approach the QPE averaged over variable grids (QAVG), that may help one to reconstruct the GF via optimization of trial parameters and modeling the probability distributions for various settings of the QPE circuits. We apply the QAVG-DMFT scheme to SrVO$_3$ to demonstrate its validity via numerical simulations.

RuO$_2$ has been regarded as a prototypical candidate for metallic altermagnet, offering a potential platform for high-speed and high-efficiency spintronics. However, the magnetic ground state of RuO$_2$ remains a topic of active debate due to conflicting experimental reports. In this work, we investigate the effect of electron correlations in RuO$_2$ using density functional theory combined with dynamical mean-field theory (DFT+DMFT). In contrast to previous DFT-based studies, DFT+DMFT captures essential dynamical correlation effects, yielding spectral functions and optical conductivities in excellent quantitative agreement with experiments, and further reveals that RuO$_2$ resides in the close vicinity of both the paramagnetic-altermagnetic phase boundary and the itinerant-localized crossover, rendering the magnetic ground state highly susceptible to external perturbations. Indeed, even a minimal compressive strain of $\sim$0.5% is sufficient to drive the system into an altermagnetic phase. These findings elucidate the origin of the conflicting experimental observations and reveal that dynamical correlation effects are the key driving force behind the highly tunable magnetic ground state of RuO$_2$.

Superconductivity in infinite-layer nickelates has drawn wide interest as a cuprate analogue, yet how the electronic structure evolves with hole doping remains unsettled. Here we map the doping- and temperature-dependent unoccupied states of the La-based infinite-layer nickelate La1-xCaxNiO2 using O K-edge and Ni L-edge x-ray absorption spectroscopy. Superconductivity occurs for 0.18<=x<=0.27. Near x~0.20-0.23, low-energy spectral weight redistributes: Ni3d-dominated states decrease while O2p-hybridized states increase, indicating an orbital-selective crossover in Ni-O covalency. This crossover coincides with a sign reversal of the Hall coefficient and precedes the reduction of the superconducting critical temperature at higher doping. By directly linking transport anomalies and the superconducting dome to a measurable Ni-O orbital reorganization, our results provide a key step toward a unified, orbital-resolved phase diagram for infinite-layer nickelates and a practical route to engineer superconductivity via hybridization control.

Recent angle-resolved photoemission measurements on La$_3$Ni$_2$O$_7$ have challenged the density-functional-theory-based picture of three Fermi surfaces by revealing that the $d_{z^2}$-derived $γ$ band can reside below the Fermi level. Motivated by this discrepancy, we investigate a realistic bilayer two-orbital Hubbard model using time-dependent variational principle (TDVP)-based cluster perturbation theory (CPT), alongside large-scale density matrix renormalization group (DMRG) calculations. Our TDVP-CPT calculations, performed on clusters of up to 16 physical sites, reveal that electronic correlations drive a pronounced orbital-selective reconstruction of the low-energy spectrum: the $d_{z^2}$ spectral weight is progressively depleted, the $γ$ band sinks below the Fermi level, and pseudogaps open on the remaining $α$ and $β$ bands, leaving Fermi arcs dominated by the $d_{x^2-y^2}$ orbital at strong coupling. Furthermore, large-scale DMRG calculations demonstrate that the leading superconducting correlations evolve consistently with this Fermi surface reconstruction, transitioning from $d_{z^2}$-dominated to $d_{x^2-y^2}$-dominated interlayer spin-singlet pairing while retaining an $s_{\pm}$ structure. Consequently, our results indicate that the disappearance of the $γ$ pocket is not detrimental to superconductivity; rather, it signals a correlation-driven shift of the pairing channel mediated by interlayer antiferromagnetism, Hund's coupling, and inter-orbital hybridization.

We investigate the emergence of nontrivial topology in a twisted cuprate bilayer described by the Hubbard model in the weak-interaction regime. Our results show that the topological character depends sensitively on the doping level. For $U/t=3.85$, the Chern number assumes a value of $\pm 8$ in the electron-doped case, whereas it vanishes (0) in the hole-doped regime. The presence of nontrivial topology is further supported by an analysis the associated edge states and their chirality in a finite-width geometry, while keeping full correlations.

The discovery of high critical-temperature $T_{\mathrm{c}}$ superconductivity in La$_3$Ni$_2$O$_7$ under high pressure has led to a rapid expansion of the $T_{\mathrm{c}}$ range through lanthanide $Ln$ substitution, and to ambient-pressure superconductivity in strained thin films, yet the exploration of new bilayer nickelates remains strongly constrained by thermodynamic stability. Beyond the difficulty of synthesis of bulk single-crystals, here we report on the pressure-induced high-$T_{\mathrm{c}}$ superconductivity in epitaxially-stabilized Pr$_3$Ni$_2$O$_7$ thin films. While the Pr$_3$Ni$_2$O$_7$ films exhibit insulating behaviour at ambient pressure regardless of ozone-annealing treatment, they show $T$-linear metallic transport and superconductivity reaching an onset $T_{\mathrm{c}}$ of 66 K and zero-resistance at nearly 40 K at 22 GPa. Furthermore, Nd$_3$Ni$_2$O$_7$, with the smaller rare-earth ion Nd, can also be stabilized, however, superconductivity is not observed in the measured pressure range. Epitaxial stabilization enables us to examine the dependence of $T_{\mathrm{c}}$ and the critical pressure $P_{\mathrm{c}}$ for superconductivity on the $Ln$ ion in $Ln_3$Ni$_2$O$_7$ ($Ln$ = La, Pr, Nd). These results suggest that a higher $P_{\mathrm{c}}$ is required for smaller $Ln$ ions, consistent with trends observed in bulk studies of $Ln$ substitution. This study demonstrates that epitaxial stabilization is a powerful technique to further expand the family of superconducting bilayer nickelates.

In cuprates, electron doping yields a much lower superconducting $T_c$ than hole doping. For recently discovered nickelate superconductors, the analogous doping strategies become more challenging. Consequently, while hole-doped Ruddlesden-Popper (RP) nickelates have been extensively studied, electron-doped RP nickelates remain rarely explored both experimentally and theoretically. Here we fill this gap by systematically investigating the two-orbital bilayer model for three representative systems: bulk La$_3$Ni$_2$O$_7$ at ambient pressure and 15\,GPa, and a heterostructure La$_3$Ni$_2$O$_7$:La$_3$Al$_2$O$_7$ that provides a feasible experimental route to electron doping. Using first-principle calculations and large-scale dynamical cluster quantum Monte Carlo simulations, we find that electron doping generically enhances $s^\pm$-wave pairing superconductivity (SC) in all three cases, with the heterostructure showing the highest $T_c$ in the underdoped regime. Furthermore, our results suggest an inter-orbital cooperative mechanism that the pairing on the $d_{x^2-y^2}$ orbital, induced by that on the $d_{z^2}$ orbital, plays a vital role in the SC. This work provides the theoretical prediction of enhanced SC in electron-doped RP nickelates and calls for future experimental verification.

Bilayer Ruddlesden-Popper nickelates display superconductivity near 80 K under high pressure, establishing a new nickelate platform for studying unconventional high-temperature superconductivity. The recent stabilization of superconducting RA3Ni2O7 (RA = rare earth or alkaline earth) films at ambient pressure has changed the experimental landscape: epitaxial strain can reproduce key structural ingredients of the high-pressure phase while making transport, spectroscopy, microscopy, and device-oriented measurements directly accessible. This Review summarizes the experimental progress on ambient-pressure superconducting bilayer nickelate films, with emphasis on synthesis routes, oxygen stoichiometry, substrate-induced strain, normal-state transport, superconducting properties, doping phase diagrams, and momentum-resolved electronic structure. We highlight several issues that remain unsettled, including the reproducibility of phase-pure ultrathin films, the microscopic origin of the two-step superconducting transition, the role of oxygen defects and substrate-derived doping, the position of the Ni 3dz2-derived γ band, and the pairing symmetry. We close by outlining experimental directions that could establish a more quantitative link among crystal structure, orbital reconstruction, and superconductivity in bilayer nickelate films.

The discovery of high transition temperature (high-Tc) superconductivity in Ruddlesden-Popper (R-P) bilayer nickelates under high pressure has stimulated extensive work to understand the underlying mechanism and search for superconductors with higher Tc. The recent realization of superconductivity in R-P bilayer nickelate thin films with onset Tc above 40 K at ambient-pressure enables the use of a wide array of powerful experimental tools to investigate the unconventional high-Tc superconductivity in bilayer nickelates. Here, using ultra-low temperature scanning tunneling microscopy/spectroscopy (STM/S) and electrical transport study, we report the first successful observation of an energy-symmetric, flat-bottom U-shaped gap with zero residual density of states around the Fermi level in the high-Tc nickelate (La,Pr)3Ni2O7 thin film grown on SrLaAlO4 substrate. Before and after STM/S studies, transport measurements on the same sample reveal consistent superconducting behaviors showing zero resistance, with an onset Tc above 40 K and zero resistance Tc above 20 K. The tunneling spectra exhibit highly unconventional temperature evolution, characterized by a rapid filling of the U-shaped energy gap to a V-shaped gap as the temperature increases. Furthermore, the U-shaped energy gap is reduced under a c-axis magnetic field of 14 T. The energy-symmetric U-shaped gap, taken together with its dependence on magnetic field and temperature, is consistent with the behavior of a superconducting gap, suggesting a nodeless gap function at ultra-low temperatures. Our findings shed new lights on the nature of high-Tc superconductivity and provide an encouraging and thought-provoking hint for a local superconductivity with Tc above liquid nitrogen boiling temperature in nickelate superconductors at ambient or zero pressure.

Ruddlesden-Popper bilayer nickelates provide an emerging platform for studying high-temperature superconductivity, yet the superconducting pairing symmetry remains under debate. Here, we use atomic-resolution scanning tunnelling microscopy and spectroscopy to investigate superconducting 1.5-unit-cell (La,Pr)3Ni2O7 films grown on SrLaAlO4. A cryogenic ultrahigh-vacuum (UHV) sample transfer preserves an ordered sqrt(2) * sqrt(2) surface and yields reproducible U-shaped spectra with two gap scales of ~14 and ~20 meV and extended flat zero-conductance bottoms. By contrast, samples exposed for a longer time in UHV without cooling during transfer show V-shaped spectra despite retaining the surface reconstruction and a transport superconducting transition onset above 40 K. Wide-energy-range spectra indicate that oxygen loss can mix density-wave-related spectral weight. Our measurements provide an atomic-scale observation of the intrinsic nodeless superconducting gap in bilayer nickelate ultrathin films.

Although spin fluctuations are central to the physics of high-Tc cuprates, their relevance to strange-metal behavior in the overdoped regime remains unclear. Here, we use high magnetic fields to suppress superconductivity and an NMR protocol tailored to electronic inhomogeneity to show that the low-energy limit of the dynamical spin susceptibility χ''(q,omega) at x=0.25 in La2-xSrxCuO4 increases continuously down to our lowest temperatures. This behavior is suggestive of quantum-critical fluctuations, a leading candidate mechanism for strange-metal transport, yet is observed well beyond the spin-stripe critical doping x=0.19. Our data further reveal that the spin dynamics are spatially inhomogeneous, suggesting that nanoscale electronic inhomogeneity may underlie this apparent paradox. These observations provide new insight into the electronic state from which strange-metal behavior emerges.

Metallic systems with delocalized valence electrons host collective charge density oscillations known as plasmons. On the other hand, conventional insulators do not have free electrons and the low energy charge degrees of freedom are pinned to the ions. The fate of the collective charge excitations in the intermediate regime is an outstanding question. This problem is especially important for strongly correlated systems such as the layered cuprates, where unconventional superconductivity and other emergent phenomena arise from valence electrons on the border between Mott localization and itinerancy. Using resonant inelastic X-ray scattering, we track this evolution in the prototypical electron-doped cuprate Nd2-xCexCuO4. We find a continuous transformation of the low-energy charge response: from an acoustic plasmon in the metallic regime, to a gapped hybrid mode at intermediate doping, and finally to a nearly dispersionless 139 meV excitation at half filling. Remarkably, the 139 meV excitation has approximately twice the energy of the oxygen breathing phonon responsible for the dispersion kink observed in angle-resolved photoemission spectroscopy, and is consistent with a putative 2-phonon excitation observed in Raman spectroscopy. These results establish a unified picture of collective charge excitations across the phase diagram of electron-doped cuprates, showing that such modes persist across the Mott transition via strong coupling to lattice degrees of freedom and revealing a missing link in the charge dynamics of carrier doped Mott insulators.

Triangular lattice antiferromagnets first attracted attention as a frustrated magnetic lattice which can serve as a platform to realize the resonating valence bond state. While the triangular lattice itself was shown to support classical 120 degree order, many theoretical phase diagrams suggest a quantum spin liquid state within a small range of parameters. One possible avenue to achieve such a state is to tune the anisotropy of the triangular lattice antiferromagnet by applying uniaxial pressure. This motivated our Raman scattering study of quasi-two-dimensional antiferromagnet $α$-SrCr$_2$O$_4$ under applied uniaxial pressure. Under ambient conditions, $α$-SrCr$_2$O$_4$ develops long-range helical magnetic order below T$_N$ = 43 K. We identify two-magnon excitations associated with this long-range antiferromagnetic order below T$_N$ at 15.5 meV and 40 meV by comparison with spin wave calculations. We observe the two features from the two-magnon excitation shift away from (towards) each other under applied tensile (compressive) pressure, indicating a decrease (increase) in anisotropy. Raman active phonons show a shift to higher (lower) frequencies under applied compressive (tensile) pressure, indicating efficient transmission of pressure and tuning of the lattice. We show spin wave and two-magnon density of states calculations under uniaxial pressure are consistent with our experimental results.

We investigate the electronic phases of an effective Hubbard model on the body-centered-cubic lattice, motivated by alkali-doped fulleride molecular solids. The model incorporates renormalized on-site interactions and an effective inverted Hund's coupling originating from electron-phonon interactions. To access complementary interaction regimes, we employ two theoretical approaches. In the intermediate-coupling regime, the on-site repulsive interaction is approximated by a long-range interaction in momentum space, yielding an exactly solvable Hatsugai-Kohmoto model supplemented by a BCS-type pairing term. Within this framework, we analyze the superconducting instability and demonstrate a first-order normal-superconducting phase transition, characterized by a discontinuous jump of the order parameter. In the strong-coupling regime, where pairing fluctuations are suppressed, we apply the spin rotationally invariant slave-boson formalism to map out the temperature-interaction phase diagram. This analysis reveals first-order transitions between a Fermi-liquid phase, an antiferromagnetic phase, and a Mott insulating phase, with a narrow intermediate region where all three phases compete. The resulting phase diagram captures the interplay of itinerancy, magnetic order, and Mott localization in three dimensions and provides a unified perspective on superconducting and correlation-driven phenomena in fulleride-inspired lattice systems.

Despite decades of intense theoretical and experimental investigation, the two-dimensional Fermi-Hubbard model still resists a complete microscopic understanding. Conventional approaches typically probe global observables and locally resolved correlation functions. Here, we develop a complementary perspective based on the measurement of entanglement. Using both an ultracold-atom quantum simulator and numerical simulations based on the dynamical vertex approximation, we find that entanglement is closely tied to the onset of the enigmatic pseudogap regime: spin-singlet entanglement emerges only as the pseudogap sets in and, in contrast to classical correlations, remains confined to nearest-neighbour sites in this regime. Our results, therefore, disfavour purely classical-fluctuation theories of the pseudogap and constrain microscopic models to those that develop nearest-neighbour spin-singlet entanglement at the pseudogap onset.

We study the normal-state, doping-driven phase diagram of the square-lattice Hubbard model using the dynamical cluster approximation combined with the numerical renormalization group as a cluster solver, which gives direct access to real-frequency dynamics at essentially zero temperature. In a parameter regime relevant for cuprates, $U=7t$ and $t'=-0.3t$, we find a critical doping $p^{\ast}$ that marks a continuous quantum phase transition between a pseudogap metal and a normal Fermi liquid. The transition is identified by a continuous collapse, from both sides, of the Fermi-liquid scale extracted from charge, spin, and $d_{x^2-y^2}$-wave pairing susceptibilities. This collapse produces a non-Fermi-liquid regime at intermediate energy scales, which appears to extend to arbitrarily low scales at $p^{\ast}$. As $p^{\ast}$ is crossed from the normal Fermi liquid at $p>p^{\ast}$ into the pseudogap metal at $p

To date, controlling the steady-state electronic band structure in high-Tc cuprate superconductors has been achieved primarily through chemical doping or magnetic fields. Here, we present that ultrafast optical excitation can instead drive the electronic band structure of Bi2Sr2CaCu2O8+δ into a photostationary, long-lived excited state. At sufficiently high pump fluences, this state undergoes a Lifshitz transition of the Fermi surface, characterized by a change in topology from hole-like to electron-like. Time- and angle-resolved photoemission spectroscopy, supported by single-band tight-binding calculations, reveals that 1.6 eV photoexcitation induces band-structure evolutions closely analogous to those produced by chemical doping. These results point to an efficient photodoping mechanism involving cooperative effects, including charge transfer, renormalization of effective electronic correlations, and defect-assisted charge trapping. Our findings raise fundamental questions regarding thermalization processes occurring on timescales comparable to the laser repetition period in cuprates. More broadly, ultrafast optical control enables access to otherwise inaccessible regions of the phase diagram by tuning the pump fluence.

Binding energy, which quantifies pair formation, is a key factor in the emergence of superconductivity. Here, we show that even when multiple spins are complexly coupled, hole-doped systems, which can be mapped onto the universal hardcore boson model in the strong-binding-energy limit, exhibit promising signatures of superconductivity. We analytically and numerically demonstrate this concept in the double Kondo lattice model. Using the density-matrix renormalization group method, we show that a pairing state is maintained via a crossover even for parameters away from the strong-coupling regime. Additionally, we find that once binding energies are sufficiently generated, pair correlations develop similarly regardless of the details of local spin correlations. Our findings suggest useful guidelines for research on superconductivity.

The exact ground state of the two-dimensional Hubbard model is critical for understanding cuprate superconductivity. Previous numerical studies on narrow cylinders found insulating, static stripes that inherently suppress superconductivity. Here, using constrained-path auxiliary-field quantum Monte Carlo on isotropic lattices up to $24 \times 24$ sites, we show static stripes are boundary artifacts. The true 2D thermodynamic limit yields a topologically melted fluid of dynamically fluctuating magnetic pockets. Furthermore, we reveal the microscopic real-space origin of cuprate particle-hole asymmetry. Hole doping actively melts the magnetic background, driving a Lifshitz transition that unleashes macroscopic $d_{x^2-y^2}$ phase coherence at an optimal $x \approx 0.150-0.200$. Conversely, electron doping preserves rigid antiferromagnetic domains, confining carriers to narrow fault lines that optimally saturate early at $x \approx 0.100$. By extracting the macroscopic off-diagonal long-range order across both regimes, we perfectly recover the skewed phenomenological superconducting dome. Our parameter-free theoretical curve aligns with empirical Uemura and Božović scaling relations, capturing the underdoped emergence, distinct optimal peaks, the 1/8 anomaly suppression, and overdoped collapse. These results prove that robust d-wave superconductivity is the intrinsic ground state of the pure Hubbard Hamiltonian.

We investigate the connection between doping-induced in-gap states and underlying spin excitations in Mott insulators by employing cluster perturbation theory on one-dimensional (1D) and quasi-1D Kitaev-Hubbard models. By manipulating Kitaev-like hopping terms ($t^{\prime}$) that selectively control spin anisotropies in the strong-coupling limit, we establish a direct correspondence between the kinetic dispersion of the in-gap states and the spin excitation spectra. Specifically, in the Z chain, in-gap states evolve from a gapless dispersion to a gapped flat band as the system transitions from the Heisenberg to the Ising model, exhibiting a gap scaling of $2t^{\prime 2}/U$ that matches the Ising spin gap. In the XY chain, the in-gap states split into a dispersive and a flat branch at the Kitaev limit, perfectly mirroring the Jordan-Wigner fermionic spectrum. For the two-leg ladder, we observe an emergent broad continuum of in-gap states that reflects the fractionalization of spin excitations, accompanied by a gap manifesting the presence of topological $Z_2$ visons. Our results establish a robust correspondence between charge and spin dynamics in doped Mott insulators and demonstrate that in-gap states can serve as a probe of exotic quantum spin phenomena, including fractionalization and topological excitations, offering a new pathway to investigate spin liquids via spectroscopic probes of charge excitations.

The microscopic mechanism by which doped holes disrupt the antiferromagnetic order is one of the fundamental questions in cuprates. In this work, we propose a unified microscopic theory in which doped holes form spatially localized Zhang-Rice singlets which actively mediate emergent spin exchange. Rather than acting as simple non-magnetic vacancies, these localized states introduce emergent next-nearest $J_2$ and third-nearest $J_3$ neighbor superexchanges. This dopant-induced exchange pathway generates significant magnetic frustration, naturally explaining the rapid collapse of the Néel AFM order and the emergence of a spin-glass phase on the hole-doped side. Our findings provide a comprehensive framework for understanding the complex doping-driven magnetic phase transitions and magnetic electron-hole asymmetry in lightly doped cuprates.

We investigate the strong correlation effect in the spinless electron-hole two-band Hubbard model using the dynamical mean-field theory. At half filling, both the renormalization factor $Z$ and the number of conduction electrons (valence holes) $n_c$ decrease with increasing the interband Coulomb interaction $U$ down to $Z\sim 0.15$ and $n_c\sim 0.1$ for $U_c \sim \mbox{bandwidth}$ at which the first-order Lifshitz transition occurs from a correlated semimetal with a large effective mass $m^*/m=Z^{-1}$ to a band insulator with a finite gap due to the Hartree shift. A slight hole doping $x$ in the band insulator with a large $U>U_c$ yields a remarkable correlated semimetal with $Z\sim 0.1$ at $x\sim 0.01$, where in-gap states emerge above the valence band top and those weights increase with increasing $x$ together with the increase in $Z$ similar to the in-gap states observed in doped Mott insulators.

Triangular-lattice systems attract a lot of attention due to various frustration-induced and strongly correlated effects. Here, we focus on the charge-ordering phenomenon by means of investigation of the extended Hubbard model with dynamical mean-field theory (DMFT). By considering the intersite nearest-neighbor interaction we have found a very rich phase diagram that contains large number of features, phases, and phase transitions. Among them are pinball-liquid (PL) phases where we distinguish charge-transfer-driven and Mott-localization-driven PLs; phase transitions that change their order as model parameters evolve (from discontinuous to continuous); very strong particle-hole asymmetry. Various features of the phase diagram are found to be better understood by means of the simple mean-field approximation (MFA). Moreover, besides helping with interpretation of the phase diagram, the MFA results together with results for the atomic-limit model are found to be able to set rather good expectations on how the DMFT phase diagram should look like. Nevertheless, a few features were not expected and are found within the DMFT, such as a small-region intermediate metallic phase on an electron-doped side of the phase diagram.

In cuprate superconductors, electron-electron repulsion results in characteristic spectroscopic features known as `waterfalls', where the sharp quasiparticle dispersion transitions into broad Hubbard bands. However, in multi-orbital systems, the additional Hund coupling results in behavior that defies the conventional Mott--Hubbard paradigm, creating qualitatively distinct `superdispersive' features in the spectral function. Here, we use tunneling spectroscopy to reveal this signature of Hund physics in Sr$_2$RuO$_4$. By combining density functional theory, dynamical mean-field theory, and continuum local density of states calculations, we show that the experimental features are in excellent agreement with theoretical predictions and intimately linked to the non-monotonous energy dependence of the real part of the self-energy in a Hund metal. Our results provide direct experimental evidence for Hund-induced spectroscopic features and open a new route to probing correlation effects in quantum materials.

We derive a discrete spectral representation of the single-particle self-energy using a discrete evaluation of Kugler's symmetric improved estimator. Our construction can be used on both the real and the complex (Matsubara) frequency axis. It is guaranteed to remain causal at the numerical level, in contrast to standard approaches that may generate unphysical negative spectral weight or require additional broadening. Our representation can be used for any Hamiltonian; here we apply it to quantum impurity models and in dynamical mean-field theory. The latter is formulated with a discrete hybridization function throughout its self-consistency loop. In both cases and across various numerical methods, we obtain significantly improved accuracy for a range of impurity properties.

We study dynamical scaling in the quantum-critical fan of the pseudogap-metal to Fermi-liquid transition of the two-dimensional Hubbard model. Using a four-patch dynamical cluster approximation with the numerical renormalization group as a cluster impurity solver, we access real-frequency dynamics over several decades at arbitrary temperatures. Close to the critical doping, the local spin and cluster-current susceptibility spectra exhibit $x=ω/T$ scaling of the form $χ''(ω,T)\sim \tanh(x/2)$, and the cluster contribution to the optical conductivity obeys $Tσ'_{\mathrm{cl}}(ω,T) \sim \tanh(x/2)/x$, implying a $1/T$ cluster dc conductivity. In the scaling regime, the vertex contribution to the cluster optical response is much larger than the bubble contribution. We further find evidence for a marginal-Fermi-liquid nodal self-energy. This, together with the $1/T$ vertex contribution to the conductivity, implies strange-metal optical transport in the quantum critical region. Our results describe several qualitative aspects of several experimental observations.

A Kondo insulator (KI) is a prototypical example of a highly entangled phase of matter, where many-body interactions between local moments and delocalized electrons engender the non-magnetic insulating ground state. Conventionally, the local moments arise from atomic multiplet states with a narrow bandwidth, limiting Kondo coherence to low temperatures. Here, we realize a new paradigm for constructing the KI state with hybridized molecular orbitals in FeSb2. Resonant inelastic X-ray scattering (RIXS) at the Fe L-edge reveals distinct signatures of band-like continuum states and localized states. Comparisons with first-principles calculations establish a mixed-configuration ground state with hybridized Fe d-Sb p molecular orbitals as basis states. By systematically investigating the RIXS momentum, temperature, and doping dependences, we find propagating collective modes commensurate with many-body charge and spin excitations. Our results pave the way for understanding the emerging class of unconventional d electron insulators and engineering high temperature Kondo many-body states.

Altermagnetism is a recently identified class of magnetic order characterized by unconventional momentum-dependent spin splitting in the absence of net magnetization, and understanding its electronic and magnetic properties is essential for revealing its fundamental physics and potential applications. In this work we investigate two intercalated transition-metal dichalcogenides, Fe$_{1/4}$NbS$_2$ and V$_{1/3}$NbS$_2$, as candidate altermagnetic materials by using effective tight-binding and spin models complemented by first-principles calculations. We show that the $g$-wave electronic spin splitting originates from bond-dependent hopping anisotropy, leading to material-dependent nodal structures. For the magnetic excitations, the emergence of chiral splitting in the magnon dispersion is controlled by single-ion anisotropy, which manifests as altermagnetic-like nodal structures when spins are oriented along an easy-axis. Conversely, this altermagnetic signature disappears when the spins are aligned in an easy-plane. Beyond linear spin-wave theory, we find that $1/S$ corrections from magnon--magnon interactions preserve the symmetry and nodal structure of the band splitting while generally reducing its magnitude, with strong antiferromagnetic exchange leading to a non-negligible renormalization of the chiral splitting. Our findings establish intercalated transition-metal dichalcogenides as promising platforms for understanding the interplay between crystal symmetry, non-relativistic spin splitting, and magnetic properties in altermagnets.

When the Fermi surface and the Brillouin zone boundary touch at a Van Hove point, mean-field analysis predicts a ferromagnetic (Stoner) instability at finite $T_{MF}$ for any coupling strength due to the divergent density of states. However, the predicted effect has not been observed experimentally. Several qualitative theoretical proposals have been put forward to explain why the mean-field prediction fails. Based on numerically exact results for the two-dimensional Hubbard model with an ordinary Van Hove singularity, we uncover the mechanisms behind the suppression of the ferromagnetic instability. We employ two diagrammatic Monte Carlo approaches: (i) the four-channel self-consistent approximation and (ii) numerically exact method of combinatorial summation of diagrams with controlled resummation of the truncated expansion. We find that the system avoids the Stoner instability down to temperatures an order of magnitude below $T_{MF}$ due to the combination of the downward renormalization of the effective coupling and the suppression of the density of states by the loss of the quasiparticle residue.

A recent article [Phys. Rev. X 15, 011047 (2025)] utilizes group-equivariant convolutional neural networks to study the ground state of the kagome Heisenberg antiferromagnet. On the largest finite-size cluster studied to date ($N=108$), the authors report variational energies significantly lower than other numerical methods, including state-of-the-art density matrix renormalization group (DMRG) calculations. In contrast to previous results suggesting a possible spin-liquid ground state, the authors observe a spinon pair-density-wave ground state. We find that: (i) the reported low energies are artifacts of broken ergodicity in the Metropolis--Hastings sampling, since the single-spin-flip update rule utilized by the authors effectively freezes the Markov chains; and (ii) when ergodic sampling is enforced via spin-exchange updates, the neural network converges to energies significantly higher than existing DMRG results, calling the paper's claims into question.

The quest for efficient devices has fueled research in thermoelectric materials. In these materials, the goal is to maximize the Figure of Merit $ZT$. One of the components of this quantity is the Seebeck coefficient, which measures the voltage generated in response to a temperature gradient. Recent studies have revealed that strong electronic correlations can enhance the Seebeck coefficient, leading to anomalous behavior near half-filling. However, the impact of interactions beyond the on-site Hubbard remains mostly unexplored. In this work, we investigate the Seebeck coefficient considering attractive interactions, nearest-neighbor interactions, sublattice potentials and electron-phonon coupling. We find that additional interaction scales can enhance the Seebeck coefficient, while also leading to multiple anomalous changes of sign as a function of doping. We also show that the anomalous behavior is connected to a gap opening in the ground state. Moreover, electron-phonon coupling also lead to a Seebeck anomaly, even without on-site repulsion. We connect these changes of sign in the Seebeck coefficient with a restructuring of the Fermi surface and a change in its topology, an effect commonly seen in cuprates.

Numerous unconventional superconductors such as cuprates, heavy-fermions, and twisted-bilayer graphene exhibit incoherent metallic transport above the superconducting critical temperature. This phenomenon cannot be described with Fermi-liquid theory and has presented a significant theoretical challenge to overcome. We utilize the two-dimensional Yukawa-SYK model of fermions with spatially random coupling to quantum-critical bosons to study transport in a manner which is non-perturbative in the coupling strength. Our work provides a microscopic model of quantum-critical incoherent metals and their concomitant properties, including a non-Boltzmann transport formula between resistivity and quasi-particle lifetime, violation of the Mott-Ioffe-Regel resistivity bound, and violation of the Kovtun-Son-Starinets shear viscosity to entropy density bound.

We find a mapping between the attractive Fermi-Hubbard model and the repulsive Bose-Hubbard model at finite temperature and at imaginary chemical potential $μ=iθ$. We show, by using a large $N$-expansion, that the partition functions of the two models are related by a simple shift $θ\to θ+ π$. This condition maps the BCS--BEC crossover of attractive fermions to a Bose--Fermi crossover (fermion-like occupation) of repulsive bosons. Central feature of this correspondence plays the thermal kernel $g(βE,φ),$ whose analytic continuation $g_B(βE,φ) = g_F(βE,φ+π)$ governs the bosonic and fermionic sectors. Interestingly, we are able to find that the special angles $φ= 2π/3,4π/3$ for fermions correspond to $φ= π/3,5π/3$ for bosons, marking the boundaries of a universal thermal window. We further argue that the present mechanism shows that fermionization can occur at finite interaction strength through a thermodynamic effect induced by the imaginary chemical potential. This suggests that it is a new way of fermionization (not a change in statistics but a fermion-like behaviour) unlike the Tonks--Girardeau limit, where fermionization arises from an infinite repulsive interaction and anyonic or Floquet-engineered systems where transmutation emerges from modified statistics or dynamics. Essentially, the phase $φ$ is a statistical parameter; by twisting the thermal phase, it generates fermion-like behaviour without hard-core constraints or infinite repulsion but only by using thermodynamics. We derive the gap equation and number equation for the bosonic model, highlighting the role of the imaginary chemical potential as a statistical regulator. Our results provide a unified framework for understanding crossovers in interacting lattice systems.

Noncollinear (NC) magnetism and spin-orbit coupling (SOC) are indispensable for predictive ab initio materials simulations with pronounced relativistic effects and magnetic frustration, yet they significantly increase the cost of cubic-scaling density functional theory (DFT) by introducing complex 2-component wavefunctions per electron and consequently much larger eigenproblems. We present a GPU-centric high-performance framework for NC-SOC DFT that combines: (i) algorithmic advances for solving finite-element (FE) discretized DFT equations; (ii) residual-based Chebyshev filtered subspace iteration (R-ChFSI), tolerant to inexact matrix-vector products, for the resulting sparse generalized eigenproblem; (iii) a matrix-free strategy for accelerating FE Poisson solver; (iv) R-ChFSI-enabled mixed-precision computation with block floating-point compressed MPI communication at compression ratios over 4x, preserving double-precision robustness while reducing compute and data movement costs; and (v) a communication efficient band-partitioning algorithm to improve scalability. Numerical results demonstrate improved time-to-solution and excellent scaling on exascale architectures, enabling fully relativistic pseudopotential DFT simulations of up to 100,000 electrons.

We present a scalable quantum simulation framework for real-time dynamics of the multi-flavor Gross-Neveu model in 1+1 dimensions. Using superconducting quantum processors at utility scale, we develop a hardware-efficient Trotterization whose per-step circuit depth scales with fermion flavor number rather than total system size, enabling simulations beyond 100 qubits. A central contribution of this work is the Localized Diagonal Operator Approximation (LDOA), which systematically reduces the overhead associated with quartic interactions. We formulate diagonal unitary synthesis as a structured least-squares problem in phase space and obtain analytic solutions via the Moore-Penrose pseudoinverse. This formulation provides a principled and quantitatively controlled approximation: in the small Trotter-step regime, the unitary error is directly linked to the phase reconstruction error and vanishes asymptotically as the Trotter step size decreases. This establishes a clear mathematical foundation for the LDOA while significantly reducing two-qubit gate counts and circuit depth, and is broadly applicable to diagonal quantum operators with long-range structure, making it particularly well suited for quantum hardware with limited qubit connectivity. Using these techniques, we run large-scale simulations on IBM superconducting processors and study real-time observables, including density-density correlators. We benchmark against exact diagonalization and tensor network-based methods, finding strong agreement across system sizes. These results show that combining hardware-aware circuit design with rigorous approximations enables practical near-term simulation of interacting fermionic field theories and provides a scalable pathway toward more complex quantum field theory simulations.

We develop an optimized continuous-field quantum Monte Carlo (QMC) algorithm to investigate the SO(5) nonlinear sigma model with a Wess-Zumino-Witten term, which describes half-filled Dirac fermions in 2+1 space-time dimensions akin to graphene and Yukawa coupled to a quintuplet of compatible mass terms. To regularize the theory, we project onto the lowest Landau level for both spherical and torus geometries. Our algorithm reduces the computational complexity to $O(βN_{\mathbf{q}} N_φ^2)$, yielding a speedup of a factor of $N_φ$ (the number of magnetic fluxes, i.e., system size) relative to prior works [1-3]. This advance enables us to simulate system sizes up to $N_φ=140$ on torus and $N_φ=49$ on sphere, far exceeding the maximum sizes accessed, and to map out the universal phase diagram of the model on both geometries. Most notably, we identify and characterize a critical point that separates an SO(5)-broken ordered phase at small coupling from an SO(5)-symmetric disordered phase at large coupling. The critical point becomes multicritical upon the inclusion of terms that break the SO(5) symmetry down to $\mathrm{U}(1) \times \mathrm{SU}(2)$, relevant for the deconfined phase transition between Néel antiferromagnetic and valence-bond-solid orders in quantum magnets. While the precise nature of the disordered phase in the thermodynamic limit remains to be determined, we argue that it is neither conformal nor trivially gapped, akin to a chiral quantum spin liquid with a small gap. Our finding of a multicritical point in the phase diagram of the SO(5) nonlinear sigma model with Wess-Zumino-Witten term resolves the long-standing open question of its global structure, and our QMC framework opens a new avenue for systematic studies of projected Hamiltonians, ranging from correlated flat bands to fractional quantum (anomalous) Hall systems.

We develop a representability-aware and interpolable neural network (NN) framework for predicting two-particle reduced density matrices (2-RDMs). The NN incorporates a subset of representability conditions through its architecture and loss function, and can operate on different momentum meshes, enabling evaluating the representability conditions across multiple meshes, which we call interpolated representability condition. The framework can be used either to predict 2-RDMs on large momentum meshes by interpolating exact results from small meshes, or as a variational 2-RDM ansatz optimized by energy minimization on arbitrary meshes. We apply this approach to the fractional Chern insulator in the one-band projected model of twisted bilayer MoTe$_2$ at twist angle $3.89^\circ$ and hole filling $2/3$. Trained on exact-diagonalization (ED) 2-RDMs from meshes with $12$ or $18$ momentum points using six different NN architectures, the best NN is the residual multilayer perceptron, which predicts the $6\times6$ 2-RDM with $97.07\%-98.18\%$ accuracy relative to the ED 2-RDM but predicts an energy $77.353$ meV above ED ground-state energy. We then variationally optimize the NN on several meshes including $6\times6$, predicting a $6\times 6$ energy of just $0.104$ meV below ED while maintaining $98.94\%-98.96\%$ accuracy. Compared with the conventional boundary-point semidefinite programming, which gives an energy $5.560$ meV below ED with $96.40\%-98.94\%$ accuracy, the NN achieves a more accurate energy and similar accuracy while using only less than 1/20 as many parameters. Eventually, we add a symmetric mesh of $48$ momentum points to the variational optimization of the NN, and provide a prediction of the many-body ground-state energy and the many-body quantum metric on that mesh.

Parafermions with non-Abelian statistics have been proposed as a promising platform for quantum computation, potentially enabling a broader set of topologically protected gates than Majorana fermions. The experimental and theoretical exploration of these exotic quasiparticles remains challenging, as their stability is linked to strong electron-electron interactions. A key step toward practical applications is the controlled shuttling of parafermionic modes, which is required for implementing geometric braiding operations. In the present work, we investigate the real-time dynamics of the elementary shuttling process by applying a combination of the density matrix renormalization group and the time-dependent variational principle approaches. We analyze the transport of $\mathbb{Z}_4$ parafermion edge states and assess the corresponding adiabatic speed limit under experimentally relevant conditions.

Recently, magnetic 2-dimensional (2D) van der Waals (vdW) materials have garnered tremendous attention. The vdW ferromagnet Fe5Ge1Te2 has a Curie temperature Tc of ~ 270 K, which is tailorable by tuning the stoichiometry and the Fe deficiency to reach room temperature. To explore the expanded compositional space, we implemented combinatorial synthesis and high-throughput characterization to investigate the structural phase distribution and ferromagnetism of a Fe-Ge-Te thin film library. The library was prepared by magnetron co-sputtering followed by annealing in vacuum or in an inert environment. Composition and structural phase distribution of the 177 pads in the library were characterized using high-throughput wavelength dispersive spectroscopy (WDS), X-ray diffraction (XRD), and two-point probe resistance measurements. We leverage unsupervised machine learning to cluster the XRD dataset into groups of compositions with similar structural phases, and further study the ferromagnetic properties via SQUID magnetometry and X-ray magnetic circular dichroism (XMCD) across different clusters. The results are compared against magnetization and structural models calculated using DFT. Our results demonstrate that the hexagonal crystal structure is a critical prerequisite for ferromagnetism in this system, and that unexplored materials adopting this structure can be efficiently identified as possible ferromagnetic materials using our high-throughput, ML-assisted framework. This workflow based on the combinatorial strategy allows us to rapidly capture the composition-structure-magnetic property map across a broad compositional landscape of novel magnetic materials.

Designing novel inorganic materials through generative models remains an important challenge for material science, driven by the complexity and diversity of inorganic structures across expansive chemical compositions and structural landscape. The vast combinatorial space of inorganic compounds demands innovative, AI-driven approaches to overcome limitations in generative accuracy and efficiency. To address this, we introduce a novel method that redefines the encoding and generation of inorganic materials by utilizing domain-specific symmetry-aware representation. Our approach not only refines the representation of intricate inorganic structures but also contributes to the field of material discovery by enhancing the precision and stability of generated candidates. Central to our methodology is a novel padding technique that exploits crystal symmetry information to enhance the encoding process. By integrating Wyckoff position length-aware padding into an encoder architecture, we achieve a more robust informed representation of inorganic materials. This symmetry-driven enhancement improves deep learning models to generate stable, previously unexplored inorganic structures with superior accuracy and computational efficiency. Furthermore, we introduce an end-to-end system that leverages the machine learning potential models to seamlessly generate novel, even those unseen in the training data, and stable inorganic materials from initial data to validated output. This pipeline integrates advanced generative models with stability analysis, marking a significant leap forward in the automated exploration and design of next-generation inorganic materials. Our method improved reconstruction accuracy 5.3% in proton conductor data, and generated 63.5% more novel stable inorganic material to baseline model on the perov-5 dataset.

The discovery of high-${T_c}$ superconductivity in pressurized bilayer nickelate La$_3$Ni$_2$O$_7$ presents a new arena for exploring unconventional pairing mechanisms. A pivotal yet unresolved issue is the specific role of the $3d_{z^{2}}$ orbital of Ni. While its inter-layer super-exchange antiferromagnetic coupling is widely considered crucial for superconductivity, the role of its itinerancy remains undetermined. Early studies showed that the superconductivity is accompanied by the emergence of a small Fermi pocket of the $3d_{z^{2}}$ orbitals. However, recent experiments show controversial results on the role of the $3d_{z^{2}}$ Fermi pocket on superconductivity. Motivated by these experimental results, we investigate an effective bilayer two-orbital model for La$_3$Ni$_2$O$_7$ using density-matrix renormalization group (DMRG) on a minimal one-dimensional geometry. By systematically varying the $3d_{z^{2}}$ orbital filling from $1/12$ doping to half-filling, we observe a pronounced suppression of superconducting correlations near half-filling. Our results demonstrate the itinerancy of $3d_{z^{2}}$ orbital is favorable for the pairing in the bilayer two-orbital model for La$_3$Ni$_2$O$_7$. Moreover, we observe that the pairing correlation is enhanced in regions where charge fluctuations are large, suggesting a competition between charge order and superconductivity in the model.

The recent discovery of high-temperature superconductivity in pressurized and thin film nickelates has generated intense interest, yet the nature of magnetism in their ambient-pressure parent phases remains poorly understood, despite its potentially crucial role in pairing. Here we use neutron scattering to resolve the spin order and dynamics of single-crystalline La3Ni2O7, an ambient-pressure parent of this class. Well defined spin excitations are observed at Q = (0, 0.5, 2.5), featuring a~5 meV spin gap and anisotropic in-plane dispersions, with zone-boundary softening along the transverse direction indicative of competing exchange interactions. The excitations exhibit pronounced out-of-plane modulations with bilayer periodicity, providing direct evidence for antiferromagnetic interlayer coupling. Their dispersion is well described by a bilayer Heisenberg Hamiltonian with strong interlayer exchange and competing in-plane couplings within a stripe-type magnetic order. Normalization of the spectra to absolute units reveals that, although the spin-wave bandwidth is only about 25% of that in cuprates, the local dynamic susceptibility at comparable energies is significantly enhanced, yielding a total fluctuating moment of comparable magnitude. These results highlight intense mid-energy spin excitations rooted in substantial electronic correlations as a defining feature of this family, establishing a magnetic framework distinct from cuprates and directly relevant to understanding superconductivity in this system.

Understanding how quasiparticles are generated following a rapid quench of superconductivity in high-Tc cuprates is a key unresolved problem in nonequilibrium superconductivity. Here we resolve these processes in optimally doped YBCO [YBa2Cu3Oy(y=6.94, Tc=92 K)] using broadband (0.16 -4.1 eV, ca. 100 fs) and nearly single-cycle (6 fs) transient reflectivity spectroscopy. We show that within a few femtosecond, enhanced electron-electron Umklapp scattering dominates, signaling a transient modulation of long-range Coulomb interactions on the eV scale. This regime is followed by a rapid suppression of the scattering rate of the mid-infrared absorption associated with carriers dressed by spin fluctuations. We attribute this observation to an ultrafast decoupling of quasiparticles from the spin-fluctuation background, occurring on a 90 fs timescale set by the inverse optical gap. These findings reveal the correlated many-body dynamics underlying quasiparticle generation in cuprates and provide further clues for unconventional pairing mechanism.

We theoretically investigate a freestanding membrane of infinite-layer nickelate Nd$_{0.85}$Sr$_{0.15}$NiO$_2$ under pressure by constructing a seven-orbital effective model based on first-principles calculations. By performing the fluctuation exchange (FLEX) approximation, we demonstrate that the seven-orbital model explains a monotonic increase in $T_c$ reported in a recent experiment. This enhancement of superconductivity is attributed to the mitigation of excessively strong electron correlations caused by exceptionally low valence of Ni atom. Furthermore, we examine the dynamical stability of the crystal structure under pressure through phonon calculation.

We theoretically investigate a freestanding membrane of infinite-layer nickelate Nd$_{0.85}$Sr$_{0.15}$NiO$_2$ under pressure by constructing a seven-orbital effective model based on first-principles calculations. By performing the fluctuation exchange (FLEX) approximation, we demonstrate that the seven-orbital model explains a monotonic increase in $T_c$ reported in a recent experiment. This enhancement of superconductivity is attributed to the mitigation of excessively strong electron correlations caused by exceptionally low valence of Ni atom. Furthermore, we examine the dynamical stability of the crystal structure under pressure through phonon calculation.

We propose a three-dimensional multi-orbital tight-binding model for rare-earth nickelates RNiO$_3$ that treats charge, spin, orbital, and lattice degrees of freedom on equal footing. All model parameters, including the on-site interactions $U$ and $J$ and the electron-phonon (el-ph) coupling to the breathing mode, are extracted from hybrid-functional DFT calculations for the small-bandwidth nickelate LuNiO$_3$. The model describes three competing insulating phases governed by the interplay of $U{-}3J$ and el-ph coupling to the breathing and Jahn--Teller (JT) modes. For large $U{-}3J$, the insulating state is stabilized by local JT distortions on high-spin Ni$^{3+}$ sites. For smaller $U{-}3J$, the system undergoes charge disproportionation, $2\mathrm{Ni}^{3+}\rightarrow\mathrm{Ni}^{2+}+\mathrm{Ni}^{4+}$, resulting in the spin-polarized charge-ordered state observed experimentally below the Néel temperature in small-bandwidth RNiO$_3$. When the JT energy on the Ni$^{2+}$ site exceeds Hund's exchange $3J$, a distinct charge- and orbital-ordered insulating phase emerges in which the two $e_g$-electrons occupy the same orbital with opposite spin. The stability of this phase is further confirmed by self-consistent calculations within the full three-dimensional tight-binding model. This newly predicted metastable state, characterized by JT distortions in a nonmagnetic charge-ordered RNiO$_3$ phase, shows that the onset of magnetic order is not required for the metal-insulator transition in RNiO$_3$.

Geometric frustration in quantum systems can stabilize unconventional phases of matter that avoid traditional magnetic ordering at low temperatures. Here, we observe this phenomenon while mapping out the finite temperature phase diagram of the spin-1/2 Heisenberg antiferromagnet on the ruby lattice with next-nearest-neighbor interactions. Using an infinite tensor network state (iTNS) optimized and measured with belief propagation (BP) and corrections to BP, we observe the low temperature formation of stable magnetic plateaus at various magnetic field strengths. We find these plateaus host a novel `simplex liquid state' -- a disordered phase involving strongly paired spin simplices that retains non-zero residual entropy due to an exponentially large subspace of crystalline configurations. We accurately quantify the energy gap associated with these states and show that, as the temperature of the system is lowered, it does not go through a phase transition to reach them: the heat capacity remains finite and continuous at all observed temperatures. Our work demonstrates how BP-based tensor network techniques provide a powerful route to understanding frustrated quantum magnets at finite temperature.

We study ground-state properties and phase diagram of the $J_{1}$-$J_{2}$ antiferromagnetic XY model on the honeycomb lattice by means of the developed corner transfer matrix renormalization group algorithm with the two-site unit cell and the infinite spiral projected entangled pair states ansatz. We identify the main phases: Néel, Ising, collinear, and incommensurate spiral phases, as well as the transitions between them, as functions of the ratio $J_{2}/J_{1}$. In the regime of competing types of ordering, we show that the energies of the dimerized states are systematically higher than the energies in the collinear phase. This collinear phase transforms to the incommensurate spiral phase through the second-order phase transition upon a further increase of $J_2/J_1$.

We report numerical ground states for the dipolar XY spin model, which describes extended antiferromagnetic interactions in two-dimensional arrays of polar molecules and two-level Rydberg atoms. Carrying out large-scale density matrix renormalization group (DMRG) calculations, we compute ground state properties on nine of the eleven Archimedean lattices--tilings of the plane by regular polygons. Four of these host trivial paramagnets, while another four develop collinear Neel magnetic order, as was found previously for the square lattice. For the ordered states, we calculate the hydrodynamic parameters (magnetization, susceptibility, and stiffness) and compare to linear spin wave theory. We also investigate the triangular lattice, for which we find several competing phases including coplanar magnetism, stripe density wave order, and a possible spin liquid; their relative stability is sensitive to the long-range couplings present in our dipolar model. Finally, the Archimedean classification is completed by the kagome lattice, which we argue in a companion work is likely to be a Dirac spin liquid.

SrCu2(BO3)2 (SCBO) is a paradigmatic realization of the Shastry-Sutherland model, hosting geometrically frustrated spin dimers and a variety of quantum magnetic phases and phenomena. Although its magnetic properties have been extensively studied, the high-energy electronic excitations that determine the crystal-field environment and Cu-O hybridization have remained largely unexplored. Here we combine Cu L3-edge resonant inelastic x-ray scattering (RIXS), broadband optical spectroscopy, and electronic-structure calculations to determine the relevant local and interband excitation energy scales in SCBO. RIXS resolves a well-defined manifold of localized Cu2+ d-d excitations between 1.8 and 2.4 eV, whose energies and polarization dependence are well reproduced by multireference quantum-chemistry calculations. In contrast, optical spectroscopy identifies charge-transfer excitations with an absorption onset near 1.2-1.6 eV and a broader higher-energy structure around 4.5 eV, which are qualitatively captured by DFT+U calculations. Taken together, these results define the characteristic energy scales of d-d and CT excitations, offering quantitative benchmarks for computational frameworks and providing essential input for refining superexchange-based magnetic models of this prototypical frustrated quantum antiferromagnet.

We propose a family of layered quantum spin-orbital models as a platform to study fractionalization, unconventional forms of symmetry-breaking order, and their possible coexistence. The models are built by stacking $N$ layers of a square-lattice system in which Kitaev-type interactions promote the formation of a $\mathbb{Z}_2$ quantum spin-orbital liquid and coupling the different layers via Ising spin interactions. Using a parton construction, we show how, at low energies, these Hamiltonians can be mapped to $N$-component Fermi Hubbard models on a $π$-flux square lattice at half filling. We also demonstrate that the models acquire an emergent SU($N$) symmetry in the limit of equal all-to-all interlayer couplings and argue that, for $N>2$, the proximity to this limit offers the potential to realize an array of competing phases. To illustrate this point, we compute the zero-temperature phase diagram of the effective $N=3$ Hubbard model within mean-field theory and uncover rich phenomena, including intertwined orders and flavor-selective localization. Mapping back to the original degrees of freedom reveals that the ground states realize distinct forms of magnetic fragmentation, wherein the orbitals remain in a quantum liquid state whereas the spins can present conventional long-range order or nonlocal order characterized by a nontrivial string order parameters. We highlight possible extensions of our construction as well as its potential to provide concrete microscopic models for different fractionalized quantum critical points.

While the interplay of correlations and geometric frustration in doped Mott insulators provides a fertile ground for exotic quantum phases, the nature of the metallic state emerging upon particle doping remains poorly understood. In this work, we investigate the triangular-lattice Hubbard model with particle doping and provide compelling evidence for an intrinsic, interaction-driven quantum state, which we term the Nagaoka supermetal. This state is characterized by a sublinear temperature dependence in the DC resistivity, along with singular behaviors in the charge compressibility and zero-frequency spectral weight. To understand the origin of these singular properties, we derive an effective low-energy model and demonstrate that a higher-order Van Hove singularity emerges from the reconstructed dispersion. This singularity gives rise to a power-law divergence in the density of states, capturing the anomalous properties observed in the supermetallic regime. Our findings offer a new perspective on non-Fermi liquid states in geometrically frustrated systems and are directly accessible in current ultracold atom experiments.

Strange metallicity, characterized by a linear temperature dependence of the resistivity, is observed in a broad range of correlated materials, including heavy-fermion compounds and cuprate superconductors. It has also recently been reported for the Kagome metal Ni$_3$In, where almost localized and itinerant electronic degrees of freedom coexist as a result of a partially flat band. We investigate the correlated electronic structure and transport properties of Ni$_3$In with dynamical mean field theory (DMFT) calculations performed on a minimal single-band Hubbard model, constructed from compact molecular orbitals. Despite the large band filling, even for moderate Hubbard repulsion, we observe a non-Fermi-liquid like frequency dependence of the self-energy, as well as the formation of local magnetic moments. With increased hole doping, a crossover to a heavy Fermi-liquid regime is found. We interpret these results in terms of an effective model for the partially filled narrow band near $k_z=0$.

Since the discovery of high-temperature superconductivity in nickelate superconductors, it is an open question how closely the superconducting state resembles that of cuprate superconductors. One salient feature of the phase diagram of the high-temperature cuprate superconductors is stripe order. Despite their prevalence, real-space imaging has been limited to the charge sector. Here we use spin-polarised scanning tunnelling microscopy to visualize the local magnetic and charge distribution emerging due to a stripe order in the trilayer nickelate La$_4$Ni$_3$O$_{10}$. The stripe order exhibits a four unit cell periodicity, closely resembling that seen in cuprates, and opens a near-complete $\sim66\mathrm{meV}$ gap at the Fermi level. Crucially, discrete phase slips can be triggered by tunneling electrons above a $\sim 20\mathrm{meV}$ threshold, allowing imaging of stripe dynamics at the atomic scale. These results highlight the importance of correlation physics driving stripe-like orders in lanthanum nickelates with striking similarities to the cuprates.

Confined ultracold atoms in optical lattices provide a versatile platform for simulating lattice models of strongly correlated quantum systems, where pairing phenomena and superfluid phases can be explored under controlled conditions. While the crossover between the Bardeen-Cooper-Schrieffer (BCS) phase and the Bose-Einstein condensation (BEC) is well understood in homogeneous systems, spatial confinement breaks translational symmetry and reshapes correlation patterns, making the BCS-BEC identification in trapped geometries challenging and allowing unconventional phases to emerge with no direct analog in homogeneous systems. Here we present a characterization of the BCS-BEC crossover in harmonically confined one-dimensional Fermi-Hubbard chains. Our analysis combines Density Matrix Renormalization Group (DMRG) simulations and entanglement-based diagnostics with effective models describing the formation of tightly bound fermion pairs. This combined approach enables a detailed understanding of how the interplay between interactions and confinement reshapes the crossover, leading to insulating regions coexisting with persistent superfluid correlations. Within this framework, we further introduce conditioned correlation functions whose power-law decay allows a clear distinction between BCS-like and BEC-like regimes. The consistency between the effective descriptions and the numerical DMRG results yields a unified picture of the crossover in harmonically confined geometries.

Moiré transition metal dichalcogenides have served as a versatile platform for simulating Hubbard physics. Recent experiments have identified robust superconductivity in moiré bilayer WSe$_2$ for certain twist angles. Here, we propose the gossamer nature of the superconductivity recently discovered at half-filling and zero displacement field in twisted WSe$_2$. By mapping the moiré continuum system to an effective extended single-orbital Hubbard model on the triangular lattice, we employ renormalized mean-field theory to investigate the strong-coupling phase diagram. We find that a moderate Coulomb repulsion partially suppresses charge fluctuations while preserving a finite density of mobile doublons and holes. In this regime, the interplay between extended kinetic hoppings and antiferromagnetic superexchange stabilizes a chiral $d+id$ superconducting phase. Our results naturally account for the twist-angle-dependent evolution from a Mott insulator to a superconductor and eventually to a correlated metal. Furthermore, the model demonstrates that this half-filled pairing state vanishes rapidly upon density doping, consistent with experimental observations.

Developing non-perturbative methods to reveal exotic properties of strongly correlated fermionic systems remains one of the most essential tasks of theoretical physics. Tensor network methods with Grassmann algebra offer powerful numerical tools for fermionic many-body systems in the coherent-state path-integral representation. Despite their vast potential for both condensed-matter and particle-physics communities, Grassmann tensor network methods are somewhat underexploited in practical simulations. In this work, we provide a detailed, self-contained introduction to Grassmann tensor network methods, from the basics of the Grassmann tensor operations to the Grassmannization of typical tensor network algorithms. Furthermore, the resulting Grassmann tensor network methods are validated in several interesting models in both particle physics and condensed matter physics.

Although the density functional theory plus Hubbard $U$ correction method (DFT+U) is broadly used to study electronic structure of strongly correlated materials, the extension of this method to electron-phonon $g$ matrices has received limited attention. Here, we implement an algorithm that integrates DFT+U method with the finite-displacement method for the calculations of phonons and electron-phonon $g$ matrices. The Hubbard $U$ corrections are applied not only to electronic and phonon structures, but, more importantly, also to electron-phonon $g$ matrices. We demonstrate our algorithm in two prototypical correlated materials: infinite-layer nickelates LaNiO$_2$ and ruthenium dioxide RuO$_2$. We find that: i) While the Hubbard $U$ corrections weakly increase the electron-phonon interaction of 20% hole-doped LaNiO$_2$, its total electron-phonon coupling remains small and is insufficient to account for the observed superconducting transition temperature of about 10-30 K. Our results contrast with the recent work showing that the full GW corrections yield an elevated electron-phonon coupling of 20% hole-doped LaNiO$_2$ five times larger than its DFT value. We attribute this discrepancy to the differences in the Fermi surface topology between DFT+$U$ and GW methods. ii) The inclusion of Hubbard $U$ corrections eliminates the imaginary phonon modes of RuO$_2$ under strain on the TiO$_2$ substrate and substantially reduces the electron-phonon coupling. Our results alleviate the discrepancy between the reported large theoretical electron-phonon coupling and the low superconducting transition temperature observed experimentally. Our work provides an algorithm that fully includes the Hubbard $U$ corrections on electron-phonon properties of correlated materials, and highlights the importance of Fermi surface shape and correlation effects on phonon spectrum and electron-phonon $g$ matrices.

Conventional approaches to simulating quantum many-body dynamics produce a single trajectory: if the Hamiltonian or the initial state is changed, the computation must be re-performed. Recent efforts toward foundation models have begun to address this limitation, yet existing methods transfer across either Hamiltonians or initial states, but not both. In this work, we introduce the Universal Neural Propagator (UNP), a single, unified model that learns the functional mapping from driving protocols to time-evolution propagators. Trained in an entirely self-supervised way, a single UNP model predicts dynamics across a function space of driving protocols and an exponentially large Hilbert space of initial states simultaneously. We benchmark on a two-dimensional driven Ising model and demonstrate the UNP's accuracy and transferability across product and entangled initial states, as well as for both in- and out-of-distribution driving protocols. The UNP remains accurate at system sizes beyond exact diagonalization, and can be efficiently fine-tuned across all initial states using observable data. By shifting the object of learning from quantum states to operators, this work opens a route toward transferable simulation of driven quantum matter.

The honeycomb cobaltate Na$_2$Co$_2$TeO$_6$ has recently been proposed as a candidate material for hosting field-induced quantum spin liquid (QSL) behavior. Here, we present a comprehensive thermodynamic study of its low-temperature, high-field phase diagram using magnetization, specific heat, and magnetocaloric-effect measurements down to 1 K. In zero field, we observe a weak residual moment that provides further insight into the nature of the magnetic ground state. For in-plane magnetic fields ($B \parallel a^*$), we identify three field-induced transitions at $B_{c1} \simeq 6$ T, $B_{c2} \simeq 7.8$ T, and $B_{c3} \simeq 10.4$ T. The magnetic Grüneisen parameter and specific heat reveal clear thermodynamic signatures of these successive phase transitions enclosing two intermediate phases. Contrary to expectations for a field-induced QSL, the phase between $B_{c2}$ and $B_{c3}$ lacks enhanced magnetic entropy but instead shows behavior consistent with a distinct ordered state. Above $B_{c3}$, the absence of additional anomalies indicates a crossover to a conventional spin-polarized regime. Our results place stringent thermodynamic constraints on the proposed QSL scenario in Na$_2$Co$_2$TeO$_6$, calling for further microscopic investigations to establish the precise nature of the field-induced phases.

We study the emergence of superconductivity in the one-band Hubbard model using the spin-rotation-invariant Kotliar-Ruckenstein slave-boson (SB) approach. Motivated by its intrinsically renormalized mean-field ground state, we construct an effective pairing vertex from dynamical fluctuations about the saddle point. Solving the anisotropic, frequency-dependent gap equation on the square lattice, we map the pairing instabilities across doping, interaction, temperature and real-frequency gap structure that qualitatively match experimental cuprate observations. This framework merges strong-correlation SB-type renormalizations with RPA-type pairing transparency, providing a scalable route to modeling multi-orbital superconductivity at strong coupling.

Although quantum computing offers a promising solution for strongly correlated system simulation, existing algorithms face significant bottlenecks on current noisy intermediate-scale quantum (NISQ) devices. Here, we introduce QiankunNet-QSCI, a hybrid quantum-classical framework that addresses this challenge by combining efficient quantum-sampling with a transformer neural network. An efficient unitary selected configuration Interaction (USCI) ansatz especially designed for quantum sampling is proposed to identify the most chemically significant electronic configurations on the Zuchongzhi 3.1 quantum processor. Subsequently, the transformer model QiankunNet learns from these sparse yet critical quantum data to infer and reconstruct the complete electronic wavefunction with high fidelity. Simulation of the challenging 40-qubit [2Fe-2S] ferredoxin active center achieves chemical accuracy. Simulation of the nitrogenase P-cluster in a 114-electron 73-orbital active space also reaches 12 milli-Hartree-level agreement with the best density matrix renormalization group (DMRG) result. QiankunNet-QSCI thus offers a practical route to accurate quantum-assisted electronic structure calculations on current devices.

We develop a self-consistent spectral quadrature (sc-SQ) framework for the calculation of many-body Green functions. The method approximates the Källén--Lehmann spectral measure by Gauss--Christoffel (GC) quadrature, yielding a rational Green function representation with guaranteed spectral positivity that exactly reproduces the first $2N$ spectral moments at pole order $N$. A key component is an SVD-based rank-selection criterion on the Hankel matrix, which identifies the numerically resolvable pole rank $N^*$ from the singular-value gap and acts as a precision-guided diagnostic of correlation complexity. The scheme is made self-consistent by requiring that the spectral function used to evaluate expectation values coincides with the spectral function generated by the quadrature reconstruction. This defines a fixed-point hierarchy that connects systematically to established approximations, including Hartree--Fock and Hubbard-I, and incorporates non-perturbative features such as multi-peak spectral structure. We benchmark the approach for the Anderson impurity model against numerical renormalization group (NRG) results and apply it within dynamical mean-field theory for the Hubbard model on the Bethe lattice. The method captures the three-peak Anderson impurity spectrum and the suppression of quasiparticle weight in the half-filled Hubbard model on the Bethe lattice, including Mott-gap formation on the insulating branch for $N\geqslant 5$, in qualitative agreement with NRG references.

We performed high-precision magnetization measurements up to 120~T on three compositions of the newly discovered kagome antiferromagnet YCu$_3$(OD)$_{7-x}$Br$_{2+x}$ (YCOB), revealing a previously unobserved 5/9 fractional magnetization plateau. All YCOB samples with different Br$^-$ concentrations exhibit nearly identical magnetization curves below 60~T, whereas the 5/9 plateau appears at markedly different fields in the ultrahigh-field regime. By modeling the experimental data using tensor-network calculations, we derive the effective spin Hamiltonians for the YCOB family with three spatially anisotropic Heisenberg couplings (the 3$J$-type model), which quantitatively reproduces the measured magnetization processes and captures the composition-dependent evolution of the 5/9 plateau. Furthermore, our theoretical analysis suggests the emergence of a spin supersolid phase in the field window between the 1/3 and 5/9 plateaus, which is sensitive to spin exchange parameters and accounts for the significant variation in the critical fields of the 5/9 plateau observed among different YCOB compositions.

Motivated by the growing interest in spin liquids and topological phases, as well as the rise of deformation engineering, we study the combined effects of deformation and magnetic fields on the honeycomb Kitaev model. The Kitaev model, as one of the prototypical and exactly solvable spin liquid-hosting models, serves as a simple platform that demonstrates the rich physics one can expect at the intersection of deformation physics and quantum spin liquids. Our work builds on a simplified solution to the undeformed base model that we present. This simplified solution allows for a straightforward extension of our analysis to the deformed case. After incorporating periodic deformations into the Kitaev model (chosen for its similarity to moiré physics), we investigate the effects of a hexagonally symmetric deformation on the band structure. We find that deformation leads to a smaller Brillouin zone with new band gaps at the edges, indicating the potential for topological transitions. Finally, we introduce a magnetic field to break time-reversal symmetry and thereby allow for non-trivial topology. We find that, under specific parameter conditions, the magnetic field leads to multiple band-gap closings and openings. An investigation into topological properties reveals nontrivial Chern numbers and a plethora of topological transitions. Our results suggest possible thermal Hall or Nernst-type responses. We also suggest a potential bulk measurement approach for he Chern numbers and possible path to physical realization. Most importantly, our results serve as a demonstration of the rich phenomenology that can arise due to the interplay between deformation and spin-liquid physics.

Altermagnets (AMs) have recently emerged as a distinct magnetic class bridging central features of ferromagnets (FMs) and antiferromagnets (AFMs), offering new opportunities for spin-based electronics. While they possess zero net magnetization like collinear AFMs, they simultaneously exhibit momentum-dependent spin splitting long thought exclusive to FMs. Despite intense theoretical interest, experimentally accessible materials hosting both altermagnetism and nontrivial band topology remain scarce. EuAgAs, crystallizing in space group $P6_3/mmc$, was previously identified via density functional theory (DFT) as a bulk altermagnetic Dirac semimetal. Contrary to these predictions, our neutron diffraction experiments reveal that the bulk ground state adopts a $\mathbf{q} = (0,0,\tfrac{1}{2})$ AFM structure with an in-plane $\uparrow\uparrow\downarrow\downarrow$ spin sequence. Systematic DFT calculations, however, uncover a remarkable near-degeneracy among competing magnetic orders: the FM and AM configurations lie only $0.11$ and $0.40~\text{meV/f.u.}$ above the AFM ground state, respectively. We further show that while a simple Heisenberg model favors a spin-spiral ground state, the inclusion of non-Heisenberg biquadratic coupling stabilizes the observed commensurate AFM phase. This near-degeneracy renders the magnetic state highly tunable, with DFT predicting a transition to the altermagnetic phase under hydrostatic pressure at approximately $14 \text{ GPa}$, establishing EuAgAs as a controllable platform for accessing topological altermagnetism.

Orbital order describes a quantum state where occupied orbitals line up in a periodic pattern. While orbital physics plays a fundamental and universal role in strongly correlated electron systems, the existence and particularly the band structure fingerprint of orbital order remain a long-standing mystery. Here, we report the discovery of rare earth 5d-orbital order developed by the surface states of intermetallic compound Tb2CoAl4Ge2. Angle-resolved photoemission spectroscopy reveals characteristic nematic features like Fermi surface deformation and band split. These experimental observations can be described by a ferro-orbital order term in the mean-field Hamiltonian. The structural and magnetic origin of such order is excluded by systematic high-resolution neutron powder diffraction and scanning tunnelling microscopy measurements. Our results provide strong evidence for a pure surface orbital order scenario avoiding complications from structural distortion as in colossal magnetoresistance manganites, magnetic order as in iron-based superconductors, and charge transfer p-orbital order in cuprates.

Since its original prediction 25 years ago, room-temperature out-of-plane ferroelectricity in compressively strained SrTiO3 remains an ongoing pursuit. In this work, we investigate the structural, electrical and electromechanical properties of highly strained epitaxial SrTiO3 capacitors with rare earth nickelate electrodes. The SrTiO3 layers experience compressive strains up to -3% and exhibit pronounced tetragonality, comparable to that of bulk PbTiO3. Variable-temperature electrical measurements and room-temperature piezoresponse force microscopy reveal butterfly-shaped capacitance-voltage hysteresis and domain-like electromechanical response typical of ferroelectric materials. However, the overall behavior is inconsistent with a stable ferroelectric state. We therefore propose an alternative mechanism for the observed polarity in our samples based on spatially inhomogeneous internal fields. Our first-principles calculations show that such fields may arise from charge discontinuities between the formally charged NdNiO3 layers and charge-neutral SrTiO3 layers.

We here explore the emergence of orbital and spin Hall effects, originating beyond the $L$-$S$ coupling, and investigate the interplay between inter-orbit hybridization, relativistic Rashba spin-orbit coupling (SOC), and non-relativistic SOC, namely altermagnetic (AM) order, in a two-dimensional model Hamiltonian. The orbital (spin) Hall responses are remarkably found to be quantized within a window of Fermi energy when the strength of AM order (Rashba SOC) exceeds (falls below) the scale set by $sp$-hybridization. Importantly, orbital and spin Hall quantizations are independent of Rashba SOC and AM order, respectively, while the uniform profiles of finite orbital and vanishingly small spin moments of bands around the Fermi energy. The microscopic origin of such quantization comes from the Fermi surface-activated orbital and spin Berry curvatures. The extent of the quantized regime is strongly controlled by the intra-orbital coupling strength. As the temperature increases, the quantization is significantly compromised in the spin Hall case. We extend our analysis to the orbital and spin Nernst coefficients where the pronounced dip-peak structures signal the existence of the quantization leading to experimental relevance.

We demonstrate that altermagnetism provides a field-free mechanism for stabilizing finite-momentum superconductivity in two dimensions. Using a non-perturbative static path approximation Monte Carlo approach, we show that a d-wave altermagnet supports a robust pair-density-wave (PDW) phase that persists over a finite temperature window despite strong thermal fluctuations. The underlying mechanism originates from momentum-dependent spin splitting, which effectively enhances pairing instabilities at finite center-of-mass momentum without Zeeman fields. We identify distinct thermal scales associated with phase coherence, gap closing, and pseudogap formation, and establish characteristic spectroscopic and real-space signatures of the PDW state. Our results reveal altermagnetism as a robust route to thermally stable finite-momentum superconductivity and provide experimentally testable signatures for altermagnetic materials.

The recent discovery of altermagnetism has sparked growing interest in compensated magnetic systems as promising platforms for highly scalable spintronics. Altermagnetism is a distinct magnetic order where opposite spin sublattices are connected by rotation, yielding zero net magnetization but momentum-dependent spin splitting. To date, experimental verification of altermagnetic order has been achieved predominantly through bulk-sensitive techniques, including spin-dependent electronic spectra and transport responses. However, direct atomic-scale evidence that explicitly correlates crystal symmetry, local structural distortions, and magnetic ordering has remained unexplored. Here, we report the direct atomic-scale observation of coexisting polar distortions and altermagnetic order in MnTe, combining atomic resolution scanning transmission electron microscopy (STEM) imaging with electron magnetic chiral dichroism (EMCD) measurements. We reveal that MnTe is not an ideal uniform P63/mmc g-wave altermagnet at the atomic scale. Instead, it hosts ubiquitous inversion-symmetry-breaking distortions that lower the spin-space-group (SSG) symmetry, admits d-wave altermagnetic components, and in lower-symmetry regimes, even allow s-wave spin splitting (net magnetization). The coexistence of ferroelectric signatures and altermagnetic order establishes local lattice symmetry in MnTe as a control knob for altermagnetic spin splitting, spin current generation, and multiferroic memory applications.

We explicitly account for electron-electron interactions when modeling low-dimensional helical organic molecules. We show that competition between various hopping channels, together with interaction-induced double- and superexchange mechanisms, can stabilize non-collinear helical magnetic order. The resulting single-electron bands exhibit partial spin polarization, a manifestation of $p$-wave magnetism. Using density-matrix renormalization group, cluster perturbation theory, and Monte Carlo methods, we find that even vanishingly small spin-orbit coupling triggers strong spin selectivity at temperatures significantly above the spin-orbit scale. While strong correlations are essential for this mechanism, long-range spin ordering is not required. We thus propose non-collinear spin correlations driven by Coulomb interactions as an explanation of chirality-induced spin selectivity and discuss connections to experiments.

The intriguing physics and rich application potential of strongly correlated first-row transition metal oxide compounds result from the complex interplay of several factors that influence the electronic structure. To shed light on the effect of composition, structure, and correlation strength, we apply a well-established charge self-consistent combination of density functional theory and dynamical mean field theory, which has proven to give electron binding energies in good agreement to experimentally derived excitation spectra. For paramagnetic NiO and CoO, we analyze the effect of rock-salt and zincblende structures and their different ligand fields on the spectral functions. By varying the value of the interaction parameter U, different correlation strengths among the transition-metal 3d electrons are considered, as well as the effect of additionally accounting for correlations in the oxygen 2p orbitals by a self-interaction-correction pseudopotential scheme.

Altermagnetism is characterized by even-parity spin-momentum locking in spin-split bands despite zero net magnetization and negligible spin-orbit coupling. Here, we formulate a microscopic framework that links altermagnetic splitting in two-dimensional (2D) square lattices to orbital character. Using tight-binding models and symmetry analysis, we show that, within the minimal antiferromagnetic square-lattice model, single-orbital lattices remain spin-degenerate, whereas interwoven dual-orbital configurations lift Kramers degeneracy and generate d-wave or g-wave altermagnetic states. The spin-splitting originates from orbital anisotropy in the same-spin hopping channels. Guided by this framework, we identify M-TCNX (M = Cr, Mn, Fe; TCNX = TCNE, TCNQ) metal-organic framework monolayers with mcm topology as candidate g-wave altermagnets. Our work provides a symmetry-explicit wavefunction-level design framework for orbital-controlled altermagnetism in 2D square lattices.

Moiré superlattices formed by transition metal dichalcogenide (TMD) heterobilayers provide a versatile platform for studying strongly correlated electronic, excitonic, and topological phenomena in solids. In particular, angle-aligned MoSe$_2$/WS$_2$ heterobilayers, which have a Type-I band alignment at zero vertical electric field, host rich correlated spin and charge physics. Here, combining large-scale first-principles calculations and optical reflection spectroscopy, we report a thorough study of the emergent moiré excitonic states and interlayer charge-transfer states in angle-aligned electron-doped MoSe$_2$/WS$_2$ moiré superlattices. The moiré excitonic states serve as sensitive optical probes to the localization profile of doped electrons. We observe a series of interlayer charge-transfer transitions from n/n$_0$ = 1 to 4 (where n$_0$ denotes the moiré density) when the vertical electric field switches the heterostructure band alignment from Type-I to Type-II. By tuning the vertical electric field, we can precisely control the interlayer electron localization, realizing a Fermi-Hubbard model with a tunable charge-transfer band on an effective honeycomb lattice. Furthermore, Monte Carlo simulation of the doping dependence of the electric-field susceptibility predicts that multiple correlated charge-ordered states appear at both integer and fractional fillings. Our results provide a holistic understanding of the emergent optical excitations and the correlated charge-transfer states in electron-doped MoSe$_2$/WS$_2$ moiré superlattices.

Spin-band splitting is a hallmark of altermagnetism, intrinsically linked to magnetic ordering driven by electron correlations. However, recent inconsistencies in the detection of altermagnetism in strongly correlated altermagnet candidates have cast doubt on the robustness of this phenomenon and its dependence on many-body effects. Here, using state-of-the-art quantum many-body frameworks, we dissect the electronic origins of altermagnetism in three prototypical candidates: MnF$_2$, MnTe, and RuO$_2$. In MnF$_2$, we identify pronounced local electron correlations within Mn-3$d$ states and uncover a distinct Mott gap in the visible range, rooted in nonlocal screening effects. The strong correlations markedly localize the Mn-3$d$ electrons, leading to a narrowing of the spin-resolved bandwidth and, consequently, a suppression of spin-band splitting. By contrast, MnTe provides an ideal platform for altermagnetism, exhibiting substantial local Mn-3$d$ magnetic moments due to the strong correlations and pronounced spin-band splitting, enabled by robust Mn 3$d$--Te-5$p$ orbital hybridization. RuO$_2$ manifests as a Pauli paramagnet with vanishing local moments, even in its antiferromagnetic phase. Nonetheless, it exhibits significant spin-band splitting, indicative of itinerant altermagnetic behavior. Our results reveal that both strong local electron correlations and judicious ligand selection to promote orbital hybridization are key prerequisites to realizing altermagnetism in strongly correlated systems. These insights pave the way for the rational design and discovery of novel altermagnetic materials.

We show that locality provides a natural principle to hierarchically organize backflow wavefunctions. This leads us to propose a family of variational fermionic states, termed hierarchical backflow (HB) wavefunctions. The expressive power of HB is systematically improvable, controlled by a path depth $K$ which reflects the range of backflow correlations. At half-filling, the HB with $K=1$ already achieves high energy precision, with an accuracy around $0.5\%$ for system sizes from $4\times 4$ to $10\times 10$. At hole doping $n_h=0.125$, the method scales efficiently to $12\times16$ and $16\times16$ systems, and the energy systematically achieves higher accuracy with $K$ increasing, yielding a clear stripe phase. The HB further enables a local-nonlocal decomposition, naturally bridging to neural quantum states, while featuring compact representations and efficient optimization. Our work reveals locality as a natural organizing principle of backflow wavefunctions, opening a new framework with systematic improvability and interpretability for large-scale simulations of correlated fermion systems.

The superconducting diode effect (SDE), characterized by a nonreciprocal critical current in superconductors, has recently been observed in strongly correlated electron systems and near quantum criticality, pointing to unconventional mechanisms beyond weak-coupling theories. Here we investigate the SDE in the Rashba-Zeeman-Hubbard model, which captures $d$-wave superconductivity in an antiferromagnetic quantum critical regime, using the Dyson-Gor'kov equation with the fluctuation exchange approximation. We show that electron correlations suppress the conventional intrinsic SDE arising from depairing currents. More importantly, a supercurrent nonreciprocally induces antiferromagnetic order, which fundamentally governs the critical current and enables perfect diode efficiency. Our results reveal a previously unrecognized correlation-driven mechanism of the SDE and establish strongly correlated superconductors as a platform for superconducting diode physics.