arXiv Curator

We present the combination of a complex-time tensor-network impurity solver with an analytic continuation scheme based on exponential fitting as an efficient framework for single and multi-orbital dynamical mean-field calculations. By performing time-evolution along a complex-time contour, the approach balances computational cost with the difficulty of spectral recovery, offering greater flexibility than methods confined to the real or imaginary axis. By complementing the complex-time evolution with an exponential fitting scheme, we faithfully extract real-time information at negligible cost. The resulting method obtains high-resolution spectra at a significantly lower computational cost than real-time evolution, offering a promising tool for ab initio studies of strongly correlated materials.

Cuprate superconductors remain central to condensed matter physics due to their technological relevance and unconventional, incompletely understood electronic behavior. While the canonical phase diagram and low-energy models have been shaped largely by studies of underdoped and moderately doped cuprates, the overdoped regime has received comparatively limited attention.Here, we track the evolution of the electronic structure from optimal to heavy overdoping in La2-xSrxCuO4(LSCO) using broadband optical spectroscopy across x=0.15-0.60. The measured spectral changes--including the redistribution of Zhang-Rice-related spectral weigh--are in qualitative agreement with determinant quantum Monte Carlo simulations of the three-orbital Emery model, which together indicate a pronounced reconstruction of the electronic structure beyond hole concentrations x>0.2. Guided by these observations, we propose a spontaneous checkerboard-type Zhang-Rice electronic configuration that captures the coexistence of itinerant and localized carriers characteristic of the heavily overdoped state. Our results refine the doping-dependent Zhang-Rice-based framework for cuprates, illuminate how correlations persist deep into the overdoped regime, and provide new constraints on microscopic mechanisms of high-temperature superconductivity, with broader implications for correlated transition-metal oxides.

The Kondo effect, a hallmark of many-body physics, emerges from the antiferromagnetic coupling between localized spins and conduction fermions, leading to a correlated many-body singlet state. Here we propose to use the mixed-dimensional (mixD) bilayer Hubbard geometry as a platform to study Kondo lattice physics with current ultracold atom experiments. At experimentally feasible temperatures, we predict that key features of the Kondo effect can be observed, including formation of the Kondo cloud around a single impurity and the competition of singlet formation with Ruderman-Kittel-Kasuya-Yosida (RKKY) interactions for multiple impurities, summarized in the Doniach phase diagram. Moreover, we show that the mixD platform provides a natural bridge between the Doniach phase diagram of the Kondo lattice model, relevant to heavy-fermion materials, and the phase diagram of cuprate superconductors as described by a single-layer Zhang-Rice type $t$-$J$ model: It is possible to continuously tune between the two regimes by changing the interlayer Kondo coupling. Our findings demonstrate that the direct connection between high-temperature superconductivity and heavy-fermion physics can be experimentally studied using currently available quantum simulation platforms.

Reliable constraints on iron's melting curve at Earth's inner-core boundary require accurate finite-temperature electronic correlations, yet DFT+DMFT calculations remain too costly for large-scale thermodynamic sampling. Here, we develop a machine-learning accelerator for charge self-consistent DFT+DMFT by training E(3)-equivariant graph neural networks to predict the local self-energy and Fermi level from atomic environments, providing an efficient warm start to the DMFT self-consistency loop. Using high-throughput data for Fe, FeO, and NiO, we obtain a 2-4 times reuduction in DMFT iterations. Leveraging this improvement, we generate correlated energies and forces for Fe at core pressures, train a neural-network interatomic potential, and determine the melting curve via two-phase coexistence simulations. We obtain a predicted melting temperature of 6225 K at 330 GPa.

In this review, we first present compelling evidence from resonant inelastic x-ray scattering data that highlights the significance of the long-range Coulomb interaction in cuprate charge dynamics, particularly around the in-plane momentum q=(0,0). We show that these experimental observations are well-captured by the layered t-J-V model, which extends the standard t-J framework to include the long-range Coulomb interaction V and the layered structure. This new perspective elucidates how charge dynamics renormalizes one-particle excitation properties, leading to several profound and often counterintuitive consequences. We demonstrate that the electron dispersion does not exhibit a sharp kink, and Landau quasiparticles persist in the low-energy limit despite a significant suppression of their spectral weight. We further show that while charge fluctuations alone cannot fully account for the pseudogap, they are a crucial component for understanding its formation. Additionally, we reveal that optical plasmon excitations generate fermionic quasiparticles, known as plasmarons, which give rise to a distinct, incoherent replica band. We argue that accurately describing these plasmonic effects requires a three-dimensional theoretical approach. This perspective on plasmon excitations may offer a critically new clue to a long-standing puzzle: why multi-layer cuprate superconductors, containing more than two CuO2 layers per unit cell, consistently exhibit a higher critical temperature Tc than their single-layer counterparts. Finally, we review the spin-fluctuation mechanism of superconductivity suffers from the "self-restraint effect" and show how important the screened Coulomb interaction is in the spin-fluctuation mechanism to realize high-Tc superconductivity.

Using the DFT+dynamical mean-field theory approach we study the effects of electronic correlations and doping on the normal state electronic structure of the double-layer nickelate superconductor La$_3$Ni$_2$O$_7$ under pressure. In agreement with experiments, we obtain significant orbital-dependent quasiparticle renormalizations of the Ni $x^2-y^2$ and $3z^2-r^2$ bands, accompanied by incoherence (bad metal behavior) of the $3z^2-r^2$ states, caused by the proximity of the Ni $3d$ states to orbital-dependent localization. Our results demonstrate a sensitive, non-monotonic dependence of $m^*/m$ on doping, with a remarkable, by about 20\%, increase for the Ni $x^2-y^2$ orbitals upon electron doping $x \sim 0.2$ (per Ni ion), implying a significant enhancement of orbital-dependent correlations with oxygen deficiency in LNO. We observe a reconstruction of the low-energy electronic structure of LNO upon doping above $x\sim -0.3$ and 0.2. It is associated with the Lifshitz transition, with a crossover to a self-doping regime characterized by partial occupation of the La $5d$ bands (upon an electron doping $x>0.2$). Our analysis of the static magnetic susceptibility $χ({\bf q})$ obtained within DFT+DMFT suggests the possible formation of the spin and charge (or bond) density wave stripes, implying strong spin and charge correlations in LNO. We show that this behavor is associated with suppression of the Néel $G$-type antiferromagnetic ordering of the Ni$^{2+}$ ions upon hole doping. Interestingly, upon a moderate electron doping of the Ni$^{2.5+}$ ions (e.g., with oxygen deficiency), we find a significant enhancement of the strength of in-plane spin and charge fluctuations. We note a close resembles of our results to those for the bilayer Hubbard model, which shows the boosting of superconductivity as one of the two electron bands approaches the Lifshitz transition (e.g., upon doping).

Using the recently proposed bond-operator technique (BOT), we discuss spin dynamics of the Heisenberg spin-$\frac12$ antiferromagnet with the ring exchange and small interactions between the second- and the third-neighbor spins on the square lattice at $T=0$. This model was suggested before for description of parent compounds of high-temperature superconducting layered cuprates. BOT describes accurately short-range spin correlations in quantum systems and provides a quantitative description of elementary excitations which appear in other approaches as bound states of conventional low-energy quasiparticles. We demonstrate that besides well-known magnons (spin-1 excitations) there are three well-defined spin-0 quasiparticles in the considered model whose energies lie near the magnon spectrum. Two of them, the amplitude (Higgs) mode and the quasiparticle which we named singlon, produce pronounced anomalies observed experimentally in the Raman scattering, resonant inelastic x-ray scattering, and infrared optical absorption. We find sets of the model parameters which describe quantitatively experimental data obtained in $\rm La_2CuO_4$ and $\rm Sr_2CuO_2Cl_2$.

We review recent theoretical progress on the charge dynamics of doped carriers in high-temperature cuprate superconductors. Advances in this field have clarified that doped charges in cuprates exhibit remarkably rich collective behavior, governed by the combined effects of strong electronic correlations, the intrinsic layered crystal structure, and long-range Coulomb interaction. First, the emergence of acousticlike plasmons has been firmly established through quantitative analyses of resonant inelastic x-ray scattering (RIXS) spectra based on the t-J-V model -- an extension of the conventional t-J model that incorporates the layered crystal structure and the long-range Coulomb interaction V. These acousticlike plasmons arise near the in-plane momentum q=(0,0) and possess characteristic energies far below the well-known ~ 1 eV optical plasmon. This behavior is found to be universal across both hole- and electron-doped cuprates, including multilayer systems. Second, in electron-doped cuprates, a pronounced tendency toward d-wave bond-charge order develops near q=(0.5pi, 0), as revealed by resonant x-ray scattering (RXS) and RIXS. As a result, the charge dynamics acquires a dual structure, in which low-energy bond-charge excitations coexist with relatively high-energy plasmons. Third, analogous signatures of charge-order tendency have also been reported in hole-doped cuprates. However, a direct application of the d-wave bond-charge-order framework fails to account for experimental observations. Similarly, the charge-stripe order in La-based cuprates remains unresolved within existing theoretical approaches. Assuming that mobile carriers behave in a largely universal manner across electron- and hole-doped systems, we discuss a possible scenario that may reconcile these diverse experimental findings.

Understanding the magnetic ground state of Ruddlesden-Popper nickelates is crucial, as these materials exhibit superconductivity under high pressure and host competing electronic orders that may play a key role in the pairing mechanism. In this work, we investigate the magnetic properties of the alternating monolayer-trilayer phase of La$_3$Ni$_2$O$_7$ (1313-La$_3$Ni$_2$O$_7$) using muon-spin rotation/relaxation ($μ$SR) under both ambient and hydrostatic pressure conditions. The monolayer-trilayer phase develops incommensurate magnetic order below approximately 150 K, with a mean ordering temperature of $T_{SDW} \simeq 123$ K and a transition width of $ΔT_{SDW} \simeq 15$ K. The abrupt onset of the internal magnetic field indicates a first-order-like transition. Hydrostatic pressure ($p$) suppresses the magnetic ordering temperature at a rate of $dT_{SDW}/d p \simeq -3.9$ K/GPa, demonstrating a progressive destabilization of the ordered state. By comparison with the bilayer 2222-La$_3$Ni$_2$O$_7$ and the trilayer 3333-La$_4$Ni$_3$O$_{10}$ systems, and within a unified phenomenological framework, systematic trends are identified linking the pressure dependence of $T_{SDW}$, the (in)commensurability of the magnetic order, and the character of the magnetic transition. These trends consistently indicate a gradual reduction of electronic correlation strength from the bilayer to the monolayer-trilayer and trilayer nickelates. This hierarchy suggests that the higher superconducting transition temperature observed in the 2222 phase may be closely connected to its more strongly correlated electronic nature. These results position the alternating monolayer-trilayer 1313-La$_3$Ni$_2$O$_7$ as an intermediate member linking the magnetic behavior of the bilayer 2222-La$_3$Ni$_2$O$_7$ and the trilayer 3333-La$_4$Ni$_3$O$_{10}$ Ruddlesden-Popper compounds.

The discovery of high-temperature superconductivity in La$_3$Ni$_2$O$_7$ under high pressure has sparked a surge of research into Ruddlesden-Popper (RP) nickelates. Currently, stabilizing the bilayer RP phases with smaller $A$-site ions remains a significant challenge. In this work, we have successfully synthesized medium- and high-entropy bilayer nickelates, La$_{1.2}$Pr$_{0.6}$Nd$_{0.6}$Sm$_{0.6}$Ni$_2$O$_{7-δ}$ and La$_{0.67}$Pr$_{0.67}$Nd$_{0.67}$Sm$_{0.33}$Eu$_{0.33}$Gd$_{0.33}$Ni$_2$O$_{7-δ}$, by utilizing the concept of configuration entropy stabilization. The high-entropy nickelate exhibits the smallest unit-cell volume and the largest orthorhombic distortion reported to date. The chemical pressure induced by the smaller A-site ions significantly enhances the NiO$_6$ octahedral rotation/distortion and shortens the interlayer Ni-Ni interatomic spacing. Physical property measurements reveal bad electrical conductivity alongside a markedly elevated density-wave transition temperature. Notably, the superconducting transition temperature extrapolated from structural correlations is projected to exceed 100 K. Our work not only demonstrates entropy stabilization of bilayer nickelates, but also reveals the effect of $A$-site-ion size on the crystal structure and physical properties, opening a new pathway for developing nickelate superconductors and tuning their electronic properties.

