We provide a simple algorithm for constructing Hamiltonian graph cycles (visiting every vertex exactly once) on a set of arbitrarily large finite subgraphs of aperiodic two-dimensional Ammann-Beenker (AB) tilings. Using this result, and the discrete scale symmetry of AB tilings, we find exact solutions to a range of other problems which lie in the complexity class NP-complete for general graphs. These include the equal-weight traveling salesperson problem, providing, for example, the most efficient route a scanning tunneling microscope tip could take to image the atoms of physical quasicrystals with AB symmetries; the longest path problem, whose solution demonstrates that collections of flexible molecules of any length can adsorb onto AB quasicrystal surfaces at density one, with possible applications to catalysis; and the three-coloring problem, giving ground states for the $q$-state Potts model ($q\ge 3$) of magnetic interactions defined on the planar dual to AB, which may provide useful models for protein folding.

In thermal equilibrium the dynamics of phase transitions is largely controlled by fluctuation-dissipation relations: On the one hand, friction suppresses fluctuations, while on the other hand, the thermal noise is proportional to friction constants. Out of equilibrium, this balance dissolves and one can have situations where friction vanishes due to antidamping in the presence of a finite noise level. We study a wide class of $\mathrm{O}(N)$ field theories where this situation is realized at a phase transition, which we identify as a critical exceptional point. In the ordered phase, antidamping induces a continuous limit cycle rotation of the order parameter with an enhanced number of $2N-3$ Goldstone modes. Close to the critical exceptional point, however, fluctuations diverge so strongly due to the suppression of friction that in dimensions $d<4$ they universally either destroy a preexisting static order or give rise to a fluctuation-induced first-order transition. This is demonstrated within a full resummation of loop corrections via Dyson-Schwinger equations for $N=2$, and a generalization for arbitrary $N$, which can be solved in the long wavelength limit. We show that in order to realize this physics it is not necessary to drive a system far out of equilibrium: Using the peculiar protection of Goldstone modes, the transition from an $xy$ magnet to a ferrimagnet is governed by an exceptional critical point once weakly perturbed away from thermal equilibrium.

Two-dimensional (2D) transition metal dichalcogenides (TMDs) is a versatile class of quantum materials of interest to various fields including, e.g., nanoelectronics, optical devices, and topological and correlated quantum matter. Tailoring the electronic properties of TMDs is essential to their applications in many directions. Here, we report that a highly controllable and uniform on-chip 2D metallization process converts a class of atomically thin TMDs into robust superconductors, a property belonging to none of the starting materials. As examples, we demonstrate the introduction of superconductivity into a class of 2D air-sensitive topological TMDs, including monolayers of ${\mathrm{T}}_{\mathrm{d}}\text{−}{\mathrm{WTe}}_2$, $1{\mathrm{T}}^{\prime }\text{−}{\mathrm{MoTe}}_2$, and $2\mathrm{H}\text{−}{\mathrm{MoTe}}_2$, as well as their natural and twisted bilayers, metallized with an ultrathin layer of palladium. This class of TMDs is known to exhibit intriguing topological phases ranging from topological insulator, Weyl semimetal to fractional Chern insulator. The unique, high-quality two-dimensional metallization process is based on our recent findings of the long-distance, non-Fickian in-plane mass transport and chemistry in 2D that occur at relatively low temperatures and in devices fully encapsulated with inert insulating layers. Highly compatible with existing nanofabrication techniques for van der Waals stacks, our results offer a route to designing and engineering superconductivity and topological phases in a class of correlated 2D materials.

