arXiv Curator

Ambient-pressure superconductivity in compressively strained bilayer nickelate films enables direct spectroscopic tests of pairing scenarios, yet how magnetism evolves with carrier doping remains largely unexplored. Here we use Ni $L_3$-edge resonant inelastic x-ray scattering (RIXS) to track electronic and spin excitations in coherently strained La$_{3-x}$Sr$_x$Ni$_2$O$_7$/SrLaAlO$_4$ thin films ($x=0$, $0.09$, $0.21$ and $0.38$), spanning superconducting and overdoped non-superconducting regimes at essentially fixed epitaxial strain. Transport confirms superconductivity for $x\le0.21$ and a weakly insulating normal state at $x=0.38$. The $dd$-excitation manifold evolves weakly up to $x=0.21$, whereas the $\sim0.4$ eV and $\sim1.6$ eV features broaden and lose intensity at $x=0.38$. In the superconducting films, dispersive spin excitations persist along both $[H, H]$ and $[H, 0]$ with nearly doping-independent undamped dispersions and only a modest reduction of spectral weight, consistent with robust double-stripe correlations. By contrast, at $x=0.38$ the magnetic response becomes strongly broadened and weakened, with enhanced damping and $\sim50\%$ lower spectral weight, indicating a collapse of coherent double-stripe spin excitations. The concomitant suppression of magnetic coherence and superconductivity establishes a direct doping-controlled link between magnetism and superconductivity in bilayer nickelate films.

The discovery of superconductivity in bilayer nickel-oxides has revived an intense effort to understand the potential of high-temperature superconductivity in these materials and their relation to cuprate superconductors. In this work, we investigate the growth and properties of bilayer La$_3$Ni$_2$O$_7$ thin films as a function of substrate, oxygen treatment and applied pressure in order to study the evolution of transport properties. We report epitaxial growth of La$_3$Ni$_2$O$_7$ thin films on LaAlO$_3$ (LAO) (001) and SrLaAlO$_4$ (SLAO) (001) substrates, and the effects of ex-situ annealing in a high pressure furnace under an oxygen-rich environment. Transport measurements show that the La$_3$Ni$_2$O$_7$ thin films on LAO(001) exhibit Fermi liquid-like metallic behavior with a slight Kondo-like upturn at low temperatures, which evolves with the application of modest hydrostatic pressures toward non-Fermi liquid behavior with a temperature dependence of resistance approaching $\sim$ T$^{1.4}$ at 1.41 GPa. The ability to tune the normal state resistivity of La$_3$Ni$_2$O$_7$ films to display non-Fermi liquid behavior under such a modest hydrostatic pressure range - only 6 - 8 % of that typically applied via diamond anvil cell (DAC) in La$_3$Ni$_2$O$_7$ single crystals to achieve comparable effects - is both noteworthy and unexpected. These findings imply the strong tunability of La$_3$Ni$_2$O$_7$ in thin film form and the likely proximity of a strongly fluctuating ordered state leading to non-Fermi liquid behavior under even modest applied pressures.

Quantum impurity solvers are the computational bottleneck of quantum embedding approaches to correlated materials, such as dynamical mean-field theory (DMFT). We show that neural networks trained on synthetic, material-agnostic data learn the impurity mapping from hybridization functions and local interactions to Green's functions with quantitative accuracy for both model systems and real materials, providing fast solvers for single- and multi-orbital models. Benchmarks against numerically controlled quantum Monte Carlo show that the method reproduces the Mott transition, multi-orbital phase diagrams of Hubbard-Kanamori models, and the electronic properties of SrVO$_3$ and SrMnO$_3$. The learned solvers achieve orders-of-magnitude speedup and can initialize controlled calculations, dramatically accelerating DMFT while preserving accuracy.

Charge density waves (CDW) are ubiquitous in the complex phase diagram of cuprate superconductors and exhibit both short- and long-range correlations. Using time-resolved resonant X-ray scattering, we investigate the photo-induced dynamics of CDW in YBa$_2$Cu$_3$O$_{6.67}$. We discover an excitation threshold ($Φ$$_\mathrm{C}$ $\approx$ 65 $μ$J/cm$^2$) above which long-range CDW disappear, revealing a persistent CDW peak with short-range correlation length. Ultrafast photo-excitation promptly uncovers this residual short-range CDW correlations, appearing within $\approx$ 0.2 ps. Long-range CDW coherence recovers within $\approx$ 0.6 ps, while the peak intensity remains partially suppressed. We rationalize the dichotomic behavior in the fluence and temporal dependencies as the signature of two coexisting CDW peaks, arising from short- and long-range correlations, which we disentangle through their distinct response to photo-excitation. We provide evidence that the collapse of long-range correlations is driven by an electronic process, while short-range correlations are characterized by distinct timescales and stiffness against photo-excitation. This approach establishes ultrafast X-ray scattering as an effective tool for disentangling coexisting density waves and correlations in quantum materials.

Competition and coexistence of charge orders and superconductivity are hallmarks in many strongly correlated electron systems. Here, we unravel the precise role of charge fluctuations on the superconducting state in the $t$-$t'$-$J$ model of the high-temperature cuprate superconductors. Using finite-temperature tensor network simulations, we investigate thermal snapshots in the underdoped regime where the ground state features a superconducting stripe phase. At intermediate temperatures, where stripes have melted and hole clustering is observed, we find that pairing correlations are tightly localized on the hole clusters. Upon entering the stripe regime at lower temperatures, pairing increasingly delocalizes across different hole clusters to ultimately become coherent across the full system in the ground state. This pair-charge locking gives rise to an intuitive picture of the parent state of the superconducting stripe phase: pairing is localized on hole clusters formed via hole attraction due to the onset of magnetic correlations at intermediate temperature. We discuss how this microscopic picture is consistent with a broad range of experimental observations in cuprate superconductors, including scanning tunneling microscopy (STM) evidence for local pairing above $T_c$ and nuclear magnetic resonance (NMR) signatures of charge clustering in the underdoped regime.

Dynamical and variational frameworks have long been viewed as distinct paradigms. In particular, in quantum embedding (QE) frameworks, dynamical mean-field theory (DMFT) captures nonperturbative dynamical correlations through a frequency-dependent self-energy, while the Gutzwiller approximation (GA) is formulated in terms of a variationally optimized ground-state wavefunction. Here we bridge these perspectives, proving that the ghost-Gutzwiller approximation (ghost-GA), which also admits a density-matrix-matching QE formulation known as ghost density matrix embedding theory (ghost-DMET), becomes strictly equivalent to DMFT in the limit of infinitely many auxiliary bath modes. This formal unification has immediate consequences. In particular, it yields a rigorous finite-temperature extension of ghost-GA and shows that the physical Green's function can be determined from static expectation values of the embedding Hamiltonians, providing a route to computational studies of competing phases in strongly correlated matter with DMFT-level accuracy, while bypassing the need to calculate dynamical spectra with conventional impurity solvers. More broadly, it shows that the variational ghost-GA, the density-matrix-matching ghost-DMET formulation, and the dynamical DMFT description are not separate constructions, but complementary formulations of the same QE structure, thereby providing a concrete formal basis for future controlled extensions beyond DMFT.

A muon-spin rotation/relaxation ($μ$SR) study of the trilayer Ruddlesden--Popper nickelate Pr$_4$Ni$_3$O$_{10}$ was performed at ambient pressure and under hydrostatic pressure up to 2.2 GPa. Three magnetic transitions were identified at ambient pressure: the onset of spin-density-wave (SDW) order at $T_{\rm SDW} \simeq 158$ K, an intermediate-temperature transition at $T^{\ast} \simeq 90$--100 K, and a low-temperature transition at $T_{\rm SDW}^{\rm Pr} \simeq 25$--27 K. While the intermediate transition at $T^{\ast}$ induces only minor changes in the internal-field distribution, the transition at $T_{\rm SDW}^{\rm Pr}$ is accompanied by a pronounced reconstruction of the magnetic structure, consistent with previous reports attributing enhanced interlayer coherence to the ordering of the Pr sublattice. The high-temperature transition at $T_{\rm SDW}$ is characterized by the sharp development of static internal magnetic fields with a narrow transition width of $0.65(4)$ K. Weak-transverse-field measurements reveal a finite thermal hysteresis of $0.27(6)$ K, with $T_{\rm SDW}^{\rm warming} > T_{\rm SDW}^{\rm cooling}$, indicating weakly first-order-like behavior. Hydrostatic pressure suppresses $T_{\rm SDW}$ linearly and reduces the ordered Ni magnetic moment $M$, with corresponding rates of ${\rm d}T_{\rm SDW}/{\rm d}p = -4.9(1)$ K/GPa and ${\rm d}\ln M/{\rm d}p = -2.0(5)\times10^{-2}$ GPa$^{-1}$, respectively, thereby demonstrating a gradual weakening of the spin-density-wave instability under compression.

The emergent behavior of quantum materials is governed by their electronic structure, which can be experimentally probed by photoemission spectroscopy techniques that generate a four-dimensional dataset of energy and momentum. However, the quantitative extraction of Hamiltonian parameters from these high-dimensional spectra remains a significant challenge, currently relying on labor-intensive, expert-dependent analysis rather than standardized workflows. Here, we introduce a deep learning framework based on implicit neural representations to accelerate the retrieval of Hamiltonian parameters in two types of transition-metal oxides: perovskite nickelates and manganites. Our approach outperforms traditional analytical fitting procedures, yielding superior agreement with experimental Fermi surface topologies and energy-momentum dispersions. This work highlights the potential of deep learning tools to bridge the gap between theory and experiment, paving the way for high-throughput, autonomous discovery pipelines in quantum materials.

Engineering artificial systems by twisting and stacking van der Waals materials has proven to be an excellent platform for exploring emergent quantum phenomena that can be significantly different from the constituents. Recent advances in the fabrication of high-quality twisted interfaces provide a unique opportunity to study the little-explored interfacial superconducting order in twisted cuprate superconductors. In our work, we fabricate superconducting quantum interference devices (SQUID) that utilize the twisted interface of $\mathrm{Bi_2Sr_2CaCu_2O_{8+δ}}$, a high-Tc cuprate superconductor. By measuring the magnetic field modulation of switching current and differential resistance, we find a $\mathrmπ$ phase difference between the two Josephson junction arms of the SQUID reflecting chiral superconducting order -- a crucial aspect inaccessible to single Josephson junction devices of the past. Our observations also indicate co-tunneling of the Cooper pairs and a time-reversal symmetry-broken emergent superconducting order. Additionally, these SQUIDs are well suited for use as state-of-the-art flux sensors close to 77 K, achieving a flux noise sensitivity of $\sim$1.5 $\mathrm{μΦ_0/\sqrt{Hz}}$. Stabilizing new superconducting orders using twisted interfaces and probing them using quantum interference opens new avenues to understanding the microscopic origin of unconventional superconductors. Our SQUID architecture is suitable for investigating the charge transport mechanisms and the symmetry of superconducting order at the interfaces of other systems, reflecting the broad applicability beyond cuprate superconductors.

