Hydroplastification offers an eco-friendly route to plasticize hydrophilic polymers by enhancing chain mobility, yet it remains ineffective for rigid crystalline nanomaterials. Here, we report a novel plasticity mechanism governed by mesopores acting as structural free volume between rigid hydrophilic nanocrystals, cellulose nanocrystals (CNCs). During the mechanosorptive vapor training process to change the mesopores, the plasticity arises from a dynamic competition between external tensile stress and internal capillary forces. We establish a semiquantitative constitutive relationship revealing that the network’s wet-state yield strength is critically regulated by the aspect ratio and surface chemistry of the nanocrystals. Specifically, the weak acidity of carboxyl groups ($pKa\sim 3.2$) combined with high aspect ratios creates a defect-dominated soft skeleton; here, the internal Laplace pressure ($\sim 26\mathrm{MPa}$) exceeds the yield strength, driving the splitting of bottlelike mesopores into smaller domains until a geometric locking limit defined by the nanocrystal diameter is reached. In contrast, strong electrostatic repulsion from sulfonates ($pKa\sim 1.0$) or spatial jamming in short nanorods results in rigid networks that resist this capillary-driven densification. This novel hydroplastification mechanism by capitalizing on mesopores with competing capillary force, opposing network stress, and external stress enables shaping CNC films into stable geometries, rather than lattice disruption. Moreover, such findings extend hydroplasticity concepts to rigid matter, establishing mesopore engineering as a physical pathway to regulate the plasticity of crystalline nanomaterials.

One of the roadblocks toward the implementation of a fault-tolerant superconducting quantum processor is impacts of ionizing radiation with the qubit substrate. Such impacts temporarily elevate the density of quasiparticles (QPs) across the device, leading to correlated qubit error bursts. The most damaging errors, $T_1$ errors, stem from QP tunneling across the qubit Josephson junctions (JJs). Recently, we demonstrated [M. McEwen, K. Miao, et al., Phys. Rev. Lett. 133, 240601 (2024)] that this type of error can be strongly suppressed by engineering the profile of superconducting gap at the JJs in a way that prevents QP tunneling. In this work, we identify a new type of impact-induced correlated error that persists in the presence of gap engineering. We observe that impacts shift the frequencies of the affected qubits, and thus lead to correlated phase errors. The frequency shifts are systematically negative, reach values up to 3 MHz, and last for $\sim 1\mathrm{ms}$. We provide evidence that the shifts originate from QP-qubit interactions in the JJ region. Further, we demonstrate that the shift-induced phase errors can be detrimental to the performance of quantum error correction protocols.

The thermodynamic limit is foundational to statistical mechanics, underlying our understanding of many-body phases. It assumes that, as the system size grows infinitely at fixed density of particles, unambiguous macroscopic phases emerge that are independent of the system’s boundary shape. We present explicit quantum spin and dimer Hamiltonians whose ground states violate this principle. Our construction relies on the previous mathematical work on classical dimers on the Aztec diamond and the square-octagon fortress, where geometry-dependent phase behaviors are observed in the infinite-size limit. We reverse engineer quantum spin Hamiltonians on the square and the square-octagon lattices whose ground states at the Rokhsar-Kivelson points are described by classical dimer coverings. On diamond-shaped domains, we find macroscopic boundary regions exhibiting distinct quantum phases from those on square-shaped domains. We study the nature of these phases by calculating the dimer-dimer and vison correlators and adapt Kasteleyn matrix based analytical and numerical methods for computing the vison correlator, which are significantly more efficient than standard Monte Carlo techniques. Our results show that the square-octagon lattice supports a single gapped short-range entangled phase, with exponentially decaying dimer correlators and a constant vison correlator. When the same model is considered on a diamond-shaped domain, two additional macroscopic regions emerge, with one near the corners and exhibiting staggered dimer order, and another exhibiting critical correlations.

A swimming microorganism stirs the surrounding fluid, creating a flow field that governs not only its locomotion and nutrient uptake, but also its interactions with other microorganisms and the environment. Despite its fundamental importance, capturing this flow field and unraveling its biological implications remains a challenge. Here, we report direct, time-resolved measurements of the three-dimensional (3D) flow field generated by a single, free-swimming microalga, Chlamydomonas reinhardtii, a model organism for microbial locomotion and flagellar dynamics. Supported by hydrodynamic modeling and simulations, our measurements resolve how established two-dimensional (2D) flow features such as in-plane vortices and the stagnation point emerge from and shape the full algal flow in 3D. Moreover, we reveal unexpected low-Reynolds-number flow phenomena including micron-sized vortex rings and periodically recurring translating vortices and uncover topological changes in the underlying flow structure associated with the puller-to-pusher transition of an alga. Biologically, access to the 3D flow field enables rigorous quantification of the alga’s energy expenditure, as well as its swimming and feeding efficiency, improving the precision of these physiological metrics. Taken together, our study demonstrates rich vortex dynamics in inertialess flows and shows their influence on microbial motility. The work also introduces an experimental method for mapping the fluid environment sculpted by beating flagella.

