arXiv Curator

In unconventional superconductivity, the pairing mechanism remains one of the central unresolved problems in solid-state physics. Moreover, the normal state shows strong deviations from conventional Fermi-liquid behavior, most notably a scattering rate that is linear in energy and temperature. Here we study the momentum dependence of the energy and damping of collective excitations in optimally doped cuprates in the normal state. In particular, we analyze existing Electron Energy Loss Spectroscopy (EELS) data of Bi2Sr2CaCu2O8 [Phys. Rev. B 44, 7155 (1991)] and new analogous optical plasmon data of Bi2Sr2Ca2Cu3O10. Furthermore, we reanalyze similar Resonant Inelastic Xray Scattering (RIXS) acoustic plasmon data [Phys. Rev. R 6, 043184 (2024)] of La1.84Sr0.16CuO4 and La1.84Ce0.16CuO4. We contrast these two-particle excitation results with one-particle excitations measured by Angle-Resolved Photoemission Spectroscopy (ARPES). Most of the plasmon behavior can be described well by combined self-energies of a hole and an electron. The dispersion of optical plasmons at high energy is only weakly determined by correlation effects, indicating resilient quasiparticles. The damping of the optical plasmons at low energy, however, is strongly influenced by the reduced scattering rates described in the marginal Fermi liquid model for "strange-metal" behavior. In addition, the same model describes the linear-in-energy damping of acoustic plasmons. Contrary to some theoretical predictions, we see no overdamped optical plasmons below the classical Lindhard critical wave vector.

Superconductivity in the Ruddlesden-Popper (RP) nickelates has so far been confined either to multilayer compounds under high pressure or to square-planar infinite-layer-type films stabilized by epitaxial strain. Here we report the high-pressure high-temperature (HPHT) growth of bulk single crystals of the quadruple-layer RP nickelate La$_5$Ni$_4$O$_{13}$ ($n=4$) and reproduce the quintuple-layer La$_6$Ni$_5$O$_{16}$ ($n=5$) synthesis, two members of the series that had previously been accessible only as thin films. Topotactic reduction of La$_5$Ni$_4$O$_{13}$ single crystals with CaH$_2$ yields the square-planar phase La$_5$Ni$_4$O$_{10}$, which we show to be a bulk, ambient-pressure superconductor with $T_\mathrm{c}^\mathrm{onset} \approx 18\,\mathrm{K}$, a shielding fraction above 60\%, and a specific-heat jump that confirms a bulk pairing transition. All crystals show besides superconductivity a coexistence with a density-wave or metal-insulator transition. The upper critical field shows a small anisotropy, $γ\approx 2.4$, which remains finite as $T \to T_c$, indicating three-dimensional superconductivity with substantial coupling between the four-layer blocks. These results identify the square-planar nickelates $R_{n+1}$Ni$_n$O$_{2n+2}$ as a bulk platform in which the layer number alone tunes the nickel valence into the superconducting regime, without chemical substitution, pressure, or epitaxial constraint.

The predominantly $d_{x^2-y^2}$-wave pairing symmetry of cuprate superconductors is well established, yet its microscopic manifestation in local Cooper-pair tunneling remains less explored. Here, we investigate Bi$_{1.2}$Pb$_{0.8}$Sr$_2$CaCu$_2$O$_{8+δ}$, using Josephson scanning tunneling spectroscopy with a conventional $s$-wave superconducting Pb tip. We observe a zero-bias conductance enhancement whose approximately quadratic dependence on junction conductance is consistent with phase-diffusive Josephson tunneling. Atomic-resolution conductance maps and site-resolved spectra reveal a pronounced intra-unit-cell modulation, with enhanced Josephson response near oxygen sites and strong suppression near copper sites and along diagonal directions. Microscopic calculations incorporating the bilayer electronic structure and realistic Wannier-based tunneling matrix elements reproduce the observed spatial pattern through interference between pairing contributions of opposite signs. Our results reveal how the sign-changing $d$-wave superconducting order parameter manifests itself in atomic-scale Josephson tunneling, providing complementary evidence for unconventional pairing symmetry in cuprates.

Superconductivity arises from Cooper pairing, but conventional photoemission probes electrons individually rather than as pairs. Here we report the observation of Cooper-pair emission from an optimally doped cuprate superconductor using two-electron angle-resolved photoemission (2e-ARPES). Below the superconducting transition, we observe a distinct two-electron emission peak at zero pair binding energy and zero pair momentum that disappears in the normal state, revealing the characteristic fingerprint of a Cooper pair. Unexpectedly, the joint energy distribution broadens markedly along the relative-energy axis as temperature decreases, while the total energy width remains nearly unchanged. This anomalous behavior cannot be explained by instrumental resolution, static disorder or thermal broadening, and instead points to scattering processes involving collective excitations of the superconductor. Our results establish 2e-ARPES as a probe of both Cooper-pair emission and the many-body dynamics of the superconducting state.

We extend the scalability of variational infinite two-dimensional tensor network methods with application to investigate frustrated magnetism of the square-lattice $J_1$-$J_2$ model and kagome-lattice Heisenberg model. In the sign-problem-free limit of the square-lattice model, we achieve high ground-state energy precision comparable to quantum Monte Carlo methods. With no obstruction beyond this limit, our method obtains direct estimates of order parameters. For the square-lattice $J_1$-$J_2$ model, we determine isotropic plaquette order in the valence-bond solid region. For the model near $J_2/J_1=0.5$, and the kagome lattice Heisenberg model, our data indicate linear-dispersion gapless states. To probe the candidate spin-liquid nature of the highly frustrated cases, we go beyond inferring the absence of symmetry-breaking orders and analyze susceptibilities that characterize gapless spin liquids. By estimating scaling dimensions, our data indicate a stable gapless spin liquid on the kagome lattice, and competing scenarios for the intermediate regime of the square-lattice model. Our developments include parallelization and randomized low-rank methods and largely expand the frontier of high-precision variational numerics for quantum phases of matter.

Whether $d$-wave or $s_\pm$-wave pairing prevails in the bilayer nickelate \ce{La3Ni2O7} remains unresolved, and how its observed spin order and magnetic excitations relate to the pairing mechanism is unclear. Here, we use finite-temperature determinant quantum Monte Carlo simulations of a minimal bilayer Hubbard model to investigate pairing competition at strong interlayer coupling, where intralayer and interlayer antiferromagnetic correlations have comparable strengths. We find that spin stripes develop spontaneously with increasing hole doping, accompanied by a crossover of the dominant pairing symmetry from $d$-wave to $s_\pm$-wave. The crossover tracks the changing relative strengths of the spin correlations, with the interlayer correlations overtaking their intralayer counterparts in magnitude. Concomitantly, the dynamical spin structure factor exhibits a momentum-dependent reorganization of the magnetic excitations, consistent with fluctuating spin stripes. Our results identify the structure of spin fluctuations as a key link between spin stripes and pairing competition, offering a minimal-model perspective for understanding \ce{La3Ni2O7}.

Rare-earth nickelates are among the most active oxide electrocatalysts for the oxygen evolution reaction (OER), owing to their tunable Ni-O covalency, near-unity eg occupancy, and sensitivity to A-site chemistry, oxygen stoichiometry, and epitaxial strain. Here, we use epitaxial A-site high-entropy (La0.2Pr0.2Nd0.2Sm0.2Eu0.2)NiO3 thin films to reveal how epitaxial strain influences defect formation, Ni valence, electronic structure, and OER activity in a chemically disordered rare-earth (RE) sublattice. Tensile strain promotes Ruddlesden-Popper fault formation, interrupts the corner-sharing NiO6 network, and lowers the average Ni oxidation state. Although the occupied density of states near the Fermi level remains largely unchanged, O K-edge X-ray absorption spectroscopy reveals increased low-energy unoccupied states under tensile strain. First-principles calculations suggest that tensile strain shifts RE 4f-derived states toward the Fermi level, enhancing RE 4f-O 2p-Ni 3d hybridization. These coupled structural and electronic changes lead to earlier OER onset, stronger pre-catalytic Ni redox features, and higher activity than compressively strained films and representative single-A-site nickelates. The high-entropy nickelate also shows enhanced resistance to OER-induced structural degradation compared with LaNiO3. These results suggest that A-site configurational disorder and epitaxial strain can jointly tune defect formation and rare-earth-derived electronic states, providing a route for improving perovskite oxide electrocatalysts.

The pseudogap state, an anomalous momentum dependent partial suppression of low energy electronic excitations persisting well above the superconducting critical temperature in underdoped cuprates, remains one of the most profound enigmas in condensed physics today. This review summarizes over three decades of experimental and theoretical progress. It appears that pseudogap is not a mere precursor to superconductivity, but a distinct, intertwined state of matter characterized by a complex interplay of competing and cooperating electronic orders. We critically examine the bulk probes like resistivity and magnetic susceptibility and the surface hallmark spectroscopic signatures observed via angle resolved photoemission spectroscopy (ARPES) in this review. ARPES reveals a Fermi surface truncated into disconnected segments and scanning tunneling microscopy (STM) visualizes a spatially heterogeneous energy gap around the Fermi level. These observations are contextualized alongside transport anomalies, including resistivity and bulk magnetic susceptibility, to construct a phenomenological portrait of a phase that violates conventional Fermi liquid theory. By mapping the electronic phase diagram, we highlight the abrupt transition at above optimal doping. We conclude with open questions still to be answered. This review underscores the role of the pseudogap as a key to unlocking the unconventional pairing mechanism of hole doped cuprates and the broader physics of strongly correlated electron systems.

Superconducting pairing symmetry provides crucial insight into the microscopic mechanisms of Cooper pairing. However, probing the local manifestations of sign-changing superconducting order parameters remains experimentally challenging. Here, we investigate the superconducting condensate in optimally doped (OPT) Bi2Sr2Ca2Cu3O10+δ (Bi-2223) using scanning Josephson tunneling microscopy (SJTM), equipped with an s-wave superconducting tip. Measurements at ultra-low temperatures reveal a pronounced modulation of the Josephson tunneling signal at different locations inside the crystallographic unit cell. A distinct Josephson tunneling signal is detected above O sites in the underlying CuO2 plane, whereas the signal nearly vanishes above the Cu-plaquette diagonals, including Cu sites and plaquette centers. This unique intra-unit-cell spatial pattern agrees with the theoretical prediction for the symmetry-dependent cancellation of Josephson tunneling between an s-wave superconducting tip and a d-wave superconducting sample. Our findings provide atomic-scale, real-space evidence for d-wave pairing symmetry in OPT high-Tc Bi-2223 and establish SJTM as a powerful approach for investigating the pairing symmetry of unconventional superconductors.

