arXiv Curator

The discovery of superconductivity in Ruddlesden-Popper nickelates has established a new frontier in the study of high-temperature superconductors. However, the underlying pairing mechanism and its relationship to the material's electronic and magnetic ground states remain elusive. Since unconventional superconductivity often emerges from a complex interplay of magnetic correlations, elucidating the magnetic ground state of the nickelates at ambient pressure is crucial for understanding the emergence of superconductivity under high pressure. Here, we combine high-resolution angle-resolved photoemission spectroscopy with tight-binding model simulation to investigate the electronic structure of the representative trilayer Ruddlesden-Popper nickelate La$_4$Ni$_3$O$_{10}$. We provide the first experimental evidence of band splitting induced by interlayer coupling and further resolve the momentum-dependent density wave gap structures along all the Fermi surfaces. Our findings identify the mirror-selective Fermi surface nesting as the origin of the interlayer antiferromagnetic spin density wave and demonstrate the dominant role of Ni-3d$_{z^2}$ orbitals in the low-energy physics of La$_4$Ni$_3$O$_{10}$. These results provide a fundamental framework for understanding the magnetic interactions and high-temperature superconductivity mechanism in the Ruddlesden-Popper nickelate family.

Quantum embedding methods, such as dynamical mean-field theory (DMFT), provide a powerful framework for investigating strongly correlated materials. A central computational bottleneck in DMFT is in solving the Anderson impurity model (AIM), whose exact solution is classically intractable for large bath sizes. In this work, we develop and benchmark a quantum-classical hybrid solver tailored for DMFT applications, using the variational quantum eigensolver (VQE) to prepare the ground state of the AIM with shallow quantum circuits. The solver uses a unified ansatz framework to prepare the particle and hole excitations of the ground-state from parameter-shifted circuits, enabling the reconstruction of the impurity Green's function through a continued-fraction expansion. We evaluate the performance of this approach across a few bath sizes and interaction strengths under noisy, shot-limited conditions. We compare three optimization routines (COBYLA, Adam, and L-BFGS-B) in terms of convergence and fidelity, assess the benefits of estimating a quantum-computed moment (QCM) correction to the variational energies, and benchmark the approach by comparing the reconstructed density of states (DOS) against that obtained using a classical pipeline. Our results demonstrate the feasibility of Green's function reconstruction on near-term devices and establish practical benchmarks for quantum impurity solvers embedded within self-consistent DMFT loops.

We introduce a tree tensor network (TTN) impurity solver that enables highly efficient and accurate real-time simulations of quantum impurity models. By decomposing a noninteracting bath Hamiltonian into a Cayley tree, the method provides a tensor network representation that naturally captures the multiscale entanglement structure intrinsic to impurity-bath systems. This geometry differs from conventional chain-based mappings and yields a substantial reduction of entanglement, allowing accurate ground-state properties and long-time dynamics to be captured at significantly lower bond dimensions. Benchmark calculations for the single-impurity Anderson model demonstrate that the TTN solver achieves markedly enhanced resolution of real-frequency spectral functions, without invoking analytic continuation. This impurity solver provides a balanced, scale-uniform description of impurity physics and offers a versatile approach for real-time dynamical mean-field theory and related applications involving quantum impurity models.

The difficulty of determining cation concentrations and oxygen stoichiometry in infinite-layer nickelate thin films has so far prevented clear experimental identification of the nickel electron configuration in the superconducting phase. We used soft x-ray absorption spectroscopy to study the successive changes in PrNiO$_x$ thin films at various intermediate stages of topotactic reduction with $x=2-3$. By comparing the Ni-$L$ edge spectra to single and double cluster ligand-field calculations, we find that none of our samples exhibit a pure $d^9$ configuration. Our quantitative analysis using the charge sum rule shows that even when films are maximally reduced, the averaged number of nickel $3d$ holes is 1.35. Superconducting samples have even higher values, calling into question the previously assumed limit of hole doping. Concomitant changes in the oxygen $K$-edge absorption spectra upon reduction indicate the presence of oxygen $2p$ holes, even in the most reduced films. Overall, our results suggest a complex interplay of hole doping mechanisms resulting from self-doping effects and oxygen non-stoichiometry.

Quantifying non-stabilizerness (``magic'') in interacting fermionic systems remains a formidable challenge, particularly for extracting high order correlations from quantum Monte Carlo simulations. In this Letter, we establish the two-point stabilizer Rényi entropy (SRE) and its mutual counterpart as robust, computationally accessible probes for detecting magic in diverse fermionic phases. By deriving local estimators suitable for advanced numerical methods, we demonstrate that these metrics effectively characterize quantum phase transitions: in the one-dimensional spinless $t$-$V$ model, they sharply identify the Luttinger liquid to charge density wave transition, while in the two-dimensional honeycomb lattice via determinant quantum Monte Carlo, they faithfully capture the critical exponents of the Gross-Neveu-Ising universality class. Furthermore, extending our analysis to the fractional quantum Hall regime, we unveil a non-trivial spatial texture of magic in the Laughlin state, revealing signatures of short-range exclusion correlations. Our results validate the two-point SRE as a versatile and sensitive diagnostic, forging a novel link between quantum resource theory, critical phenomena, and topological order in strongly correlated matter.

The discovery of nickelate superconductors has opened a new arena for studying the behavior of correlated electron liquids that give rise to unconventional superconductivity. While critical information about a material's charge dynamics is encoded in its plasmons, collective modes of the electron gas, these excitations have not yet been observed in nickelate materials. Here, we use resonant inelastic x-ray scattering (RIXS) to detect plasmons in the metallic, low-valence nickelate Pr4Ni3O8. Although qualitatively similar to those in cuprates, the nickelate plasmons are more heavily damped and have a lower velocity than those in a cuprate at comparable doping, which we attribute to reduced electronic hopping and enhanced screening of the long-range Coulomb interactions. Furthermore, the plasmons in Pr4Ni3O8 soften with increasing temperature, in contrast to the cuprate, where plasmons remain at nearly fixed energy but become more strongly damped. Taken together, these results reveal a distinct charge-screening landscape in nickelates and place quantitative constraints on analogies to cuprates.

Identifying superconducting states of matter without prior assumptions is a central challenge in strongly correlated electron systems. We introduce a canonical framework for diagnosing the formation of Cooper pair condensates based on the Penrose-Onsager criterion, in which superconducting order is encoded in the spectral properties of the two-particle reduced density matrix (2RDM). Within this formulation, the symmetry and structure of the condensate are obtained by projecting the 2RDM onto irreducible representations of the underlying symmetry group, enabling an unbiased identification of both conventional and exotic superconducting states. We demonstrate the power and versatility of the approach through applications to the two-dimensional Hubbard model, using both auxiliary-field quantum Monte Carlo (AFQMC) and the density matrix renormalization group (DMRG). For attractive interactions without a magnetic field, we reveal a clear finite-size scaling of the condensate fraction on square lattices of size up to $20\times 20$. The framework further provides direct access to the internal structure and extent of Cooper pairs, which we track across the BCS-BEC crossover. Moreover, it enables a clean diagnosis of the finite-momentum Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase in a magnetic field. Finally, we apply the approach to a supersolid phase in the repulsive Hubbard model with an additional next-nearest neighbor hopping $t^\prime$, where a charge-density wave coexists with a superconductor. We confirm the fragmented nature of the condensate and uncover substantial pairing correlations in the triplet channel with $p$-wave spatial symmetry in addition to the dominant singlet $d$-wave pairing. Our results establish the 2RDM-based Penrose-Onsager framework as a broadly applicable and unbiased tool for characterizing superconducting order in correlated quantum matter.

CsCr$_3$Sb$_5$ is a newly identified strongly correlated kagome superconductor, characterized by non-Fermi-liquid behavior at elevated temperatures and intertwined charge- and spin-density-wave order below $T_{DW}\approx 54$K. Under external pressure, this order is suppressed and a superconducting phase emerges. This phase diagram, which closely resembles that of high-$T_c$ superconductors, together with a kagome flat band near the Fermi level and possible altermagnetic order, has motivated extensive theoretical and experimental investigations. To better understand how pressure influences the ordered states, we present a systematic study of the evolution of the electronic properties under applied pressure. Performing DFT+DMFT (density functional theory combined with dynamical mean field theory) calculations, we uncover a complex interplay between the redistribution of spectral weight in the flat bands and the strength of electronic correlations under pressure. Our results further strengthen the interpretation that pressure effectively weakens electronic correlations through enhanced orbital hybridization. This, in turn, strongly suggests that superconductivity emerges as a direct consequence of the suppression of the system's ordered phase.

We present a time-domain iteration scheme for solving the Dynamical Mean-Field Theory (DMFT) self-consistent equations using retarded Green's functions in real time. Unlike conventional DMFT approaches that operate in imaginary time or frequency space, our scheme operates directly with real-time quantities. This makes it particularly suitable for near-term quantum computing hardware with limited Hilbert spaces, where real-time propagation can be efficiently implemented via Trotterization or variational quantum algorithms. We map the effective impurity problem to a finite one-dimensional chain with a small number of bath sites, solved via exact diagonalization as a proof-of-concept. The hybridization function is iteratively updated through time-domain fitting until self-consistency. We demonstrate stable convergence across a wide range of interaction strengths for the half-filled Hubbard model on a Bethe lattice, successfully capturing the metal-to-insulator transition. Despite using limited time resolution and a minimal bath discretization, the spectral functions clearly exhibit the emergence of Hubbard bands and the suppression of spectral weight at the Fermi level as interaction strength increases. This overcomes major limitations of two-site DMFT approximations by delivering detailed spectral features while preserving efficiency and compatibility with quantum computing platforms through real-time dynamics.

A polymorph of the bilayer nickelate La3Ni2O7 that displays an alternating single-layer (SL) and trilayer (TL; 1313) stacking pattern has recently been discovered. Signatures of superconductivity under pressure have been found in this phase. At ambient pressure, La3Ni2O7-1313 has been reported to crystallize in three different space-group symmetries Cmmm, Imma, and Fmmm. Unlike the commonly observed tilted NiO6 octahedra in perovskite nickelates, the Cmmm phase exhibits no NiO6 tilts, implying that this structural feature alone may be insufficient to give rise to superconductivity in Ruddlesden-Popper nickelates. Here, we employ first-principles calculations and group theory analysis to study the pressure dependence of the structural instabilities in this SL-TL La3Ni2O7 polymorph. At ambient pressure, we identify multiple unstable phonon branches in the highest symmetry (Cmmm) structure at various high-symmetry points of the Brillouin zone. Distortions associated with these instabilities lead to one of the other experimentally reported space groups (Imma) that does display octahedral tilts. The magnetic tendencies indicate that the electronic structure of La3Ni2O7-1313 at ambient pressure is dominated by the TL block, as the SL is in a Mott-insulating regime. Under pressure, a tetragonal P4/mmm structure becomes stable, in agreement with experiments.

We probe the density-wave transition of the trilayer nickelate La$_4$Ni$_3$O$_{10}$ with polarization-resolved infrared spectroscopy. The low-energy electrodynamics is strongly anisotropic, with metallic in-plane and insulating out-of-plane character. In the ordered phase, the anisotropy grows more than an order of magnitude as the out-of-plane conductivity is sharply suppressed. We interpret this enhancement as an effective electronic decoupling of the Ni-O layers, driven by a spin-density-wave-induced redistribution of Ni-$d_{z^2}$ occupation within the trilayers. This electronic response is accompanied by clear shifts and splittings of the out-of-plane phonons, compatible with a density-wave instability of electronic origin.