Recently, a series of experiments which control the T$_\text{c}$ of the bilayer nickelates La$_3$Ni$_2$O$_7$ through varying environmental conditions, including the rare-earth Sm/Nd substitution, the pressure on the bulk material, the compressive strain on the film and the hole doping through over-oxidation or alkaline earth element substitution have caught great interests. Here, we provide a unified understanding toward all these experiments based on the minimal single $d_{x^2-y^2}$-orbital bilayer $t-J_\parallel-J_\perp$ model proposed previously. With model parameters input from density-functional-theoretical calculations under varying experimental conditions, we adopt combined slave-boson-mean-field and density-matrix-renormalization-group approaches to solve the model and compare with experiments. Our results yield that, the bulk T$_\text{c}$ under pressure enhances with the Sm/Nd substitution fraction, the bulk T$_\text{c}$-pressure relation takes a dome shaped curve, the T$_\text{c}$ of the thin film enhances with compressive strain. The obtained parameters dependence of T$_\text{c}$ in these three experiments mainly originates from the variation of $J_\perp$ with experimental conditions. As for the hole doping, our results provide that T$_\text{c}$ decreases with the hole doping level $δ$, due to reduced density of state for the $d_{x^2-y^2}$-orbital. All these results are qualitatively consistent with experiments. We further conduct a comparative weak-coupling random-phase-approximation (RPA) based study on these experiments and find that our strong-coupling $t-J_\parallel-J_\perp$ model provides a more natural understanding of the experiments. We propose that electron doping implemented through substitution of La by element with higher valence, or further enhancement of the compressive strain in the film, can enhance T$_\text{c}$.

First-order structural phase transition is a common phenomenon in materials that qualitatively alters their physical properties. Yet, the abrupt first-order nature is usually unexplained by realistic computations, implying an omission of important physics in describing the electronic structure of the nearby stable phases. Using the recently discovered nickelate superconductors La$_3$Ni$_2$O$_7$ as a prototypical example, we demonstrate that such first-order nature is typically beyond intra-atomic correlation considered in state-of-the-art material computations. Instead, a full many-body treatment of low-energy active orbitals reveals a generic inter-atomic "orbital dimerization" mechanism of first-order structural phase transition, corresponding to abrupt energy reduction upon a spin-singlet bond formation. Such an inter-atomic correlation qualitatively changes not only the essential lattice bonding but also the characteristics of low-energy electronic properties across the transition. This strong mechanism and the developed computational framework are generally applicable to a wide variety of ionic materials, to produce valuable insights into atomic and electronic structures essential for their physical properties and functionalities.

This review summarizes recent progress of ambient-pressure superconductivity in bilayer nickelate La$_3$Ni$_2$O$_7$ thin films, a major advancement following the discovery of high-pressure superconductivity in bulk La$_3$Ni$_2$O$_7$. First, we explain how epitaxial strain engineering enables ambient-pressure superconductivity in La$_3$Ni$_2$O$_7$ thin films, with compressive strain from substrates like SrLaAlO$_4$ stabilizing superconductivity. Next, we review experimental characterizations of related systems, with particular emphasis on ARPES measurements that have shown conflicting Fermi surface topologies. We then discuss progress in increasing the superconducting transition temperature $T_c$. Finally, we summarize theoretical studies of the electronic structure and pairing symmetry of La$_3$Ni$_2$O$_7$ thin films. Together, these advances establish bilayer nickelate thin films as a highly tunable and promising platform for exploring high-$T_c$ superconductivity.

Competing charge and spin orders are central to uncovering the nature of unconventional superconductivity. Here we utilize synchrotron X-ray diffraction on a high-quality single crystal to reveal the charge order of La$_3$Ni$_2$O$_7$ at ambient pressure, which competes with the high-temperature superconducting phase under pressure. Enabled by the high synchrotron photon flux and a large dynamic range, we resolve faint reflections -- nearly four orders of magnitude weaker than the main Bragg reflections -- that were overlooked in prior diffraction studies. This observation evidences a broken glide-mirror symmetry, leading to a polar crystal structure, rather than the widely used centrosymmetric structure model. The polarity is induced by checkerboard charge order on nickel sites in combination with octahedral tilting, reminiscent of bilayer manganese oxides. Our results provide a foundation for understanding phase competition and the mechanism of pressure-induced superconductivity in bilayer nickelates.

High-pressure superconductivity in the bilayer nickelate La$_3$Ni$_2$O$_7$, with a transition temperature approaching 80 K, has stimulated intense debate regarding its microscopic origin. Although an $s_{\pm}$ gap symmetry has been widely proposed, the electronic degrees of freedom responsible for pairing remain unsettled. Here we investigate a bilayer two-orbital Hubbard model using the variational Monte Carlo method and reveal a hierarchical pairing structure in bilayer nickelates. The primary pairing interaction originates from the bonding--antibonding splitting of the Ni $3d_{z^2}$ orbitals, while orbital hybridization redistributes superconducting correlations to the $d_{x^2-y^2}$ channel despite its weak intrinsic pairing interaction. This distinction between the origin of pairing and resulting superconducting correlations explains why the two orbital channels exhibit comparable long-range correlations. The resulting $s_{\pm}$ state is robust against changes in Fermi-surface topology. These results reconcile apparently competing theoretical scenarios and provide a comprehensive understanding, highlighting the distinctive role of orbital hybridization in multilayer correlated superconductors.

The discovery of superconductivity in reduced square-planar nickelates marked a major advance in identifying structural and electronic analogs to the high-$T_c$ cuprates. The more recent observation of superconductivity in parent Ruddlesden-Popper (RP) octahedral nickelates with a clear difference in electron count with respect to cuprates raises new questions about the nature of superconductivity across these related but distinct nickelate families. Here, we use Ni $L_3$-edge resonant inelastic x-ray scattering (RIXS) to probe the low-energy excitations in a representative compound of both families: the parent octahedral $n=8$ RP phase Nd$_9$Ni$_8$O$_{25}$ (p-RP), which is non-superconducting, and its reduced square-planar counterpart Nd$_9$Ni$_8$O$_{18}$ (r-RP), which exhibits superconducting correlations with a $T_c \approx 5$ K. The $n=8$ p-RP develops a spin-density-wave (SDW) ground state with ordering wave vector $q_{\mathrm{SDW}} = (1/4,\, 1/4)$, analogous to the bilayer RP, while the $n=8$ r-RP shows an elastic peak at $q^\star = (1/3,\, 0)$. Polarimetric RIXS shows that the p-RP exhibits low-energy spectra dominated by weakly dispersive paramagnons along the 0$\rightarrowπ$ and $π\!\rightarrow\! π$ directions, whereas the r-RP with superconducting correlations displays dispersionless magnetic excitations. Our results comprehensively map out the spin excitations and reveal fundamental differences in the ground state between these two distinct structural families.

Quantum impurity solvers are the computational bottleneck of quantum embedding approaches to correlated materials, such as dynamical mean-field theory (DMFT). We show that neural networks trained on synthetic, material-agnostic data learn the impurity mapping from hybridization functions and local interactions to Green's functions with quantitative accuracy for both model systems and real materials, providing fast solvers for single- and multi-orbital models. Benchmarks against numerically controlled quantum Monte Carlo show that the method reproduces the Mott transition, multi-orbital phase diagrams of Hubbard-Kanamori models, and the electronic properties of SrVO$_3$ and SrMnO$_3$. The learned solvers achieve orders-of-magnitude speedup and can initialize controlled calculations, dramatically accelerating DMFT while preserving accuracy.

It is generally believed that the spin-$\tfrac{1}{2}$ triangular-lattice $J_1$-$J_2$ Heisenberg model hosts a quantum spin liquid in the intermediate regime between the $120^\circ$ and stripe ordered phases. Density matrix renormalization group studies on cylinders have consistently found two nearly degenerate ground states, commonly interpreted as distinct topological sectors. Using state-of-the-art matrix product state simulations on YC6 cylinders, we compare the static and dynamical properties of these two sectors at $J_2/J_1 = 0.125$. Noticeable differences appear already in static correlations; moreover, high-resolution dynamical structure factors reveal qualitatively distinct low-energy excitations. These results suggest that the two ground states cannot be understood as merely topologically distinct sectors of a gapped $\mathbb{Z}_2$ spin liquid.

It is generally believed that the spin-$\tfrac{1}{2}$ triangular-lattice $J_1$-$J_2$ Heisenberg model hosts a quantum spin liquid in the intermediate regime between the $120^\circ$ and stripe ordered phases. Density matrix renormalization group studies on cylinders have consistently found two nearly degenerate ground states, commonly interpreted as distinct topological sectors. Using state-of-the-art matrix product state simulations on YC6 cylinders, we compare the static and dynamical properties of these two sectors at $J_2/J_1 = 0.125$. Noticeable differences appear already in static correlations; moreover, high-resolution dynamical structure factors reveal qualitatively distinct low-energy excitations. These results suggest that the two ground states cannot be understood as merely topologically distinct sectors of a gapped $\mathbb{Z}_2$ spin liquid.

The cuprate system La$_{2-x}$Ba$_{x}$CuO$_{4}$ (LBCO) exhibits a pronounced sensitivity to in-plane uniaxial stress, particularly near the 1/8 doping anomaly, where stripe order strongly suppresses bulk superconductivity. While previous studies have focused on compositions close to 0.125, the commensurate $x$=0.125 phase remains largely unexplored under symmetry-selective lattice tuning. Here, we combine muon-spin rotation ($μ$SR), AC susceptibility, and electrical resistivity to investigate superconductivity, spin-stripe order, and structural response in LBCO-0.125 under in-plane uniaxial stress applied 45$^\circ$ to the Cu-O bond direction. Complementary resistivity measurements on $x$=0.115 and 0.135 track the evolution across both sides of the anomaly. We observe a giant enhancement of the bulk superconducting transition temperature in LBCO-0.125, increasing from 5 K to 37 K under 0.5 GPa. While the onset temperature of spin-stripe order decreases only modestly, the magnetic volume fraction is reduced by about a factor of two, with local order preserved. Simultaneously, the resistivity peak associated with the LTT phase is fully suppressed across all dopings. These results demonstrate that suppression of the LTT phase and reduction of the static spin-stripe-ordered volume fraction are crucial for the development of optimal three-dimensional superconductivity. Strikingly, the composition $x$=0.125, with the most robust stripe stability and the lowest ambient-pressure $T_{\rm c}$ develops the highest $T_{\rm c}$ under stress, reaching a zero-resistance state at 37 K and an onset of the superconducting transition as high as 46 K. This behavior indicates that stripe-related interactions enhance pairing strength, while static stripe order competes with superconductivity primarily at the level of phase coherence rather than pairing itself.

Motivated by recent angle-resolved photoemission spectroscopy (ARPES) experiments on infinite-layer (IL) nickelates, we employ determinant quantum Monte Carlo (DQMC) to study the three-orbital Emery model ($d$-$p$ model) coupled to an additional interstitial $s$ orbital retaining the three-dimensional dispersion. Our large-scale simulations reveal that: (1) the interstitial $s$-orbital-derived electron pocket is significantly reduced by the strong interaction but persists upon 20\% hole doping, reaching a size comparable to experimental observations; (2) the $d_{x^2-y^2}$-orbital dispersion is strongly renormalized by interactions, leading to a weak $k_z$ dependence consistent with ARPES measurements. Furthermore, compared with the conventional three-orbital $d$-$p$ model, the $d$-$p$-$s$ model exhibits enhanced short-range antiferromagnetic correlations. These results highlight the crucial role of strong correlations and multi-orbital effects in shaping the low-energy electronic structure and many-body correlations in IL nickelates, and demonstrate the necessity of treating interaction-driven many-body physics within a realistic multi-orbital framework.

