Recent developments in stochastic thermodynamics have elucidated various relations between the entropy production rate (thermodynamic dissipation) and the physical limits of information processing in nonequilibrium dynamical systems. These findings have been actively utilized and have opened new perspectives in the analysis of real biological systems. In neuroscience also, the importance of quantifying entropy production has attracted increasing attention as a means to understand the properties of information processing in the brain. However, the relationship between entropy production rate and oscillations, which is prevalent in many biological systems, is unclear. For example, neural oscillations, such as delta, theta, and alpha waves, play crucial roles in brain information processing. Here, we derive a novel decomposition of the entropy production rate of linear Langevin systems. We show that one of the components of the entropy production rate, called the housekeeping entropy production rate, can be decomposed into independent positive contributions from oscillatory modes. Our decomposition enables us to calculate the contribution of oscillatory modes to the housekeeping entropy production rate. In addition, when the noise matrix of the Langevin equation is diagonal, the contribution of each oscillatory mode is further decomposed into the contribution of each element of the system. To demonstrate the utility of our decomposition, we apply our decomposition to an electrocorticography dataset recorded during awake and anesthetized conditions in monkeys, wherein the properties of oscillations change drastically. We show the consistent trends across different monkeys; i.e., the contribution of oscillatory modes from the delta band are larger in the anesthetized condition than in the awake condition, while those from the higher-frequency bands, such as the theta band, are smaller. These results allow us to interpret the change in neural oscillation in terms of stochastic thermodynamics and the physical limit of information processing.

Active nematic liquid crystals have the remarkable ability to spontaneously deform and flow in the absence of any external driving force. While living materials with orientational order, such as the mitotic spindle, can self-assemble in quiescent active phases, reconstituted active systems often display chaotic, periodic, or circulating flows under confinement. Quiescent in vitro active nematics are, therefore, quite rare, despite the prediction from active hydrodynamic theory that confinement between two parallel plates can suppress flows. This spontaneous flow transition—named the active Fréedericksz transition by analogy with the conventional Fréedericksz transition in passive nematic liquid crystals under a magnetic field—has been a cornerstone of the field of active matter. Here, we report experimental evidence that confinement in spherical droplets can stabilize the otherwise chaotic dynamics of a 3D extensile active nematics, giving rise to a quiescent—yet still out-of-equilibrium—nematic liquid crystal. The active nematics spontaneously flow when confined in larger droplets. The composite nature of our model system composed of extensile bundles of microtubules and molecular motors dispersed in a passive colloidal liquid crystal allows us to demonstrate how the interplay of activity, nematic elasticity, and confinement impacts the spontaneous flow transition. The critical diameter increases when motor concentration decreases or nematic elasticity increases. Experiments and simulations also demonstrate that the critical confinement depends on the confining geometry, with the critical diameter in droplets being larger than the critical width in channels. Biochemical assays reveal that neither confinement nor nematic elasticity impacts the energy-consumption rate, confirming that the quiescent active phase is the stable out-of-equilibrium phase predicted theoretically. Further experiments in dense arrays of monodisperse droplets show that fluctuations in the droplet composition can smooth the flow transition close to the critical diameter. In conclusion, our work provides experimental validation of the active Fréedericksz transition in 3D active nematics, with potential applications in human health, ecology, and soft robotics.

In chemical reactions, the nuclear motion of the molecules plays a crucial role in determining the reaction rates and outcomes. Employing the cold target recoil ion momentum spectroscopy and femtosecond pump-probe techniques, we perform a molecular-level study into the influence of nuclear vibrations on light-induced bimolecular reactions within ${\mathrm{H}}_2\text{−}{\mathrm{D}}_2$ dimers. The study focuses on the formation dynamics of ${\mathrm{D}}_2{\mathrm{H}}^{+}$ and ${\mathrm{H}}_2{\mathrm{D}}^{+}$ cations, shedding light on the interplay between translational and vibrational motions of the nuclei steering the bimolecular reactions. Our observations reveal a notable yield ratio of $1:1.6$ between ${\mathrm{H}}_2{\mathrm{D}}^{+}$ and ${\mathrm{D}}_2{\mathrm{H}}^{+}$ channels, accompanied with a faster formation of ${\mathrm{D}}_2{\mathrm{H}}^{+}$ compared to ${\mathrm{H}}_2{\mathrm{D}}^{+}$. Molecular dynamics simulations unveil that the faster vibrational motion of ${\mathrm{H}}_2^{+}$ than that of ${\mathrm{D}}_2^{+}$ upon single ionization within the dimer accounts for these differences. Our findings provide new insight into the time-resolved kinetic isotope effect on the bimolecular reactions, highlighting the critical relationship between nuclear vibrational motions and reaction dynamics.

