arXiv Curator

We investigate the superconducting gap structure in bilayer nickelates within a model where conduction bands of dx2-y2 symmetry coexist with localized d3z2-r2 spins. Strong interlayer coupling drives the local moments into a singlet ground state, whose virtual singlet-triplet excitations ("triplons") mediate the pairing interaction between conduction electrons. This yields interband s+- pairing, with opposite signs of the order parameter on the bonding beta and antibonding alpha bands. Our theory naturally explains two key experimental features: a larger gap on the alpha band despite its smaller density of states, and pronounced gap anisotropy arising from momentum-dependent nonlocal Kondo coupling. These results support triplon-mediated pairing as the microscopic origin of superconductivity in bilayer nickelates.

We present a variational quantum eigensolver (VQE) approach for solving the Anderson Impurity Model (AIM) arising in Dynamical Mean-Field Theory (DMFT). Recognizing that the minimal two-site approximation often fails to resolve essential spectral features, we investigate the efficacy of VQE for larger bath discretizations while adhering to near-term hardware constraints. We employ a symmetry-adapted ansatz enforcing conservation of particle number $(N)$, spin projection $(S_z=0)$, and total spin $(S^2=0)$ symmetry, benchmarking the performance against exact diagonalization across different interaction strengths using bath parameters extracted from the DMFT self-consistency loop. For a four-site model, the relative error in the ground state energy remains well below $0.01%$ with a compact parameter set $(N_p \le 30)$. Crucially, we demonstrate that the single-particle Green's function-the central quantity for DMFT-can be accurately extracted from VQE-prepared ground states via real-time evolution in the intermediate to strong interaction regimes. However, in the weak interaction regime, the Green's function exhibits noticeable deviations from the exact benchmark, particularly in resolving low-energy spectral features, despite the ground state energy showing excellent agreement. These findings demonstrate that VQE combined with real-time evolution can effectively extend quantum-classical hybrid DMFT beyond the two-site approximation, particularly for describing insulating phases. While this approach offers a viable pathway for simulating strongly correlated materials on near-term devices, the observation that accurate ground state energy does not guarantee accurate dynamical properties highlights a key challenge for applying such approaches to correlated metals.

We study quasi-one-dimensional strongly correlated materials using a multi-step approach based on density functional theory, downfolding techniques, and tensor-network simulations. The downfolding procedure yields effective multiband Hubbard models that capture the competition between electron hopping and local Coulomb interactions relevant to the system's low-energy properties. The resulting multiband Hubbard models are solved using matrix product states. Applied to Sr$_2$CuO$_3$, SrBaCuO$_3$, and Ba$_2$CuO$_3$, this purely one-dimensional treatment yields no long-range magnetic order, in contrast to the magnetic ordering observed experimentally. To account for this behavior, we extend the multi-step approach by incorporating interchain couplings through a self-consistent mean-field scheme. This combined approach stabilizes finite staggered magnetizations, providing a consistent description of magnetic order in agreement with experiment. For Sr$_2$CuO$_{3.5}$ and SrCuO$_2$, we also tested an approach proposed for ladder materials, however, we find that these materials are not well suited for this approach due to the small magnitude of the intraladder hopping parameters.

Integrating real-space topological spin textures with momentum-space topological electronic states within a single altermagnetic system has remained a persistent challenge. Here, we introduce a symmetry-locked bilayer altermagnet that concurrently hosts d-wave altermagnetism, momentum-space topology, and stable antiskyrmions. In momentum-space, it enables strain-triggered transitions to an antiferromagnetic Weyl semimetal phase, where the Néel vector acts as a switch for spin-layer-polarized quantum anomalous Hall and Weyl states, alongside coupled topological states and valley polarization effects. In real-space, the formation of interlayer co-directional and locked in-plane Dzyaloshinskii-Moriya interactions facilitates the creation of coupled antiskyrmion pairs with compensated topological charges. This locking symmetry fully cancels the transverse Magnus forces, resulting in current-driven, strictly longitudinal motion of antiskyrmions without any Hall-like deflection. Our work establishes a robust platform for dual-space topological magnetism and offers a definitive solution to the vanishing Hall angle in charge-neutral antiskyrmions, opening pathways toward high-density, low-power topological spintronics.

We study the proximity effect in the Hubbard model coupled to BCS superconductors describing a single-layer strongly correlated electron system in a phase-biased Josephson junction. We find two distinct gapped solutions, one Mott-like insulating (M-phase) and one proximitized superconducting phase (S-phase), separated by first-order transition with hysteresis. In the M-phase the large correlation charge gap strongly suppresses the critical current, while the S-phase behaves as a $0$-junction, with a proximitized gap that closes for $φ=π$ to yield a correlated metal. Phase bias and junction transparency can thus serve as tuning knobs to switch between conducting and insulating regimes. Working within the dynamical mean field theory using the numerical renormalization group as the impurity solver, we associate M- and S-phase solutions with the doublet and singlet fixed points of the underlying superconducting Anderson impurity problem. We obtain detailed insight into the spectral structure on all energy scales. In the M-phase, the self-energy has sub-gap resonances symmetrically located around the Fermi level resulting from the splitting of the ''mid-gap pole'' found in Mott insulators; this structure accounts for phase insensitivity.

Dimensionality of a physical system, conventionally an invariant geometric characteristic, fundamentally governs the universality class of phase transitions and the landscape of emergent collective phenomena. In low-dimensional or layered high-temperature superconductors, the macroscopic phase coherence of superconducting orders is typically confined in two dimensions, underscoring the critical role of phase fluctuations in determining the overall phase diagrams. Here, we strategically enhance the phase fluctuations by fabricating periodically arranged nano-holes in the infinite-layer nickelate superconducting films, effectively constructing Josephson junction arrays. In the nano-patterned films, the weakening of macroscopic phase coherence drives a two-stage superconducting transition towards an anomalous metallic ground state with saturated resistance. The emergence of charge-2e quantum oscillations manifests the coherence across the array, while an anomalous zero-field magnetoresistance peak signifies the extreme quantum phase fluctuations persisting to ultralow temperatures. Remarkably, with quantum fluctuations enhanced synergistically by nano-patterning and magnetic fields, an anomalous reversal of superconducting anisotropy is observed in Nd-nickelates, where in-plane critical fields fall below out-of-plane values. The evolution of anisotropy may unmask an internal exchange-Zeeman field coupled to the collective electronic states. Our results unveil how superconductivity evolves in response to phase fluctuations, establishing nano-patterning as a powerful paradigm to uncover hidden intertwined orders in strongly correlated systems.

MnTe has recently attracted exceptional attention due to its well-established altermagnetism, prompting a thorough reexamination of its properties. In particular, it was found that a Raman-active excitation at ~175 cm$^{-1}$, routinely assigned to the E2g phonon, is incompatible with this interpretation. It was further hypothesized that this mode is a "leakage", due to symmetry lowering, of an otherwise forbidden phonon. Here, using first-principles calculations, we decisively rule out this hypothesis and propose an alternative interpretation that the "mystery mode" is an electronic excitation, i.e., a plasmon, enabled by hole self-doping. The resolution of this mystery will require additional experiments and shed new light on the nature of electronic transport in MnTe.

Consecutive topological phase transitions (TPTs) between strongly correlated electronic phases that differ simultaneously in symmetry breaking and topological order are of fundamental interest in condensed matter physics, yet are rarely realized experimentally. We report two consecutive electric-field-driven TPTs at half filling (nu = 1) in angle-aligned MoTe2/WSe2 moire heterobilayers. With increasing out-of-plane displacement field, a geometrically frustrated Mott insulator evolves into a ferromagnetic quantum anomalous Hall (QAH) Mott insulator, i.e., a spin-polarized topological Mott insulator without an observable charge-gap closure, and subsequently into an antiferromagnetic, valley-coherent Mott insulator (VC-AFM) accompanied by a continuous charge-gap collapse and the emergence of a critical metallic state. Layer-resolved magnetic circular dichroism (MCD), magneto-transport, and compressibility measurements jointly determine the phase diagram. The high-field evolution of the antiferromagnetic state reveals a metamagnetic-like transition at a critical field B*, above which a Chern insulating transport response reappears. Our results establish the MoTe2/WSe2 moire platform as a tunable realization of an extended Kane-Mele-Hubbard model hosting sequential correlation-topology-intertwined transitions.

Altermagnetism arises from composite real-space and spin-space symmetries, combining zero net magnetization with pronounced momentum-dependent spin splitting. This review highlights the pivotal role of angle-resolved photoemission spectroscopy (ARPES), along with its spin-resolved (SARPES) and circular-dichroism (CD-ARPES) variants, in directly visualizing nonrelativistic band splitting and spin textures in altermagnets. Within the spin-group framework, we distinguish ferromagnetic, antiferromagnetic, and altermagnetic orders and elucidate the symmetry origin of spin polarization. We then systematically review representative systems: the debated d-wave prototype RuO2, layered d-wave altermagnets KV2Se2O and Rb1-delta V2Te2O, and a series of g-wave compounds including MnTe (domain-tunable) and CrSb (topological), together with the noncoplanar antiferromagnet MnTe2 and other emerging candidate platforms. Overall, ARPES has become a key probe for resolving symmetry-driven spin splitting. Future advances in micro/nano-beam and in-situ spectroscopies, combined with strain and domain engineering, heterostructure design, and exploration of broader unconventional magnetic states, are expected to drive the joint evolution of altermagnetism and photoemission spectroscopy, paving the way for spintronic and correlated quantum research.

We study superconductivity and magnetism in bilayer nickelates from an itinerant perspective. Starting from a tight binding fit to recent ARPES measurements on compressively strained thin films, we incorporate the standard set of onsite repulsive interactions among partially filled $e_g$ orbitals: intra-orbital $U$, inter-orbital $U'$, Hund's coupling $J_H$ and a pair hopping $J_P$. We obtain the effective pairing interaction by dressing these bare interactions with particle-hole fluctuations via the RPA. In the strong Hund's coupling regime, we find that $s$-wave superconductivity and $(π/2, π/2)$ SDW order are the favored ground states. With weaker Hund's coupling, we find that $d$-wave pairing and $(π, π)$ SDW are the leading ground states. Our results are qualitatively consistent with earlier DMRG studies, and point to the key role played by Hund's coupling in determining the nature of superconductivity and magnetism in this system.

We study superconductivity and magnetism in bilayer nickelates from an itinerant perspective. Starting from a tight binding fit to recent ARPES measurements on compressively strained thin films, we incorporate the standard set of onsite repulsive interactions among partially filled $e_g$ orbitals: intra-orbital $U$, inter-orbital $U'$, Hund's coupling $J_H$ and a pair hopping $J_P$. We obtain the effective pairing interaction by dressing these bare interactions with particle-hole fluctuations via the RPA. In the strong Hund's coupling regime, we find that $s$-wave superconductivity and $(π/2, π/2)$ SDW order are the favored ground states. With weaker Hund's coupling, we find that $d$-wave pairing and $(π, π)$ SDW are the leading ground states. Our results are qualitatively consistent with earlier DMRG studies, and point to the key role played by Hund's coupling in determining the nature of superconductivity and magnetism in this system.

Dynamical mean-field theory (DMFT) provides an optimal local approximation for correlated lattice systems by mapping the lattice onto a self-consistent effective impurity model. To account for the missing long-range correlations, we propose a novel extended approach, which we term fluctuating dynamical mean-field theory (fDMFT). It incorporates collective fluctuations of auxiliary impurity models across different sites via functional integration. Technically, this method involves obtaining a family of DMFT solutions on a grid for a self-consistent auxiliary classical field applied to the lattice. While the result can, in principle, be improved diagrammatically, we find that the minimal version of the theory already yields accurate results, with lowest-order diagrammatic corrections offering only minor improvements. This consistent framework, based on our fluctuating local field concept, demonstrates superior performance for the nearly half-filled Hubbard model compared to other known diagrammatic extensions of DMFT.