Altermagnetism simultaneously possesses nonrelativistic spin responses and zero net magnetization, thus combining advantages of ferromagnetism and antiferromagnetism. This superiority originates from its unique dual feature, i.e., opposite-magnetic sublattices in real space and alternating spin polarization in momentum space enforced by the same crystal symmetry. Therefore, the determination of an altermagnetic order and its unique spin response inherently necessitates atomic-scale spin-resolved measurements in real and momentum spaces, an experimental milestone yet to be achieved. Here, via utilizing the helical edge (hinge) modes of a higher order topological insulator as the spin sensor, we realize spin-resolved scanning tunneling microscopy which enables us to pin down the dual-space feature of a layered $d$-wave altermagnet, KV$_2$Se$_2$O. In real space, atomic-registered mapping demonstrates the checkerboard antiferromagnetic order together with density-wave lattice modulation, and in momentum space, spin-resolved spectroscopic imaging provides a direct visualization of d-wave spin splitting of the band structure. Critically, using this new topology-guaranteed spin filter we directly reveal the unidirectional, spin-polarized quasiparticle excitations originating from the crystal symmetry-paired X and Y valleys around opposite magnetic sublattices simultaneously --the unique spin response for $d$-wave altermagnetism. Our experiments establish a solid basis for the exploration and utilization of altermagnetism in layered materials and further facilitate access to atomic-scale spin sensing and manipulating of 2D quantum materials.

Although the Kadanoff-Baym equations are typically solved using time-stepping methods, iterative global-in-time solvers offer potential algorithmic advantages, particularly when combined with compressed representations of two-time objects. We examine the computational complexity and stability of several global-in-time iterative methods, including multiple variants of fixed point iteration, Jacobian-free methods, and a Newton-Krylov method using automatic differentiation. We consider the ramped and periodically-driven Falicov-Kimball and Hubbard models within time-dependent dynamical mean-field theory. Although we observe that several iterative methods yield stable convergence at large propagation times, a standard forward fixed point iteration does not. We find that the number of iterations required to converge to a given accuracy with a fixed time step size scales roughly linearly with the number of time steps. This scaling is associated with the formation of a propagating front in the residual error, whose velocity is method-dependent. We identify key challenges which must be addressed in order to make global solvers competitive with time-stepping methods.

We present a diffuse neutron and x-ray scattering study of structural, spin- and charge-density-wave fluctuations in the electrical insulator La$_{2-x}$Sr$_x$NiO$_4$. This lamellar nickelate is an isostructural analogue of the high-temperature cuprate superconductor La$_{2-x}$Sr$_x$CuO$_4$, for which recent experiments uncovered evidence for unusual structural and superconducting fluctuations indicative of rare-region physics due to inherent inhomogeneity unrelated to common point disorder effects. We find closely analogous nanoscale orthorhombic fluctuation behavior in La$_{2-x}$Sr$_x$NiO$_4$, including exponential scaling of the diffuse scattering intensity and power-law scaling of the characteristic length with relative temperature. Moreover, our neutron and x-ray scattering data reveal similar behavior for short-range magnetic and charge fluctuations above the respective ordering temperatures. These observations indicate that rare-region effects are a generic feature of perovskite-related structures and lead to universal fluctuations of both structural and electronic degrees of freedom over extended temperature ranges.

The discovery of high-temperature superconductivity in bilayer nickelate La$_{3}$Ni$_{2}$O$_{7}$ necessitates a minimal theoretical model that unifies the superconducting phase with the spin-density-wave (SDW) phase without external pressure or strain. We propose a model where half-filled $d_{z^{2}}$ local moments interact with itinerant $d_{x^{2}-y^{2}}$ electrons via strong Hund's coupling $J_H$, which reduces to a bilayer type-II t-J model in the large $J_H$ limit. Using iDMRG calculations on an $L_y=4, L_z=2$ cylinder, we demonstrate that the competition between double-exchange ferromagnetism and in-plane superexchange generates period-4 stripe-like SDW order-a feature absent in one-orbital t-J model with only $d_{x^2-y^2}$ orbital. Furthermore, increasing the interlayer exchange coupling suppresses magnetic order and stabilizes interlayer s-wave superconductivity. These results identify the type-II t-J model as a minimal framework for capturing the interplay of magnetism and superconductivity in bilayer nickelates.

The discovery of high-T$_c$ superconductivity in La$_3$Ni$_2$O$_7$ has opened the door to a new route to high temperature superconductivity, distinct from that in cuprates and iron-based materials. Yet, despite intense recent activity, we lack experimentally testable protocols for distinguishing between different pairing scenarios. In this Letter, we construct a minimal two-band model that reproduces the Fermi-surface topology observed in recent ARPES measurements and DFT calculations, and we analyze superconductivity arising from two distinct pairing mechanisms. We show that these mechanisms yield sharply different responses to an applied perpendicular electric field. Thus, La$_3$Ni$_2$O$_7$ offers the unique opportunity to cleanly distinguish between different pairing scenarios. Finally, we propose three concrete experimental proposals designed to distinguish these scenarios and thereby identify the pairing mechanism most relevant to the real material.

The recent experimental observations of high temperature superconductivity in bilayer nickelate have attracted lots of attentions. Previous studies have assumed a mirror symmetry $\mathcal M$ between the two layers and focused on uniform and clean superconducting states. Here, we show that breaking this mirror symmetry via an applied displacement field can stabilize a pair-density-wave (PDW) superconductor, which is similar to the Fulde--Ferrell--Larkin--Ovchinnikov (FFLO) state, but at zero magnetic field. Based on a mean-field analysis of a model of $d_{x^2-y^2}$ orbital with an effective inter-layer attraction, we demonstrate that the PDW phase is robust over a wide range of displacement field, interlayer hopping strengths, and electron fillings. Finally, we analyze disorder effects on interlayer superconductivity within the first Born approximation. Based on symmetry considerations, we show that pairing is weaken by disorders which break the mirror symmetry, even with unbroken time reversal symmetry. Our results establish bilayer nickelate as a tunable platform for realizing finite-momentum pairing and for exploring generalized disorder effects.

To understand the intricate exchange between electrons of different bands in strongly correlated materials, it is essential to treat multi-orbital models accurately. For this purpose, dynamical mean-field theory (DMFT) provides an established framework, whose scope crucially hinges on the availability of efficient quantum impurity solvers. Here we present a real-frequency impurity solver based on neural quantum states (NQS) combined with an operator-Lanczos construction. NQS are an asymptotically unbiased variational ground-state ansatz that employs neural networks to capture long-range correlations on complicated graph structures. We leverage this ability to solve multi-orbital impurity problems using a systematically improvable Segmented Commutator Operator-Lanczos (SCOL) construction. Our benchmarks on both the single-orbital Anderson model and the multi-orbital Hubbard-Kanamori impurity Hamiltonian reveal excellent ground-state precision and the capacity to accurately resolve zero temperature spectral functions and self-energies. These results open avenues for extending DMFT to more challenging problems.

The recent discovery of high-temperature superconductivity in pressurized nickelates has renewed interest in the broken-symmetry states of their ambient-pressure parent phases, where a density-wave (DW) order emerges and competes with superconductivity, but its microscopic origin remains unresolved. Using ultrafast optical spectroscopy, we track quasiparticle relaxation dynamics across the DW transition at $T_{\rm DW} \approx$ 136 K in trilayer nickelate La$_4$Ni$_3$O$_{10}$ single crystals, revealing the opening of an energy gap of $\sim$ 52 meV. Multiple coherent phonons, including $A_g$ modes near 3.88, 5.28, and 2.09 THz, display pronounced mode-selective anomalies across the transition, demonstrating that the DW is coupled with lattice degree of freedom stabilized through electron-phonon coupling. At higher excitation densities, the DW is nonthermally suppressed, producing a temperature-fluence phase diagram that parallels pressure-tuned behavior. These results establish the DW in La$_4$Ni$_3$O$_{10}$ as a lattice-entangled instability involving multiple phonon modes, and highlight ultrafast optical excitation as a powerful tool to manipulate competing orders in nickelates.

We study anisotropic superconductivity in the nearly quarter-filled $α$-$\mathcal{T}_3$ lattice. We analyze an extended Hubbard model with off-site attractive interactions within the mean-field framework and find two distinct chiral $d+id'$-wave superconducting phases characterized by different Chern numbers. We further investigate the superconducting mechanism mediated by spin fluctuations arising from purely repulsive interactions by applying the fluctuation-exchange (FLEX) approximation to the Hubbard model. The gap symmetry obtained by solving the linearized Eliashberg equation is $d$-wave, which corresponds to a $d+id'$-wave superconducting state with a Chern number of $8$, including the spin degree of freedom. The $\mathbf{q}=\mathbf{0}$ antiferromagnetic spin fluctuation, which possesses the largest spectral weight at finite energies arising from the incipient flat band, gives rise to an effective spin-singlet pairing glue between rim sites.

We develop a microscopic theory for pseudospin-$\frac{1}{2}$ screw-chain compounds with spin-orbit coupling that goes beyond the phenomenological site-dependent $g$-tensor description traditionally used for XXZ-like BaCo$_2$V$_2$O$_8$ and related materials. Starting from the symmetry-allowed $JKΓ$ Hamiltonian with Heisenberg $J$, Kitaev $K$, and off-diagonal $Γ$ interactions, we show that the $Γ$ interaction naturally generates the four-sublattice pattern associated with the crystal's screw symmetry. Using the density matrix renormalization group, we identify two distinct field-induced transitions. The first is a continuous transition into an intermediate phase, where the symmetry responsible for the two-fold ground-state degeneracy is broken. The second is a first-order transition into the high-field phase, characterized by a discontinuous jump in the spin-spin correlator. Entanglement-entropy scaling confirms that the first transition belongs to the Ising critical point with the central charge $1/2$. These results establish a microscopic framework for pseudospin-$\frac{1}{2}$ screw-chain systems such as Co$^{2+}$ materials, uncover an intermediate phase whose width increases with $Γ$, and provide guidance for systematic exploration of additional field orientations and structural distortions.

Ruddlesden-Popper (RP) nickelates have attracted intense interest following the discovery of superconductivity in several members of the series, including bilayer (BL) La$_3$Ni$_2$O$_7$, trilayer (TL) La$_4$Ni$_3$O$_{10}$, and structural polymorphs composed of monolayer-bilayer or monolayer-trilayer (ML-TL) units. However, an inherent propensity of the RP series to form intergrown phases during single-crystal synthesis, together with spatial variations in oxygen stoichiometry, has complicated the determination of their intrinsic material properties. As a consequence, conflicting reports have emerged on both their electronic phase transitions and lattice dynamics. In this work, we perform a comparative study of the phononic and electronic Raman responses of high-quality ML-TL single crystals and contrast them with those of other RP nickelates, using samples with optimized oxygen content. We establish several Raman spectral features that enable unambiguous phase identification across the series. Moreover, we uncover characteristics in the phononic and electronic Raman response of ML-TL that are not reflected in the pure ML and TL compounds. We attribute these differences to a distinctive electronic structure arising from self-doping and confinement effects induced by the ML unit within the ML-TL lattice architecture.

Understanding superconductivity emerging from repulsive fermions remains a major challenge in condensed matter physics. In this paper, we investigate the pairing tendencies in a one-dimensional, three component repulsive Hubbard model, using the density matrix renormalization group method. At half-filling, the system exhibits density wave ground state due to strong Hubbard repulsions. Upon doping, we find that Cooper pairs can emerge, whose fluctuations predominate the long-range physics in the system across a wide parameter range. The effective attractions between Cooper pairs are mediated by the particle-hole fluctuations in the third non-pairing component, resembling an excitonic mechanism of superconductivity. The coexistence of multiple density waves and superconductivity at different fermion fillings is explored. We also present an analytical study of the pairing mechanism in both weak and strong coupling limits. Our results provide a new perspective for understanding and exploring unconventional superconductivities in strongly correlated fermionic systems.