We introduce Neural Tensor Network States ($ν$TNS), a variational many-body wave-function ansatz that integrates deep neural networks with tensor-network architectures. In the $ν$TNS framework, a neural network serves as a disentangler of the wave-function, transforming the physical degrees of freedom into renormalized variables with much less entanglement. The renormalized state is then efficiently encoded by a back-flow tensor network. This construction yields a compact yet highly expressive representation of strongly correlated quantum states. Using convolutional neural networks combined with matrix product states as a concrete implementation, we obtain state-of-the-art variational energies for the spin-$1/2$ $J_1$-$J_2$ Heisenberg model on the square lattice at the highly frustrated point $J_2/J_1=0.5$, for systems up to $20\times 20$ with periodic boundary conditions. Finite-size scaling of spin, dimer, and plaquette correlations exhibits power-law decay without magnetic or valence-bond long-range order, consistent with a gapless quantum spin-liquid ground state at that point.This $ν$TNS framework is flexible and naturally extensible to other neural and tensor-network structures, offering a general platform for investigating strongly correlated quantum many-body systems.

This thesis investigates the magnetic, spectral, and transport properties of strongly correlated electronic systems, with a primary focus on the Hubbard model and its extensions relevant for real materials. Within the dynamical mean-field theory (DMFT) framework, different regimes of interaction strength, temperature, doping, and magnetic field are explored, highlighting the central role of local electronic correlations in shaping spectral reconstruction and nontrivial transport responses. For the antiferromagnetic Hubbard model under a Zeeman field, magnetoresistance and local metamagnetism are characterized, revealing the coexistence of distinct energy scales associated with charge and spin degrees of freedom. A minimal, purely correlation-driven mechanism for generating spin-polarized charge transport in structurally conventional collinear antiferromagnets is identified, controlled by the simultaneous breaking of particle--hole symmetry and antiferromagnetic sublattice equivalence. Finally, these concepts are applied to correlated materials with strong spin--orbit coupling, such as Sr$_2$IrO$_4$ and Sr$_3$Ir$_2$O$_7$, and to nanoparticle solids dominated by Coulomb blockade and disorder. The results show how ideas developed in correlated lattice models provide a unified interpretation of metal--insulator transitions and spectral reconstruction in complex systems.

Quantum spin liquids (QSLs) are exotic phases of matter characterized by long-range entanglement and the absence of magnetic order even at zero temperature. Here, we present a comprehensive theoretical study of the frustrated magnet Y$_{3}$Cu$_{2}$Sb$_{3}$O$_{14}$ to elucidate its electronic and magnetic properties. We uncover two inequivalent Cu sites with fundamentally distinct oxygen coordination, i.e. octahedral CuO$_6$ and axially compressed CuO$_8$, that give rise to completely opposite crystal-field splittings. This inversion places the unpaired hole in the $d_{z^2}$ orbital at the Cu-2 site, while Cu-1 maintains conventional $d_{x^2-y^2}/d_{xy}$ character, producing narrow bands at commensurate filling. Dynamical mean-field theory (DMFT) calculation reveals a selective band-renormalization of orbitals from the two Cu ions with Cu-1 undergoing a Mott transition while Cu-2 remaining metallic under electronic correlations. Crucially, FLEX calculations demonstrate that multiple magnetic instabilities compete with nearly equal strength: the spin susceptibility lacks dominant peaks, and the leading eigenvalues approach unity simultaneously across all wavevectors with increasing interaction. This competitive interplay, originating from the distinct local environments and geometric frustration on the triangular lattice, agrees well with the absence of long-range magnetic order in experiment. Our results suggest Y$_{3}$Cu$_{2}$Sb$_{3}$O$_{14}$ as a promising QSL candidate where the unique combination of disparate crystal-field environments, strong correlations, and competing exchange interactions conspire to stabilize an exotic quantum ground state.

The square-planar layered nickelates R$_{n+1}$Ni$_n$O$_{2n+2}$ (R= Nd, $n=4-7$) have been recently shown to be superconducting without the need for chemical doping or pressure. Here, we examine the electronic structure of the analog higher-order square-planar palladates --that have not yet been synthesized-- via \textit{ab initio} calculations. These layered palladates exhibit larger bandwidths, an increased $p-d$ hybridization, and less interference from R-$d$ bands at the Fermi level. These characteristics make them closer cuprate analogs and promising candidates to pursue in the context of unconventional superconductivity.

Because the normal state of underdoped cuprate superconductors is an enigmatic Fermi-arc metal, it is valuable to analyze an exactly solvable model that exhibits both Fermi arcs and $d$-wave superconductivity. Here, we focus on a recently proposed solvable model in which the emergence of Fermi arcs is especially transparent. Upon incorporating a $d$-wave pairing interaction, the model produces an asymptotically exact solution for the superconducting transition temperature $T_c$ that traces out a superconductivity dome as a function of hole doping, in qualitative agreement with experimental observations in cuprates. Crucially, we show analytically that the Fermi arcs generate an additional many-body effect that suppresses $T_c$ beyond the simple reduction expected from a shrinking Fermi surface. The many-body nature of the Fermi arcs further introduces the gap-to-$T_c$ ratio greatly surpassing the mean-field limit. These findings provide an analytic benchmark for understanding how Fermi-arc physics competes with $d$-wave superconductivity in high-$T_c$ superconductors.

The microscopic origin of high-temperature superconductivity in cuprates remains one of the central open questions in condensed matter physics. Growing experimental and theoretical evidence suggests that the bare single-band Fermi-Hubbard model may not fully capture properties of cuprates such as superconductivity, motivating us to revisit the canonical three-band model of the copper-oxide planes - the Emery model - from which the single-band counterpart was originally derived. Here, we propose and analyze a quantum simulation scheme for realizing the Emery model in regimes relevant to cuprates and infinite-layer nickelates with today's ultracold atom quantum simulation platforms, enabling the exploration of the three-band physics on system sizes that are challenging for current numerical methods. Specifically, we show that a two-dimensional optical lattice with a superimposed pattern of repulsive potentials can be designed to study low-temperature properties for variable parameter regimes of the Emery model relevant to cuprates as well as infinite-layer nickelates. Our results pave the way for real material simulations with ultracold atom quantum simulators and a better understanding of the physics of unconventional superconductors.

Antiferromagnets offer an attractive platform for spintronics due to their absence of net magnetization and ultrafast spin dynamics, yet their intrinsically spin-compensated electronic structure has traditionally limited their active role in spin transport. Here we identify a minimal, correlation-driven route to spin-polarized charge transport in collinear antiferromagnets. Using the doped antiferromagnetic Hubbard model within dynamical mean-field theory, we show that electronic correlations generate strong spin-dependent scattering upon doping away from half filling, while a uniform magnetic field lifts the residual symmetries that enforce spin-degenerate transport. Only the combined breaking of particle--hole symmetry by doping and of the antiferromagnetic sublattice equivalence by the applied magnetic field converts these dynamical asymmetries into a finite spin polarization of the charge current. Our results establish electronic correlations as an active ingredient for antiferromagnetic spintronics and reveal a correlated analogue of the symmetry-breaking mechanism underlying altermagnetic spin-polarized transport in structurally conventional, collinear antiferromagnets.

Antiferromagnets offer an attractive platform for spintronics due to their absence of net magnetization and ultrafast spin dynamics, yet their intrinsically spin-compensated electronic structure has traditionally limited their active role in spin transport. Here we identify a minimal, correlation-driven route to spin-polarized charge transport in collinear antiferromagnets. Using the doped antiferromagnetic Hubbard model within dynamical mean-field theory, we show that electronic correlations generate strong spin-dependent scattering upon doping away from half filling, while a uniform magnetic field lifts the residual symmetries that enforce spin-degenerate transport. Only the combined breaking of particle--hole symmetry by doping and of the antiferromagnetic sublattice equivalence by the applied magnetic field converts these dynamical asymmetries into a finite spin polarization of the charge current. Our results establish electronic correlations as an active ingredient for antiferromagnetic spintronics and reveal a correlated analogue of the symmetry-breaking mechanism underlying altermagnetic spin-polarized transport in structurally conventional, collinear antiferromagnets.

Machine learning and high-throughput screening approaches to superconductor discovery require physically meaningful descriptors that capture essential physics while remaining computationally tractable. The superfluid weight is an ideal descriptor as it is a prerequisite for superconductivity, determines the magnetic penetration depth and the Berezinskii-Kosterlitz-Thouless transition temperature in two-dimensional materials, may limit the critical temperature in unconventional superconductors through phase coherence, and reveals quantum geometric contributions to supercurrent transport. We develop a computationally efficient framework for calculating the zero-temperature, mean-field superfluid weight for uniform pairing from density functional theory band structures and Bloch wavefunctions. We separately evaluate the conventional contribution from band curvature and the geometric contribution from quantum geometry. To validate the method, we calculate London penetration depths for a few conventional superconductors (Al, Pb, Nb, MgB$_2$, LuRu$_3$B$_2$ and YRu$_3$B$_2$) and find good agreement with experiment after accounting for nonlocal corrections, strong-coupling effects, and sample quality. The conventional contribution dominates by orders of magnitude in these wide-band materials, as expected. This framework provides a foundation for large-scale screening of superconducting candidates and exploring quantum geometric effects in unconventional superconductors.

Understanding pairing in the strong-coupling regime of doped Mott insulators remains an open problem in the context of cuprate superconductors. We perform ultra-high resolution numerical simulations of spectral functions in the highly underdoped $t-J$ model and discover two coupled branches of hole pairs emerging at low energies in the largely unexplored two-particle spectrum. As spin anisotropy is tuned from the Ising limit to the $SU(2)$-symmetric Heisenberg regime, the lowest $d$-wave pair evolves from a single bipolaronic branch into two hybridized branches separated by an avoided crossing. We explain this behaviour using an effective two-channel model involving a tightly bound bipolaronic state and a second channel associated with two magnetic polarons. The model reproduces the qualitative low-energy spectra and implies near-resonant $d$-wave interactions in the $SU(2)$-symmetric $t-J$ model, consistent with proximity to an emergent Feshbach-type resonance. To probe these predictions experimentally, we propose a Raman spectroscopy scheme for the attractive Hubbard model that can be directly implemented using ultracold atoms in optical lattices. Our work establishes two-particle spectroscopy, beyond single-particle Green's functions, as a powerful tool for revealing the microscopic origins of unconventional superconductivity.

Order parameter analysis of the t/U series reveals a uniform altermagnet endemic to the doped Mott insulator, driven by kinetic interactions, occupying a position between the antiferromagnet and hole-doped d-wave superconductor that is normally reserved for the pseudogap. The metastable boundary of the altermagnet punctures and divides the superconductor into underdoped and overdoped regions, reminiscent of the $T^\ast$ crossover or transition in the cuprates. Similarly, the $T_{pair}$ boundary of the superconductor divides the altermagnet, leading to a low temperature phase susceptible to Cooper fluctuations. The altermagnet is unstable to inhomegeneous spin and charge order of sites, bonds, and currents. Its leading instability is to the $π$-flux state, suggesting the possible emergence of spin-charge liquids and quantum ordered states from a physically realistic microscopic model.