Quantum Monte Carlo methods are widely used to solve quantum many-body problems. However, rigorous guarantees for their performance are poorly understood, especially for simulating interacting fermions. In this work, we develop a quantum Monte Carlo algorithm for a broad family of fermionic systems exhibiting Kramers time-reversal symmetry, and prove that their partition functions can be estimated in polynomial time. This family, which contains many important classes of Hubbard-type models, is known to both avoid the sign problem and obey a Lee-Yang zero-freeness theorem. Thus our results suggest that these two properties together may enable efficient classical algorithms for quantum many-body systems more broadly.

We show that the boundary fractionalization familiar from gapped topological insulators can enforce bulk gaplessness in a dipole-conserving system. In a one-dimensional fermion model with charge conservation, dipole conservation, and particle-hole symmetry, we realize an intrinsically gapless symmetry-protected topological (igSPT) phase combining a dipolar Luttinger liquid with half-quantized boundary charges. This boundary structure obstructs a unique gapped bulk preserving these symmetries and lattice translation. Remarkably, the two charge-neutral edge configurations remain exactly degenerate even at finite size due to dipole conservation and particle-hole symmetry. The boundary fractionalization also manifests as shifted conformal towers in the dipole-resolved entanglement spectrum, reflecting distinct conformal boundary conditions of the Luttinger liquid. We support our analysis of the igSPT by studying exactly-solvable limits and performing DMRG computations, establishing its stability for a wide range of tuning parameters.

Quantum impurity calculations become challenging when multiple correlated orbitals are treated explicitly, particularly in nonequilibrium steady states. We introduce Coupled-Cluster Coupled Lindblad (C$^3$L), a deterministic approach that combines the recently developed coupled-Lindblad compact bath representation with coupled-cluster theory to compute steady states and Green's functions directly in imaginary and real frequency. C$^3$L has favorable polynomial computational complexity and avoids the sign problem of quantum Monte Carlo (QMC). In a real-time multiorbital benchmark, it agrees closely with state-of-the-art QMC while reducing computational cost by four to five orders of magnitude, with a further 10--100-fold speedup on a GPU. Beyond the practical reach of established real-time QMC solvers, four-orbital nonequilibrium impurities with general interactions and biased reservoirs become accessible with C$^3$L. We obtain resonance-resolved transport signatures relevant to scanning tunneling microscopy experiments on molecular junctions.

Rigorously establishing spontaneous symmetry breaking (SSB) in the thermodynamic limit is generally challenging. This is exemplified by the two-dimensional (2D) $S = 1/2$ antiferromagnetic Heisenberg model where the existence of Néel order is not rigorously established despite overwhelming evidence from a variety of methods. We show that a bootstrap approach can be used to rigorously establish SSB and use it to prove the existence of Néel order in the 2D $S = 1/2$ XXZ model by deriving a strictly positive lower bound on the Néel order parameter for a set of discrete values of the anisotropy parameter $Δ$ in the range $[0,2]$ including the isotropic (Heisenberg) point $Δ= 1$. Our approach employs a generalization of the Kennedy--Lieb--Shastry inequalities, which use an infrared bound to reduce the existence of long-range order to bounds on certain finite-distance correlators. We compute these bounds for all sufficiently large systems with periodic boundary conditions using the bootstrap. An important ingredient in our proof is a novel method for obtaining exactly certified variational upper bounds on the energy density in the thermodynamic limit from matrix product states. Our approach can be readily generalized to other local lattice models for which infrared susceptibility bounds can be obtained.

We present an implementation of the bare and interacting longitudinal susceptibility in a full-potential linear muffin-tin orbital code. The bare susceptibility is built in real space from Matsubara sums of Kohn-Sham Green's functions and is resolved by site, spin and orbital degrees of freedom; the spin and charge susceptibilities follow from a Dyson equation, with an on-site exchange-correlation kernel. After testing the solver on the Stoner-like instability of bcc Fe, we apply it to four density-wave materials. The bare susceptibility peaks at H in bcc Cr, where the nonmagnetic state is close to a Stoner instability, and at the stripe vector M in FeSe. In LaTe3 the response of the Te nets peaks close to the observed charge-density-wave vector, at its image in the neighbouring Brillouin zone; which of the two equivalent vectors carries the response is set by the phases of the sites in the cell. In 2H-NbSe2 the response is weakly modulated, and the observed vector lies on the flank of a broad maximum, consistent with the known importance of electron-phonon coupling there. Antiferromagnetic order reduces the bare longitudinal response at the ordering vector by ~8% in Cr and about half in FeSe. The orbital decomposition of the bare susceptibility shows which orbitals contribute the most for the nearest-neighbor bonds.

Sample-based quantum diagonalization (SQD) accuracy remains limited for strongly correlated systems like iron-sulfur clusters. This accuracy heavily depends on the mean-field reference state, which determines the sampling ansatz and configuration-recovery target. We extend the qcsc-prefect workflow to support a broken-symmetry unrestricted Hartree-Fock (BS-UHF) starting point, derived from atom-localized antiferromagnetic initial guesses. For the $[Fe_4S_4(SCH_3)_4]^{2-}$ cubane (54 electrons, 36 orbitals, 72 data qubits; 88 total qubits), this lowers the mean-field reference by $\approx 534$ mHa relative to the Restricted Hartree-Fock (RHF) solution (-327.0809 vs. -326.5474 Ha). Combined with deep multi-step configuration recovery on ROQUO GPU nodes at subspace dimension $d = 3 \times 10^{10}$, we reach $E = -327.235076$ Ha. This is 4.3 mHa above the finite-$D = 4000$ DMRG value (-327.2393952 Ha) and $\sim 512-517$ mHa below a previous same-active-space RHF-reference SQD calculation. Extrapolating to the $d \to \infty$ limit via a Gaussian-in-log $d$ fitting ansatz gives an asymptotic estimate $E_\infty = -327.244468$ Ha $\pm 1.05$ mHa, just 0.33 mHa above our discarded-weight-extrapolated DMRG estimate ($E_{D\to\infty} = -327.244797$ Ha $\pm 0.46$ mHa), well within their combined $1σ$ uncertainty. We characterize the determinant-subspace composition driving convergence, show matching validation on $[Fe_2S_2(SCH_3)_4]^{2-}$, and report multi-HPC execution on Fugaku, Miyabi-G, and ROQUO. We conclude that reference-state preparation, rather than subspace size alone, is a key driver of SQD energy accuracy under strong electronic correlation.

We propose a weakly coupled description of the long-range transverse-field Ising chain near its short-range crossover at $s_*=7/4$, where the couplings decay as $r^{-1-s}$. The Ising conformal field theory is coupled to an instantaneous Gaussian field through $σχ$, extending the classical long-range Ising duality to quantum chains. The equal-time interaction generates a stress-tensor deformation that causes the velocity to run and renormalizes the dynamical critical exponent $z$. Conformal perturbation theory gives $z=1-0.720\,δ+O(δ^2)$, $1/ν=(1+z)/2+O(δ^2)$, and $ω=δ+O(δ^2)$, with $δ=s_*-s$. DMRG calculations on periodic chains of up to $128$ sites provide a consistency test using theoretical inputs for $1/ν$ and $ω$. After accounting for slow finite-size corrections, static and dynamical estimates of $z$ agree, supporting the crossover expansion at $s=1.5$ and a two-sided Padé interpolation at $s=1$, $1.25$, and $1.5$.

Many-body physics faces an apparent paradox concerning violations of Luttinger's theorem in strongly correlated systems. The theorem relates the electron number to the Luttinger count -- the monodromy of the logarithm of the imaginary-frequency Green's function $G(iε)$ -- and, at finite magnetic field, connects the Hall conductivity to the magnetic-field derivative of the Luttinger count via the Streda formula. Here, we provide a non-perturbative solution to this paradox. We demonstrate that the reported violations of Luttinger's theorem do not require a breakdown of the skeleton diagram representation of physical quantities, thereby resolving a long-standing point of contention in many-body theory. Furthermore, we show that the Luttinger integral, which accounts for the mismatch between the particle number and the Luttinger count, is itself a genuine topological invariant directly tied to the Luttinger-Ward functional $Φ[G]$ rather than to the Green's function alone. We demonstrate that this invariant, like the monodromy of $\ln G(iε)$ itself, stems exclusively from the zero-frequency discontinuity of the Matsubara trajectory.

The spin-1/2 kagome Heisenberg antiferromagnet is a central model in the study of quantum spin liquids, yet its ground-state properties remain unresolved. We perform large-scale variational Monte-Carlo calculations using neural-network quantum states (NQSs) with kagome-geometry-aware message-passing convolution on two-dimensional tori containing up to 300 spins, and achieve exceptionally low variational energies for the kagome model. The spin and dimer structure factors show no evidence of magnetic or valence-bond long-range order. Moreover, the dominant symmetry-resolved spin and dimer correlations exhibit a common scaling behavior consistent with the emergent $SO(6)$ monopole multiplet of the Dirac spin liquid. Finally, the fidelity of the optimized NQSs with Gutzwiller-projected mean-field states increasingly favors the Dirac state over the competing $Z_{2}$ candidates as the system size grows. Together, these complementary results provide evidence for a $U(1)$ Dirac spin-liquid ground state in the kagome Heisenberg antiferromagnet.

We investigate the interplay between antiferromagnetism and topological superconductivity in twisted cuprate bilayers. We demonstrate how topological $s+d+id$ superconductivity is induced by a twist of $θ\sim 45^\circ$ in a microscopic $t$-$J$ type model of cuprate bilayers, where it robustly coexists with antiferromagnetic order at low hole doping. This novel pairing state becomes physically accessible at realistic interlayer couplings in the overdoped regime. Crucially, it displays a characteristic $2$-fold rotational symmetry in angle-resolved photoemission spectroscopy (ARPES) not observed in the competing time-reversal symmetry breaking $d+id$ or $s+id$ states previously found. A finite Josephson effect at $45^\circ$ twist together with signatures in scanning tunneling spectroscopy of its two Majorana edge modes identifies this topological $s+d+id$ state unambiguously.