Progress toward a complete understanding of cuprate superconductors has been hindered by their intricate phase diagram, potentially linked to a quantum critical point (QCP). However, conclusive evidence for the QCP is lacking, as the presumed QCP is buried under the superconducting dome, disguising its presence. Here, we use high-resolution resonant inelastic X-ray scattering to examine the dynamical charge-charge correlation in La$_{2-x}$Sr$_x$CuO$_4$ and uncover the quantum critical scaling, a key feature required for a QCP. Specifically, \djh{we observed that the inverse correlation lengths for various dopings and temperatures collapsed onto a universal scaling curve, yielding a critical exponent $ν$ of $0.74 \pm 0.08$. The non-negativity of this exponent confirms the presence of a QCP. Remarkably, the value of $ν$ suggests that while the QCP is manifested through the charge-density wave, other orders also participate, such that the QCP appears to belong to the universality class characterized by the O(4) symmetry, reminiscent of the microscopic SO(4) symmetry in the Hubbard model at half-filling. Further analysis indicates that the QCP is highly dissipative with a short quasi-particle lifetime, reflecting the intertwined quantum fluctuations due to its being buried inside the superconducting state.

The interplay between magnetic fluctuations and superconductivity is fundamental for understanding unconventional high-temperature superconductors. In the recently discovered Ruddlesden-Popper phase nickelates-which achieve superconducting transition temperatures up to ~100 K-this connection has been theoretically predicted but experimentally unverified. Using compressively strained La2SmNi2O7 thin films, we report direct evidence for this interplay. We observe a characteristic 'Mexican hat'-shaped magnetoresistance, a signature of superconductivity coexisting with emergent antiferromagnetic (AFM) fluctuations. This distinctive feature arises from a competition between the magnetic field's suppression of AFM fluctuation-driven scattering and its suppression of superconducting fluctuations. We quantify these AFM fluctuations by a crossover field, B*, whose magnitude decreases with increasing temperature and vanishes near the superconducting onset-behavior contrasting sharply with that in high-Tc cuprates. Our results establish not merely coexistence, but also a direct and novel correlation between AFM instability and superconductivity in nickelates, offering crucial insight into their pairing mechanism.

Electron-phonon coupling (EPC) governs lattice dynamics, charge transport, and collective electronic phases in quantum materials. In several families of unconventional superconductors, including transition-metal dichalcogenides and kagome metals, growing evidence points to a cooperative role of EPC and dynamic charge-density fluctuations (CDF) in stabilizing superconductivity. However, how the EPC strength evolves across phase diagrams and relates to superconducting properties in strongly correlated systems remains an open question. Here we investigate the interplay between phonons and the CDF recently identified in cuprate superconductors. Using resonant inelastic x-ray scattering, we track the dispersion and intensity of bond-stretching phonons in YBa$_2$Cu$_3$O$_{7-δ}$ over wide ranges of doping, temperature, and momentum. We find that both the phonon softening at the CDF wave vector and the EPC strength, extracted from a pronounced phonon intensity anomaly, are maximized near $p = 0.19$, where superconducting properties are optimal and CDF intensity is strongest. These results identify dynamic charge-density fluctuations, rather than quasi-static charge density waves, as the dominant source of phonon renormalization in cuprates, and establish a direct correlation between EPC strength and the superconducting dome. More broadly, our measurements highlight EPC as a doping-dependent property of correlated materials, shaped by the electronic environment in which lattice vibrations are embedded.

We perform a renormalized mean-field study of the two-dimensional repulsive Hubbard model, focusing on the intricate interplay and possible coexistence of magnetic, charge, and superconducting orders. We improve on conventional mean-field theory by utilizing a renormalization group framework that captures high-energy fluctuations. This method generates effective magnetic and $d$-wave pairing interactions, and allows for an unbiased exploration of coexisting phases at weak and moderate interaction strengths. Unrestricted mean-field calculations of the effective Hamiltonian on large finite lattices are combined with analyses in the thermodynamic limit, revealing a rich phase diagram with extensive regions of coexisting orders. We find that $d$-wave superconductivity coexists with Néel order on the electron-doped side. On the hole-doped side, superconductivity is found to coexist with spiral or stripe magnetic orders. Within the stripe ordered region, the superconducting order parameter is spatially modulated, with a period that follows the charge modulation of the stripes. Below van Hove filling, pairing provides the primary energy gain, while the stripe order yields only a small, and hence fragile, additional energy lowering.

Motivated by the recent discovery of superconductivity in the bilayer nickelate La$_3$Ni$_2$O$_7$ (LNO) under pressure, we study a mixed-dimensional (mixD) bilayer $t_\parallel$-$J_\parallel$-$J_\perp$ model, which has been proposed as an effective low-energy description of LNO. Using neural quantum states (NQS), and in particular Gutzwiller-projected Hidden Fermion Pfaffian State, we access the ground-state properties on large lattices up to $8\times 8\times 2$ sites. We show that this model exhibits superconductivity across a wide range of dopings and couplings, and analyze the pairing behavior in detail. We identify a crossover from tightly bound, Bose-Einstein-condensed interlayer pairs at strong interlayer exchange to more spatially extended Bardeen-Cooper-Schrieffer-like pairs as the interlayer exchange is decreased. Furthermore, upon tuning the intralayer exchange, we observe a sharp transition from interlayer $s$-wave pairing to intralayer $d$-wave pairing, consistent with a first-order change in the pairing symmetry. We verify that our simulations are accurate by comparing with matrix product state simulations on coupled ladders. Our results represent the first simulation of a fermionic multi-orbital system with NQS, and provide the first evidence for superconductivity in two-dimensonal bilayers using high-precision numerics. These findings provide insight into superconductivity in bilayer nickelates and cold atom quantum simulation platforms.

Simulating quantum many-body systems (QMBS) is one of the long-standing, highly non-trivial challenges in condensed matter physics and quantum information due to the exponentially growing size of the system's Hilbert space. To date, tensor networks have been an essential tool for studying such quantum systems, owing to their ability to efficiently capture the entanglement properties of the systems they represent. One of the well-known tensor network architectures, namely matrix product states (MPS), is the standard method for simulating one-dimensional QMBS. Here, we propose a simple, yet efficient, method to augment the already available MPS algorithms to simulate the dynamics of time-dependent Hamiltonians with better accuracy and a faster convergence rate, giving a second-order convergence compared to the first-order convergence of the standard method. We apply our proposed method to simulate the dynamics of a chain of single spins associated with nitrogen-vacancy color centers in diamonds, which has potential applications for practical and scalable quantum technologies, and find that our method improves the average error for a system of few NV centers by a factor of about 1000 for moderate step sizes. Our work paves the way for efficient simulation of QMBS under the influence of time-dependent Hamiltonians.

The evolution of magnetic correlation in strongly correlated electron systems with altermagentic spin splitting remains largely unexplored. Here we investigate how spin splitting generated by spin-dependent next-nearest-neighbor hopping $t'$ reshapes the Fermi surface nesting and van Hove singularities in the two-dimensional square-lattice Hubbard model, leading evolution of magnetic instabilities. Using the constrained-path quantum Monte Carlo method, we find the dominant magnetic correlation as functions of the filling and $t'/t$ by computing the momentum-resolved spin structure factor. The analysis reveals a transition from antiferromagnetic $(π,π)$ order in the isotropic, half-filled system to non-collinear spiral $(π,q)$ order upon increasing the spin-dependent anisotropy or doping away from half-filling, ultimately entering a short-range correlation regime where stripe and spiral correlation coexist. These findings highlight a possible route to realizing spiral correlation in altermagnetic systems, potentially providing a platform for spintronic devices that exploit non-collinear spin textures.

The interplay among orbital-selective Mott physics, Hund's coupling, tunable structural motifs, and Kondo-like scattering establishes a compelling paradigm for understanding and engineering correlated multi-orbital systems, as vividly exemplified by nickelate superconductors. Here, using high-resolution angle-resolved photoemission spectroscopy combined with theoretical calculations, we systematically investigate the electronic properties of trilayer nickelates. In La4Ni3O10, we observe pronounced interorbital hybridization, whereas in Pr4Ni3O10, the flat d_(z^2 ) band becomes markedly incoherent and diminishes in spectral weight. By contrast, the dispersive d_(x^2-y^2 ) bands retain coherence in both compounds. This striking incoherence/coherence dichotomy identifies an orbital-selective Mott phase modulated by the interlayer Ni-O-Ni bonding angle. The depletion of the d_(z^2 ) orbitals further frustrates the interorbital hybridization and influences the density-wave transition in Pr4Ni3O10. Moreover, the density-wave gap is substantially reduced in Pr4Ni3O10, likely due to extra scattering channels provided by the local moments of Pr3+ cations. Our findings elucidate the intricate interplay among lattice, orbital, spin, and electronic degrees of freedom and reveal a feasible structural control parameter for the multi-orbital correlated state in trilayer nickelates, which provide a concrete framework for understanding the emergence of superconductivity under high pressure.

Unraveling the interplay between density-wave (DW) instabilities and multi-orbital physics is critical for understanding superconductivity in Ruddlesden-Popper nickelates, yet intrinsic electronic features have been persistently obscured by material inhomogeneity and thus the multi-domain averaging effect. Here, we employ micro-focused angle-resolved photoemission spectroscopy ($μ$-ARPES) on single-domain Pr$_4$Ni$_3$O$_{10}$ to disentangle the complex hierarchy of intrinsic and back-folded bands, explicitly identifying the electronic states driving the DW phase transition. We provide decisive spectroscopic evidence that the low-energy reconstruction is governed by inter-orbital nesting between the $α$ and $β$ bands. Specifically, we resolve a orbital-dependent gap of $\sim44$ meV on the $α$ pocket, a value quantitatively consistent with prior measurements, unifying previously conflicting experimental reports regarding the locus and magnitude of the DW gap. Furthermore, we reveal strong orbital-selective mass renormalization in the $d_{z^2}$ states and successfully resolve the long-sought intrinsic trilayer $β$-band splitting, establishing a critical lower bound for the outer-layer hopping. These results define a coherent microscopic fingerprint for the trilayer nickelates, identifying the specific nesting channels and correlation effects that underpin the phase diagram.

Inspired by recent discovery of metallic $d$-wave altermagnetism in KV$_2$Se$_2$O, we develop a self-consistent microscopic many-body calculation of density-wave order for an itinerant altermagnetic metal. We show that the strong $d$-wave spin-momentum locking inherent to the altermagnetic band structure reconstructs the Fermi surface into spin-selective quasi-1D open sheets. This unique topology of Fermi surface drives an instability toward spin-resolved electronic charge density waves (CDWs), in which the ordering wave vectors for spin-up and spin-down electrons condense along two mutually orthogonal directions, forming spin-resolved stripe phases. As a consequence, this results in pronounced gap openings near the Fermi surface, and the superposition of these spin-resolved stripe orders leads to a checkerboard CDW in the charge channel and an antiphase spin-density-wave modulation in the spin channel. Upon increasing temperature, the density-wave order melts at $T_c$ due to thermal phase fluctuation while the gap opening persists, giving rise to a robust pseudogap regime, which eventually closes at a higher temperature $T_g$. The resulting simulations quantitatively reproduce the key features observed in the spectroscopic measurements, offering a consistent and generic understanding of the reported phenomena in KV$_2$Se$_2$O and, more broadly, in metallic altermagnets with strong spin-momentum locking.

Spin-orbital liquids provide an exactly solvable route to three-dimensional Z2 quantum spin liquids beyond the original Kitaev setting. Built from higher-dimensional Clifford-algebra representations, spin-orbital Hamiltonians can be realized on both three- and four-coordinated lattices, giving rise to phases with 3 and 2 itinerant Majorana flavors. We demonstrate that these models host a rich set of gapless Majorana metals, characterized, in particular, by topological Fermi surfaces, nodal lines, and Weyl semimetal phases. We analyze the stability of these structures under physically motivated perturbations and identify generic splitting patterns and topological transitions driven by symmetry breaking and flavor mixing. This yields a unified organizing framework for three-dimensional Majorana metals in fractionalized spin liquids.

The La2-xSrxCuO4+delta (LSCO) system provides a unique experimental setting for exploring how magnetism, superconductivity, and disorder jointly shape charge transport in a doped Mott insulator. Transport measurements in lightly doped and oxygen-enriched LSCO reveal a strongly insulating normal state governed by variable-range hopping, accompanied by pronounced nonlinear current-voltage characteristics and, at low temperatures, current-induced negative differential resistance. With increasing carrier concentration, these features evolve into regimes characterized by granular and percolative superconductivity near the threshold of bulk superconductivity and, eventually, into a homogeneous strange-metal state close to optimal doping. Throughout this evolution, the transport response shows marked sensitivity to disorder, electronic inhomogeneity, and external control parameters, such as bias current and magnetic field. Rather than reflecting a sequence of sharply distinct phases, the observed transport regimes form a continuous crossover from a localization-dominated insulating state to granular superconductivity and further to a coherent metallic state. This crossover is driven primarily by the progressive enhancement of electronic screening, inter-region coupling, and superconducting connectivity, rather than by abrupt changes in the underlying microscopic scattering mechanisms. Taken together, the available transport data provide a coherent experimental basis for understanding how disorder and mesoscale electronic inhomogeneity organize charge transport and superconductivity across the LSCO phase diagram, underscoring the central role of percolation and nonequilibrium effects in underdoped cuprates.

Altermagnets emerge as a novel platform for realizing unconventional superconductivity through their exotic momentum-dependent spin-splitting of electronic band structures. Recent experiments have uncovered a novel form of altermagnetism with distinctive $g$-wave symmetry in CrSb. However, the potential for unconventional superconductivity arising from $g$-wave altermagnetism in such systems remains largely unexplored. In this study, we discover the emergence of chiral superconducting states in three-dimensional $g$-wave altermagnetic metals. Through systematic self-consistent mean-field analysis on the extended attractive Hubbard model combined with $g$-wave altermagnetic exchange fields in a three-dimensional hexagonal lattice, as observed in CrSb, we find that the altermagnetic spin splitting of Fermi surfaces favors chiral $p$-wave states as the dominant pairing channel under strong altermagnetic fields and high electron densities, while chiral $d$-wave states become predominant under weak altermagnetic fields and intermediate electron densities. Conversely, at weak altermagnetic fields and typical electron densities, non-chiral $s$-, extended $s$-, or $f$-wave states become stabilized. We also showcase the possible experimental detection using the quasiparticle energy dispersions and the density of states to distinguish different pairing symmetries. These findings underscore the potential of $g$-wave altermagnets to host sought-after chiral and gapless superconductivity.

Quantum circuits have been shown to be a fertile ground for realizing long-range entangled phases of matter. While various quantum double models with non-chiral topological order have been theoretically investigated and experimentally implemented, the realization and characterization of chiral topological phases have remained less explored. Here we show that chiral topological phases in spin systems, i.e., chiral spin liquids, can be prepared in quantum circuits using the variational quantum eigensolver (VQE) framework. On top of the VQE ground state, signatures of the chiral topological order are revealed using the recently proposed tangent space excitation ansatz for quantum circuits. We show that, both topological ground state degeneracy and the chiral edge mode can be faithfully captured by this approach. We demonstrate our approach using the Kitaev honeycomb model, finding excellent agreement of low-energy excitation spectrum on quantum circuits with exact solution in all topological sectors. Further applying this approach to a non-exactly solvable chiral spin liquid model on square lattice, the results suggest this approach works well even when the topological sectors are not exactly known.

One of the long standing puzzles in strongly correlated materials is the microscopic origin of the quantum critical Planckian strange metal phase with universal linear in temperature scattering rate from which unconventional superconductivity directly emerges by lowering temperatures. Recently, the linear in temperature and linear in field resistivity have been simultaneously observed in high temperature cuprate superconductors, manifested by the universal field to temperature scaling in magnetoresistivity. To date, there has been a lack of coherent and unified understanding of these coexisting linear behaviors and their possible link to quantum criticality. In this work, we establish the universality in linear in temperature and linear in field Planckian behaviors in underdoped LSCO near optimal doping. Experimentally, we observe the linear in field Planckian scattering rate and its relation to its linear in temperature counterpart. Theoretically, we propose a spin based common microscopic mechanism based on Kondo-like charge fluctuations near local quantum criticality of heavy fermion formulated tJ model subject to a Zeeman term. Similar to frequency to temperature scaling near quantum criticality, we find the magnetic field here effectively introduces a Zeeman energy, reminiscent of an external energy in the quantum critical regime, leading to field to temperature scaling. Our analytically predicted universal field to temperature scaling in isotropic scattering rate and the relation between the linear in temperature and linear in field Planckian coefficients, unifies these two phenomena over an extended doping range, pointing toward a unified quantum-critical origin of Planckian transport in cuprates.

Despite being well established in cuprates, an intrinsic translational symmetry-breaking charge order has not been clearly identified in iron-based superconductors. Using spectroscopic-imaging scanning tunneling microscopy on epitaxial Ca(Fe1-xCox)2As2 (x = 0 ~ 0.055) thin films, we observe smectic, near-commensurate charge-stripe order in the underdoped regime that intervenes between the nematic parent phase and optimally doped superconductivity. Distinct from the bidirectional checkerboard-like order in cuprates, these charge stripes are unidirectional along the antiferromagnetic Fe-Fe bond direction and are accompanied by a van Hove singularity near the Fermi level, inherited from the Fermi surface reconstruction driven by intertwined antiferromagnetic and nematic correlations. Both local and global suppression of the charge-stripe instability enhance superconductivity, tunable via epitaxial strain and Co doping. These results establish charge-stripe order as an intermediate electronic phase in iron pnictides and reveal a coherent pathway from nematicity to superconductivity. Our findings highlight charge ordering as a unifying element across different families of high-temperature superconductors.

We investigate superconductivity emerging in the photodoped Mott insulating Hubbard model using steady-state dynamical mean-field theory implemented on the real-frequency axis. By employing high-order strong-coupling impurity solvers, we obtain the nonequilibrium phase diagram for photoinduced $η$-pairing superconductivity with a remarkably high effective critical temperature. We further identify a superconducting gap in the momentum-resolved spectral function and optical conductivity, providing spectroscopic signatures accessible to experiments. Our results highlight a route to a controllable form of high-temperature superconductivity in nonequilibrium strongly correlated systems, fundamentally distinct from the equilibrium $s$-wave pairing state in the attractive Hubbard model or cuprate-like $d$-wave superconductors.

Quantum geometry quantifies how the single-particle Bloch wavefunction changes in phase and amplitude across the Brillouin Zone. In multi-orbital systems where bands have strongly mixed orbital composition, quantum geometry plays a vital role in determining the ground state and low-energy properties of interacting electronic systems. In this work, we show that Mott metal-insulator transitions, as well as transitions between different magnetic orders within the Mott insulating phase, can be driven by the quantum geometry of the underlying Bloch band, thereby providing a mechanism complementary to conventional bandwidth-tuned Mott transitions. By studying the Kane-Mele-Hubbard model using exact diagonalization, we demonstrate that in in half-filled and topologically-trivial bands, quantum geometric properties of the Bloch states alone can act as a tuning knob for Mott metal-to-insulator and affect the competition between ferromagnetism and antiferromagnetism. We show that both transitions may be heuristically understood via non-local Coulomb scattering in a basis of exponentially localized Wannier functions. These results highlight the role of quantum geometry beyond topological settings as a governing principle for conventional Mott and magnetic physics in multi-orbital and moiré materials.

Exact formulas for the optical conductivity and the Hall conductivity of the two dimensional Hubbard model are derived in terms of an effective theory description of the local moment fluctuation in the system. In this framework, the quantum Monte Carlo simulation of the electromagnetic response of such a strongly correlated electron system becomes sign-problem-free in many physically relevant cases. In particular, it is sign-problem-free when we assume the widely used Millis-Monien-Pines form for the phenomenological susceptibility in the effective action of the fluctuating local moment, even though these local moments are now subjected to Landau damping as a result of their coupling to the itinerant quasiparticle on the fermi surface. This is true more generally when a $\varphi^{4}$ term is included in the effective action and is thus not restricted to the Gaussian limit. Here we demonstrate the power of this framework by studying the effect of thermal fluctuation of the local moment on the optical conductivity $σ^{xx}(ω)$ and the Hall conductivity $σ^{xy}(ω)$ of the cuprate superconductors. Both $σ^{xx}(ω)$ and $σ^{xy}(ω)$ calculated are found to exhibit a two-component structure, with a Drude component at low energy and a mid-infrared component at higher energy. Depending on the relative importance of the hole pocket and the electron pocket on the reconstructed fermi surface and the coupling strength to the local moment, the Drude component in $\mathrm{Im}σ^{xy}(ω)$ can be either positive or negative.(full-length abstract can be found in the main text.)

Superconductivity is usually understood as a phase in which charge-$2e$ Cooper pairs are condensed. Charge-$4e$ superconductivity has largely been discussed as a vestigial order at finite temperature emerging from charge-$2e$ states. Primary charge-$4e$ superconducting phases at zero temperature remain scarce in both experiments and microscopic models. Here we argue that a doped featureless Mott insulator with $SU(4)$ symmetry provides a natural platform for primary charge-$4e$ superconductivity, based on perturbative renormalization group arguments and group theoretic considerations. As a concrete realization, we construct a bilayer Hubbard model with tunable onsite $SU(4)$ and $Sp(4)$ symmetries that exhibits a featureless Mott insulating phase at half filling. Its low energy physics is captured by a generalized ESD model, featuring an effective Hamiltonian that is purely kinetic within the constrained Hilbert space. Using density matrix renormalization group (DMRG) simulations, we find a primary charge-$4e$ superconducting phase in the $SU(4)$ ESD model and a conventional primary charge-$2e$ phase in the $Sp(4)$ case. We further characterize the corresponding normal states and discuss the resulting finite temperature phase diagram.

Altermagnetism (AM) and its associated spin-transport phenomena are typically linked to spin-split electronic band structures in bulk materials. However, the crystal surface has a reduced symmetry with respect to the bulk, which can induce AM at the surface of conventional antiferromagnets (AFMs) -- a local effect which cannot be detected using bulk properties. In this work we define the symmetry conditions necessary for surface AM and show how it can be topologically protected, rendering it a robust effect. We provide a minimal model for one trivial and two topological examples of surface AM. We show that the spin spectral density, accessible by spin- and angle-resolved photoemission spectroscopy, can exhibit a $d$-wave-like altermagnetic character at the surface, even when the full band structure is completely spin degenerate. Our topological model describes the Dirac semimetal CuMnAs, which provides an existing realization of our theory. Our results identify crystal surfaces as a platform to realize robust, topology- and symmetry-driven unconventional magnetism beyond the bulk classification of magnetic materials.

The pressurized Shastry-Sutherland Mott insulator SrCu2(BO3)2 has been found to host a plaquette-singlet phase and an antiferromagnetic phase that break different symmetries spontaneously.The recent experiment showed that their transition is of a first order nature, which seems against the pursuit of exotic and deconfined degrees of freedom in this famous frustrated quantum magnet. We found a new direction in this study. By applying a magnetic field to the material, we discover that SrCu2(BO3)2 exhibits a universal and metallic T-linear specific heat behavior in a large magnetitic field range close to the pressure of zero-field first order transition between plaquette-singlet and antiferromagnetic phases. Such an unexpected gapless response from an electronically gapped Mott insulator could be attributed to magnetized Dirac spinons liberated by the combined effect of magnetic field and pressure, consistently seen from our quantum many-body thermal tensor network computation of the Shastry-Sutherland model under magnetic field. Such a robust and universal T-linear specific heat phase points out the richness of the phase diagram of the material expanded by the axes of pressure and magnetic field and is calling for new theoretical frameworks to its full explanation.

The strange metal remains one of the great unsolved problems for 21st century science. Since the early development of the marginal Fermi liquid phenomenology, it has been clear that progress requires detailed knowledge of the momentum- and frequency-dependent charge susceptibility, $χ(\mathbf{q},ω)$, particularly at large momenta. Electron energy-loss spectroscopy (EELS), performed in either reflection or transmission geometry, provides the most direct probe of $χ(\mathbf{q},ω)$. However, measurements over the past four decades have yielded conflicting results, with some studies reporting a dispersing RPA-like plasmon and others observing a strongly overdamped, incoherent response. Here we report a transmission EELS study of Bi$_2$Sr$_2$CaCu$_2$O$_{8+x}$ (Bi-2212) that simultaneously achieves high energy resolution ($ΔE \approx 30$ meV) and high momentum resolution ($Δq \approx 0.01$ Å$^{-1}$). To address issues of reproducibility, measurements were repeated ten times on five different Bi-2212 flakes, benchmarked against aluminum, a well-characterized Fermi liquid, and quantitatively compared with prior studies spanning four decades. At momenta $q < 0.15$ Å$^{-1}$, we observe a highly damped plasmon whose linewidth is comparable to its energy. At larger momenta, $q > 0.15$ Å$^{-1}$, this excitation does not disperse but instead evolves into an incoherent continuum, with no evidence for the RPA-like dispersion reported in some earlier works. Comparison with recent RIXS measurements on Bi-based cuprates supports the view that Bi-2212 is an incoherent metal with strongly damped charge excitations.

Motivated by the recent experimental discovery of the $T$-linear specific heat in pressurized and magnetized Shastry-Sutherland Mott insulator SrCu$_2$(BO$_3$)$_2$, we perform the state-of-the-art thermal tensor-network computation on the Shastry-Sutherland model under a magnetic field. Our simulation results suggest the existence of a symmetric intermediate phase with $T$-linear specific heat at low temperature, occupying a large parameter space and separating the plaquette-singlet phase and antiferromagnetic phase at low fields and other symmetry-breaking phases at high fields before the system is fully polarized. Such an unexpected novel state bears an astonishing similarity to the experimental findings in the material. It opens the door to further investigations of the possible liberation of deconfined magnetized Dirac spinons by the competing interactions in this highly frustrated quantum magnet model, and by the combined effects of magnetic field and pressure in the the associated Shastry-Sutherland Mott insulator SrCu$_2$(BO$_3$)$_2$.

We have investigated the electronic structure of Ba4Ir3O10 within the density-functional theory (DFT) using the generalized gradient approximation while considering strong Coulomb correlations (GGA+U) in the framework of the fully relativistic spin-polarized Dirac linear muffin-tin orbital band-structure method. Ba4Ir3O10 has a quasi-2D structure composed of buckled sheets, which constitute corner-connected Ir3O12 trimers containing three distorted face-sharing IrO6 octahedra. The Ir atoms are distributed over two symmetrically inequivalent sites: the center of the trimer (Ir1) and its two tips (Ir2). The Ir1 - Ir2 distance within the trimer is quite small and equals to 2.58 A at low temperature. As a result, the clear formation of bonding and antibonding states at the Ir1 site occurs. The large bonding-antibonding splitting stabilizes the dyz-orbital-dominant antibonding state of t2g holes and produces a wide energy gap at the Fermi level. However, the energy gap opens up only with taking into account strong Coulomb correlations at the Ir2 site. Therefore, we have quite a unique situation when the insulating state is driven by both the dimerization at the Ir1 site and Mott insulating behavior at the Ir2 one. We have investigated resonant inelastic x-ray scattering (RIXS) spectra at the Ir L3 edge. The calculated results are in good agreement with experimental data. The RIXS spectrum possesses several sharp features below 2.1 eV corresponding to transitions within the Ir t2g levels. The excitation located from 2.1 to 4.6 eV is due to t2g to eg and O2p to t2g transitions. The wide structure situated at 6.2-12 eV appears due to charge transfer and O2p to eg transitions. We have also presented comprehensive theoretical calculations of the RIXS spectrum at the oxygen K edge.

Determining the Fermi surface (FS) volume in underdoped cuprates is crucial for understanding the nature of the strongly correlated pseudogap phase. Conventional quantum oscillation techniques, typically used for this purpose, are inapplicable in this high-temperature regime due to thermal and disorder-induced smearing of Landau levels. We propose Sondheimer oscillations (SO), semiclassical oscillations of in-plane magnetoresistivity in thin films, as a robust alternative probe of FS reconstruction. SO arise from the commensuration between the cyclotron radius and film thickness, do not rely on Landau quantization, and remain observable at moderate fields and elevated temperatures where quantum oscillations are suppressed. Their frequencies depend solely on the FS parameters (e.g., curvature), and not on specific details of scattering mechanisms. SO are also sensitive to the coherence of inter-layer tunneling, allow contributions from individual FS pockets to be distinguished in the frequency domain, and naturally include the Yamaji angle effect (if present in the system) as a prominent feature in the frequency spectrum. We compute SO spectra as a function of the magnetic field orientation for three representative scenarios: (i) an unreconstructed large FS, (ii) a spin density wave reconstructed FS with volume $p/4$, and (iii) a fractionalized Fermi liquid (FL$^*$) with pocket volume $p/8$ (here $p$ is the hole doping). We show that the SO spectrum offers a wealth of universal features that could be used to differentiate between these scenarios. In particular, we highlight a FS geometry-dependent phase shift between oscillations in longitudinal and transverse conductivities, characterize how the FS curvature can be extracted from SO if the film orientation is perpendicular to the crystallographic $c$-axis, and analyze the evolution of the SO spectrum with doping.

Two widely-used non-local extensions of dynamical mean field theory (DMFT), cellular DMFT (CDMFT) and the dynamical cluster approximation (DCA), both yield self-energies marred by having some unphysical properties: CDMFT yields real-space self-energies that are not translationally invariant, and DCA yields momentum-space self-energies with discontinuities in their momentum dependence. It is often desirable to remove these flaws by post-processing cluster DMFT results, using strategies called periodization for CDMFT and interpolation for DCA -- for brevity, we refer to both cases as interpolation. However, traditional interpolation approaches struggle to capture intricate structures such as hole pockets in the hole-doped square-lattice Hubbard model, as highlighted in Phys. Rev. B 105, 35117 (2022). Further, these approaches interpolate frequency-dependent functions, which may lead to causality violations. Here, we propose Liouvillian interpolation, a novel, intuitive, and robust scheme for interpolating cluster DMFT results. Our key idea is to interpolate frequency-independent matrix elements of the single-particle irreducible part of the Liouvillian, obtained from a continued-fraction expansion of the cDMFT self-energy. We demonstrate that the ingredients of such an expansion possess a more local Fourier expansion than the functions involved in traditional interpolation schemes, and that Liouvillian interpolation inherently conserves causality. We illustrate our method for the one-dimensional Hubbard model using CDMFT, and for the two-dimensional Hubbard model using four-patch DCA. For the latter, we find that L-interpolation can (depending on doping) yield Fermi and Luttinger arcs which together form a closed surface.

Li$_2$MnO$_3$ remains a crucial component of the Li-excess layered cathode family, $(1-x)\,\mathrm{LiMO_2} + x\,\mathrm{Li_2MnO_3}$ ($M$ = Mn, Ni, Co, \dots), but its role in limiting Li-ion mobility remains under debate. Here we combine DFT+$U$, finite-temperature DMFT with a continuous-time quantum Monte Carlo impurity solver, and nudged-elastic-band (NEB) calculations to investigate Li$^+$ migration for a single Li vacancy in paramagnetic Li$_2$MnO$_3$. Dynamical electronic correlations within DMFT substantially reduce the activation energies of the lowest-barrier pathways, yielding $E_a = 0.18$ eV for the shortest-range Li jump and $E_a = 0.50$ eV for the next-lowest pathway. The 0.18 eV barrier quantitatively reproduces the short-range activation energy extracted from $μ^+$SR measurements, whereas the 0.50 eV barrier is consistent with the long-range transport activation energy obtained from ac-impedance measurements. This single-vacancy, paramagnetic DMFT description therefore provides a coherent explanation of both local and macroscopic probes without requiring highly clustered vacancy configurations or strong extrinsic disorder, an assumption compatible with nearly stoichiometric Li$_2$MnO$_3$ powders. Our results highlight the importance of finite-temperature dynamical correlations for Li-ion migration in correlated oxides and provide a framework for incorporating strong Coulomb interactions in future studies of transition-metal oxide battery materials.

Non-Hermitian systems have attracted significant theoretical interest due to their extreme properties. However, realizations have mostly been limited to classical applications or artificial setups. In this study, we focus on the quantum nature inherent in bosonic Bogoliubov-de Gennes (BdG) systems, which from the perspective of spectral theory corresponds to non-Hermiticity. Based on this insight, we propose a steady-state skin effect in quantum condensed matter utilizing such BdG non-Hermiticity. Specifically, we introduce BdG quantum terms arising from parametric pumping to the edge states of an underlying bosonic Hermitian Chern insulator, thereby realizing non-Hermiticity without dissipation. This system design has the advantage of being largely independent of microscopic model details. Through analysis using non-equilibrium Green's functions, we find that under open boundary conditions, a steady state exhibiting the non-Hermitian skin effect is realized. The pronounced corner particle accumulation observed in this steady state shows quadrature anisotropy, which manifests the bosonic quantum nature. Our results bridge the gap between the fascinating mathematics of non-Hermitian matrices and practical quantum physical systems.

Phenomena concerning altermagnets have opened up a window for unconventional analysis of the momentum space spin polarization (MSSP) of antiferromagnetic materials. Taking the example of one of the widely investigated altermagnets, CrSb, we explore the underlying mechanisms leading to the formation or breaking of altermagnetism. With the aid of DFT calculation and symmetry analysis, we study the behavior of MSSP in the altermagnetic bands of pristine CrSb, along with a few model structures designed from the pristine one by hypothetical vacancy engineering and interstitial doping. We show that the six-fold rotational symmetry of the pristine CrSb can be reduced to a two-fold rotational symmetry via vacancy and doping engineering. We discover the formation of fragmented nodal curves (FNCs) across the Brillouin zone when in an altermagnetic material when the symmetry is restricted to two-fold rotation. Unlike the typical nodal planes and axes, the location of the FNCs in the momentum space is found to be band-specific. The formation of FNCs is further validated by introducing uniaxial strain to CrSb and by examining the band structure of RbMnPO$_4$, as they both exhibit a two-fold rotational symmetry responsible for altermagnetism. We observe that, unlike the pristine case, these FNCs have the potential to manifest anomalous Hall conductivities (AHC), while the Néel vector orients along both in-plane and out-of-plane directions. This flexibility of the AHC will pave the way for the application of altermagnets in the futuristic quantum devices.

Exciton-phonon interactions govern the optical response of semiconductors, yet disentangling multiple coupling channels in lead halide perovskites remains challenging. We investigate CsPbBr3 microcrystals using photoluminescence, Raman and reflectance spectroscopy at low temperature, revealing the simultaneous presence of high-energy and Rashba excitons, each accompanied by distinct phonon replica series. High-energy exciton replicas are uniquely spaced by approximately 9 meV, whereas Rashba exciton replicas exhibit a characteristic approximately 6 meV spacing, indicating the specificity of the exciton-phonon coupling. Unsupervised machine learning applied to a large low-temperature photoluminescence dataset reveals these replica features are prevalent. With increasing temperature, replica features broaden and merge, evolving into a dominant longitudinal optical phonon coupling regime at room temperature. This work establishes direct spectroscopic evidence for concurrent, exciton-specific phonon coupling within a single material, offering new pathways to engineer light-matter interactions for optoelectronic and phonon-photon-based quantum device applications.

Dynamical mean-field theory (DMFT) provides an optimal local approximation for correlated lattice systems by mapping the lattice onto a self-consistent effective impurity model. To account for the missing long-range correlations, we propose a novel extended approach, which we term fluctuating dynamical mean-field theory (fDMFT). It incorporates collective fluctuations of auxiliary impurity models across different sites via functional integration. Technically, this method involves obtaining a family of DMFT solutions on a grid for a self-consistent auxiliary classical field applied to the lattice. While the result can, in principle, be improved diagrammatically, we find that the minimal version of the theory already yields accurate results, with lowest-order diagrammatic corrections offering only minor improvements. This consistent framework, based on our fluctuating local field concept, demonstrates superior performance for the nearly half-filled Hubbard model compared to other known diagrammatic extensions of DMFT.

We investigate the emergence of altermagnetism -- a collinear magnetic phase characterized by large non-relativistic spin splitting and zero net magnetization -- driven by electronic correlations on 3x3 geometrically frustrated square lattice. Using exact diagonalization of the simple and the extended Hubbard model, we analyze the interplay between on-site repulsion U , nearest-neighbor (NN) Coulomb interaction V and geometric frustration across various filling factors (N=8,9,10). Unlike traditional models that rely on single-particle anisotropy or specific sublattice geometries, our results demonstrate that altermagnetic signatures arise from many-body fluctuations in doped Mott insulators. We find that while geometric frustration suppresses altermagnetism at half-filling (N=9), carrier doping (hole (N=8) or electron (N=10)) stabilizes robust, rotationally symmetric altermagnetic correlations. Furthermore, we identify a critical threshold for the NN interaction V at intermediate coupling (U=4) where the altermagnetic state in the electron-doped sector undergoes a first-order-like transition, whereas strong coupling (U=10) stabilizes well-formed local moments and preserves the anisotropic spin texture. This work establishes a fluctuation-mediated route to altermagnetism on symmetric geometrically frustrated lattices and identify carrier concentration and the NN Coulomb interaction as a critical tunable parameter for controlling magnetic anisotropy in strongly correlated systems.

Angle-resolved photoemission spectroscopy (ARPES) reveals ubiquitous quasiparticle ``kinks'' near $\sim$70 meV and $\sim$40 meV across cuprate superconductors, often accompanied by peak--dip--hump (PDH) structures. These features point to strong coupling between electrons and low-energy bosonic excitations, but the microscopic origin has remained elusive due to the limitations of conventional density-functional theory (DFT) and the high cost of beyond-DFT methods. Here, we systematically study the electron--phonon coupling (EPC) in hole-doped infinite-layer CaCuO$_2$ using the Strongly Constrained and Appropriately Normed (SCAN) density functional, explicitly including magnetic effects. We find a substantial EPC strength $λ$ of $\sim$0.5 in the magnetic phase, producing kinks and PDH structures in the 40-80~meV window in excellent agreement with experiments. The dominant contribution arises from rotational oxygen phonons, while breathing modes contribute little. Our results establish strong EPC in cuprates, highlight the key role of rotational phonons, and provide a framework for understanding spectral anomalies in cuprates and beyond.

The discovery of superconductivity in square-planar nickelates has offered a rich materials platform to explore the origins of cuprate-like superconductivity. Experimental investigations however have largely been limited to the infinite-layer $R$NiO$_2$ ($R$=rare-earth) nickelates. Here, we construct a phase diagram of multi-layer square-planar Nd$_{n+1}$Ni$_n$O$_{2n+2}$ compounds and discover signatures of superconductivity for $n$ = 4 - 8. Upon decreasing the dimensionality $n$, the superconducting anisotropy evolves due to 4$f$ electron effects, and electronic structure characteristics approach cuprate-like behavior. Magnetic fluctuations persist from within the superconducting regime and into the over-doped, non-superconducting regime. Remarkably, the superconducting regime overlaps with that of chemically-doped infinite-layer nickelates, demonstrating underlying commonalities and distinct differences across varying structural realizations of square-planar nickelates. Our work establishes this layered template for creating new nickel-based superconductors.

Previous studies of the $S=1/2$ triangular-lattice $J_1$--$J_2$ Heisenberg antiferromagnet have inferred the existence of a non-magnetic ground-state phase for an intermediate range of $J_2$, but disagree concerning whether it is a gapped $\mathbb{Z}_2$ quantum spin liquid (QSL), a gapless (Dirac) QSL, or a weakly symmetry-broken phase. Using an improved dynamical density-matrix renormalization group method, we investigate the relevant intermediate $J_2$ regime for cylinders with circumferences from 6 to 9. Depending on the initial state and boundary conditions, we find two {\it distinct} variational states. The higher energy state is consistent with a Dirac QSL. In the lower-energy state, both the static and dynamical properties are qualitatively similar to the magnetically ordered state at $J_2=0$, suggestive of either a weakly magnetically ordered non-QSL or a gapped QSL proximate to a continuous transition to such an ordered state.

A comprehensive spectroscopic map of the electronic, magnetic, and lattice excitations is presented for the bilayer nickelate La$_3$Ni$_2$O$_7$ using Raman scattering at ambient pressure. Upon entering the spin density wave state below 153 K, the $A_{1g}$ channel exhibits an abrupt electronic spectral gap with a clear isosbestic point. In contrast, the $B_{1g}$ and $B_{2g}$ channels are dominated by pronounced two-magnon (2M) excitations, representing an unambiguous signature of incipient Mottness. These 2M signals in both channels constitute direct evidence for two distinct in-plane spin exchange interactions along the Ni-O bonding and its diagonal directions. Calculations based on the spin wave theory further reveal that the 2M mode in the $B_{2g}$ channel arises from the competition between two bond-diagonal antiferromagnetic interactions mediated by nickel $d_{x^2-y^2}$ orbitals. Furthermore, emergent low-energy 2M excitations below 10 meV are found to originate from distinct, weaker spin moments, strongly supporting spin disproportionation. Simultaneously, an anomalous softening of $B_{1g}$ phonons from 280 down to 4.5 K is uncovered, suggesting the presence of an incipient lattice instability leading to checkerboard-type breathing modulations. Collectively, these findings identify a ground state of the bilayer nickelate characterized by competing bond-diagonal interactions, spin disproportionation, and an incipient lattice instability, establishing key ingredients for understanding the mechanism of nickelate superconductivity.

Identifying whether a many-body ground state is gapped or gapless is a fundamental yet challenging problem, especially in quantum spin liquids. In this work, we develop a gap-diagnostic scheme based on the polarization amplitude defined via a twist operator, evaluated within the infinite density-matrix renormalization group (iDMRG) framework. As a benchmark, analysis of the spin-$1/2$ XXZ chain demonstrates that the polarization amplitude clearly distinguishes the gapless Tomonaga-Luttinger liquid from the gapped Néel phase. We then extend this framework to infinite cylinders of the spin-$1/2$ XY-$J_χ$ model on the square lattice. We find that the polarization amplitude sharply detects the transition between the gapless XY phase and the gapped chiral spin liquid phase. These results show that polarization amplitudes provide a strong energy-gap diagnostic in two-dimensional frustrated quantum magnets, including quantum spin liquids.

The spin-$1/2$ $J_1$-$J_2$ Heisenberg model on the square lattice exhibits a debated intermediate phase between Néel antiferromagnetic and stripe ordered regimes, with competing theories proposing plaquette valence bond, nematic, and quantum spin liquid ground states. We apply the Prometheus variational autoencoder framework -- previously validated on classical (2D, 3D Ising) and quantum (disordered transverse field Ising) phase transitions -- to systematically explore the $J_1$-$J_2$ phase diagram using a multi-scale approach. For $L=4$, we employ exact diagonalization with full wavefunction analysis via quantum-aware VAE. For larger systems ($L=6, 8$), we introduce a reduced density matrix (RDM) based methodology using DMRG ground states, enabling scaling beyond the exponential barrier of full Hilbert space representation. Through dense parameter scans of $J_2/J_1 \in [0, 1]$ and comprehensive latent space analysis, we identify the structure factor $S(π,π)$ and $S(π,0)$ as the dominant order parameters discovered by the VAE, with correlations exceeding $|r| > 0.97$. The RDM-VAE approach successfully captures the Néel-to-stripe crossover near $J_2/J_1 \approx 0.5$--$0.6$, demonstrating that local quantum correlations encoded in reduced density matrices contain sufficient information for unsupervised phase discovery. This work establishes a scalable pathway for applying machine learning to frustrated quantum systems where full wavefunction access is computationally prohibitive.

The spin-$1/2$ $J_1$-$J_2$ Heisenberg model on the square lattice exhibits a debated intermediate phase between Néel antiferromagnetic and stripe ordered regimes, with competing theories proposing plaquette valence bond, nematic, and quantum spin liquid ground states. We apply the Prometheus variational autoencoder framework -- previously validated on classical (2D, 3D Ising) and quantum (disordered transverse field Ising) phase transitions -- to systematically explore the $J_1$-$J_2$ phase diagram via unsupervised analysis of exact diagonalization ground states for a $4 \times 4$ lattice. Through dense parameter scans of $J_2/J_1 \in [0.3, 0.7]$ with step size 0.01 and comprehensive latent space analysis, we investigate the nature of the intermediate regime using unsupervised order parameter discovery and critical point detection via multiple independent methods. This work demonstrates the application of rigorously validated machine learning methods to open questions in frustrated quantum magnetism, where traditional order parameter identification is challenged by competing interactions and limited accessible system sizes.

Qualitative changes in thermodynamic and single-particle properties characterize the transition between the pseudogapped electronic liquid and the Fermi liquid. Recent cold-atom experiments on a simulator of the Hubbard model with nearest-neighbor hoppings \cite{kendrick2025pseudogap} showed that the isothermal compressibility $κ(δ)$ has a maximum as a function of doping $δ$. Here we use the two-particle self-consistent plus (TPSC+) approach to explain these experiments and connect the maximum in $κ(δ)$ to the transformation of the single-particle spectrum from the pseudogapped to the metallic regime. This elucidates the nature of the pseudogap (PG). Specifically, the maximum in $κ(δ)$ practically coincides with the doping at which the precursor of the lower $(π,π)$ spin density wave (SDW) band at the antinodal point crosses the zero-frequency $ω=0$. The Knight shift, $χ_{sp}(0,0)(δ)$, as a function of doping, should also have a maximum. The maxima in both quantities should exist, at sufficiently low temperatures ($T$), in both the intermediate $U \approx U_{Mott}$ and weak $U < U_{Mott}$ interaction limits. In both limits, the mechanism is critical thermal SDW fluctuations. At the antinodal pseudogap, the correlation length at $δ_{max}(T)$ can be small, controlled not by static but by dynamic critical thermal fluctuations. We also find that the SDW fluctuations are incommensurate at $δ=δ_{max}$. We predict that, at low $T$, the multiple peaks in the spin susceptibility in the incommensurate case lead to more than two SDW precursor peaks in the spectral function and density of states. By allowing access to parameter regimes relevant to cuprates-including further-neighbor hopping ($t', t''$) and low temperatures, our work provides a high-impact tool for further studies by the broader community.

Quantum phase transitions driven by electronic correlations are central to understanding the physics of graphene and related two-dimensional materials. A paradigmatic example is the semimetal-to-Mott-insulator transition on the honeycomb lattice, governed by the Gross-Neveu-Heisenberg universality class, yet consensus on its critical exponents has remained elusive for over a decade due to severe finite-size effects and the absence of rigorous conformal bootstrap benchmarks. Here we try to resolve this long-standing controversy by performing projector determinant quantum Monte Carlo simulations on lattices of unprecedented size, reaching 10,368 sites. By developing a novel projected submatrix update algorithm, we achieve a significant algorithmic speedup that enables us to access the thermodynamic limit with high precision. We observe that the fermion anomalous dimension and the correlation length exponent converge rapidly, while the boson anomalous dimension exhibits a systematic size dependence that we resolve via linear extrapolation. To validate our analysis, we perform parallel large-scale simulations of the spinless $t$-$V$ model on the honeycomb lattice, which belongs to the Gross-Neveu-Ising class. Our results for the $t$-$V$ model show agreement with conformal bootstrap predictions, thereby corroborating the robustness of our methodology. Our work provides state-of-the-art critical exponents for the honeycomb Hubbard model and establishes a systematic finite-size scaling workflow applicable to a broad class of strongly correlated quantum systems, paving the way for resolving other challenging fermionic quantum critical phenomena.

The interplay between Dirac fermions and electronic correlations on the honeycomb lattice hosts a fundamental quantum phase transition from a semimetal to a Mott insulator, governed by the Gross-Neveu-Heisenberg (GNH) universality class. Despite its importance, consensus on the precise critical exponents remains elusive due to severe finite-size effects in numerical simulations and the lack of conformal bootstrap benchmarks. Here we try to resolve this long-standing controversy by performing projector determinant quantum Monte Carlo (QMC) simulations on lattices of unprecedented size, reaching 10,368 sites. By developing a novel projected submatrix update algorithm, we achieve a significant algorithmic speedup that enables us to access the thermodynamic limit with high precision. We observe that the fermion anomalous dimension and the correlation length exponent converge rapidly, while the boson anomalous dimension exhibits a systematic size dependence that we resolve via linear extrapolation. To validate our analysis, we perform parallel large-scale simulations of the spinless $t$-$V$ model on the honeycomb lattice, which belongs to the Gross-Neveu-Ising class. Our results for the $t$-$V$ model, including the first QMC determination of the fermion anomalous dimension, show agreement with conformal bootstrap predictions, thereby corroborating the robustness of our methodology. Our work provides state-of-the-art critical exponents for the honeycomb Hubbard model and establishes a systematic finite-size scaling workflow applicable to a broad class of strongly correlated quantum systems, paving the way for resolving other challenging fermionic quantum critical phenomena.

Topology without quasiparticles has emerged as a key framework for understanding Mott insulators, where Green's-function zeros encode nontrivial topological structure. Yet, experimental detection of these zeros represents a challenge. Using exact diagonalization of the one-dimensional Hubbard model with an impurity and Zeeman field, supported by exact analytic results, we show that Green's-function zeros manifest as an in-gap spectral weight in the unitary scattering regime. In this limit, we map the impurity problem onto a doped Mott insulator and identify the resulting in-gap state as a "zeron" excitation which is a localized doublon (holon) for an attractive (repulsive) potential. The zeron spectral weight and its associated zero vanish above a critical Zeeman field. Our results imply that Green's function zeros have in fact already been observed in experiments, and establish impurity and magnetic-field tuning as practical tools for controlling their topology.

The spatially-resolved laser-based high resolution angle resolved photoemission spectroscopy (ARPES) measurements have been performed on the optimally-doped HgBa$_2$Ca$_2$Cu$_3$O$_{8+δ}$ (Hg1223) superconductor with a $T_c$ at 133 K. Two distinct regions are identified on the cleaved surface: the single Fermi surface region where only one Fermi surface is observed, and the double Fermi surface region where two Fermi surface sheets are resolved coming from both the inner (IP) and outer (OP) CuO$_2$ planes. The electronic structure and superconducting gap are measured on both of these two regions. In both cases, the observed electronic states are mainly concentrated near the nodal region. The momentum dependence of superconducting gap deviates from the standard d-wave form. These results indicate that the surface electronic structure of Hg1223 behaves more like that of underdoped cuprates.

Altermagnets provide a natural platform for studying and exploiting piezomagnetism. In this paper, we introduce a class of insulating altermagnets in two dimensions (2D) referred to as Dirac quadrupole altermagnets, and show based on microscopic minimal models that the orbital piezomagnetic polarizability of such altermagnets has a topological contribution described by topological response theory. The essential low-energy electronic structure of Dirac quadrupole altermagnets can be understood from a gapless parent phase (i.e., the Dirac quadrupole semimetal), which has important implications for their response to external fields. Focusing on the strain-induced response, here we demonstrate that the topological piezomagnetic effect is a consequence of the way in which strain affects the Dirac points forming a quadrupole. We consider two microscopic models: a spinless two-band model describing a band inversion of $s$ and $d$ states, and a Lieb lattice model with collinear Néel order. The latter is a prototypical minimal model for altermagnetism in 2D and is realized in a number of recently proposed material compounds, which are discussed.

The understanding of the mechanisms responsible for superconductivity in strongly correlated systems is an interesting and important subject in condensed matter physics. Several theoretical proposals were considered for these systems. The Coulomb interaction between electrons allow a new approach to study this problem. In this paper, we use a usual Hubbard model with a local repulsive interaction to describe a 2D system. The system of equations are solved using the Green's functions method, within a Hubbard-I mean field approximation, which allows to treat the strong interaction limit. We consider both cases of attractive and repulsive interactions and obtain the zero temperature phase diagram of the model. Our results show, in the repulsive case, the existence of a superconducting ground state mediated by the kinetic electronic energy and described by a non-local order parameter. A minimum value of the repulsive interaction $U_{min}$ is required to create a pairing state. At finite temperatures, for strong interactions, the critical temperature $T_c$ shows a saturation similar to the Bose-Einstein condensation observed for strong attractive interactions.

The undoped Kondo necklace in 1D is a paradigmatic and well understood model of a Kondo insulator. This work performs the first large-scale study of the 1D Anderson-lattice underlying the Kondo necklace with quasi-exact numerical methods, comparing this with the perturbative effective 1D Kondo-necklace model derived from the former. This study is based on an exact mapping of the Anderson model to one of a superconducting pairing layer connected to a metallic reservoir which is valid in arbitrary spatial dimensions, thereby linking the previously disparate areas of reservoir-enhanced superconductivity, following Kivelson's pioneering proposals, and that of periodic Kondo-systems. Our work reveals that below the length-scales on which the insulating state sets in, which can be very large, superconducting and density-density correlations are degenerate and may both appear to approach an almost ordered state, to a degree that far exceeds that of any isolated 1D pairing layer with short-range interactions. We trace these effects to the effective extended-range coupling that the metallic layer mediates within the pairing layer. These results translate directly to the appearance of near-long-range magnetic order at intermediate scales in the Kondo-systems, and explain the strong renormalization of the RKKY-coupling that we effectively observe, in terms of the back-action of the pairing layer onto the metallic layer. The effects we predict could be tested either by local probes of quasi-1D heavy fermion compounds such as CeCo$_2$Ga$_8$, in engineered chains of ad-atoms or in ultracold atomic gases.

Strontium titanate ($\mathrm{SrTiO_3, STO}$) dome-shaped superconducting transition temperature as a function of chemical potential, consistent with STO experiments, and shows that tunable s-wave and d-wave symmetries are modulated by doping. Superconducting fluctuations validate the mean-field approximation at low temperatures but destroy pairing at higher temperatures. The charge-density-wave order competes with superconductivity, enhances the effective electron mass inversely with the chemical potential, and increases with the interaction strength $U$ and the temperature $T$. SDW order is rare and fragile, while an additional magnetic term induces subtle band splitting. These findings suggest e-e contributions to STO's transport anomalies and provide criteria to distinguish e-e from e-ph origins, offering insights for engineering higher $T_c$ in dilute systems.

We present the mixed-configuration approximation (MCA) based on the auxiliary master equation approach impurity solver to study multiorbital correlated systems under equilibrium and nonequilibrium conditions within dynamical mean-field theory (DMFT). We benchmark the method for bulk and layered SrVO$_3$ in equilibrium and apply it to a prototypical nonequilibrium geometry in which a voltage bias is applied perpendicular to the layer via reservoirs held at different chemical potentials. For bulk SrVO$_3$, MCA reproduces the metallic state at moderate interaction strengths, but it overestimates the weight of the lower band relative to quantum Monte Carlo (QMC) and fork tensor product state (FTPS) solvers. With respect to QMC and FTPS, MCA yields an earlier metal-to-insulator transition as the electron-electron interaction is increased. In layered SrVO$_3$ at equilibrium, MCA partially captures the orbital polarization in favor of the in-plane $xy$ orbital, although not as strong as in the DMFT-converged results obtained with QMC. However, when performing a one-shot impurity calculation initialized with the DFMT-QMC results, MCA yields orbital occupations which show a stronger charge polarization in favor of orbital $xy$. This suggests that our approach can be used to study multiorbital impurity problems when the focus is to assess properties without performing the full DMFT self-consistent loop. Finally, under applied bias, we observe a pronounced redistribution of orbital occupations, demonstrating that the method captures bias-driven orbital charge transfer in realistic materials in nonequilibrium conditions.

Lattice gauge theories (LGTs) provide a powerful framework for studying confinement, topological order, and exotic quantum matter. In particular, the paradigmatic phenomenon of confinement, where dynamical matter is coupled to gauge fields and forms bound states, remains an open problem. In addition, LGTs can provide low-energy descriptions of quantum spin liquids, which is the focus of ongoing experimental research. However, the study of LGTs is often limited theoretically by their numerical complexity and experimentally in implementing challenging multi-body interactions, such as the plaquette terms crucial for the realization of many exotic phases of matter. Here we investigate a $(2+1)$D $\mathbb{Z}_2$ LGT coupled to hard-core bosonic matter featuring a global U(1) symmetry, and show that dynamical matter naturally induces sizable plaquette interactions even in the absence of explicit plaquette terms in the Hamiltonian. Using a combination of density matrix renormalization group simulations and neural quantum state calculations up to a system size of $20 \times 20$, we analyze the model across different fillings and electric field strengths. At small coupling strength, we find a large plaquette expectation value, independent of system size, for a wide range of fillings, which decreases in the presence of stronger electric fields. Furthermore, we observe signatures of a confinement-deconfinement transition at weak coupling strengths. Our results demonstrate that dynamical U(1) matter can induce complex multi-body interactions, suggesting a natural route to the realization of strong plaquette terms and paving the way for realizing a topological quantum spin liquid protected by a large gap.

We explore the effects of non-magnetic impurities (vacancy disorder) on the quantum paramagnetic phases stabilized by SU($N$) designer Hamiltonians on bipartite lattices. Using the results of our quantum Monte Carlo simulations, we demonstrate that isolated vacancies seed emergent spin $S=1/2$ moments in their vicinity when the low-temperature state has valence bond solid order. Indeed, our quantum Monte Carlo results for the low-temperature susceptibility in such regimes shows clear evidence of the vacancy-induced Curie tails associated with these emergent moments, and our zero-temperature projector Monte Carlo results on the ground-state wavefunction in the valence bond basis provide additional evidence in support of this picture. Further, for such designer Hamiltonians on the Lieb lattice with two additional sites on each bond of a square lattice, we identify a low-temperature spin liquid-like regime with no sign of spin or valence bond order. This liquid-like regime serves as a test bed for validating a recently-developed argument concerning the effects of vacancy disorder in such low temperature regimes. Consistent with this argument, we find that isolated vacancies do not seed emergent local moments in such spin liquids. Instead, in the presence of vacancy disorder, emergent local moments are associated with the presence of monomers in maximum-density dimer packings of the corresponding diluted lattice.

Characterizing the ground-state properties of disordered systems, such as spin glasses and combinatorial optimization problems, is fundamental to science and engineering. However, computing exact ground states and counting their degeneracies are generally NP-hard and #P-hard problems, respectively, posing a formidable challenge for exact algorithms. Recently, Tensor Networks methods, which utilize high-dimensional linear algebra and achieve massive hardware parallelization, have emerged as a rapidly developing paradigm for efficiently solving these tasks. Despite their success, these methods are fundamentally constrained by the exponential growth of space complexity, which severely limits their scalability. To address this bottleneck, we introduce the Branch-and-Bound Tensor Network (BBTN) method, which seamlessly integrates the adaptive search framework of branch-and-bound with the efficient contraction of tropical tensor networks, significantly extending the reach of exact algorithms. We show that BBTN significantly surpasses existing state-of-the-art solvers, setting new benchmarks for exact computation. It pushes the boundaries of tractability to previously unreachable scales, enabling exact ground-state counting for $\pm J$ spin glasses up to $64 \times 64$ and solving Maximum Independent Set problems on King's subgraphs up to $100 \times 100$. For hard instances, BBTN dramatically reduces the computational cost of standard Tropical Tensor Networks, compressing years of runtime into minutes. Furthermore, it outperforms leading integer-programming solvers by over 30$\times$, establishing a versatile and scalable framework for solving hard problems in statistical physics and combinatorial optimization.

Spectroscopic measurements in model one-dimensional (1D) correlated systems offer insights for understanding their two-dimensional counterparts, which include the cuprate and iron pnictide/chalcogenide superconductors. A major challenge is the identification of such correlated systems with dominantly 1D physics. In this work, inelastic neutron scattering measurements on NaFe$_{0.53}$Cu$_{0.47}$As single crystal directly reveal quasi-1D spin excitations, resulting from atomic order that lead to magnetic Fe and nonmagnetic Cu chains. The dominant exchange interaction is antiferromagnetic along the chain ($SJ_{\rm \parallel}\approx90.1(3)$~meV), whereas the inter-chain couplings are much weaker ($SJ_{\rm \perp}\approx-2.4(1)$~meV and $SJ_{\rm c}\approx0.15(5)$~meV). The quasi-1D spin excitations in NaFe$_{0.53}$Cu$_{0.47}$As stem from both the Néel and stripe vectors, with Néel excitations sensitive to Fe impurities on the Cu site. The spin excitations in quasi-1D NaFe$_{0.53}$Cu$_{0.47}$As and quasi-2D FeSe exhibit a striking resemblance, suggesting a common origin for their coexistent stripe and Néel excitations. Our findings demonstrate magnetic dilution in NaFeAs leads to dimension reduction of its magnetic degree of freedom, presenting a strategy for discovering low-dimensional quantum materials.

We experimentally investigate a frustrated spin-exchange antiferromagnet in a quantum simulator, composed of N = 114 dipolar Rydberg atoms arranged into a kagome array. Motivated by a recent theoretical proposal of a gapless U(1) Dirac spin liquid ground state, we use local addressing to adiabatically prepare low-energy states. We measure the local polarization and spin-spin correlations over this adiabatic protocol, and observe our system move from a staggered product state, through an intermediate magnetic crystal, and finally into a disordered, correlated liquid. We estimate the entropy density of this atomic liquid to be similar to that of frustrated magnetic insulators at liquid nitrogen temperatures. We compare the correlations in our liquid to those of a simple, parameter-free ansatz for the Dirac spin liquid, and find good agreement in the sign structure and spatial decay. Finally, we probe the static susceptibility of our system to a local field perturbation and to a geometrical distortion. Our results establish Rydberg atom arrays as a promising platform for the preparation and microscopic characterization of quantum spin liquid candidates.

Altermagnets have recently emerged as a new class of magnetic materials that combine compensated magnetic order with spin-split electronic band structures. In this work, we employ molecular beam epitaxy to grow MnTe thin films with controlled thickness on InP(111)A substrates. By performing angle-resolved photoemission spectroscopy measurements, we observe a large spin splitting of ~230 meV for bulk bands well below the Fermi level and identify surface states that cross the Fermi level. Electrical transport measurements reveal that a robust anomalous Hall (AH) effect persists down to 2 K and an AH sign reversal occurs near 175 K. By systematically tuning film thickness, growth conditions, and interfacial structure, we demonstrate that the AH response in MnTe films originates from the Berry curvature of surface states rather than from bulk bands. Our first-principles calculations reveal that this surface-state-driven AH effect is imprinted by the bulk altermagnetic order and remains unchanged for terminations with opposite Mn magnetic orientations. Our results establish a unique surface transport probe of bulk altermagnetism, demonstrate interface engineering as an effective route to generate and control the AH effect in altermagnets, and provide a unified understanding of the AH response in altermagentic MnTe films.

For the $t$-$J$ model, variational wave functions can generally be constructed based on an accurate description of antiferromagnetism (AFM) at half-filling and an exact phase-string sign structure under doping. The single-hole-doped and two-hole-doped states, as determined by variational Monte Carlo (VMC) simulations, display sharply contrasting behaviors. The single-hole state constitutes a ``cat state'' that resonates strongly between a quasiparticle component and a local loop-current component, with approximately equal weights. In the ground state, the quasiparticle spectral weight $Z_{\mathbf{k}}$ peaks at momenta $\mathbf{k}_0 \equiv (\pm\fracπ{2},\pm\fracπ{2})$. The total-energy dispersion versus $\mathbf{k}$ agrees remarkably well with the Green function Monte Carlo results. However, Landau's one-to-one correspondence hypothesis for quasiparticles breaks down here with the incoherent component exhibiting intrinsic magnetization originating from a minimal $2\times2$ loop current that forms a $4\times4$ pattern on the square lattice--a finding in excellent agreement with density matrix renormalization group (DMRG) calculations. In the two-hole ground state, a new pairing mechanism is revealed: the two holes are automatically fused into a tightly bound object consisting of an incoherent $d_{xy}$ pairing along the diagonal direction by compensating the local loop currents. This hole pair is again a ``cat state'' that resonates strongly between the incoherent $d_{xy}$ and a coherent $d_{x^2-y^2}$ Cooper channel to gain substantial hopping energy. Its size extends over an area of about $4\times 4$ lattice spacings, much smaller than the divergent AFM correlation length, implying that it should survive as a minimal superconducting building block even in the dilute doping regime. Experimental implications and the generalization to the finite-doping case are briefly addressed.

We study Rokhsar-Kivelson (RK) dimer and spin ice models realizing $U(1)$-lattice gauge theories in a wide class of quasi-one-dimensional settings, which define a setup for the study of few quantum strings (closed electric field lines) interacting with themselves and each other. We discover a large collection of mappings of these models onto three quantum chains: the spin-1/2 XXZ chain, a spin-1 chain, and a kinetically constrained fermion chain whose configurations are best described in terms of tilings of a rectangular strip. We show that the twist of boundary conditions in the chains maps onto the transverse momentum of the electric field string, and their Drude weight to the inverse of the string mass per unit length. We numerically determine the phase diagrams for these spin chains, employing DMRG simulations and find global similarities but also many interesting new features in comparison to the full 2D problems. For example, the spin-1 chain we obtain features a continuous family of degenerate ground states at its RK point analogous to a Bloch sphere, but without an underlying microscopic global $SU(2)$ symmetry. We also argue for the existence of a (stable) Landau-forbidden gapless critical point away from the RK point in one of the models we study using bosonization and numerics. This is surprising given that the full 2D problem is generically gapped away from the RK point. The same model also displays extensively many local conserved quantities which fragment the Hilbert space, arising as a consequence of destructive resonances between the electric field lines. Our findings highlight spin-chain mappings as a potent technique for the exploration of unusual dynamics, exotic criticality, and low-energy physics in lattice gauge theories.

We study when a density-matching Kohn--Sham (KS) description can reproduce a many-body Berry phase in a correlated insulator, despite the fact that geometric phases are functionals of the wave function. Focusing on the one-dimensional SSH--Hubbard chain on a ring as a controlled interacting topological model, we introduce a $U(1)$ twist $θ$ (flux insertion). The many-body ground state along the full twist cycle is computed by the density-matrix renormalization group (DMRG), while the onsite interaction $U$ is tuned from the noninteracting to the strong-coupling regime. At half filling in the inversion-symmetric gapped regime, our DMRG calculations show that the density remains constant within numerical accuracy over the entire $(θ,U)$ range studied. Thus, the density has no dependence on either the flux $θ$ or the interaction strength $U$. Accordingly, the symmetry-preserving density constraint collapses the KS reference to an SSH-type quadratic representative with $U$-independent geometric diagnostics. Nevertheless, the many-body wave function exhibits a nontrivial geometric response: the quantum metric associated with the $θ$-parametrized ground-state manifold depends on $θ$ at intermediate $U$ and is strongly suppressed at large $U$, consistent with the charge fluctuation freezing. Intriguingly, the KS and many-body Berry phases coincide throughout the gapped regime as $U$ is tuned from weak to strong coupling. We show that this agreement is best understood as symmetry-enforced $\mathbb{Z}_2$ sector matching, rather than as evidence that the density encodes the many-body Berry connection.

We study when a density-matching Kohn--Sham (KS) description can reproduce a many-body Berry phase in a correlated insulator, despite the fact that geometric phases are functionals of the wave function. Focusing on the one-dimensional SSH--Hubbard chain on a ring as a controlled interacting topological model, we introduce a $U(1)$ twist $θ$ (flux insertion). The many-body ground state along the full twist cycle is computed by the density-matrix renormalization group (DMRG), while the onsite interaction $U$ is tuned from the noninteracting to the strong-coupling regime. At half filling in the inversion-symmetric gapped regime, our DMRG calculations show that the density remains constant within numerical accuracy over the entire $(θ,U)$ range studied. Thus, the density has no dependence on either the flux $θ$ or the interaction strength $U$. Accordingly, the symmetry-preserving density constraint collapses the KS reference to an SSH-type quadratic representative with $U$-independent geometric diagnostics. Nevertheless, the many-body wave function exhibits a nontrivial geometric response: the quantum metric associated with the $θ$-parametrized ground-state manifold depends on $θ$ at intermediate $U$ and is strongly suppressed at large $U$, consistent with the charge fluctuation freezing. Intriguingly, the KS and many-body Berry phases coincide throughout the gapped regime as $U$ is tuned from weak to strong coupling. We show that this agreement is best understood as symmetry-enforced $\mathbb{Z}_2$ sector matching, rather than as evidence that the density encodes the many-body Berry connection.

An important problem in many-body physics is to reconstruct the spectral density from the imaginary-time domain Green's function. Typically, the imaginary-time Green's function is generated by Monte Carlo methods. As the one-point fermionic kernel diverges exponentially for large frequencies, numerical noise generically causes instabilities. We use a convolutional neural network to obtain the spectral density for a given imaginary time Green's function. The network is trained by data which we generate using random Gaussians. We improve the training data set available by including collision centers for the Gaussians rather than employing uniformly distributed Gaussians. Our network is constructed in such a way that its output fulfills positive semidefiniteness. We compare the results of our network with results of the Maximum Entropy method (MaxEnt), a standard method for the same reconstruction problem for the spectral density. This comparison is performed for three different cases, namely our Gaussian based test data as well as two physical models, the 1d Hubbard model showing spin-charge separation, and the two-dimensional SSH model in the self-consistent Born approximation. We find that the network outperforms MaxEnt when presented data close to the training set. For the physical models considered, MaxEnt recognizes physical features more precisely as compared to our network prediction. While it is hard to improve MaxEnt, the quality of the network depends on the training data set which can be systematically enhanced and improved.

This text reviews, hopefully in a pedagogical manner, a series of work on the automatic calculations of Feynman diagrams in the context of quantum nanoelectronics (Keldysh formalism) with an application to the Kondo effect in the out-of-equilibrium single impurity Anderson model. It includes a discussion of (A) how to deal with the proliferation of diagrams, (B) how to calculate them using the Tensor Cross Interpolation algorithm instead of Monte-Carlo and (C) how to resum the obtained series. These notes correspond to a lecture given at the Autumn School on Correlated Electrons 2025 in Jullich, Germany. The book with all the lectures of the school (edited by Eva Pavarini, Erik Koch, Alexander Lichtenstein, and Dieter Vollhardt) is available in open access.

Spin-orbit coupling (SOC) is incorporated into the phaseless plane-wave-based auxiliary-field quantum Monte Carlo (pw-AFQMC) method. This integration is implemented using optimized multiple-projector norm-conserving pseudopotentials, which are derived from the fully-relativistic (FR) atomic all-electron Dirac-like equation. The inclusion of SOC enables accurate phaseless pw-AFQMC calculations that capture both electronic correlation and SOC effects concurrently, greatly improving the method's applicability for studying systems containing heavy atoms. We discuss the form of FR pseudopotentials and detail the corresponding formulations of phaseless pw-AFQMC with a two-component Hamiltonian in the spinor basis. The accuracy of our approach is demonstrated by computing the dissociation energy of molecule I2 and the cohesive energy of bulk Pb, highlighting the large influence of SOC in both. Subsequently, we determine the transition pressure of the III-V compound InP from its zinc-blende to rock-salt phase by constructing and analyzing their respective equations of state.

We demonstrate the emergence of altermagnetism at the surfaces of antiferromagnets, vastly expanding the number of material candidates with altermagnetic characteristics and establishing a route to two-dimensional altermagnetism through surface-induced symmetry breaking. We do so by developing a surface spin group formalism that fully classifies all surface magnetic states and identifies altermagnetic surface spin groups that can arise at the surfaces of antiferromagnets. We use this formalism to identify over 140 antiferromagnetic entries from the MAGNDATA database with at least one altermagnetic surface, often times with multiple such surfaces in the same material. We illustrate this emergent phenomenon in a realistic Lieb lattice-based minimal model and present ab initio calculations on two representative material candidates, NaMnP and FeGe$_2$, exhibiting $d$-wave and $g$-wave surface altermagnetism, respectively. Our theory naturally resolves the contradiction of recent experimental reports of $d$-wave ARPES measurements on metallic Lieb lattice compounds that have been shown to be antiferromagnetic in the bulk. Hence, we establish a new paradigm for generating two-dimensional altermagnetism by functionalizing the abundant material class of collinear antiferromagnets as viable platforms for controlled surface altermagnetism, creating natural materials for future hybrid device implementation.

Ab initio auxiliary-field quantum Monte Carlo (AFQMC) is a systematically improvable many-body method, but its application to extended solids has been severely limited by unfavorable computational scaling and memory requirements that obstruct direct access to the thermodynamic and complete-basis-set limits. By combining tensor hypercontraction with $\mathbf{k}$-point symmetry, we reduce the computational and memory scaling of ab initio AFQMC for solids to $\mathcal O(N^3)$ and $\mathcal O(N^2)$, respectively, with an arbitrary basis, comparable to diffusion Monte Carlo. This enables direct and simultaneous thermodynamic-limit and complete-basis-set AFQMC calculations across insulating, metallic, and strongly correlated solids, without embedding, local approximations, empirical finite-size corrections, or composite schemes. Our results establish AFQMC as a general-purpose, systematically improvable alternative to diffusion Monte Carlo and coupled-cluster methods for predictive ab initio simulations of solids, enabling accurate energies and magnetic observables within a unified framework.

High-temperature spin transport in integrable quantum spin chains exhibits a rich dynamical phase diagram, including ballistic, superdiffusive, and diffusive regimes. While integrability is known to survive in static and periodically driven systems, its fate in the complete absence of time-translation symmetry, particularly in interacting random quantum circuits, has remained unclear. Here we construct integrable random quantum circuits built from inhomogeneous XXZ R-matrices. Remarkably, integrability is preserved for arbitrary sequences of gate layers, ranging from quasiperiodic to fully random, thereby explicitly breaking both continuous and discrete time-translation symmetry. Using large-scale time-dependent density-matrix renormalization group simulations at infinite temperature and half filling, we map out the resulting spin-transport phase diagram and identify ballistic, superdiffusive, and diffusive regimes controlled by the spectral parameters of the R-matrices. The spatiotemporal structure of spin correlations within each regime depends sensitively on the inhomogeneity, exhibiting spatial asymmetry and sharp peak structures tied to near-degenerate quasiparticle velocities. To account for these findings, we develop a generalized hydrodynamics framework adapted to time-dependent integrable circuits, yielding Euler-scale predictions for correlation functions, Drude weights, and diffusion bounds. This approach identifies the quasiparticles governing transport and quantitatively captures both the scaling exponents and fine structures of the correlation profiles observed numerically. Our results demonstrate that exact Yang-Baxter integrability is compatible with stochastic quantum dynamics and establish generalized hydrodynamics as a predictive framework for transport in time-dependent integrable systems.

We analyze the magnetization curve of the spin-1/2 kagome Heisenberg model in a magnetic field. Using state-of-the-art variational wavefunctions based on neural networks, we confirm the presence of robust magnetization plateaus at $m=1/3$, $5/9$ and $7/9$ of the saturation value, stabilized by a spontaneous symmetry breaking of lattice translations with a $\sqrt{3}\times \sqrt{3}$ unit cell. Regarding the more challenging $m=1/9$ plateau, we find two competing valence bond crystals depending on the system size, both breaking translation as well as point group symmetries and with a larger $3\times 3$ unit cell. Such quantum states with local modulations of the magnetization average values could be observed experimentally in the near future.

High entropy oxides (HEO) hold the potential to revolutionize the conventional material paradigms by leveraging high order of chemical disorder that induces highly desirable exotic phases for advanced applications. Here, we devise a methodology to enhance the efficiency of an artificial photonic synapse using a high entropy rare earth nickelate. Combined with epitaxial strain, we show that high entropy can further manipulate the phase of these locally disordered materials. Using time-averaged and time resolved Terahertz (THz) spectroscopy as dynamic probe, for the first time we show a rare combination of i) crystal axis dependent insulator to metal THz electronic phase transition and ii) coexistence of negative and positive THz photoconductivity at room temperature. Detailed analysis within theoretical models, including density functional theory (DFT)-based band structure calculations, suggest origin of these properties as disproportionate ordering of oxygen vacancies. Based on these findings, a conceptual THz-based artificial photonic synapse is proposed. This work underlines the pivotal role of HEO in advancing diverse THz functionalities, representing a critical step toward futuristic applications like THz-based high-speed computing and communication with an emphasis in THz frequency domain.

