2026-08-14 | | Total: 18
Combined with numerical analytic continuation techniques, quantum Monte Carlo (QMC) methods enable the extraction of real-frequency dynamical properties from imaginary-time correlation functions. However, the efficient computation of imaginary-time correlation functions by QMC simulations can (depending on the particular model used) be challenging, particularly for operators that are off-diagonal in the computational basis. In this work, we present an efficient and general algorithm within the stochastic series expansion (SSE) framework for evaluating imaginary-time correlation functions of both diagonal and off-diagonal operators. The algorithm builds on a discrete imaginary-time slicing of the SSE operator string, which provides correlation functions on a grid of well-defined imaginary-time points with no discretization error. For off-diagonal operators, we derive estimators that integrate directly into the existing SSE directed-loop or cluster updating schemes, introducing only minimal computational overhead. We benchmark the method on the one-dimensional transverse-field Ising model (sampling with cluster updates) and XXZ spin chain (using directed-loop sampling), demonstrating excellent agreement (with only statistical errors) with exact diagonalization of small systems. We also study larger systems to demonstrate efficiency.
A main theme in modern condensed matter physics is the emergence of fractionalized excitations and gauge structures from quantum spin systems. However, understanding how these exotic degrees of freedom reconfine into new phases of matter remains a fundamental challenge. In this work, we demonstrate that Polyakov's confinement mechanism-a foundation of compact gauge theory-can serve as a dynamical engine to transform a Dirac spin liquid into a symmetry-protected topological (SPT) state. Starting with a Dirac spin liquid with spin-dependent gauge fluxes, we show that the single-occupancy constraint inherent to the physical Hilbert space triggers monopole condensation, dynamically confining bulk spinons while preserving gapless edge modes-realizing a spinon analog of the quantum spin Hall effect. Using a large-scale variational Monte Carlo simulation on a triangular antiferromagnet with an additional Dzyaloshinskii-Moriya interaction, we provide microscopic evidence for this confined SPT phase, including a characteristic area law for the Wilson loop and vanishing topological entanglement entropy. Furthermore, we identify a measurable spin-pump response under magnetic fields, which directly encodes the Berry curvature of the parent Dirac cones. Our results reveal a previously unexplored pathway to SPT physics, bridging the fields of gauge theory, quantum magnetism, and topological matter.
Including higher order diagrammatic corrections to the strong-coupling expansion is mainly limited by the evaluation of high-dimensional, time-ordered integrals. In this work we present and compare four different parametrizations of the integrands in order to obtain a low-rank (quantics) tensor-train representation using tensor cross interpolation. Particular emphasis is placed on a quantics time-difference formulation in which the required retarded convolutions are performed directly in quantics tensor-train form. Using controlled Gaussian benchmarks, we analyze the accuracy, bond dimensions, and computational scaling of the different approaches. We then validate the most promising formulations in self-consistent equilibrium and nonequilibrium DMFT calculations and demonstrate calculations up to the third order in the strong-coupling expansion. Finally, we extend the solver to impurity models with retarded density-density interactions and apply it within nonequilibrium extended DMFT. Our results show that tensor cross interpolation substantially reduces the cost of evaluating higher-order diagrams and provides a controlled, systematically improvable framework for nonequilibrium quantum impurity calculations.
We report the magnetic properties of a series of lithium-rich antiperovskites (Li$_2M$)$Ch$O ($M$ = Fe, Co, Mn and $Ch$ = Se, S) where transition metal and lithium ions are randomly distributed on the X-sites of the X$_3$BA structure, thereby forming a strongly diluted magnetic sublattice. Our study hence enables us to investigate the evolution of magnetic order at fixed 1/3-filling -- which is in the vicinity but slightly above the percolation threshold -- upon variation of the spin size, the magnetic anisotropy, and the orbital configuration. The data imply the absence of a distinct Curie-Weiss behavior up to 350~K but show rather large and weakly temperature-dependent magnetic susceptibility. We observe clear signatures of long-range antiferromagnetic order evolving in the 1/3-filled and strongly diluted magnetic X-site lattice with increasing Néel temperatures from $T_{\rm{N}}\simeq 30$~K in (Li$_2$Mn)$Ch$O to $\simeq 50$~K in (Li$_2$Fe)$Ch$O and $70-90$~K in (Li$_2$Co)$Ch$O. Except for $M$ = Co, the chalcogenide has no sizable effect on $T_{\rm N}$. We conclude significant magnetic coupling and short-range magnetic correlations at well above $T_{\rm N}$ which is in line with the observation of a broad electron spin resonance signal at room temperature. The actual ordering temperatures are strongly diminished by magnetic dilution. While structural parameters such as the tolerance factor and bonding angles do not strongly affect $T_{\rm N}$, a key parameter is the magnetic anisotropy of the transition metals.
The gapped Kondo effect describes the screening of an $S=\frac12$ impurity spin locally coupled via an antiferromagnetic exchange interaction to a conduction-electron system exhibiting a finite hard gap. Using a combination of a Lanczos transformation and a self-consistent configuration-interaction scheme, we numerically investigate the local phase diagram. Furthermore, we show that the different phases can be characterized by several topological invariants: the conventional momentum-space Chern number of the underlying two-dimensional host system, corresponding to a Chern insulator; the B-space Chern number, defined by coupling the impurity spin to a fictitious local magnetic field $\boldsymbol B$, in the limit $B \to 0$; and the S-space Chern number, defined for a classical impurity spin, i.e., a vector of fixed length. The classical-spin limit is obtained for $B \to \infty$. By varying the field strength, we can therefore continuously interpolate between quantum-impurity-spin and classical-impurity-spin Hamiltonians and investigate whether the corresponding phase diagrams are likewise continuously connected. The gapped underscreened Kondo effect is studied for an impurity spin $S>\frac12$ as well as in the classical-spin limit approached via $B \to \infty$ or $S \to \infty$. Different variants of scattering theory are employed to interpret the resulting phases. Finally, the gapped two-channel overscreened Kondo effect, realized by coupling a quantum impurity spin equally to the local electron spins of both orbitals within a unit cell, is shown to be characterized by spontaneous particle-hole symmetry breaking. This leads to a highly nontrivial quantum-classical phase diagram.
Correlation functions $\langle O_1(x) O_2(0) \rangle$ reveal scaling dimensions through spatial decay. We instead consider static susceptibility, the change in $\langle O_1(x) \rangle$ from perturbing the Hamiltonian by $O_2(0)$, which we term causation for short. In a conformal field theory (CFT), dimensional analysis predicts decay of $|x|^{-2Δ}$ for correlation and $|x|^{-2Δ+1}$ for causation. Yet we find causation can decay up to fifteen additional orders in $x$ through a general mechanism, which we trace to time-derivative fields being unable to contribute to static response. In higher-dimensional CFTs, this mechanism ensures leading causation arises from primaries, even when descendants dominate correlation, which we leverage with DMRG to identify a previously unresolved corner primary of $Δ\approx 8.8$ and a heavy magnetic line defect primary of $Δ\approx 4.6$ in the $(2+1)$D critical Ising model. Moreover, the same mechanism governs edge-mode localization in $(1+1)$D gapless symmetry-protected topological phases, explaining previously observed anomalously small edge-mode splittings and guiding our construction of spin chains with splittings as small as $1/L^{18}$ and $1/L^{25}$.
We report mean-field critical behavior in Eu$_2$MnSi$_2$O$_7$, a melilite-type ferrimagnet with spin-only Eu$^{2+}$ and Mn$^{2+}$ moments and negligible orbital contributions. Magnetization measurements combined with neutron powder diffraction reveal critical exponents close to the mean-field values, indicating that long-range dipolar interactions govern the asymptotic critical behavior in this insulating ferrimagnet. The refined magnetic structure, described by the magnetic space group $P2_12_1^\prime2^\prime$, exhibits a tilted ferrimagnetic configuration driven by the Dzyaloshinskii-Moriya interaction, reflecting the noncentrosymmetric nature of the lattice. These results extend the applicability of mean-field theory to complex insulating magnets and establish Eu$_2$MnSi$_2$O$_7$ as a platform for exploring ferrimagnetism and long-range interactions. To our knowledge, this is the first insulating ferrimagnet in which dipolar interactions drive mean-field criticality.
Magnetic defect engineering offers a route to manipulate and control magnetic and electronic states in quantum materials. In the Mn(Bi,Sb)2Te4 family of topological magnetic insulators, Sb substitution facilitates site mixing between Mn and Sb atoms that affects magnetic order and band topology. Here we directly probe these defect-induced magnetic interactions using inelastic neutron scattering on single crystals of MnSb2Te4. We find that antisite mixing generates inequivalent and disordered magnetic sublattices where strong defect-induced antiferromagnetic coupling produces an optical magnon with a large gap. Semi-classical spin-dynamics simulations accurately capture magnetic excitations in the ordered and paramagnetic states and identify that linear Mn-Te-Mn bonds mediate coupling to antisite magnetic defects.
Despite its fundamental importance in the theory of strongly correlated electrons, the nature of the ground state of the two-dimensional doped Hubbard model remains intensely debated. Variational approaches provide a powerful route to this problem, but their conclusions can depend sensitively on the chosen wave-function parameterization, the mean-field initialization, or the pinning fields used to guide the optimization, as well as on boundary conditions. This can favor one type of symmetry breaking over another, making it difficult to distinguish the genuine interplay of intertwined or competing orders from biases induced by the variational parameterization. Here, we introduce the Symmetry-Preserving Backflow Pairing (SBP) ansatz, a neural-network wave function that respects translational symmetry by construction and thereby avoids these broken-symmetry minima. The SBP ansatz reaches state-of-the-art variational energies for the $t$-$t'$ Hubbard model on lattices up to $24\times24$ with $504$ electrons, below those of competing pure stripe solutions. By extrapolating to the thermodynamic limit, we find robust evidence for $d$-wave superconducting order, resolving a long-standing question about the $1/8$-doped model at $t'/t=-0.2$ and $U/t=8.0$. Built on general principles of symmetry and locality, the SBP wave function provides a broadly applicable variational representation for challenging interacting fermionic systems.
A comprehensive analysis of pairing and bound states of magnons in the Heisenberg $S\!=\!1/2$ square-lattice $J_1\!-\!J_2$ ferro-antiferromagnetic model is presented. We highlight the similarities and differences between the bound states of spin flips on a lattice and those of particles in the continuum. Magnon bound states at finite pair momentum are studied throughout the Brillouin zone and a convenient lattice partial-wave nomenclature is advocated. The mechanisms of enhanced stability or fragility of these bound states for the high-symmetry pair momenta are identified and quantified as relating to the effective dimensional reduction or enhancement, respectively. These effects are also shown to control the evolution of the bound states with the model parameters, providing a transparent framework for understanding magnon pairing in a more general setting. A method to determine the bound state phase boundaries is presented.
We develop a lattice characterization of fermionic 't Hooft anomalies of finite internal symmetries in (1+1)D and (2+1)D, formulated in terms of obstructions to symmetric short-range-entangled (SRE) states. We consider lattice systems formed by tensor product of onsite fermionic and bosonic Hilbert spaces, and finite internal symmetry given by a central extension $\mathbb Z_2^F\to G_f\to G_b$. We extract a hierarchy of fermionic anomaly indices for a given symmetry operator. In (1+1)D, an exact lattice symmetry is characterized by a pair of cohomological data $(n_2,ν_3)$. For $G_f=G_b\times\mathbb Z_2^F$, we show that a symmetry with trivial anomaly indices $(n_2,ν_3)$ is onsiteable and hence admits a symmetric SRE state, establishing that these indices faithfully detect the lattice anomaly. Comparing with continuum QFT, we find that exact lattice symmetries do not realize the additional $H^1(BG_b,\mathbb Z_2)$ anomaly layer in continuum QFT. In particular, for $G_b=\mathbb Z_2$, exact lattice symmetries realize only the even $\mathbb Z_4$ subgroup of the continuum $\mathbb Z_8$ classification. In (2+1)D, we identify three successive anomaly layers of cohomological data $(n_2,n_3,ν_4)$. We show that a nontrivial value of any layer obstructs symmetric SRE. For $G_f=G_b\times\mathbb Z_2^F$, it also forbids a symmetric invertible state. When the bosonic group $G_b$ is non-trivially extended by fermion parity, the lattice obstruction to SRE states does not generally coincide with the continuum 't Hooft anomaly. We explicitly construct a $\mathbb Z_4^F$ lattice symmetry in (2+1)D with nontrivial lattice anomaly index that forbids symmetric SRE states, even though its continuum anomaly is trivial. Our results highlight a mismatch between lattice and continuum anomalies and motivate a systematic study of which continuum anomalies admit exact microscopic lattice realizations.
We study a spin-$\frac12$ impurity coupled to the boundary of a strongly correlated spin-$s$ Takhtajan--Babujian chain, an integrable model whose low-energy physics is described by a perturbed $SU(2)_{2s}$ Wess--Zumino--Witten conformal field theory. While boundary conformal field theory determines the low-energy universality class of the weak-coupling regime, exact Bethe Ansatz methods reveal a sequence of boundary quantum phase transitions in which impurity-bound states emerge and reorganize the Hilbert space into multiple excitation towers built on distinct boundary configurations. This tower restructuring provides the organizing principle for a rich boundary phase diagram extending beyond the conventional Kondo regime. Weak antiferromagnetic coupling realizes the overscreened $2s$-channel Kondo universality class, whereas stronger couplings generate localized boundary modes and qualitatively new screening mechanisms. To describe the resulting thermodynamics, we develop a generalized thermodynamic Bethe Ansatz framework that captures the multi-tower structure across all regimes. The impurity entropy reproduces the boundary conformal field theory prediction in the overscreened Kondo regime but develops pronounced nonmonotonic temperature dependence once boundary-bound states appear, in quantitative agreement with large-scale finite-temperature matrix-product-operator simulations. Complementary dynamical calculations reveal sharp threshold features in the impurity spectral function that directly track the underlying tower structure. Together, boundary conformal field theory, exact Bethe Ansatz, generalized thermodynamic Bethe Ansatz, and tensor-network simulations provide a unified description of impurity screening, boundary-bound-state formation, and excitation-tower reconstruction in a correlated spin-$s$ chain.
While the necessity of a moiré potential for fractional Chern insulators (FCIs) in rhombohedral graphene-hBN superlattices, first predicted in Ref. 1, is now grounded in experiments, a theory of its origin and importance remains at large. We present a mechanism---the moiré capacitor effect---that enhances the moiré potential by electrostatically imprinting the valence charge density onto the conduction bands. We derive the analytical form of this term and reveal its crucial role in stabilizing a parent state with Chern number $C=1$ at filling $ν=1$. We propose a parent state theory which posits that stability of the Chern insulator and flatness of its hole excitations are necessary for obtaining FCIs upon doping. We then perform multi-band exact diagonalization calculations to confirm the emergence of FCIs at $ν= 2/3$ in the presence of the moiré capacitor effect. Our FCI state is stabilized by inter-band fluctuations, unlike in the moire-free case which collapses with band-mixing. We provide the first consistent theory for this state in aligned samples and explain its absence in unaligned ones.
Unlike elementary particles, which obey either bosonic or fermionic exchange statistics, certain quasiparticles, known as anyons, are predicted to exhibit Abelian or non-Abelian braiding statistics. While braiding Abelian anyons modifies the wavefunction by a 'statistical phase', braiding non-Abelian anyons implements a unitary transformation of the state within a degenerate subspace of states. Experimental evidence of non-Abelian braiding has thus far remained elusive. Here, we report a 'time-domain braiding' signature of non-Abelian anyons in the $ν= 5/2$ fractional quantum Hall state, by extending our previously demonstrated approach with Abelian anyons at $ν= 1/3$. Our approach is based on measurements of the current fluctuations arising from weak partitioning of a highly dilute one-dimensional edge mode. We independently probe the partition noise of the downstream charged mode and also that of the upstream neutral mode. These independent measurements agree with our theoretical predictions for 'time-domain braiding' of the downstream Abelian and the upstream non-Abelian anyons, respectively, in the 'particle-hole Pfaffian' topological order. Together, these results provide evidence for the presence of non-Abelian anyons.
A chiral superconductor carries chiral Majorana modes along its edges, and a single integer, the Bogoliubov--de Gennes Chern number, counts them. Thirty years of candidate materials have not yielded a measurement of that integer, because the magnetic signatures usually invoked are not topologically protected. Rhombohedral graphene makes the question both urgent and answerable: magnetic imaging resolves rewritable time-reversal-breaking domains inside the superconducting phase, while quantum oscillations reveal a normal state too intricate to reconstruct pocket by pocket. We show that the low-temperature thermal Hall conductance returns the integer directly, with no such reconstruction. For band-projected pairing it equals the pairing-vortex winding enclosed by the occupied regions of momentum space. Splitting the intravalley Hamiltonian into symmetric and antisymmetric parts isolates the trigonal warping and finite Cooper pair momentum of the real material: the antisymmetric part is topologically inert, direct Chern calculations across $525$ parameter points show the invariant preserved, and one inequality marks where a Bogoliubov Fermi surface removes quantization. The plateau $κ_{xy}/T=(π^2k_B^2/6h)\,C_{\rm BdG}$ then reads out the integer, its sign reverses with the imaged domain, a written domain wall should carry $2|C_{\rm BdG}|$ Majorana channels, and the thermometry required already resolves single thermal quanta in encapsulated graphene at millikelvin temperatures.
Understanding and coherently controlling the properties of interacting quantum many-body systems is a central challenge in non-equilibrium quantum physics. While in the past decades, a wide range of many-body Hamiltonians have been introduced to study quantum chaos, atypical eigenstates, and quantum resources, systematically engineering and continuously tuning these properties within a single microscopic model remains largely unexplored. Here, we employ a pulse engineering scheme to construct an effective Hamiltonian that continuously interpolates between a chaotic Heisenberg (XXX) chain with a local impurity and the ZX Hamiltonian. Along this interpolation, we identify several families of analytically tractable atypical eigenstates embedded in the excited-state spectrum with distinct entanglement and nonstabilizerness properties. In the XXX limit, these states exhibit exact plateaus in both entanglement and stabilizer Rényi entropy and correspond to coherent superpositions of long-range valence-bond solid (VBS) states. As the pulse strength increases, the effective Hamiltonians exhibit a hierarchy of new set of approximate entanglement plateaus in the low-energy spectrum. Interestingly, in the fully pulse-engineered ZX limit, we uncover a distinct pair of long-range entangled stabilizer eigenstates, corresponding to rainbow scar states. We further show that pulse engineering preserves the distinct chaotic and non-chaotic regimes of the original model, while that is largely absent in the dynamical generation of entanglement and nonstabilizerness. The pulse-engineered models generate nearly identical quantum resources in both regimes, revealing a partial decoupling between quantum chaos and quantum-resource generation. Our results establish pulse engineering as a versatile framework for generating many-body Hamiltonians with structured eigenstates and tunable quantum resources.
Moiré engineering in twisted two-dimensional (2D) materials radically alters low-energy bands, interactions and topological quantum states. Despite extensive studies, optical spectroscopy of interacting moiré bands in the characteristic far-infrared (FIR) regime has remained largely unexplored due to extreme experimental challenges. Using a newly developed millikelvin FIR platform, we report the observation of the long-sought-after characteristic FIR resonances of flat-band electrons in magic-angle twisted bilayer graphene (MATBG). We observe highly tunable spectroscopic signatures of interacting light and heavy fermions that constitute the flat bands in MATBG. Using the topological heavy-fermion model (THF), we show that itinerant topological electrons act as an "antenna" that couples strongly to the optical field, with resonant frequencies renormalized by the hybridization with localized heavy electrons. We establish optical selection rules of MATBG which uncovers the key symmetry governing light-heavy fermion hybridization. At charge neutrality, we observe pronounced resonances at energies below the on-site Coulomb energy, implying the emergence of new many-body modes. Our experiments and modeling provide a fundamental understanding of light-matter interactions in MATBG and enable resonant optical spectroscopy of moiré bands down to millikelvin temperatures.
Disordered thin films undergoing a superconductor-insulator transition provide a controlled setting for studying pseudogap physics in the absence of competing electronic orders. Although theory predicts that local Cooper pairing can survive deep into the insulating phase, direct spectroscopic confirmation has remained experimentally inaccessible because tunneling measurements in highly insulating films require ultra-high-resistance junctions and picoampere current sensitivity. Here we combine ultra-high-resistance planar tunneling spectroscopy on amorphous indium oxide films with quantum Monte Carlo simulations of the attractive Hubbard model to probe the single-particle excitation spectrum deep in the insulating regime. We find striking agreement between experiment and theory: the single-particle gap not only survives across the superconductor-insulator transition, but increases substantially with disorder, reaching values more than twice those observed on the superconducting side. With increasing temperature, the gap fills rather than closes, while coherence peaks are suppressed and spectral weight redistributes to energies far exceeding the gap scale. Our results provide direct quantitative experimental confirmation of the theoretically predicted Cooper-pair insulating state with localized Cooper pairs and establish a unified connection between the pseudogap above $T_c$ and the insulating gap as manifestations of pairing without global phase coherence.