Excitons, the correlated electron-hole pairs governing optical and transport properties in organic semiconductors, have long resisted direct experimental access to their full quantum-mechanical wave functions. Here, we use femtosecond time-resolved photoemission orbital tomography (trPOT) combining high-harmonic probe pulses with time- and momentum-resolved photoelectron spectroscopy to directly image the momentum-space distribution and ultrafast dynamics of excitons in $\alpha$-sexithiophene thin films. We introduce a model that enables reconstruction of the exciton wave function in real space, including both its spatial extent and its internal phase structure. The reconstructed wave function reveals coherent delocalization across approximately three molecular units and exhibits a characteristic phase modulation, consistent with ab initio calculations within the framework of many-body perturbation theory. Time-resolved measurements further indicate an approximately 25% contraction of the exciton radius within 400 fs, suggesting self-trapping driven by exciton-phonon coupling. These results establish trPOT as a general and experimentally accessible approach for resolving exciton wave functions—with spatial, phase, and temporal sensitivity—in a broad class of molecular and low-dimensional materials.

Preparation of quantum thermal states of many-body systems is a key computational challenge for quantum processors, with applications in physics, chemistry, and classical optimization. We provide a simple and efficient algorithm for thermal state preparation, combining engineered bath resetting and modulated system-bath coupling to derive a quantum channel approximately satisfying quantum detailed balance relations. We show that the fixed point $\widehat{\sigma }$ of the channel approximates the Gibbs state as $\parallel \widehat{\sigma }-{\widehat{\sigma }}_{\beta }\parallel \sim {\theta }^2$, where $\theta$ is the system-bath coupling and ${\widehat{\sigma }}_{\beta }\propto e^{-\beta {\widehat{H}}_S}$. We provide extensive numerics, for the example of the 2D quantum Ising model, confirming that the protocol successfully prepares the thermal state throughout the finite-temperature phase diagram, including near the quantum phase transition. Simulations for free-fermion systems provide further evidence for the accuracy of the protocol for large system sizes in the weak-coupling limit. Our algorithm provides a path to efficient quantum simulation of quantum-correlated states at finite temperature with current and near-term quantum processors.

Quantifying the irreversibility and dissipation of nonequilibrium processes is crucial to understanding their behavior, assessing their possible capabilities, and characterizing their efficiency. We introduce a physical quantity that quantifies the irreversibility of stochastic Langevin systems from the observation of individual molecules’ displacements. Categorizing these displacements into a few groups based on their initial and final positions allows us to measure irreversibility precisely without the need to know the forces and magnitude of the fluctuations acting on the system. For short times, our model-free estimate of irreversibility is related to entropy production by a conditional fluctuation theorem. For short times and, in general, for stationary protocols, our estimate provides a lower bound to the average entropy production. We validate the method on single-molecule force spectroscopy experiments of proteins subject to force ramps. We show that irreversibility is sensitive to detailed features of the energy landscape underlying the protein folding dynamics and suggest how our methods can be employed to unveil key properties of protein folding processes.

Quantum routers (QRouters) are essential components of bucket-brigade quantum random access memory (QRAM), enabling quantum applications such as Grover’s search and quantum machine learning. Despite significant theoretical advances, achieving scalable and coherent QRouters experimentally remains challenging. Here, we demonstrate coherent quantum routers using a superconducting quantum processor, laying a practical foundation for scalable QRAM systems. The quantum router at the core of our implementation utilizes the transition composite gate scheme, wherein auxiliary energy levels temporarily mediate conditional interactions. By leveraging two-qubit primitives in its construction, the quantum router preserves straightforward initialization and calibration while significantly reducing circuit depth and duration relative to traditional gate-based decompositions. Moreover, by encoding routing addresses in the nonadjacent qutrit states $|0⟩$ and $|2⟩$, our design inherently enables erasure-detection capability, providing efficient postselection to mitigate routing errors. Experimentally, we achieve an average fidelity of 94.8% across three individual QRouters, and we validate scalability through a two-layer quantum routing network achieving a fidelity of 82.4%. Our results represent a significant advancement in quantum routing technology, providing enhanced fidelity, built-in error resilience, and practical scalability crucial for the development of future QRAM and large-scale quantum computing architectures.