We study charge transport in the two-dimensional infinite-$U$ Hubbard model with hopping amplitude $t$ using full exact diagonalization (ED) on $N=16$- and $18$-site square clusters with a single doped hole. We first provide a proof, treating degeneracies explicitly, that a constant dc $f_x$ function, defined as the sum of the squared current matrix elements connecting eigenstates at a given energy to eigenstates of nearly the same energy, averaged over degenerate initial states, is necessary and sufficient for strictly linear-in-$T$ resistivity: $ρ_x(T)=αT$ for all $T>0$ with a constant $α$. The numerically extrapolated dc $f_x$ is not exactly constant, but it is approximately flat over a broad region in the interior of the many-body spectrum and varies strongly within a scale $\varepsilon_{f_x}$ of the spectral edges. This scale $\varepsilon_{f_x}$ is related to a characteristic temperature $T_{\mathrm{cr}}$, which quantifies the leading departure of the resistivity from high-temperature linearity upon cooling. Within the studied temperature range $T/t \geq 0.05$, our ED resistivity results reveal a gradual crossover between high- and low-temperature approximately linear-in-$T$ regimes, with characteristic crossover temperatures $T_{\mathrm{cr}}/t\approx0.414$ and $0.467$ for $N=16$ and $18$, respectively. We also find a broad specific-heat peak associated with this gradual crossover, at $T/t \approx 0.70$ and $0.68$ for $N=16$ and $18$, respectively.

We present a first-principles study of the formal-charge polar discontinuity in superlattices of BaTiO$_3$ interfaced with SmNiO$_3$ or BiNiO$_3$, comparing the results for the insulating and metallic nickelate phases. Polar discontinuities arising from differently charged layers in perovskites are well-studied in heterostructures such as LaAlO$_3$/SrTiO$_3$, and recently, signatures consistent with such a polar discontinuity have been observed in systems where the formally charged layers would reside in a metal. We first identify signatures of the polar discontinuity in the structure and the density of states of superlattices with BaTiO$_3$ and insulating BiNiO$_3$, and find partial screening by modification of the BiNiO$_3$ charge disproportionation pattern analogous to our previous study on BaTiO$_3$/SmNiO$_3$. We then examine superlattices constructed from the corresponding metallic nickelate phases and find that, despite the presence of screening charges, the qualitative interfacial response is remarkably similar to that of the insulating counterparts. We show that the metallic screening charges have little influence on the apparent formal-charge polar discontinuity and discuss practical limitations of the polar-discontinuity picture when interpreting and quantifying screening at metallic nickelate/ferroelectric interfaces.

We study the effect of electron-phonon coupling on spin-fluctuation-mediated pairing in bulk La$_6$Ni$_5$O$_{12}$. A five-layer Hubbard model treated within the random-phase approximation yields a leading $d_{x^2-y^2}$ pairing state. Two phenomenological couplings motivated by oxygen breathing modes favor $s$-wave pairing when considered alone, while their attractive scattering between opposite-sign regions of the spin-fluctuation-induced gap suppresses $d$-wave pairing. Electronic charge screening strongly reduces this suppression. Comparing the $s$- and $d$-wave pairing strengths with and without screening shows that screening shifts their boundary in favor of $d$-wave pairing. Momentum-resolved projections locate the reduction of the unfavorable scattering processes. At the representative electronic interaction, $d$-wave pairing remains dominant.

The complex-temperature phase structure of the two-dimensional Potts model has remained largely unresolved despite decades of study. Here we determine the complex-temperature phase diagram of the 3-state Potts model, and uncover three gapless phases beyond the familiar ferromagnetic and paramagnetic phases. This is done by combining tensor-network free energies with the Poisson equation to compute the Fisher zero density in the thermodynamic limit. The resulting zero distribution unifies and revises long-standing conjectures on the zero loci, and outlines each phase. One is a non-conformal gapless phase, with an exponentially growing number of low-lying transfer-matrix levels distributed across all momentum sectors. The other two are Luttinger-liquid phases with complex-valued Luttinger parameters related by the Kramers-Wannier duality. At the antiferromagnetic critical point, the Luttinger parameter equals 3, while its limit from the adjacent Luttinger-liquid phase is 6. These findings reveal rich phase structures and unconventional critical behavior in a classical textbook model.

We show that a microscopic model of twisted MoTe$_2$ with screened Coulomb interactions produces topological superconductivity. Using the density matrix renormalization group (DMRG) on infinite cylinders, we find that the narrow Chern band of twisted MoTe$_2$ has Cooper pair correlations that strengthen with cylinder circumference, consistent with two dimensional superconductivity. The superconductor has chiral $f$-wave pairing around the valley center, and its topology produces chiral Majorana edge modes and non-Abelian vortex cores. At zero displacement field, the superconductor spans fillings $-0.78\lesssimν\lesssim-0.72$, close to the experimental range, and extends in a dome to finite displacement field. We then predict how its stiffness evolves with the filling and the tunable displacement field. The interplay with screening length is also examined. Nearby, we find a fractional Chern insulator, a topological charge density wave, a metal, and a charge-ordered state. The last is a candidate for the reentrant integer Hall regime in experiments. These results establish a route from repulsive interactions in microscopic models of moiré materials to chiral topological uperconductors with chiral Majorana edge modes.

The sign problem in path integrals severely limits the application of quantum Monte Carlo methods to strongly correlated systems, whereas the Lefschetz-thimble method offers a rigorous framework to mitigate this issue. In this paper, we establish a numerical framework based on the Lefschetz thimbles for the Sachdev-Ye-Kitaev (SYK) model. By introducing continuous auxiliary fields via the Hubbard-Stratonovich transformation, the partition function is formulated as a high-dimensional complex integral. We employ the complex Langevin equation to explore the saddle-point distribution of the effective action in the complex plane. Monte Carlo integration is then performed to compute the partition function and Green functions. Our numerical results exhibit excellent agreement with high-temperature expansions, infrared-limit solutions, and determinantal quantum Monte Carlo (DQMC) simulations, confirming the validity of our framework. Furthermore, we demonstrate that our complex Langevin search successfully captures the relevant dominant saddle points governing the SYK model in the large-N limit. Our work deepens the understanding of the sign problem in the SYK model and provides a viable non-perturbative avenue for addressing sign problems in more complex strongly correlated systems.

The spin splitting of an altermagnet arises from a rotation that relates its two magnetic sublattices. In known cases it is fixed by the crystal structure. Here we show that guest-host frameworks separate the two roles. In the boron-carbon clathrates $M$B$_3$C$_3$ the ordered cage removes the translation between the two guest sites and enforces cubic $i$-wave altermagnetism, with six symmetry-enforced nodal planes, for any magnetic guest in the compensated state; strain retunes the class to $d$- or $g$-wave. We demonstrate the platform with GdB$_3$C$_3$: its compensated antiferromagnetic ground state is a narrow-gap semiconductor with spin splitting of 0.19 eV at the conduction-band edge, and its weakly pinned Néel vector switches the helicity-resolved optical response on and off according to the magnetic point group. Thermodynamic stability at high pressure and dynamical stability at ambient pressure indicate that it can be made by the route established for this family.

The complete statistical characterization of defects generated in a quantum phase transition has mostly been studied in exactly solvable models, and remains little explored in nonintegrable systems such as long-range interacting spin chains. We characterize the full defect counting statistics in a nonintegrable long-range spin chain, the fractional Ising model, by time-dependent variational-principle simulations within the matrix-product-state and matrix-product-operator formalism. Fractional interactions capture emergent multiscale correlations and long-range effects in nature, governed by the order of the underlying fractional derivative. Such couplings can be approximated on existing quantum simulation platforms, where the full defect statistics can be measured. The defect number cumulants follow power laws in the driving time, with exponents that vary continuously with the fractional order. At low fractional order, deep in the long-range regime, the mean defect number decays more slowly than predicted by the Kibble-Zurek mechanism (KZM), while at intermediate order the KZM holds with the numerically extracted many-body critical exponents. At large fractional order, the exactly solvable free-fermion limit of the model becomes progressively more accurate towards the short-range limit. We derive exact bounds on the kink-number cumulants and on the kink-kink correlations, implying an exact Gaussian limiting distribution for large system sizes. In the many-body simulations, however, kinks are more correlated, consistent with the confinement of domain walls, and these bounds are violated in the long-range regime, also in the bulk of the chain. Nevertheless, the kink distribution still approaches a Gaussian for large defect numbers, with proportionality partially broken between the cumulants.

The hierarchical equations of motion (HEOM) are a powerful, numerically exact approach for simulating the time evolution of open quantum systems. Recently, several tensor-train- and tensor-network-based approaches have been suggested and realized to extend the range of applicability of the HEOM approach to larger and more realistic systems. In this paper, we build upon the existing twin-space formulation of the HEOM and present a novel approach to propagate the fermionic HEOM efficiently using the multi-configurational time-dependent Hartree method (MCTDH). We benchmark the resulting HEOM+MCTDH method with electron transport calculations for a Holstein-type nanojunction model. Furthermore, we demonstrate the capabilities of the HEOM+MCTDH method by applying it to a more complex vibronic coupling nanojunction model with strong molecule-electrode coupling for which fully quantum simulations have not been feasible previously.

We investigate light-driven inversion of orbital selectivity in a two-orbital Hubbard-Kanamori model initially in an orbital-selective Mott phase. Using the matrix-product-state we show that tailored pulse excitations drive interorbital charge transfer that nearly reverses the orbital occupations. Single-particle spectral snapshots exhibit changes consistent with an interchange of the localized and itinerant characters of the orbitals. We analyze the relevant one- and two-photon excitation channels and show that chirping around the lower-frequency resonance enhances the transfer. These results establish a microscopic route to reversing orbital selectivity and accessing long-lived nonequilibrium states in multiorbital systems.

Quantum computers have advanced to the stage that interesting many-body quantum phenomena can be observed. A striking example is the magnetic monopole paradox, a problem that was only solved recently using the language of generalized symmetries. Here we demonstrate a quantum spin chain analogue of magnetic monopole scattering on Quantinuum System Model H2. We implement a novel protocol for preparing a localized quasiparticle wavepacket on top of the ground state. By combining this protocol with time evolution across a duality defect, we efficiently convert local particles into non-local string-like excitations. We demonstrate this particle-to-string conversion in the transverse-field Ising model and a symmetry-protected topological (SPT) phase of the cluster model. Our results establish duality defects and generalized symmetries as practical resources for preparing and probing non-local excitations on quantum hardware.

