We use resonant inelastic x-ray scattering (RIXS) at the $\mathrm{Fe}\text{−}{\mathrm{L}}_3$ edge to study the spin excitations of uniaxial-strained and unstrained ${\mathrm{FeSe}}_{1-x}{\mathrm{S}}_x$ ($0\le x\le 0.21$) samples. The measurements on unstrained samples reveal dispersive spin excitations in all doping levels, which show only minor doping dependence in energy dispersion, lifetime, and intensity, indicating that high-energy spin excitations are only marginally affected by sulfur doping. RIXS measurements on uniaxial-strained samples reveal that the high-energy spin-excitation anisotropy observed previously in FeSe is also present in the doping range $0<x\le 0.21$ of ${\mathrm{FeSe}}_{1-x}{\mathrm{S}}_x$. The spin-excitation anisotropy persists to a high temperature up to $T>200\mathrm{K}$ in $x=0.18$ and reaches a maximum around the nematic quantum critical doping ($x_c\approx 0.17$). Since the spin-excitation anisotropy directly reflects the existence of nematic spin correlations, our results indicate that high-energy nematic spin correlations pervade the regime of nematicity in the phase diagram and are enhanced by the nematic quantum criticality. These results emphasize the essential role of spin fluctuations in driving electronic nematicity and highlight the capability of uniaxial strain in tuning spin excitations in quantum materials hosting strong magnetoelastic coupling and electronic nematicity.

A $(2+1)\mathrm{D}$ topologically ordered phase may or may not have a gappable edge, even if its chiral central charge $c_{-}$ is vanishing. Recently, it was discovered that a quantity regarded as a “higher” version of chiral central charge gives a further obstruction beyond $c_{-}$ to gapping out the edge. In this Letter, we show that the higher central charges can be characterized by the expectation value of the partial rotation operator acting on the wave function of the topologically ordered state. This allows us to extract the higher central charge from a single wave function, which can be evaluated on a quantum computer. Our characterization of the higher central charge is analytically derived from the modular properties of edge conformal field theory, as well as the numerical results with the $\nu =1/2$ bosonic Laughlin state and the non-Abelian gapped phase of the Kitaev honeycomb model, which corresponds to $\mathrm{U}(1{)}_2$ and Ising topological order, respectively. The Letter establishes a numerical method to obtain a set of obstructions to the gappable edge of $(2+1)\mathrm{D}$ bosonic topological order beyond $c_{-}$, which enables us to completely determine if a $(2+1)\mathrm{D}$ bosonic Abelian topological order has a gappable edge or not. We also point out that the expectation values of the partial rotation on a single wave function put a constraint on the low-energy spectrum of the bulk-boundary system of $(2+1)\mathrm{D}$ bosonic topological order, reminiscent of the Lieb-Schultz-Mattis-type theorems.

We study a model of a hole-doped collinear Ising antiferromagnet on the honeycomb lattice as a route toward realization of subsystem symmetry. We find nearly exact conservation of dipole symmetry verified both numerically with exact diagonalization on finite clusters and analytically with perturbation theory. The emergent symmetry forbids the motion of single holes—or fractons—but allows hole pairs—or dipoles—to move freely along a one-dimensional line, the antiferromagnetic direction, of the system; in the transverse direction both fractons and dipoles are completely localized. This presents a realization of a “unidirectional” subsystem symmetry. By studying interactions between dipoles, we argue that the subsystem symmetry is likely to continue to persist up to finite (but probably small) hole concentrations.

We numerically model a two-dimensional active nematic confined by a periodic array of fixed obstacles. Even in the passive nematic, the appearance of topological defects is unavoidable due to planar anchoring by the obstacle surfaces. We show that a vortex lattice state emerges as activity is increased, and that this lattice may be tuned from “ferromagnetic” to “antiferromagnetic” by varying the gap size between obstacles. We map the rich variety of states exhibited by the system as a function of distance between obstacles and activity, including a pinned defect state, motile defects, the vortex lattice, and active turbulence. We demonstrate that the flows in the active turbulent phase can be tuned by the presence of obstacles, and explore the effects of a frustrated lattice geometry on the vortex lattice phase.