The insulator-quasi-metal (bad metal) transition observed in Jahn-Teller (JT) magnets orthonickelates RNiO$_3$ (R = rare earth, or yttrium Y) is considered a canonical example of the Mott transition, traditionally described in the framework of Hubbard's $U-t$ model. However, in reality, the insulating phase of nickelates is the result of charge disproportionation (CD) with the formation of a system of spin-triplet ($S = 1$) electron [NiO$_6$]$^{10-}$ and spinless ($S = 0$) hole [NiO$_6$]$^{8-}$ centers, equivalent to a system of effective spin-triplet composite bosons moving in a nonmagnetic lattice. The effective CD-phase Hamiltonian takes into account local ($U$) and nonlocal ($V$) correlations, and the transfer of composite bosons ($t_b$). Within the framework of the effective field approximation, we have shown the existence of two types of CD phases: the high-temperature classical paramagnetic CO-phase of charge ordering of electron and hole centers, and the low-temperature magnetic quantum CDq phase with charge and spin density transfer between electron and hole centers, with ''uncertain valence'' [NiO$_{6}$]$^{(9\pmδ)-}$ ($0 \le δ\le 1$) and spin density $(1 \pm δ)/2$ NiO$_6$-centers. In the classical CO phase, spin-triplet electron centers are surrounded by the nearest nonmagnetic hole centers, which ''turns off'' the strong superexchange interaction of the nearest neighbors. The magnetic ordering in the quantum CDq phase is determined by a strong traditional superexchange and an unusual bosonic double exchange mechanism.

We study the quantum phase diagram of the spin-$1/2$ nearest-neighbor Heisenberg model on the maple-leaf lattice using infinite projected entangled pair states (iPEPS) combined with a corner transfer matrix renormalization group scheme adapted to $C_3$-symmetric lattices. Focusing on the fully antiferromagnetic $J$-$J_d$ model with $J_h = J_t := J$, we map out the ground-state phase diagram as a function of the dimer coupling $J_d$. Our results show that the system hosts only two phases: a magnetically ordered canted-$120^\circ$ phase and an exact dimer singlet product phase. We identify a first-order transition between these two phases at $J_d/J \approx 1.45$. Within the magnetically ordered phase, we observe small but finite magnetic moments. We also resolve the quantum renormalization of the canting angle, which deviates from the classical prediction over almost the entire magnetically ordered phase.

The functional renormalization group (fRG) is an established tool in the treatment of correlated electron systems, notably for the description of competing instabilities. In recent years, methodological advancements led to the multiloop extension of the fRG, which systematically includes loop corrections beyond the conventional one-loop truncation and yields a quantitatively accurate description of two-dimensional lattice systems. At the same time, the single-boson exchange (SBE) decomposition of the two-particle vertex has been shown to offer both computational and interpretative advantages paving the way to more affordable approximation schemes. We here apply their combination coined as multiloop SBE fRG to the two-dimensional Hubbard model at weak coupling. After providing a detailed account of the underlying formalism in physical channels, we analyze the results for the frequency- and momentum-dependent vertex functions. We find that the SBE approximation, i.e., without calculating explicitly multi-boson exchange contributions, accurately reproduces the parquet approximation at loop convergence. The presented algorithmic improvement opens the route for the treatment of more challenging parameter regimes and more realistic models.