A topological “Thouless” pump represents the quantized motion of particles in response to a slow, cyclic modulation of external control parameters. The Thouless pump, like the quantum Hall effect, is of fundamental interest in physics, because it links physically measurable quantities, such as particle currents, to geometric properties of the experimental system, which can be robust against perturbations and, thus, technologically useful. So far, experiments probing the interplay between topology and interparticle interactions have remained relatively scarce. Here, we observe a Thouless-type charge pump in which the particle current and its directionality inherently rely on the presence of strong interactions. Experimentally, we utilize a two-component Fermi gas in a dynamical superlattice which does not exhibit a sliding motion and remains trivial in the single-particle regime. However, when tuning interparticle interactions from zero to positive values, the system undergoes a transition from being stationary to drifting in one direction, consistent with quantized pumping in the first cycle. Remarkably, the topology of the interacting pump trajectory cannot be adiabatically connected to a noninteracting limit, highlighted by the fact that only one atom is transferred per cycle. Our experiments suggest that Thouless charge pumps are promising platforms to gain insights into interaction-driven topological transitions and topological quantum matter.

Quantum many-body systems, particularly, the ones with large near-$E_F$ density states, are well known for exhibiting rich phase diagrams as a result of enhanced electron correlations. The recently discovered locally noncentrosymmetric heavy fermion superconductor ${\mathrm{CeRh}}_2{\mathrm{As}}_2$ has stimulated extensive attention due to its unusual $H\text{−}T$ phase diagram consisting of two-phase superconductivity, antiferromagnetic order, and possible quadrupole-density wave orders. However, the critical near-$E_F$ electronic structure remains experimentally elusive. Here, we provide this key information by combining soft-x-ray and vacuum ultraviolet (VUV) angle-resolved-photoemission-spectroscopy measurements and atom-resolved density-functional-theory $(\mathrm{DFT})+U$ calculations. With bulk-sensitive soft x ray, we reveal quasi-2D hole and electron pockets near the $E_F$. On the other hand, under VUV light, the Ce flat bands are resolved with the $c\text{−}f$ hybridization persisting up to well above the Kondo temperature. Most importantly, we observe a symmetry-protected fourfold Van Hove singularity (VHS) coexisting with the $\mathrm{Ce}4f_{5/2}^1$ flat bands at the $X$ point, which, to the best of our knowledge, has never been reported before. Such a rare coexistence is expected to lead to a large density of states at the zone edge, a large upper critical field of the odd-parity phase, as well as spin and/or charge instabilities with a vector of ($1/2$, $1/2$, 0). Uniquely, it will also result in a new type of $f$-VHS hybridization that alters the order and fine electronic structure of the VHS and flat bands. Our findings provide not only key insights into the nature of multiple phases in ${\mathrm{CeRh}}_2{\mathrm{As}}_2$ but also open up new prospects for exploring the novelties of many-body systems with $f$-VHS hybridization.

We propose a new electrical breakdown mechanism for exciton insulators in the BCS limit, which differs fundamentally from the Zener breakdown mechanism observed in traditional band insulators. Our new mechanism results from the instability of the many-body ground state for exciton condensation, caused by the strong competition between the polarization and condensation energies in the presence of an electric field. We refer to this mechanism as “many-body breakdown.” To investigate this new mechanism, we propose a BCS-type trial wave function under finite electric fields and use it to study the many-body breakdown numerically. Our results reveal two different types of electric breakdown behavior. If the system size is larger than a critical value, the Zener tunneling process is first turned on when an electrical field is applied, but the excitonic gap remains until the field strength reaches the critical value of the many-body breakdown, after which the excitonic gap disappears and the system becomes a highly conductive metallic state. However, if the system size is much smaller than the critical value, the intermediate tunneling phase disappears since the many-body breakdown happens before the onset of Zener tunneling. The sudden disappearance of the local gap leads to an “off-on” feature in the current-voltage ($I\text{−}V$) curve, providing a straightforward way to distinguish excitonic insulators from normal insulators.

