Among various proteinaceous nanocontainers and nanoparticles, the most promising ones for various applications in nano- and medical science appear to be those whose structures differ fundamentally from icosahedral viral capsids described by the paradigmatic Caspar-Klug model. By analyzing such anomalous assemblies represented in the Protein Data Bank, we identify a series of shells with square-triangular local order and find that most of them originate from short-period approximants of a dodecagonal tiling consisting of square and triangular tiles. Examining the nonequilibrium assembly of such packings, we propose a new method for obtaining periodic square-triangle approximants and then construct the simplest models of tetragonal, octahedral, and icosahedral shells based on cubic and icosahedral nets cut from the approximant structures. Since gluing the nets can change the distances between adjacent vertices of the resulting shell, we introduce an effective energy, the minimization of which equalizes these distances. While the obtained spherical polyhedra reproduce the structures of experimentally observed protein shells and nanoparticles, the principles of protein organization that we lay out, and the ensuing structural models, can help to discover and investigate similar systems in the future.

First-order phase transitions produce abrupt changes to the character of both ground and excited electronic states. Here we conduct electronic compressibility measurements to map the spin phase diagram and Landau level (LL) energies of monolayer ${\mathrm{WSe}}_2$ in a magnetic field. We resolve a sequence of first-order phase transitions between completely spin-polarized LLs and states with LLs of both spins. Unexpectedly, the LL gaps are roughly constant over a wide range of magnetic fields below the transitions, which we show reflects spin-polarized ground states with opposite spin excitations. These transitions also extend into compressible regimes, with a sawtooth boundary between full and partial spin polarization. We link these observations to the important influence of LL filling on the exchange energy beyond a smooth density-dependent contribution. Our results show that ${\mathrm{WSe}}_2$ realizes a unique hierarchy of energy scales where such effects induce reentrant magnetic phase transitions tuned by density and magnetic field.

Dynamical decoupling and Hamiltonian engineering are well-established techniques that have been used to control qubit systems. However, designing the corresponding methods for qudit systems has been challenging due to the lack of a Bloch sphere representation, more complex interactions, and additional control constraints. By identifying several general structures associated with such problems, we develop a formalism for the robust dynamical decoupling and Hamiltonian engineering of strongly interacting qudit systems. Our formalism significantly simplifies qudit pulse-sequence design while naturally incorporating robustness conditions necessary for experimental practicality. We experimentally demonstrate these techniques in a strongly interacting, disordered ensemble of spin-1 nitrogen-vacancy centers, achieving more than an order-of-magnitude improvement in coherence time over existing pulse sequences. We further describe how our techniques enable the engineering of exotic many-body phenomena such as quantum many-body scars, and open up new opportunities for quantum metrology with enhanced sensitivities. These results enable wide-reaching new applications for dynamical decoupling and Hamiltonian engineering in many-body physics and quantum metrology.

Fault-tolerant quantum computation with bosonic qubits often necessitates the use of noisy discrete-variable ancillae. In this work, we establish a comprehensive and practical fault-tolerance framework for such a hybrid system and synthesize it with fault-tolerant protocols by combining bosonic quantum error correction (QEC) and advanced quantum control techniques. We introduce essential building blocks of error-corrected gadgets by leveraging ancilla-assisted bosonic operations using a generalized variant of path-independent quantum control. Using these building blocks, we construct a universal set of error-corrected gadgets that tolerate a single-photon loss and an arbitrary ancilla fault for four-legged cat qubits. Notably, our construction requires only dispersive coupling between bosonic modes and ancillae, as well as beam-splitter coupling between bosonic modes, both of which have been experimentally demonstrated with strong strengths and high accuracy. Moreover, each error-corrected bosonic qubit is comprised of only a single bosonic mode and a three-level ancilla, featuring the hardware efficiency of bosonic QEC in the full fault-tolerant setting. We numerically demonstrate the feasibility of our schemes using current experimental parameters in the circuit-QED platform. Finally, we present a hardware-efficient architecture for fault-tolerant quantum computing by concatenating the four-legged cat qubits with an outer qubit code utilizing only beam-splitter couplings. Our estimates suggest that the overall noise threshold can be reached using existing hardware. These developed fault-tolerant schemes extend beyond their applicability to four-legged cat qubits and can be adapted for other rotation-symmetrical codes, offering a promising avenue toward scalable and robust quantum computation with bosonic qubits.