Bounding the correlations predicted by quantum theory is an important challenge in quantum information science. Today’s leading approach is semidefinite programming relaxations, but existing methods still cannot account for many relevant types of constraints. Here, we propose a general semidefinite relaxation methodology that can incorporate a breadth of constraints needed in various quantum correlation problems, thereby generalising the seminal Navascués-Pironio-Acín hierarchy. It yields useful results at reasonable computational cost. We showcase the methodology and its features by using it to address five different quantum information problems. These are (i) entanglement witnessing from imperfect measurement devices, (ii) certifying measurements from fidelity-constrained sources, (iii) computing dimensionality in genuine multiparticle entangled states, (iv) benchmarking dimensionality for state preparation devices, and (v) finding uncertainty relations for nearly anticommuting observables. These applications reflect both the usefulness and versatility of the methodology, as well as its potential for broader relevance in the field.

Long-term stability in materials is commonly associated with resistance to external perturbations. Yet metastable defect populations may relax spontaneously even in the absence of environmental stimuli, reshaping macroscopic properties over extended timescales. Here, we show the real-space visualization of intrinsic aging of $n$-type ${\mathrm{Mg}}_3{(\mathrm{Sb},\mathrm{Bi})}_2$ under inert and room-temperature conditions after two-year storage by atomic-scale characterizations. Density-functional-theory calculations and machine-learning molecular dynamics simulations prove that the aging originates from relaxation of a metastable Mg-rich state, in which low vacancy-formation energies and migration barriers enable thermodynamically favorable and kinetically accessible Mg redistribution. Multiscale characterizations further reveal that grain boundaries (GBs) act as fast-diffusion pathways and effective sinks, establishing a hierarchical Mg redistribution process from grain interiors to GBs and surfaces. These results identify intrinsic defect physics as one of the fundamental stability constraints in functional materials and underscore the importance of controlling defect thermodynamics and transport pathways for achieving practically durable materials.

In the study of electronic structure, the wave function view dominates the current research landscape and forms the theoretical foundation of modern quantum mechanics. In contrast, valence bond (VB) theory represents chemical bonds as shared electron pairs and can provide an intuitive, particle-based insight into chemical bonding. In this work, using a newly developed periodic dynamic Voronoi Metropolis sampling (PDVMS) method, we project classical many-body electronic configurations from the neural network wave function and apply VB theory to construct a complementary particle-view paradigm. The powerful neural network wave function can achieve near-exact ab initio solutions for the ground state of both molecular and solid systems. It allows us to definitively characterize the ground state of benzene by reassessing the competition between its two VB structures. Extending PDVMS to solids, we also predict a spin-staggered VB structure in graphene, which explains emergent magnetic properties in graphene-based nanostructures. Furthermore, this particle view provides insight into the optimization process of the neural network wave function itself. Our work thus introduces a framework for analyzing many-body electronic structure in molecules and solids, opening new avenues for investigating this complex problem and its associated exotic phenomena.

The generation of non-Gaussian quantum states is a key requirement for universal continuous-variable quantum information processing. We report the experimental generation of large-amplitude squeezed coherent-state superpositions (squeezed cat states) on free-space optical pulses, reaching an amplitude of $\alpha =2.47$, which, to our knowledge, exceeds all previously reported values. Our protocol relies on the controlled mixing of the Fock states $|1⟩$ and $|2⟩$ through a tunable beam splitter, followed by heralding via homodyne detection. The resulting state displays three well-resolved negative regions in its Wigner function and achieves a fidelity of 0.53 with the target state $\propto \widehat{S}(z)(|\alpha ⟩-|-\alpha ⟩)$, with $\alpha =2.47$ and squeezing parameter $z=0.56$. These results constitute a significant milestone for temporal breeding protocols and for the iterative generation of optical Gottesman-Kitaev-Preskill states, opening perspectives for scalable and fault-tolerant photonic quantum architectures.