Berry curvature multipoles appearing in topological quantum materials have recently attracted much attention. Their presence can manifest in novel phenomena, such as nonlinear anomalous Hall effects (NLAHE). The notion of Berry curvature multipoles extends our understanding of Berry curvature effects on the material properties. Hence, research on this subject is of fundamental importance and may also enable future applications in energy harvesting and high-frequency technology. It was shown that a Berry curvature dipole can give rise to a second-order NLAHE in materials of low crystalline symmetry. Here, we demonstrate a fundamentally new mechanism for Berry curvature multipoles in antiferromagnets that are supported by the underlying magnetic symmetries. Carrying out electric transport measurements on the kagome antiferromagnet FeSn, we observe a third-order NLAHE, which appears as a transverse voltage response at the third harmonic frequency when a longitudinal ac drive is applied. Interestingly, this NLAHE is strongest at and above room temperature. We combine these measurements with a scaling law analysis, a symmetry analysis, model calculations, first-principle calculations, and magnetic Monte Carlo simulations to show that the observed NLAHE is induced by a Berry curvature quadrupole appearing in the spin-canted state of FeSn. At a practical level, our study establishes NLAHE as a sensitive probe of antiferromagnetic phase transitions in other materials—such as moiré superlattices, two-dimensional van der Waal magnets, and quantum spin liquid candidates, which remain poorly understood to date. More broadly, Berry curvature multipole effects are predicted to exist for 90 magnetic point groups. Hence, our work opens a new research area to study a variety of topological magnetic materials through nonlinear measurement protocols.

Spin current is a crucial element in spintronics, and its diffusion in materials is typically characterized by monotonic decay. However, when the material hosting the spin current is also a magnet, the spin current is expected to exhibit spatial oscillations, the observation of which remains elusive. Here, we demonstrate the spatial oscillation of a spin current in a nickel film by measuring the thickness-dependent inverse spin Hall effect in $\mathrm{Ni}/\mathrm{YIG}$ bilayers. The inverse spin Hall current in nickel is found to oscillate with its film thickness, in contrast to nonmagnetic materials, and that the oscillation period quantitatively agrees with theoretical predictions based on differences in the Fermi wave vector between majority and minority carriers. Our findings reveal a previously hidden behavior of spin-transport dynamics and identify a new degree of freedom for manipulating spin current, with potential implications for future spintronic devices.

The deconfined quantum critical point (DQCP) is an example of phase transitions beyond the Landau symmetry-breaking paradigm that attracts wide interest. However, its nature has not been settled after decades of study. In this paper, we apply the recently proposed fuzzy-sphere regularization to study the SO(5) nonlinear sigma model with a topological Wess-Zumino-Witten term, which serves as a dual description of the DQCP with an exact SO(5) symmetry. We demonstrate that the fuzzy sphere functions as a powerful microscope, magnifying and revealing a wealth of crucial information about the DQCP, ultimately paving the way toward its final answer. In particular, through exact diagonalization, we provide clear evidence that the DQCP exhibits approximate conformal symmetry. The evidence includes the existence of a conserved SO(5) symmetry current, a stress tensor, and integer-spaced levels between conformal primaries and their descendants. Most remarkably, we identify 23 primaries and 76 conformal descendants. Furthermore, by examining the renormalization group flow of the lowest symmetry singlet as well as other primaries, we provide numerical evidence in favor of DQCP being pseudo-critical, with the approximate conformal symmetry plausibly emerging from nearby complex fixed points. The primary spectrum we compute also has important implications, including the conclusion that the SO(5) DQCP cannot describe a direct transition from the Néel to valence bond solid phase on the honeycomb lattice.

${\mathrm{FeTe}}_{0.55}{\mathrm{Se}}_{0.45}$ (FTS) occupies a special spot in modern condensed matter physics at the intersections of electron correlation, topology, and unconventional superconductivity. The bulk electronic structure of FTS is predicted to be topologically nontrivial due to the band inversion between the $d_{xz}$ and $p_z$ bands along $\mathrm{\Gamma }\text{−}Z$. However, there remain debates in both the authenticity of the Dirac surface states (DSSs) and the experimental deviations of band structure from the theoretical band inversion picture. Here we resolve these debates through a comprehensive angle-resolved photoemission spectroscopy investigation. We first observe a persistent DSS independent of $k_z$. Then, by comparing FTS with FeSe, which has no band inversion along $\mathrm{\Gamma }\text{−}Z$, we identify the spectral weight fingerprint of both the presence of the $p_z$ band and the inversion between the $d_{xz}$ and $p_z$ bands. Furthermore, we propose a renormalization scheme for the band structure under the framework of a tight-binding model preserving crystal symmetry. Our results highlight the significant influence of correlation on modifying the band structure and make a strong case for the existence of topological band structure in this unconventional superconductor.