Entanglement is an extraordinary feature of quantum mechanics. Sources of entangled optical photons were essential to test the foundations of quantum physics through violations of Bell’s inequalities. More recently, entangled many-body states have been realized via strong nonlinear interactions in microwave circuits with superconducting qubits. Here, we demonstrate a chip-scale source of entangled optical and microwave photonic qubits. Our device platform integrates a piezo-optomechanical transducer with a superconducting resonator which is robust under optical illumination. We drive a photon-pair generation process and employ a dual-rail encoding intrinsic to our system to prepare entangled states of microwave and optical photons. We place a lower bound on the fidelity of the entangled state by measuring microwave and optical photons in two orthogonal bases. This entanglement source can directly interface telecom wavelength time-bin qubits and gigahertz frequency superconducting qubits, two well-established platforms for quantum communication and computation, respectively.

Given the fundamental importance of combinatorial optimization across many diverse domains, there has been widespread interest in the development of unconventional physical computing architectures that can deliver better solutions with lower resource costs. However, a theoretical understanding of their performance remains elusive. We develop such understanding for the case of the coherent Ising machine (CIM), a network of optical parametric oscillators that can be applied to any quadratic unconstrained binary optimization problem. We focus on how the CIM finds low-energy solutions of the Sherrington-Kirkpatrick spin glass. As the laser gain of this system is annealed, the CIM interpolates between gradient descent on coupled soft spins to descent on coupled binary spins. By combining the Kac-Rice formula, the replica method, and supersymmetry breaking, we develop a detailed understanding of the evolving geometry of the high-dimensional energy landscape of the CIM as the laser gain increases, finding several phase transitions in the landscape, from flat to rough to rigid. Additionally, we develop a novel cavity method that provides a geometric interpretation of supersymmetry breaking in terms of the reactivity of a rough landscape to specific external perturbations. Our energy landscape theory successfully matches numerical experiments, provides geometric insights into the principles of CIM operation, and yields optimal annealing schedules.

Stimulated Brillouin scattering (SBS) is often an unwanted loss mechanism in both active and passive fibers. Highly multimode excitation of fibers has been proposed as a novel route toward efficient SBS suppression. Here, we develop a detailed, quantitative theory which confirms this proposal and elucidates the physical mechanisms involved. Starting from the vector optical and scalar acoustic equations, we derive appropriate nonlinear coupled mode equations for the signal and Stokes modal amplitudes and an analytical formula for the SBS (Stokes) gain with applicable approximations, such as the neglect of shear effects. This allows us to calculate the exponential growth rate of the Stokes power as a function of the distribution of power in a highly multimode signal. The peak value of the gain spectrum across the excited modes determines the SBS threshold—the maximum SBS-limited power that can be sent through the fiber. The theory shows that the peak SBS gain is greatly reduced by highly multimode excitation due to gain broadening and relatively weaker intermodal SBS gain. The inclusion of exact vector optical modes in the calculation is crucial in order to capture the incomplete intermodal coupling due to mismatch of polarization patterns of higher-order modes. We demonstrate that equal excitation of the 160 modes of a commercially available, highly multimode circular step index fiber raises the SBS threshold by a factor of 6.5 and find comparable suppression of SBS in similar fibers with a $\mathsf{D}$-shaped cross section.