Low-dimensional materials exhibit exotic properties due to enhanced quantum fluctuations, making the understanding of their microscopic origin central in condensed matter physics. Analogue quantum simulators offer a powerful approach for investigating these systems at the microscopic level, particularly in large-scale regimes where quantum entanglement limits classical numerical methods. To date, analogue simulators have largely focused on universal Hamiltonians rather than material-specific quantitative comparisons. Here we use a Rydberg-based quantum simulator to study the bulk-layered frustrated quantum magnet TmMgGaO$_4$. Magnetisation measurements obtained from the quantum simulator show excellent agreement with independent measurements performed in a magnetic laboratory facility, validating the proposed effective two-dimensional microscopic Hamiltonian. Building on this quantitative correspondence, we investigate on both platforms the antiferromagnetic phase transition. We further probe the role of quantum fluctuations via snapshot analysis, connecting our results to integrated inelastic neutron scattering data. Finally, we access, with the simulator, non-equilibrium dynamics on picosecond material timescales, including frequency response and thermalisation of observables.

We present an SU(2) gauge theory of fluctuating stripe order in the two-dimensional Hubbard model. The theory is based on a fractionalization of the electron operators in fermionic chargons with a pseudospin degree of freedom, and charge neutral spinons capturing fluctuations of the spin orientation. The chargons are treated in a renormalized mean-field theory. We focus on regions of the phase diagram where they undergo stripe order. The spinons are described by a non-linear sigma model with pseudospin stiffnesses determined by the chargons. They prevent breaking of the physical SU(2) spin symmetry at any finite temperature, resulting in a charge ordered pseudogap phase with a reconstructed Fermi surface and a spin gap. The spectral function for single-particle excitations exhibits a collection of Fermi arcs and other structures. The arcs appear in various regions of the Brillouin zone, but never exclusively around the Brillouin zone diagonals.

The $t$-model represents the Hubbard model in the limit $U \to \infty$ and is one of the basic models of strongly correlated electrons. On a one-dimensional chain, the model is integrable, and the charge dynamics corresponds to that of free spinless fermions. However, the sequence of spins is frozen, leading to the Hilbert space fragmentation and nontrivial spin dynamics. We consider integrable and perturbed models with perturbations that break integrability while preserving fragmentation, and show that they exhibit various types of spin dynamics, from ballistic transport to anomalous diffusion in the integrable case, and from diffusion to subdiffusion in the perturbed case. Due to fragmentation, in all cases considered, spin transport is mediated by charge transport, with a particular magnetization dependence, most notably leading to subdiffusion in the grandcanonical average of the perturbed model, with a mechanism distinct from subdiffusion in disordered or dipole-conserving models.

Symmetric mass generation (SMG) transitions defy the conventional Landau-Ginzburg-Wilson paradigm by opening a many-body gap without spontaneous symmetry breaking or topological order, attracting intense interest across particle physics and condensed matter physics. Here, we utilize unbiased quantum Monte Carlo simulations to investigate the equilibrium and nonequilibrium critical dynamics of the SMG transition in a bilayer honeycomb lattice model. We unambiguously confirm the existence of an SMG transition at $J_{\text{c}}=2.584(8)$ that separates the Dirac semimetal phase from a symmetry-preserving SMG phase. High-precision extraction of the critical exponents reveals a novel universality class that profoundly departs from mean-field theory. We then extend our study to the nonequilibrium regime, exploring the driven dynamics of the SMG transition. Notably, despite the breakdown of the prerequisites for the celebrated Kibble-Zurek mechanism, the nonequilibrium SMG transition still follows the generalized finite-time scaling. By bridging equilibrium criticality and nonequilibrium dynamics, our work uncovers the universal critical properties of SMG transitions, providing a solid theoretical basis for future experimental studies of SMG physics.

The kagome lattice naturally hosts flat bands, Dirac fermions, and van Hove singularities, yet whether its geometry can stabilize orbital-selective phases - a hallmark of Hund's physics in multi-orbital correlated systems - has remained an open question. Here, we combine resonant inelastic X-ray scattering with density functional theory and dynamical mean-field theory to demonstrate that YMn$_6$Sn$_6$ exhibits a spontaneous orbital differentiation into coexisting itinerant and localized electrons within the same Mn $3d$ manifold. Orbitals directed along Mn-Mn bonds provide coherent quasiparticles and metallic bands, while those pointing toward ligands become strongly correlated and display non-Fermi-liquid behavior. Hund's intra-atomic exchange suppresses orbital fluctuations, stabilizing this dichotomy and providing a natural double-exchange-like mechanism for the observed ferromagnetic bilayer coupling. Our work establishes YMn$_6$Sn$_6$ as a kagome platform where orbital selectivity, flat-band topology, and Hund's metallicity converge - revealing that geometric frustration and correlation-driven orbital differentiation can cooperatively design exotic quantum phases beyond the canonical paradigms of Mott physics or band topology alone.

The prototypical Mott-Hubbard insulator LaVO3 undergoes a structural phase transition accompanied by the onset of spin and orbital ordering below 140 K. By combining ultrafast optical pump-probe spectroscopy and two-dimensional electronic spectroscopy, we investigate the interplay between fluctuations of the local spin and orbital order parameter and the lifetime of high-energy electron-hole excitations. Specifically, we demonstrate that the pump-induced perturbation of the order parameter leads to a change of the Hubbard exciton decoherence time and, consequently, of its homogeneous linewidth. Dynamical mean-field theory calculations confirm that the exciton scattering rate is crucially affected by the degree of order of the spin and orbital lattices in LaVO3. Our results demonstrate that multi-dimensional ultrafast optical spectroscopy can be used to track the dynamics of the order parameter, thus opening new routes in the study of correlated quantum materials characterized by intertwined orders.

First-principles analysis of the spectroscopic signatures of correlated quantum materials poses significant challenges due to the interplay between spin-orbit and vibronic couplings, as well as the need to describe both dynamic and static electron correlation to reach decent accuracy. In this work, we apply the equation-of-motion coupled-cluster (EOM-CC) method to derive the spin-orbit-lattice entangled vibronic states and predict the Re $L_3$ edge resonant inelastic x-ray scattering (RIXS) spectra of Ba$_2$MgReO$_6$. The EOM-CC yields interaction parameters in close agreement with those extracted from RIXS spectra, with errors of less than 5\%. In particular, the EOM-CC method allowed us to determine the weak vibronic coupling to the $T_{2g}$ vibrations, which is difficult to address experimentally. The simulated spectra indicate that vibronic coupling to the $T_{2g}$ modes gives rise to a shoulder on the elastic peak. Going beyond the conventional treatment, which focuses solely on $E_g$ modes, we show that vibronic couplings to both $T_{2g}$ and $E_g$ modes are required to account for the fine structure of the RIXS spectra. This work demonstrates that the EOM-CC method is a powerful tool for accurately predicting the complex local states at metal sites and spectroscopic signatures of correlated insulating materials.

We introduce Large Electron Model, a single neural network model that produces variational wavefunctions of interacting electrons over the entire Hamiltonian parameter manifold. Our model employs the Fermi Sets architecture, a universal representation of many-body fermionic wavefunctions, which is further conditioned on Hamiltonian parameter and particle number. On interacting electrons in a two-dimensional harmonic potential, a single trained model accurately predicts the ground state wavefunction while generalizing across unseen coupling strengths and particle-number sectors, producing both accurate real-space charge densities and ground state energies, even up to $50$ particles. Our results establish a foundation model method for material discovery that is grounded in the variational principle, while accurately treating strong electron correlation beyond the capacity of density functional theory.

Incoherent electronic states in strongly correlated materials are commonly attributed to disorder or material specific mechanisms. Here we show that incoherent spectra instead arise from self-generated dynamical disorder associated with competing fluctuations. In this regime, electron dynamics coupled to time-dependent scattering naturally produce a spectral function of the form rho (z) = exp(-z^2/4) Dnu (z), where z is a scaled energy and Dnu denotes the parabolic cylinder function. This form reflects a marginal dynamical regime characterized by non-Markovian temporal correlations. Applying this scaling function to angle resolved photoemission spectroscopy (ARPES) energy distribution curves from the cuprates Nd2-xCexCuO4 and Bi2Sr2CaCu2O8+delta, the Kagome metal CsCr3Sb5, and the double-layer nickelate La3Ni2O7, we find that incoherent spectra are quantitatively described by rho (z), differing only in non-universal amplitude and energy scales. After rescaling, the datasets collapse onto a single universal curve characterized by a fixed parabolic-cylinder order nu = -1/2. The observed spectral collapse indicates a fixed-point-like regime in which microscopic details such as lattice geometry, band structure, and chemical composition become irrelevant at low energies. These results establish a unified and quantitative framework for continuum-dominated ARPES spectra across diverse strongly correlated materials.

The strange metal behavior in cuprate superconductors - characterized by linear in temperature resistivity and anomalous Hall transport - stands in stark contrast to the expectation of conventional Fermi liquid (FL) theory. Remarkably, the similar transport behavior has also been observed in the heavy fermion metal CeCoIn$_5$, whose d-wave superconducting ground state and strong antiferromagnetic fluctuations draw parallels to the cuprates. Here we have investigated the optical conductivity of the strange metal state of CeCoIn$_5$ over a wide magnetic field range using time-domain THz spectroscopy (TDTS). Using unique high-field THz spectroscopy we have shown that the current relaxation rate scales approximately as T$^2$, giving evidence for a hidden Fermi liquid state over a large field range. This result can be reconciled with linear in T resistivity with the realization that heavy quasiparticles have an optical mass that scales roughly like 1/T. This optical mass contrasts with the mass that characterizes cyclotron motion, which does not suffer the same large temperature dependent renormalization. Although by itself anomalous, this allows one to understand a number of other phenomena in CeCoIn$_5$ that have been taken to be signatures of strange metals, including the coexistence of a conventional T$^2$ dependence of the cotangent of the Hall angle with the linear in T resistivity, which with our observation also reflects FL-like physics.

X-ray absorption spectroscopy (XAS) and resonant inelastic X-ray scattering (RIXS) are powerful probes of electronic structure owing to their chemical and orbital selectivity. For powder samples, however, interpreting RIXS spectral intensities remains challenging as the measured signal is an average over all orientations. Existing theoretical treatments rely largely on spherical-tensor formalisms, which often involve complex derivations and case-specific analyses. Meanwhile, recent advances in quantum-chemistry methods have made the evaluation of transition tensors in Cartesian coordinates both accurate and straightforward. Here, we present a general theoretical framework that translates Cartesian transition tensors into physically meaningful, orientation-averaged intensities for powder samples. The formalism allows predicting angular and polarization dependences \textit{ab initio} for both XAS and RIXS and is extendable to other spectroscopies. The resulting predictions show excellent agreement with RIXS experimental data at the Ce L$_3$ edge.

We study the conditions under which Matrix Product States (MPS) or Matrix Product Operators are exact eigenvectors of an extensive local operator, such as a Hamiltonian. By suitably choosing the local operator, this covers a wide range of settings: Exact eigenstates of Hamiltonians, including scar states, exact MPS trajectories for driven quantum systems, steady states of local Lindbladians, generalized symmetries of either Hamiltonians or density matrices, and many more. Our key result is that that a local, fixed-size equation -- namely, how a single term in the operator acts on a block of tensors -- provides a necessary and sufficient condition for exact solutions. This allows to characterize the full space of solutions in all of the aforementioned problems, and to identify them both analytically and numerically. We elaborate on the concrete application of this characterization to all of the aforementioned settings, and in particular exemplify the power of our local characterization by using it to recover the quantum group symmetries of the XXZ model. We also discuss applications to numerical algorithms with MPS and the generalization of our results to 2D, i.e., projected entangled pair states (PEPS).