In this article we show that the condensation mechanism in cuprates involves the strong coupling of the condensate to pairon excited states. We present an accessible formalism that significantly extends our previous work, providing a theoretical basis for the energy-dependent gap function $Δ(E)$. The latter is proportional to the effective spin exchange energy, $J_{eff}$, with no retardation effects, such as the case of spin-fluctuation or phonon mediated couplings. The fundamental parameters of the superconducting (SC) state are the condensation energy per pair, $β_c$, and the antinodal energy gap, $Δ_p$, which are quantitatively extracted by fitting the cuprate quasiparticle spectrum from tunneling experiments. An explicit formula for the critical temperature is also derived in the model. Valid for any doping, we find $T_c$ to be proportional to $β_c$, and not the gap $Δ_p$, in sharp contrast to conventional SC. The numerical factor $β_c/k_BT_c\simeq 2.24$ originates from pair excitations of the condensate, following Bose statistics, with a mini-gap $δ_M \simeq 1\,$meV in the excitation spectrum. These results strongly suggest that the same `all-electron' mechanism is at work all along the $T_c$-dome.

We report a polarized ultrafast pump-probe study of the normal-state electronic dynamics in bilayer La$_3$Ni$_2$O$_7$ and trilayer La$_4$Ni$_3$O$_{10}$ single crystals at ambient pressure. While both nickelates exhibit density-wave (DW) transitions accompanied by the opening of a quasiparticle relaxation bottleneck, their electronic responses display strikingly different symmetry properties. La$_4$Ni$_3$O$_{10}$ maintains an isotropic optical response across the entire temperature range. In contrast, La$_3$Ni$_2$O$_7$ exhibits a pronounced twofold ($C_2$) anisotropy in its low-temperature electronic dynamics. This electronic nematicity, evident in both the relaxation dynamics and the effective gap scales, competes with a secondary isotropic order emerging below 115 K. The presence of macroscopic electronic anisotropy in the bilayer system, and its absence in the trilayer system, suggests an intimate relation between electronic nematic fluctuations and superconducting pairing in La$_3$Ni$_2$O$_7$ that worth for deeper explorations.

Altermagnetism combines fully compensated magnetic order with a magnetic symmetry that relates inequivalent spin sublattices, offering a promising, still underexplored platform for unconventional topological phases. Here we show that both isotropic tensile strain and electron localization, controlled by an effective Hubbard interaction $U_{\text{eff}}$, can act as efficient and systematic topological control parameters in the altermagnetic Weyl semimetal CrSb. While CrSb hosts Weyl fermions at equilibrium, modest tensile strain of 4-5% stabilizes additional symmetry allowed Dirac crossings and triple-point fermions, with further strain selectively favoring the triple-point phase. We propose a 3D low-energy Hamiltonian that captures the interplay between the Hubbard interaction $U$ and the sublattice symmetry of the altermagnet, giving rise to an interaction-driven Dirac crossing. Our results establish CrSb as a model altermagnet in which either strain or electron localization can selectively access and control the distinct topologies inherent to the altermagnets.

Owing to strong electronic correlations, high-temperature superconductivity always exhibits intricate intertwinement with various competing electronic orders in phase diagrams, such as spin/charge density waves (S/CDWs). In cuprate superconductors,the intertwinement of superconductivity and CDW order could strongly affect the fundamental properties of superconductivity, such as the critical temperature(Tc) and critical magnetic field(Hc). Recent high-field transport measurements indicate that when quantum fluctuations become important at low temperatures and high magnetic fields, the CDW order also reshapes the vortex states, which leads to fragile superconductivity with extremely low critical current(Jc). Here, by performing comprehensive high-field transport measurements, the H-T phase diagram of vortex states is mapped to H = 33 T in a quasi-two-dimensional FeSe-based superconductor (TBA+)xFeSe with a zero-resistivity transition temperature above 40 K. Our results indicate that (TBA+)xFeSe is an extremely type II superconductor with significant thermal fluctuations.At low temperatures, high magnetic fields cause the vortex solid state to exhibit similar current-dependent zero-resistance behavior as the fragile superconductivity in cuprate superconductors with CDW order. When the vortex solid state is melted with increasing temperature, a superconducting regime with vortex-like phase fluctuations emerges as an intermediate state, which features finite longitudinal resistance and vanishing Hall resistance. At higher temperatures, a vortex liquid state with finite Hall resistance eventually appears due to thermal fluctuations. All these observations suggest exotic vortex states beyond the classical paradigm of vortex matter.

NaCaNi$_2$F$_7$ is a unique example of spin-1 Heisenberg antiferromagnet on the pyrochlore lattice, but the presence of Na$^{1+}$/Ca$^{2+}$ $A$-site disorder complicates the local electronic and magnetic environment of the Ni$^{2+}$ $B$-site. We utilize resonant inelastic X-ray scattering (RIXS) to study the influence of $A$-site disorder on the $B$-site electronic structure of NaCaNi$_2$F$_7$. Ni L-edge RIXS measurements reveal a Ni$^{2+}$ electronic structure in nearly ideal octahedral coordination, with only a small trigonal compression ($δ$ = -200$\;$meV) required to capture all spectral features. Within the $D_{3d}$ symmetry of the Ni local environment, we extract an anisotropic $g$-factor of $g_{\parallel} = 2.26$ and $g_{\perp} = 2.27$, and a corresponding paramagnetic moment of $μ_{\rm{eff}}=3.2\;μ_B$. To simulate disorder, RIXS spectra were calculated with realistic distributions of crystal field parameters; however, these spectra are invariant relative to a disorder-free model, demonstrating the robustness of the Ni$^{2+}$ electronic environment to the $A$-site disorder, within the resolution of our measurement.

The density matrix renormalisation group (DMRG) is one of the most powerful numerical methods for strongly correlated condensed matter systems. We extend DMRG to the case with the $\mathrm{SU}(N_c)$ symmetry with $N_c > 2$, including two-dimensional systems. As a killer application, we simulate the ground state of the $\mathrm{SU}(4)$ Heisenberg model on the honeycomb lattice, which can potentially be realised in cold atomic systems and solid state systems like $α$-ZrCl$_3$. We keep up to 12800 $\mathrm{SU}(4)$ states equivalent to more than a million $\mathrm{U}(1)$ states. This supermassive DMRG simulation reveals the quantum spin-orbital liquid ground state, which has been conjectured for more than a decade. The methodology developed here can be extended to any classical Lie groups, paving the way to a next-generation DMRG with a full symmetry implementation.

Cuprate high-temperature superconductors, such as Yttrium Barium Copper Oxide (YBCO), are extremely promising for emerging technologies such as low-power computing, data storage, quantum sensors and superconducting electronics. However, the realization of high-performance functional nanostructures presents formidable challenges due to the difficulty of applying conventional nanofabrication methods to such sensitive materials, making the search for alternative methods a key enabling factor. Since YBCO's superconducting and normal-state properties are highly dependent on oxygen stoichiometry, precise nanoscale control of the oxygen content represents a highly appealing approach for creating advanced nanoengineered devices. In this work, we demonstrate the precise fabrication of sub-micrometer, grayscale patterns over large areas in epitaxial YBCO thin films, achieving finely tuned optical and superconducting transport properties by locally controlling the stoichiometry through maskless direct laser writing under ambient conditions. Cryogenic magneto-optical imaging and transport measurements in irradiated devices directly demonstrate the spatial tuning of the critical temperature and carrier density with the patterning conditions. Correlated Raman microscopy and reflectometry of the patterned regions indicate a laser-power dependent oxygen depletion in the irradiated regions. The proposed laser-controlled stoichiometry approach provides a direct and scalable method to navigate the phase diagram of high-Tc superconducting oxides, offering new possibilities for integrating functional nanostructures into superconducting devices.

We present an extension to the two-dimensional functional renormalization group to efficiently treat interactions on the surface or at interfaces of three-dimensional systems. As an application, we consider a semi-infinite stack of two-dimensional square lattices, including a Hubbard interaction on the surface layer and an alternating interlayer coupling. We investigate how strongly correlated states of the decoupled two-dimensional Hubbard model on the surface evolve under inclusion of such an SSH-like interlayer coupling. For large parts of the phase diagram as a function of the interlayer hopping parameters, the physics of the two-dimensional system prevails, with antiferromagnetic, superconducting d-wave, and ferromagnetic correlations taking center stage. However, for intermediate interlayer couplings the superconducting state at intermediate interaction strengths separates into two regimes by a small region of incommensurate spin-density-wave and spin-bond order, enabling the potential realization of chiral spin-bond order.

Correlated physics in nearly flat topological bands is a central theme in the study of moiré materials. While ground states at integer fillings are typically identified as quantum Hall ferromagnets within a Hartree-Fock framework, we propose the existence of symmetric topological Mott insulators (STMIs) that transcend this Slater determinant picture. Focusing on half-filling of each flavor per unit cell, we demonstrate the existence of STMIs which exhibit a quantized charge or spin Hall response. We first establish this phase in a bilayer Haldane-Hubbard model with localized orbitals on the $A$ sublattice and dispersive band on the $B$ sublattice. Starting from a trivial Mott insulator on the $A$ sublattice, tuning the sublattice potential drives a Bose-Einstein-condensation (BEC) to Bardeen-Cooper-Schrieffer (BCS) transition of the associated $p-\mathrm{i}p$ exciton pairing, realizing a topological Mott insulator with $C=1$ per flavor. We further generalize this construction to a single-layer spinful model, where the resulting STMI hosts charge edge modes coexisting with bulk local moments. A Mott semimetal is identified at the quantum critical point between the STMI and the trivial Mott insulator. Finally, we discuss applications to AA-stacked MoTe$_2$/WSe$_2$, proposing a ferromagnetic Chern insulator phase as a low-temperature descendant of the symmetric Mott semimetal.

To investigate the coexistence of superconductivity and charge density wave (CDW) in a correlated regime, we employ the Green's functions formalism, as well as the Hubbard-I approximation, as a way to introduce the correlations into the problem, in the form of a repulsive Coulomb interaction $U$. In addition, we investigate the effects of second-nearest neighbor hopping $t_1$ on a pure CDW state. The analysis of the results show that, for small values of $t_1$, both CDW and superconducting gaps compete for the same region on the Fermi surface. The increase of $t_1$ decreases the competition and may lead the system to a coexistence regime. Effects of temperature in the coexistence regime, are also investigated.

Quantum spin liquids (QSLs) host a variety of fractionalized particles. In Kitaev's paradigmatic honeycomb model a spin-$\tfrac{1}{2}$ fractionalizes into $Z_2$ flux due to emergent $Z_2$ gauge field and matter Majorana fermions. Although these excitations have well-defined dynamics in the integrable limit, their direct experimental identification is notoriously challenging: realistic materials inevitably host additional symmetry-allowed interactions that break integrability and hybridize gauge and matter sectors, while magnetic fields, which are often required to suppress competing order and stabilize a putative QSL regime, further entangle the responses of different fractionalized quasiparticles and may even drive the system into field-induced spin-liquid phases that are not adiabatically connected to the integrable limit. A prominent example is the quantum Majorana metal, in which the distinct dynamics of fractionalized Majorana fermions can become directly visible in scattering. This report highlights recent progress on these related questions: in which field-stabilized QSL regimes and nearby emergent phases, and under what conditions, can the response of a specific fractionalized quasiparticle be isolated and positively understood, thereby clarifying the existence and the experimental scope of putative spin liquids? We review the progress on these questions across Abelian, non-Abelian, and an emergent quantum phases under magnetic field that are not perturbatively connected to the integrable limit. We connect these field-induced dynamical phenomena to concrete experimental observables, relevant for neutron scattering, resonant inelastic X-ray scattering, and pump-probe spectroscopy that are capable of resolving specific types of fractionalized particles, including Majoranas and $Z_2$ fluxes.

An $S=1$ Haldane chain is a one-dimensional (1D) quantum magnet where strong fluctuations result in quantum disordered singlet ground state with a gapped excitation spectrum. The gap magnitude is primarily set by the dominant intrachain interaction ($J_\text{1D}$). An applied magnetic field closes the gap at $B_\text{c1}$ and drives the system into a gapless Tomonaga-Luttinger liquid (TLL) regime, followed by, at lower temperatures, a Bose-Einstein condensate (BEC) ground state, persisting up to $B_\text{c2} \propto 4 J_\text{1D}/gμ_B$. Almost all previously studied experimental realizations of such systems were based on transition-metal complexes which typically suffer from intrinsic anisotropies or large $J_\text{1D}$ values, limiting the access to the full theoretical phase diagram. We report a comprehensive study of TLL and BEC phases in the organic Haldane chain system 3,5-bis(N-tert-butylaminoxyl)-3'-nitrobiphenyl (BoNO). The absence of anisotropy and a moderate $J_\text{1D}$ enable exploration of the complete $B-T$ phase diagram. Through $^1$H nuclear magnetic resonance, combined with theoretical analysis, we characterize the TLL properties, map the BEC phase boundary $T_c (B)$, determine the associated critical exponent $ν\approx 0.66$ at $B_\text{c2}$, and demonstrate universal quasiparticle scaling in the quantum-critical regime. These results provide full experimental validation of theoretical predictions for field-induced phases in an $S=1$ Haldane chain, made over two decades ago.

Topology forms a cornerstone in modern condensed matter and statistical physics, offering a new framework to classify the phases and phase transitions beyond the traditional Landau paradigm. However, it is widely believed that topological properties are destroyed when the bulk energy gap closes, making it highly nontrivial to consider topology in gapless quantum critical systems. To address these challenges, recent advancements have sought to generalize the notion of topology to systems without a bulk energy gap, including quantum critical points and critical phases, collectively referred to as gapless symmetry-protected topological states. Extending topology to gapless quantum critical systems challenges the traditional belief in condensed matter physics that topological edge states are typically tied to the presence of a bulk energy gap. Furthermore, it suggests that topology plays a crucial role in classifying quantum phase transitions even if they belong to the same universality class, fundamentally enriching the textbook understanding of phase transitions. Given its importance, here we give a pedagogical review of the current progress of topological physics in quantum critical systems. We introduce the topological properties of quantum critical points and generalize them to stable critical phases, both for noninteracting and interacting systems. Additionally, we discuss further generalizations and future directions, including higher dimensions, nonequilibrium phase transitions, and realizations in modern experiments.

Altermagnetism (AM), a collinear magnetic phase with momentum-dependent spin splitting, is a promising candidate for strong magnetoelectric coupling. However, realizing direct and tunable coupling between ferroelectricity (FE) and AM within a single two-dimensional (2D) material remains an outstanding challenge. Here, based on first-principles calculations, we identify the distorted phase of monolayer Janus VOClBr as an intrinsic 2D FE-AM. This phase demonstrates robust magnetoelectric coupling, as evidenced by a complete reversal of momentum-space spin polarization upon FE switching, and further supported by spin texture analysis and the magneto-optical Kerr effect. Notably, the FE properties are highly strain-tunable: biaxial compression strain of -4% reduces the FE polarization switching barrier by approximately 87%, whereas a tensile strain of +3% induces a phase transition to an antiferromagnet. Leveraging the lock-in between the electrically controlled spin texture and the magneto-optical Kerr effect signal, we propose a non-volatile, polymorphic spintronic memory device featuring all-electrical writing and optical readout. This work establishes 2D FE-AMs as a versatile platform for coupled ferroic orders and paves the way for voltage-controlled, multifunctional spin-logic devices.

Entanglement entropy (EE) is a fundamental probe of quantum phases and critical phenomena, which was thought to reflect only bulk universality for a long time. Very recently, people realized that the microscopic geometry of the entanglement cut can induce distinct entanglement-edge modes, whose coupling to bulk critical fluctuations may alter the scaling of the EE. However, this perception is very qualitative and lacks quantitative consideration. Here, we investigate this problem through high-precision quantum Monte Carlo simulations combined with the analysis of scaling theory to build a quantitative understanding. By considering three distinct bipartitions corresponding to three surface criticality types, we reveal a striking dependence of the constant term γ on the microscopic cut at the quantum critical point. Notably, cuts that generate extra gapless edge modes yield a sign reversal in γ compared to those producing gapped edges. We explain this behavior via a modified scaling form that incorporates contributions from both bulk and surface critical modes. Furthermore, we demonstrate that the derivative of EE robustly extracts the bulk critical point and exponent ν regardless of the cut geometry, providing a reliable diagnostic of bulk universality in the presence of strong surface effects. Our work for the first time establishes a direct quantitative connection between surface criticality and entanglement scaling, challenging the conventional view that EE solely reflects bulk properties and offering a refined framework for interpreting entanglement in systems with boundary-sensitive criticality.

The recent discovery of high-temperature superconductivity in pressurized Ruddlesden-Popper nickelates stimulated intense research into their correlated electron physics. Establishing the diversity of ground states across different Ruddlesden-Popper phases is crucial for elucidating the superconducting mechanisms in these nickelates. Motivated by the recent report of superconductivity in hybrid 1212-type La5Ni3O11, we synthesized and investigated the long-range-ordered hybrid 1313-type La3Ni2O7. In contrast to its bilayer counterpart, the 1313-type La3Ni2O7 exhibits characteristic semiconducting behavior at ambient pressure, displaying a distinct anomaly at 170 K. This behavior is consistently evidenced by measurements of both magnetic susceptibility and specific heat. Nuclear magnetic resonance spectroscopy unambiguously indicates a spin-density-wave transition occurring at 170 K. High-pressure electrical transport measurements demonstrate the induction of metallization under pressure, yet reveal no discernible traces of superconductivity up to 65 GPa. Our findings establish hybrid 1313-type La3Ni2O7 as a new member of the Ruddlesden-Popper nickelate family exhibiting a distinct spin-density-wave transition, and offers a new platform for investigating the interplay among crystal structure, electronic orders, and superconductivity in hybrid nickelates.

Spin fluctuations have been generally believed as the pairing glue of high-$T_c$ superconductivity. Recent inelastic neutron scattering (INS) studies have revealed a weak flat spin-fluctuation signal around 45 meV in the bilayer nickelate La$_3$Ni$_2$O$_{7-δ}$, suggesting strong interlayer and weak intralayer magnetic couplings ($SJ_{\perp}\approx$ 60 meV, $SJ_{\parallel}\leq$ 3.5 meV) in contrast to cuprate and pnictide superconductors. Here, we report further INS studies on the Pr and Nd doped La$_3$Ni$_2$O$_{7-δ}$ powder samples at ambient pressure. Besides the crystalline electric field excitations at low energies, we have found that the 45 meV flat mode splits into two modes in doped compounds, along with another weak mode at about 60 meV, where the spin fluctuations in La$_2$NdNi$_2$O$_{7-δ}$ are stronger than La$_3$Ni$_2$O$_{7-δ}$ and La$_2$PrNi$_2$O$_{7-δ}$. Based on an effective Heisenberg model by only considering the nearest-neighbor exchange couplings on the stripe-type antiferromagnetic orders, we conclude that the interlayer coupling $SJ_{\perp}$ is enhanced to about 69 meV and 73 meV for Pr and Nd doped samples, respectively. Our results highlight the crucial role of interlayer coupling in the rare-earth doped bilayer nickelates, which towards to promote high $T_c$ via interlayer $s\pm$ pairing.

We report a combined x-ray absorption (XAS) and oxygen $K$-edge resonant inelastic x-ray scattering (RIXS) study of charge order in underdoped La$_{2-x}$Sr$_x$CuO$_4$ ($x=0.125$) under uniaxial $c$-axis pressure. We find that compressive $c$-axis strain modifies the charge order only within the superconducting state, with a striking polarization dependence: suppression in the $d_{x^2-y^2}$ channel and enhancement in the $d_{z^2}$ channel. X-ray absorption spectra reveal concomitant strain-induced modifications of the oxygen pre-edge and upper Hubbard band, consistent with increased $d_{z^2}$ orbital admixture. Our results suggest that $c$-axis pressure drives an orbital separation between superconductivity, rooted in $d_{x^2-y^2}$ states, and charge order, which gradually shifts to the $d_{z^2}$ channel. This orbital separation reveals a way for superconductivity and charge order to coexist in the cuprates with minimal competition. Furthermore, it suggests that the multi-order phase diagram of La$_{2-x}$Sr$_{x}$CuO$_4$ cannot be realistically described within single band models usually used to describe cuprate physics.

Understanding unconventional superconductivity (SC) driven by strong electronic correlations is a central challenge in condensed matter physics. In this work, we employ sign-problem-free quantum Monte Carlo (QMC) simulations to systematically investigate a bilayer fermionic model featuring strong interlayer antiferromagnetic (AFM) exchange and on-site repulsive Hubbard interactions. This system serves as a prototypical model for realizing a symmetric mass generation (SMG) insulator. Our numerically exact results unambiguously demonstrate that robust superconducting pairing emerges upon doping the SMG phase. Remarkably, we find that the SC order is significantly enhanced by the repulsive Hubbard interaction. Given its potential relevance to the essential features of the high-$T_c$ superconductor $\mathrm{La}_{3}\mathrm{Ni}_{2}\mathrm{O}_{7}$ under pressure, our study establishes a new paradigm for superconductivity arising from a doped SMG parent state and provides key theoretical guidance for future experimental investigations.

Normalizing flows have recently demonstrated the ability to learn the Boltzmann distribution of the Hubbard model, opening new avenues for generative modeling in condensed matter physics. In this work, we investigate the steps required to extend such simulations to larger lattice sizes and lower temperatures, with a focus on enhancing stability and efficiency. Additionally, we present the scaling behavior of stochastic normalizing flows and non-equilibrium Markov chain Monte Carlo methods for this fermionic system.

We present comprehensive inelastic neutron spectroscopic maps of the magnetic field-induced disordered phase of the Kitaev quantum spin liquid candidate material $α$-RuCl$_3$. For fields along both in-plane high-symmetry directions we observe that the spin excitation spectrum at and above a magnetic field of 8~T is gapped. Excitation modes then sharpen for increasing field but are consistently broader than experimental resolution even at 13.5~T. The out-of-plane dispersion diminishes in the 7-10~T regime, signifying enhanced two-dimensional behavior as the in-plane liquid correlations are established. In this regime, excitations are very broad and largely flat for all accessible energy-momenta, which is kinematically at odds with a magnon-decay picture. By contrast, a continuum of fractionalized excitations naturally yields a broad continuum response, which crucially may be accompanied by sharper modes of bound states of fractionalized excitations. Their damping by the continuum accounts for the observed spectral broadening and field dependence. Our results provide strong evidence for the existence of fractionalized excitations in $α$-RuCl$_3$ in a magnetic field.

The discovery of superconducting transitions in pressurized La$_3$Ni$_2$O$_{7}$ and La$_4$Ni$_3$O$_{10}$ has highlighted the pivotal role of density wave (DW) orders in nickelate superconductors. To gain a comprehensive understanding of the superconducting state, it is essential to elucidate the nature of the DW order. In this study, we utilized $^{139}$La nuclear quadrupole resonance (NQR) to investigate the charge density wave (CDW) and spin density wave (SDW) orders in both single-crystal and polycrystalline La$_4$Ni$_3$O$_{10}$. Near $T_{\rm{DW}} \approx 133$ K, an abrupt change in both the linewidth and frequency of the La(2) site in the single-crystal sample provides compelling evidence for a first-order-like phase transition. The pronounced broadening of the NQR lines indicates the incommensurate nature of the DW order. Furthermore, the spin-lattice relaxation rate divided by temperature 1/$T_1$$T$ exhibits a strong enhancement at $T_{\rm{DW}}$, indicating the strong spin fluctuations above the first-order DW transition. These observations suggest an intricate interplay between incommensurate CDW and SDW orders. Our findings offer critical insights into the microscopic mechanisms of the DW state in La$_4$Ni$_3$O$_{10}$ and establish an essential framework for exploring the interplay between DW and superconducting phases in nickelate superconductors.

Since the discovery of high-temperature superconductivity, studying the upper critical field and its anisotropy has been crucial for understanding superconducting mechanism and guiding applications. Here we perform in situ high-pressure angular-dependent electrical transport measurements on Pr4Ni3O10 single crystals using a custom diamond anvil cell (DAC) rotator and confirming its anisotropic superconductivity. The anisotropy parameter is approximately 1.6, decreasing with increasing temperature and approaches 1 near Tc. Comparing effective mass anisotropy and inter-block distance in cuprates and iron-based superconductors (FeSCs) reveals that Pr4Ni3O10 single crystals superconductors are consistent with a two-band model, where intralayer quantum confinement within the unit cell induces interlayer coherence, thereby leading to three-dimensional (3D) superconductivity. This study not only establishes the existence of anisotropic superconductivity in bulk Ruddlesden-Popper nickelates, but also provide critical insight into the role of dimensionality in high-temperature superconductivity.

We present a theoretical study of the low-energy physics of a quarter-hole-filled two-orbital bilayer Hubbard model motivated by transition-metal bilayer systems with strong orbital-selective interlayer hybridization. By explicitly treating the strong interlayer bonding of dz2 orbitals within a molecular orbital basis and projecting out high-energy electronic states, we derive a low-energy effective Kugel-Khomskii Hamiltonian describing the interplay between electron spin and emergent layer pseudospin degrees of freedom. We map out a rich ground state phase diagram featuring diverse spin and charge ordered states. These include conventional ferromagnetic and antiferromagnetic phases with layer staggered charge densities, a layer-coherent phase characterized by spontaneous interlayer quantum coherence, and a novel maximally spin-layer-entangled phase with a hidden composite spin-layer order. We show that this exotic hidden ordered phase arises from the spontaneous breaking of an emergent O(4) symmetry down to a O(3), manifesting a unique excitation spectrum with three entangled gapless Goldstone modes. Our results uncover a geometrically driven mechanism for realizing composite entanglement in strongly correlated bilayer systems and provide a concrete theoretical framework relevant to bilayer nickelate superconductors and other multi-component correlated materials.

We present a finite-temperature canonical-ensemble determinant quantum Monte Carlo algorithm that enforces an exact fermion number and enables stable simulations of correlated lattice fermions. We propose a stabilized QR update that reduces the computational complexity from standard cubic scaling $O(βN^3)$ to linear scaling $O(βN N_e^2)$ with respect to the system size $N$, where $N_e$ is the particle number. This yields a dramatic speedup in dilute regimes ($N_e \ll N$), opening unbiased access to large-scale simulations of strongly correlated low-density phases. We validate the method on the dilute Fermi gas with onsite Hubbard interactions, observing the suppression of the fermion sign problem in the dilute limit. Furthermore, we apply this approach to an one-dimensional flat-band system, where the canonical ensemble allows for precise control over filling. We reveal a ferromagnetic instability at low temperatures in the half-filling regime. Our linear-scaling approach provides a powerful tool for investigating emergent phenomena in dilute quantum matter.

Nonequilibrium electronic forces play a central role in voltage-driven phase transitions but are notoriously expensive to evaluate in dynamical simulations. Here we develop a machine learning framework for adiabatic lattice dynamics coupled to nonequilibrium electrons, and demonstrate it for a gating induced insulator to metal transition out of a charge density wave state in the Holstein model. Although exact electronic forces can be obtained from nonequilibrium Green's function (NEGF) calculations, their high computational cost renders long time dynamical simulations prohibitively expensive. By exploiting the locality of the electronic response, we train a neural network to directly predict instantaneous local electronic forces from the lattice configuration, thereby bypassing repeated NEGF calculations during time evolution. When combined with Brownian dynamics, the resulting machine learning force field quantitatively reproduces domain wall motion and nonequilibrium phase transition dynamics obtained from full NEGF simulations, while achieving orders of magnitude gains in computational efficiency. Our results establish direct force learning as an efficient and accurate approach for simulating nonequilibrium lattice dynamics in driven quantum materials.

Interaction quenches in strongly correlated electron systems provide a powerful route to probe nonequilibrium many-body dynamics. For the Hubbard model, nonequilibrium dynamical mean-field theory has revealed coherent post-quench oscillations, dynamical crossovers, and long-lived transient regimes. However, these studies are largely restricted to spatially homogeneous dynamics and therefore neglect the role of spatial structure formation during ultrafast evolution. Here we investigate interaction quenches in the half-filled Hubbard model using a real-space time-dependent Gutzwiller framework. We show that homogeneous nonequilibrium dynamics is generically unstable: even arbitrarily weak spatial fluctuations grow dynamically and drive the system toward intrinsically inhomogeneous states. Depending on the interaction strength, the post-quench evolution exhibits spatial differentiation, nucleation, and slow coarsening of Mott-like domains. Our results establish spatial self-organization as a generic feature of far-from-equilibrium correlated matter and reveal a fundamental limitation of spatially homogeneous nonequilibrium theories.

Studying finite-temperature properties with tensor networks is notoriously difficult, especially at low temperatures, due to the rapid growth of entanglement and the complexity of thermal states. Existing methods like purification and minimally entangled typical thermal states offer partial solutions but struggle with scalability and accuracy in low-temperature regime. To overcome these limitations, we propose a new approach based on generating-function matrix product states (GFMPS). By directly computing a large set of Bloch-type excited states, we construct Gibbs states that moderate the area-law constraint, enabling accurate and efficient approximation of low-temperature thermal behavior. Our benchmark results show magnificent agreement with both exact diagonalization and experimental observations, validating the accuracy of our approach. This method offers a promising new direction for overcoming the longstanding challenges of studying low-temperature properties within the tensor network framework. We also expect that our method will facilitate the numerical simulation of quantum materials in comparison with experimental observations.

Unconventional superconductivity under strong Zeeman fields--particularly beyond the Pauli paramagnetic limit--remains a central challenge in condensed matter physics. The exotic Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state, in particular, remains in need of definitive study within fundamental electronic models. Here we employ state-of-the-art finite-temperature and ground-state tensor network approaches to systematically explore the superconducting (SC) phase diagram of the $t$-$t'$-$J$ model subjected to Zeeman fields. We find that zero-momentum $d$-wave superconductivity persists until the spin gap closes, coexisting with charge density waves. A novel $d$-wave FFLO phase emerges under a higher Zeeman field even above the Pauli limit, concomitant with a field-enhanced spin density waves. We identify these phases, characterized by the simultaneous presence of pairing condensate and density wave orders, as charge-2e supersolids. Analysis of Matsubara Green's function reveals that the FFLO pairing momentum is locked to the underlying Fermi surface. Our results provide microscopic insights into field-induced unconventional pairing mechanisms and reveal the long-sought FFLO state in a fundamental correlated electron model, offering a promising route for its realization in ultracold atom optical lattice.

We introduce SeeMPS, a Python library dedicated to implementing tensor network algorithms based on the well-known Matrix Product States (MPS) and Quantized Tensor Train (QTT) formalisms. SeeMPS is implemented as a complete finite precision linear algebra package where exponentially large vector spaces are compressed using the MPS/TT formalism. It enables both low-level operations, such as vector addition, linear transformations, and Hadamard products, as well as high-level algorithms, including the approximation of linear equations, eigenvalue computations, and exponentially efficient Fourier transforms. This library can be used for traditional quantum many-body physics applications and also for quantum-inspired numerical analysis problems, such as solving PDEs, interpolating and integrating multidimensional functions, sampling multivariate probability distributions, etc.

Pump-probe experiments performed on $\mathrm{K}_3\mathrm{C}_{60}$ have unveiled both optical and transport signatures of metastable light-induced superconductivity up to room temperature, far above $T_c$. Recent experiments have uncovered that excitation in the vicinity of $50 ~\textrm{meV}$ enables the observation of high temperature light-induced superconductivity at significantly lower fluences. Inspired by these experiments we develop a mechanism which can explain such a giant resonant enhancement of light-induced superconductivity. Within a minimal non-linear Holstein model, we show that resonantly driving optical Raman modes leads to a time-dependent electron-phonon coupling. Such a coupling then modulates the effective electron-electron attraction, with the strongest modulations occurring when the drive is resonant with the phonon frequency. These dynamical modulations of the pairing interactions lead to Floquet-BCS instabilities at temperatures far exceeding equilibrium $T_c$, as observed in experiments. We conclude by discussing the implications of our general analysis on the $\mathrm{K}_3\mathrm{C}_{60}$ experiments specifically and suggesting experimental signatures of our mechanism.

We solve the two-electron problem on a square lattice with $d$-wave altermagnetism, considering both on-site and nearest-neighbor attractive interactions. The altermagnetic spin-splitting in the single-particle dispersion naturally gives rise to a ground state of two-electron bound pairs with nonzero center-of-mass momentum. This finite-momentum pairing can be interpreted as a two-body mechanism underlying the recently proposed altermagnetism-induced Fulde--Ferrell--Larkin--Ovchinnikov (FFLO) superconducting state. Additionally, when the nearest-neighbor attraction is strong, the resulting finite-momentum bound pairs exhibit a mixture of both spin-singlet and spin-triplet characteristics, suggesting the possibility of unconventional superconductors, where spin-singlet and spin-triplet pairings coexist.

Recent research suggests the possibility of the two-dimensional breathing-Kagome magnet Nb$_3$Cl$_8$ hosting a quantum spin liquid state, warranting further study into its magnetic properties. Using ab initio calculations, we show that monolayer Nb$_3$Cl$_8$ has short-range antiferromagnetic correlations among Nb$_3$ trimers with S = 1/2, and becomes magnetically frustrated due to the underlying effective triangular lattice geometry, and is evidenced by a frustration index of f > 1. The high-temperature susceptibility shows a negative Weiss temperature from Monte Carlo calculations. Considering spin-orbit coupling, we investigate the magnetic anisotropy, including anisotropic exchange, single-ion anisotropy and the Dzyaloshinskii-Moriya interaction using the four-state energy mapping formalism. Although the elements have relatively small atomic numbers, the Dzyaloshinskii-Moriya interaction is comparable in magnitude to the anisotropic exchange. Additionally, we show that biaxial strain tunes the short-range correlations between antiferromagnetic, paramagnetic and ferromagnetic. These findings strengthen our understanding of Nb$_3$Cl$_8$ and advance its applications in current condensed matter physics and materials science research, including nanoscale mechanical and spintronics applications.

Van Hove singularities (VHSs) amplify electronic correlations, providing a crucial platform for discovering novel quantum phase transitions. Here, we show that VHSs in multi-orbital systems can stabilize a variety of competing $\bm{Q}=0$ magnetic orders, including intrinsic altermagnetism emerging from spontaneous orbital antiferromagnetism. This intrinsic phase, in which antiparallel spins reside on distinct orbitals, is realized across all four 2D Bravais lattices. It is driven by orbital-resolved spin fluctuations enhanced by inter-orbital hopping and favors suppressed Hund's coupling $J_H$, strong inter-orbital hybridization, and filling near a VHS from quadratic band touching. Through Hubbard-$U$-$J_H$ phase diagrams we map several magnetic phase transitions: (i) ferrimagnet to $d$-wave extrinsic altermagnet, (ii) $d$-wave intrinsic altermagnet to ferromagnet, and (iii) $g$-wave extrinsic altermagnet to either $d$-wave extrinsic altermagnet or ferromagnet. Our work identifies VHSs as a generic route to altermagnetism in correlated materials.

Scanning tunneling microscopy measurements exploiting the natural superstructure modulation of the cuprate superconductor Bi$_2$Sr$_2$CaCu$_2$O$_{8+x}$ (Bi-2212) have revealed a possible correlation between the Cu-apical-O distance $δ_{\mathrm{api}}$ and the superconducting order parameter $m_{\mathrm{SC}}$, as reported recently by O'Mahony et al. (Proc. Natl. Acad. Sci. 119, e2207449119 (2022)). These observations were interpreted as evidence for a direct link between superconductivity and the charge-transfer gap, and more broadly revived the long-standing question of the role of apical oxygens in cuprate superconductivity. Using a combination of density-functional theory and cluster dynamical mean-field theory, we compute from first principles the variations of $m_{\mathrm{SC}}$ induced solely by apical oxygen displacement in Bi$_2$Sr$_2$CuO$_{6+δ}$, Bi-2212, and HgBa$_2$CuO$_{4+δ}$. The quantitative agreement between our calculations and experiments allows us to unambiguously attribute the observed variations of $m_{\mathrm{SC}}$ to changes in $δ_{\mathrm{api}}$. We demonstrate, however, that these variations of $m_{\mathrm{SC}}$ originate predominantly from changes in the effective hole-doping of the CuO$_2$ planes, with negligible effect on the charge-transfer gap. The modest magnitude of the $m_{\mathrm{SC}}$ modulation induced by apical-oxygen displacement alone therefore warrants caution in interpreting correlations between $T_c$ and $δ_{\mathrm{api}}$ inferred from comparisons across different cuprate compounds.

Several ground state properties of the square-lattice $S=1/2$ Heisenberg antiferromagnet are computed (the energy, order parameter, spin stiffness, spinwave velocity, long-wavelength susceptibility, and staggered susceptibility) using extensive quantum Monte Carlo simulations with the stochastic series expansion method. Moderately sized lattices are studied at temperatures $T$ sufficiently low to realize the $T \to 0$ limit. Results for periodic $L\times L$ lattices with $L \in [6,96]$ are tabulated versus $L$ and extrapolations to infinite system size are carried out. The extrapolated ground state energy density is $e_0=-0.669441857(7)$, which represents an improvement in precision of three orders of magnitude over the previously best result. The leading and subleading finite-size corrections to $e_0$ are in full quantitative agreement with predictions from chiral perturbation theory, thus further supporting the soundness of both the extrapolations and the theory. The extrapolated sublattice magnetization is $m_s=0.307447(2)$, which agrees well with previous estimates but with a much smaller statistical error. The coefficient of the linear in $L^{-1}$ correction to $m^2_s$ agrees with the value from chiral perturbation theory and the presence of a factor $\ln^γ(L)$ in the second-order correction is also confirmed, with the previously not known value of the exponent being $γ= 0.82(4)$. The finite-size corrections to the staggered susceptibility point to logarithmic corrections also in this quantity. To facilitate benchmarking of methods for which periodic boundary conditions are challenging, results for systems with open and cylindrical boundaries are also listed and their spatially inhomogeneous order parameters are analyzed.

Solving the quantum many-body ground state problem remains a central challenge in computational physics. In this context, the Variational Monte Carlo (VMC) framework based on Projected Entangled Pair States (PEPS) has witnessed rapid development, establishing itself as a vital approach for investigating strongly correlated two-dimensional systems. However, standard PEPS-VMC algorithms predominantly rely on sequential local updates. This conventional approach often suffers from slow convergence and critical slowing down, particularly in the vicinity of phase transitions or within frustrated landscapes. To address these limitations, we propose an efficient autoregressive row-wise sampling algorithm for PEPS that enables direct, rejection-free sampling via single-layer contractions. By utilizing autoregressive single-layer row updates to generate collective, non-local configuration proposals, our method significantly reduces temporal correlations compared to local Metropolis moves. We benchmark the algorithm on the two-dimensional transverse-field Ising model and the quantum spin glass. Our results demonstrate that the row-wise scheme effectively mitigates critical slowing down near the Ising critical point. Furthermore, in the rugged landscape of the quantum spin glass, it yields improved optimization stability and lower ground-state energies. These findings indicate that single-layer autoregressive row updates provide a flexible and robust improvement to local PEPS-VMC sampling and may serve as a basis for more advanced sampling schemes.

Jahn--Teller (JT) distortions are a key driver of physical properties in many correlated oxide materials. Cooperative JT distortions, in which long-range orbital order reduces the symmetry of the average structure macroscopically, are common in JT-distorted materials at low temperatures. This long-range order will often melt on heating, \textit{via} a transition to a high-temperature state without long-range orbital order. The nature of this transition has been observed to vary with different materials depending on crystal structure; in LaMnO$_3$ the transition has generally been interpreted as order-disorder, whereas in layered nickelates $A$NiO$_2$ ($A$=Li,Na) there is a displacive transition. Alternatively, recent theoretical work has suggested that previous attributions of order-disorder may in fact be a consequence of phonon anharmonicity, rather than persistence of JT distortions, which would suggest that the displacive transition may be more common than currently believed. In this work, we run Monte Carlo simulations with a simple Hamiltonian which is modified to include terms dependent on the JT amplitude $ρ$, which is allowed to vary within the simulation \textit{via} the Metropolis algorithm. Our simulations yield distributions of JT amplitudes consistent with displacive rather than order-disorder behaviour for both perovskites and layered nickelates, suggesting that displacive-like JT transitions may be more common than previously assumed in both perovskites and layered nickelates. We also find significant differences between the transition observed for perovskites compared with layered nickelates, which we attribute to differing extensivity of configurational entropy on the two lattices, showing the crucial role of lattice geometry in determining behaviour.

In trapped-ion quantum simulators, effective spin-1/2 XY interactions can be engineered via laser-induced coupling between internal atomic states and collective phonon modes. In the simplest one-dimensional ($1d$) traps, these interactions decay as a power-law with distance $1/r^α$, with a tunable exponent $α$. For small $α$, the resulting long-range $1d$ XY model exhibits continuous symmetry breaking, in marked contrast to its nearest neighbor counterpart. In this paper, we examine this model near the phase transition at $α_c$ from the lens of the spin stiffness, or superfluid density. We develop a stochastic series expansion (SSE) quantum Monte Carlo (QMC) simulation and a generalized winding number estimator to measure the superfluid density in the presence of power-law interactions, which we test against exact diagonalization for small lattice sizes. Our results show how conventional superfluidity in the $1d$ XY model is enhanced in the long-range interacting regime. This is observed as a diverging superfluid density as $α\rightarrow 0$ in the thermodynamic limit, which we show is consistent with linear spin-wave theory. Finally, we define a normalized superfluid density estimator that clearly distinguishes the short, medium, and long-range interacting regimes, providing a novel QMC probe of the critical value $α_c$.

We present a systematic comparison of the real-time Dyson expansion (RTDE) with established non-equilibrium Green's function approaches for simulating driven, interacting quantum systems. Focusing on density matrix dynamics, time-off-diagonal Green's functions, and time-resolved photoemission spectra, we benchmark RTDE against fully self-consistent Kadanoff-Baym equation (KBE) calculations, the generalized Kadanoff-Baym ansatz (GKBA), and exact diagonalization for small systems using second order many-body perturbation theory. Using a driven two-band Hubbard model, we show that mean-field single particle density matrix trajectories provide a reliable baseline for RTDE across a broad range of interaction strengths and excited-carrier populations. Further, RTDE accurately captures correlation effects in the Green's functions, including long-lived oscillations and revivals that are strongly suppressed by the overdamping inherent to self-consistent KBE schemes. As a consequence, RTDE resolves rich non-equilibrium spectral structure in time-resolved photoemission, such as interaction- and population-dependent quasiparticle splittings and bandgap renormalization, which are largely washed out in self-consistent approaches, yet are present in the exact solutions. Our results demonstrate that RTDE bridges the gap between mean-field propagation and full two-time KBE simulations, retaining favorable linear scaling while capturing essential dynamical correlations relevant for ultrafast spectroscopy.

Neural quantum states (NQS) have proven highly effective in representing quantum many-body wavefunctions, but their loss landscape remains poorly understood and debated. Here, we demonstrate that the NQS loss landscape is more benign and similar to conventional deep learning than previously thought, exhibiting mode connectivity: independently trained NQS are connected by paths in parameter space with essentially no energy barrier. To construct these paths, we develop GeoNEB, a path optimizer integrating efficient stochastic reconfiguration with the nudged elastic band method for constructing minimum energy paths. For the strongly interacting six-electron quantum dot modeled by a $1.6$M-parameter Psiformer, we find two independent minima with expected energy barrier $\sim10^{-5}$ times smaller than the system's overall energy scale and $\sim10^{-3}$ times smaller than the linear path's barrier. The path respects physical symmetry in addition to achieving low energy, with the angular momentum remaining well quantized throughout. Our work is the first to construct optimized paths between independently trained NQS, and it suggests that the NQS loss landscape may not be as pathological as once feared.

Kitaev iridates have been proposed as candidates for realizing an elusive quantum spin liquid (QSL) state, in which strong spin-orbit coupling and bond-directional exchange generate a highly frustrated and entangled ground state. However, all physical systems proposed to host this ground state, including Li$_2$IrO$_3$, Na$_2$IrO$_3$, and RuCl$_3$, develop magnetic order at low temperatures due to competing interactions. Nonetheless, theoretical modeling of experimental data has shown that Kitaev interactions are still present, motivating the application of perturbations such as pressure, magnetic field, and chemical doping to drive the system into the QSL phase. Here we study $β$-Li$_{2}$Ru$_{x}$Ir$_{1-x}$O$_{3}$ with dilute levels of Ru, $x \lesssim 10\%$. Through a combination of magnetometry, resonant elastic X-ray scattering, ac-heat capacity, and muon spin relaxation/resonance, we show that weak magnetic disorder suppresses long-range antiferromagnetic order and stabilizes an anisotropic spin glass that retains key signatures of Kitaev exchange. This Kitaev spin glass shows pronounced directional anisotropy in its magnetic susceptibility and thermoremenant magnetization. These results demonstrate that dilute magnetic disorder can access an anisotropic Kitaev spin glass: a proximate phase that freezes the Kitaev frustration landscape. This could provide a new window into the degeneracy, anisotropy, and competing interactions underlying the Kitaev QSL.

Kitaev iridates have been proposed as candidates for realizing an elusive quantum spin liquid (QSL) state, in which strong spin-orbit coupling and bond-directional exchange generate a highly frustrated and entangled ground state. However, all physical systems proposed to host this ground state, including Li$_2$IrO$_3$, Na$_2$IrO$_3$, and RuCl$_3$, develop magnetic order at low temperatures due to competing interactions. Nonetheless, theoretical modeling of experimental data has shown that Kitaev interactions are still present, motivating the application of perturbations such as pressure, magnetic field, and chemical doping to drive the system into the QSL phase. Here we study $β$-Li$_{2}$Ru$_{x}$Ir$_{1-x}$O$_{3}$ with dilute levels of Ru, $x \lesssim 10\%$. Through a combination of magnetometry, resonant elastic X-ray scattering, ac-heat capacity, and muon spin relaxation/resonance, we show that weak magnetic disorder suppresses long-range antiferromagnetic order and stabilizes an anisotropic spin glass that retains key signatures of Kitaev exchange. This Kitaev spin glass shows pronounced directional anisotropy in its magnetic susceptibility and thermoremenant magnetization. These results demonstrate that dilute magnetic disorder can access an anisotropic Kitaev spin glass: a proximate phase that freezes the Kitaev frustration landscape. This could provide a new window into the degeneracy, anisotropy, and competing interactions underlying the Kitaev QSL.

The Kondo effect is a paradigmatic model of strongly-correlated physics, where a magnetic impurity forms a many-body singlet with a fermionic environment. Cold gases of ytterbium (Yb) atoms have been proposed to be an ideal platform to study the Kondo effect since different internal states of the atom can be used to create both the impurity and the fermionic environment. In Yb gases, however, the atomic impurity interacts with the fermionic environment both through magnetic and potential scattering. These two scattering mechanisms counteract one another, raising the question of how robust Kondo screening remains. Here, we show that potential scattering can quench the Kondo screening in one-dimensional Yb gases; yet, strikingly, Kondo physics survives this quench in well-defined regimes. Combining analytical renormalization-group theory for a Luttinger liquid with density matrix renormalization group (DMRG) simulations, we identify a transition from a strongly- to a weakly-entangled impurity as potential scattering is increased. The two approaches show excellent agreement concerning the stability of Kondo physics throughout the different parameter regimes considered. Our results provide a quantitative criterion for the emergence of Kondo screening in one-dimensional Yb gases and delineate experimentally accessible regimes for its realization in cold-atom platforms.

For quantum spin models defined on a two-dimensional lattice, we look for the best numbering of the lattice sites (a layout) that, at fixed bond dimension and other parameters of the density matrix renormalization group (DMRG) algorithm, gives the lowest value of the variational energy, maximum entropy and truncation error. We consider the conjecture that the optimal layout is a Hamiltonian path, and that it optimizes a simply computable geometric cost function. Finding the minimum of such a function, which is a variant of the minimum linear arrangement problem, provides the DMRG with an efficient layout of the lattice and improves both accuracy and convergence time. We present applications to the antiferromagnetic and spin glass spin-1/2 models on the square and triangular lattices.

We consider the two-band Hubbard model, where electrons from different bands interact through an on-site one- and two-particle hybridization. The proposed Hamiltonian makes it possible to construct an effective theory and answer the question of the nature of pairing: conduction electrons form pairs due to two-particle hybridization of electrons from different bands, compensating for the direct Hubbard repulsion between conduction electrons. It is shown that an effective attraction between conduction electrons leads to $η$-pairing. The electron-electron pairing mechanism explains the presence of superconductivity at high temperatures experimentally observed in hydrogen-rich materials at high pressure.

Theoretical spectroscopy, and more generally, electronic-structure theory, are powerful concepts for describing the complex many-body interactions in materials. They comprise a variety of methods that can capture all aspects, from ground-state properties to lattice excitations to different types of light-matter interaction, including time-resolved variants. Modern electronic-structure codes implement either a few or several of these methods. Among them, exciting is an all-electron full-potential package that has a very rich portfolio of all levels of theory, with a particular focus on excitations. It implements the linearized augmented planewave plus local orbital basis, which is known as the gold standard for solving the Kohn-Sham equations of density-functional theory. Based on this, it also offers benchmark-quality results for a wide range of excited-state methods. In this review, we provide a comprehensive overview of the features implemented in exciting in recent years, accompanied by short summaries on the state of the art of the underlying methodologies. They comprise density-functional theory and time-dependent density-functional theory, density-functional perturbation theory for phonons and electron-phonon coupling, many-body perturbation theory in terms of the $GW$ approach and the Bethe-Salpeter equation. Moreover, we capture resonant inelastic x-ray scattering, pump-probe spectroscopy as well as exciton-phonon coupling. Finally, we cover workflows and a view on data and machine learning. All aspects are demonstrated with examples for scientific relevant materials.

We extend a recently developed "tangent fermion" method to discretize the Hamiltonian of a helical Luttinger liquid on a one-dimensional lattice, including two-particle backscattering processes that may open a gap in the spectrum. The fermion-doubling obstruction of the sine dispersion is avoided by working with a tangent dispersion, preserving the time-reversal symmetry of the Hamiltonian. The numerical results from a tensor network calculation on a finite lattice confirm the expectation from infinite-system analytics, that a gapped phase with spontaneously broken time-reversal symmetry emerges when the Fermi level is tuned to the Dirac point and the Luttinger parameter crosses a critical value.

Strongly-correlated fermion systems on a lattice have been a subject of intense focus in the field of condensed-matter physics. These systems are notoriously difficult to solve, even with state-of-the-art numerical methods, especially in regimes of parameters where degrees of freedom compete or cooperate at similar energy and length scales. Here, we introduce a spin-1 slave-particle technique to approximately treat the t-U-V-J fermionic model at arbitrary electron dopings in an economical manner. This formalism respectively maps the original charge and spin degrees of freedom into effective pseudo-spin and pseudo-fermion sectors, which are treated using a self-consistent cluster mean-field method. We study the phase diagram of the model under various conditions and report the appearance of charge and spin stripes within this formalism. These stripes are a consequence of the cluster mean-field treatment of the pseudo-particle sectors and have not been detected in previous slave-particle studies. The results obtained agree qualitatively well with what more reliable numerical methods capture.

We combine $GW$ plus extended dynamical mean field theory ($GW$+EDMFT) with the density functional theory for superconductors (SCDFT) framework to study the electronic properties of LiTi$_2$O$_4$. Excellent agreement with experiment is obtained for the density of states, mass enhancement, Sommerfeld coefficient and superconducting $T_c$, if the dynamical nature of the screened Coulomb interaction is taken into account. Our results show that the coupling to collective charge fluctuations (plasmons) plays an important role in the pairing mechanism and explains the remarkably high $T_c$ of this moderately correlated spinel compound.

We study the Lanczos coefficients in a quadratic model given by an impurity interacting with a multi-mode field of fermions, also known as single impurity Anderson model. We analytically derive closed expressions for the Lanczos coefficients of Majorana fermion operators of the impurity for different structures of the coupling to the hybridisation band at zero temperature. While the model remains quadratic, we find that the growth of the Lanczos coefficients structurally depends strongly on the chosen coupling. Concretely, we find $(i)$ approximately constant, $(ii)$ exactly constant, $(iii)$ square root-like as well $(iv)$ linear growth in the same model. We further argue that in fact through suitably chosen couplings, essentially arbitrary Lanczos coefficients can be obtained in this model. These altogether evince the inadequacy of the Lanczos coefficients as a reliable criterion for classifying the integrability or chaoticity of the systems. Eventually, in the wide-band limit, we find exponential decay of autocorrelation functions in all the settings $(i)-(iv)$, which demonstrates the different structures of the Lanczos coefficients not being indicative of different physical behavior.

We perform a systematic study on the effects of electric fields in the Luttinger liquid phase of the one-dimensional spin-$1/2$ $K_1J_1Γ_1Γ_1^\prime K_2J_2$ model in the region of ferromagnetic nearest-neighboring Kitaev coupling. We find that while electric fields along $(1,1,1)$-direction maintain the Luttinger liquid behavior, fields along other directions drive the system to a dimerized state. An estimation is made on how effective a $(1,1,1)$-field is for tuning the Luttinger parameter in real materials. Our work is useful for understanding the effects of electric fields in one-dimensional generalized Kitaev spin models, and provides a starting point for exploring the electric-field-related physics in two dimensions based on a quasi-one-dimensional approach.

Utilizing a combination of dynamical mean field theory and density functional theory (DMFT/DFT), it has been theoretically proposed that PuB$_6$ is a strongly correlated topological insulator characterized by nontrivial $\mathbf{Z}_{2}$ topological invariants and metallic surface states (\textit{X. Deng et al., Phys. Rev. Lett. 111, 176404 (2013)}). Here, we demonstrate through low-temperature magneto-transport measurements and first-principles calculations that PuB$6$ exhibits characteristics of a topological Kondo insulating state. These features include a transition in electrical resistivity from high-temperature, thermally activated behavior with a narrow gap at the Fermi level ($Δρ \sim$ 20 meV) to a distinctive low-temperature plateau, as well as a surface-to-volume dependence of electrical resistivity at low temperatures. The topological nature of PuB$_6$ is further supported by the theoretical calculations, which show that GGA+$U$ is capable of capturing electronic, topological, and lattice properties of PuB$_6$ with much lower computational cost than DMFT.

We use determinant quantum Monte Carlo to study the half-filled `bond-Holstein' model on a square lattice. We find that the model exhibits a charge-density-wave (CDW) phase transition with a critical temperature $T_\mathrm{cdw}$ considerably higher than that of the canonical `site-Holstein' model. Using a finite-size scaling analysis of the charge structure factor $S_{\rm cdw}$, we obtain $T_\mathrm{cdw}$ to greater than one percent accuracy. At the same time, local observables also show clear signatures consistent with the transition temperatures inferred from our scaling analysis. We attribute the enhanced CDW tendencies to a phonon-mediated nearest-neighbor electron repulsion that is directly proportional to the dimensionless electron-phonon coupling $λ$ in the atomic ($t\rightarrow 0$) limit. This behavior contrasts with the site-Holstein case, where the same limit yields only an on-site attraction. We supplement our analysis with results from several unsupervised machine learning methods, which not only confirm our estimates of $T_\mathrm{cdw}$ but also provide insight into the high-temperature crossover between a metallic and bipolaron liquid regime.