Open quantum systems host mixed-state phases that go beyond the symmetry-protected topological and spontaneous symmetry-breaking paradigms established for closed, pure-state systems. Developing a unified and physically transparent classification of such phases remains a central challenge. In this work, we introduce a purification-based framework that systematically characterizes all mixed-state phases in one-dimensional systems with $\mathbb{Z}_2^σ \times \mathbb{Z}_2^τ$ symmetry. By introducing an ancillary $κ$ chain and employing decorated domain-wall constructions, we derive eight purified fixed-point Hamiltonians labeled by topological indices $(μ_{στ},μ_{τκ},μ_{κσ}) \in \{\pm1\}^3$. Tracing out the ancilla recovers the full structure of mixed-state phases, including symmetric, strong-to-weak spontaneous symmetry breaking, average symmetry-protected topological phases, and their nontrivial combinations. Interpolations between the eight fixed points naturally define a three-dimensional phase diagram with a cube geometry. The edges correspond to elementary transitions associated with single topological indices, while the faces host intermediate phases arising from competing domain-wall decorations. Along the edges, we identify a class of critical behavior that connects distinct strong-to-weak symmetry-breaking patterns associated with distinct strong subgroups, highlighting a mechanism unique to mixed-state settings. Large-scale tensor-network simulations reveal a rich phase structure, including pyramid-shaped symmetry-breaking regions and a fully symmetry-broken phase at the cube center. Overall, our purification approach provides a geometrically transparent and physically complete classification of mixed-state phases, unified with a single $\mathbb{Z}_2^σ \times \mathbb{Z}_2^τ \times \mathbb{Z}_2^κ$ model.