Superionic conductors are fascinating systems representing a state of matter intermediate between solids and liquids, yet their underlying physics remain poorly understood. Correlations in jump dynamics have been theoretically proposed as a key mechanism underlying superionic diffusion, but a direct experimental demonstration remains elusive. Here, we report the direct observation of collective cation hopping in the argyrodite compound ${\mathrm{Cu}}_7{\mathrm{PS}}_6$, using single-crystal neutron scattering to access the full four-dimensional space-time correlations of ionic motion. By leveraging the rarely used coherent scattering channel in quasielastic neutron scattering, together with diffuse scattering, we reveal pronounced short-range dynamical correlations between diffusing Cu ions, providing direct experimental evidence of collective hopping in a superionic phase. Machine-learned molecular dynamics simulations quantitatively reproduce all our measurements and pinpoint the dominant collective hopping mechanism. Furthermore, we establish through both experiments and simulations that the collective diffusive hops originate from the breakdown of low-energy phonons dominated by the Cu motions in the parent ordered phase, demonstrating a key link between phonon precursors and the activation of facile collective hops. This work presents the first direct experimental evidence of collective hopping in superionic argyrodites and identifies wave-vector-resolved soft phonons as its physical origin, advancing the microscopic understanding of ionic transport in superionic solids.

Divalent atoms have emerged as powerful alternatives to alkalis in ultracold atom platforms, offering unique advantages arising from their two-electron structure. Among these species, ytterbium (Yb) is especially promising, yet its anionic properties and its Rydberg spectrum remain comparatively unexplored. In this work, we perform a comprehensive experimental and theoretical investigation of ultralong-range Rydberg molecules (ULRMs) of ${}^{174}\mathrm{Yb}$ in $6sns{{}^{1}S}_0$ Rydberg states across nearly two decades in principal quantum number $n$ and 3 orders of magnitude in molecular binding energy. Using the Coulomb Green’s function formalism, we compute Born-Oppenheimer molecular potentials describing the Rydberg atom in the presence of a ground-state perturber and achieve quantitative agreement with high-resolution molecular spectra. This enables the extraction of low-energy electron-Yb scattering phase shifts, including the zero-energy $s$-wave scattering length and the positions of two spin-orbit split $p$-wave shape resonances. Our results provide strong evidence that the ${\mathrm{Yb}}^{-}$ anion exists only as a metastable resonance. We additionally show the sensitivity of ULRM spectra to the atomic quantum defects, using this to determine the quantum defect of the $6s23f{{}^{1}F}_3$ state. Together, these findings establish Yb ULRMs as a powerful probe of electron-Yb interactions and lay essential groundwork for future Rydberg experiments with divalent atoms.

Ab initio auxiliary-field quantum Monte Carlo (AFQMC) is a systematically improvable many-body method, but its application to extended solids has been severely limited by unfavorable computational scaling and memory requirements that obstruct rigorous access to the thermodynamic and complete-basis-set limits. By combining tensor hypercontraction with $\mathbf{k}$-point symmetry, we reduce the computational and memory scaling of ab initio AFQMC for solids to $\mathcal{O}(N^3)$ and $\mathcal{O}(N^2)$, respectively, with an arbitrary basis, comparable to diffusion Monte Carlo. This enables simultaneous, controlled extrapolation of full-system AFQMC calculations to the thermodynamic and complete-basis-set limits across insulating, metallic, and strongly correlated solids, without embedding, local approximations, empirical finite-size corrections, or composite schemes. Our results establish AFQMC as a general-purpose, systematically improvable alternative to diffusion Monte Carlo and coupled-cluster methods for predictive ab initio simulations of solids, enabling accurate energies and magnetic observables within a unified framework.

