2026-04-17 | | Total: 22
Magnetic correlations of doped Mott insulators hold the key to the unusual characteristics of many quantum materials. Recent experiments with ultracold atoms in optical lattices have provided new information about the magnetic properties of the Fermi-Hubbard model on a square lattice. We demonstrate that recent measurements indicate that a single doping-dependent energy scale determines both static correlations and dynamical response of these systems. To understand these experimental findings, we employ a self-consistent formalism to describe the coupling between antiferromagnetic magnons and doped holes, and we uncover the emergence of a universal magnetic energy scale at finite doping, which we denote by $J^*$. We present the single- and two-magnon spectral properties at finite doping and discuss the appearance of a bimagnon peak in lattice-modulation spectroscopy, at frequencies set by $J^*$. Furthermore, we argue that this same energy scale sets the onset of pseudogap phenomena, leading to the hypothesis $k_BT^* = c J^*$, with $c$ an order one number. We identify another low-energy scale emerging from our analysis of magnetic excitations, and argue that it controls the stability of Néel order at the lowest temperatures, ultimately driving a transition to an incommensurate spin-density-wave at finite doping. We discuss the relation between this low-energy scale and the nature of fermionic quasiparticles. Our analysis suggests that stability of the commensurate antiferromagentic phase at finite doping can be controlled experimentally by introducing additional quasiparticle broadening via disorder or low-frequency noise.
The condensation of electron quartets, known as charge-4e superconductivity (SC), represents a novel quantum state of matter beyond the standard paradigm of Cooper pairing. However, concrete microscopic models realizing this phase in two dimensions remain a central challenge. Here, we introduce a non-engineered and sign-problem-free model, unambiguously demonstrating the emergence of a robust and high-temperature charge-4e SC phase using unbiased quantum Monte Carlo simulations. At zero temperature, the phase diagram reveals that charge-4e SC is the primary ground state in the strong-coupling regime. At finite temperature in the absence of charge-2e SC, we identify charge-4e SC through a Berezinskii-Kosterlitz-Thouless transition, marked by a universal jump in the superfluid stiffness consistent with a condensate of charge 4e. Remarkably, the transition temperature Tc increases nearly linearly with interaction strength, providing a robust mechanism for high-Tc quartet superconductivity. Furthermore, spectral analysis reveals a prominent pseudogap above Tc arising from strong phase fluctuations. Our results establish a canonical and numerically exact model system for charge-4e superconductivity, offering crucial guidance for its realization in experimental platforms such as moiré materials and ultracold atomic systems.
Plasmon modes, their dispersion, and the onset of damping when approaching the electron-hole continuum are well understood when electron correlations are weak. However, we know little about how this picture is modified and what additional features emerge in strongly correlated materials. Here, we present a fully ab initio approach to plasmon excitations that combines density functional theory with dynamical mean-field theory, and we use it to reconcile controversial electron energy-loss spectroscopy results in Sr$_2$RuO$_4$. In particular, we show that electronic correlations reproduce the plasmon dispersion, while generating a large intrinsic width already below the electron-hole continuum. An additional high-energy peak reflecting transitions between incoherent features and a sharp increase of the plasmon's energy-momentum dispersion, akin to waterfalls in photoemission spectroscopy, are identified as genuine correlation effects.
We investigate the ground-state magnetic structure of the spin-$1/2$ XXZ antiferromagnet on the triangular lattice in the easy-axis regime using the density-matrix renormalization group. By applying spiral boundary conditions, we exactly map finite $L\times L$ clusters onto one-dimensional chains while avoiding the spatial anisotropy inherent in cylindrical geometries. From symmetry-broken local magnetization profiles, we extract the three-sublattice moments and track their evolution with anisotropy. At the isotropic point, we obtain a positive sublattice moment of $0.21671$, consistent with previous numerical estimates. In the easy-axis regime, the ordered moments remain close to a zero-magnetization three-sublattice structure of the form $(2m,-m,-m)$ over a broad range of $Δ$. Extrapolation in $1/Δ$ shows that the positive sublattice moment stays well below the classical saturation value $1/2$, approaching $0.41873$ as $Δ\to\infty$, while the magnitude of the negative sublattice moment approaches $0.20832$. We further compare the energies of the Y state and the up-down-down state and find that the Y state is favored at zero field. Independent thermodynamic-limit energy calculations, performed without assuming any particular ordered pattern, yield an energy consistent with the Y-state solution. These results show that the easy-axis ground state does not simply cross over to a trivially saturated collinear Ising state, but instead remains a nontrivial three-sublattice ordered state selected from the macroscopically degenerate Ising manifold by quantum fluctuations.
Understanding how the interaction range and various types of disorder affect the level statistics of many-body quantum systems and lead to the emergence of many-body localization (MBL) is a challenging open frontier. We study the level statistics of a variant of the spin-$1/2$ Haldane-Shastry model with $1/r^α$ interactions, where $α{\geq}0$ parametrizes the range of the interactions, in the presence of position disorder and/or random magnetic fields. We find that neither position disorder nor random magnetic fields alone yields pristine Poisson statistics in this long-range interacting system; however, Poisson statistics emerge in their combined presence, suggesting the emergence of MBL when both types of disorder coexist. Interestingly, once random magnetic fields break the $SU(2)$ symmetry, the strength of the position disorder, $δ$, appears to play an important role, as evidenced by an approximate scaling collapse of the disorder-averaged gap ratios that is parametrized in terms of a single parameter, $αδ$.
The spontaneous breaking of symmetries is a cornerstone of physics, defining the phases of matter from the cosmological scale to the quantum realm. In condensed matter, electronic orders are classified by their behavior under fundamental symmetries like spatial inversion (parity). While even-parity orders, such as conventional superconductivity and charge density waves, are ubiquitous, their odd-parity counterparts--predicted to host exotic phenomena such as gapless quasiparticle excitations and novel collective modes--are comparatively elusive states of quantum matter. Here, using high-resolution scanning tunneling microscopy and angle-resolved photoemission spectroscopy on the kagome metal CsV$_3$Sb$_5$, we report the discovery of an inversion symmetry-breaking $f$-wave charge bond order. We show that this phase, which preserves translation symmetry, is stabilized by the spontaneous opening of a spectral gap at a previously overlooked Dirac point, providing a textbook condensed-matter realization of the Gross-Neveu model for dynamical mass generation and parity breaking. Intriguingly, this $f$-wave order is itself a intervening phase, vanishing abruptly below a temperature of 10\,K and pointing to a subsequent transition into a `hidden' electronic state that is invisible to local STM probes. Our findings establish odd-parity charge order as a novel phase of matter, here, embedded within the intricate hierarchy of correlated electronic orders on the kagome lattice.
The single orbital, one-dimensional, Hatano-Nelson Hamiltonian provides deep insight into the physics of non-Hermiticity, resulting from asymmetric left/right hopping, and its connections to localization. In the absence of disorder, its single particle eigenvalues $E_α$ lie on an ellipse in the complex plane whose extent in the imaginary direction is controlled by the degree of asymmetry. When randomness is introduced, two sets of real eigenvalues emerge at the extremes of the largest and smallest real part of $E_α$. These real eigenvalues are associated with localized eigenvectors. For spinless fermions, increasing near-neighbor interactions first cause a transition to a charge density wave phase, and ultimately, on finite lattices, a collapse of all eigenvalues to the real axis. In this paper, we explore the presence of real eigenvalues in the interacting, two-particle sector for the spinful case (Hubbard model) in a two-chain (two-band) geometry with a Hermitian interchain hopping. Our key results are to obtain the ``phase" diagrams for the existence of a purely real spectrum, as a function of the interaction strength, degree of non-Hermiticity, and interchain hopping. We study the sensitivity to boundary conditions of the spectral properties of our two-chain model with winding number analysis and explore the relationship between PBC doublon states and OBC skin modes. To address the question of stability in such non-equilibrium systems, we solve the dynamics at low filling according to Lindbladian evolution and find that the non-Hermitian description is able to qualitatively describe such systems.
Correlated transition metal oxides present exciting prospects as switches or memory and storage devices owing to the possibility to control electronic properties using various external stimuli. While their complex behaviour is known to stem from interplay between electronic correlations, atomic structure and orbital physics, they remain poorly understood on the microscopic level. Here, we investigate such origins as a function of temperature and pressure in the transition metal oxide Ti3O5. We find that the insulating room-temperature phase is characterized by one-dimensional zig-zag chains composed by two types of titanium dimers forming orbital selective valence bonds. At the thermal phase transition, one type of titanium dimer breaks up, resulting in an insulator to metal transition with a large orbital repopulation between the two states. Moreover, optical spectroscopy reveals that an additional pressure-driven insulator to metal transition occurs in Ti3O5 at room temperature. The phenomenology of this novel pressure-induced metallic transition is completely different from the insofar studied transitions and results from a competition between intra- and inter-dimer hopping. Our combined results suggest that Ti3O5 is a prototypical correlated transition metal oxide, where both correlations as well as orbital interactions need to be considered to fully understand the evolution of the electronic states.
Twisted bilayer graphene provides a paradigmatic platform for exploring the interplay between electronic topology and strong correlations. Within the topological heavy fermion model [Song and Bernevig, Phys. Rev. Lett. 129, 047601 (2022)], topology and electron interactions are brought together by including a weak hybridization between the bands of itinerant $c$- and localized $f$-electrons. Hybridization infuses concentrated Berry curvature into the $f$-band, while leaving it flat. These band features have motivated recent proposals of a Mott semimetal phase above the flavor-ordering temperature at charge neutrality. In this work, we develop an analytic theory of the quasiparticle dispersion and lifetime in the Mott semimetal. We reformulate the interacting flat-band Hamiltonian as an on-site Hubbard interaction defined on a set of non-orthogonal orbitals, and compute the electron Green's function using the equation-of-motion method, in close analogy with the Hubbard-III approximation. Unlike the conventional Hubbard model, in our case this approximation is controlled by a well-defined small parameter in the theory. We evaluate the electron self-energy and demonstrate the emergence of well-defined low-energy quasiparticles with the dispersion and relaxation rate proportional to the interaction strength. The quasiparticle spectrum is well-resolved in energy and in momentum down to the very vicinity of the Fermi level. Our results illustrate unconventional spectral properties arising from strong correlations and nontrivial quantum geometry, and have direct relevance for spectroscopic probes such as quantum twisting microscope experiments.
We develop a general framework to calculate the many-body density of states (DOS) of isolated and interacting quantum systems. Based on the generalized coherent state formalism and the Simon-Lieb bounds for a quantum partition function, our method provides a general method of calculation for the DOS in high-dimensional irreducible sectors. This framework further provides rigorous bounds for the ground state energy in each sector and enables the calculation of microcanonical observables across the entire spectrum. Using the Lipkin-Meshkov-Glick (LMG) model as a test bed, we validate our framework by successfully identifying quantum phase transitions (QPTs) and excited-state quantum phase transitions (ESQPTs) across spin sectors. Unlike existing model-specific numerical or analytical techniques, our formalism relies on general underlying symmetries, making it broadly applicable. Applying our method to the ferromagnetic transverse field Ising chain with power law interactions, we demonstrate that its highest-spin-sector DOS is qualitatively identical to that of LMG-type Hamiltonians. Our work establishes a versatile and computationally efficient bridge between algebraic structure and many-body thermodynamics.
Twisted bilayer graphene near integer fillings hosts correlated single-particle excitations whose dispersion and linewidth are increasingly accessible experimentally. We study these excitations using the topological heavy-fermion model, which captures both strong correlations and band topology of twisted bilayer graphene. In the decoupled limit, where both the single-particle fc hybridization and the Hund coupling between f and c electrons are absent, the model admits exact solutions in which free Dirac fermions coexist with interacting f electrons that form zero-width Hubbard bands. By treating the fc hybridization and Hund coupling perturbatively around this solvable limit, we obtain analytical results for the single-particle self-energy. From the resulting self-energy, we derive explicit expressions for both dispersion renormalization and scattering rates of both Hubbard-band excitations and low-energy Dirac modes, thereby establishing an analytical framework for understanding correlated excitations in twisted bilayer graphene. We analyze the scattering of the two kinds, Gamma3 and Gamma1,2, of Dirac electrons and find that they arise from different mechanisms. We also briefly investigate the effect of strain. Finally, we compare these analytical expressions with DMFT results for the same model.
In many strongly correlated electron systems, non-Fermi liquid behavior and unconventional superconductivity can be viewed as emerging from an effective 4-fermion interaction with a singular frequency dependence. A pairing instability in such a system is qualitatively different from that in a Fermi liquid and generally gives rise to multiple pairing states with topologically distinct gap functions. However, in the systems studied so far, a topologically trivial solution has the lowest energy. Here we show that a repulsive Hubbard-type interaction with a finite cutoff added to a model with a singular dynamical interaction selects, in some parameter range, the theretofore subleading, topologically nontrivial solution. We consider a minimal model that displays this behavior and show that the transformation between the topologically trivial and nontrivial gap functions necessarily occurs via an intermediate phase with topology-induced breaking of time-reversal symmetry.
The Migdal-Eliashberg theory (MET) describes electrons interacting with phonons in the adiabatic limit when the phonon Debye frequency is much smaller than the Fermi energy. A conventional belief is that MET holds even at strong coupling, when electron self-energy is large, and breaks down only near the point where the dressed phonon spectrum softens to near zero. We analyze numerically and analytically a different option -- collapse to a polaronic/bipolaronic ground state. The last scenario has never been analyzed in precise quantitative terms for a generic electron density. Using variational considerations, we establish rigorous upper bounds on the coupling $λ$, at which a FL state transforms into the bipolaron/polaron state. We show that at small and near-maximum densities, this happens well before a dressed phonon softens. This is true both in 2D and 3D systems; in the latter the upper bound on $λ$ tends to zero in the limit of small or near-full density. We present analytical reasoning for this behavior based on hints extracted from exact diagrammatic treatment of the on-site Holstein model for the spin polarized case and argue that polarons are produced by fermions with energies comparable to the bandwidth; i.e., polaron formation is outside the realm of MET. Closer to half-filling, the leading instability upon increasing $λ$ is towards a charge-density-wave state (CDW), and there exists a strong coupling regime of MET near this instability, while the polaron/bipolaron state develops at larger $λ$ out of a CDW-ordered state and inherits a CDW order over some range of coupling.
It is widely believed that in an adiabatic limit a Fermi liquid state of an electron-phonon system described by Migdal-Eliashberg theory remains stable before a dressed phonon softens. Using Holstein model as a prototypical example and variational/analytic considerations we demonstrate that in a wide range of fillings both in 3D and 2D, a polaronic/bi-polaronic state emerges before phonon softening; at small filling in 3D this happens already at weak coupling. We show that a polaronic/bi-polaronic state emerges, upon increasing coupling, via an intermediate pseudogap-type mixed state, in which some fermions regain Fermi liquid behavior, yet Luttinger theorem is broken. At even larger couplings the density of states gradually approaches its form in the atomic limit.
Unlike conventional charge-2e superconductors, a charge-4e superconductor exhibits long-range coherence of electron quartets rather than Cooper pairs. Clear zero-temperature realizations of charge-4e superconductivity remain rare. Here, we investigate the zero-temperature phase diagram of the attractive SU(4) Hubbard model with numerically exact, large-scale quantum Monte Carlo (QMC) simulations overcoming major technical hurdles. We identify both charge-2e and charge-4e superconducting phases. Upon increasing interaction, charge-2e correlations are suppressed and eventually vanish, while the charge-4e correlations remain robust and converge with system size, signaling the onset of a quartet-condensed phase. Interestingly, across the charge-2e--charge-4e transition, single electrons remain gapped, while charge-2e correlations exhibit a scaling behavior inconsistent with a conventional Landau description. These features are naturally captured by a fractionalized framework in which the physical charge-2e order parameter is a composite field coupled to an emergent non-Abelian gauge structure. We formulate an Sp(4) gauge-Higgs theory that realizes deconfined quantum pseudocriticality between the Higgs (charge-2e) phase and the confined (charge-4e) phase. The Sp(4) gauge-Higgs theory yields pseudocriticality through a fixed-point collision, and its one-loop collision-point exponents quantitatively track the QMC results. Our results establish charge-4e superconductivity as a bona fide zero-temperature phase, provide a simple model for future studies in a numerically exact framework, and reveal an unconventional route to superconducting criticality.
We show that the spin conductivity of a metal approaching a ferromagnetic quantum critical point exhibits divergent fluctuation corrections. This effect arises from critical spin fluctuations and constitutes a spin analog of the Aslamazov-Larkin theory of paraconductivity in superconductors. The spin current is derived in linear response within a Gaussian-level treatment of the effective action for a system with easy-plane magnetic anisotropy. We demonstrate the consistency of our spin transport theory by showing that it (i) fulfills the Ward identity and (ii) yields vanishing spin stiffness in the normal state. The critical enhancement of the spin conductivity is interpreted as incipient spin superfluidity in the quantum critical region. This is further supported by an intuitive picture based on the current-loop representation of the easy-plane ferromagnet.
We develop a controlled theoretical framework for the topological heavy fermion model relevant to magic-angle twisted bilayer graphene, where low density conduction electrons hybridize with a lattice of strongly interacting f-sites. By tracing out the localized electrons, we derive an effective action for the conduction electrons with long-range in time effective interactions, built from correlators of the single f-site problem. We identify a small hybridization-phase-space parameter resulting in a controlled loop expansion, enabling the derivation of nonperturbative results in either the interaction or the hybridization strength. To tree-level, the results are equivalent to the Hubbard I approximation. At higher loop order, we derive two key results applicable to temperatures above the flavor ordering temperature and below the on-site charging energy: 1) the quasi-particle lifetime, 2) the flavor susceptibility of the system. Remarkably, despite being strongly interacting, we find the susceptibility to accurately obey a Curie-Weiss law parametrically close to the Curie temperature.
Hydrodynamic electron transport, namely, the electric behaviors in solid materials at the macroscopic level are similar to the fluid hydrodynamics when the momentum-conserving electron-electron scattering plays the leading role, has got much attention in the past ten years. However, most of previous studies mainly focus on the electric properties. In this work, the thermal behaviors of hydrodynamic electron transport in a homogeneous 2D Corbino disk geometry is studied by the electron Boltzmann transport equation (eBTE) coupled with the Poisson equation under the magnetic field perpendicular to disk plane. Results show that in the electron hydrodynamic regime, the heat flux deflection phenomenon appears under the radial electric field or temperature gradient, namely, the heat flux no longer flows only along the radial direction and there is heat flux in the tangential direction of the radius. Heat flux deflection phenomenon is suppressed by momentum-relaxing scattering process and promoted by momentum-conserving scattering process. When an electric potential gradient or temperature gradient in the same direction is applied separately, the direction of heat flux is reversed in the electron hydrodynamic regime.
The question of a possible quantum critical point lying inside of a superconducting phase is central for understanding unconventional superconductivity. In various unconventional superconductors, non-Fermi-liquid/'strange-metal' $T^{n}$ variations, with $n<2$, of the electrical resistivity have been identified as the signature of magnetic quantum criticality. However, a difficulty is to prove experimentally that a non-Fermi-liquid/'strange-metal' law identified at temperatures above the superconducting temperature is the signature of an intrinsic zero-temperature quantum critical regime. In the heavy-fermion paramagnet UTe$_2$, unconventional superconductivity develops in the vicinity of a metamagnetic quantum phase transition induced by a magnetic field, and the quantum critical magnetic properties are suspected to play a role for the superconducting mechanism. In this work, we present a comparative analysis of electrical resistivity data collected on two UTe$_2$ samples of different qualities, in magnetic fields tilted by angles $θ\simeq35-40$~$^\circ$ from $\mathbf{b}$ to $\mathbf{c}$. Fits to the data have been performed either with a Fermi-liquid function $ρ=ρ_0+AT^{2}$ or with a non-Fermi-liquid/'strange-metal' function $ρ=ρ_0+A_nT^n$. Near to a superconducting phase induced beyond 40~T, non-physical residual resistivities $ρ_0<0$ are extracted from the $T^n$ fits, revealing that a 'hidden' Fermi-liquid $T^2$ regime may be ultimately recovered at low temperature. The results obtained here highlight the importance to investigate high-quality samples with low residual resistivity to confirm - or not - the presence of a suspected 'hidden' quantum critical behavior masked by superconductivity.
Altermagnets are a new class of magnets accompanying global time-reversal symmetry breaking (TRSB) without net magnetization. The TRSB results in formation of novel altermagnetic domains. Features of altermagnetic domains, in particular their responses to external stimuli, are essentially important but yet unexplored. Here, we report visualization of bulk altermagnetic domains in MnTe based on scanning magneto-optical Kerr-effect microscopy using telecom infrared wavelength. We found two distinct TRSB domains with large Kerr rotations that do not scale with its tiny bulk magnetization. We also revealed controllability and stability of domains against magnetic or thermal perturbations. Our first observation of altermagnetic domains using a laboratory-scale simple optical technique showing their movable nature provide firm bases for future fundamental and application studies of altermagnets.
Ruddlesden-Popper nickelates offer a new route to high-temperature superconductivity beyond the cuprates and iron-pnictides. However, the electronic reorganization that enables superconductivity in bilayer nickelates remain unresolved, largely due to the difficulty of directly probing the superconducting phase. Here, we overcome this limitation by stabilizing superconducting (La,Pr)$_3$Ni$_2$O$_7$ thin films with a protective capping layer, thereby enabling direct spectroscopic access via X-ray absorption and resonant inelastic X-ray scattering. We resolve the evolution of in-plane and out-of-plane electronic states, spin and orbital excitations, and spin-density-waves across insulating, superconducting, and metallic regimes. Combining experimental results with theoretical analysis, we show that the in-plane $d_{x^2-y^2}$ states form an itinerant backbone, while superconductivity emerges only when coherent $d_{z^2}$-$p_z$-$d_{z^2}$ interlayer hybridization develops, accompanied by suppressed static spin order and strongly damped spin excitations. Oxygen stoichiometry and epitaxial strain both act on this interlayer channel, placing superconductivity within a narrow window of interlayer coherence and correlation strength. These findings identify the microscopic ingredients required for superconductivity in bilayer nickelates and provide a multiorbital picture of its emergence.
We derive the reduced Dirac Hamiltonian on the non-rotating BTZ background and use its redshift structure to construct a gauge-covariant single-band lattice model on the constant-time BTZ cylinder. In equal-area coordinates the AdS radius $L$ fixes the local Gaussian curvature, while the horizon radius $r_h$ fixes the throat size and the strength of the near-horizon redshift. The lattice model therefore has a direct geometric interpretation and is not presented as an unshown reduction of the two-component Dirac lattice. Its angular Fourier transform yields an exact curved Harper equation with BTZ-dependent hopping amplitudes and a consistent dimensionless angular quasi-momentum. We then supplement global parameter scans with state-resolved diagnostics: spectra color-coded by mean radius, local density of states, direct flux-response versus radius correlations, and Aharonov--Bohm spectral flow and persistent current on the BTZ cycle. These results show that weaker curvature sharpens the butterfly-like fragmentation, whereas larger horizons suppress both magnetic and Aharonov--Bohm response by creating weakly dispersing near-horizon states.