A remarkable feat of active matter physics is that systems as diverse as collections of self-propelled particles, nematics mixed with molecular motors, and interacting robots can all be described by symmetry-based continuum theories. These descriptions rely on reducing complex effects of individual motors to a few key active parameters, which increase with activity. Here we observe a striking anomaly in the continuum description of nonreciprocal active solids, a ubiquitous class of active materials. Using a combination of metamaterial experiments and coarse-graining theory we find that as microscopic activity increases, macroscale active response can vanish: more is less. In this highly active regime, nonaffine and localized modes prevail and destroy the large-scale signature of microscopic activity. These modes exist in any dilute periodic structure and emerge in random lattices below a percolation transition. Our results unveil a counterintuitive facet of active matter, offering new principles for engineering materials far from equilibrium.

Collisions between ultracold calcium monohydroxide (CaOH) molecules are realized and studied. Inelastic collision rate constants are measured for CaOH prepared in ground and excited vibrational states, and the electric field dependence of these rates is measured for molecules in single quantum states of the parity-doubled bending mode. Theoretical calculations of collision rate coefficients are performed and found to agree with measured values. The lowest collisional loss rates are for states with repulsive long-range potentials that shield ultracold molecules from loss channels at short distance. These results unveil the collisional behavior of parity-doublet molecules in the ultracold regime and lay the foundation for future experiments to evaporatively cool polyatomic molecules to quantum degeneracy.

Driving a quantum system out of equilibrium while preserving its subtle quantum mechanical correlations on large scales presents a major challenge, both fundamentally and for technological applications. At its core, this challenge is pinpointed by the question of how quantum effects can persist at asymptotic scales, analogous to quantum critical points in equilibrium. In this work, we construct such a scenario using fermions as building blocks. These fermions undergo an absorbing-to-absorbing state transition between two topologically distinct and quantum-correlated dark states. Starting from a microscopic, interacting Lindbladian, we derive an effective Lindblad-Keldysh field theory in which critical fermions couple to a bosonic bath with hydrodynamic fluctuations associated with particle number conservation. A key feature of this field theory is an emergent symmetry that protects the purity of the fermions’ state even in the presence of the thermal bath. We quantitatively characterize the critical point using a leading-order expansion around the upper critical dimension, thereby establishing the nonequilibrium universality class of fermions. The symmetry protection mechanism, which exhibits parallels to the problem of directed percolation, suggests a pathway toward a broader class of robust, universal quantum phenomena in fermionic systems.

Many quantum networking applications require efficient photonic interfaces to quantum memories which can be produced at scale and with high yield. Synthetic diamond offers unique potential for the implementation of this technology as it hosts color centers which retain coherent optical interfaces and long spin coherence times in nanophotonic structures. Here, we report a technique enabling wafer-scale processing of thin-film diamond that combines ion implantation and membrane liftoff, high-quality overgrowth, targeted color center implantation, and serial, high-throughput thermocompression bonding with yields approaching unity. The deterministic deposition of thin diamond membranes onto semiconductor substrates facilitates consistent integration of photonic crystal cavities with silicon-vacancy (SiV) quantum memories. We demonstrate reliable, strong coupling of SiVs to photons with cooperativities approaching 50. Furthermore, we show that photonic crystal cavities can be reliably fabricated across several membranes bonded to the same handling chip. Our platform enables modular fabrication where the photonic layer can be integrated with functionalized substrates featuring electronic control lines such as coplanar waveguides for microwave delivery. Finally, we implement passive optical packaging with sub-decibel insertion loss. Together, these advances pave the way to the scalable assembly of optically addressable quantum memory arrays which are a key building block for modular photonic quantum interconnects.