We describe coupled nonequilibrium electron-phonon systems semiclassically—Ehrenfest dynamics for the phonons and quantum mechanics for the electrons—using a classical Monte Carlo approach that determines the nonequilibrium response to a large pump field. The semiclassical approach is expected to be accurate, because the phonons are excited to average energies much higher than the phonon frequency, eliminating the need for a quantum description. The numerical efficiency of this method allows us to perform a self-consistent time evolution out to very long times (tens of picoseconds), enabling us to model pump-probe experiments of a charge-density-wave (CDW) material. Our system is a half-filled, one-dimensional (1D) Holstein chain that exhibits CDW ordering due to a Peierls transition. The chain is subjected to a time-dependent electromagnetic pump field that excites it out of equilibrium, and then a second probe pulse is applied after a time delay. By evolving the system to long times, we capture the complete process of lattice excitation and subsequent relaxation to a new equilibrium, due to an exchange of energy between the electrons and the lattice, leading to lattice relaxation at finite temperatures. We employ an indirect (impulsive) driving mechanism of the lattice by the pump pulse due to the direct driving of the electrons. We identify two driving regimes, where the pump can either cause small perturbations or completely invert the initial CDW order. Our work successfully describes the ringing of the amplitude mode in CDW systems that has long been seen in experiment but never successfully explained by microscopic theory. We also describe the fluence-dependent crossover that inverts the CDW order parameter and changes the phonon dynamics. Finally, we illustrate how this method can examine a number of different types of experiments including photoemission, x-ray diffraction, and two-dimensional (2D) spectroscopy.

Deuterons are atomic nuclei composed of a neutron and a proton held together by the strong interaction. Unbound ensembles composed of a deuteron and a third nucleon have been investigated in the past using scattering experiments, and they constitute a fundamental reference in nuclear physics to constrain nuclear interactions and the properties of nuclei. In this work, $K^{+}\text{−}d$ and $p\text{−}d$ femtoscopic correlations measured by the ALICE Collaboration in proton-proton ($pp$) collisions at $\sqrt{s}=13\mathrm{TeV}$ at the Large Hadron Collider (LHC) are presented. It is demonstrated that correlations in momentum space between deuterons and kaons or protons allow us to study three-hadron systems at distances comparable with the proton radius. The analysis of the $K^{+}\text{−}d$ correlation shows that the relative distances at which deuterons and protons or kaons are produced are around 2 fm. The analysis of the $p\text{−}d$ correlation shows that only a full three-body calculation that accounts for the internal structure of the deuteron can explain the data. In particular, the sensitivity of the observable to the short-range part of the interaction is demonstrated. These results indicate that correlations involving light nuclei in $pp$ collisions at the LHC will also provide access to any three-body system in the strange and charm sectors.

Brain activity is characterized by brainwide spatiotemporal patterns that emerge from synapse-mediated interactions between individual neurons. Calcium imaging provides access to in vivo recordings of whole-brain activity at single-neuron resolution and, therefore, allows the study of how large-scale brain dynamics emerge from local activity. In this study, we use a statistical mechanics approach—the pairwise maximum entropy model—to infer microscopic network features from collective patterns of activity in the larval zebrafish brain and relate these features to the emergence of observed whole-brain dynamics. Our findings indicate that the pairwise interactions between neural populations and their intrinsic activity states are sufficient to explain observed whole-brain dynamics. In fact, the pairwise relationships between neuronal populations estimated with the maximum entropy model strongly correspond to observed structural connectivity patterns. Model simulations also demonstrated how tuning pairwise neuronal interactions drives transitions between observed physiological regimes and pathologically hyperexcitable whole-brain regimes. Finally, we use virtual resection to identify the brain structures that are important for maintaining the brain in a physiological dynamic regime. Together, our results indicate that whole-brain activity emerges from a complex dynamical system that transitions between basins of attraction whose strength and topology depend on the connectivity between brain areas.

We investigate the steady state of an ellipsoidal active nematic shell using experiments and numerical simulations. We create the shells by coating microsized ellipsoidal droplets with a protein-based active cytoskeletal gel, thus obtaining ellipsoidal core-shell structures. This system provides the appropriate conditions of confinement and geometry to investigate the impact of nonuniform curvature on an orderly active nematic fluid that features the minimum number of defects required by topology. We identify new time-dependent states where topological defects periodically oscillate between translational and rotational regimes, resulting in the spontaneous emergence of chirality. Our simulations of active nematohydrodynamics demonstrate that, beyond topology and activity, the dynamics of the active material are profoundly influenced by the local curvature and viscous anisotropy of the underlying droplet, as well as by external hydrodynamic forces stemming from the self-sustained rotational motion of defects. These results illustrate how the incorporation of curvature gradients into active nematic shells orchestrates remarkable spatiotemporal patterns, offering new insights into biological processes and providing compelling prospects for designing bioinspired micromachines.