Proofs generated by large language models are often difficult to understand, generalize, or formalize directly. Here we verify, digest, and formalize OpenAI's computer-assisted proof of a uniform gap in the periodic spin-1 Heisenberg chain using tensor-network diagrams and Lean 4. Reading the partition-function network along space and imaginary time shows how finite-ring bounds propagate across doubled lengths and inverse temperatures to all longer even rings, sharpening the spin-1 bound from the same finite inputs by a factor of 19 to $Δ_L > 0.0729$ for all even $L \geq 36$. We identify reflection positivity as the general condition behind the proof and certify uniform gaps in seven additional spin chains and ladders. More broadly, digesting machine-generated proofs into transparent physical mechanisms yields sharper bounds, broader theorems, and reusable formal verifications.

We present a preconditioning pipeline that accelerates Tensor Train (TT) linear solvers for symmetric positive-definite problems. Extending the Factorized Sparse Approximate Inverse framework to tensor networks, we construct a preconditioner as the product of a TT operator and its transpose. This formulation is positive semidefinite for any bond dimension and truncation, and definite whenever the factor is non-singular, satisfying the mathematical requirements of the Conjugate Gradient method. We detail the construction of this preconditioner and its integration into Density Matrix Renormalization Group (DMRG)-like linear solvers. Performance evaluations show consistent iteration reductions for anisotropic diffusion and Poisson problems across various configurations, leading to net speedups of up to 25% in ill-conditioned settings of time stepping schemes that can reuse the approximate inverses. These benchmarks were performed with the open-source Julia library DifferenTT.jl, which we also introduce in this work to provide out-of-the-box access to these preconditioned pipelines.

We establish a pressure-correlation phase diagram that yields multiple pressure-tunable quantum phase transitions in iron oxyhydroxide ($ε$-FeOOH) across three critical pressures, $P_{c1}>P_{c2}>P_{c3}$. The transitions arise from the competition among electron-correlation, crystal-field, and Hund's coupling. First-principles calculations reveal that the earlier reported high-spin (HS) to low-spin (LS) crossover at correlation-sensitive $P_{c1}$ is accompanied by a transition from charge-transfer to Mott insulating (CTI-MI) phase, and a change in the spin-quantization axis (SQA) from $\hat{y}$ to $\hat{x}$. The HS-LS crossover and CTI-MI transition is driven by the simultaneous increase in crystal-field splitting and decrease in Hund's coupling. The critical pressure $P_{c2}$ governing the structural transition from centrosymmetric $Pnnm$ phase to noncentrosymmetric $P2_1nm$ phase is found to be independent of electron-correlation. When correlation is strong, a third critical pressure ($P_{c3}$) emerges in the lower-pressure regime of the phase diagram. Validating the recent experimental observation, $P_{c3}$ is responsible for changing the SQA from $\hat{z}$ to $\hat{y}$. By exhibiting altermagnetism in both $P2_1nm$ and $Pnnm$ phases and possessing finite anomalous Hall conductivity (AHC) only for SQA along $\hat{x}$ and $\hat{y}$, $ε$-FeOOH enables a unique piezo-anomalous Hall effect where one can switch on the AHC at a lower pressure and change its direction at a higher pressure.

Quantum spin ice attracts considerable attention due to its remarkable low-energy properties, described by the effective gauge theory, which mimics the quantum electrodynamics (QED). The ``effective" QED in quantum spin ice is beyond the ``real" QED in our universe in many respects. Above all, the effective QED hosts two kinds of gapped particles, the electric monopole, in addition to the magnetic monopole. Among the two monopoles, the electric monopole is of purely quantum origin, and its character has not been well resolved. In this work, we focus on the Rokhsar Kivelson point of the quantum dimer model, as the low-energy description of the quantum XXZ model on a pyrochlore lattice, and address the excitation spectrum of the electric monopoles, by combining the quantum Monte Carlo simulation with the stochastic analytical continuation. Most strikingly, the electric monopole spectrum exhibits clear double-peak structure, corresponding to the bound states and scattering states of electric monopoles. The lower-energy peak is attributed to the exciton-like bound state of an electric monopole pair, a signal of positronium in quantum spin ice. We found that the flippable loop operator serves as a local physical observable to excite a pair of electric monopoles, implying that monopoles can be detected through lattice or optical response.

A recent remarkable paper by OpenAI showed that the spin-one Heisenberg chain with periodic boundary conditions has a spectral gap in the thermodynamic limit. The proof elegantly deduces this gap from the Hamiltonian's spectrum on a finite system. This basic idea is an inequality involving the purity of the original Gibbs state and a dual Gibbs state that is defined via space-time rotation. Here we consider a few generalizations of this purity bootstrap argument. We first extend it to arbitrary spatial dimensions. Our main result uses the space-time duality between symmetry generators and symmetry twists to obtain upper bounds on the splitting of nearly degenerate low-energy states, together with a nonzero gap above them, from finite-system information. As applications, we show that finite-system calculations for the 2+1-D toric code in a field and the 2+1-D transverse-field Ising model establish parameter regimes with four and two low-lying states whose splittings are bounded from above by $e^{-cL}$ and $e^{-cL^2}$ respectively for arbitrarily large system size $L$, while the gap above these states remains nonzero in the thermodynamic limit. Next, we extend it to fermionic systems and use a finite calculation to establish a symmetric mass generation phase in a $2+1$-D fermionic model. Finally, we apply the bootstrap to tensor-network states and establish exponential decay of connected correlations of all local operators in the symmetric square-lattice spin-one state introduced by Jian and Zaletel.

We present a general proof of different power-law decays versus distance of the eigenvalues of two-point covariance matrices in critical lattice quantum models, allowing extraction of multiple scaling dimensions. We demonstrate the practical utility of this method in the context of a spin-1/2 chain with critical ground state described by the SU(2)_1 Wess--Zumino--Witten field theory. Using quantum Monte Carlo simulations to compute covariance matrices in different symmetry sectors, we resolve the expected relevant and marginal primary operators as well as some descendants. The covariance method can be integrated into many existing numerical schemes and promises to significantly enhance computational studies of quantum criticality.

We introduce SqDMRG, which combines sample-based quantum diagonalization (SQD) with the density matrix renormalization group (DMRG). SqDMRG expands the sampled determinant space into a product of fragment spaces constructed from the samples, and solves the Hamiltonian in this product space using DMRG. We demonstrate the method on benzene pentamer and decamer clusters using classical reference samples, and on the pentamer using quantum samples obtained from the IBM Quantum Nighthawk r2 processor ibm phoenix. It achieves sub-kcal/mol accuracy with far fewer configurations than direct diagonalization.

Long-range hoppings offer a powerful mechanism to engineer exotic topological phases beyond conventional tight-binding models. Here, we experimentally demonstrate a twofold topological phase transition driven by third-nearest-neighbor couplings in a one-dimensional phononic chain that emulates cis-polyacetylene. Tunable hopping parameters in highly structured aluminum plates allow us to observe the emergence of Dirac points, a valley pseudospin structure, and topological edge states. The experimental results are supported by non-equilibrium Green's function transport calculations. Our findings confirm the predicted long-range phase diagram and establish phononic metamaterials as a versatile platform for exploring topological phenomena beyond the Su-Schrieffer-Heeger paradigm.

Preparing ground and thermal states of quantum many-body systems is a central task in quantum simulation. Dissipative protocols have emerged as a promising solution for both near-term and fault-tolerant hardware, but most constructions rely on weak system-bath coupling, leading to slow cooling. Here, we construct cooling protocols at strong system-bath coupling, and demonstrate accurate state preparation within even a single cooling cycle. Our protocols combine dynamically ramped ancilla energy levels with time-dependent couplings to optimize the trade-off between cooling error and physical time. For a general class of free-fermion systems, we reduce the dynamics exactly to a two-level scattering problem: ground-state cooling becomes a control problem of maximizing resonant excitation transfer while suppressing non-adiabatic heating. We then construct a family of ``Möbius-Ingham'' protocols, which reach optimal error scaling for gapless spectra and almost-exponential suppression of the leading heating error for gapped spectra, as a function of the cycle period. In our free-fermion benchmarks, Möbius-Ingham cooling lies on the Pareto frontier between steady-state error and cooling time among the controls tested. We further demonstrate its effectiveness in interacting Ising and $t$-$V$ chains and the two-dimensional Fermi-Hubbard model. Finally, we extend our construction to Gibbs state preparation, showing how free-fermion thermal states can be prepared in a single cycle.

We introduce a class of two-dimensional tensor network states obtained by categorifying matrix product states. This procedure leads us to consider bond spaces provided by a 2-vector space associated with the 2-cells of the underlying lattice. Tensor network representations are then defined as module categories over the free monoidal category on a monoidal signature encoding the lattice geometry. We prove a fundamental theorem stating that two representations generating the same states with suitable open boundary conditions are equivalent as module categories. Finally, we introduce a notion of 2-injectivity, a categorified analogue of the injectivity condition for matrix product states.

We consider a class of quantum circuits, known as quantum many-body spin ratchets, which act on alternating spin chains with two distinct spin types on even and odd sites. This setup breaks parity symmetry and exhibits chiral transport. It has been shown that the drift velocity in the case of integrable, SU(2)-symmetric quantum spin ratchets is universal, and conjectured that the same holds in the non-integrable case. Here, we prove that this universal drift is exact at all times, requiring only translation invariance and U(1) symmetry. Furthermore, we prove that it is sufficient for the circuit to be translation-invariant on average by deriving the same result in the random case. Finally, we demonstrate that breaking translation invariance immediately results in deviations from this universal prediction.

Tensors are higher dimensional generalisations of vectors and matrices. Tensor networks are a diagrammatic representation of complicated relations between tensors, and tensor contraction is the evaluation of a tensor network analogous to matrix multiplication. Tensor networks and tensor contraction are ubiquitous in computer science, physics, and applied mathematics. Most existing tensor contraction algorithms are focused on finding an optimal contraction order of the tensors' bonds. However, finding the optimal order is an NP-hard problem. It is equivalent to another possibly more straightforward but still NP-hard problem of finding the treewidth of the line graph of the tensor network. While in general finding the optimal bond-by-bond contraction order is hard, in many practical cases there is no need to find an optimal order. This observation leads to the core innovation of this thesis, which is a new tensor contraction paradigm that uses graph theoretic operations to produce an equivalent tensor network whose line graph has lower treewidth than the original tensor network, before applying the bond-by-bond contraction algorithm. We call these graph theoretical operations \textit{rewrite preprocessing}, and we name the algorithms that apply rewrite preprocessing before regular bond-by-bond contraction \textit{Tensor Contraction with Graph Operation}. There are many possible ways to perform rewrite preprocessing; in this thesis, we focus on the graph operations of ZX calculus, a diagrammatic language for quantum circuits. Ample numerical simulations have been performed to demonstrate the effectiveness of our proposed method.