Magic-angle twisted bilayer graphene is the best-studied physical platform featuring moiré potential-induced narrow bands with nontrivial topology and strong electronic correlations. Despite their significance, the Chern insulating states observed at a finite magnetic field—and extrapolating to a band filling $s$ at zero field—remain poorly understood. Unraveling their nature is among the most important open problems in the province of moiré materials. Here, we present the first comprehensive study of interacting electrons in finite magnetic field while varying the electron density, twist angle, and heterostrain. Within a panoply of correlated Chern phases emerging at a range of twist angles, we uncover a unified description for the ubiquitous sequence of states with the Chern number $t$ for $(s,t)=\pm (0,4)$, $\pm (1,3)$, $\pm (2,2)$, and $\pm (3,1)$. We also find correlated Chern insulators at unconventional sequences with $s+t\ne \pm 4$, as well as with fractional $s$, and elucidate their nature.

A fundamental distinction between many-body quantum states are those with short- and long-range entanglement (SRE and LRE). The latter cannot be created by finite-depth circuits, underscoring the nonlocal nature of Schrödinger cat states, topological order, and quantum criticality. Remarkably, examples are known where LRE is obtained by performing single-site measurements on SRE, such as the toric code from measuring a sublattice of a 2D cluster state. However, a systematic understanding of when and how measurements of SRE give rise to LRE is still lacking. Here, we establish that LRE appears upon performing measurements on symmetry-protected topological (SPT) phases—of which the cluster state is one example. For instance, we show how to implement the Kramers-Wannier transformation by adding a cluster SPT to an input state followed by measurement. This transformation naturally relates states with SRE and LRE. An application is the realization of double-semion order when the input state is the ${\mathbb{Z}}_2$ Levin-Gu SPT. Similarly, the addition of fermionic SPTs and measurement leads to an implementation of the Jordan-Wigner transformation of a general state. More generally, we argue that a large class of SPT phases protected by $G\times H$ symmetry gives rise to anomalous LRE upon measuring $G$-charges, and we prove that this persists for generic points in the SPT phase under certain conditions. Our work introduces a new practical tool for using SPT phases as resources for creating LRE, and we uncover the classification result that all states related by sequentially gauging Abelian groups or by Jordan-Wigner transformation are in the same equivalence class, once we augment finite-depth circuits with single-site measurements. In particular, any topological or fracton order with a solvable finite gauge group can be obtained from a product state in this way.

The value of fundamental physical constants is affected by the coupling of matter to the electromagnetic vacuum state, as predicted and explained by quantum electrodynamics. In this work, we present a millikelvin magnetotransport experiment in the quantum Hall regime that assesses the possibility of the von Klitzing constant being modified by strong cavity vacuum fields. By employing a Wheatstone bridge, we measure the difference between the quantized Hall resistance of a cavity-embedded Hall bar and the resistance standard, achieving an accuracy down to one part in ${10}^5$ for the lowest Landau level. While our results do not suggest any deviation that could imply a modified Hall resistance, our work represents pioneering efforts in exploring the fundamental implications of vacuum fields in solid-state systems.