Ever since the discovery of the charge density wave (CDW) transition in the kagome metal ${\mathrm{CsV}}_3{\mathrm{Sb}}_5$, the nature of its symmetry breaking has been under intense debate. While evidence suggests that the rotational symmetry is already broken at the CDW transition temperature ($T_{\mathrm{CDW}}$), an additional electronic nematic instability well below $T_{\mathrm{CDW}}$ has been reported based on the diverging elastoresistivity coefficient in the anisotropic channel ($m_{E_{2g}}$). Verifying the existence of a nematic transition below $T_{\mathrm{CDW}}$ is not only critical for establishing the correct description of the CDW order parameter, but also important for understanding low-temperature superconductivity. Here, we report elastoresistivity measurements of ${\mathrm{CsV}}_3{\mathrm{Sb}}_5$ using three different techniques probing both isotropic and anisotropic symmetry channels. Contrary to previous reports, we find the anisotropic elastoresistivity coefficient $m_{E_{2g}}$ is temperature independent, except for a step jump at $T_{\mathrm{CDW}}$. The absence of nematic fluctuations is further substantiated by measurements of the elastocaloric effect, which show no enhancement associated with nematic susceptibility. On the other hand, the symmetric elastoresistivity coefficient $m_{A_{1g}}$ increases below $T_{\mathrm{CDW}}$, reaching a peak value of 90 at $T^{*}=20\mathrm{K}$. Our results strongly indicate that the phase transition at $T^{*}$ is not nematic in nature and the previously reported diverging elastoresistivity is due to the contamination from the $A_{1g}$ channel.

A characteristic feature of “quantum chaotic” systems is that their eigenspectra and eigenstates display universal statistical properties described by random matrix theory (RMT). However, eigenstates of local systems also encode structure beyond RMT. To capture this feature, we introduce a framework that allows us to compare the ensemble properties of eigenstates in local systems with those of pure random states. In particular, our framework defines a notion of distance between quantum state ensembles that utilizes the Kullback-Leibler divergence to compare the microcanonical distribution of entanglement entropy (EE) of eigenstates with a reference RMT distribution generated by pure random states (with appropriate constraints). This notion gives rise to a quantitative metric for quantum chaos that not only accounts for averages of the distributions but also higher moments. The differences in moments are compared on a highly resolved scale set by the standard deviation of the RMT distribution, which is exponentially small in system size. As a result, the metric can distinguish between chaotic and integrable behaviors and, in addition, quantify and compare the degree of chaos (in terms of proximity to RMT behavior) between two systems that are assumed to be chaotic. We implement our framework in local, minimally structured, Floquet random circuits, as well as a canonical family of many-body Hamiltonians, the mixed-field Ising model (MFIM). Importantly, for Hamiltonian systems, we find that the reference random distribution must be appropriately constrained to incorporate the effect of energy conservation in order to describe the ensemble properties of midspectrum eigenstates. The metric captures deviations from RMT across all models and parameters, including those that have been previously identified as strongly chaotic, and for which other diagnostics of chaos such as level spacing statistics look strongly thermal. In Floquet circuits, the dominant source of deviations is the second moment of the distribution, and this persists for all system sizes. For the MFIM, we find significant variation of the KL divergence in parameter space. Notably, we find a small region where deviations from RMT are minimized, suggesting that “maximally chaotic” Hamiltonians may exist in fine-tuned pockets of parameter space.