Neural quantum states provide accurate variational representations for interacting lattice fermions and spins. Turning a new ansatz into a reproducible calculation, however, requires coordinated implementations of the Hamiltonian, sampling, optimization, observables, and accelerator parallelism. We present LaQX (Lattice Quantum simulation based on jaX), an open-source neural quantum state library based on JAX that organizes these components through common functional interfaces. Beyond ground-state variational Monte Carlo, LaQX supports excited-state calculations, Green's Function Monte Carlo refinement, adiabatic evolution, and real-time evolution with time-dependent variational Monte Carlo. It consolidates methods developed in recent high-accuracy neural quantum state studies of the Hubbard model and frustrated magnets, and the package has been used to investigate the triangular-lattice Hofstadter-Hubbard model. We document these capabilities with reproducible calculations for Hubbard, $J_1$-$J_2$ Heisenberg, altermagnetic Hubbard, and interacting Hofstadter models. Function-oriented modules separate the scientific components and limit the source context needed for extensions by human developers and large language models. Source code, documentation, and executable examples are available at https://github.com/bytedance/laqx.

Electronic band structures govern how crystals conduct electrons and absorb light. Density-functional theory (DFT) computes them through costly self-consistent iterations, while existing machine-learning models either predict isolated scalar properties or output Hamiltonians that still require DFT solvers. Here we present EpNet, a learned empirical pseudopotential for crystals spanning 63 elements across the \textit{s}, \textit{p} and \textit{d} blocks. An equivariant neural network builds its Hamiltonian from physics-motivated local potentials and nonlocal projectors. Trained solely on band energies, EpNet predicts complete band structures directly from atomic coordinates. On unseen chemical systems, its median errors are \winAll\ eV for band energies and \gapAll\ eV for band gaps. On 600,000 structures from an independent database, its band-gap errors are overall less than half those of existing models, at about 1\% of DFT cost. EpNet thus enables fast screening of electronic properties for materials discovery and electronic-structure simulations of large systems.

Designing a superconducting quantum processor rests on predicting the low-lying spectrum of its coupled circuit. The Hilbert-space dimension of the circuit's lumped-element model grows exponentially with the number of circuit nodes. The usual methods split the circuit into subsystems, but that split is hard to choose when the elements are strongly coupled. We present a tensor-network method that chooses no subsystems and whose memory and run time per solve are set by the number of modes and two truncations, the Fock cutoff of each normal mode and the bond dimension of a matrix product state (MPS), not by the Hilbert-space dimension. From the capacitance matrix, the inverse-inductance matrix and the junction list alone, the linear circuit is diagonalized into normal modes. Each junction cosine is then a sum of two products of single-mode operators, so the complete lumped-element Hamiltonian is a matrix product operator (MPO) of small bond dimension with one site per mode, solved by the density matrix renormalization group (DMRG). The normal modes make the method generic, and the tensor network makes it affordable. For two fluxonia joined by a tunable coupler, the method agrees with exact diagonalization of the same truncated space to within 1 MHz. Our circuit of three fluxonia joined by two couplers has 11 normal modes and a truncated space of $10^9$ states, out of reach of exact diagonalization in dense or sparse storage. At a bond dimension of 32, its lowest eight levels are converged to within 1 MHz in both truncations, and each state is stored in a small fraction of the numbers a dense vector needs. One solve at that bond dimension takes 81 min and 1.3 GB of memory on one CPU node.

Machine learning from local atomic environments offers an efficient route to electronic structure, but evaluating global electronic properties in large systems remains computationally demanding. Here we combine real space charge density learning with stochastic propagation using Kohn-Sham Hamiltonians constructed from the predicted densities. The framework avoids target system self-consistent iterations and explicit diagonalization, while hierarchical blocking enables dense grid inference within limited GPU memory. We demonstrate density of states calculations for one million Si atoms and optical and transport calculations for over one hundred thousand atoms, without training separate models for each property. This framework extends the practical scale of electronic structure calculations within Kohn-Sham DFT, opening opportunities to investigate how atomic structure governs electronic, optical and transport properties at length scales relevant to real materials and devices.

We extend the exact fixed-point tensor construction of Cheng et al. to fermionic minimal models, providing an analytic construction of fixed-point tensors for fermionic rational conformal field theories. The fermionic fixed-point tensor is defined as a correlation function of boundary-changing operators on a triangle, with careful treatment of fermionic parity and spin structure. Boundary tensors are also constructed and combined with triangle tensors to form strips, enabling the extraction of scaling dimensions for boundary-changing operators. The construction incorporates modified weights for shrinkable boundary conditions, systematic depiction of spin structures, recovering both the Neveu-Schwarz and Ramond sectors. As a benchmark, our tensor with a minimal truncation of descendant levels yields convergence to the exact bulk scaling dimensions of many low energy states for the free Majorana fermion ($c=1/2$), demonstrating the efficiency of approximate fixed point tensors from boundary changing operators. We also make use of a more general description of fermionic CFTs based on fermionisation of rational bosonic CFTs, and make a generic prescription for discretisation making use of their connections to modular tensor categories. We show that the discretised spin structures on a spin surface naturally recover the signs needed to recover the R/NS sectors. We propose that the collection of OPE coefficients of the discretised 2d CFT can be naturally interpreted as a state-sum of spin-TQFTs, analogous to the bosonic case. This work provides analytically computed starting points for tensor network simulations of fermionic critical phenomena.

Rutile RuO2 has emerged as a prominent candidate altermagnet, yet the microscopic origin of its reported charge-angular-momentum conversion remains unsettled. Recent theoretical and experimental advances have highlighted a potentially dominant role of orbital transport, raising a fundamental ambiguity between spin- and orbital-mediated conversion in RuO2. Here, we report the first experimental disentanglement of the spin- and orbital-to-charge conversion channels in RuO2 heterostructures using terahertz emission spectroscopy. We identify Fe as a predominantly spin-current injector, whereas Ni and permalloy predominantly inject orbital currents, revealing a substantial orbital contribution to ultrafast terahertz emission. Under controlled spin-current injection, we find no evidence for inverse altermagnetic spin-splitting effect. Instead, the observed anisotropy is consistent with an anisotropic inverse spin Hall response that is strongly conditioned by the near-interface environment. Under orbital-current injection, the anisotropy is markedly weaker, indicating that orbital-to-charge conversion is largely isotropic in RuO2. Our results establish orbital transport and interface-conditioned spin Hall conversion as key to terahertz emission from RuO2, underscoring the need to disentangle spin, orbital, and interfacial contributions before assigning anisotropy to altermagnetism.

Symmetry-protected topological (SPT) phases are the simplest example of quantum phases that cannot be detected by local order parameters. Their classification is by now well-established from multiple perspectives, a prominent one being via tensor network states. The practical question, however, is not mere classification, but detection: given the ground state of an interacting Hamiltonian and a symmetry group $G$ along with its representation, is there an efficient way to tell if the state is in a SPT phase with respect to $G$? While there exist such practical tests in one spatial dimension (1D), methods in 2D have so far been restricted to finite-size systems or require an explicit form of virtual symmetries of the tensor networks representing the ground state. Here, we develop an agnostic and efficient method to detect interacting $\mathbb{Z}_2$ 2D SPT phases, using an infinite projected entangled-pair (PEPS) representation of the ground state. The method is agnostic to the symmetry representation in the virtual space of the PEPS tensor and applicable independent of the bond dimension of the tensors. Moreover, it is efficient as it does not require finite-size scalings: it works directly in the thermodynamic limit. We numerically demonstrate the procedure by mapping out phase boundaries away from fixed-point wavefunctions. We extract the SPT label cocycle data explicitly from the PEPS tensor.

Efficient encoding of data into quantum states remains a fundamental bottleneck in many quantum computing protocols, which would otherwise offer advantages over their classical counterparts. Here, we introduce deterministic and scalable amplitude-encoding methods with low two-qubit gate costs and tunable trade-offs between circuit depth and encoding fidelity. Our approaches require no amplitude amplification or mid-circuit measurement and are capable of encoding arbitrary complex data or queryable functions by using tensor-networks to go beyond what can be held in classical memory. Our principle contributions are two complementary algorithms for simplifying fully controlled multiplexer circuits to one- and two-qubit gates. In the first, we formally map the elimination of redundant gate controls to a Boolean logic minimization problem, while in the second we utilize a single ancilla qubit to map the multi-controlled rotation gates to a sequence of one- and two-qubit operations via the Walsh-Hadamard spectrum of the gate angles. We benchmark these methods across a broad range of function classes, characterizing their total depth in an all-to-all connectivity graph against system size and target fidelity. Hybridizing our methods gives up to orders of magnitude improvement in circuit depth at equivalent state fidelity compared to the default state-preparation method in the Qiskit package. We furthermore demonstrate our method on IQM hardware with a square lattice connectivity, finding excellent fidelity, even with minimal noise mitigation. By significantly reducing gate depth for encoding, this work enables the experimental execution of applications on current and near-term hardware which would otherwise be reduced to noise before the first calculation step, shortening the time horizon to robust quantum utility.

Altermagnetism combines vanishing net magnetization with momentum-dependent spin splitting, enabling spin-transport functionalities without ferromagnetism. Although rutile RuO$_2$ has been proposed as a d-wave altermagnet, its magnetic ground state remains highly controversial. Moreover, existing transport verifications suffer from severe ambiguity due to the entanglement of spin-splitting effects and time-reversal-even spin-Hall contributions. Here, we introduce a two-pronged strategy to break the current verification impasse. First, we establish a definitive, symmetry-based transport criterion: by identifying specific crystallographic geometries where the spin-Hall background is strictly forbidden, we provide a zero-background window for unambiguous altermagnetic detection. Concurrently, to bypass the uncertain bulk energetics, we propose the RuO$_2$/MnO$_2$ superlattice as a prototypical platform to robustly enforce the d-wave altermagnetic order via interfacial magnetic proximity. Within this platform, symmetry-resolved spin-transport analysis reveals unconventional responses, including collinear spin currents enabling field-free magnetization switching and sizable spin-polarized charge currents. By combining a universal transport criterion with a robustly engineered material platform, our work transcends the debates over bulk properties, providing an actionable, definitive experimental roadmap for verifying d-wave altermagnetism.

