Friction causes mechanical energy dissipation and material degradation in machinery and devices. While phononic friction is well understood via anharmonic lattice dynamics, the physics of electronic friction remains unclear due to challenges in separating electronic degrees of freedom from phononic ones in experiments and analyzing the nonequilibrium interactions between ionic movement and electronic dynamics in theory. To tackle this problem, we construct a sliding device featuring 2D crystalline interfaces that possess ultrasmooth and minimally interacting surfaces, achieving the state of structural superlubricity with no wear and minimal friction. Using electrical and mechanical controls, we tuned the nature of interfacial electronic coupling and charge densities in materials in an on-device manner, which allows us to disentangle the electron and phonon contributions to friction. Our experimental data and theoretical analysis supported by first-principles calculations demonstrate that electronic friction can well surpass phononic contributions and dominate energy dissipation at structural superlubricity contacts. These findings offer fresh insights into the mechanism of electronic friction and promising opportunities for friction control in device applications.

Noncollinear magnets are notoriously difficult to describe within first-principles approaches based on density-functional theory because of the presence of low-lying spin excitations. At the level of ground-state calculations, several methods exist to constrain the magnetic moments to a predetermined configuration and thereby accelerate convergence toward self-consistency. Their use in a perturbative context, however, remains very limited. Here, we present a general methodological framework to achieve parametric control over the local spin moments at the linear-response level. Our strategy builds on the concept of Legendre transform to switch between various flavors of magnetic functionals and to relate their second derivatives via simple linear-algebra operations. Thereby, we can address an arbitrary response function at the time-dependent density-functional level of theory with optimal accuracy and minimal computational effort. In the low-frequency limit, we identify the leading correction to the existing adiabatic formulation of the problem [S. Ren et al., Phys. Rev. X 14, 011041 (2024)], consisting in a renormalization of the phonon and magnon masses due to electron inertia. As a demonstration, we apply our methodology to the terahertz optical response of bulk ${\mathrm{CrI}}_3$ and ${\mathrm{Cr}}_2{\mathrm{O}}_3$, where we identify contributions from hybrid (electro)magnons with mixed spin-lattice character.

Bosonic excitations like phonons and magnons dominate the low-temperature transport of magnetic insulators. Similar to electronic Hall responses, the thermal Hall effect (THE) of charge-neutral bosons has been proposed as a powerful tool for probing topological properties of their wave functions. For example, the intrinsic contribution of the THE of a perfectly clean system is directly governed by the distribution of Berry curvature, and many experiments on topological magnon and phonon insulators have been interpreted in this way. However, disorder is inevitably present in any material, and its contribution to the THE has remained poorly understood. Here, we develop a rigorous kinetic theory of the extrinsic side-jump contribution to the THE of bosons. We show that the extrinsic THE is always relevant for bosonic systems and can be of the same order as the intrinsic one but sensitively depends on the type of local imperfection. We study different types of impurities and show that a THE can even arise as a pure impurity-induced effect in a system with a vanishing intrinsic contribution. As a side product, we also generalize existing results for the electronic anomalous Hall effect to general types of impurities beyond the standard assumption of local potential scattering. We discuss the importance of our results for the correct interpretation of THE measurements, for example, in the Kitaev magnet $\alpha \text{−}{\mathrm{RuCl}}_3$, and provide a ready-to-use formula for comparison to experimental data.

We present an experimental protocol for the fabrication and characterization of scalable microarchitected elastic waveguides. Using silicon microfabrication techniques, we develop free-standing 2D truss-based architected waveguides with a maximum diameter of 80 mm, unit cells size of $100\mathrm{\mu }\mathrm{m}$, and minimum beam width of $5\mathrm{\mu }\mathrm{m}$, thus achieving scale separation. To characterize elastic wave propagation, we introduce a custom-built scanning optical pump-probe experiment that enables contactless excitation of elastic wave modes and full spatiotemporal reconstruction of wave propagation across hundreds of unit cells with subunit cell resolution. Results on periodic architectures show excellent agreement with finite element simulations and equivalent experimental data at larger length scales. Motivated by scalable computational inverse design, we fabricate a specific example of a spatially graded waveguide and demonstrate its ability to guide elastic waves along an arbitrary predesigned path.