The composition of cellular metabolism is different across species. Empirical data reveal that bacterial species contain similar numbers of metabolic reactions but that the cross-species popularity of reactions is so heterogenous that some reactions are found in all the species while others are in just few species, characterized by a power-law distribution with the exponent one. Introducing an evolutionary model concretizing the stochastic recruitment of chemical reactions into the metabolism of different species at different times and their inheritance to descendants, we demonstrate that the exponential growth of the number of species containing a reaction and the saturated recruitment rate of brand-new reactions lead to the empirically identified power-law popularity distribution. Furthermore, the structural characteristics of metabolic networks and the species’ phylogeny in our simulations agree well with empirical observations.

Dual-unitary circuits are a class of quantum systems for which exact calculations of various quantities are possible, even for circuits that are nonintegrable. The array of known exact results paints a compelling picture of dual-unitary circuits as rapidly thermalizing systems. However, in this Letter, we present a method to construct dual-unitary circuits for which some simple initial states fail to thermalize, despite the circuits being “maximally chaotic,” ergodic, and mixing. This is achieved by embedding quantum many-body scars in a circuit of arbitrary size and local Hilbert space dimension. We support our analytic results with numerical simulations showing the stark contrast in the rate of entanglement growth from an initial scar state compared to nonscar initial states. Our results are well suited to an experimental test, due to the compatibility of the circuit layout with the native structure of current digital quantum simulators.

Measurements are able to fundamentally affect quantum dynamics. We here show that a continuously measured quantum many-body system can undergo a spontaneous transition from asynchronous stochastic dynamics to noise-free stable synchronization at the level of single trajectories. We formulate general criteria for this quantum phenomenon to occur and demonstrate that the number of synchronized realizations can be controlled from none to all. We additionally find that ergodicity is typically broken, since time and ensemble averages may exhibit radically different synchronization behavior. We further introduce a quantum type of multiplexing that involves individual trajectories with distinct synchronization frequencies. Measurement-induced synchronization appears as a genuine nonclassical form of synchrony that exploits quantum superpositions.

We study experimentally and numerically the noisy evolution of multipartite entangled states, focusing on superconducting qubit devices accessible via the cloud. We find that a valid modeling of the dynamics requires one to properly account for coherent frequency shifts, caused by stochastic charge-parity fluctuations. We introduce an approach modeling the charge-parity splitting using an extended Markovian environment. This approach is numerically scalable to tens of qubits, allowing us to simulate efficiently the dissipative dynamics of some large multiqubit states. Probing the continuous-time dynamics of increasingly larger and more complex initial states with up to 12 coupled qubits in a ring-graph state, we obtain a good agreement of the experiments and simulations. We show that the underlying many-body dynamics generate decays and revivals of stabilizers, which are used extensively in the context of quantum error correction. Furthermore, we demonstrate the mitigation of 2-qubit coherent interactions (crosstalk) using tailored dynamical decoupling sequences. Our noise model and the numerical approach can be valuable to advance the understanding of error correction and mitigation and invite further investigations of their dynamics.

Quantum control is a ubiquitous research field that has enabled physicists to delve into the dynamics and features of quantum systems, delivering powerful applications for various atomic, optical, mechanical, and solid-state systems. In recent years, traditional control techniques based on optimization processes have been translated into efficient artificial intelligence algorithms. Here, we introduce a computational method for optimal quantum control problems via physics-informed neural networks (PINNs). We apply our methodology to open quantum systems by efficiently solving the state-to-state transfer problem with high probabilities, short-time evolution, and using low-energy consumption controls. Furthermore, we illustrate the flexibility of PINNs to solve the same problem under changes in physical parameters and initial conditions, showing advantages in comparison with standard control techniques.