The appearance of fractional Chern insulators in moiré systems can be rationalized by the presence of a fictitious magnetic field associated with the spatial texture of layer-resolved electronic wave functions. Here, we present a systematic study of real-space topology and the associated fictitious magnetic fields in moiré systems. We first show that at the level of individual Bloch wave functions, the real-space Chern number, akin to a Pontryagin index, is a fragile marker. It generically vanishes except for specific limits where the Bloch functions exhibit fine-tuned zeros within the unit cell, such as the chiral limit of twisted bilayer graphene (TBG) or the adiabatic regime of twisted homobilayer transition metal dichalcogenides (TMDs). We then show that these limitations do not apply to textures associated with ensembles of Bloch wave functions, such as entire bands or the ensemble of states at a given energy. The Chern number of these textures defines a robust topological index protected by a spectral gap. We find that symmetries constrain it to be nonzero for both twisted TMDs and TBG across all twist angles and levels of corrugation, which can be verified in scanning tunneling microscopy measurements. By projection to the band ensemble texture, a single-component Hamiltonian under fictitious magnetic field emerges in broad regime, enabling a direct comparison of multicomponent bands to Landau-level wave functions. We also study real-space topology within the topological heavy fermion model of TBG, finding that the real-space topological features are supported only by the light $c$ electrons.

The eigenstate thermalization hypothesis (ETH) posits how isolated quantum many-body systems thermalize, assuming that individual eigenstates at the same energy density have identical expectation values of local observables in the limit of large systems. While the ETH apparently holds across a wide range of interacting quantum systems, in this work, we show that it may require generalization in the presence of thermal first-order phase transitions. We introduce a class of all-to-all spin models, featuring first-order thermal phase transitions that stem from two distinct local maxima of entropy (two mean-field solutions that we dub “branches”) that exchange dominance in the many-body density of states as the energy is varied. We argue that, for energies in the vicinity of the thermal phase transition, eigenstate expectation values do not need to converge to the same thermal value. The system has a regime with coexistence of two classes of eigenstates corresponding to the two branches with distinct expectation values at the same energy density and another regime with Schrödinger-cat-like eigenstates that are interbranch superpositions; these two regimes are separated by an eigenstate phase transition. We propose a more general form of the ETH Ansatz, support our results by semiclassical calculations and an exact diagonalization study of a microscopic spin model, and argue that the structure of eigenstates in the vicinity of thermal first-order phase transitions can be experimentally probed via nonequilibrium dynamics.

Ferroelectric domain walls are described as interfacial regions between oppositely polarized domains, where the original bulk distortions are suppressed or reconstructed at energetic cost. Here, we demonstrate that $\mathrm{Pca}{2}_1$ ${\mathrm{HfO}}_2$ exhibits a distinct ferroelectric domain wall where the polarization-reversal boundaries remain structurally bulk compatible and exceptionally low in energy. The central mechanism is a pair-condensation—the simultaneous condensation of the two trilinearly coupled phonon modes located at opposite termini of a single phonon branch. We show that this pair condensation not only establishes the microscopic origin of layer-selective ferroelectricity in bulk but also governs the exceptionally low energy of the 180° domain walls. Specifically, the simultaneous reversal and complementary energetic responses of the paired condensed modes compensate the domain-wall energy cost and facilitates a wall-centered structure where bulk motifs are locally recovered rather than suppressed. This bulklike domain wall is corroborated by the agreement between density-functional-theory calculations and annular-bright-field-STEM measurements of oxygen displacements. This logic extends to suboptimal and higher-energy domain walls which emerge as accessible structural bridges via the selective, energy-minimizing reversal of individual condensed phonon modes. By unifying these structural and energetic observations under the principle of pair condensation, this work provides a generalized predictive framework for engineering domain-wall physics in multimode functional materials.