Polymers with active segments constitute prospective future materials and are used as a model for some biological systems such as chromatin. The directions of the active forces are typically introduced with temporal or spatial correlations to establish directional motion of the chain and corresponding active dynamics. Instead, here we consider an active-passive copolymer, where the two segments differ only by the magnitude of their fluctuations and feature no artificial correlations. Here we show that although the model itself does not possess directional dynamics, if the chains are concentrated, directional persistent motion spontaneously arises as a consequence of the broken translational symmetry owing to the topological constraints. Using scaling arguments and simulations, we explain the phenomenon and describe the ensuing dynamics. Our work has thus far-reaching consequences for the mechanical properties of all dense active polymeric systems with heterogeneous fluctuations and in particular for chromatin conformation and dynamics that are crucial for biological functionality.

In the past decades, many efforts have focused on analyzing typical-case hardness in optimization and inference problems. Some recent work has pointed out that polynomial algorithms exist, running with a time that grows more than linearly with the system size, which can do better than linear algorithms, finding solutions to random problems in a wider range of parameters. However, a theory for polynomial and superlinear algorithms is in general lacking. Here, we examine the performance of the simulated annealing algorithm, a standard, versatile, and robust choice for solving optimization and inference problems, in the prototypical random $K$-SAT problem. For the first time, we show that the algorithmic thresholds depend on the time scaling of the algorithm with the size of the system. Indeed, one can identify not just one but different thresholds for linear, quadratic, and cubic regimes (and so on). This observation opens new directions in studying the typical case hardness in optimization problems.

The local structure and chemical environment critically govern transport properties in quantum and energy materials, but their atomic-scale manifestations remain elusive. Here we uncovered the atomic characteristics, intrinsic origin, and thermoelectric implications of chemical inhomogeneity in ${\text{Bi}}_2{\text{Te}}_3$-based materials with varying Se alloying. Using aberration-corrected electron microscopy, chemical bonding analysis, and thermodynamic modeling, we revealed pronounced Se enrichment in the central layer of a quintuple structure and constructed precise atomic models for this intrinsic inhomogeneity. Atomic-scale charge-density and valence-state measurements provided direct experimental evidence for chemical bonding disparity in ${\text{Bi}}_2{\text{Te}}_3$, with more covalent bonding in central layers and greater ionic bonding in outer layers. We also showed that this bonding disparity, coupled with entropy-enthalpy competition, drove the observed site preference and local chemical segregation in ${\text{Bi}}_2{\text{Te}}_{3\text{−}x}{\text{Se}}_x$ ($x=0.5$, 1.0, 1.5, 2.0, and 2.5). The site-resolved bonding analyses further demonstrated that selective Se occupation tuned both ionic and covalent components of Bi–$X$ bonds in ${\text{Bi}}_2{\text{Te}}_{3\text{−}x}{\text{Se}}_x$. Therefore, this intrinsic chemical inhomogeneity modulated bonding polarity, the electronic band gap, and the Seebeck coefficient, serving as a fundamental determinant of thermoelectric performance. Our findings highlight atomic-scale chemical inhomogeneity as a key factor in tailoring functional properties of layered quantum materials beyond thermoelectrics.

Anderson localization is a quantum phenomenon in which disorder localizes electronic wave functions. In this work, we propose a new approach to study Anderson localization based on the density matrix formalism. Drawing an analogy to the standard transfer matrix method, we extract the localization length from the modular density matrix in quasi-one-dimensional systems. This approach successfully captures the metal-insulator transition in the three-dimensional Anderson model and in the two-dimensional Anderson model with spin-orbit coupling. It can be also readily extended to multi-orbital systems. We further generalize the formalism to interacting systems, showing that the one-dimensional spinless attractive model exhibits the expected metallic phase, consistent with previous studies. More importantly, we demonstrate the existence of a two-dimensional metallic tendency in the presence of Hubbard interactions and disorder. This method offers a new perspective on Anderson localization and its interplay with interactions.