Chiral Floquet phases support quantized chiral flow along their boundaries even when the bulk is fully localized. In this Letter, we establish a complementary phenomenon, which we call Hilbert-space pumping: under a boundary condition twisted by time translation, a part of the boundary Hilbert space, of quantized dimension, is transferred between the opposite boundaries. The twist is the temporal analogue of Laughlin's flux insertion in the quantum Hall effect. For noninteracting fermions, the number of transferred boundary modes is determined by the winding number, as we show analytically for the boundary dynamics and numerically for the anomalous Floquet Anderson insulator. For interacting systems whose boundary dynamics is a quantum cellular automaton, the transferred dimension is set by the Gross-Nesme-Vogts-Werner (GNVW) index: the twisted boundary Floquet operator is the Floquet operator on a tilted time slice, whose Hilbert-space dimension is multiplied by the index for each period of twist. We derive this relation for matrix product unitaries and test it numerically in random unitary circuits. Applied to the matrix product operators of non-invertible symmetries, the same twisted rank returns the quantum dimension of the symmetry operator in addition to its flow index.

Heat transport in solids is usually described in a phenomenological way via Fourier's law. However, at small time and length scales in low thermal conductivity materials, a microscopic description is needed. This is traditionally done by invoking the phonon Boltzmann Transport Equation (BTE), which has been successful in describing transport in high thermal conductivity materials. Recently, a generalization of the BTE, the Wigner Transport Equation, was developed to be able to microscopically describe transport in low thermal conductivity materials by taking into account not only phonon scattering processes but also phonon tunneling processes, so-called coherences. We present a Python package that solves the Wigner Transport Equation (WTE) in spatial and temporal Fourier space for arbitrary source terms, thereby enabling the study of heat transport at small space and time scales. This allows one to compute both static and dynamic thermal conductivities from bulk to nanoscale. Beyond that, it can be used to study the response of phonon systems in materials to arbitrary heat sources.

We present a systematic framework for gauging continuous symmetries, and show that the Schur-Weyl duality arises naturally as a realization of this construction with open boundary conditions. Restricting to compact Lie group symmetries, we construct dual models for one-dimensional quantum lattice models exhibiting those symmetries. We demonstrate that matrix product operators enable the explicit construction of intertwiners between the dual theories and their (non-invertible) symmetries. With closed boundary conditions, we can relate the full spectrum of these dual theories by considering all possible twisted boundary conditions.We illustrate how the Peter-Weyl transform appears when the duality acts on the boundary Hilbert space of a $G$-symmetric model. The presented framework generalizes previous work on dualities for unitary fusion category symmetric models. It exemplifies how dualities for both finite and continuous symmetries rely on the relation between a bond algebra and a representation of the symmetry in the lattice as double centralizers. We argue that dual non-invertible spontaneous symmetry breaking is the mechanism that underlies symmetric tensor network algorithms. We illustrate these results for $U(1)$- and $SU(2)$-symmetric quantum lattice models.

We present a systematic framework for gauging continuous symmetries, and show that the Schur-Weyl duality arises naturally as a realization of this construction with open boundary conditions. Restricting to compact Lie group symmetries, we construct dual models for one-dimensional quantum lattice models exhibiting those symmetries. We demonstrate that matrix product operators enable the explicit construction of intertwiners between the dual theories and their (non-invertible) symmetries. With closed boundary conditions, we can relate the full spectrum of these dual theories by considering all possible twisted boundary conditions. We illustrate how the Peter-Weyl transform appears when the duality acts on the boundary Hilbert space of a $G$-symmetric model. The presented framework generalizes previous work on dualities for unitary fusion category symmetric models. It exemplifies how dualities for both finite and continuous symmetries rely on the relation between a bond algebra and a representation of the symmetry in the lattice as double centralizers. We argue that dual non-invertible spontaneous symmetry breaking is the mechanism that underlies symmetric tensor network algorithms. We illustrate these results for $U(1)$- and $SU(2)$-symmetric quantum lattice models.

We study a Bayesian inverse problem in EUV reflectometry, where optical constants and layer thicknesses are inferred from angle- and wavelength-resolved reflection measurements. The physical process is modeled and evaluated using an efficient forward model based on the transfer matrix method. As an alternative to computationally expensive sampling-based inference, we employ an Eulerian, entropically regularized Jordan-Kinderlehrer-Otto (JKO) scheme that minimizes the Kullback--Leibler divergence via proximal steps in the Wasserstein space. Using the Benamou--Brenier dynamic formulation and a Hopf-Cole transform, each JKO step is reformulated as a pair of heat equations coupled only through their initial and terminal conditions, enabling an efficient fixed-point solver. The resulting high-dimensional densities are discretized on tensor-product grids and represented in tensor-train format, exploiting Kronecker structure for scalable computation. Preconditioning based on estimated posterior moments further improves tractability. We test the method with numerical experiments on EUV reflectometry data for a ruthenium--silicon ($\ce{Ru}$ on $\ce{Si}$) specimen and compare the learned tensor-train posterior to ensemble Markov chain Monte Carlo (MCMC) reference samples. We demonstrate a use case in which we are able to reuse the tensor-train posterior to obtain samples for shifted target distributions with less computational effort than a comparable MCMC method.

Determining the entanglement Hamiltonian (EH) of an interacting quantum many-body system remains a fundamental challenge beyond special analytically tractable cases. Here we develop an observable-matching framework for learning effective EHs, which we demonstrate using response estimates from quantum Monte Carlo. As applications of this framework, we first calibrate the overall prefactor of the lattice Bisognano--Wichmann (LBW) construction using equal-time correlations and then learn spatially resolved couplings beyond its prescribed profile, without assuming a specific functional form for their spatial dependence. Trained solely on local correlations, the learned EHs reproduce benchmarked von Neumann and second Rényi entropies in $1$D and second Rényi entropies in $2$D. For a strong-bond cut associated with special surface criticality at the $(2+1)$D $O(3)$ Wilson--Fisher critical point, free-form learning outperforms the calibrated LBW ansatz in predicting the second Rényi entropy and reveals spatially nonuniform coupling renormalizations. Importantly, the learned EH also provides access to the von Neumann entropy in this setting, revealing a negative subleading correction and extending previous Rényi entropy findings to the von Neumann case. Directly evaluating the von Neumann entropy in such $(2+1)$D interacting quantum many-body systems remains challenging for conventional numerical methods. These results establish observable-based learning as a route to constructing effective entanglement Hamiltonians beyond prescribed spatial profiles in higher-dimensional interacting systems.

Krylov subspace methods efficiently approximate extremal eigenstates and eigenvalues of many-body Hamiltonians, but non-orthogonal bases can suffer from ill-conditioning. This challenge is amplified in quantum computing, where direct state access and thus straightforward orthogonalization are unavailable. We investigate a variational Krylov-like framework that iteratively optimizes a parametrized sequence of states to mimic Lanczos-like couplings while penalizing linear dependencies. Analyzed against exact Krylov and Lanczos on a finite Anderson impurity model, ansatz expressivity strongly influences the performance: Hamiltonian-inspired ansätze preserve Krylov-like convergence while improving the conditioning, whereas shallow ansätze remain stable but miss low-energy directions.

The anyon-Hubbard model is the paradigmatic lattice model of one-dimensional (1D) anyons, in which chiral interactions make bosons equivalent to particles with fractional statistics. Its continuum limit, the chiral BF theory, emerges as the low-energy description of a Raman-coupled Bose gas with imbalanced interactions. Here, we show that loading such a gas into a 1D optical lattice realizes the lattice counterpart of this theory, which we name the Raman anyon-Hubbard model. The Raman-dressed bosons acquire density-dependent Peierls phases when tunneling through occupied sites. For a small statistical phase, the model reduces to the anyon-Hubbard model, while beyond this limit it defines a distinct lattice anyon model with unbounded rather than periodic statistical interactions. Numerically, we show that the dressed gas displays the defining properties of 1D anyons, namely a density-dependent asymmetry of the ground-state momentum distribution, interaction-induced asymmetric transport and chiral bound states. Remarkably, these features are readily accessible with Raman-coupled $^{39}$K atoms, opening up the experimental exploration of 1D anyons beyond the few-particle regime.

We present a unified theoretical study of the interplay between Coulomb disorder and Coulomb interactions in two-dimensional electron systems, covering the Wigner crystallization transition, the metal-insulator transition, and the crossover from a Wigner crystal to an amorphous localized solid in both zero magnetic field and the lowest Landau level. Using Hartree-Fock theory, variational Monte Carlo, and percolation and strong-localization arguments, we show that quenched charged disorder destroys true long-range Wigner crystalline order in the thermodynamic limit at any nonzero charged-impurity density, within the Imry-Ma-Larkin framework; the finite-system density image is that of an amorphous solid with crystalline coherence over at most a few lattice spacings. For the two-dimensional metal-insulator transition we show, in a controlled nonlinear-Thomas-Fermi model of long-range Coulomb disorder at setback d, that the half-depletion density of the semiclassical landscape and the Ioffe-Regel-Mott density of Born-RPA transport track each other to within a factor of order unity over the physically relevant range of impurity density and setback, with a slow n_i^(-1/10) d^(-1/5) spread; the agreement reflects the common energy scale E_dis ~ E_F at the transition rather than a theorem, and modest independent short-range transport channels can shift the ratio by factors of order two. The article is a theoretical study under well-defined approximations; it takes the recent scanning tunneling microscopy imaging experiments on two-dimensional electron systems and the longstanding disorder-versus-interaction debate on the two-dimensional metal-insulator transition as motivation, and leaves to the reader the comparison between its calculated signatures and the published experimental analyses.