Geometrical frustration is a powerful route to realize exotic phases such as quantum spin liquids. Despite extensive efforts, systematic searches targeting specific frustration motifs and their potential to host unconventional magnetic ground states remain rare, thus highlighting the need for a more focused and predictive materials discovery approach. Here we present a new strategy combining high-throughput first-principles calculations, magnetic force theory, and spin Hamiltonian analysis. Starting from the 150,000 material database, we catalogue candidate materials that may host competing exchange interactions and new types of magnetic states with the focus on kagome or triangular lattices. Our workflow not only reproduces the majority of known frustrated magnetic materials, validating our approach, but also predicts novel candidate compounds with targeted frustration profiles that have not yet been experimentally synthesized. Among these, we identify six promising new materials: one triangular lattice compound, KMgNiIO6, and five kagome lattice compounds; Li4Fe3WO8, Li2V3F8, Li5VP2(O4F)2, and Li2MgCo3O8 (P2/m and C2/m). For each candidate, we identify detailed magnetic properties and further propose their potential magnetic ground states, revealing that some of them may host entirely new magnetic phases driven by their distinct frustration characteristics.

Topological flat bands (FBs) offer an ideal platform for realizing exotic topological phases, such as fractional Chern insulators, yet their realization with both exact flatness and stable topology in local lattice models has been long hindered by fundamental no-go theorems. Here, we overcome this barrier by demonstrating the existence of critical topological FBs (CTFBs) in finite-range hopping models. They saturate the no-go theorems via a unique structure of Bloch wavefunctions: While continuous over the whole Brillouin zone, the wavefunctions are non-analytic at isolated band touching points, thereby relaxing the inherent restrictions on the coexistence of exact flatness and stable topology. We establish a general principle to construct CTFBs, as well as their parent Hamiltonians, that carry desired topological invariants in given space groups. Explicit examples exhibiting Chern numbers, strong $\mathbb{Z}_2$ index, and crystalline-symmetry-protected invariants in two and three dimensions are provided. Furthermore, an automated algorithm identifies more than 50,000 robust, symmetry-indicated CTFBs. Filling such CTFBs yields short-range entangled topological states that exhibit power-law correlations. Crucially, all filled CTFB states admit exact tensor-network representations with finite bond dimensions, providing a tractable starting point for exploring strongly correlated topological matter.

The simulation of strongly correlated electron systems remains a formidable challenge. Certain experimentally relevant dynamical response functions are especially difficult to calculate, due to issues of finite-size effects and the ill posed analytic continuation. To address this we propose the quantum cluster embedding theory, an embedded cluster method aimed at computing the response of the system following an external perturbation; the frequency dependent dynamical susceptibility is obtained subsequently by means of inverse linear response theory. The embedded clusters, used within the method as representative of short range correlations, are open quantum systems governed by the Lindblad equation. The short-range correlations extracted from the clusters are used to close the equations of motion for the fermionic bilinear and the local double occupancy on the lattice. In turn, the clusters' Markovian baths are tuned to keep the bilinear and the double occupancy expectation values on the clusters and the lattice identical, throughout the concomitant evolution of the two sets of equations. The theory becomes numerically exact in the non-interacting, atomic and infinite cluster size limits, and it respects the total charge and energy conservation laws. We show that our approach can treat very large lattices while avoiding analytic continuation through the explicit time evolution. Finally we compute the charge-charge correlation function in the square lattice Hubbard model and compare with a recent cold atom experiment, finding good qualitative agreement.

Trilayered Bi-2223 superconductor features the highest critical temperature $T_c$ among the bismuth-based cuprate collection and symbolizes an ideal prototype for studying intrinsic superconducting properties. The previous solid-state reaction method substantiated the growth of the high-quality Bi-2223 compounds but was accompanied by excessively laborious time and effort in terms of multiple grinding, pressing, as well as calcining stages, %causing risk of constituent loss, so finding a less tedious synthesis path is imperative. Here, we present an advanced sol-gel synthesis for assembling the multicomponent complexity of Bi1.4Pb0.6Sr2Ca2(Cu1-xLix)3O10 superconductors (Li-doped Bi-2223), with $x$ = 0.0--0.20, utilizing metallic cationic molecular mixing within the chemical Pechini polyesterization route followed by single-step pyrolysis and sintering stages. Although monovalent cations such as Li$^+$ pose limitations in establishing a perplex crosslinking network or chelating mechanism in the Pechini method, they represent a unique probe to elucidate the major chemical process during polymerization. We observe that a 5 molar~\% Li-doped sample pronounces the highest $T_c$ = 111.4 K among the series of samples, as confirmed by both ac susceptibility and dc resistivity measurements, and is equivalent to the value obtained by our preceding solid state fabrication. In addition, we showcase a rare observation of layer-by-layer crystalline phase growth under microstructure probe. Through analyzing the reliable ac susceptibility data at low magnetic fields in a wide range of frequency, we provide the quantum flux formation and flux creep mechanism by Anderson-Müller's model and Cole-Cole plot.

Chromium antimonide has emerged as a key material platform for studying altermagnetism because of its simple binary composition, high Néel temperature, and semimetallic electronic structure. Here, we investigate electrical and thermal magnetotransport in single-crystalline CrSb using steady-and pulsed-magnetic fields up to 65 T, and complement these measurements with neutron diffraction and magnetization data. We confirm the compensated magnetic structure and observe a large nonsaturating magnetoresistance together with a pronounced nonlinear Hall response at low temperatures. Multicarrier modeling, supported by mobility-spectrum analysis, reveals coexisting electron- and hole-like charge carriers with mobilities up to ~3000 cm2/Vs and shows that the number of transport channels that can be resolved strongly depends on the accessible magnetic-field range. Thermal-transport measurements further reveal a nonlinear thermal Hall response and a thermal conductivity substantially exceeding a simple Wiedemann-Franz law. The broadly similar field and temperature evolution of electrical and thermal transport point to a dominant electronic contribution, while the remaining deviations indicate additional heat-carrying channels.

Luttinger's theorem constrains the particle density of interacting fermions through global properties of the single-particle Green's function, and its violation signals a breakdown of the identification between the quantized Hall response and the Green-function-based Ishikawa-Matsuyama invariant. This phenomenon becomes especially compelling in strongly correlated topological phases, such as fractional Chern insulators, where fractionalized quasiparticles lack an adiabatic connection to electrons, raising the question of how Green's-function-based topological invariants manifest in such phases. Using exact diagonalization of the fermionic Harper-Hofstadter-Hubbard model, we compute bulk single-particle Green's functions deep inside a fractional Chern insulating phase and directly evaluate the Luttinger count, its possible correction (the Luttinger integral), and the Ishikawa-Matsuyama invariant $N_3[\mathrm{G}]$. We demonstrate a clear violation of Luttinger's theorem and show that the fractional nature of the many-body Chern number is encoded in the Středa response of the Luttinger integral, while the integer invariant $N_3[\mathrm{G}]$ arises from the Středa response of the Luttinger count. We also analytically prove that $N_3[\mathrm{G}]$ is fully determined by the Luttinger count together with the Chern number of the occupied Bloch band, upon neglecting Bloch-band mixing. Finally, we propose an experimental protocol to extract all Green-function-based topological invariants from local density-of-states measurements, experimentally accessible in fractional quantum Hall systems.

This paper shows how data-driven machine learning approaches can improve growth control, reproducibility, and physical insight in the pulsed laser deposition (PLD) growth of correlated oxides. Despite well-known relationships between growth conditions and material properties, consistently producing high-quality films of complex materials like LaVO$_3$ remains difficult due to the highly non-equilibrium nature of PLD and the defects and competing phases that accumulate during growth. Here, we use an active learning framework based on Gaussian process Bayesian optimization that incorporates measured bulk and surface lattice properties along with impurity phase information to efficiently map the multidimensional growth space of LaVO$_3$ by PLD. By tuning the relative weighting of these properties, the model identifies an optimized region where phase-pure films of LaVO$_3$ exhibit two-dimensional surfaces, near-ideal lattice parameters, and minimal sub-band gap optical absorption. The trained model reveals clear competition among different defect formation mechanisms that are connected to unseen parameters like supersaturation and surface mobility, thus giving insight into the highly non-equilibrium process of PLD growth. Together, this demonstrates that property-guided machine learning can accelerate materials optimization while providing a new way to address fundamental growth mechanisms in PLD that enable understanding and utilization of quantum phenomena found in complex oxides.

A correspondence exists between the multi-orbital Hubbard model and the bilayer Hubbard model, in which superconductivity is optimized in an incipient-band regime in both cases. In the multi-orbital system, the orbital level offset $ΔE$ plays a role analogous to the interlayer hopping in bilayer systems, and superconductivity is enhanced for large $ΔE$. We refer to such a multi-orbital model as an orbital-space bilayer model (OSBM). In this study, we theoretically propose that a reduced bilayer nickelate La$_3$Ni$_2$O$_6$ can be a candidate for a superconductor described by OSBM when an appropriate amount of holes is doped. By constructing a tight-binding model based on first-principles calculations, a large $ΔE$ between the Ni $d_{x^2-y^2}$ and the other $d$ orbitals is obtained due to the absence of outer apical oxygens. Furthermore, our fluctuation exchange approximation calculations indicate the emergence of $s\pm$-wave superconductivity driven by interorbital interactions in an incipient-band situation, where the superconducting gap function changes its sign between the $d_{x^2-y^2}$ and other $d$ orbital bands. We also investigate the energetic and dynamical stability of the crystal structure under atomic substitution and pressure. Although La$_3$Ni$_2$O$_7$ and La$_3$Ni$_2$O$_6$ share a similar chemical formula, our study shows that an entirely different pairing mechanism can take place in the latter.

Quantum spin liquids are commonly thought to be highly sensitive to lattice geometry, symmetry, and microscopic exchange patterns, leading to a proliferation of seemingly distinct phases across frustrated magnets. Here, we provide a framework that unifies phases that appear distinct from the viewpoint of this intuition. We postulate that the spin-$\tfrac{1}{2}$ Heisenberg antiferromagnets on the square, Shastry-Sutherland, and checkerboard lattices can realize a single unified quantum phase: a gapless $\mathbb{Z}_2$ Dirac quantum spin liquid, despite their markedly different lattice symmetries. Using a systematic projective symmetry group analysis, we identify a checkerboard spin-liquid state that completes a closed set of adiabatically connected phases linking the well-established square-lattice and Shastry-Sutherland spin liquids. Crucially, we show that this lattice-level unification is mirrored exactly in the continuum description. In all three cases, the spin liquids descend from a common SU(2) $π$-flux parent state and are governed by the same gauge theory, QED$_3$ with two Dirac fermion flavors coupled to two adjoint Higgs fields. As a result, we postulate that the surrounding Néel and valence-bond-solid phases and their confinement transitions admit a unified interpretation within the framework of deconfined quantum criticality. More broadly, our results suggest that quantum spin liquids are most fundamentally classified not by lattice geometry or microscopic couplings, but by the emergent gauge theory and its Higgs structure. Distinct frustrated lattices can thus host the same quantum phase and exhibit the same confinement mechanisms, despite substantial differences in their microscopic symmetries.