Understanding and manipulating emergent phases, which are themes at the forefront of quantum-materials research, rely on identifying their underlying symmetries. This general principle has been particularly prominent in materials with coupled electronic and magnetic degrees of freedom, in which magnetic order influences the electronic band structure and can lead to exotic topological effects. However, identifying symmetry of a magnetically ordered phase can pose a challenge, particularly in the presence of small domains. Here we introduce a multimodal approach for determining magnetic structures, which combines symmetry-sensitive optical probes, scattering, and group-theoretical analysis. We apply it to ${\mathrm{EuIn}}_2{\mathrm{As}}_2$, a material that has received attention as a candidate axion insulator. While first-principles calculations predict this state on the assumption of a simple collinear antiferromagnetic structure, subsequent neutron-scattering measurements reveal a much more intricate magnetic ground state characterized by two coexisting magnetic wave vectors reached by successive thermal phase transitions. The proposed high- and low-temperature phases are a spin helix and a state with interpenetrating helical and Néel antiferromagnetic order termed a “broken helix,” respectively. Employing a multimodal approach, we identify the magnetic structure associated with these two phases of ${\mathrm{EuIn}}_2{\mathrm{As}}_2$. We find that the higher-temperature phase is characterized by a variation of the magnetic moment amplitude from layer to layer, with the moment vanishing entirely in every third Eu layer. The lower-temperature structure is similar to the broken helix, with one important difference: Because of local strain, the relative orientation of the magnetic structure and the lattice is not fixed. Consequently, the symmetry required to protect the axion phase is not generically protected in ${\mathrm{EuIn}}_2{\mathrm{As}}_2$, but we show that it can be restored if the magnetic structure is tuned with uniaxial strain. Finally, we present a spin Hamiltonian that identifies the spin interactions that account for the complex magnetic order in ${\mathrm{EuIn}}_2{\mathrm{As}}_2$. Our work highlights the importance of a multimodal approach in determining the symmetry of complex order parameters.

Convincing evidence of domain growth in the heating of ultrastable glasses suggests that the equilibration dynamics of supercooled liquids could be driven by a nucleation and growth mechanism. We investigate this possibility by simulating the equilibration dynamics of a model glass during both heating and cooling between poorly and well-annealed states. Though we do observe the growth of domains during heating, we find that domains are absent during cooling. This absence is inconsistent with classical nucleation theory. By comparing the equilibration dynamics of our glass with that of two models with kinetic constraints, we demonstrate that dynamical facilitation generically leads to heating driven by domain growth and cooling without domains. Our results provide strong evidence that dynamical facilitation, not nucleation and interfacial-tension-driven domain growth, is the driving mechanism for the equilibration dynamics of glass formers.

Cells contain multiple condensates which spontaneously form due to the heterotypic interactions between their components. Although the proteins and disordered region sequences that are responsible for condensate formation have been extensively studied, the rule of interactions between the components that allow demixing, i.e., the coexistence of multiple condensates, is yet to be elucidated. Here, we construct an effective theory of the interaction between heteropolymers by fitting it to the molecular dynamics simulation results obtained for more than 200 sequences sampled from the disordered regions of human proteins. We find that the sum of amino acid pair interactions across two heteropolymers predicts the Boyle temperature qualitatively well, which can be quantitatively improved by the dimer pair approximation, where we incorporate the effect of neighboring amino acids in the sequences. The improved theory, combined with the finding of a metric that captures the effective interaction strength between distinct sequences, allowed the selection of up to three disordered region sequences that demix with each other in multicomponent simulations, as well as the generation of artificial sequences that demix with a given sequence. The theory points to a generic sequence design strategy to demix or hypermix thanks to the low-dimensional nature of the space of the interactions that we identify. As a consequence of the geometric arguments in the space of interactions, we find that the number of distinct sequences that can demix with each other is strongly constrained, irrespective of the choice of the coarse-grained model. Altogether, we construct a theoretical basis for methods to estimate the effective interaction between heteropolymers, which can be utilized in predicting phase separation properties as well as rules of assignment in the localization and functions of disordered proteins.