Hyperuniform particle arrangements are characterized by a local number variance within a spherical window of radius $R$ that grows more slowly than the volume of the window, i.e., $R^d$, in $d$-dimensional Euclidean space. We generalize this concept to describe the large-scale behavior of particle systems in which particles carry weights: internal degrees of freedom such as scalars (charges and masses), vectors (electric dipole moments, velocities, and torques), pseudovectors (spins and angular momenta), directors (bond orientations), tensors (quadrupole moments), or extrinsic local attributes (Voronoi-cell characteristics). The underlying hyperuniform arrangement may be ordered (crystals and quasicrystals) or disordered, the latter of which has been extensively studied for its novel properties. Our generalization extends hyperuniformity from fluctuations in particle positions to fluctuations in the spatial distribution of weights and examines how weighted fluctuations compare to their unweighted counterparts. We derive generalized weighted pair correlation functions, autocovariance functions, and spectral functions and show their relation to formulas for the local variance in weighted many-particle systems. Then, we apply our formalism to determine the hyperuniformity or nonhyperuniformity of bond-orientational ordered phases, dipolar liquid water, Voronoi-cell volumes, and certain ionic liquids in various Euclidean space dimensions. We demonstrate that hyperuniformity in the particle system does not necessarily translate to hyperuniformity of the weighted system. In fact, cases exist where a hyperuniform particle system becomes antihyperuniform when weighted and others where nonhyperuniform or antihyperuniform particle systems yield hyperuniform weighted systems. This theoretical framework provides a road map for quantifying large-scale fluctuations in weighted many-particle systems, offering a powerful tool for identifying systems with novel physical properties.

Periodically driven Floquet quantum many-body systems have revealed new insights into the rich interplay of thermalization and growth of entanglement. The phenomenology of dynamical freezing, whereby a translationally invariant many-body system exhibits emergent conservation laws and a slow growth of entanglement entropy at certain fixed ratios of a drive amplitude and frequency, presents a novel paradigm for retaining memory of an initial state up to late times. Previous studies of dynamical freezing have largely been restricted to a high-frequency Floquet-Magnus expansion and numerical exact diagonalization. Both techniques are unable to capture the slow approach to thermalization, or lack thereof, in a systematic fashion. By employing Floquet flow renormalization, where the time-dependent part of the Hamiltonian is gradually decoupled from the effective Hamiltonian using a sequence of unitary transformations, we unveil the universal approach to dynamical freezing and beyond, at asymptotically late times. We analyze the fixed-point behavior associated with the flow renormalization at and near freezing using both exact-diagonalization and tensor-network-based methods and contrast the results with the conventional prethermal phenomenon. For a generic nonintegrable spin Hamiltonian with a periodic cosine wave drive, the flow approaches an unstable fixed point with an approximate emergent symmetry. We observe that at freezing the thermalization timescales are delayed compared to away from freezing, and the flow trajectory undergoes a series of instanton events. Our numerical results are supported by analytical solutions to the flow equations.

Starting from a microscopic description of weak system-bath interactions, we derive from first principles a quantum master equation that does not rely on the well-known rotating wave approximation. This includes generic many-body systems, with Hamiltonians with vanishingly small energy spacings that forbid that approximation. The equation satisfies a general form of detailed balance, called KMS (Kubo-Martin-Schwinger), which ensures exact convergence to the many-body Gibbs state. Unlike the more common notion of GNS (Gelfand-Naimark-Segal) detailed balance, this notion is compatible with the absence of the rotating wave approximation. We show that the resulting Lindbladian dynamics not only reproduces the thermal equilibrium point up to a small renormalization of the system Hamiltonian, but it also approximates the true system evolution with an error that grows at most linearly in time, giving an exponential improvement upon previous estimates. This master equation has quasilocal jump operators, can be efficiently simulated on a quantum computer, and reduces to the usual Davies dynamics in the limit of a coarse-graining time much larger than the inverse of the smallest frequency difference. With it, we provide a rigorous model of many-body thermalization relevant to both open quantum systems and quantum algorithms.

Topological and crystalline orders of electrons both benefit from enhanced Coulomb interactions in partially filled Landau levels. In bilayer graphene (BLG), the competition between fractional quantum Hall liquids and electronic crystals can be tuned electrostatically. Applying a displacement field leads to Landau-level crossings, where the interaction potential is strongly modified due to changes in the orbital wave functions. Here, we leverage this control to investigate phase transitions between topological and crystalline orders at constant filling factors in the lowest Landau level of BLG. Using transport measurements in high-quality hBN-encapsulated devices, we study transitions as a function of displacement field near crossings of $N=0$ and $N=1$ orbitals. The enhanced Landau-level mixing near the crossing stabilizes electronic crystals at all fractional fillings, including a resistive state at $\nu =\frac{1}{3}$ and a reentrant integer quantum Hall state at $\nu =\frac{7}{3}$. On the $N=0$ side, the activation energies of the crystal and fractional quantum Hall liquid vanish smoothly and symmetrically at the transition, while the $N=1$ transitions out of the crystal appear discontinuous. Additionally, we observe quantized plateaus forming near the crystal transition at half filling of the $N=0$ levels, suggesting a paired composite fermion state stabilized by Landau level mixing.