Quantum magnets in two dimensions display strong quantum interaction effects even when magnetically ordered. Using the metal-organic framework material CuF$_2$(D$_2$O)$_2$(pyz), we investigate the field-dependent spin dynamics of the $S = 1/2$ square-lattice Heisenberg antiferromagnet by high-resolution inelastic neutron scattering to applied fields beyond one third of saturation. We discover an anomalously sharp, dispersive ``shadow mode'' residing within the two-magnon continuum, which shadows the dispersion of the transverse one-magnon branches across the Brillouin zone at an offset equal to the Larmor energy. We perform cylinder matrix-product-state (MPS) calculations that reproduce the field-induced spectrum quantitatively and apply a spectrally consistent $1/S$ spin-wave theory to deduce that the ``Larmor-shadow mode'' is a composite two-magnon resonance: a dispersing magnon at wavevector ${\bf Q}$ couples to the uniform Larmor precession at $Γ$, its small intrinsic linewidth indicating a non-perturbative effect of attractive magnon-magnon interactions. Another quantum-fluctuation phenomenon, the zero-field $(π,0)$ anomaly, is lost at increasing fields, which tighten the spectral weight into the one-magnon and Larmor-shadow modes. To our knowledge, these results constitute the first observation of a sharp quasi-bound state embedded in the continuum of a gapless two-dimensional antiferromagnet.

We benchmark constrained-path Monte Carlo (CPMC) on the triangular-lattice Hubbard model for several fillings and $U$ values and show that symmetry-adapted trial wave functions are essential for quantitative accuracy. Away from half-filling, simple free-electron-based trials that preserve the ground state symmetry yield energy deviations $\lesssim 1\%$ from exact diagonalization and density matrix renormalization group results. At half-filling, strong frustration in the intermediate to large $U$ regimes necessitates symmetry-projected trials to reach comparable accuracy, where both free-electron and symmetry-broken Hartree-Fock trials incur substantial constraint bias. Since the computational cost of CPMC with symmetry projection scales polynomially with system size, our results motivate its use as a practical route for studying competing ground states in strongly correlated, frustrated systems.

We introduce a functional renormalization group framework formulated directly in the Floquet steady-state that systematically incorporates frequency-dependent interaction effects. By retaining the frequency structure of the two-particle vertex up to second order in interaction strength, our approach provides controlled access to dynamical response functions and nonequilibrium transport in driven, interacting systems. Using the periodically driven single-impurity Anderson model as a paradigmatic example, we benchmark our results against state-of-the-art Floquet Green's function methods and find quantitative agreement for finite-frequency observables up to intermediate interaction strengths. Remarkably, we also show that static properties are often captured reliably by much simpler approximations, suggesting practical pathways for modeling driven quantum materials. Finally, we demonstrate that although periodic driving of the dot strongly broadens the Kondo resonance through inelastic scattering, it leaves the many-body Kondo cloud largely intact. This robustness suppresses Floquet replicas of the Kondo peak and leads to a partial persistence of Kondo pinning, highlighting the resilience of emergent many-body correlations under local periodic driving.

Identifying materials hosting an excitonic insulator ground state has been one of the major pursuits in condensed matter physics in recent years. Promising candidates in transition metal chalcogenide compounds (TMC), including $1T-\mathrm{TiSe_2}$, $\mathrm{Ta_2Pd_3Te_5}$, and $\mathrm{Ta_2NiSe_5}$, share a crucial common characteristic: their low-energy physics is governed by electrons in $d-$ orbitals subject to strong on-site Coulomb interactions. In this work, we investigate spatially indirect exciton condensation in two-dimensional semimetals on triangular lattice. Using a combination of dynamical mean-field theory and the determinant quantum Monte Carlo method, we study two- and three-orbital Hubbard models incorporating strong on-site ($U$) and inter-orbital interactions ($V$). Our results demonstrate that on-site Hubbard $U$ can strongly suppress the condensation temperature $T_c$, an effect that is particularly pronounced at higher electron-hole pair densities. This behavior contrasts sharply with the case without on-site $U$, where $T_c$ grows with pair density at fixed $V$. Moreover, we uncover competition among multiple electron-hole pairing channels in the three-orbital model, which also acts to suppress $T_c$ of exciton condensation. An orbital-selective electron-hole pairing state is identified. These findings may help explain the large discrepancy between strong binding-energy and relative low transition temperature for indirect excitons in TMCs materials, offering important insights for understanding and engineering exciton condensation in materials with strongly correlated $d-$ shell electrons.

Long-range interactions can qualitatively reorganize correlated-electron ground states. In the square-lattice Hubbard model, on-site repulsion produces antiferromagnetic spin and charge stripes upon doping. We show that including a nearest-neighbor repulsion $V$ can dramatically alter this behavior. Using auxiliary-field quantum Monte Carlo and density matrix renormalization group methods, we find that, above a critical ratio $V/U$ ($\sim 0.25$), the system develops a modulated ferrimagnetic order intertwined with checkerboard charge-density-wave. Inside the ferrimagnetic domains, spin density alternates between positive (or negative) and nearly zero values. When the total spin is fixed to zero, positive and negative domains alternate in space; when spins are unconstrained, a ferrimagnetic state emerges with finite magnetization. Including a next-nearest-neighbor hopping $t'$ changes the modulation wavelength but leaves the order robust. Our results demonstrate that even short-range nonlocal interactions can stabilize qualitatively new magnetic textures, with implications for cuprate materials and programmable quantum simulators.

We use first-principles calculations to investigate the crystal and electronic structure of the hybrid bilayer-trilayer Ruddlesden-Popper (RP) nickelate La$_7$Ni$_5$O$_{17}$ under hydrostatic pressure and biaxial compressive strain. By analyzing the irreducible representations of the dynamically unstable phonon modes in the high-symmetry $P4/mmm$ structure, we identify a dynamically stable lower-symmetry $C2/c$ structure containing octahedral tilts. The application of both pressure and compressive strain tends to suppress the octahedral tilts, effectively tetragonalizing the structure, in analogy with the conventional RPs. The electronic structure under hydrostatic pressure and strain has similarities, but it differs in the position of the $d_{z^2}$ bonding band from the trilayer block. This band crosses the Fermi level at a pressure of 30 GPa, but it remains below it for any level of compressive strain. This strain-induced modification mirrors the electronic structure changes observed in the conventional bilayer nickelate.

Accurate exchange-correlation (XC) potentials are essential for density functional theory, yet reliable approximations remain challenging for strongly correlated systems. In this work, we present a quantum algorithmic framework to determine spin-resolved XC potentials using a variational quantum eigensolver. Using the Hubbard model as a prototypical strongly correlated lattice system, we prepare ground states in fixed spin sectors through a Hamiltonian variational ansatz combined with a continuation strategy that gradually increases the interaction strength. From the resulting many-body ground states, we extract the XC energy and compute the corresponding spin-resolved XC potentials via finite differences. The accuracy of the approach is benchmarked against exact diagonalization for one- and two-dimensional Hubbard systems of various lattice sizes. We demonstrate that the variational ansatz reproduces the ground-state energies and densities with high fidelity, enabling accurate construction of both magnetic and non-magnetic XC potentials. We analyzed the dependence of the XC potentials on the interaction strength, charge, spin densities, and magnetization. We also present an empirical complexity scaling relation for the computational cost of the method at a fixed fidelity. These results illustrate how quantum simulations can be used to construct spin-resolved XC functionals for correlated lattice models, providing a potential pathway for improving density functional approximations in strongly correlated materials.

We have developed a semi-analytical framework formulated in the canonical fermion representation to investigate strongly correlated electron systems. We consider the U=$\infty$ Hubbard model and used the equation of motion method to calculate the fermion self-energy which has two parts: single and two-boson exchange processes. The emergent bosons here are self-generated local charge and spin-density fluctuations which become strongly time-dependent due to extreme correlations. The computed boson spectral density is a diffusive damped mode with a long tail. The electron self-energy at $d=\infty$ is computed self-consistently. The corresponding fermionic spectral density displays a pronounced coherence peak at $ω=0$, while its frequency derivative develops a two-peak structure at finite $ω$. The resistivity shows a linear temperature dependence over a broad range, crossing over to coherent Fermi-liquid behavior at extremely low temperatures.

The kagome-lattice Hubbard model attracts widespread interest due to its flat-band and Van Hove singularity features, which can give rise to unconventional magnetism. We employ determinant quantum Monte Carlo simulations to systematically investigate the uniform magnetic susceptibility across a range of on-site interactions and electron fillings on a two-dimensional kagome lattice. Beyond the Van Hove singularity, dominant ferromagnetic fluctuations emerge. Magnetic susceptibility grows markedly with increasing interaction strength and decreasing temperature, indicating that the Van Hove singularity acts as a critical point for the crossover of dominant magnetic fluctuations. Finite-size analysis further suggests the potential stabilization of a finite-temperature ferromagnetic phase. We also examine the sign problem to identify numerically reliable parameter regimes. These results provide valuable insights into controlling magnetic fluctuations in kagome systems and establish a computational framework for exploring flat-band physics in regimes characterized by novel quantum phases and competing orders.

Because the normal state of underdoped cuprate superconductors is an enigmatic Fermi-arc metal, it is valuable to analyze an exactly solvable model that exhibits both Fermi arcs and $d$-wave superconductivity. Here, we focus on a recently proposed solvable model in which the emergence of Fermi arcs is especially transparent. Upon incorporating a $d$-wave pairing interaction, the model produces an asymptotically exact solution for the superconducting transition temperature $T_c$ that traces out a superconductivity dome as a function of hole doping, in qualitative agreement with experimental observations in cuprates. Crucially, we show analytically that the Fermi arcs generate an additional many-body effect that suppresses $T_c$ beyond the simple reduction expected from a shrinking Fermi surface. The many-body nature of the Fermi arcs further introduces the gap-to-$T_c$ ratio greatly surpassing the mean-field limit. These findings provide an analytic benchmark for understanding how Fermi-arc physics competes with $d$-wave superconductivity in high-$T_c$ superconductors.

A promising application of ML is in creating low-cost surrogate models to replace computational bottlenecks in quantum many-body simulations. Here, we explore whether a NN can be trained in the low-data regime, with one to two orders of magnitude fewer training examples than previous works, as an efficient substitute for the impurity solver in DMFT simulations of correlated electron models. We show that the NN solver achieves accuracy comparable to popular CTQMC impurity solvers when interpolating between samples within the training set. While the NN's performance decreases notably when extrapolating to lower temperatures outside the training distribution, its output still provides an excellent initial guess for input to more accurate CTQMC impurity solvers, thus accelerating the time to solution up to a factor of five. We discuss our results in the context of rapid phase-space exploration.

Despite extensive research, the precise spin Hamiltonian of the van der Waals antiferromagnet NiPS$_3$ -- which hosts a zigzag-ordered ground state -- remains debated. While consensus has emerged on ferromagnetic nearest-neighbor ($J_1$) and antiferromagnetic third-nearest-neighbor ($J_3$) Heisenberg interactions, recent studies suggest a biquadratic ($B$) exchange term may also play a role, though its estimated magnitude varies widely. To address this controversy, we perform density functional theory calculations and extract a positive biquadratic interaction with $B/J_3 \approx 0.44$. Within the minimal $J_1$-$J_3$-$B$ model, we show that these parameters naturally stabilize zigzag ordering using minimally augmented spin-wave theory. Density-matrix renormalization group calculations further validate our extracted parameters as a reasonable description of the ground state. Although fully resolving the spin Hamiltonian of NiPS$_3$ requires further investigation, our findings provide new insights into its biquadratic interaction.