We investigate the role of strong Coulomb interactions beyond the standard Hubbard model in two distinct physical contexts. First, we analyze the superconducting phase transition occurring near the Mott metal-insulator transition. Second, we study transport properties of artificial nano-scale structures containing quantum dots coupled to external electrodes. In both cases, we focus on the impact of the correlated (assisted) hopping (CH) interaction. For superconductors, CH acts as a driving mechanism for the phase transition and modifies the spectral properties of the system. We present the evolution of the spectral function as the system approaches the Mott-type transition under varying model parameters. In quantum-dot-based devices, CH influences the tunneling amplitude between the dot and metallic leads. We demonstrate that the characteristic changes in the conductance of a normal metal-quantum dot-normal metal structure provide a clear signature of the presence and sign of CH interaction.

Large-scale tensor network simulations are crucial for developing robust complexity-theoretic bounds on classical quantum simulation, enabling circuit cutting approaches, and optimizing circuit compilation, all of which aid efficient quantum computation on limited quantum resources. Modern exascale high-performance computing platforms offer significant potential for advancing tensor network quantum circuit simulation capabilities. We implement TN-Sim, a tensor network simulator backend within the NWQ-Sim software package that utilizes the Tensor Algebra for Many-body Methods (TAMM) framework to support both distributed HPC-scale computations and local simulations with ITensor. To optimize the scale up in computation across multiple nodes we implement a task based parallelization scheme to demonstrate parallelized gate contraction for wide quantum circuits with many gates per layer. Through the integration of the TAMM framework with Matrix Product State (MPS) tensor network approaches, we deliver a simulation environment that can scale from local systems to HPC clusters. We demonstrate an MPS tensor network simulator running on the state-of-the-art Perlmutter (NVIDIA) supercomputer and discuss the potential portability of this software to HPC clusters such as Frontier (AMD) and Aurora (Intel). We also discuss future improvements including support for different tensor network topologies and enhanced computational efficiency.