Hamiltonian formulations of lattice field theories provide access to real-time dynamics, but their simulation is difficult to implement efficiently. Trotter-Suzuki decompositions are at the center of time evolution computation, either on quantum hardware or classically, for instance with the use of tensor networks. While low-order Trotterizations remain the standard choice due to their simplicity, higher-order schemes offer the potential for improved efficiency. In this work we outline a short guide to Trotter-Suzuki schemes and their implementations in general. To help with this, we highlight new efficient schemes found by our optimization framework, and demonstrate their performance on the Heisenberg model.

Computational discovery of magnetic materials remains challenging because magnetism arises from the competition between kinetic energy and Coulomb interaction that is often beyond the reach of standard electronic-structure methods. Here we tackle this challenge by directly solving the many-electron Schrödinger equation with neural-network variational Monte Carlo, which provides a highly expressive variational wavefunction for strongly correlated systems. Applying this technique to transition metal dichalcogenide moiré semicondutors, we predict itinerant ferromagnetism in WSe$_2$/WS$_2$ and an antiferromagnetic insulator in twisted $Γ$-valley homobilayer, using the same neural network without any physics input beyond the microscopic Hamiltonian. Crucially, both types of magnetic states are obtained from a single calculation within the $S_z=0$ sector, removing the need to compute and compare multiple $S_z$ sectors. This significantly reduces computational cost and paves the way for faster and more reliable magnetic material design.