The atomic cluster expansion provides local, complete basis functions that enable efficient parametrization of many-atom interactions. We extend the atomic cluster expansion to incorporate graph basis functions. This naturally leads to representations that enable the efficient description of semilocal interactions in physically and chemically transparent form. Simplification of the graph expansion by tensor decomposition results in an iterative procedure that comprises current message-passing machine learning interatomic potentials. We demonstrate the accuracy and efficiency of the graph atomic cluster expansion for a number of small molecules, clusters, and a general-purpose model for carbon. We further show that the graph atomic cluster expansion scales linearly with the number of neighbors and layer depth of the graph basis functions.

We use symmetry analysis and density-functional theory to determine and characterize surface terminations that have a finite equilibrium magnetization density in antiferromagnetic materials. A nonzero magnetic dipole moment per unit area or “surface magnetization” can arise on particular surfaces of many antiferromagnets due to the bulk magnetic symmetries. Such surface magnetization underlies intriguing physical phenomena like interfacial magnetic coupling and can be used as a readout method of antiferromagnetic domains. However, a universal description of antiferromagnetic surface magnetization is lacking. We first introduce a classification system based on whether the surface magnetization is either sensitive or robust to roughness and on whether the magnetic dipoles at surface of interest are compensated or uncompensated when the bulk magnetic order is retained at the surface. We show that roughness-sensitive categories can be identified by a simple extension of a previously established group-theory formalism for identifying roughness-robust surface magnetization. We then map the group-theory method of identifying surface magnetization to a novel description in terms of bulk magnetic multipoles, which are already established as symmetry indicators for bulk magnetoelectric responses at both linear and higher orders. We use density-functional calculations to illustrate that nominally compensated surfaces in magnetoelectric ${\mathrm{Cr}}_2{\mathrm{O}}_3$ and centrosymmetric altermagnetic ${\mathrm{FeF}}_2$ develop a finite magnetization density at the surface, in agreement with our predictions based on both group theory and the ordering of the bulk multipoles. Our analysis provides a comprehensive basis for understanding the surface magnetic properties and their intimate correspondence to bulk magnetoelectric effects in antiferromagnets and has important implications for technologically relevant phenomena such as exchange-bias coupling.

Understanding superfluidity remains a major goal of condensed matter physics. Here, we tackle this challenge utilizing the recently developed fermionic neural network (FermiNet) wave function Ansatz [D. Pfau et al., Phys. Rev. Res. 2, 033429 (2020).] for variational Monte Carlo calculations. We study the unitary Fermi gas, a system with strong, short-range, two-body interactions known to possess a superfluid ground state but difficult to describe quantitatively. We demonstrate key limitations of the FermiNet Ansatz in studying the unitary Fermi gas and propose a simple modification based on the idea of an antisymmetric geminal power singlet (AGPs) wave function. The new AGPs FermiNet outperforms the original FermiNet significantly in paired systems, giving results which are more accurate than fixed-node diffusion Monte Carlo and are consistent with experiment. We prove mathematically that the new Ansatz, which differs from the original Ansatz only by the method of antisymmetrization, is a strict generalization of the original FermiNet architecture, despite the use of fewer parameters. Our approach shares several advantages with the original FermiNet: The use of a neural network removes the need for an underlying basis set; sand the flexibility of the network yields extremely accurate results within a variational quantum Monte Carlo framework that provides access to unbiased estimates of arbitrary ground-state expectation values. We discuss how the method can be extended to study other superfluid.

We derive an exact solution for the steady state of a setup where two $XX$-coupled $N$-qubit spin chains (with possibly nonuniform couplings) are subject to boundary Rabi drives and common boundary loss generated by a waveguide (either bidirectional or unidirectional). For a wide range of parameters, this system has a pure entangled steady state, providing a means for stabilizing remote multiqubit entanglement without the use of squeezed light. Our solution also provides insights into a single boundary-driven dissipative $XX$ spin chain that maps to an interacting fermionic model. The nonequilibrium steady state exhibits surprising correlation effects, including an emergent pairing of hole excitations that arises from dynamically constrained hopping. Our system could be implemented in a number of experimental platforms, including circuit QED.