We construct asymptotic quantum many-body scars (AQMBS) in one-dimensional SU($N$) Hubbard chains ($N\geq 3$) by embedding the scar subspace into an auxiliary Hilbert subspace $\mathcal{H}_P$ and identifying a parent Hamiltonian within it, together with a corresponding extension of the restricted spectrum-generating algebra to the multi-ladder case. Unlike previous applications of the parent-Hamiltonian scheme, we show that the parent Hamiltonian becomes the SU($N$) ferromagnetic Heisenberg model rather than the spin-1/2 case, so that its gapless magnons realize explicit AQMBS of the original model. Working in the doublon-holon subspace, we derive this mapping, obtain the one-magnon dispersion for periodic and open boundaries, and prove (i) orthogonality to the tower of scar states, (ii) vanishing energy variance in the thermodynamic limit, and (iii) subvolume entanglement entropy with rigorous MPS/MPO bounds. Our results broaden the parent-Hamiltonian family for AQMBS beyond spin-1/2 and provide analytic, low-entanglement excitations in SU($N$)-symmetric systems.

We have studied the magnetic properties of clinoatacamite Cu$_2$Cl(OH)$_3$, the parent compound of the quantum spin liquid candidate herbertsmithite and a longstanding puzzle among frustrated quantum magnets. As we reveal using density-functional theory, clinoatacamite belongs to the class of distorted kagome antiferromagnets with the kagome plane being embedded into a low-symmetry crystal structure. By means of thermodynamic measurements, muon spin rotation/relaxation as well as neutron diffraction on single crystals, we find a complex sequence of phases/regions below 18.1 K in zero magnetic field. We propose this complexity in multicritical clinoatacamite to arise from the competition of antiferromagnetic ordering modes from the underconstrained manifold of modes, which can lead to a metamagnetic texture in zero magnetic field.