The origin of superconductivity in magic-angle twisted bilayer graphene has been a subject of intense debate. While some experimental evidence indicated an unconventional pairing mechanism which should be sensitive to Coulomb screening, experimental attempts to tune the critical temperature by screening Coulomb interactions so far have remained unsuccessful, possibly indicating a conventional phonon-mediated pairing. Here we study a double-layer electronic system consisting of two twisted graphene bilayers in immediate proximity of each other but remaining electronically decoupled. By increasing the carrier density in one bilayer, we completely suppressed both the superconductivity and the correlated-insulator state in the adjacent magic-angle graphene. The observation of such an effect from screening offers support for an unconventional mechanism of Cooper pairing in magic-angle twisted bilayer graphene, shedding new light on the underlying physics governing their properties.

We present a general, systematic, and quantitative theory of electromagnetism and chirality in crystalline solids. Symmetry is its basic guiding principle, enabling us to consider macroscopic multipole densities without reference to any specific microscopic configurations. We use a formal analogy between space inversion $i$ and time inversion $\theta$ to identify two complementary, comprehensive classifications of crystals, based on five categories of electric and magnetic multipole order—called polarizations—and five categories of chirality. The five categories of polarizations (parapolar, electropolar, magnetopolar, antimagnetopolar, and multipolar) embody the variety of ways in which electromagnetic multipole order can be realized in solids, thus expanding the familiar notion of electric dipolarization in ferroelectrics and magnetization in ferromagnets to higher-rank multipole densities characterizing, e.g., altermagnets and alterelectrics. Group theory also permits the absolute, gauge-invariant quantification of these multipole densities. The five categories of chirality (parachiral, electrochiral, magnetochiral, antimagnetochiral, and multichiral) extend the notion of enantiomorphism—conventionally associated with the lack of spatial mirror symmetries—to include all possibilities for creating nonsuperposable images by applying the inversions $i$, $\theta$, and $i\theta$. In particular, multichiral systems lack all inversion symmetries and therefore have four different enantiomorphs. Each category of chirality is shown to arise from particular superpositions of electric and magnetic multipole densities. Our theory has great predictive power by relating electric and magnetic multipolar order with, e.g., distinctive features in the crystal structure, the electronic band structure, property tensors, and response tensors of quantum materials. These findings are critical for a large community working on fundamental and applied aspects of novel quantum materials.

Spin space groups (SSGs) impose fundamentally different constraints on magnetic configurations in real and reciprocal spaces. As a consequence, the correspondence between real-space and momentum-space spin arrangements is far richer than traditionally assumed. Building on the complete enumeration of SSGs, we develop a systematic, symmetry-based framework that classifies all possible spin arrangements allowed by these groups. This unified approach naturally incorporates conventional magnetic orders, altermagnetism, and coplanar odd-wave magnetism (including the previously proposed $p$-wave case) as distinct symmetry classes. Moreover, by determining the symmetry-allowed wave patterns of spin-polarized phases, our framework classifies coplanar odd-wave magnets beyond the $p$-wave case into higher-wave types ($f$, $h$, $j$, and $l$ wave). Crucially, our classification predicts a variety of novel magnetic phases, highlighted by the discovery of the coplanar even-wave magnet(including $d$, $g$, and $i$ wave): A state that is noncollinear in real space but hosts a collinear even-wave spin polarization in $\mathbf{k}$ space. Analysis of a minimal model reveals that this phase is characterized by nonquantized spin polarization and exhibits a novel mechanism for symmetry-enforced zero polarization on nondegenerate bands. Extending the framework from bulk crystals to layer SSGs appropriate for two-dimensional systems, we further predict layered counterparts and provide symmetry guidelines for designing bilayer coplanar odd-wave and even-wave magnets. We further validate this finding through first-principles calculations and propose ${\mathrm{CoCrO}}_4$ as a promising candidate for its experimental realization, thereby demonstrating the completeness and predictive power of the SSG-based classification of magnetic orders.