The flux quantization is a key indication of electron pairing in superconductors. For example, the well-known $h/2e$ flux quantization is considered strong evidence for the existence of charge-$2e$, two-electron Cooper pairs. Here we report evidence for multicharge flux quantization in mesoscopic ring devices fabricated using the transition-metal kagome superconductor ${\mathrm{CsV}}_3{\mathrm{Sb}}_5$. We perform systematic magnetotransport measurements and observe unprecedented quantization of magnetic flux in units of $h/4e$ and $h/6e$ in magnetoresistance oscillations. Specifically, at low temperatures, magnetoresistance oscillations with period $h/2e$ are detected, as expected from the flux quantization for charge-$2e$ superconductivity. We find that the $h/2e$ oscillations are suppressed and replaced by resistance oscillations with $h/4e$ periodicity when the temperature is increased. Increasing the temperature further suppresses the $h/4e$ oscillations, and robust resistance oscillations with $h/6e$ periodicity emerge as evidence for charge-$6e$ flux quantization. Our observations provide the first experimental evidence for the existence of multicharge flux quanta and emergent quantum matter exhibiting higher-charge superconductivity in the strongly fluctuating region above the charge-$2e$ Cooper pair condensate, revealing new insights into the intertwined and vestigial electronic order in kagome superconductors.

First-principles calculations of electron interactions in materials have seen rapid progress in recent years, with electron-phonon ($e\text{−}\mathrm{ph}$) interactions being a prime example. However, these techniques use large matrices encoding the interactions on dense momentum grids, which reduces computational efficiency and obscures interpretability. For $e\text{−}\mathrm{ph}$ interactions, existing interpolation techniques leverage locality in real space, but the high dimensionality of the data remains a bottleneck to balance cost and accuracy. Here we show an efficient way to compress $e\text{−}\mathrm{ph}$ interactions based on singular value decomposition (SVD), a widely used matrix and image compression technique. Leveraging (un)constrained SVD methods, we accurately predict material properties related to $e\text{−}\mathrm{ph}$ interactions—including charge mobility, spin relaxation times, band renormalization, and superconducting critical temperature—while using only a small fraction (1%–2%) of the interaction data. These findings unveil the hidden low-dimensional nature of $e\text{−}\mathrm{ph}$ interactions. Furthermore, they accelerate state-of-the-art first-principles $e\text{−}\mathrm{ph}$ calculations by about 2 orders of magnitude without sacrificing accuracy. Our Pareto-optimal parametrization of $e\text{−}\mathrm{ph}$ interactions can be readily generalized to electron-electron and electron-defect interactions, as well as to other couplings, advancing quantitative studies of condensed matter.

A direct current through a metal-insulator-metal tunneling junction emits light when surface-plasmon polaritons (SPPs) are excited. Two distinct processes are believed to coexist in this light emission mediated by surface plasmons: inelastic tunneling, where electrons excite SPPs in the insulator gap, and hot-electron radiative decay, which occurs in the electrodes after elastic tunneling. Previous theoretical approaches to study light emission by inelastic tunneling have relied on Bardeen’s approximation where the electronic wave functions are considered only in the barrier of the junction. In this work, we introduce an extension to models of inelastic tunneling by incorporating the full quantum device solution of the Schrödinger equation, which can also account for processes in the metallic electrodes. The extension unveils the existence of long-range correlations of the current density across the barrier and enables us to establish the equivalence between two models widely used in the past: (i) a calculation of the inelastic transition rate between two states across the barrier based on Fermi’s golden rule and (ii) a calculation of the power transferred to plasmons by current fluctuations. Importantly, the new model accounts for processes that take place in the metallic electrodes and that could not be described within Bardeen’s approximation. Hence, it is no longer necessary to invoke a hot-electron mechanism to obtain a dependence on the geometry of metallic electrodes. The new framework enables to discuss the role of surface plasmons localized in different metal-insulator interfaces and to include possible nonlocal effects at the interfaces.