Can superconducting order be engineered as a robust attractor of open-system dynamics in strongly correlated systems? We demonstrate this possibility by proposing a minimal dissipation-engineering protocol for the particle-hole symmetric Hubbard model. By applying a rotated quantum jump operator, specifically a locally transformed $η$-pair lowering operator, on a single lattice site only, we show that the Lindblad evolution autonomously pumps the system from the vacuum into a nonequilibrium steady state (NESS) exhibiting macroscopic $η$-pair off-diagonal long-range order (ODLRO). Crucially, this local-to-global synchronization stands in stark contrast to schemes reliant on spatially extensive reservoirs: here, a single local dissipative seed suffices to establish long-range coherence across the entire interacting lattice system. We elucidate the underlying mechanism via three core features: local dark-state selection, the controlled elimination of off-manifold excursions induced by hopping, and a Liouvillian invariant-subspace structure that yields an attractive fixed point with a finite dissipative gap. Furthermore, we systematically classify the stability of this NESS with respect to static disorder, and identify a wide regime in which the superconducting attractor remains robust against Hamiltonian perturbations that preserve the effective subspace structure. We also pinpoint specific perturbations that directly cause dephasing of the $% η$-pseudospin coherence and suppress ODLRO. These findings open up a disorder-tolerant pathway for stabilizing superconducting order as a non-thermal attractor through minimal local quantum-jump control.