We study the phase diagram of interacting spinless fermions on the honeycomb bilayer at charge neutrality using large-scale quantum Monte Carlo simulations. In the noninteracting limit, the low-energy spectrum features quadratically dispersing bands that touch at the corners of the hexagonal Brillouin zone. Weak to intermediate interactions induce a splitting of each of the quadratic band touching points into four Dirac points, located along high-symmetry directions of the reciprocal lattice. Strong interactions lead to the formation of a layer-polarized charge density wave, which spontaneously breaks the $\mathbb Z_2$ layer inversion symmetry and opens an insulating gap in the spectrum. We show that the semimetal-to-insulator quantum phase transition as a function of interaction is continuous and characterized by emergent relativistic symmetry. Our results for the values of the correlation-length exponent $ν$, the order-parameter anomalous dimension $η_φ$, and the fermion anomalous dimension $η_ψ$ agree with those of the theoretically predicted 2+1D Gross-Neveu-Ising universality class with eight two-component Dirac fermions within less than 5\%\ deviation. We also determine the crossover scale as a function of interaction strength between the nonrelativistic semimetal state at high temperatures, characterized by dynamical critical exponent $z = 2$, and the Dirac semimetal state at intermediates temperatures, characterized by $z=1$. Further reducing the temperature below the crossover scale at a fixed value of the interaction strength above the quantum critical point results in a classical ordering transition in the 2D Ising universality class.

Strongly correlated fermionic systems are of great interest in condensed matter physics and numerical methods are indispensable tools for their study. However, existing approaches such as exact diagonalization (ED) and stochastic quantum Monte Carlo methods each suffer from fundamental limitations: ED is hindered by exponential scaling in system size, while Monte Carlo methods are plagued by sign problems and long autocorrelation times. These limitations restrict the accessible parameter space and developing algorithms that efficiently alleviate them remains a central challenge in computational physics. In this work, we propose a hybrid algorithm that combines ED and Hamiltonian Monte Carlo (HMC) to simulate 2D arrays of coupled quantum wires, modeled as interacting fermionic Hubbard chains. We demonstrate how our hybrid implementation of HMC, which we dub H$^2$MC, outperforms either method alone across several key simulation facets. When compared to pure ED, H$^2$MC has a much more favorable computational scaling, which allows us to push simulations to much larger 2D arrays. H$^2$MC also greatly alleviates the sign problem and reduces autocorrelation times when compared to pure HMC formulations utilizing either real or imaginary auxiliary fields. Our formalism demonstrates how complementary strengths of seemingly disparate methods can be leveraged to enable feasible simulations in an extended parameter space.

Continuous control over lattice geometry, when combined with long-range interactions, offers a powerful yet underexplored tool to generate highly frustrated quantum spin systems. By considering long-range dipolar antiferromagnetic interactions on a breathed Kagome lattice, we demonstrate how these tools can be leveraged to stabilize a chiral spin liquid. We support this prediction with large-scale density-matrix renormalization group calculations and explore the surrounding phase diagram, identifying a route to adiabatic preparation via a locally varying magnetic field. At the same time, we identify the relevant low-energy degrees of freedom in each unit cell, providing a complementary language to study the chiral spin liquid. Finally, we carefully analyze its stability and signatures in finite-sized clusters, proposing direct, experimentally viable measurements of the chiral edge mode in both Rydberg atom and ultracold polar molecule arrays.

For the Fermi-Hubbard model in the Mott insulator phase, we employ the hierarchy of correlations to study how doublon and holon quasi-particle excitations are affected by adding disorder to the system. We study two types of disorder: charge disorder, in the form of on-site potential randomness; and spin disorder, in the form of a fixed, randomly generated background spin arrangement. By analysing the quasi-particle eigen-spectra and quantifying the degree to which the corresponding eigen-states localise, we find both an energetic and spatial separation between localised and delocalised states in the charge disorder. In contrast, the spin disorder results in localised states throughout the quasi-particle bands. Finally, we repeat our calculations using strong-coupling perturbation theory, and compare the results obtained from both methods.