Triangular magnetic lattices have attracted significant attention as a model platform for exploring frustrated magnetism and quantum spin liquids, particularly in rare-earth-based systems. These materials also hold considerable promise for cryogenic magnetic refrigeration. In this work, we investigate the magnetic and magnetocaloric properties of the rare-earth (RE) based compounds $\text{Ba}_3\text{RE}(\text{BO}_3)_3$ ($\text{RE} = \text{Dy}$ and $\text{Ho}$). Polycrystalline samples synthesized via a solid-state reaction method crystallize in a hexagonal structure with the space group $P6_3mc$, featuring triangular nets separated by approximately $8.8~\textÅ$ along the $c$-axis. Detailed magnetic investigation reveals a deviation from pure paramagnetic behavior at low temperatures ($T < 20\text{ K}$). This is further corroborated by isothermal magnetization measurements, which exhibit non-overlapping $M$ vs. $H/T$ plots and reduced magnetization moments, suggesting significant crystal electric field (CEF) splitting of localized $f$-states. For an applied magnetic field change of $70\text{ kOe}$, the maximum entropy change (Relative Cooling Power) reaches $4.2\text{ J kg}^{-1}\text{ K}^{-1}$ ($204\text{ J kg}^{-1}$) and $7.3\text{ J kg}^{-1}\text{ K}^{-1}$ ($230\text{ J kg}^{-1}$) for the Dy- and Ho-based compounds, respectively. Overall, the positive slopes of Arrott plots, the universal scaling behaviour of $ΔS_{\text{m}}/ΔS_{\text{m}}^{\text{max}}$ versus reduced temperature ($t$), and the field exponent $c$ of the RCP versus $H$ curves indicate a second-order nature of entropy distribution, making them promising candidates for magnetocalorics in the millikelvin range.

We demonstrate that a single machine-learning correction, applied to bandgaps from density functional theory calculations with two of the most frequently used energy functionals, results in bandgaps that are close to the experimental values. This solves a classical problem of computationally efficient electronic structure methods, to give reliable bandgaps that can be compared to experiments. The employed approach relies on materials specific feature engineering of machine learning models, and is shown to have close accuracy to that of much more expensive calculations based on the GW-approximation or hybrid functionals. The findings reported here open the possibility to use efficient calculations in materials design workflows that need reliable bandgaps, with implications for investigations targeting large computational databases or for investigations that are entirely predictive in nature.

The correlation between two subsystems $A$ and $B$ of a many-body system is usually measured through a pair of observables chosen in advance. The Hirschfeld-Gebelein-Rényi maximal correlation removes the choice: it is the largest correlation coefficient between any function of the configuration of $A$ and any function of the configuration of $B$, and the two functions that reach it, the optimal transformations, are the observables that carry the correlation. Both follow from a singular value decomposition (SVD) of the joint distribution normalized by its marginals, a matrix whose dimension is the number of configurations. We show that this SVD can be computed variationally for a classical distribution of discrete variables, such as the Boltzmann distribution of a spin model, given in the tensor-train (matrix-product) format, without ever forming the matrix. The optimal transformations are sought as tensor trains; each step updates a single core in closed form, and the inverse of a small matrix contracted from the distribution and the remaining cores performs the division by the marginals core by core. Every iterate is a feasible pair, so the objective is a lower bound on the maximal correlation that increases monotonically, and the cost grows linearly with the number of variables. The algorithm reproduces the exact SVD on small systems, recovers the known nonlinear optimal transformations of synthetic distributions of up to ten angular variables ($10^{15}$ configurations), and, for the bipartite antiferromagnetic model with $2^{32}$ configurations, returns an optimal transformation that is a function of the group magnetization, linear above the transition and sign-like below it: the staggered order parameter is read off from the distribution instead of being postulated.

Machine learning supports materials design, atomistic simulation and experimental characterization, yet models and workflows remain fragmented. We present PaddleMaterials, an open-source framework integrating heterogeneous materials-learning tasks while preserving task-specific representations and scientific semantics. It unifies data processing, training, evaluation, pretrained-model management and workflows across six task classes: property prediction, structure generation, machine-learning interatomic potentials, electronic-structure fields, spectrum elucidation and microscopy image processing. Task interfaces record physical quantities, units, aggregation conventions and references, making cross-task composition verifiable. A shared compiler backend enables switching from dynamic-graph to compiled execution through one setting, using graph optimization, operator fusion and hardware-targeted code generation without changing model definitions. Across 65 model instances in 11 groups, compiled execution accelerated every instance in eight groups, with group median speedups of 11%-224% and group maximum speedups of 11%-336%. Case studies on LiMnO2 and toluene link structure generation or spectrum elucidation, interatomic-potential relaxation, property prediction and electron-density prediction, using same-geometry density functional theory (DFT) calculations as references. For LiMnO2, predicted formation energies for generated and relaxed structures differed from same-geometry PBE+U references by at most 7 meV per atom, while the computational cost of model steps was less than one-twentieth of a DFT single-point calculation. PaddleMaterials extends unified materials workflows beyond conventional property, potential and structure-generation tasks while supporting verifiable cross-task operations and compiled execution.

We propose a Thermodynamic Bethe Ansatz (TBA) approach for describing thermalization in integrable quantum field theories. The exponential damping and oscillations when approaching thermal equilibrium can be understood as the analogs of quasi-normal mode (QNM) frequencies of black hole perturbations within the framework of holography. Holographic intuitions and some general properties of thermal Green's functions motivate the existence of such a TBA formulation, as well as lead to the imposition of a nonstandard analytical structure on the QNM TBA, while comparison with integrable form-factor expansions fixes the remaining ingredients. We conjecture the QNM TBA to be equivalent to performing an all-order resummation of secular terms in a thermal form factor expansion of retarded 2-point correlation functions. We compare the predictions with tensor network simulations of the E8 deformation of Ising CFT at nonzero temperature. We conjecture that the QNM TBA framework should also be applicable for describing out-of-equilibrium approach to Generalized Hydrodynamics (GHD).

In this paper we propose an ensemble model based on a swarm of dense neural networks, able to predict Heisenberg exchange interactions within less than 1% error. The model uses the inputs of the ab-initio Schrödinger equation, i.e. the compound structure and chemical composition. In this first work, we trained the model on a few compounds of the hexagonal manganite family and were able to infer, with similar accuracy, the J values for other compounds of the same family at randomly chosen geometries, thus showing that we accurately learned the effect of the chemical composition and structure on effective magnetic integrals J.

We study 3-tensors satisfying the ``zipper condition,'' which play an important role in tensor network approaches to two-dimensional topological order. In our previous work, we identified the 4-tensors appearing in these approaches with bi-unitary connections in the Jones theory of subfactors in operator algebras. These bi-unitary connections give rise to subfactors of finite index. We prove that, when the associated subfactor has finite depth, the 3-tensors satisfying the zipper condition correspond precisely to the quantum $6j$-symbols, also known as $F$-symbols, of the unitary multi-fusion category arising from the subfactor. This realization of quantum $6j$-symbols differs from their realization through a flat bi-unitary connection arising from the higher relative commutants of a subfactor.

Symmetry governs the organization of physical structures, while symmetry breaking reveals their transformations. In geometric deep learning, symmetry imposes rank constraints on equivariant representations, limiting the expressive capacity of equivariant graph neural networks (GNNs). Here we show that these constraints can instead be exploited: the singular-value spectra of equivariant representations provide a unified fingerprint of structural symmetry. We introduce spectral symmetry fingerprints (SymFP), which enable discrete point-group identification through rank-template matching and continuous symmetry analysis by quantifying spectral energy beyond target rank constraints. SymFP is orientation invariant and architecture agnostic, enabling a parameter-free GNN implementation. Its spectra can be reused across multiple target groups, and the resulting leakage scores satisfy subgroup monotonicity. Across synthetic structures, molecules, coordination complexes and large biological assemblies, SymFP accurately identifies point groups, tracks controlled symmetry breaking and reveals local and degree-resolved symmetry changes that can be obscured by global geometric measures. These results establish the spectral structure of equivariant representations as a general framework for characterizing symmetry across physical systems.

The absence of magnetic Bragg peaks and a broad excitation continuum are commonly recognized experimental features of a quantum spin liquid. The Ce-based triangular-lattice hexaaluminate CeMgAl$_{11}$O$_{19}$ exhibits both signatures, yet its effective model and ground state remain under debate. Experiments have variously suggested a quantum spin liquid, proximate umbrella point, and disorder-suppressed magnetic order. Here we trace these contradictions to an overlooked multipolar degree of freedom: the transverse pseudospin of the $|\pm5/2\rangle$ Kramers doublet is a rank-5 triakontadipole rather than an ordinary dipole. Analyzing the magnetization, neutron spectra, and thermodynamics on equal footing, we establish a dipole-triakontadipole XXZ model with $J_{\perp} \simeq 0.65$ K and $J_c \simeq -0.13$ K. The resulting planar triakontadipolar superfluid order is hidden from dipolar probes, explaining the absence of a magnetic Bragg peak in neutron scattering. The zero-field excitations separate into a low-energy triakontadipole wave and a broad higher-energy continuum. We further propose nuclear spin-lattice relaxation and spin-current experiments to probe the gapless excitations of this hidden triakontadipole superfluid. Our findings thus open a new avenue for high-rank multipolar physics in rare-earth triangular-lattice magnets.

Predicting magnetic ordering temperatures requires accurate exchange interactions and an appropriate treatment of spin thermodynamics. Here we assess these requirements through complementary quantum Monte Carlo benchmarks. For five insulating antiferromagnets, we compare ferromagnetic-antiferromagnetic energy differences from diffusion Monte Carlo (DMC) and density functional theory with references reconstructed from inelastic neutron scattering (INS) exchanges. DMC underestimates these references but achieves a mean absolute error comparable to the best-performing functionals tested. Separately, we directly test the previously proposed (S+1)/S rescaling of classical Monte Carlo temperatures against quantum stochastic series expansion simulations for 17 compounds using corresponding INS-derived Heisenberg Hamiltonians. Rescaling reduces the mean absolute classical-quantum temperature difference from 81.7 to 8.2 K, with a mean absolute percentage difference of 5.9 %. These results establish a direct quantum benchmark for the rescaling within the tested models and emphasize that improvements in electronic magnetic energy scales should be accompanied by attention to quantum spin thermodynamics.

Weak first-order transitions can exhibit extended pseudocritical scaling regimes that mimic genuine criticality. We show that the long-debated apparent criticality of the dilute Baxter--Wu model arises from a marginally relevant renormalization-group flow, which ultimately drives the system into a first-order regime with phase coexistence. By integrating conformal perturbation theory with symmetry-resolved tensor-network level spectroscopy, we parameterize a slowly running coupling in a numerically accessible form and construct a generic flow-based diagnostic. The renormalization-group trajectories of the dilute and self-dual Baxter--Wu models under marginal operator deformations collapse onto the multicritical phase diagram of the dilute four-state Potts model, exposing the flow structure on both the marginally relevant and irrelevant sides of the Potts point.