We investigate the ground-state phase diagram of the $S$=1/2 XXZ model on the two-dimensional Shastry-Sutherland lattice using exact diagonalization (ED), density-matrix renormalization group (DMRG), and cluster mean-field theory (CMFT) with DMRG as a solver. In the isotropic case ($Δ=1$), CMFT results reveal an intermediate empty plaquette (EP) phase that has a lower energy than the full plaquette (FP) phase. However, due to mean-field artifacts, CMFT alone is not suitable for accurately determining phase boundaries. Therefore, we combined three methods to map out the reliable phase diagram. Our calculations show that the EP phase narrows as $Δ$ deviates from unity and eventually vanishes. More importantly, we identify a spiral phase at small $Δ$, which has not been reported in previous studies. This phase is clearly captured by DMRG simulations on long cylindrical geometries. The competition between the EP, spiral, and $xy$-AFM phases near their boundaries provides a plausible explanation for the emergent spin-liquid-like behavior in RE$_2$Be$_2$GeO$_2$, while shedding new light on the role of XXZ anisotropy in the Shastry-Sutherland XXZ model.

Can superconducting order be made a robust attractor of open-system dynamics in strongly correlated lattices? We demonstrate that it can by proposing a minimal dissipation-engineering protocol for the particle--hole symmetric Hubbard model. By applying a rotated quantum jump operator, a locally transformed $η$-pair lowering operator, on as little as a single lattice site, we show that the Lindblad evolution autonomously pumps the system from the vacuum into a nonequilibrium steady state (NESS) with macroscopic $η$-pair off-diagonal long-range order (ODLRO). Crucially, this local-to-global synchronization contrasts with schemes requiring spatially extensive reservoirs: here, a strictly local dissipative seed suffices to establish coherence across the interacting system. We elucidate the mechanism via local dark-state selection, controlled elimination of off-manifold excursions induced by hopping, and a Liouvillian invariant-subspace structure that yields an attractive fixed point with a finite dissipative gap. Furthermore, we classify the stability of this NESS against static disorder, identifying a broad regime where the superconducting attractor is resilient to Hamiltonian perturbations that leave the effective subspace structure intact, while pinpointing specific perturbations that directly dephase the $η$-pseudospin coherence and suppress ODLRO. Our results establish a disorder-tolerant route to stabilizing superconducting order as a non-thermal attractor via minimal local quantum-jump control.