We present TensorCircuit-NG, a next-generation quantum software platform designed to bridge the gap between quantum physics, artificial intelligence, and high-performance computing. Moving beyond the scope of traditional circuit simulators, TensorCircuit-NG establishes a unified, tensor-native programming paradigm where quantum circuits, tensor networks, and neural networks fuse into a single, end-to-end differentiable computational graph. Built upon industry-standard machine learning backends (JAX, TensorFlow, PyTorch), the framework introduces comprehensive capabilities for approximate circuit simulation, analog dynamics, fermion Gaussian states, qudit systems, and scalable noise modeling. To tackle the exponential complexity of deep quantum circuits, TensorCircuit-NG implements advanced distributed computing strategies, including automated data parallelism and model-parallel tensor network slicing. We validate these capabilities on GPU clusters, demonstrating a near-linear speedup in distributed variational quantum algorithms. TensorCircuit-NG enables flagship applications, including end-to-end QML for CIFAR-100 computer vision, efficient pipelines from quantum states to neural networks via classical shadows, and differentiable optimization of tensor network states for many-body physics.

Accurate many-body treatments of condensed-phase systems are challenging because correlated solvers such as full configuration interaction (FCI) and the density matrix renormalization group (DMRG) scale exponentially with system size. Downfolding and embedding approaches mitigate this cost but typically require prior selection of a correlated subspace, which can be difficult to determine in heterogeneous or extended systems. Here, we introduce a stochastic cluster expansion framework for efficiently recovering the total correlation energy of large systems with near-DMRG accuracy, without the need to select an active space a priori. By combining correlation contributions from randomly sampled environment orbitals with an exactly treated subspace of interest, the method reproduces total energies for non-reacting and reactive systems while drastically reducing computational cost. The approach also provides a quantitative diagnostic for molecule-solvent correlation, guiding principled embedding decisions. This framework enables systematically improvable many-body calculations in extended systems, opening the door to high-accuracy studies of chemical processes in condensed phase environments.