Anyons are two-dimensional particles with fractional exchange statistics that emerge as elementary excitations of fractional quantum Hall phases. Experimentally, their exchange statistics can be measured in the edge-state Fabry-Pérot interferometer, wherein the presence of $N_{qp}$ localized anyons contributes a phase $N_{qp}{\theta }_a$ to the interference pattern where ${\theta }_a$ is twice the exchange phase. Here we report the observation of large, hysteretic phase jumps in a monolayer graphene Fabry-Pérot interferometer at $\nu =1/3$. When the filling factor is increased from $\nu <1/3$ toward the center of the plateau, we observe phase slips with magnitude $\mathrm{\Delta }\theta \approx 2\pi /3$, consistent with the addition of individual quasiparticles to the interferometer bulk. These phase slips occur as instantaneous jumps in the interference signal, with intervals between the jumps indicating quasiparticle equilibration times exceeding 20 min. We use this long timescale to investigate the effect of changes in interferometer area $A_I$ and $N_{qp}$ independently at fixed magnetic field, revealing a striking memory effect in the phase slip magnitude. In particular, as the $\nu =1/3$ plateau is approached from higher filling, we observed phase slips with $\mathrm{\Delta }\theta$ significantly larger than $2\pi /3$ over the same range of gate voltage where quantized jumps are seen for increasing $\nu$. We discuss this asymmetry in terms of bulk-edge coupling of quasiparticles localized near the edge or in the bulk, and argue that this effect can be qualitatively reconciled with theoretical expectations for strongly interacting quasiparticles in the presence of weak disorder and strongly nonequilibrium charge dynamics. Besides providing a replication of interferometric measurements sensitive to ${\theta }_a$, our results highlight the key role played by charge dynamics on signatures of the anyon phase, and demonstrate that fractional quasiparticles can be indefinitely localized in nonequilibrium configurations.

Mitigating ancilla-mediated error channels is a critical challenge in controlling high-quality superconducting cavities using circuit quantum electrodynamics (cQED). We address this by weakening the dispersive coupling while demonstrating fast, high-fidelity multimode control through transmon-mediated sideband interactions. We implement transmon-cavity SWAP gates with speeds up to 30 times larger than the bare dispersive coupling. Combined with transmon rotations, this enables universal state preparation in a single mode, though achieving unitary gates and extending control to multiple modes remains a challenge. In this work, we overcome this limitation by introducing two control strategies: (i) a shelving technique that stores populations in sideband-transparent states, and (ii) a method that exploits the dispersive shift to implement photon-number-selective transmon-cavity SWAP gates. We use these protocols to prepare Fock and binomial code states across any of the ten modes of a multimode cavity with millisecond coherence times—serving as a multimode quantum memory. We demonstrate unitaries that encode and decode an arbitrary qubit state from the transmon into corresponding vacuum and Fock state superpositions, as well as entangled NOON states of cavity mode pairs—a scheme extendable to arbitrary multimode Fock encodings in the cavity modes. Furthermore, we implement a new binomial encoding gate that converts arbitrary transmon superpositions into binomial code states in any cavity mode at a rate exceeding the dispersive shifts in our system, achieving an average post-selected state fidelity of 96.3% in a $4\mathrm{\mu }\mathrm{s}$ gate time. By using precalibrated transmon and sideband pulses, our work demonstrates multimode control with significantly reduced calibration overhead, enabling efficient unitary operations using sideband interactions in multimode cQED systems.

Superconductors are famously capable of supporting persistent electrical currents, that is, currents that flow without any measurable decay as long as the material is kept in the superconducting state. We introduce here a class of materials—superconducting altermagnets—that can both generate and carry persistent spin currents. This includes spin-polarized electrical supercurrent as well as pure spin supercurrent that facilitates spin transport in the absence of any charge transport. A key to this remarkable property is the realization that the leading superconducting instability of altermagnetic metals consists of two independent condensates formed of spin-up and spin-down electrons. In the nonrelativistic limit, the two condensates are decoupled and can, thus, naturally support persistent currents with any spin polarization, including pure spin supercurrents realized in the charge counterflow regime. We describe a novel “spin-current dynamo effect” that can be used to generate pure spin supercurrent in such systems by driving a charge current along certain crystallographic directions. Away from the nonrelativistic limit, when spin-orbit interactions and magnetic disorder are present, we find that the spin current generically develops spatial oscillations but, importantly, no dissipation or decay. This is in stark contrast to spin currents in normal diffusive metals which tend to decay on relatively short length scales. We illustrate the above properties by performing model calculations relevant to two distinct classes of altermagnets and various device geometries.