Recent numerical and theoretical studies on the two-dimensional $J$-$Q_3$ model suggests that the deconfined quantum critical point is actually a $SO(5)$-symmetry-enhanced first-order phase transition that is spontaneously broken to $O(4)$. However, this conclusion has mainly relied on finite-size scaling of the entanglement entropy, lacking direct evidence from physical observables.} Here, we investigate the dynamical spectra of spin and bond operators at the deconfined critical point of the $J$-$Q_3$ model using large-scale quantum Monte Carlo simulations, and contrasting them with the well-established $\mathrm{O(3)}$ Wilson-Fisher criticality in the $J_1$-$J_2$ Heisenberg model. Although both models exhibit two gapless magnon modes in the Néel phase, their critical behaviors diverge strikingly. At the $J_1$-$J_2$ critical point, the Higgs mode becomes gapless, yielding three gapless modes that reflect the full restoration of the $\mathrm{O(3)}$ symmetry. {In the $J$-$Q_3$ model, we instead observe four gapless transverse modes at the either side of the transition. This spectral feature, together with the entanglement entropy results, provides direct evidence for the weakly first-order scenario that the deconfined quantum critical point exhibits an emergent $\mathrm{SO(5)}$ symmetry that spontaneously breaks to $\mathrm{O(4)}$.

Murunskite K$_2$Cu$_3$FeS$_4$ is a representative sulfosalt, isostructural to the pnictides, but with electronic properties more similar to the insulating parent compounds of the cuprates. We use it as a bridge to compare the chemical and physical roles of metal and ligand orbitals in cuprates and pnictides. In cuprates, ionicity, covalency, and metallicity are tightly interwoven to give rise to high-temperature superconductivity (SC). Their most remarkable property is the interaction of an ionically localized hole on the copper (Cu) with a Fermi liquid (FL) on the oxygens (O), which is critically important for understanding all key properties of these materials. The localization is due to strong correlations on the Cu $3d$ orbital. We describe a scenario in which the localized hole gives rise both to SC by Cooper scattering of O holes, and to Fermi arcs, as observed in cuprate spectroscopy, the latter by a purely kinematic projection of the static local disorder, without invoking any residual interactions between the mobile O FL carriers. In the pnictides, the orbitals responsible for binding and metallic conduction appear to be separate. The Fe $3d$ $e_{g}$ orbitals hybridized with the ligands set the lattice spacing. The $3d$ $t_{2g}$ orbitals overlap directly between the Fe atoms, resulting in several electronic bands appearing at the Fermi level. The ensuing Fermi liquid exhibits both charge and magnetic correlations. We argue that a similar SC scenario as in the cuprates is plausible in the pnictides, except that a light FL scatters on a slow nearly-antiferromagnetic (AF) one, rather than on localized holes as in the cuprates.

The experimental discovery of fractional quantum anomalous Hall (FQAH) states in tunable moiré superlattices has sparked intense interest in exploring the interplay between topological order and symmetry breaking phases. In this paper, we present a comprehensive numerical study of this interplay through large-scale density matrix renormalization group (DMRG) simulations on a minimal two-band lattice model of twisted bilayer MoTe$_2$ at filling $ν=-2/3$. We find robust FQAH ground states and provide clear numerical evidences for anyon excitations with fractional charge and pronounced real-space density modulations, directly supporting the recently proposed anyon density-wave halo picture. We also map out the displacement field dependent phase diagram, uncovering a rich landscape of charge ordered states emerging from the FQAH, including a quantum anomalous Hall crystal (QAHC) with an integer quantized Hall conductance. We expect our work to inspire further research interest of intertwined correlated topological phases in moiré systems.

Two of the iconic phases of the hole-doped cuprate materials are the intermediate temperature pseudogap metal and the lower temperature $d$-wave superconductor. Following the prescient suggestion of P.~W.~Anderson, there were numerous early theories of these phases as doped quantum spin liquids. However, these theories have had difficulties with two prominent observations: (i) angle-dependent magnetoresistance measurements (ADMR) in the pseudogap metal, including observation of the Yamaji effect, present convincing evidence of small hole pockets which can tunnel coherently between square lattice layers (Fang et al., Nature Physics 18, 558 (2022); Chan et al., Nature Physics 21, 1753 (2025)) and (ii) the velocities of the nodal Bogoliubov quasiparticles in the $d$-wave superconductor are highly anisotropic, with $v_F \gg v_Δ$ (Chiao et al., Phys. Rev. B 62, 3554 (2000)). These lecture notes review how the fractionalized Fermi Liquid (FL*) state, which dopes quantum spin liquids with gauge-neutral electron-like quasiparticles, resolves both difficulties. Theories of insulating quantum spin liquids employing fractionalization of the electron spin into bosonic or fermionic partons are discussed. Doping the bosonic parton theory leads to a holon metal theory: while not appropriate for the cuprate pseudogap, this theory is argued to apply to the Lieb lattice. Doping the fermionic parton theory leads to a $d$-wave superconductor with nearly isotropic quasiparticle velocities. The construction of the FL* state in a single band Hubbard-type model is described using a quantum dimer model, followed by a more realistic description using the Ancilla Layer Model (ALM), which is then used to obtain the theory of the pseudogap and the $d$-wave superconductor.

It is well known that the electronic structure of hole-doped cuprate superconductors is tunable through both chemical doping and external pressure, which frequently offer us new insights of understanding on the high-Tc superconducting mechanism. While, for electron-doped cuprate superconductors, although the chemical doping effects have been systematically and thoroughly investigated, there is still a notable lack of experimental evidence regarding the pressure-driven coevolution of Tc and electronic structure. In this study, we report the first observation on the signatures of pressure-induced Lifshitz transition in Pr0.87LaCe0.13CuO4+delta (PLCCO) single crystal, a typical electron-doped cuprate superconductor, through the comprehensive high-pressure measurements of electrical resistance, Hall coefficient (RH) and synchrotron X-ray diffraction (XRD). Our results reveal that, at 40 K, the ambient-pressure RH with a significantly negative value decreases with increasing pressure until it reaches zero at a critical pressure (Pc ~ 10 GPa). Meanwhile, the corresponding Tc exhibits a slight variation within this pressure range. As pressure is further increased beyond Pc, RH changes its sign from negative to positive and then shows a slight increase, while Tc displays a continuous decrease. Our XRD measurements at 40 K demonstrate that no crystal structure phase transition occurs across the Pc. These results reveal that applying pressure to PLCCO can induce a Lifshitz transition at Pc, manifesting the reconstruction of the Fermi surface (FS), which turns the superconductivity toward fading out. Our calculation further reinforces the Fermi surface reconstruction from electron-dominated to hole-dominated ones at around Pc. These findings provide new evidence that highlights the strong correlation between the superconductivity and the Fermi surface topology in the electron-doped cuprates.

We investigate the stabilization of topological Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phases, with a specific emphasis on the intraband FFLO phase, in a one-dimensional (1D) Fermi gas subjected to an external magnetic field. This research highlights the crucial role of the interplay between Rashba spin-orbit coupling (RSOC) and Dresselhaus spin-orbit coupling (DSOC). Employing a Fermi-Hubbard model alongside the density matrix renormalization group (DMRG) method, we examine the combined effects of RSOC and DSOC on these exotic superfluid phases, taking into account attractive fermionic interactions. Our principal finding reveals that while RSOC primarily stabilizes conventional zero-momentum pairing, DSOC performs a distinct and crucial role in selectively stabilizing the intraband FFLO phase. This stabilization is achieved by enhancing spin polarization within a single helicity band and suppressing interband coherence, thereby facilitating the formation of finite-momentum FFLO pairs within the same band and resulting in the emergence of a topologically nontrivial superfluid. This targeted control of intraband FFLO pairing paves the way for new strategies in the manipulation of superfluid phases in spin-orbit coupled systems and offers essential insights for experimental realizations in ultracold atomic gases, with implications for topological quantum computing and Majorana fermions.

Two of the iconic phases of the hole-doped cuprate materials are the intermediate temperature pseudogap metal and the lower temperature $d$-wave superconductor. Following the prescient suggestion of P.W. Anderson, there were numerous early theories of these phases as doped quantum spin liquids. However, these theories have had difficulties with two prominent observations: (i) angle-dependent magnetoresistance measurements (ADMR), including observation of the Yamaji effect, present convincing evidence of small hole pockets which can tunnel coherently between square lattice layers, and (ii) the velocities of the nodal Bogoliubov quasiparticles in the $d$-wave superconductor are highly anisotropic, with $v_F \gg v_Δ$. These lecture notes review how the fractionalized Fermi Liquid (FL*) state, which dopes quantum spin liquids with gauge-neutral electron-like quasiparticles, resolves both difficulties. Theories of insulating quantum spin liquids employing fractionalization of the electron spin into bosonic or fermionic partons are discussed. Doping the bosonic parton theory leads to a holon metal theory: while not appropriate for the cuprate pseudogap, this theory is argued to apply to the Lieb lattice. Doping the fermionic parton theory leads to a $d$-wave superconductor with nearly isotropic quasiparticle velocities. The construction of the FL* state is described using a quantum dimer model, followed by a more realistic description using the Ancilla Layer Model (ALM), which is then used to obtain the theory of the pseudogap and the $d$-wave superconductor.

Hematite $α$-Fe$_2$O$_3$ is a $g$-wave altermagnetic material, which has an easy-axis phase and easy-plane weak ferromagnetic phase below and above Morin transition temperature, respectively. The presence of these phases renders it a good candidate to study the characteristic spin splitting in altermagnets under the impacts of relativistic effect and finite temperature. In this regard, we have calculated the band structure of $α$-Fe$_2$O$_3$ based on density functional theory (DFT) which also considers the Hubbard-U correction and spin-orbit coupling (SOC) effects. Additionally, the charge self-consistent DFT + dynamical mean-field theory (DMFT) calculations have been performed at finite temperatures. It is found that the altermagnetic spin splitting in $α$-Fe$_2$O$_3$ preserves with either SOC or temperature effect taken into account. Furthermore, we present a numerical simulation of the x-ray magnetic circular dichroism (XMCD) of the L$_{2,3}$ edge of Fe using a combination of DFT with multiplet ligand-field theory (MLFT). In terms of the different Néel vectors present in $α$-Fe$_2$O$_3$, we calculate the x-ray absorption spectroscopy (XAS) of the L$_{2,3}$ edge of Fe in the form of conductivity tensor and analyze the XMCD response from a perspective of symmetry. A characteristic XMCD line shape is expected when the Néel vector is along [010] direction (magnetic point group $2^\prime/m^\prime$) and the light propagation vector is perpendicular to the Néel vector, which can be further distinguished from the XMCD response originated from weak ferromagnetism with the light propagation vector parallel to the Néel vector.

We introduce a method for constructing global approximations to correlation functions of strongly interacting quantum field theories, starting from perturbative results. The key idea is to employ interpolation method, such as the two-point Padé expansion, to interpolate the weak and strong coupling expansions of correlation function. We benchmark this many-body interpolation approach on two prototypical models: the lattice $φ^4$ field theory and the 2D Hubbard model. For the $φ^4$ theory, the resulting two point Padé approximants exhibit uniform and global convergence to the exact correlation function. For the Hubbard model, we show that even at second order, the Padé appproximant already provides reasonable characterization of the Matsubara Green's function for a wide range of parameters. Finally, we offer a heuristic explanation for these convergence properties based on analytic function theory.

Band flattening has been identified as key ingredient to correlation phenomena in Moiré materials and beyond. Here, we examine strongly repulsive fermions on a ladder -- a minimal platform for unconventional $d$-wave pairing -- and show that flattening of the lower band through an additional diagonal hopping term produces non-Fermi liquid behavior, evidenced by the violation of Luttinger's theorem, as well as axial $d$-wave pairing correlations. Alternatively, plaquette ring exchange can also generate pairing, albeit with a distinct diagonal $d$-wave pairing symmetry. Hence, our finding showcases a competition of different unconventional pairing channels, and demonstrates via a simple model how band geometry can induce fermionic pairing. This offers broadly relevant insights for correlated flat-band systems, ranging from ultracold atoms to strongly interacting electrons in solids.

Reliably simulating two-dimensional many-body quantum dynamics with projected entangled pair states (PEPS) has long been a difficult challenge. In this work, we overcome this barrier for low-energy quantum dynamics by developing a stable and efficient time-dependent variational Monte Carlo (tVMC) framework for PEPS. By analytically removing all gauge redundancies of the PEPS manifold and exploiting tensor locality, we obtain a numerically well-conditioned stochastic reconfiguration (SR) equation amenable to robust solution using the efficient Cholesky decomposition, enabling long-time evolution in previously inaccessible regimes. We demonstrate the power and generality of the method through four representative real-time problems in two dimensions: (I) chiral edge propagation in a free-fermion Chern insulator; (II) fractionalized charge transport in a fractional Chern insulator; (III) vison confinement dynamics in the Higgs phase of a Z2 lattice gauge theory; and (IV) superfluidity and critical velocity in interacting bosons. All simulations are performed on 12x12 or 13x13 lattices with evolution times T = 10 to 12 using modest computational resources (1 to 5 days on a single GPU card). Where exact benchmarks exist (case I), PEPS-tVMC matches free-fermion dynamics with high accuracy up to T = 12. These results establish PEPS-tVMC as a practical and versatile tool for real-time quantum dynamics in two dimensions. The method extends the reach of classical tensor-network simulations for studying elementary excitations in quantum many-body systems and provides a valuable computational counterpart to emerging quantum simulators.

We investigate the electronic structure of AgF2, AgFBF4, AgF and Ag2O using X-ray absorption spectroscopy (XAS) and resonant inelastic X-ray scattering (RIXS) at the Ag L3 edge. XAS results were compared with density functional theory computations of the spectra, allowing an identification of main features and an assessment of the theoretical approximations. Our RIXS measurements reveal that AgF2 exhibits charge transfer excitations and dd excitations, analogous to those observed in La2CuO4. We propose to use the ratio of dd to CT spectral weight as a measure of the covalence of the compounds and provide explicit equations for the weights as a function of the scattering geometry for crystals and powders. The measurements at the metal site L3 edge and previous measurements at the ligand K edge reveal a striking similarity between the fluorides and cuprates materials, with fluorides somewhat more covalent than cuprates. These findings support the hypothesis that silver fluorides are an excellent platform to mimic the physics of cuprates, providing a promising avenue for exploring high-Tc superconductivity and exotic magnetism in quasi-two-dimensional (AgF2) and quasi-one-dimensional (AgFBF4) materials.

We present a symbolic implementation of recursion method for the dynamics of strongly correlated fermions on one-, two- and three-dimensional lattices. Focusing on two paradigmatic models, interacting spinless fermions and the Hubbard model, we first directly confirm that the universal operator growth hypothesis holds for interacting fermionic systems, manifested by the linear growth of Lanczos coefficients. Equipped with symbolically computed Lanczos coefficients and knowledge of their asymptotics, we are able to compute infinite-temperature autocorrelation functions up to times long enough for thermalization to occur. In turn, the knowledge of autocorrelation functions unlocks transport properties. We compute with high precision the charge diffusion constant over a broad range of interaction strengths, $V$. Surprisingly, we observe that these results are well described by a simple functional dependence $\sim 1/V^2$ universally valid both for small and large $V$. All results are obtained directly in the thermodynamic limit. Our results highlight the promise of symbolic computational paradigm where the most costly step is performed once and outputs symbolic results that can further be used multiple times to easily compute physical quantities for specific values of model parameters.

Infinite projected entangled-pair states (iPEPS) provide a powerful variational framework for two-dimensional quantum matter and have been widely used to capture bosonic topological order, including chiral spin liquids. Here we extend this approach to \emph{fermionic} topological order by variationally optimizing $U(1)$-symmetric fermionic iPEPS for a fractional Chern insulator (FCI), with bond dimensions up to $D=9$. We find evidence for a critical bond dimension, above which the ansatz faithfully represents the FCI phase. The FCI state is characterized using bulk observables, including the equal-time single-particle Green's function and the pair-correlation function, as well as the momentum-resolved edge entanglement spectrum. To enable entanglement-spectrum calculations for large iPEPS unit cells, we introduce a compression scheme and show that the low-lying part of the spectrum is already well converged at relatively small cutoff dimensions.

Computing spectral functions in large, non-periodic super-moiré systems remains an open problem due to the exceptionally large system size that must be considered. Here, we establish a tensor network methodology that allows computing momentum-resolved spectral functions of non-interacting and interacting super-moiré systems at an atomistic level. Our methodology relies on encoding an exponentially large tight-binding problem as an auxiliary quantum many-body problem, solved with a many-body kernel polynomial tensor network algorithm combined with a quantum Fourier transform tensor network. We demonstrate the method for one and two-dimensional super-moiré systems, including super-moiré with non-uniform strain, interactions treated at the mean-field level, and quasicrystalline super-moiré patterns. Furthermore, we demonstrate that our methodology allows us to compute momentum-resolved spectral functions restricted to selected regions of a super-moiré, enabling direct imaging of position-dependent electronic structure and minigaps in super-moiré systems with non-uniform strain. Our results establish a powerful methodology to compute momentum-resolved spectral functions in exceptionally large super-moiré systems, providing a tool to directly model scanning twisting microscope tunneling experiments in twisted van der Waals heterostructures.

Spin-charge stripes belong to the most prominent low-temperature orders besides superconductivity in high-temperature superconductors. This phase is particularly challenging to study numerically due to finite-size effects. By investigating the formation of long, isolated stripes, we offer a perspective complementary to typical finite-doping phase diagrams. We use the density-matrix renormalization group algorithm to extract the ground states of an 18-leg cylindrical strip geometry, making the diameter significantly wider than in previous works. This approach allows us to map out the range of possible stripe filling fractions on the electron versus hole-doped side. We find good agreement with established results, suggesting that the spread of filling fractions observed in the literature is governed by the physics of a single stripe. Taking a microscopic look at stripe formation, we reveal two separate regimes - a high-filling regime captured by a simplified squeezed-space model and a low-filling regime characterized by the structure of individual pairs of dopants. Thereby, we trace back the phenomenology of the striped phase to its microscopic constituents and highlight the different challenges for observing the two regimes in quantum simulation experiments.

Fractional quantum Hall states host emergent anyons with exotic exchange statistics, but obtaining direct access to their topological properties in real systems remains a challenge. Neural-network wavefunctions provide a flexible computational approach, as they can represent highly correlated states without requiring a tailored basis. Here we use the neural-network variational Monte Carlo method to study the fractional quantum Hall effect on the torus and find the three degenerate ground states at filling factor nu=1/3. From these, we extract the modular S matrix via entanglement interferometry, a technique previously only applied to lattice models. The resulting S matrix encodes the quantum dimensions, fusion rules, and exchange statistics of the emergent anyons, providing a direct numerical demonstration of the topological order. The calculated anyon properties match the well-known theoretical and experimental results. Our work establishes neural-network wavefunctions as a powerful new tool for investigating anyonic properties.

Fractional quantum anomalous Hall (FQAH) effect, a lattice analogue of fractional quantum Hall effect, offers a unique pathway toward fault-tolerant quantum computation and deep insights into the interplay of topology and strong correlations. The exploration has been successfully guided by the paradigm of ideal flat Chern bands, which mimic Landau levels in both band topology and local quantum geometry. Yet, given the near-infinite possibilities for Bloch bands in lattices, it remains a major open question whether FQAH states can emerge in scenarios fundamentally different from this paradigm. Here we turn to a class of gapless flat bands, featuring divergent quantum geometry at singular band touching, non-integer Berry flux threading the Brillouin zone (BZ), and ill-defined band topology. Our exact diagonalization and density matrix renormalization group calculations unambiguously demonstrate FQAH phase that is virtually independent of the interaction strength, persisting from the weak-interaction to the strong-interaction limit. We find the stability of the FQAH states does not uniquely correlate with the singularity strength or the BZ-averaged quantum geometric fluctuations. Instead, the many-body topological order can adapt to the singular and fluctuating quantum geometric landscape by spontaneously developing an inhomogeneous carrier distribution, while its quenching accompanies the drop in the occupation-weighted Berry flux. Our work reveals a profound interplay between quantum geometry and many-body correlation, and significantly expands the design space for exploring FQAH effect and flat-band correlation phenomena in general.

We investigate the influence of correlations and plasmonic excitation on valence band-satellite of $α$-Fe, along with magnetic and spectral properties as function of temperature. Coulomb interaction parameters are obtained by systematically employing various schemes in constrained random phase approximation (cRPA). This study identifies the presence of valence band satellite in Fe at $\sim$6 eV binding energy supported by (i) substantial incoherent spectral weight in the valence band spectra obtained from Density Functional Theory plus Dynamical Mean Field Theory (DFT+DMFT) and (ii) plasmonic excitations in the frequency range $\sim$6-8 eV suggested by $G_0W_0$ calculations. We note presence of significant contribution of temperature-dependent Pauli-spin susceptibility indicating competing degree of itinerancy. $e_g$ state shows a strong temperature driven non-Fermi-liquid behavior emerging near $T_c$. Our results reveal a high-temperature orbital-selective loss of coherence eventually leads to a orbital selective collapse of magnetization at $T_c$, suggesting a ferromagnetic phase characterized by strong correlation- and temperature- dependent spectral features.

Entanglement plays a central role in quantum technologies, yet its characterization and control in materials remain challenging. Recent developments in spectrum-based entanglement witnesses have enabled new strategies for quantifying many-body entanglement in macroscopic materials. Here, we develop a protocol for detecting spin--orbital entanglement using experiment-accessible resonant inelastic x-ray scattering (RIXS). Central to our approach is the construction of a Hermitian generator from experimentally measurable spectra, which allows us to compute the quantum Fisher information (QFI) available in spin--orbital systems. The resulting QFI provides upper bounds for $k$-producible states and thus serves as a robust witness of spin--orbital entanglement. To account for realistic experimental limitations, we further extend our framework to include relaxed QFI bounds applicable to measurements lacking full polarization resolution.

Magnetic exchange interactions govern the macroscopic magnetic behavior of solids and underpin both fundamental spin phenomena and emerging technologies. The accurate and efficient determination of these interactions is therefore critical for predictive modeling of magnetic materials. Here we present a systematic first-principles comparison of three widely used approaches-the Least-Squares Total Energy (LSTE), the Four-State Total Energy (FSTE), and the Green's function-based Liechtenstein \textit{et al.} (LKAG) methods-applied to thirteen antiferromagnetic compounds. We introduce an framework for identifying the minimal supercells required for an accurate exchange parameter extraction in the FSTE method, significantly reducing computational cost while preserving precision. Our results show that LSTE and FSTE yield nearly identical exchange parameters, whereas the LKAG method reproduces the dominant exchange interactions but exhibits quantitative deviations. A detailed analysis of computational efficiency versus accuracy reveals that the LSTE scheme offers the most favorable balance, establishing a general, reproducible, and scalable workflow for Heisenberg mapping, while the FSTE approach remains the most straightforward for extracting specific exchange interactions.

Topological materials exhibit unique electronic structures that underpin both fundamental quantum phenomena and next-generation technologies, yet their discovery remains constrained by the high computational cost of first-principles calculations and the slow, resource-intensive nature of experimental synthesis. Recent machine-learning approaches, such as the heuristic topogivity rule, offer data-driven alternatives by quantifying each element's intrinsic tendency toward topological behavior. Here, we develop a quantum-classical hybrid artificial neural network (QANN) that extends this rule into a quantum-inspired formulation. Within this framework, the QANN maps compositional descriptors to quantum probability amplitudes, naturally introducing pairwise inter-element correlations inaccessible to classical heuristics. The physical validity of these correlations is substantiated by constructing an equivalent complex-valued neural network (CVNN), confirming both the consistency and interpretability of the formulation. Retaining the simplicity of chemical reasoning while embedding quantum-native features, our quantum-inspired rule enables efficient and generalizable topological classification. High-throughput screening combined with first-principles (DFT) validation reveals five previously unreported topological compounds, demonstrating the enhanced predictive power and physical insight afforded by quantum-inspired heuristics.

Using determinant Quantum Monte Carlo, we investigate the interplay between doping, inter-layer tunneling and onsite Hund's coupling in stabilizing superconductivity (SC) in a two-orbital model for the bilayer Nickelate $\text{La}_3\text{Ni}_2\text{O}_7$. With realistic dispersion and for certain values of the interaction parameters, the auxiliary-field-decoupled fermion Hamiltonian has Kramers anti-unitary symmetries which guarantee the absence of a sign problem. The same anti-unitary symmetries can also be used to show there is a second instability towards $(π,π)$ exciton condensation in the strong interaction limit. We indicate the possible connection between this exciton order and the enigmatic density wave state observed in experiment, and clarify the decisive role played by the inter-layer tunneling in the competition between SC and exciton condensation. Finally, possible directions on how to enhance the SC transition temperature and stabilize the SC phase are also discussed.

The recent experimental observations of high temperature superconductivity in bilayer nickelate have attracted lots of attentions. Previous studies have assumed a mirror symmetry $\mathcal M$ between the two layers and focused on uniform and clean superconducting states. Here, we show that breaking this mirror symmetry via an applied displacement field can stabilize a pair-density-wave (PDW) superconductor, which is similar to the Fulde--Ferrell--Larkin--Ovchinnikov (FFLO) state, but at zero magnetic field. Based on a mean-field analysis of a model of $d_{x^2-y^2}$ orbital with an effective inter-layer attraction, we demonstrate that the PDW phase is robust over a wide range of displacement field, interlayer hopping strengths, and electron fillings. Finally, we analyze disorder effects on interlayer superconductivity within the first Born approximation. Based on symmetry considerations, we show that pairing is weaken by disorders which break the mirror symmetry, even with unbroken time reversal symmetry. Our results establish bilayer nickelate as a tunable platform for realizing finite-momentum pairing and for exploring generalized disorder effects.

We propose an analytical method for the construction of Hartree-Fock phase diagrams for the (fermion) Hubbard model and various generalizations thereof. Such phase diagrams are traditionally constructed numerically, but we argue that, by using asymptotic techniques, it is possible to obtain analytic formulas approximating the curves separating the different phases to very high accuracy. To illustrate the new method, we apply it to the two-dimensional Hubbard model on the square lattice at zero temperature. This yields formulas for the Hartree-Fock phase boundaries that agree with, but also improve on, earlier numerical results. In particular, our results provide the first rigorous proof of the existence of mixed phases in this model.

Quantum spin liquids (QSLs) host exotic fractionalized magnetic and gauge-field excitations whose microscopic origins and experimental verification remain frustratingly elusive. In the absence of static magnetic order, the spin excitation spectrum constitutes the crucial probe of QSL behavior, but its theoretical computation remains a serious challenge. Here we employ state-of-the-art tensor-network methods to obtain the full dynamical spectral function of the $J_1$-$J_2$ kagome Heisenberg model and benchmark our results by tracking their evolution across the magnetically ordered and QSL phases. Reducing $|J_2|/J_1$ causes increasingly strong spin-wave renormalization, flattening these modes then merging them into a continuum characteristic of deconfined spinons at all finite energies in the QSL. The low-energy continuum and the occurrence of gap closure at multiple high-symmetry points identify this gapless QSL as the U(1) Dirac spin liquid. These results establish a unified understanding of spin excitations in highly frustrated quantum magnets and provide clear spectral fingerprints for experimental detection in candidate kagome QSL materials.

In a recent experiment [Nature 642, 58 (2025)], a new record for the superconducting critical temperature $T_c$ among infinite-layer nickelates has been reported in doped SmNiO$_2$. Here, we use the cutting-edge dynamical vertex approximation (D$Γ$A), and qualitatively as well as quantitatively reproduce the $T_c$ vs. doping dome for this compound. Encouraged by this, we go further and identify a path towards realizing even higher $T_c$'s by changing the cation along the line Nd$\rightarrow$Sm$\rightarrow$Y$\rightarrow$Lu with matching substrates. The successively smaller cation radius allows for smaller lattice constants of the substrate. This in turn increases the in-plane hopping and thus eventually $T_c$.

Ultracold multicomponent fermions (atoms/molecules) loaded in optical lattices provide an ideal platform for simulating SU($N$) Hubbard models that host unconventional many-body quantum states beyond SU(2). A prime example is the attractive three-color Hubbard model, in which trion states emerge at strong coupling. Nevertheless, much of its trion ordering on two-dimensional lattices remains uncertain. Here, we employ the determinant quantum Monte Carlo (DQMC) method to simulate the attractive three-color Hubbard model on a $π$-flux square lattice at half filling. We show that color-dependent attractive interaction can induce coexisting charge density wave (CDW) and Néel ordered states in the three-color $π$-flux Hubbard model. In particular, enhanced charge fluctuations (cf. honeycomb lattice) cause much stronger Néel ordering on the $π$-flux square lattice. The coexisting charge and Néel orders survive up to a melting temperature, at which they vanish simultaneously. The Ginzburg-Landau (GL) analysis on the coexistence of CDW and Néel orders demonstrates how color-dependent Hubbard interactions stabilize coexisting orders from the perspective of GL free energy principle.

By means of cellular dynamical mean-field theory (CDMFT) we study how short-range correlations drive the breakdown of the self-consistent perturbation theory in two-dimensional systems and the most relevant physical consequences associated to it. To this aim, we first derive in a structured and consistent way the Bethe-Salpeter equation (BSE) formalism at the CDMFT level in all physical channels, explicitly addressing the important aspect of the related Ward identities. In this context, we perform systematic calculations of the BSE for the two-dimensional Hubbard model at half-filling at intermediate coupling. Our study illustrates how the divergence of a fundamental building block of the BSE in the charge channel, the two-particle irreducible vertex, systematically occurs at lower interactions than in the (purely local) DMFT case, due to short-range antiferromagnetic fluctuations. Further, the change of sign of the eigenvalues of the generalized charge susceptibility associated to the vertex divergences is identified as the essential prerequisite to drive, at larger interaction values, the physics of the Mott transition in two dimensions, as well as of the adjacent phase-separation instabilities.

We analyze the impact of the lattice geometry on the thermodynamic transition to magnetically ordered phases in strongly interacting electron systems for various Bravais lattices in three and four dimensions, including both local and nonlocal correlation effects. In a first step we use the dynamical mean field theory (DMFT), which takes into account purely local correlations, to calculate the magnetic susceptibilities of the Hubbard model on three (3d-sc) and four dimensional (4d-sc) simple cubic/hypercubic, as well as on three dimensional body- (bcc) and face-centered (fcc) cubic lattices, and determine the transition temperature to the corresponding magnetically-ordered state. In a second step, we exploit the dynamical vertex approximation (D$Γ$A), a diagrammatic extension of DMFT, to include the effect of nonlocal correlations which are particularly important in the vicinity of the corresponding phase transition. For the bipartite 3d-sc, 4d-sc and bcc lattices nonlocal fluctuations lead to a substantial reduction of the DMFT transition temperature consistent to the overall tendency of mean-field approaches to overestimate the stability of ordered phases. As expected, the magnitude of the difference between the DMFT, being exact in the limit of large connectivity/dimensions, and D$Γ$A transition temperatures decreases with increasing coordination number. On a more practical perspective, these results also provide a reasonable guidance to evaluate the expected overestimation of the DMFT ordering temperature for different material geometries. For the fcc lattice, on the other hand, the ordered phase observed in DMFT vanishes completely within D$Γ$A which is consistent with the existence of strong geometric frustration in this lattice.

At strong on-site repulsion $ U $, the fermionic Hubbard model realizes an extremely correlated electron system. In this regime, it is natural to derive the low-energy physics with the help of non-canonical operators acting on a projected Hilbert space without double occupancies. Using a strong-coupling functional renormalization group technique, we study the physics of such extreme correlations in the strict $ U = \infty $ limit, where only kinematic interactions due to the Hilbert space projection remain. For nearest-neighbor hopping on a square lattice, we find that the electronic spectrum is significantly renormalized, with bandwidth and quasi-particle residue strongly decreasing with increasing electron density. On the other hand, damping and particle-hole asymmetry increase, while a polaronic continuum forms in the hole sector, below the single-particle band. Fermi liquid phenomenology applies only at low densities, where the system remains paramagnetic. At higher densities, we find a bad metal with strong magnetic correlations, indicating that the ground state is the Nagaoka ferromagnet at high densities and a stripe antiferromagnet at intermediate densities. Both in the paramagnetic and the ferromagnetic regimes, we observe a violation of Luttinger's theorem.

High-temperature and high-magnetic-field-induced re-entrant superconductivity has been discovered in the infinite-layer nickelate $\mathrm{Sm_{1-x-y} Eu_x Ca_y Ni O_2}$ (SECNO). Infinite-layer nickelates are the closest known analogues of high-$\mathrm{T}_c$ cuprate superconductors, yet they host distinct magnetic ground states. Using low-energy muon spin relaxation and polarized neutron reflectometry, we reveal the magnetic order in SECNO. We find that magnetic freezing occurs at a higher-temperature than in other nickelate compounds, and that a substantial net magnetization of 55 $\,\mathrm{kA}\,\mathrm{m}^{-1}$ $\pm10 \,\mathrm{kA}\,\mathrm{m}^{-1}$ emerges and remains largely unchanged across the superconducting transition. The magnetism in SECNO is disordered and nonuniform.

Preparing ground states of strongly correlated quantum systems is a central goal in quantum simulation and optimization. The feedback-based quantum algorithm (FALQON) provides an attractive alternative to variational methods with a fully quantum feedback rule, but it fails in the presence of spectral degeneracies, where the feedback signal collapses and the evolution cannot reach the ground state. Using the Fermi-Hubbard model on lattices up to 3x3, we show that this breakdown appears at half-filling on the 2x2 lattice and extends to both half-filled and doped configurations on the 3x3 lattice. We then introduce an imaginary-time-enhanced FALQON (ITE-FALQON) scheme, which inserts short imaginary-time evolution steps into the feedback loop. The hybrid method suppresses excited-state components, escapes degenerate subspaces, and restores monotonic energy descent. The ITE-FALQON achieves a reliable ground-state convergence across all fillings, providing a practical route to scalable ground-state preparation in strongly correlated quantum systems.

We investigate the binding of holes and the emergence of competing spin-charge order in the simple and extended Hubbard model using exact diagonalization on the 3x4 cylindrical lattice. For the simple Hubbard model (V=0), we find weakly bound hole pairing mediated by magnetic correlations at intermediate repulsive U, without any evidence of phase separation. Introducing nearest-neighbor interaction V reveals a rich phase diagram: attractive V drives multi-hole clustering and phase separation with localized magnetic quenching, while repulsive V stabilizes charge-density-wave (CDW) order that coexists with bound hole pairs within a modulated magnetic background. At strong coupling (U=10), the competition sharpens, with attractive V overcoming on-site repulsion to form magnetically quenched clusters and repulsive V producing robust CDW order that constrains pairing. Real-space analysis of spin and charge correlations provides microscopic evidence of distinct binding mechanisms -- phase separation versus correlation-mediated pairing -- depending on the sign and strength of intersite interaction V . Our results establish a comprehensive picture of how nonlocal Coulomb interactions reshape the landscape of hole-binding and collective order in correlated electron systems.

Electron-phonon coupling is believed to be responsible for many spectral anomalies in the cuprate superconductors. In particular, the $B_{1g}$ buckling mode of the oxygen ion in the $CuO_{2}$ plane has been proposed to be responsible for the dramatic peak-dip-hump(PDH) structure in the anti-nodal spectrum. The recent observation of the exceptional flat quasiparticle dispersion in the anti-nodal region and the sudden suppression of the PDH structure around the pseudogap end point cast doubts on such a scenario. Instead, a scenario involving the coupling to the antiferromagnetic spin fluctuation seems to resolve both puzzles naturally. Here we present a systematic study on the interplay between antiferromagnetic spin fluctuation and electron-phonon coupling in the cuprate superconductors. We show that the coupling strength to the $B_{1g}$ buckling mode is strongly suppressed by the vertex correction caused by the antiferromagnetic spin fluctuation in the $\mathbf{q}\rightarrow 0$ limit as a result of the destructive interference between electron-phonon coupling at electron momentum differ by the antiferromagnetic wave vector. Counterintuitively, we find that the same vertex correction enhances the phonon contribution to the PDH structure. We also find that while the coupling to either the antiferromagnetic spin fluctuation or the $B_{1g}$ buckling mode can generate a PDH structure in the anti-nodal spectrum with similar phenomenologies, the sudden suppression of such a structure around the pseudogap end point should be mainly attributed to the dramatic change in the nature of the spin fluctuation at such a critical doping. We suggest to take the PDH structure in the anti-nodal spectrum as a spectral signature for the emergence of fluctuating local moment in the pseudogap phase and the entrance of a doped Mott insulating state.

We propose a single-layer tensor network framework for the variational determination of ground states in two-dimensional quantum lattice models. By combining the nested tensor network method [Phys. Rev. B 96, 045128 (2017)] with the automatic differentiation technique, our approach can reduce the computational cost by three orders of magnitude in bond dimension, and therefore enables highly efficient variational ground-state calculations. We demonstrate the capability of this framework through two quantum spin models: the antiferromagnetic Heisenberg model on a square lattice and the frustrated Shastry-Sutherland model. Even without GPU acceleration or symmetry implimention, we have achieved the bond dimension of 9 and obtained accurate ground-state energy and consistent order parameters compared to prior studies. In particular, we confirm the existence of an intermediate empty-plaquette valence bond solid ground state in the Shastry-Sutherland model. We have further discussed the convergence of the algorithm and its potential improvements. Our work provides a promising route for large-scale tensor network calculations of two-dimensional quantum systems.

The ionic Hubbard model is a paradigmatic setup for studying the competition between band and Mott insulating behavior. Within the variationally exact in infinite dimensions Gutzwiller approximation, we derive a compact analytic expression for the phase boundary between Mott and band insulator. While the method reproduces the expected band-Mott insulator phenomenology, it does not capture the correlated metallic state at finite staggered potential found for example in dynamical mean-field theory. This absence highlights that the metallic phase originates from incoherent Hubbard-band physics rather than Fermi-liquid behavior well captured by Gutzwiller approximation. Our formulation establishes a concise variational framework to ionic Hubbard model, with natural extensions to nonequilibrium setups and spin-exchange dynamics.

We introduce an unsupervised machine-learning framework that discovers optimally compressed representations of quantum many-body ground states. Using an autoencoder neural network architecture on data from $L$-site Fermi-Hubbard models, we identify minimal latent spaces with a sharp reconstruction quality threshold at $L-1$ latent dimensions, matching the system's intrinsic degrees of freedom. We demonstrate the use of the trained decoder as a differentiable variational ansatz to minimize energy directly within the latent space. Crucially, this approach circumvents the $N$-representability problem, as the learned manifold implicitly restricts the optimization to physically valid quantum states.

We study superconductivity and superfluid weight of the two-dimensional $α$-$\mathcal{T}_3$ lattice with on-site asymmetries, hosting an isolated quasi-flat band with tunable bandwidth via a parameter $α$. Within a mean-field approximation of the attractive Hubbard model, we obtain the superconducting order parameters on the three inequivalent sublattices and show their strong dependence on $α$, interaction strength, and electron filling. At quasi-flat band filling, a superconducting gap opens and grows power-law fast with interaction strength, instead of the usual slow exponential growth, due to diverging density of states. We calculate the superfluid weight from linear response theory and study its band dispersion and geometric contributions. While the conventional part proportional to band derivatives is suppressed in the quasi-flat band regime, the contribution dominated by the quantum metric grows linearly for small interaction strength. We further demonstrate how tuning $α$ enhances the quantum metric and thus the geometric superfluid weight especially near half-filling, while increasing on-site asymmetries increases the conventional contribution by broadening the quasi-flat band. We obtain the Berezinskii-Kosterlitz-Thouless transition temperature and demonstrate its strong dependence and enhancement with the parameter $α$. Our results establish a tunable flat band system, the $α$-$\mathcal{T}_3$ lattice model, as a candidate for tunable quantum geometry and superfluid weight and as a prototype of related behavior in tunable quantum materials.

High-temperature superconductivity in cuprate materials remains a major challenge in physics due to the complexity of their strongly correlated electronic states. Interfacial strain is a powerful lever for tuning electronic correlations in complex oxides, offering new pathways to control emergent quantum phases. Here, we report the discovery of interfacial strain modulated correlated plasmons observed exclusively in superconducting La1.85Sr0.15CuO4 (LSCO) through spectroscopic ellipsometry. This form of plasmons is absent in the non-superconducting LSCO counterparts. Detailed analysis reveals that these correlated plasmons, arising from the collective excitations within Mott-correlated bands, are driven by long-range electronic correlations in the Cu-O planes. Furthermore, long-range electronic correlations, intricately modulated by interfacial strain, may play a crucial role in the emergence of superconductivity and in tuning the transition temperature. Dynamical cluster approximation (DCA) with quantum Monte Carlo (QMC) calculations of the extended Hubbard model suggest that long-range Coulomb interactions play an important role in LSCO, showing good agreement with our experimental findings. The collective evidence from both the experimental results and theoretical findings provides new insights into the nature of collective excitations and their pivotal role in the emergence of high-temperature superconductivity.

We investigate the impact of quantum and thermal phase fluctuations on the suppression of superconducting order in two-dimensional systems. Within the two-dimensional quantum XY model in the phase representation, where on-site interaction terms govern quantum phase fluctuations, we perform extensive path-integral quantum Monte Carlo simulations. The resulting temperature-interaction phase diagram establishes the presence of a well-defined critical line ending at a quantum critical point at vanishing temperature with no indication of reentrant behavior. We further demonstrate that the resistance above the critical line reproduces the two expected different critical behaviors. For stronger interactions, above the quantum critical point, the system exhibits a crossover to an insulating regime at low temperatures. Finally, Monte Carlo calculations of current-current correlation functions enable us to extract the frequency-dependent conductivity in both superconducting and normal regimes, revealing a finite-frequency response that we attribute to quantum phase fluctuations.

We investigate an S=3/2 quantum spin model on a two-dimensional honeycomb lattice that continuously interpolates between two paradigmatic quantum disordered states with distinct entanglement structures: the Kitaev quantum spin liquid and the Affleck-Kennedy-Lieb-Tasaki (AKLT) valence bond solid. Combining classical, semi-classical, and exact diagonalization approaches, we map out the ground-state phase diagram and elucidate the role of quantum fluctuations across the entire parameter range. While classical and semi-classical frameworks predict noncoplanar orders competing with a collinear Néel state, we find these phases to be fragile: once full quantum fluctuations are included, they melt into a quantum-entangled state characterized by suppressed spin correlations and enhanced entanglement entropy. Our findings highlight how competition between qualitatively different quantum disordered phases provides a fertile playground for unconventional phases emerging from their interplay and quantum fluctuations.

The breaking of spatial inversion symmetry in condensed matter gives rise to intriguing physical properties, such as ferroelectricity, piezoelectricity, spin-momentum locking, and nonreciprocal responses. Here we propose that a concentration gradient, which often persists as a quasi-stable nonequilibrium state with long relaxation times in solids, can serve as a general platform for inversion symmetry breaking. We demonstrate this concept in an epitaxial thin film of the hydrogen-doped SmFeAsO (Sm1111:H) superconductor with a depthwise hydrogen-concentration gradient introduced via an optimized topotactic reaction. This film exhibits nonreciprocal charge transport, meaning that the electrical resistance depends on the direction of the applied current, which serves as a key signature of broken inversion symmetry. A pronounced nonreciprocal signal emerges in the vicinity of the superconducting transition, which we attribute to vortex-motion nonreciprocity arising from an asymmetric pinning landscape created by the hydrogen-concentration gradient. Owing to the high critical temperature of Sm1111:H, vortex-origin nonreciprocity is observed above 40 K, representing the highest temperature reported to date among single bulk materials without an artificially hetero-layered structure. Our findings establish concentration-gradient engineering as a versatile and broadly applicable route for realizing inversion-broken states in otherwise centrosymmetric hosts, opening pathways toward a broader landscape of odd-parity-driven functionalities.

Seniority is a useful way of organizing Hilbert space for strongly correlated systems. The exact zero-seniority wave function, doubly-occupied configuration interaction (DOCI), provides accurate results (given the right orbitals) for many strongly-correlated electronic systems, but has combinatorial computational cost. In many cases, pair coupled cluster doubles provides a polynomial-cost approximation that closely reproduces the energies of DOCI, but it breaks down in some cases and, as shown herein, it does not provide particularly good density matrices. In this work, we demonstrate that by using the Jordan-Wigner approximation to turn the seniority zero problem back into a fermionic one, we can provide variational results of DOCI quality for the Hubbard model and a few small molecular dissociation examples, with polynomial cost, both for the energies and for density matrices, all while being protected from collapse.

Ruthenates, materials with a single magnetic Ruthenium (Ru) atom, often display an exotic array of ground states ranging from superconductivity to altermagnetism. In this work, we investigated the magnetic ground state of a least explored member of the 4d3 double perovskite ruthenate series A2ScRuO6 (A = Ca, Sr, Ba): Ca2ScRuO6. Interestingly, temperature-dependent bulk susceptibility curve shows ferrimagnetic-like behaviour above the magnetic ordering at around 40 K, which were corroborated by the identification of the mixed valence states, Ru5+ and Ru4+ via X-ray absorption spectroscopy. Structural analysis further revealed atomic-site exchange between the Ru and Sc sites, which results in the Ru mixed valence states. Neutron powder diffraction measurements detected the presence of magnetic Bragg peaks at a low temperature near 4 K and a moderate magnetoelastic coupling near the ordering temperature of 40 K. However, the corresponding symmetry analysis shows a weak Type I antiferromagnetic ground state with a reduced magnetic moment of 1.1μB/Ru atom. Our findings establish an unusual magnetic ground state in the Mott insulating Ca2ScRuO6, where a long range ordered antiferromagnet coexists with small magnetic clusters, which manifests a ferrimagnetic-like high temperature inverse magnetic susceptibility. This system presents a unique platform to study long-range magnetic order in the presence of antisite disorder.

We propose an alternative to the infinite density-matrix renormalization approach for accessing quantum many-body states within a finite-size calculation that faithfully mimics the thermodynamic limit. Our method constructs environment matrix product operators (MPOs) representing the Hamiltonian of semi-infinite regions surrounding the target system. Starting from the finite-size ground-state MPS, we contract its Hamiltonian representation to generate effective environment MPOs, which are then attached to a renewed finite system in a recursive manner. This iterative embedding drives the system toward a bulk-like state with negligible finite-size effects. The scheme requires no assumption of homogeneity and achieves unprecedentedly long real-time dynamics free from boundary reflections.

Recent experiments have revealed signatures of small Fermi pockets in the pseudogap phase of cuprate superconductors, most notably the Yamaji effect observed in $\mathrm{HgBa}_2\mathrm{CuO}_{4+δ}$. The Yang-Rice-Zhang (YRZ) model provides a successful phenomenological description of the pseudogap state and naturally predicts such small pockets. In this work, we use a microscopic framework to calculate angle-dependent magnetoresistance and quantum oscillation within the YRZ model. Our calculations simultaneously reproduce the experimentally observed Yamaji oscillations and the Shubnikov-de Haas oscillation corresponding to a pocket area of about $p/8$, with $p$ the hole density. By further testing the effect of Green's-function zeros, we confirm that isolated zeros leave the oscillation period unchanged, whereas an extended zero segment suppresses and modifies the oscillation. Our findings demonstrate that the YRZ model captures essential features of the pseudogap regime and provides a general quantum approach that can be applied to more complex electronic structures.

We evaluate the Sample-based Krylov Quantum Diagonalization (SKQD) algorithm on one- and two-dimensional Heisenberg models, including strongly correlated regimes in which the ground state is dense. Using problem-informed initial states and magnetization-sector sweeps, SKQD accurately reproduces ground-state energies and field-dependent magnetization across a range of anisotropies. Benchmarks against DMRG and exact diagonalization show consistent qualitative agreement, with accuracy improving systematically in more anisotropic regimes. We further demonstrate SKQD on quantum hardware by implementing 18- and 30-qubit Heisenberg chains, obtaining magnetization curves that match theoretical expectations. Simulations on small 2D square-lattice systems further demonstrate that the method applies effectively beyond 1D geometries.

Kitaev materials are of great interest due to their potential in realizing quantum spin liquid (QSL) states and applications in topological quantum computing. In the pursuit of realizing Kitaev QSL, a Mott insulator with strong bond-dependent frustration and weak geometric frustration is highly desirable. Here we explore Kitaev physics in the van der Waals triangular antiferromagnet (AF) CoI$_2$, through the spin-orbital states and Wannier function analyses, exact diagonalization and density matrix renormalization group study of the electronic structure and magnetic properties. We find that the high-spin Co$^{2+}$ ion is in the $J_\mathrm{eff}=1/2$ state because of strong spin-orbit coupling, and the weak trigonal elongation and crystal field contribute to the observed weak in-plane magnetic anisotropy. The strong $t_{2g}$-$e_g$ hopping via the strong Co 3$d$-I 5$p$ hybridization gives rise to a strong Kitaev interaction ($K_1$) at the first nearest neighbors (1NN), and the long Co-Co distance and the weak $t_{2g}$-$t_{2g}$ hoppings determine a weak Heisenberg interaction $J_1$. The resultant $|K_1/J_1|$ = 6.63 confirms a strong bond-dependent frustration, while the geometric frustration due to the 3NN Heisenberg interaction $J_3$ gets involved, and they all together result in the experimental helical AF order in CoI$_2$. We then propose to suppress the $J_3$ using a partial Mg substitution for Co, and indeed we find that Co$_{2/3}$Mg$_{1/3}$I$_2$ has the much reduced geometric frustration but hosts the robust bond-dependent frustration, and thus it would be a promising Kitaev material being so far closest to the QSL state.

Extending corresponding results for matrix product states [Verstraete and Cirac, PRB 73, 094423 (2006); Schuch et al. PRL 100, 030504 (2008)], it is shown how the approximation error of tree tensor network states (TTNS) can be bounded using Schmidt spectra or Rényi entanglement entropies of the target quantum state. Conversely, one obtains bounds on TTNS bond dimensions needed to achieve a specific approximation accuracy. For tree lattices, the result implies that efficient TTNS approximations exist if $α<1$ Rényi entanglement entropies for single-branch cuts obey an area law, as in ground and low-energy states of certain gapped systems.

Reliably simulating two-dimensional many-body quantum dynamics with projected entangled pair states (PEPS) has long been a difficult challenge. In this work, we overcome this barrier for low-energy quantum dynamics by developing a stable and efficient time-dependent variational Monte Carlo (tVMC) framework for PEPS. By analytically removing all gauge redundancies of the PEPS manifold and exploiting tensor locality, we obtain a numerically well-conditioned stochastic reconfiguration (SR) equation amenable to robust solution using the efficient Cholesky decomposition, enabling long-time evolution in previously inaccessible regimes. We demonstrate the power and generality of the method through four representative real-time problems in two dimensions: (I) chiral edge propagation in a free-fermion Chern insulator; (II) fractionalized charge transport in a fractional Chern insulator; (III) vison confinement dynamics in the Higgs phase of a Z2 lattice gauge theory; and (IV) superfluidity and critical velocity in interacting bosons. All simulations are performed on 12x12 or 13x13 lattices with evolution times T = 10 to 12 using modest computational resources (1 to 5 days on a single GPU card). Where exact benchmarks exist (case I), PEPS-tVMC matches free-fermion dynamics with high accuracy up to T = 12. These results establish PEPS-tVMC as a practical and versatile tool for real-time quantum dynamics in two dimensions. The method extends the reach of classical tensor-network simulations for studying elementary excitations in quantum many-body systems and provides a valuable computational counterpart to emerging quantum simulators.

Simulations of quantum matter rely mainly on Kohn-Sham density functional theory (DFT), which often fails for strongly correlated systems. Quantum embedding (QE) theories address this limitation by mapping the system onto an auxiliary embedding Hamiltonian (EH) describing fragment-environment interactions, but the EH is typically large and its iterative solution is the primary computational bottleneck. We introduce a linear foundation model for QE that utilizes principal component analysis (PCA) to compress the space of quantum states needed to solve the EH within a small variational subspace. Using a data-driven active-learning scheme, we learn this subspace from EH ground states and reduce each embedding solve to a deterministic ground-state eigenvalue problem in the reduced space. Within the ghost Gutzwiller approximation (ghost-GA), we show for a three-orbital Hubbard model that a variational space learned on a Bethe lattice is transferable to square and cubic lattices without additional training, while substantially reducing the cost of the EH step. We further validate the approach on plutonium, where a single variational space reproduces the energetics of all six crystalline phases while reducing the cost of the EH solution by orders of magnitude. This provides a practical route to overcome the main computational bottleneck of QE frameworks, paving the way for high-throughput ab initio simulations of strongly correlated materials at a near-DFT cost.

Tensor networks provide a powerful new framework for classifying and simulating correlated and topological phases of quantum matter. Their central premise is that strongly correlated matter can only be understood by studying the underlying entanglement structure and its associated (generalised) symmetries. In essence, tensor networks provide a compressed, holographic description of the complicated vacuum fluctuations in strongly correlated systems, and as such they break down the infamous many-body exponential wall. These lecture notes provide a concise overview of the most important conceptual, computational and mathematical aspects of this theory.

Cr$_{1+δ}$Te$_2$, a two-dimensional van der Waals ferromagnet, displays a contested magnetic structure, poised between collinear and non-collinear spin configurations. In this work, we investigate the magnetic structure of Cr$_{1.33}$Te$_2$ at the microscopic level by combining single-crystal neutron diffraction, X-ray absorption spectroscopy, and first-principles calculations. Neutron diffraction measurements reveal a distinct collinear spin alignment, whereas spectroscopic analyses reveal inherent structural vacancies at both Cr and Te sites. These vacancies lead to local symmetry breaking that elevates the orbital degeneracy of the Cr 3$d$ states, as demonstrated by our first-principles analysis. The resulting modification of magnetocrystalline anisotropy emerges as the key mechanism stabilising the collinear magnetic ground state over the non-collinear one in the presence of vacancies. Our findings uncover a vacancy-driven route to control spin anisotropy and magnetic ordering in layered ferromagnets, offering new insights into the design of tunable 2D magnetic materials.

The string hypothesis for Bethe roots represents a cornerstone in the study of quantum integrable systems, providing access to physical quantities such as the ground-state energy and the finite-temperature free energy. While the $t-W$ scheme and the inhomogeneous $T-Q$ relation have enabled significant methodological advances for systems with broken $U(1)$ symmetry, the underlying physics induced by symmetry breaking remains largely unexplored, due to the previously unknown distributions of the transfer-matrix roots. In this paper, we propose a new approach to determining the patterns of zero roots and Bethe roots for the $Λ-θ$ and inhomogeneous Bethe ansatz equations using tensor-network algorithms. As an explicit example, we consider the isotropic Heisenberg spin chain with non-diagonal boundary conditions. The exact structures of both zero roots and Bethe roots are obtained in the ground state for large system sizes, up to ($N\simeq 60$ and $100$). We find that even in the absence of $U(1)$ symmetry, the Bethe and zero roots still exhibit a highly structured pattern. The zero roots organize into bulk strings, boundary strings, and additional roots, forming two dominant lines with boundary-string attachments. Correspondingly, the Bethe roots can be classified into four distinct types: regular roots, line roots, arc roots, and paired-line roots. These structures are associated with a real-axis line, a vertical line, characteristic arcs in the complex plane, and boundary-induced conjugate pairs. Comparative analysis reveals that the $t-W$ scheme generates significantly simpler root topologies than those obtained via off-diagonal Bethe Ansatz.

The kagome lattice, with its unique geometric structure, has emerged as a leading platform for exploring quantum many-body physics, particularly in the study of quantum spin liquids (QSLs) and unconventional superconductivity. This review highlights recent advancements in the investigations of QSLs, fractional magnetization plateau phases in kagome antiferromagnets, and unconventional superconductivity in vanadium-based kagome superconductors. We begin by examining the classical ground-state properties of the nearest-neighbor kagome antiferromagnetic Heisenberg model and introducing recent experimental progress in the study of QSLs and fractional magnetization plateau phases. Next, we discuss the fermionic description of the QSL states, along with related gauge theory and the variational Monte Carlo (VMC) method. We then focus on discussing the VMC studies of QSLs and magnetization plateau phases in kagome antiferromagnets. For superconductivity in kagome systems, we first analyze the characteristics of the electronic structure and the possible associated electronic instabilities. Finally, we review recent experimental advances in unconventional superconductivity in AV$_3$Sb$_5$ (A = K, Rb, Cs), with a particular focus on chiral superconductivity and pairing density waves.

We report a spin-(1/2, 5/2) three-leg ladder realized in a radical-Mn polymer, exhibiting an antiferromagnetic transition and magnetization curves accurately described by classical mean-field theory. Although the underlying spin model intrinsically supports strong quantum fluctuations, as confirmed by quantum Monte Carlo simulations, the real system shows an anomalously complete suppression of quantum behavior. These findings provide a key experimental benchmark for the quantum-classical crossover and suggest that lattice topology can play a crucial role in tuning the balance between quantum and classical physics in strongly correlated systems.

Two-dimensional (2D) materials are promising candidates for next-generation spintronic devices due to their tunable properties and potential for efficient spin-charge interconversion. However, discovering materials with intrinsically high spin Hall conductivity (SHC) is hindered by the vast chemical space and expensive nature of conventional experimental and first-principles methods. In this work, we employ an active learning framework to accelerate the discovery of high-SHC 2D materials. Machine learning (ML) models were trained on SHC values computed from density functional theory calculations, incorporating the Kubo formalism via tight-binding Hamiltonians constructed from maximally localized Wannier functions, with explicit treatment of spin-orbit coupling. Starting from random but chemically diverse 24 2D systems, the dataset was expanded to 41 cases (from an overall pool of around 2000 materials) over three active learning loops using an expected improvement acquisition strategy. The ML technique successfully identified several high SHC candidates with the best candidate exhibiting a SHC of 271.52 (hbar/e) Ohm^-1, nearly 23 times higher than the top performer in the initial round. Beyond candidate discovery, several features such as orbital symmetry near the Fermi energy, types of atomic species, material composition, covalent radii, and electronegativity of constituent atoms were found to play critical role in shaping the spin Hall response in 2D systems. The data generated is made publicly available to facilitate further advances in 2D spintronics.

Relaxation processes in topological phases such as quantum spin liquids are controlled by the dynamics and interaction of fractionalized excitations. In layered materials hosting two-dimensional topological phases, elementary quasiparticles can diffuse freely within the layer, whereas only pairs (or more) can hop between layers - a fundamental consequence of topological order. Using exact solutions of emergent nonlinear diffusion equations and particle-based stochastic simulations, we explore how pump-probe experiments can provide unique signatures of the presence of $2d$ topological excitations in a $3d$ material. Here we show that the characteristic time scale of such experiments is inversely proportional to the initial excitation density, set by the pump intensity. A uniform excitation density created on the surface of a sample spreads subdiffusively into the bulk with a mean depth $\bar z$ scaling as $\sim t^{1/3}$ when annihilation processes are absent. The propagation becomes logarithmic, $\bar z \sim \log t$, when pair-annihilation is allowed. Furthermore, pair-diffusion between layers leads to a new decay law for the total density, $n(t) \sim (\log^2 t)/t$ - slower than in a purely $2d$ system. We discuss possible experimental implications for pump-probe experiments in samples of finite width.

In the local-moment regime, heavy fermions are most economically described by a compact U(1) gauge theory. With this formulation of the Kondo lattice, we study a spin chain coupled to two-dimensional Dirac conduction electrons. The spin chain is described by fermionic partons carrying spin and U(1) gauge charge. The heavy-fermion quasiparticle is a bound state of a U(1) matter field carrying unit electric and U(1) gauge charge, and the fermionic parton. Using sign-problem-free determinant quantum Monte Carlo simulations, we identify two symmetry-equivalent regimes: a heavy-fermion metal with a sharp composite-fermion resonance and robust low-frequency transport, and a Kondo-breakdown metal with an incoherent resonance and vanishing low-frequency transport. For any finite lattice extent in the direction perpendicular to the chain, the Luttinger volume of the heavy-fermion phase counts both composite and conduction electrons, while in the Kondo-breakdown phase it counts only the conduction electrons. The evolution of the composite-fermion spectrum, dynamical spin structure factor, and optical conductivity provides a nonperturbative demonstration of gauge-mediated Kondo breakdown and establishes transport fingerprints of an orbital-selective Mott transition in the context of U(1) gauge theories of heavy fermions.

Decisive experimental confirmation of the $U(1)$ quantum spin liquid phase in quantum spin ice remains an outstanding challenge. In this work, we propose stray-field magnetometry as a direct probe of the emergent photons -- the gapless excitation of the emergent electrodynamics in quantum spin ice. The emergent photons are transverse magnetization waves, which, in a finite sample, form discrete modes governed by one of two sets of natural boundary conditions: ``insulating'' or ``superconducting''. Considering cavity and thin film geometries, we find that the spectrum and spatial structure of the stray magnetic noise provide a sharp qualitative signature of the underlying electrodynamics. The predicted stray-field noise power lies comfortably within the detection range of present-day solid-state defect magnetometry.

We identify an unconventional form of the non-Hermitian skin effect that occurs not in position space but in many-body Fock space, which we call the Fock space skin effect (FSSE). Using quantum dimer models, we characterize FSSE analytically and numerically, and propose a concrete route toward its realization in Rydberg atom arrays. The dimer constraint is enforced through Rydberg gadgets employing the blockade mechanism, while directional reservoirs generate non-Hermitian flipping amplitudes. We show that FSSE enables the preparation of gapped spin liquid states, and in particular, we demonstrate how a Rydberg geometry realizing a square lattice quantum dimer model with next-nearest neighbor dimers can be driven by non-Hermiticity into an exact spin liquid ground state. Our results establish Fock-space non-Hermiticity as a powerful principle for engineering exotic quantum phases and dynamical state-preparation protocols.

The homogeneous electron gas is a cornerstone of quantum condensed matter physics, providing the foundation for developing density functional theory and understanding electronic phases in semiconductors. However, theoretical understanding of strongly-correlated electrons in realistic semiconductor systems remains limited. In this work, we develop a neural network based variational approach to study quantum wells in three dimensional geometry for a variety of electron densities and well thicknesses. Starting from first principles, our unbiased AI-powered method reveals metallic and crystalline phases with both monolayer and bilayer charge distributions. In the emergent bilayer, we discover a new quantum phase of matter: the electronic quasicrystal.

Spin space groups (SSGs) impose fundamentally different constraints on magnetic configurations in real and reciprocal spaces. As a consequence, the correspondence between real-space and momentum-space spin arrangements is far richer than traditionally assumed. Building on the complete enumeration of SSGs, we develop a systematic, symmetry-based framework that classifies all possible spin arrangements allowed by these groups. This unified approach naturally incorporates conventional magnetic orders, altermagnetism, and p-wave magnetism as distinct symmetry classes. Crucially, our classification predicts a variety of novel magnetic phases, highlighted by the discovery of the coplanar even-wave magnet: a state that is non-collinear in real space but hosts a collinear even-wave spin polarization in k-space. Analysis of a minimal model reveals that this phase is characterized by non-quantized spin polarization and exhibits a novel mechanism for symmetry-enforced zero polarization on non-degenerate bands. Extending the framework from bulk crystals to layer SSGs appropriate for two-dimensional systems, we further predict layered counterparts and provide symmetry guidelines for designing bilayer coplanar p-wave and even-wave magnets. We further validate this finding through first-principles calculations and propose CoCrO4 as a promising candidate for its experimental realization, thereby demonstrating the completeness and predictive power of the SSG-based classification of magnetic orders.

Phase separated manganites provide a unique platform to study the dynamics of competing electronic and structural orders in correlated systems. In $La_{0.275}Pr_{0.35}Ca_{0.375}MnO_{3}$ (LPCMO), we use temperature cycling Raman spectroscopy to uncover a previously unidentified regime of structural irreversibility, emerging from the interplay between lattice distortions and phase competition across the phase separation and charge and orbital ordering temperatures. This irreversible behavior encodes a thermal memory effect reflecting the system's history dependent energy landscape. Correlated magnetic and transport responses confirm the coupling between lattice and electronic degrees of freedom, revealing a nem form of nonequilibrium phase dynamics in mixed valence oxides. These results advance the understanding of metastability and memory phenomena in strongly correlated materials, opening pathways toward adaptive and neuromorphic functionalities in quantum materials.

Correlated flat bands and altermagnetism are two important directions in quantum materials, centred respectively on interaction-dominated phases and symmetry-enforced spin-textured states, yet both derive from lattice symmetry and orbital hybridization. This common origin implies that extreme crystal distortion, by narrowing bandwidths, enhancing correlations and reshaping the symmetries of altermagnetic spin splittings, could unify flat-band and altermagnetic physics in a single material; in practice, however, achieving such large distortions in a crystalline altermagnet is a formidable challenge. Here we combine a dedicated strain device with a tailored single-crystal mounting scheme to impose a highly tensile strain gradient in bulk CrSb, a prototypical altermagnet, creating a near-surface layer in which the in-plane lattice is strongly distorted relative to the weakly strained bulk, while the average bulk distortion remains small. Angle-resolved photoemission reveals a reversible regime at moderate strain, where a deeper flat-band feature, attributed to a strain-gradient-driven suppression of Cr-Sb hybridization, coexists with a correlation-enhanced Cr 3d flat band, and an irreversible regime at larger strain where partial bond decoupling drives a predominantly insulating spectral response. Density-functional calculations show that an orbital-selective altermagnetic spin texture persists across this correlated regime despite strong bandwidth renormalisation. These results define a strain-symmetry-correlation map for CrSb and establish extreme tensile strain as a route to co-engineer flat-band tendencies and spin-textured, zero-net-moment correlated states in altermagnets, pointing toward strain-adaptive, spin-selective Mott filtering and related device concepts.

Establishing symmetry-protected topological (SPT) phases with interactions in chemically realistic systems remains an open challenge. We show that a single, synthetically plausible organic one-dimensional chain, tunable via chemical modification of its radical sites, hosts two such phases: an odd-Haldane phase of a dimerized $S=\tfrac{1}{2}$ Heisenberg chain and a Haldane phase of an $S=1$ chain realized when Hund coupling locks two $S=\tfrac{1}{2}$ spins per monomer into $S=1$. Density-functional theory places the active manifold deep in the Mott regime ($U/t\!\approx\!126$), justifying a spin-only Heisenberg description; a compact $(t,U)\!\to\!J$ mapping then fixes exchange couplings. Exact diagonalization and DMRG reveal a consistent SPT fingerprint across both phases, including a quantized many-body Zak phase, even-degenerate entanglement spectrum, protected edge spins, and characteristic triplon/Haldane features in $S^{+-}(q,ω)$. Our results identify a chemically programmable molecular platform for interacting SPT physics in one dimension and suggest concrete spectroscopic routes to organic Haldane spin chains for nanoscale quantum devices.

Research on strongly correlated electron systems faces a fundamental challenge due to the complex nature of intrinsic many-body correlations. A key strategy to address this challenge lies in advancing experimental methods that can directly probe and elucidate the underlying many-body correlations. In this perspective article, we discuss the theoretically proposed coincidence detection techniques, which are designed to directly measure two-body correlations in various particle-particle and particle-hole channels, with momentum, energy, and/or spatial resolution. We also explore the prospects of these coincidence detection techniques for future theoretical and experimental developments. The successful implementation and refinement of these coincidence detection techniques promise to deliver powerful new approaches for unraveling long-standing puzzles in strongly correlated electron systems, such as the enigmatic mechanism of unconventional superconductivity and the long-sought quantum spin liquids. Furthermore, these coincidence detection techniques will offer powerful new methods to investigate novel phenomena like itinerant magnetism and electronic nematicity in quantum materials.

We study energy minimization problems in quantum chemistry through the lens of computational algebraic geometry. We focus on minimizing the Rayleigh quotient of a Hamiltonian over a tensor train variety. The complex critical points of this problem approximate eigenstates of the quantum system, with the global minimum approximating the ground state. We call the number of critical points the Rayleigh-Ritz degree. After introducing tensor train varieties, we identify instances when they are Segre products of projective spaces. We also report what we know about the defining ideals of tensor trains. We present a birational parametrization of them from products of Grassmannians. Along the way, we study the Rayleigh-Ritz degree, and we introduce the Rayleigh-Ritz discriminant, which describes Hamiltonians that lead to deficient number of critical points. We use homotopy continuation to compute all critical points of this optimization problem over various tensor train and determinantal varieties. Finally, we use these results to benchmark state-of-the-art methods, the Alternating Linear Scheme and Density Matrix Renormalization Group.