A relevant problem in the theory of open quantum systems is the lack of complete positivity of dynamical maps obtained after weak-coupling approximations, a famous example being the Redfield master equation. A number of approaches exist to recover well-defined evolutions under additional Markovian assumptions, but much less is known beyond this regime. Here, we propose a numerical method to cure the complete-positivity violation issue while preserving the non-Markovian features of an arbitrary original dynamical map. The idea is to replace its unphysical Choi operator with its closest physical one, mimicking recent work on quantum process tomography. We also show that the regularized dynamics is more accurate in terms of reproducing the exact dynamics, which allows us to heuristically push the utilization of these master equations in moderate coupling regimes, where the loss of positivity can have a relevant impact.

Numerous natural and human-made systems exhibit critical transitions whereby slow changes in environmental conditions spark abrupt shifts to a qualitatively distinct state. These shifts very often entail severe consequences; therefore, it is imperative to devise robust and informative approaches for anticipating the onset of critical transitions. Real-world complex systems can comprise hundreds or thousands of interacting entities, and implementing prevention or management strategies for critical transitions requires knowledge of the exact condition in which they will manifest. However, most research so far has focused on low-dimensional systems and small networks containing fewer than ten nodes or has not provided an estimate of the location where the transition will occur. We address these weaknesses by developing a deep-learning framework which can predict the specific location where critical transitions happen in networked systems with size up to hundreds of nodes. These predictions do not rely on the network topology, the edge weights, or the knowledge of system dynamics. We validate the effectiveness of our machine-learning-based framework by considering a diverse selection of systems representing both smooth (second-order) and explosive (first-order) transitions: the synchronization transition in coupled Kuramoto oscillators; the sharp decline in the resource biomass present in an ecosystem; and the abrupt collapse of a Wilson-Cowan neuronal system. We show that our method provides accurate predictions for the onset of critical transitions well in advance of their occurrences, is robust to noise and transient data, and relies only on observations of a small fraction of nodes. Finally, we demonstrate the applicability of our approach to real-world systems by considering empirical vegetated ecosystems in Africa.

Flocking, as paradigmatically exemplified by birds, is the coherent collective motion of active agents. As originally conceived, flocking emerges through alignment interactions between the agents. Here, we report that flocking can also emerge through interactions that turn agents away from each other. Combining simulations, kinetic theory, and experiments, we demonstrate this mechanism of flocking in self-propelled Janus colloids with stronger repulsion on the front than on the rear. The polar state is stable because particles achieve a compromise between turning away from left and right neighbors. Unlike for alignment interactions, the emergence of polar order from turn-away interactions requires particle repulsion. At high concentration, repulsion produces flocking Wigner crystals. Whereas repulsion often leads to motility-induced phase separation of active particles, here it combines with turn-away torques to produce flocking. Therefore, our findings bridge the classes of aligning and nonaligning active matter. Our results could help to reconcile the observations that cells can flock despite turning away from each other via contact inhibition of locomotion. Overall, our work shows that flocking is a very robust phenomenon that arises even when the orientational interactions would seem to prevent it.