Triboelectric charging of insulating particles through contact is critical in diverse physical and engineering processes, from dust storms and volcanic eruptions to industrial powder handling. However, many experiments over the years have consistently revealed counterintuitive charging patterns, including variable impact charge under identical conditions, charge sign reversal with repeated impacts, and bipolar charging of differently sized particles. Existing computational models cannot predict these patterns; they either rely on oversimplified heuristics or require inaccessible detailed surface properties. We present a stochastic scaling model (SSM) for particle charging that unifies same-material (particle-particle) and different-material (particle-wall) charging in a single theoretical framework. The model is grounded in a physics-based stochastic closure by the mean, variance, skewness, and minimum impact charge measured in a highly controlled reference experiment. To test the SSM, we implemented it in an open-source Lagrangian-Eulerian computational fluid dynamics solver. When simulating 300 000 insulating particles transported by turbulent wall-bounded flows, the SSM takes less than 0.01% of the CPU time. By scaling the statistical parameters of the reference impact to each collision, the new model reproduces the complex charging patterns observed in experiments without requiring surface-level first-principles inputs. The SSM offers a physically grounded route to large-scale simulations of electrostatic effects across many fields of particle-laden flows.

We report measurements of phase-matched nonlinear x-ray and optical mixing from single-crystal silicon using subresonant 0.95 eV laser pulses and 9.5 keV hard x-ray pulses from the Linear Coherent Light Source free-electron laser. The mixing signal appears as energy and momentum sidebands to the elastic Bragg peak. It is proportional to the magnitude squared of the relevant temporal and spatial Fourier components of the optically induced microscopic charges and currents. We measure the first- and second-order sideband to the 220 Bragg peak and find that the efficiency is maximized when the applied field is along the reciprocal lattice vector. For an optical intensity of approximately ${10}^{12}\mathrm{W}/{\mathrm{cm}}^2$, we measure peak efficiencies of $3\times {10}^{-7}$ and $3\times {10}^{-10}$ for the first- and second-order sideband, respectively (relative to the elastic Bragg peak). The first-order sideband is consistent with induced microscopic currents along the applied electric field and an isotropic response. The second-order sideband depends nontrivially on the optical field orientation and is consistent with an anisotropic response originating from induced charges along the bonds with ${\mathrm{C}}_{3v}$ site symmetry. The results agree well with first-principles Floquet-Bloch calculations.

Biological and robotic systems often operate in confined environments where materials must be gathered without centralized control. Inspired by the effective collection strategies of aquatic worms (Lumbriculus variegatus and Tubifex tubifex), we investigate how active filaments autonomously aggregate dispersed non-Brownian passive particles. We study this process across four platforms: living worms, a robotic chain, Brownian dynamics simulations of active polymers, and a coarse-grained toy model. We show that aggregation emerges from repeated contact and body deformation—effectively, a sweeping or brooming motion—and demonstrate that clustering dynamics are governed by filament length and bending stiffness. Across systems, particle gathering follows a shared aggregation-fragmentation process, where the steady-state cluster size approximately follows a quadratic increase with the effective width of the path cleared by the filament $W$ relative to the domain size $D$. This scaling provides a minimal geometric baseline for comparison across active filamentous systems. We find that filament flexibility modulates $W$, enabling more flexible filaments to sweep larger areas and collect more particles. These results establish a unifying framework for understanding how shape and flexibility influence transport and organization in active filament systems and filamentous robots.

Stabilizing and shaping autonomous flows of active fluids is a fundamental challenge and a prerequisite for applications. We embed a light-responsive microtubule-based nematic in a proportional-integral control loop that adjusts the applied light intensity in response to real-time measurements of the spatially averaged flow speed. The self-regulating hardware-software-wetware system maintains a target flow speed against external or internal perturbations, including protein aging and aggregation, sample-to-sample variability, and temperature variation. Varying the controller’s gains reveals antagonistic roles between feedback and intrinsic processes, leading to nontrivial dynamics observed in fluctuation spectra. In particular, oscillations emerge from the interplay between the controller, motor binding kinetics, and active hydrodynamic relaxation. Accounting for the underlying binding timescale, our coarse-grained model and nematohydrodynamics simulations corroborate these observations. This work provides insight into the coupled dynamics of controlled active matter, laying the foundation for spatiotemporal patterning of active stress to generate and stabilize new dynamical configurations.