We study the statistics of matrix elements of local operators in the basis of energy eigenstates in a paradigmatic, integrable, many-particle quantum theory, the Lieb-Liniger model of bosons with repulsive delta-function interactions. Using methods of quantum integrability, we determine the scaling of matrix elements with system size. As a consequence of the extensive number of conservation laws, the structure of matrix elements is fundamentally different from, and much more intricate than, the predictions of the eigenstate thermalization hypothesis for generic models. We uncover an interesting connection between this structure for local operators in interacting integrable models and the one for local operators that are not local with respect to the elementary excitations in free theories. We find that typical off-diagonal matrix elements $⟨\mathit{\mu }|\mathcal{O}|\mathit{\lambda }⟩$ in the same macrostate scale as $\exp \mathbf{(}-c^{\mathcal{O}}L\ln (L)-L{M}_{\mathit{\mu },\mathit{\lambda }}^{\mathcal{O}}\mathbf{)}$, where the probability distribution function for $M_{\mathit{\mu },\mathit{\lambda }}^{\mathcal{O}}$ is well described by Fréchet distributions and $c^{\mathcal{O}}$ depends only on macrostate information. In contrast, typical off-diagonal matrix elements between two different macrostates scale as $\exp (-d^{\mathcal{O}}{L}^2)$, where $d^{\mathcal{O}}$ depends only on macrostate information. Diagonal matrix elements depend only on macrostate information up to finite-size corrections.

Understanding the ability of particles to maneuver through disordered environments is a central problem in innumerable settings, from active matter and biology to electronics. Macroscopic particles ultimately exhibit diffusive motion when their energy exceeds the characteristic potential barrier of the random landscape. In stark contrast, wave-particle duality causes electrons in disordered media to come to rest even when the potential is weak—a remarkable phenomenon known as Anderson localization. Here, we present a hydrodynamic active system with wave-particle features, a millimetric droplet self-guided by its own wave field over a submerged random topography, whose dynamics exhibits localized statistics analogous to those of electronic systems. Consideration of an ensemble of particle trajectories reveals a suppression of diffusion when the guiding wave field extends over the disordered topography. We rationalize mechanistically the emergent statistics by virtue of the wave-mediated resonant coupling between the droplet and topography, which produces an attractive wave potential about the localization region. This hydrodynamic analog, which demonstrates how a classical particle may localize like a wave, suggests new directions for future research in various areas, including active matter, wave localization, many-body localization, and topological matter.

In an ideal perfectly straight multimode fiber with a circular core, the symmetry ensures that rotating the input wave front leads to a corresponding rotation of the output wave front. This invariant property, known as the rotational memory effect (RME), remains independent of the typically unknown output profile. The RME thus offers significant potential for imaging and telecommunication applications. However, in real-life fibers, this effect is degraded by intrinsic imperfections and external perturbations, and is challenging to observe because of its acute sensitivity to misalignments and aberrations in the optical setup. Building on a previously established method for precisely estimating fiber transmission properties, we demonstrate an accurate extraction of RME properties. Additionally, we introduce a comprehensive theoretical framework for both qualitative and quantitative analysis, which specifically links the angular-dependent correlation of the RME to the core deformation’s geometrical properties and the fiber’s mode characteristics. Our theoretical predictions align well with experimental data and simulations for various amounts of fiber distorsion. Finally, we demonstrate the ability to engineer wave fronts with significantly enhanced correlation across all rotation angles. This work enables accurate characterization of distributed disorder from the fabrication process and facilitates calibration-free imaging in multimode fibers.