One of the hallmarks of topological systems is the robust quantization of particle transport. It is the origin of the integer-valued quantum Hall conductivity and a potential tool for quantum information technology. Recent experiments on topological pumps constructed by using arrays of photonic waveguides and described by the (lattice-translational invariant) Aubry-André-Harper model, have demonstrated both integer and fractional transport of lattice solitons. In these systems, a background medium mediates interactions between photons via a Kerr nonlinearity and leads to the formation of self-bound multiphoton states. Upon increasing the interaction strength, a sequence of transitions was observed from a phase with integer transport in a pump cycle through different phases of fractional transport to a phase with no transport. We here present a quantum description of topological pumps of self-bound many-particle states in terms of an effective Hamiltonian of their center-of-mass (c.m.) motion, which allows one to introduce an effective band structure $E_{\mu }(K)$ with $K$ being the c.m. momentum and to classify topological phases in terms of generalized symmetries. We provide an explicit analytic expression of the effective Hamiltonian for few particles in the strong interaction limit and present numerical results in the more general case. We identify a topological invariant, an effective single-particle Chern number, which fully governs the soliton transport. Increasing the interaction strength in the Aubry-André-Harper model leads to a successive merging of c.m. bands, which is the origin of the observed sequence of topological phase transitions and also the potential breakdown of topological quantization for some interaction strength.

Understanding universal aspects of many-body systems is one of the central themes in modern physics. Recently, the stabilizer Rényi entropy (SRE) has emerged as a computationally tractable measure of nonstabilizerness, a crucial resource for fault-tolerant universal quantum computation. While numerical results suggested that the SRE in critical states can exhibit universal behavior, its comprehensive theoretical understanding has remained elusive. In this work, we develop a field-theoretical framework for the SRE in a ($1+1$)-dimensional many-body system and elucidate its universal aspects by using boundary conformal field theory. We demonstrate that the SRE is equivalent to a participation entropy in the Bell basis of a doubled Hilbert space, which can be calculated from the partition function of a replicated field theory with the interlayer line defect created by the Bell-state measurements. This identification allows us to characterize the universal contributions to the SRE on the basis of the data of conformal boundary conditions imposed on the replicated theory. We find that the SRE of the entire system contains a universal size-independent term determined by the noninteger ground-state degeneracy known as the $g$ factor. In contrast, we show that the mutual SRE exhibits the logarithmic scaling with a universal coefficient given by the scaling dimension of a boundary-condition-changing operator, which elucidates the origin of universality previously observed in numerical results. As a concrete demonstration, we present a detailed analysis of the Ising criticality, where we analytically derive the universal quantities at arbitrary Rényi indices and numerically validate them with high accuracy by employing tensor-network methods. These results establish a field-theoretical approach to understanding the universal features of nonstabilizerness in quantum many-body systems.

Light-induced quantum phases offer the potential for simple and powerful tuning of material properties. For example, simply illuminating 2D materials in the integer quantum Hall regime with polarized light is predicted to drive quantum phase transitions. Such phenomena are largely beyond the current frontier of solid state experiments due to technical limitations on laser intensity and material purity. However, the Harper-Hofstadter mapping which relates a two-dimensional integer quantum Hall system to a 1D quasicrystal enables the same polarization-dependent light-induced phase transitions to be observed using a quantum gas in a driven quasiperiodic optical lattice. We report results of such an experiment. We observe an interlaced phase diagram of localization-delocalization phase transitions as a function of drive polarization and amplitude. Elliptically polarized driving can stabilize an extended critical phase featuring multifractal wave functions; we observe signatures of this phenomenon in anomalous polarization-dependent subdiffusive transport. In this regime, increasing the strength of the quasiperiodic potential can enhance rather than suppress transport. These experiments demonstrate a simple method for synthesizing exotic multifractal states and exploring light-induced quantum phases across different dimensionalities.

We propose a novel computational strategy to study the glass transition of molecular fluids. Our approach combines the construction of simple yet realistic models with the development of Monte Carlo algorithms to accelerate equilibration and sampling. Inspired by the well-studied orthoterphenyl glass former, we construct a molecular model with an analogous triangular geometry and construct a “flip” Monte Carlo algorithm. We demonstrate that the flip Monte Carlo algorithm achieves a sampling speedup of about ${10}^9$ at the experimental glass transition temperature $T_g$. This allows us to systematically analyze the equilibrium structure and molecular dynamics of the model over a temperature regime previously inaccessible. We carefully compare the observed physical behavior to earlier studies that used atomistic models. In particular, we find that the glass fragility and the departure from the Stokes-Einstein relation are much closer to experimental observations. We characterize the development and temperature evolution of spatial correlations in the relaxation dynamics using both orientational and translational degrees of freedom. Excess wings emerge at intermediate frequencies in dynamic rotational spectra, and we directly visualize the corresponding molecular motion near $T_g$. Our approach can be generalized to a broad range of molecular geometries and paves the way to a deeper understanding of how molecular details may affect more universal physical aspects characterizing molecular liquids approaching their glass transition.