The fundamental optical excitations of correlated transition-metal compounds are typically identified with multielectronic transitions localized at the transition-metal site, such as $dd$ transitions. In this vein, intense interest has surrounded the appearance of sharp, below-band-gap optical transitions, i.e., excitons, within the magnetic phase of correlated ${\mathrm{Ni}}^{2+}$ van der Waals magnets. The interplay of magnetic and charge-transfer insulating ground states in ${\mathrm{Ni}}^{2+}$ systems raises intriguing questions on the roles of long-range magnetic order and of metal-ligand charge transfer in the exciton nature, which inspired microscopic descriptions beyond typical $dd$ excitations. Here we study the impact of charge transfer and magnetic order on the excitation spectrum of the nickel dihalides ($\mathrm{Ni}{X}_2$, $X=\mathrm{Cl}$, Br, and I) using Ni-$L_3$ edge resonant inelastic x-ray scattering (RIXS). In all compounds, we detect sharp excitations, analogous to the recently reported excitons, and assign them to spin-singlet multiplets of octahedrally coordinated ${\mathrm{Ni}}^{2+}$ stabilized by intra-atomic Hund’s exchange. Additionally, we demonstrate that these excitons are dispersive using momentum-resolved RIXS. Our data evidence a ligand-mediated multiplet dispersion, which is tuned by the charge-transfer gap and independent of the presence of long-range magnetic order. This reveals the mechanisms governing nonlocal interactions of on-site $dd$ excitations with the surrounding crystal or magnetic structure, in analogy to ground-state superexchange. These measurements thus establish the roles of magnetic order, self-doped ligand holes, and intersite-coupling mechanisms for the properties of $dd$ excitations in charge-transfer insulators.

A ubiquitous problem in quantum physics is to understand the ground-state properties of many-body systems. Confronted with the fact that exact diagonalization quickly becomes impossible when increasing the system size, variational approaches are typically employed as a scalable alternative: Energy is minimized over a subset of all possible states and then different physical quantities are computed over the solution state. Despite remarkable success, rigorously speaking, all that variational methods offer are upper bounds on the ground-state energy. On the other hand, so-called relaxations of the ground-state problem based on semidefinite programming represent a complementary approach, providing lower bounds to the ground-state energy. However, in their current implementation, neither variational nor relaxation methods offer provable bound on other observables in the ground state beyond the energy. In this work, we show that the combination of the two classes of approaches can be used to derive certifiable bounds on the value of any observable in the ground state, such as correlation functions of arbitrary order, structure factors, or order parameters. We illustrate the power of this approach in paradigmatic examples of 1D and 2D spin-$1/2$ Heisenberg models. To improve the scalability of the method, we exploit the symmetries and sparsity of the considered systems to reach sizes of hundreds of particles at much higher precision than previous works. Our analysis therefore shows how to obtain certifiable bounds on many-body ground-state properties beyond energy in a scalable way.

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

Interfacing cold atoms with integrated nanophotonic devices could offer new paradigms for engineering atom-light interactions and provide a potentially scalable route for quantum sensing, metrology, and quantum information processing. However, it remains a challenging task to efficiently trap a large ensemble of cold atoms on an integrated nanophotonic circuit. Here, we demonstrate direct loading of an ensemble of up to 70 atoms into an optical microtrap on a nanophotonic microring circuit. Efficient trap loading is achieved by employing degenerate Raman-sideband cooling in the microtrap, where a built-in spin-motion coupling arises directly from the vector light shift of the evanescent-field potential on a microring. Atoms are cooled into the trap via optical pumping with a single free space beam. We have achieved a trap lifetime approaching 700 ms under continuous cooling. We show that the trapped atoms display large cooperative coupling and superradiant decay into a whispering-gallery mode of the microring resonator, holding promise for explorations of new collective effects. Our technique can be extended to trapping a large ensemble of cold atoms on nanophotonic circuits for various quantum applications.

An important outstanding challenge that must be overcome in order to fully utilize the XY surface code for correcting biased Pauli noise is the phenomenon of fragile temporal boundaries that arises during the standard logical state-preparation and measurement protocols. To address this challenge we propose a new logical state-preparation protocol based on locally entangling qubits into small Greenberger-Horne-Zeilinger-like states prior to making the stabilizer measurements that place them in the XY-code state. We prove that in this new procedure $O(\sqrt{n})$ high-rate errors along a single lattice boundary can cause a logical failure, leading to an almost quadratic reduction in the number of fault configurations compared to the standard state-preparation approach. Moreover, the code becomes equivalent to a repetition code for high-rate errors, guaranteeing a 50% code-capacity threshold during state preparation for infinitely biased noise. With a simple matching decoder we confirm that our preparation protocol outperforms the standard protocol in terms of both threshold and logical error rate in the fault-tolerant regime where measurements are unreliable and at experimentally realistic biases. We also discuss how our state-preparation protocol can be inverted for similar fragile-boundary-mitigated logical-state measurement.