A major challenge of AI plus science lies in its inherent incompatibility: Today’s AI is primarily based on connectionism, while science depends on symbolism. To bridge the two worlds, we propose a framework to seamlessly synergize Kolmogorov-Arnold networks (KANs) and science. The framework highlights KANs’ usage for three aspects of scientific discovery: identifying relevant features, revealing modular structures, and discovering symbolic formulas. The synergy is bidirectional: science to KAN (incorporating scientific knowledge into KANs), and KAN to science (extracting scientific insights from KANs). We highlight major new functionalities in pykan: (1) MultKAN, KANs with multiplication nodes, (2) kanpiler, a KAN compiler that compiles symbolic formulas into KANs; (3) tree converter, convert KANs (or any neural networks) into tree graphs. Based on these tools, we demonstrate KANs’ capability to discover various types of physical laws, including conserved quantities, Lagrangians, symmetries, and constitutive laws.

Understanding how the vibrational and thermal properties of solids are influenced by atomistic structural disorder is of fundamental scientific interest and paramount to designing materials for next-generation energy technologies. While several studies indicate that structural disorder strongly influences the thermal conductivity, the fundamental physics governing the disorder-conductivity relation remains elusive. Here we show that order-of-magnitude, disorder-induced variations of conductivity in network solids can be predicted from a “bond-network” entropy, an atomistic structural descriptor that quantifies heterogeneity in the topology of the atomic-bond network. We employ the Wigner formulation of thermal transport to demonstrate the existence of a relation between the bond-network entropy and observables such as smoothness of the vibrational density of states and macroscopic conductivity. We also show that the smoothing of the vibrational density of states encodes information about the thermal resistance induced by disorder and can be directly related to phenomenological models for phonon-disorder scattering based on the semiclassical Peierls-Boltzmann equation. Our findings rationalize the conductivity variations of disordered carbon polymorphs ranging from nanoporous electrodes to defective graphite used as a moderator in nuclear reactors.

Quantum processors based on neutral atoms trapped in arrays of optical tweezers have appealing properties, including relatively easy qubit number scaling and the ability to engineer arbitrary gate connectivity with atom movement. However, these platforms are inherently prone to atom loss, and the ability to replace lost atoms during a quantum computation is an important but previously elusive capability. Here, we demonstrate the ability to measure and reinitialize, and if necessary replace, a subset of atoms while maintaining coherence in other atoms. This allows us to perform logical circuits that include single- and two-qubit gates as well as repeated midcircuit measurement while compensating for atom loss. We highlight this capability by performing up to 41 rounds of syndrome extraction in a repetition code, and combine midcircuit measurement and atom replacement with real-time conditional branching to demonstrate heralded state preparation of a logically encoded Bell state. Finally, we demonstrate the ability to replenish atoms in a tweezer array from an atomic beam while maintaining coherence of existing atoms—a key step toward execution of logical computations that last longer than the lifetime of an atom in the system.

We provide a polynomial-time classical algorithm for noisy quantum circuits. The algorithm computes the expectation value of any observable for any circuit, with a small average error over input states drawn from an ensemble (e.g., the computational basis). Our approach is based upon the intuition that noise exponentially damps nonlocal correlations relative to local correlations. This enables one to classically simulate a noisy quantum circuit by keeping track of only the dynamics of local quantum information. Our algorithm also enables sampling from the output distribution of a circuit in quasipolynomial time, so long as the distribution anticoncentrates. A number of implications are discussed, including a fundamental limit on the efficacy of noise mitigation strategies: For constant noise rates, any quantum circuit for which error mitigation succeeds in polynomial-time on most input states can also be classically simulated in polynomial-time on most input states. Our algorithms scale exponentially in the inverse noise rate, which is fundamental and makes them impractical for current quantum devices.

In equilibrium, the physical properties of matter are set by the interactions between the constituents. In contrast, the energy input of the individual components controls the behavior of synthetic or living active matter. Great progress has been made in understanding the emergent phenomena in active fluids, though their inability to resist shear forces hinders their practical use. This motivates the exploration of active solids as shape-shifting materials, yet, we lack controlled synthetic systems to devise active solids with unconventional properties. Here we build active elastic beams from dozens of active colloids and unveil complex emergent behaviors such as self-oscillations or persistent rotations. Developing tensile tests at the microscale, we show that the active beams are ultrasoft materials, with large (nonequilibrium) fluctuations. Combining experiments, theory, and stochastic inference, we show that the dynamics of the active beams can be mapped on different phase transitions which are tuned by boundary conditions. More quantitatively, we assess all relevant parameters by independent measurements or first-principles calculations, and find that our theoretical description agrees with the experimental observations. Our results demonstrate that the simple addition of activity to an elastic beam unveils novel physics and can inspire design strategies for active solids and functional microscopic machines.