We show that the concept of coherent phonon generation by second order response to incident electric laser fields, which is a hallmark of pump-probe spectroscopy on conventional solids, can be expanded to include frustrated quantum magnets. For that purpose, we analyze the Raman force on the shear phonons of a frustrated magnetoelectric bilayer spin system. The bilayer is a stacked triangular magnet, motivated by recently emerging type-II van der Waals multiferroic transition metal dihalides and comprises a spin system which allows for incommensurate spiral order. The magnon excitations are treated by linear spin wave theory. In the spiral state, a finite electric polarization is obtained from the spin-current interaction which induces a coupling of the magnons to the electric field. Scattering of the bilayer shear phonons from the magnons is derived from a magnetoelastic energy. In this scenario, a mixed three-point response function for the Raman force is evaluated. We find it to be strongly anisotropic and very sensitive to the magnon lifetime.

Determining atomic structure from spectroscopic data is central to materials science but remains restricted to a limited set of techniques and material classes, largely due to the computational cost and complexity of structural refinement. Here we introduce ActiveStructOpt, a general framework that integrates graph neural network surrogate models with active learning to efficiently determine candidate structures that reproduce target spectra with minimal computational expenditure. Benchmarking with X-ray pair distribution function data, and with the more computationally demanding simulations of X-ray absorption near-edge spectra (XANES) and extended X-ray absorption fine structure (EXAFS), demonstrate that ActiveStructOpt reliably determines structures that match closely in spectra across diverse materials classes. Under equivalent computational budgets, ActiveStructOpt outperforms existing structure determination methods. By enabling data-efficient, multi-objective structural refinement across a broad range of computable spectroscopic techniques, ActiveStructOpt provides a flexible and extensible approach to atomic structure determination in complex materials.