Ultracold polar molecules are promising for quantum information processing and searches for physics beyond the standard model. Laser cooling to ultracold temperatures is an established technique for trapped diatomic and triatomic molecules. Further cooling of the molecules to near the motional ground state is crucial for reducing various dephasings in quantum and precision applications. In this work, we demonstrate Raman sideband cooling (RSC) of CaF molecules in optical tweezers to near their motional ground state, with average motional occupation quantum numbers of ${\overline{n}}_x=0.16(12)$, ${\overline{n}}_y=0.17(17)$ (radial directions), and ${\overline{n}}_z=0.22(16)$ (axial direction), and a 3-D motional-ground-state probability of $54\pm 18\%$ of the molecules that survive the RSC. This process paves the way to increase molecular coherence times in optical tweezers for robust quantum computation and simulation applications.

Deep learning algorithms demonstrate a surprising ability to learn high-dimensional tasks from limited examples. This is commonly attributed to the depth of neural networks, enabling them to build a hierarchy of abstract, low-dimensional data representations. However, how many training examples are required to learn such representations remains unknown. To quantitatively study this question, we introduce the random hierarchy model: a family of synthetic tasks inspired by the hierarchical structure of language and images. The model is a classification task where each class corresponds to a group of high-level features, chosen among several equivalent groups associated with the same class. In turn, each feature corresponds to a group of subfeatures chosen among several equivalent groups and so on, following a hierarchy of composition rules. We find that deep networks learn the task by developing internal representations invariant to exchanging equivalent groups. Moreover, the number of data required corresponds to the point where correlations between low-level features and classes become detectable. Overall, our results indicate how deep networks overcome the curse of dimensionality by building invariant representations and provide an estimate of the number of data required to learn a hierarchical task.

We present a systematic framework to classify symmetry-enriched topological quantum spin liquids in two spatial dimensions. This framework can deal with all topological quantum spin liquids, which may be either Abelian or non-Abelian and chiral or nonchiral. It can systematically treat a general symmetry, which may include both lattice symmetry and internal symmetry, may contain antiunitary symmetry, and may permute anyons. The framework applies to all types of lattices and can systematically distinguish different lattice systems with the same symmetry group using their quantum anomalies, which are sometimes known as Lieb-Schultz-Mattis anomalies. We apply this framework to classify $\mathrm{U}(1{)}_{2N}$ chiral states and non-Abelian ${\mathrm{Ising}}^{(\nu )}$ states enriched by a $p6\times \mathrm{SO}(3)$ or $p4\times \mathrm{SO}(3)$ symmetry and ${\mathbb{Z}}_N$ topological orders and $\mathrm{U}(1{)}_{2N}\times \mathrm{U}(1{)}_{-2N}$ topological orders enriched by a $p6m\times \mathrm{SO}(3)\times {\mathbb{Z}}_2^T$, $p4m\times \mathrm{SO}(3)\times {\mathbb{Z}}_2^T$, $p6m\times {\mathbb{Z}}_2^T$, or $p4m\times {\mathbb{Z}}_2^T$ symmetry, where $p6$, $p4$, $p6m$, and $p4m$ are lattice symmetries while SO(3) and ${\mathbb{Z}}_2^T$ are spin rotation and time-reversal symmetries, respectively. In particular, we identify symmetry-enriched topological quantum spin liquids that are not easily captured by the usual parton-mean-field approach, including examples with the familiar ${\mathbb{Z}}_2$ topological order.