We investigate the detectability of gravitational waves that have been lensed by a spinless stellar-mass black hole, with respect to the advanced LIGO. By solving the full relativistic linear wave equations in the spacetime of a Schwarzschild black hole, we find that the strong gravity can create unique signals in the lensed waveform, particularly during the merger and ring-down stages. The differences in terms of fitting factor between the lensed waveform and best-fitted unlensed general relativity template with spin precession and higher-order multipoles are greater than 5% for the lens black hole mass within $70M_{\odot }<M_{\mathrm{lens}}<133.33M_{\odot }$ under advanced LIGO’s sensitivity. This is up to 5 times more detectable than the previous analysis based on the weak field approximation for a point mass and covers most part of the black hole mass gap predicted by stellar evolution theory. Based on Bayesian inference, the lensing feature can be distinguished with a signal-to-noise ratio of 12.5 for $M_{\mathrm{lens}}=70M_{\odot }$ and 19.2 for $M_{\mathrm{lens}}=250M_{\odot }$, which is attainable for advanced LIGO.

Surface polaritons have proven to be uniquely capable of controlling light-matter interactions. Here we explore surface magnon polaritons in low-loss ferrimagnetic semiconductors, with a focus on their topological phases. We propose several surface magnon polariton devices, including microwave resonators that can strongly enhance magnetic fields and low-loss interconnects joining waveguides with vastly different impedances. Our work can facilitate the exploration of topological phases in polaritons and the development of topological microwave devices for quantum sensing and information processing.

Topological photonic states provide intriguing strategies for robust light manipulations, however, it remains challenging to perfectly excite these topological eigenstates due to their complicated mode profiles. In this work, we propose to realize the exact eigenmode of the topological edge states by supersymmetric (SUSY) structures. By adiabatically transforming the SUSY partner to its main topological structure, the edge modes can be perfectly excited with simple single-site input. We experimentally verify our strategy in integrated silicon waveguides in telecommunication wavelength, showing a broad working bandwidth. Moreover, a shortcut-to-adiabaticity strategy is further applied to speed up the adiabatic pump process by inverse-design approaches, thus enabling fast mode evolutions and leading to reduced device size. Our method is universal and beneficial to the topology-based or complex eigenmodes systems, ranging from photonics and microwaves to cold atoms and acoustics.

Recent studies have found that fluctuations of magnetization transfer in integrable spin chains violate the central limit property. Here, we revisit the problem of anomalous counting statistics in the Landau-Lifshitz field theory by specializing to two distinct anomalous regimes featuring a dynamical critical point. By performing optimized numerical simulations using an integrable space-time discretization, we extract the algebraic growth exponents of time-dependent cumulants which attain their threshold values. The distinctly non-Gaussian statistics of magnetization transfer in the easy-axis regime is found to converge toward the universal distribution of charged single-file systems. At the isotropic point, we infer a weakly non-Gaussian distribution, corroborating the view that superdiffusive spin transport in integrable spin chains does not belong to any known dynamical universality class.

Classical and quantum systems are used to simulate the Ising Hamiltonian, an essential component in large-scale optimization and machine learning. However, as the system size increases, devices like quantum annealers and coherent Ising machines face an exponential drop in their success rate. Here, we introduce a novel approach involving high-dimensional embeddings of the Ising Hamiltonian and a technique called “dimensional annealing” to counteract the decrease in performance. This approach leads to an exponential improvement in the success rate and other performance metrics, slowing down the decline in performance as the system size grows. A thorough examination of convergence dynamics in high-performance computing validates the new methodology. Additionally, we suggest practical implementations using technologies like coherent Ising machines, all-optical systems, and hybrid digital systems. The proposed hyperscaling heuristics can also be applied to other quantum or classical Ising devices by adjusting parameters such as nonlinear gain, loss, and nonlocal couplings.

Energy-energy correlators (EECs) are promising observables to study the dynamics of jet evolution in the quark-gluon plasma (QGP) through its imprint on angular scales in the energy flux of final-state particles. We carry out the first complete calculation of EECs using realistic simulations of high-energy heavy-ion collisions and dissect the different dynamics underlying the final distribution through analyses of jet propagation in a uniform medium. The EECs of $\gamma$-jets in heavy-ion collisions are found to be enhanced by the medium response from elastic scatterings instead of induced gluon radiation at large angles. In the meantime, EECs are suppressed at small angles due to energy loss and transverse momentum broadening of jet shower partons. These modifications are further shown to be sensitive to the angular scale of the in-medium interaction, as characterized by the Debye screening mass. Experimental verification and measurement of such modifications will shed light on this scale and the short-distance structure of the QGP in heavy-ion collisions.