Motivated by recent experimental progress on high-temperature superconductivity in bilayer nickelates, we investigate the phase diagram of the normal state in a bilayer Kondo lattice model using single-site dynamical mean-field theory (DMFT). When the interlayer tunneling $t_\perp$ is absent, we identify a non-Fermi-liquid (NFL) critical point tuned by the interlayer spin coupling $J_\perp$ or hole doping $x$, which separates a standard Fermi liquid in the overdoped region from a distinct pseudogap (PG) metal in the underdoped regime. This PG phase, which we term the `second Fermi liquid' (sFL), exhibits small hole pockets and violates the perturbative Luttinger theorem despite the absence of symmetry breaking or fractionalization. The PG metal behaves like a heavy Fermi liquid, with small quasi-particle residue and large effective mass. We also provide an intuitive analytical description of the pseudogap and the ground-state wave function based on an ancilla-fermion framework. Inside the PG phase, we interpret the ancilla fermion as a spin-polaron and demonstrate a Kondo-like resonance peak in the spectral function of this composite fermion directly in DMFT calculation. Extending the analysis to finite $t_\perp$, we apply this framework to the bilayer nickelate $\mathrm{La}_3\mathrm{Ni}_2\mathrm{O}_7$. We propose that current experimental samples ($x \approx 0.5$) reside in the overdoped FL regime, suggesting that the pseudogap phase and the NFL criticality may be accessed via electron doping.

We study the hole-doped Kitaev-Heisenberg chain using the density-matrix renormalization group. In the Kitaev-only limit, the bond-directional exchange itself promotes pairing, favoring spin-singlet and spin-triplet superconducting tendencies for antiferromagnetic and ferromagnetic Kitaev couplings, respectively, together with finite-size Majorana edge correlations suggestive of topological superconductivity. In the full Kitaev-Heisenberg chain, cooperative $J$ and $K$ exchanges broadly stabilize superconductivity, while competition between $J$ and $K$ induces a strong filling dependence and enables superconductivity even when both $J$ and $K$ are weak. At quarter filling, this competition produces a Mott insulator with spontaneous hopping dimerization. These results identify magnetic frustration as a common mechanism underlying superconducting and interaction-driven insulating phases in doped Kitaev systems.

We develop a highly efficient framework for computing the thermal entropy in the doped Fermi-Hubbard model within the grand-canonical ensemble. The framework comprises four calculation schemes that express the entropy as path integrals in the parameter space of temperature, interaction strength, and chemical potential. The integrands involve only fundamental observables, including the total energy, fermion density, and double occupancy, which are readily accessible in a wide range of theoretical and numerical methods. We further derive useful Maxwell relations connecting the entropy to other quantities, and present practical formulas for directly evaluating the grand potential. As an application, we compute the entropy of the doped Hubbard model in two and three dimensions, using the numerically unbiased auxiliary-field quantum Monte Carlo method. The test results show excellent agreement across the different schemes and quantitatively verify the Maxwell relations, confirming the reliability of the framework. In two dimensions, we further benchmark our entropy results in physically relevant parameter regimes against diagrammatic Monte Carlo calculations and observe excellent quantitative consistency between the two approaches. By providing an efficient and broadly applicable route for entropy evaluation, our work facilitates the thermodynamic characterization of complex correlated states in the doped Hubbard model.

We construct a Hubbard model with a nearly flat band whose quantum geometry can be tuned independently of the energy dispersion and the Coulomb interaction. We show that, when the nearly flat band is half-filled, the exact ground state of the model exhibits ferromagnetism and that this ferromagnetism is stabilized by the quantum metric through the spin stiffness. Furthermore, we demonstrate that tuning the quantum geometry alone drives a magnetic phase transition. Our nonperturbative results without resorting to mean-field approximations reveal the quantum-geometric origin of ferromagnetism and the underlying many-body physics in dispersive-band systems.

We present the x-ray magnetic circular dichroism (XMCD) at the Ni $L_{2,3}$-edge as an evidence of the $d$-wave altermagnetism in rutile-structure NiF$_2$. Sizable XMCD signal is observed in excellent agreement with theoretical simulations. Owing to a considerable net magnetization due to spin canting, the XMCD spectrum consists of an altermagnetic signal as well as a non-negligible ferromagnetic contribution. We verify experimentally that the XMCD spectrum can be written as a sum of contributions from altermagnetism and weak ferromagnetism. Two experimental methods to isolate the ferromagnetic contribution are shown to yield essentially the same result. These are dependence of XMCD on applied magnetic fields below the Néel temperature and the XMCD measured in applied field above the Néel temperature. Our results demonstrate the utility of XMCD as a probe for altermagnetic materials with the coexisting weak ferromagnetism induced by the relativistic spin-orbit coupling.

Novel superconducting phases have been found in various moiré heterostructures based on hexagonal lattices. However, the archetypal high-temperature superconductors (cuprates, iron-based and nickelate families) all share a square lattice foundation. These materials host a rich landscape of correlated phenomena, such as charge and spin stripes, pseudogap behavior, and unconventional metallicity, which continue to challenge our fundamental understanding of strongly correlated electrons. In this work, we investigate the possibility of simulating the effective models governing these high-$T_c$ superconductors using twisted homobilayers of $Γ$-valley square-lattice systems. We develop a universal theoretical framework and carry out a detailed analysis of a promising candidate material ZnF$_2$. We find that the first moiré band realizes a single-orbital square-lattice Hubbard model, widely believed to capture cuprate physics, while the second and third moiré bands map to a $p_x,p_y$ two-orbital square-lattice Hubbard model, which shares common physics to the minimal $d_{xz}, d_{yz}$ models proposed for iron pnictides. Our study combines continuum Hamiltonian modeling, first-principle calculations, and Hartree-Fock mean field theory. The latter focuses on the quarter-filling regime of the two-orbital model and in particular leads to, among others, a stable antiferro-orbital, ferromagnetic insulating phase. These results highlight $Γ$-valley square-lattice moiré systems as a new and important generation of van der Waals heterostructures to realize interesting strongly correlated phases of matter.

We present an efficient protocol leveraging classical computation to support Initial State Preparation for strongly correlated fermionic systems, a critical bottleneck for fault-tolerant quantum simulation. Focusing on nuclear shell model eigenstates, we first demonstrate that the Density Matrix Renormalization Group algorithm can efficiently approximate target states as Matrix Product States, capitalizing on the favourable entanglement structure of these fermionic systems. These high-fidelity approximations are then leveraged as a classical resource in a variational circuit optimization scheme to compile shallow quantum circuits. We establish concrete resource estimates by decomposing the resulting circuits into the industry-standard Clifford$+T$ gateset, exploring the benefits of specialized $U3$ synthesis techniques. For all nuclear systems tested, on up to 76 qubit Hamiltonians, we consistently find low $T$-count circuits preparing the nuclear eigenstates to high fidelity with $\sim 2\times 10^4$ total $T$ gates. This low number gives confidence these eigenstates can be prepared on early fault-tolerant quantum computers. Our work establishes a viable path toward practical ground state preparation for nuclear structure and other fermionic applications.

Spin-orbit coupling (SOC) gives rise to complex magnetic states such as spin liquids, skyrmion crystals, and topological spin-wave excitations. We consider exchange interactions in multi-orbital Mott insulators where SOC is strong on ligand ions. SOC on the ligands enables electron hopping accompanied by spin flips and fluctuations in the orbital state of the ligand hole. These processes generate anisotropic exchange interactions and greatly increase the number of possible exchange paths. The number grows further with the inclusion of hopping between ligands, which mediates interactions between more distant spins. We propose an effective method to calculate exchange interactions at arbitrary separations between spins. Applying it to monolayer CrI$_3$, we obtain anisotropic interactions between nearest-neighbor and next-nearest-neighbor Cr spins, as well as single-ion anisotropy induced by long-range hopping. In this material, magnetic anisotropy stabilizes long-range ferromagnetic order and opens a magnon gap at the Dirac points, which defines a nontrivial magnon band topology. Using Hubbard model parameters from first-principles calculations, the resulting spectrum agrees well with the spin-wave dispersion observed experimentally in bulk CrI$_3$, except that the calculated Dirac gap is much smaller.

We introduce Bi$_2$CuO$_3$(SO$_4$) as a rare example of a spin-ladder magnet with ferromagnetic interactions on the rungs. Its magnetic response is studied through measurements of heat capacity, temperature-dependent magnetic susceptibility, and field-dependent magnetization, as well as electron spin resonance spectroscopy. These experiments are complemented by density-functional-theory calculations combined with the construction of maximally localized Wannier functions and an analysis of the relevant superexchange pathways. Quantum Monte Carlo simulations are employed to model thermodynamic properties and to quantitatively determine the magnetic exchange parameters. Our combined approach identifies Bi$_2$CuO$_3$(SO$_4$) as a two-leg spin-ladder system with ferromagnetic rungs ($J'$ $\approx -208$ K) and antiferromagnetic legs ($J$ $\approx 258$ K). These interactions of similar magnitude arise from remarkably different superexchange pathways, with the Cu--Cu distance along the leg being almost twice as long than the respective distance along the rung. The antiferromagnetic leg coupling represents the strongest oxygen-mediated long-range superexchange in a Cu$^{2+}$ compound reported to date and sets the benchmark for the role of complex superexchange pathways in quantum magnets.

We investigated the pressure dependence of the superconducting transition temperature ($T_{\rm c}$) in optimally doped Bi$_2$Sr$_2$Ca$_2$Cu$_3$O$_{10+δ}$ (Bi-2223) single crystals using different pressure-transmitting media. Previous high-pressure studies have reported conflicting behaviors, ranging from a resurgence of $T_{\rm c}$ of optimally doped Bi-2223 in fluid media to an insulating-like transition in solid media. However, a direct comparison of the effects of different pressure-transmitting media is lacking. Here, we employed wide-field quantum microscopy based on nitrogen-vacancy centers to probe the magnetic response under high pressure, utilizing cBN and KBr as media. We observed that a diamagnetic response near 70 K, indicative of the superconducting transition, persisted up to 23 GPa in KBr, whereas it disappeared above 11 GPa and 70 K in cBN. These results demonstrate the high sensitivity of Bi-2223 to the pressure environment and highlight the critical role of hydrostatic pressure in cuprate superconductors.

Ultracold Bose gases in one-dimensional optical lattices constitute an important benchmark problem in the study of strongly interacting many-body quantum phases. Here we present a combined experimental and theoretical study of their phase-coherence properties over a wide range of lattice depths. Experimentally, we extract the single-particle correlation function directly from the measured momentum distribution. Theoretically, we perform tensor-network simulations of the Bose-Hubbard model that incorporate all relevant experimental parameters. For deep lattices well within the Mott insulator regime, the experimental results are in good agreement with the expected zero-temperature behavior, with only small temperature-dependent corrections. As the lattice depth is reduced, finite-temperature effects become increasingly important. We find that the experimental data are quantitatively described by an effective temperature extracted from the tensor-network simulations, and that this effective temperature decreases markedly with increasing lattice depth. Rather than indicating actual cooling, we interpret this behavior as evidence of inhibition of thermalization caused by the formation of Mott domains that suppress heat transport. Counterintuitively, the inhibition of thermalization favors the preparation of an effectively low-entropy quantum gas in the trap center for large lattice depths.

Building upon the spin-1/2 Kitaev model on a honeycomb lattice, the Yao-Lee spin-orbital model provides exactly solvable quantum spin liquids with potentially better stability against perturbations due to the additional degree of freedom. Recently, the microscopic mechanism underlying the Yao-Lee interaction in honeycomb materials has been uncovered, leading to an extended Kitaev-Yao-Lee spin-orbital model when the celebrated Kugel-Khomskii interaction is included. Numerical studies of this model have identified various disordered phases, including a broad region of the nematic phase that is reminiscent of a spin-orbital liquid. Here, we investigate the origin and stability of this nematic phase via thermal and quantum fluctuations using classical Monte Carlo simulations and a generalized spin wave theory appropriate for the spin-orbital model. We demonstrate that the additional spin-orbital degree of freedom gives rise to strong thermal and quantum fluctuations in spin-orbital models, providing insight into the emergence of disordered phases.

The low-energy spin- and charge-sensitive thermodynamic properties of a broad range of strongly correlated 4f-electron systems follow Kondo scaling, with a characteristic Kondo temperature, $T_K$. While the theory is known for thermodynamic properties and high-energy spectroscopies of Kondo materials, the surface sensitivity of electron spectroscopy limits the extent to which Kondo scaling can be quantitatively verified. In this study, bulk-sensitive photon-in photon-out temperature-dependent resonant inelastic X-ray scattering (RIXS), in combination with single-impurity Anderson model (SIAM) calculations, is used to provide quantitative evidence of low- and high-energy Kondo scaling in CeSi$_2$. RIXS Ce M$_5$-edge spectra show a clear decrease in the occupancy of the $f^0$ state as temperature increases accompanied by an increase of the spectral weight of the $f^1\underline L^1$ state, in good agreement with the SIAM calculations. The results demonstrate the breakdown of the Kondo singlet state, coupled with thermal occupation of the low-lying first-excited magnetic states. The RIXS data reveal a temperature evolution of the $f^n$ spectral weights, which is in stark contrast to that extracted from photoemission and inverse photoemission spectroscopies. This study provides an accurate spectroscopic method to determine the Kondo energy $k_B$$T_K$ that is consistent with thermodynamic measurements, and highlights soft X-ray RIXS as a quantitative bulk probe of low- and high-energy-scale hybridization effects in strongly correlated materials.

Low-dimensional materials exhibit exotic properties due to enhanced quantum fluctuations, making the understanding of their microscopic origin central in condensed matter physics. Analogue quantum simulators offer a powerful approach for investigating these systems at the microscopic level, particularly in large-scale regimes where quantum entanglement limits classical numerical methods. To date, analogue simulators have largely focused on universal Hamiltonians rather than material-specific quantitative comparisons. Here we use a Rydberg-based quantum simulator to study the bulk-layered frustrated quantum magnet TmMgGaO$_4$. Magnetisation measurements obtained from the quantum simulator show excellent agreement with independent measurements performed in a magnetic laboratory facility, validating the proposed effective two-dimensional microscopic Hamiltonian. Building on this quantitative correspondence, we investigate on both platforms the antiferromagnetic phase transition. We further probe the role of quantum fluctuations via snapshot analysis, connecting our results to integrated inelastic neutron scattering data. Finally, we access, with the simulator, non-equilibrium dynamics on picosecond material timescales, including frequency response and thermalisation of observables.

Moiré superlattices have emerged as a premier platform for simulating the Hubbard model, yet achieving high tunability in square-lattice systems remains a key challenge. We demonstrate that $Γ$-valley twisted square homobilayers provide a faithful and highly tunable realization of $t-t'-U$ Hubbard model, extending the recent proposal in M-valley systems. We show that at small twist angles, an emergent layer-exchange symmetry decouples electronic states into flat bands residing on two nested square sublattices. An interlayer displacement field breaks this symmetry to induce controllable inter-sublattice hybridization, enabling wide-range experimental tuning of the effective hopping ratio $t'/t$. By establishing a direct correspondence between $Γ$- and M-valley systems, we provide a unified framework for understanding displacement-field tunability in square moiré physics. These findings establish $Γ$-valley twisted bilayers as a versatile platform for simulating the square-lattice Hubbard model and exploring its rich landscape of correlated phenomena.

Following the Bethe ansatz we determine the dynamical spectra of the one-dimensional supersymmetric t-J model. A series of fractionalized excitations are identified through two sets of Bethe numbers. Typical patterns in each set are found to yield wavefunctions containing elementary spin and charge carriers, manifested as distinct boundaries of the collective excitations in the spectra of single electron Green functions. In spin channels, gapless excitations fractionalized into two spin and a pair of postive and negative charge carriers, extending to finite energy as multiple continua. These patterns connect to the half-filling limit where only fractionalized spinons survive. In particle density channel, apart from spin-charge fractionalization, excitations involving only charge fluctuations are observed. Furthermore, nontrivial Bethe strings encoding bound state structure appear in channels of reducing or conserving magnetization, where spin and charge constituents can also be identified. These string states contribute significantly even to the low-energy sector in the limit of vanishing magnetization.

The recent discovery of altermagnetism has demonstrated that spin-split electronic band structures can emerge in magnetic systems with zero net magnetization. In contrast, fully compensated ferrimagnetic (fFIM) systems remain far less explored, despite exhibiting similar characteristics such as vanishing magnetization and spin-split bands. Here, based on first-principles calculations combined with theoretical analysis, we demonstrate that monolayer CoS and CoSe can be driven into fFIM states by an external electric field. These materials possess collinear antiferromagnetic ground states with out-of-plane Néel vectors, and their electronic bands are spin degenerate due to $\mathcal{PT}$ symmetry. When an out-of-plane electric field is applied, $\mathcal{PT}$ symmetry is broken, inducing fFIM states with pronounced spin splitting. Moreover, we show that the resulting fFIM states host fully spin-polarized currents, anomalous Hall effects, and magneto-optical Kerr and Faraday effects. Our results establish monolayer CoS and CoSe as promising platforms for electric-field-controlled fFIM states and spintronic applications.

We report emergent magnetic behavior in heterostructures composed of (111)-oriented La$_{0.7}$Sr$_{0.3}$MnO$_3$ (LSMO) and (00$l$)-oriented Bi$_2$Te$_3$ (BT), controlled by interfacial reconstructions. When BT is deposited directly onto LSMO, an intermediate interfacial layer forms between the two materials. Polarized Neutron Reflectometry modeling reveals that this reconstructed region stabilizes a secondary magnetically ordered phase that is coupled to the underlying ferromagnetic LSMO layer. As a consequence, the heterostructures exhibit unconventional self-crossing magnetic hysteresis loops at room temperature, characterized by a reversal of the net magnetization at low applied fields. In contrast, the introduction of a tellurium seed layer results in a sharper LSMO-BT interface, while preserving the anomalous hysteresis behavior and enhancing the saturation magnetization. Element-specific X-ray absorption spectroscopy suggests that the emergent magnetic phase originates from the chemical reconstruction of manganese species. These results demonstrate that interface engineering in magnetic oxide-topological insulator heterostructures provides a pathway to control emergent magnetic coupling and emergent magnetic states in oxide-topological insulator heterostructures.