Although electrostatics can be incorporated into machine-learned interatomic potentials, existing approaches are computationally very demanding, limiting large-scale, long-time simulations of electrostatics-driven phenomena such as dielectric response, infrared activity, and field-matter coupling. Here, we extend the neuroevolution potential (NEP), a highly efficient machine-learned interatomic potential, to a charge-aware framework (qNEP) by introducing explicit, environment-dependent partial charges. Each ionic partial charge is represented by a neural network as a function of the local descriptor vector, analogous to the NEP site-energy model. This formulation enables the direct prediction of the Born effective charge tensor for each ion and, consequently, the polarization. As a result, dielectric properties, infrared spectra, and coupling to external electric fields can be evaluated within a unified framework. We derive consistent expressions for the forces and virials that explicitly account for the position dependence of the partial charges. The qNEP method has been implemented in the free-and-open-source GPUMD package, with support for both Ewald summation and particle-particle particle-mesh treatments of electrostatics. We demonstrate the accuracy and efficiency of the qNEP approach through representative applications to water, Li7La3Zr2O12, BaTiO3, and a magnesium-water interface. These results show that qNEP enables accurate atomistic simulations with explicit long-range electrostatics, scalable to million-atom systems on nanosecond time scales using consumer-grade GPUs.

We investigate operator delocalization in disordered one-dimensional spin chains by introducing -- besides the already known operator mass -- a complementary measure of operator complexity: the operator length. Like the operator nonstabilizerness, both these quantities are defined from the expansion of time-evolved operators in the Pauli basis. They characterize, respectively, the number of sites on which an operator acts nontrivially and the spatial extent of its support. We show that both the operator mass and length can be computed efficiently and exactly within a matrix-product-operator (MPS) framework, providing direct access to their full probability distributions, without resorting to stochastic sampling. Applying this approach to the disordered XXZ spin-1/2 chain, we find sharply distinct behaviors in non-interacting and interacting regimes. In the Anderson-localized case, operator mass, length, and operator entanglement entropy rapidly saturate, signaling the absence of scrambling. By contrast, in the many-body localized (MBL) regime, for arbitrarily weak interactions, all quantities exhibit a robust logarithmic growth in time, consistent with the known logarithmic light cone of quantum-correlation propagation in MBL. We demonstrate that this behavior is quantitatively captured by an effective $\ell$-bit model and persists across system sizes accessible via tensor-network simulations.

Integer and fractional Chern insulators exhibit a nonzero quantized anomalous Hall conductivity due to a spontaneous breaking of time reversal symmetry. To identify nontrivial topology in their time-reversal symmetric many-body spectra, we introduce many-body Euler numbers as a counterpart to many-body Chern numbers. Exemplarily, we perform calculations in a topological Hubbard model that can realize Chern and fractional Chern insulating phases. Furthermore, we lay out a classification scheme to realize different topological phases in interacting systems using symmetry indicators in analogy to topological band theory.

Recent advances in X-ray free-electron laser diagnostics have enabled direct probing of the electronic structure under extreme pressures and temperatures, such as those encountered in stellar interiors and inertial confinement fusion experiments, challenging theoretical models for interpreting experimental data. Kohn-Sham density functional theory (KSDFT) has been successfully applied to analyze experimental X-ray scattering measurements, but its high computational cost renders routine application impractical. Orbital-free DFT (OFDFT) is a substantially more efficient alternative, with computational cost scaling linearly with system size and a weak temperature dependence, yet it often lacks the accuracy required for electronic structure description. Overcoming this limitation, we present a non-empirical Kohn-Sham-assisted orbital-free density functional framework for calculations at extreme conditions, which enables efficient OFDFT simulations with KSDFT-level accuracy for electron densities, electron-ion structure factors, and equations of state across a broad range of conditions. Benchmark comparisons with quantum Monte Carlo data for dense hydrogen and validation against Rayleigh weight measurements of hot dense beryllium demonstrate the reliability of the framework and speedups of up to several hundred times compared with KSDFT. We further show that even at temperatures on the order of 100 eV, quantum non-locality remains essential for correctly describing the electronic structure of dense hydrogen.