The interplay of dynamical correlations and electronic ordering is pivotal in shaping phase diagrams of correlated quantum materials. In magic-angle twisted bilayer graphene, transport, thermodynamic, and spectroscopic experiments pinpoint a competition between distinct low-energy states with and without electronic order, as well as between localized and delocalized charge carriers. In this study, we utilize dynamical mean-field theory on the topological heavy fermion model of twisted bilayer graphene to investigate the emergence of electronic correlations and long-range order in the absence of strain. We contrast moment formation, Kondo screening, and ordering on a temperature basis and explain the nature of emergent correlated states based on three central phenomena: (i) the formation of local spin and valley isospin moments around 100 K, (ii) the ordering of the local isospin moments around 10 K preempting Kondo screening, and (iii) a cascadic redistribution of charge between localized and delocalized electronic states upon doping. At integer fillings, we find that low-energy spectral weight is depleted in the symmetric phase, while we find insulating states with gaps enhanced by exchange coupling in the zero-strain ordered phases. Doping away from integer filling results in distinct metallic states: a “bad metal” above the ordering temperature, where scattering off the disordered local moments suppresses electronic coherence, and a “good metal” in the ordered states with coherence of quasiparticles facilitated by isospin order. This finding reveals coherence from order as the microscopic mechanism behind the Pomeranchuk effect observed experimentally by Rozen et al. [Nature (London) 592, 214 (2021)] and by Saito et al. [Nature (London) 592, 220 (2021)]. Upon doping, there is a periodic charge reshuffling between localized and delocalized electronic orbitals leading to cascades of doping-induced Lifshitz transitions, local spectral weight redistributions, and periodic variations of the electronic compressibility ranging from nearly incompressible to negative. Our findings highlight the essential role of charge transfer, hybridization, and ordering in shaping the electronic excitations and thermodynamic properties in twisted bilayer graphene and provide a unified understanding of the most puzzling aspects of scanning tunneling spectroscopy, transport, and compressibility experiments.

Open system quantum dynamics can generate a variety of long-range entangled mixed states, yet it has been unclear in what sense they constitute phases of matter. To establish that two mixed states are in the same phase, as defined by their two-way connectivity via local quantum channels, we use the renormalization group (RG) and decoders of quantum error correcting codes. We introduce a real-space RG scheme for mixed states based on local channels which ideally preserve correlations with the complementary system, and we prove this is equivalent to the reversibility of the channel’s action. As an application, we demonstrate an exact RG flow of finite temperature toric code in two dimensions to infinite temperature, thus proving it is in the trivial phase. In contrast, for toric code subject to local dephasing, we establish a mixed-state toric code phase using local channels obtained by truncating an RG-type decoder and the minimum weight perfect matching decoder. We also discover a precise relation between mixed-state phase and decodability, by proving that local noise acting on toric code cannot destroy logical information without bringing the state out of the toric code phase.

In this paper, we develop a novel method to solve problems involving quantum optical systems coupled to coherent quantum feedback loops featuring time delays. Our method is based on exact mappings of such non-Markovian problems to equivalent Markovian driven dissipative quantum many-body problems. In this work, we show that the resulting Markovian quantum many-body problems can be solved (numerically) exactly and efficiently using tensor network methods for a series of paradigmatic examples, consisting of driven quantum systems coupled to waveguides at several distant points. In particular, we show that our method allows solving problems in so far inaccessible regimes, including problems with arbitrary long time delays and arbitrary numbers of excitations in the delay lines. We obtain solutions for the full real-time dynamics as well as the steady state in all these regimes. Finally, motivated by our results, we develop a novel mean-field approach, which allows us to find the solution semianalytically, and we identify parameter regimes where this approximation is in excellent agreement with our tensor network results.

In solution, DNA, the “most important molecule of life,” is a highly charged macromolecule that bears a unit of negative charge on each phosphate of its sugar-phosphate backbone. Although partially compensated by counterions (cations of the solution) adsorbed at or condensed near it, DNA still produces a substantial electric field in its vicinity, which is screened by buffer electrolytes at longer distances from the DNA. This electric field is experienced by any charged or dipolar species approaching and interacting with the DNA. So far, such a field has been explored predominantly within the scope of a primitive model of the electrolytic solution, not considering more complicated structural effects of the water solvent. In this paper, we investigate the distribution of electric field around DNA using linear response nonlocal electrostatic theory, applied here for helix-specific charge distributions, and compare the predictions of such a theory with specially performed, fully atomistic, large-scale, molecular dynamics simulations. Both approaches are applied to unravel the role of the structure of water at close distances to and within the grooves of a DNA molecule in the formation of the electric field. As predicted by the theory and reported by the simulations, the main finding of this study is that oscillations in the electrostatic potential distribution are present around DNA, caused by the overscreening effect of structured water. Surprisingly, electrolyte ions at physiological concentrations do not strongly disrupt these oscillations and are rather distributed according to these oscillating patterns, indicating that water structural effects dominate the short-range electrostatics. We also show that (i) structured water adsorbed in the grooves of DNA leads to a positive electrostatic potential core relative to the bulk, (ii) the Debye length some 10 Å away from the DNA surface is reduced, effectively renormalized by the helical pitch of the DNA molecule, and (iii) Lorentzian contributions to the nonlocal dielectric function of water, effectively reducing the dielectric constant close to the DNA surface, enhance the overall electric field. The impressive agreement between the atomistic simulations and the developed theory substantiates the use of nonlocal electrostatics when considering solvent effects in molecular processes in biology.