First-principles investigations of electron-phonon interactions (EPIs) play a crucial role in understanding a wide range of phenomena in physics and materials science. Among various approaches, the finite difference method offers a direct route to capture higher-order EPIs and is compatible with diverse electronic structure solvers. However, its considerable computational cost limits its broader application. To overcome this bottleneck, we present a machine learning electron-phonon interaction (MLEPI) pipeline that predicts force constants and electronic Hamiltonians for modeling EPIs from finite difference calculations, improving efficiency by orders of magnitude without compromising accuracy. The performance of MLEPI is validated by studying the temperature dependence of the electronic band properties in bilayer graphene, where both first- and second order EPIs are treated on an equal footing. Using a heterogeneous edge network, the pipeline integrates both interlayer and intralayer interactions, making it particularly suitable for studying multilayer materials. With its inherent adaptability and ease of transfer to other applications, our methodology provides a robust tool with a very favorable accuracy/efficiency balance for investigating EPIs in large-scale material systems.

Coherent phonon excitation offers a nonthermal route to control quantum phases of condensed matter. In this work, we employ first-principles calculations to investigate the phonon landscape of (TaSe4)2I in its charge-density-wave (CDW) phase. We identify nine symmetry-preserving Raman-active modes that can suppress the Gamma-Z direct gap to the meV scale and render the system globally gapless by generating Weyl nodes at generic k points. Among them, the 2.51 THz CDW amplitude mode A(18) directly weakens the Ta-chain tetramerization, approaching a transient restoration of the uniform-chain geometry. It is also the most efficient mode owing to its low frequency and a relatively small critical displacement dominated by Ta motions. Other Raman modes, dominated by Se vibrations, require significantly larger displacements to reach the Weyl-semimetallic regime and are generally less effective than A(18) at reducing the Ta-chain tetramerization. Furthermore, we explore nonlinear phonon-phonon interactions and find that the low-frequency infrared-active mode B3(7) (1.14 THz) exhibits strong anharmonic coupling with A(18), providing an indirect pathway to drive the system toward a Weyl-semimetallic regime. Our results provide predictive insight for ultrafast pump-probe experiments and present a generalizable framework for lattice-driven topological switching in quasi-one-dimensional quantum materials.