The interplay among electronic correlation, topology, and time-reversal-symmetry (TRS) often leads to exotic quantum states of matter. Primary examples include the recently realized fractional Chern insulators (FCIs) in twisted MoTe2 bilayers (tMoTe2) and multilayer graphene aligned with hBN, where TRS is broken in partially filled flat moire Chern bands. Among the FCIs in tMoTe2, the most robust is at a hole filling of v = -2/3 per moire unit cell. Interestingly, transient optical sensing and more recent transport measurements revealed a correlated state at v = -4/3, twice the filling factor for the v = -2/3 FCI. Here, employing pump-probe circular dichroism (CD) measurements on tMoTe2 with twist angles = 3.9 degree and 3.7 degree, we find that the v = -4/3 state exhibits vanishing magnetization (m = 0) in finite windows of out-of-plane magnetic field less than ~2-4 mT, and a first order phase transition to + - m states at higher fields. This out-of-plane antiferromagnetic (AFM) like behavior is notably absent for all other correlated states and disappears for the v = -4/3 state at higher or lower twist angles = 4.0 degree and 3.3 degree. The observed magnetic signature at v = -4/3 is consistent with a predicted fractional topological insulator (FTI) with TRS, consisting of two copies of -2/3 FCIs with opposite chiralities. We support these findings with calculations in the interacting continuum model of tMoTe2. Our work presents a candidate for fractional topological insulators with TRS.

Polaron formation in pump-probe experiments is an inherently non-equilibrium phenomenon, driven by the ultrafast coupled dynamics of electrons and phonons, and culminating in the emergence of a localized quasiparticle state. In this work, we present a first-principles quantum-kinetic theory of polaron formation that captures the real-time evolution of electronic and lattice degrees of freedom in presence of electron-phonon coupling. We implement this framework to investigate the ultrafast polaron formation in the prototypical polar insulator MgO. This approach allows us to determine the characteristic timescales of polaron localization and to identify its distinctive dynamical fingerprint. Our results establish clear and experimentally accessible criteria for identifying polaron formation in pump-probe experiments.

Angle-resolved photoemission spectroscopy (ARPES) provides a direct access to the electronic band structure of solid and molecular systems. The momentum range accessible by this technique depends directly on the photon energy used, and low-photon-energy sources are insufficient to photoemit electrons over the full Brillouin zone of most quantum materials. In addition, while electrons are emitted over a 2$π$ solid angle, conventional hemispherical analyzers only collect a small subset of those electrons. A previous work [RSI 92, 123907 (2021)] demonstrated that electrons emitted over a larger field-of-view can be acquired in one fixed configuration by accelerating them towards the analyzer with a bias voltage. Here, we extend this work by leveraging the deflector technology of novel ARPES hemispherical analyzers. We demonstrate the ability to detect all $2π$ photoemitted electrons in a fixed configuration for various materials such as gold, cuprates and transition-metal dichalcogenides. This approach is especially advantageous for time-resolved ARPES, as electron dynamics over a large momentum range can be accessed with identical measurement conditions.

Coherent interconnects between gate-defined silicon quantum processing units are essential for scalable quantum computation and long-range entanglement. We argue that one-dimensional electron channels formed in the silicon quantum well of a Si/SiGe heterostructure exhibit strong Coulomb interactions and realize strongly interacting Luttinger liquid physics. At low electron densities, the system enters a Wigner regime characterized by dominant 4kF correlations; increasing the electron density leads to a crossover from the Wigner regime to a Friedel regime with dominant 2kF correlations. We support these results through large-scale density matrix renormalization group (DMRG) simulations of the interacting ground state under both screened and unscreened Coulomb potentials. We propose experimental signatures of the Wigner-Friedel crossover via charge transport and charge sensing in both zero- and high-magnetic field limits. We also analyze the impact of short-range correlated disorder - including random alloy fluctuations and valley splitting variations - and identify that the Wigner-Friedel crossover remains robust until disorder levels of about 400 micro eV. Finally, we show that the Wigner regime enables long-range capacitive coupling between quantum dots across the interconnect, suggesting a route to create long-range entanglement between solid-state qubits. Our results position silicon interconnects as a platform for studying Luttinger liquid physics and for enabling architectures supporting nonlocal quantum error correction and quantum simulation.

Over-parameterization is commonly used to increase the expressivity of variational quantum circuits (VQCs), yet deeper and more highly parameterized circuits often exhibit poor trainability and limited generalization. In this work, we provide a theoretical explanation for this phenomenon from a function-class perspective. We show that sufficiently expressive, unstructured variational ansatze enter a Haar-like universality class in which both observable expectation values and parameter gradients concentrate exponentially with system size. As a consequence, the hypothesis class induced by such circuits collapses with high probability to a narrow family of near-constant functions, a phenomenon we term simplicity bias, with barren plateaus arising as a consequence rather than the root cause. Using tools from random matrix theory and concentration of measure, we rigorously characterize this universality class and establish uniform hypothesis-class collapse over finite datasets. We further show that this collapse is not unavoidable: tensor-structured VQCs, including tensor-network-based and tensor-hypernetwork parameterizations, lie outside the Haar-like universality class. By restricting the accessible unitary ensemble through bounded tensor rank or bond dimension, these architectures prevent concentration of measure, preserve output variability for local observables, and retain non-degenerate gradient signals even in over-parameterized regimes. Together, our results unify barren plateaus, expressivity limits, and generalization collapse under a single structural mechanism rooted in random-matrix universality, highlighting the central role of architectural inductive bias in variational quantum algorithms.

A state-of-the-art Resonant Inelastic X-ray Scattering (RIXS) facility, composed of a dedicated beamline and a 2D-RIXS spectrometer, has been constructed and commissioned at BL02U in NanoTerasu, Japan. This paper reports the optical design and optimization of this spectrometer, aiming for an ultrahigh energy resolution of <10 meV in the soft X-ray range. Conventional RIXS spectrometers using monochromatic incident X-rays suffer from a severe trade-off between energy resolution and measurement efficiency, which makes achieving resolutions of <10 meV practically unfeasible. To overcome this limitation, we adopted the $hν^2$ concept using dispersive incident X-rays, based on a comparative study. This approach decouples the energy resolution from the incident bandwidth. We estimate that our 2D-RIXS spectrometer improves efficiency by more than a factor of 10 compared to conventional spectrometers. The 2D-RIXS spectrometer has been successfully commissioned, demonstrating that sub-10 meV resolution measurement is achievable with practical throughput. While development continues towards the ultimate design resolution, the successful implementation validates the design strategy and offers a pathway for future high-performance RIXS instrumentation.

The study of quantum and classical correlations between subsystems is fundamental to understanding many-body physics. In quantum information theory, the quantum mutual information, $I(A;B)$, is a measure of correlation between the subsystems $A,B$ in a quantum state, and is defined by the means of the von Neumann entropy: $I\left(A;B\right)=S\left(ρ_{A}\right)+S\left(ρ_{B}\right)-S\left(ρ_{AB}\right)$. However, such a computation requires an exponential amount of resources. This is a defining feature of quantum systems, the infamous ``curse of dimensionality'' . Other measures, which are based on Rényi divergences instead of von Neumann entropy, were suggested as alternatives in a recent paper showing them to possess important theoretical features, and making them leading candidates as mutual information measures. In this work, we concentrate on the maximal Rényi divergence. This measure can be shown to be the solution of a generalized eigenvalue problem. To calculate it efficiently for a 1D state represented as a matrix product state, we develop a generalized eigenvalue version of the density matrix renormalization group algorithm. We benchmark our method for the paradigmatic XXZ chain, and show that the maximal Rényi divergence may exhibit different trends than the von Neumann mutual information.

Developing new metal hydrides is a critical step toward efficient hydrogen storage in carbon-neutral energy systems. However, existing materials databases, such as the Materials Project, contain a limited number of well-characterized hydrides, which constrains the discovery of optimal candidates. This work presents a framework that integrates causal discovery with a lightweight generative machine learning model to generate novel metal hydride candidates that may not exist in current databases. Using a dataset of 450 samples (270 training, 90 validation, and 90 testing), the model generates 1,000 candidates. After ranking and filtering, six previously unreported chemical formulas and crystal structures are identified, four of which are validated by density functional theory simulations and show strong potential for future experimental investigation. Overall, the proposed framework provides a scalable and time-efficient approach for expanding hydrogen storage datasets and accelerating materials discovery.

Accurate prediction of dielectric tensors is essential for accelerating the discovery of next-generation inorganic dielectric materials. Existing machine learning approaches, such as equivariant graph neural networks, typically rely on specially-designed network architectures to enforce O(3) equivariance. However, to preserve equivariance, these specially-designed models restrict the update of equivariant features during message passing to linear transformations or gated equivariant nonlinearities. The inability to implicitly characterize more complex nonlinear structures may reduce the predictive accuracy of the model. In this study, we introduce a frame-averaging-based approach to achieve equivariant dielectric tensor prediction. We propose GoeCTP, an O(3)-equivariant framework that predicts dielectric tensors without imposing any structural restrictions on the backbone network. We benchmark its performance against several state-of-the-art models and further employ it for large-scale virtual screening of thermodynamically stable materials from the Materials Project database. GoeCTP successfully identifies various promising candidates, such as Zr(InBr$_3$)$_2$ (band gap $E_g = 2.41$ eV, dielectric constant $\overline{\varepsilon} = 194.72$) and SeI$_2$ (anisotropy ratio $α_r = 96.763$), demonstrating its accuracy and efficiency in accelerating the discovery of advanced inorganic dielectric materials.

We introduce a family of generalized Rokhsar-Kivelson (RK) Hamiltonians, which are reverse-engineered to have an arbitrary edge-weighted superposition of dimer coverings as their exact ground state at the RK point. We then focus on a quantum dimer model on the triangular lattice, with doubly-periodic edge weights. For simplicity we consider a $2\times1$ periodic model in which all weights are set to one except for a tunable horizontal edge weight labeled $α$. We analytically show that the model exhibits a continuous quantum phase transition at $α=3$, changing from a topological $\mathbb{Z}_2$ quantum spin liquid ($α<3$) to a columnar ordered state ($α>3$). The dimer-dimer correlator decays exponentially on both sides of the transition with the correlation length $ξ\propto1/|α-3|$ and as a power-law at criticality. The vison correlator exhibits an exponential decay in the spin liquid phase, but becomes a constant in the ordered phase. We explain the constant vison correlator in terms of loops statistics of the double-dimer model. Using finite-size scaling of the vison correlator, we extract critical exponents consistent with the 2D Ising universality class.

We investigate the ground-state and finite-temperature properties of the $J_1$-$J_2$ Heisenberg antiferromagnet on the square lattice in the presence of an external magnetic field. We focus on the highly frustrated regime around $J_2 \approx J_1/2$. The $h$-$T$ phase diagram is investigated with particular emphasis on the finite-temperature transition into the "up-up-up-down" state that is stabilized by thermal and quantum fluctuations and manifests itself as a plateau at one half of the saturation magnetization in the quantum case. We also discuss the enhanced magnetocaloric effect associated to the ground-state degeneracy that arises at the saturation field for $J_2=J_1/2$. For reference, we first study the classical case by classical Monte Carlo simulations. Then we turn to the extreme quantum limit of spin-1/2 where we perform zero- and finite-temperature Lanczos calculations.