X-ray free-electron lasers (XFELs) provide the temporal resolution needed to extend total scattering and pair distribution function (PDF) analysis to investigate structural dynamics on femtosecond-to-picosecond timescales. We recently demonstrated that high-quality total scattering data and PDFs can be obtained from individual XFEL pulses; here, we take the next step towards time-resolved studies of structural dynamics in solution. Using the Materials Imaging and Dynamics (MID) instrument at the European XFEL, we obtain total scattering data up to $Q_{\max} = 16.4\,\mathring{\mathrm{A}}^{-1}$ from both conventional capillary samples and, critically, from samples delivered by a liquid jet. Despite the large air and solvent backgrounds inherent to this setup, PDFs suitable for quantitative structural analysis are recovered from liquid dispersions and molecular clusters in solution. We further demonstrate pump-probe total scattering from $\mathrm{CsPbBr}_3$ nanocrystals in a liquid jet, resolving a fluence-dependent photoinduced isotropic lattice strain. Additionally, empirical potential structure refinement (EPSR) of XFEL total scattering data from an aqueous Keggin cluster reveals site-specific hydration environments and enables three-dimensional mapping of the surrounding solvent structure. Together, these advances open the door to tracking not only ultrafast structural changes in solution, but also the evolving solvent environments that govern molecular behaviour and reactivity.

Ground states of Sachdev-Ye-Kitaev (SYK)-type models are volume-law entangled, which makes them costly both for tensor networks and for quantum hardware simulations. Thus, it remains unclear where quantum hardware could outperform classical simulation of such models. In this work, we study a spin SYK-dot chain, in which q dense SYK dots of p qubits are joined by single ZZ bonds of strength K. The entanglement across any cut is then bounded by floor(p/2) ln 2 plus the inter-dot entanglement, and the two-qubit depth of one Trotter step is set by p alone. A resource estimate shows that a quantum advantage for the ground state could appear only for dots of about 15 to 21 qubits or more. There, a single Trotter step already requires two-qubit error rates of order 10^{-6}, beyond the reach of noisy intermediate-scale quantum (NISQ) hardware. We test the assumptions of this estimate with an imaginary-time matrix product state (MPS) that applies each dot as one block. It agrees with exact diagonalisation at (p, q) = (6, 2). At (5, 2), a variational circuit with the same dot-plus-link layout comes within 3.6% of the exact energy, and noisy simulations show that its fidelity decays as assumed in the resource estimate. At (8, 5), i.e., N = 40 qubits, the energy converges at bond dimension chi = 768, far below the exact value 2^20, and the inter-dot entanglement stays small and is set mainly by K. Out-of-time-order correlators across the links decay on a time scale that shrinks rapidly with |K|. Thus, the dot size p and the inter-dot coupling K together set the entanglement of the chain, and with it where a quantum advantage could appear.

We address a central open question in driven topological matter by investigating the nonequilibrium response of a two-dimensional magnet with Dzyaloshinskii-Moriya interactions subjected to a periodically oscillating magnetic field. Using extensive Monte Carlo simulations, we uncover a dynamically reentrant skyrmion phase that emerges within a finite window of field amplitudes and bias fields, sandwiched between the conventional skyrmion and polarized ferromagnetic phases. This phase has no equilibrium counterpart and originates from an imbalance between skyrmion creation and annihilation over a single field cycle, leading to a reversal of the period-averaged topological charge. We show that the re-entrant phase exists only in the dynamic response regime where the system can follow the drive, and is suppressed at higher driving frequencies. Our results establish periodic driving as a nonequilibrium control parameter for stabilizing and reversing skyrmion topology.

Single-particle spectroscopy has long served as a criterion for the strength of electronic correlations: strong collective motion is expected to shift and broaden spectral peaks and to renormalize the dispersion. Here we show that this correspondence can fail on the same microscopic orbital. The mechanism relies on local interference between two resonant side orbitals: one suppresses low-energy real transitions into the interacting environment, while the other suppresses spectral-weight backflow to the source orbital; virtual mixing nevertheless transfers the environmental spin polarization to the source. We prove a rigorous relative spectral bound that is uniform in environment size and interaction strength and holds for arbitrary stationary states; the protection becomes stronger as the relative frequency resolution improves. In a fully quantum two-cell device, a collective mode softens by 60%, whereas the particle intensity at the free peak changes by only 3.6 parts per million. Quantum Monte Carlo shows that the transferred spin response grows with spatial correlations and is quantitatively predicted by the free-band polarization. These results recast the meaning of a simple particle spectrum and provide a design principle for quantum probes that combine predictable particle spectra with strong collective sensitivity.

Phase identification is a computational task involving the classification of material structure using the positions of atoms. Using data from Molecular Dynamics (MD) simulations, our task is to identify the material phase of a set of atoms, based on how the simulations predict their positioning relative to each other. By leveraging existing research on mechanistic neural network theory, particularly Neural Tangent Kernels (NTK), we develop a new and efficient method for classifying simulated material samples. Our method involves defining a set of "atom detectors" that are given positions around the origin in cartesian coordinates. The sum of Gaussians centered at each atom is calculated at each atom detector, and is used as the input to a neural network. We present theoretical and experimental evidence to show that this method is highly effective and efficient at the task. Our method obtains better or similar accuracy scores to commonly used methods while having a major computational speedup.

We develop a hybrid tensor-train/quantized tensor-train (TT/QTT) upwind discontinuous Galerkin method for three-dimensional transient multigroup neutron transport. The central contribution is a fixed-matrix upwind construction: two spatial DG matrices serve all directions and energy groups, with signed and absolute direction-cosine factors preserving Kronecker structure. For homogeneous materials on uniform Cartesian grids, operator rank bounds are independent of mesh size, group count and direction count at fixed scattering order. We derive a conditional global step-defect bound and enforce group-integrated balance of the compressed system. Crank-Nicolson tests show second-order convergence on a separable manufactured solution whose finest final state occupies 13 kB instead of 35 TB, a compression factor of $2.7\times10^9$. A modified rank-adaptive basis-update and Galerkin integrator (TT-BUG) advances passive and active bursts with approximately half a billion unknowns. Selected runs are 15-20 times faster than AMEn at comparable verified scalar-flux errors. In the reported anisotropic active-burst configurations, TT-BUG takes 14.3 s on one processor versus 716 s for OpenSn on 24 processors: approximately 50 times less wall time and 1200 times less allocated processor-time cost. This compares reported configurations rather than matched physical errors. OpenSn provides an independent discretization check. The modified TT-BUG source-onset and saturation treatment is tested numerically; its order, convergence and general tolerance-to-error relation remain unproved.

In the iron oxide Fe$_3$BO$_6$ weak ferromagnetism stems from competing contributions from two-site Dzyaloshinskii-Moriya interactions (DMI) and from on-site anisotropies at the two $4c$ and $8d$ types of nonequivalent Fe$^{3+}$ sites. As a result, spin reorientation (SR) transition $Γ_2$ $\rightarrow$ $Γ_4$ in Fe$_3$BO$_6$ uniquely comprises 90$^\circ$ spin switching and a 35% change of weak ferromagnetic moment. We explore ultrafast kinetics of such a transition as a proxy for weak ferromagnetism contributions modified by femtosecond laser pulses. Using magneto-optical pump-probe technique we reveal that the SR transition is triggered by an abrupt decrease of the weak ferromagnetic moment in the $Γ_2$ phase, and relate it to the laser-induced modification of the DMI and anisotropy contributions to the spin canting. By tuning the pump pulse phonon energy between different optical absorption ranges, we show that laser-driven impact on spin canting possesses a nonthermal character and is achieved through the charge-transfer transitions $2p\rightarrow3d$ within Fe$^{3+}$O$^{2-}_6$ octahedra rather than by localized $3d\rightarrow3d$ transitions at Fe$^{3+}$. The decreased weak ferromagnetic moment controls the threshold and the latency of the SR transition, as well as a final state. The latter is achieved via the precessional switching to the $Γ_4$ phase activated by the laser-induced heating. Final state is found to be a coexistence of $Γ_4$ and $Γ_2$ phase with modified spin canting persisting even at hundreds of picoseconds.

Quantum-inspired numerical methods for partial differential equations represent discretized fields and operators as quantics tensor trains (QTT), so that their cost depends on how compressible the problem is rather than on the number of unknowns. Global spectral formulations in this setting have so far relied on Fourier bases, restricted to periodic problems, or on collocation, which does not preserve the variational structure exploited by tensor-network solvers. We introduce a QTT Chebyshev--Galerkin method that assembles and solves global spectral systems directly in compressed form, without forming the underlying arrays. Starting from the weak formulation, the method handles coupled elliptic systems with variable coefficients, smoothly mapped geometries, and mixed Dirichlet--Robin boundary conditions, and yields symmetric positive-definite systems suited to variational tensor-network solvers. Parameters and random inputs enter as additional coordinates, so the same construction returns the solution over an entire parameter space in a single solve. In linear elasticity, the method retains spectral convergence and reaches errors below $10^{-8}$ on a three-dimensional benchmark, beyond those attained by QTT finite differences and Q2/Q4 FEniCSx finite elements over the tested range. A six-parameter cantilever design study covering $1.7\times10^{7}$ designs is solved in a single 92~s run, and a stochastic problem with ten random coordinates is represented with $2.4\times10^{6}$ times less memory than its full coefficient array.

The free energy of a crystal is an integral over the positions of all its atoms, in 3N dimensions. We evaluate it without random sampling for the ordered alloys B2 NiAl and gamma-TiAl with embedded-atom potentials, using a self-consistent harmonic approximation (SCHA) corrected to third order in temperature, calculated via thermal averages of tensor-train (TT) integrals built from one atom and its neighbors, which we call TT-SCHA. The result is for an infinite crystal without sampling error; its numerical error is set by the starting pivots of the interpolation (up to 9 MPa in the pressure, reduced by averaging three runs; SM Table S26). Against the molecular dynamics of 65,536 atoms (Frenkel-Ladd free energy at 600 K, integrated energy above), TT-SCHA free energy agrees to 0.13 meV per atom at 600-1000 K, using 420-500 times fewer core-seconds for one run. At zero pressure, results agree with published calculations that use the same embedded-atom potential: the published elastic constants of gamma-TiAl (to 1 GPa, by a static lattice calculation) and its thermal expansion from Monte Carlo sampling (to 0.012 percentage points).