Quantifying how difficult it is to control an emerging infectious disease is crucial to public health decision-making, providing valuable evidence on if targeted interventions, e.g., quarantine and isolation, can contain spread or when population wide controls, e.g., lockdowns, are warranted. The disease reproduction number $R$ or growth rate $r$ are universally assumed to measure controllability because $R=1$ and $r=0$ define when infections stop growing and hence the state of critical stability. Outbreaks with larger $R$ or $r$ are therefore interpreted as less controllable and requiring more stringent interventions. We prove this common interpretation is impractical and incomplete. We identify a positive feedback loop among infections intrinsically underlying disease transmission and evaluate controllability from how interventions disrupt this loop. The epidemic gain and delay margins, which, respectively, define how much we can scale infections (this scaling is known as gain) or delay interventions on this loop before stability is lost, provide rigorous measures of controllability. Outbreaks with smaller margins necessitate more control effort. Using these margins, we quantify how presymptomatic spread, surveillance limitations, variant dynamics, and superspreading shape controllability and demonstrate that $R$ and $r$ measure controllability only when interventions do not alter timings between the infections and are implemented without delay. Our margins are easily computed, interpreted, and reflect complex relationships among interventions, their implementation, and epidemiological dynamics.

By applying a microwave drive to a specially designed Josephson circuit, we have realized a double-well model system: a Kerr oscillator submitted to a squeezing force. We have observed, for the first time, the spectroscopic fingerprint of a quantum double-well Hamiltonian when its barrier height is increased: a pairwise level kissing (coalescence) corresponding to the exponential reduction of tunnel splitting in the excited states as they sink under the barrier. The discrete levels in the wells also manifest themselves in the activation time across the barrier which, instead of increasing smoothly as a function of the barrier height, presents steps each time a pair of excited states is captured by the wells. This experiment illustrates the quantum regime of Arrhenius’s law, whose observation is made possible here by the unprecedented combination of low dissipation, time-resolved state control, 98.5% quantum nondemolition single shot measurement fidelity, and complete microwave control over all Hamiltonian parameters in the quantum regime. Direct applications to quantum computation and simulation are discussed.

Symmetries of three-dimensional periodic scalar fields are described by 230 space groups (SGs). Symmetries of three-dimensional periodic (pseudo)vector fields, however, are described by the spin-space groups (SSGs), which were initially used to describe the symmetries of magnetic orders. In SSGs, the real-space and spin degrees of freedom are unlocked in the sense that an operation could have different spatial and spin rotations. SSGs give a complete symmetry description of magnetic structures and have natural applications in the band theory of itinerary electrons in magnetically ordered systems with weak spin-orbit coupling. Altermagnetism, a concept raised recently that belongs to the symmetry-compensated collinear magnetic orders but has nonrelativistic spin plitting, is well described by SSGs. Because of the vast number and complicated group structures, SSGs have not yet been systematically enumerated. In this work, we exhaust SSGs based on the invariant subgroups of SGs, with spin operations constructed from three-dimensional (3D) real representations of the quotient groups for the invariant subgroups. For collinear and coplanar magnetic orders, the spin operations can be reduced into lower-dimensional real representations. As the number of SSGs is infinite, we consider only SSGs that describe magnetic unit cells up to 12 times crystal unit cells. We obtain 157 289 noncoplanar, 24 788 coplanar-noncollinear, and 1421 collinear SSGs. The enumerated SSGs are stored in an online database with a user-friendly interface. We develop an algorithm to identify SSGs for realistic materials and find SSGs for 1626 magnetic materials. We also discuss several potential applications of SSGs, including the representation theory, topological states protected by SSGs, structures of spin textures, and refinement of magnetic neutron diffraction patterns using SSGs. Our results serve as a solid starting point for further studies of symmetry and topology in magnetically ordered materials.