Analog quantum simulators offer a route to exploring strongly correlated many-body dynamics beyond classical computation, but their predictive power remains limited by the absence of quantitative error estimation. Establishing rigorous uncertainty bounds is essential for elevating such devices from qualitative demonstrations to quantitative scientific tools. Here we introduce a general framework for bounded-error quantum simulation, which provides predictions for many-body observables with experimentally quantifiable uncertainties. The approach combines Hamiltonian and Lindbladian Learning—a statistically rigorous inference of the coherent and dissipative generators governing the dynamics—with the propagation of their uncertainties into the simulated observables, yielding confidence bounds directly derived from experimental data. We demonstrate this framework on trapped-ion quantum simulators implementing long-range Ising interactions with up to 51 ions. We analyze error bounds on two levels. First, we learn an open-system model from experimental data collected in an initial time window of quench dynamics, simulate the corresponding master equation, and quantitatively verify consistency between theoretical predictions and measured dynamics at long times. Second, we discuss short-time error bounds directly from experimental measurements alone, without relying on classical simulation—crucial for entering regimes of quantum advantage. In both cases, the learned models reproduce the experimental evolution within the predicted bounds, demonstrating quantitative reliability and internal consistency. By integrating statistical learning, open-system modeling, and precise experimental control, bounded-error quantum simulation provides a scalable foundation for trusted analog quantum computation, bridging the gap between experimental quantum platforms and predictive many-body physics. The techniques presented here directly extend to digital quantum simulation.

We investigate how polycrystallinity evolves in homogeneously nucleated crystals in an experimental hard-sphere-like colloidal system using three-dimensional confocal microscopy. We capture the whole process of nucleation and crystal growth and use a recently developed method to resolve grain structure and misorientation between adjacent grains and nearest-neighbor clusters. Closely following nucleation, crystals are shown to give rise to grains through the emergence of stacking faults and twinning. The birth rate of new grains from a single nucleation event are measured, observed to increase, and then plateau. We visualize the spatial structure of a nucleus which comprises a highly ordered core and a radial increase in misorientation angles. Finally, we use local misorientation to formulate a measure for positional order which reveals a sharp decoupling between the development of orientational order and positional order at times significantly after nucleation.

P-type MgAgSb with excellent room-temperature thermoelectric performance holds great promise for cooling and power generation. However, its practical application has been significantly limited by phase transitions, which are long considered a polymorphic transformation where the compositions remain identical but the crystal structures differ. Here we report the non-polymorphic nature of the complex phase transitions in Mg-Ag-Sb during thermal cycling, yielding four compositionally distinct ternary compounds that crystallize in three different structures. During heating, the tetragonal $\alpha$-phase MgAgSb decomposes into the tetragonal $\beta$-phase ${\mathrm{Mg}}_5{\mathrm{Ag}}_{4.32}{\mathrm{Sb}}_{4.5}$, and subsequently into the cubic ${\gamma }_{\mathrm{H}}$-phase ${\mathrm{Mg}}_5{\mathrm{Ag}}_{2.11}{\mathrm{Sb}}_4$, with the progressive precipitation of ${\mathrm{Ag}}_3\mathrm{Sb}$ and Sb phases. Upon cooling, the ${\gamma }_{\mathrm{H}}$-phase ${\mathrm{Mg}}_5{\mathrm{Ag}}_{2.11}{\mathrm{Sb}}_4$ converts to the cubic ${\gamma }_{\mathrm{L}}$-phase ${\mathrm{Mg}}_5{\mathrm{Ag}}_{1.18}{\mathrm{Sb}}_{3.63}$, accompanied by the $\alpha$ phase with substantial ${\mathrm{Ag}}_3\mathrm{Sb}$ and Sb phases. The electron localization function and bonding analysis indicate that the Mg-Ag-Sb compound is characterized by ionic interactions, wherein the Ag-Sb bond is notably weaker compared to the Mg-Sb bond. Molecular dynamics simulations reveal that significant Mg and Ag migration within the Sb sublattices underlies the phase evolution. Surprisingly, the rapid atomic diffusion also enables the restoration of the $\alpha$-phase MgAgSb upon annealing at low temperatures. Successful recovery of the $\alpha$ phase can be realized even after 500 thermal cycles between 423–673 K, demonstrating unprecedented thermal recoverability.