Spin waves (magnons) in two-dimensional (2D) materials have received increasing interest due to their unique states and potential for tunability. However, many interesting features of these systems, including Dirac points and topological states, occur at high frequencies, where experimental probes are limited. Here, we study a crystal formed by patterning a hexagonal array of holes in a perpendicularly magnetized thin film. Through simulation, we find that the magnonic band structure imitates that of graphene, but additionally has some kagomelike character and includes a few flat bands. Surprisingly, its nature can be understood using a nine-band tight-binding Hamiltonian. This clear analogy to 2D materials enables band-gap engineering in 2D, topological magnons along 1D phase boundaries, and spectrally isolated modes at 0D point defects. Interestingly, the 1D phase boundaries allow access to the valley degree of freedom through a magnonic analog of the quantum valley Hall insulator. These approaches can be extended to other magnonic systems, but are potentially more general due to the simplicity of the model, which resembles existing results from electron, phonon, photon, and cold-atom systems. This finding brings the physics of spin waves in 2D materials to more experimentally accessible scales, augments it, and outlines a few principles for controlling magnonic states.

The high complexity of many-body quantum dynamics means that essentially all analytical or numerical approaches either exploit special structure or are approximate in nature. One such approach—the memory function formalism—involves a carefully chosen split into “fast” and “slow” modes. An approximate model for the fast modes can then be used to solve for Green’s functions $G(z)$ of the slow modes, and the success of this approach depends on the accuracy of the fast space approximation. Using a formulation in operator Krylov space known as the recursion method, we prove the emergence of a universal random matrix description of the fast mode dynamics. This is captured by the “level-$n$ Green’s function” $G_n(z)$, which we show approaches universal scaling forms in the “fast limit” $n\to \infty$. Notably, this emergent universality can occur in both chaotic and nonchaotic systems, provided their spectral functions are sufficiently smooth. This universality of $G_n(z)$ turns out to be precisely analogous to the universality of eigenvalue correlations in random matrix theory (RMT), even though there is no explicit randomness present in the Hamiltonian. Concretely, at finite $z$ we show that $G_n(z)$ approaches the Wigner semicircle law, while if $G(z)$ is the Green’s function of certain hydrodynamical variables, we show that at low frequencies $G_n(z)$ is instead governed by the Bessel universality class from RMT. As an application of this universality, we give a new numerical method, the spectral bootstrap, for approximating spectral functions, including hydrodynamic transport data, from a finite number of Lanczos coefficients. Our proof involves a map to a Riemann-Hilbert problem which we solve using a steepest-descent-type method, rigorously controlled in the $n\to \infty$ limit. Via the steepest-descent procedure, we are led to a related Coulomb gas optimization problem, and we discuss how a recent conjecture—the “Operator Growth Hypothesis”—implies that chaotic operator dynamics can generically be identified with the critical point of a confinement transition in this Coulomb gas. These results elevate the recursion method from a useful numerical technique to a theoretically principled framework with universal content.

Over the past decade, autonomous stabilization of bosonic qubits has emerged as a promising approach for hardware-efficient protection of quantum information. However, applying these techniques to more complex encodings than the Schrödinger cat code requires exquisite control of high-order wave mixing processes. The challenge is to enable specific multiphotonic dissipation channels while avoiding unintended nonlinear interactions. In this work, we leverage a genuine six-wave mixing process enabled by a near Kerr-free Josephson element to enforce dissipation of quartets of excitations in a high-impedance superconducting resonator. Owing to residual nonlinearities stemming from stray inductances in our circuit, this dissipation channel is only effective when the resonator holds a specific number of photons. Applying it to the fourth excited state of the resonator, we show an order-of-magnitude enhancement of the state decay rate while only marginally impacting the relaxation and coherence of lower-energy states. Given that stray inductances could be strongly reduced through simple modifications in circuit design and that our methods can be adapted to activate even higher-order dissipation channels, these results pave the way toward the dynamical stabilization of four-component Schrödinger cat qubits and even more complex bosonic qubits.