Explore the vast range of titles available.

Simple models of interacting spins have an important role in physics. They capture the properties of many magnetic materials, but also extend to other systems, such as bosons and fermions in a lattice, gauge theories, high-temperature superconductors, quantum spin liquids, and systems with exotic particles such as anyons and Majorana fermions 1 , 2 . To study and compare these models, a versatile platform is needed. Realizing such systems has been a long-standing goal in the field of ultracold atoms. So far, spin transport has only been studied in systems with isotropic spin–spin interactions 3 – 12 . Here we realize the Heisenberg model describing spins on a lattice, with fully adjustable anisotropy of the nearest-neighbour spin–spin couplings (called the XXZ model). In this model we study spin transport far from equilibrium after quantum quenches from imprinted spin-helix patterns. When spins are coupled only along two of three possible orientations (the XX model), we find ballistic behaviour of spin dynamics, whereas for isotropic interactions (the XXX model), we find diffusive behaviour. More generally, for positive anisotropies, the dynamics ranges from anomalous superdiffusion to subdiffusion, whereas for negative anisotropies, we observe a crossover in the time domain from ballistic to diffusive transport. This behaviour is in contrast with expectations from the linear-response regime and raises new questions in understanding quantum many-body dynamics far away from equilibrium. Spin transport far from equilibrium is studied in a Heisenberg model with adjustable anisotropy realized with coupled ultracold 7 Li atoms, and different dynamical regimes are found for positive and negative anisotropies.

The Lie-Trotter formula, together with its higher-order generalizations, provides a direct approach to decomposing the exponential of a sum of operators. Despite significant effort, the error scaling of such product formulas remains poorly understood. We develop a theory of Trotter error that overcomes the limitations of prior approaches based on truncating the Baker-Campbell-Hausdorff expansion. Our analysis directly exploits the commutativity of operator summands, producing tighter error bounds for both real- and imaginary-time evolutions. Whereas previous work achieves similar goals for systems with geometric locality or Lie-algebraic structure, our approach holds, in general. We give a host of improved algorithms for digital quantum simulation and quantum Monte Carlo methods, including simulations of second-quantized plane-wave electronic structure,k-local Hamiltonians, rapidly decaying power-law interactions, clustered Hamiltonians, the transverse field Ising model, and quantum ferromagnets, nearly matching or even outperforming the best previous results. We obtain further speedups using the fact that product formulas can preserve the locality of the simulated system. Specifically, we show that local observables can be simulated with complexity independent of the system size for power-law interacting systems, which implies a Lieb-Robinson bound as a by-product. Our analysis reproduces known tight bounds for first- and second-order formulas. Our higher-order bound overestimates the complexity of simulating a one-dimensional Heisenberg model with an even-odd ordering of terms by only a factor of 5, and it is close to tight for power-law interactions and other orderings of terms. This result suggests that our theory can accurately characterize Trotter error in terms of both asymptotic scaling and constant prefactor.

Some measures of genuine multipartite entanglement in the three-qubit XYZ Heisenberg model are reported for a system that is embedded in a uniform magnetic field. A quantitative characterization of these measures is accomplished to determine the optimal parameters that provide maximal entanglement.

Supersolid, an exotic quantum state of matter that consists of particles forming an incompressible solid structure while simultaneously showing superfluidity of zero viscosity 1 , is one of the long-standing pursuits in fundamental research 2 , 3 . Although the initial report of 4 He supersolid turned out to be an artefact 4 , this intriguing quantum matter has inspired enthusiastic investigations into ultracold quantum gases 5 , 6 , 7 – 8 . Nevertheless, the realization of supersolidity in condensed matter remains elusive. Here we find evidence for a quantum magnetic analogue of supersolid—the spin supersolid—in the recently synthesized triangular-lattice antiferromagnet Na 2 BaCo(PO 4 ) 2 (ref. 9 ). Notably, a giant magnetocaloric effect related to the spin supersolidity is observed in the demagnetization cooling process, manifesting itself as two prominent valley-like regimes, with the lowest temperature attaining below 100 mK. Not only is there an experimentally determined series of critical fields but the demagnetization cooling profile also shows excellent agreement with the theoretical simulations with an easy-axis Heisenberg model. Neutron diffractions also successfully locate the proposed spin supersolid phases by revealing the coexistence of three-sublattice spin solid order and interlayer incommensurability indicative of the spin superfluidity. Thus, our results reveal a strong entropic effect of the spin supersolid phase in a frustrated quantum magnet and open up a viable and promising avenue for applications in sub-kelvin refrigeration, especially in the context of persistent concerns about helium shortages 10 , 11 . Evidence for a quantum magnetic analogue of a supersolid appears in a recently synthesized antiferromagnet showing a strong magnetocaloric effect of the spin supersolid phase with potential for applications in sub-kelvin refrigeration.

Topology in quantum many-body systems has profoundly changed our understanding of quantum phases of matter. The model that has played an instrumental role in elucidating these effects is the antiferromagnetic spin-1 Haldane chain 1 , 2 . Its ground state is a disordered state, with symmetry-protected fourfold-degenerate edge states due to fractional spin excitations. In the bulk, it is characterized by vanishing two-point spin correlations, gapped excitations and a characteristic non-local order parameter 3 , 4 . More recently it has been understood that the Haldane chain forms a specific example of a more general classification scheme of symmetry-protected topological phases of matter, which is based on ideas connected to quantum information and entanglement 5 – 7 . Here, we realize a finite-temperature version of such a topological Haldane phase with Fermi–Hubbard ladders in an ultracold-atom quantum simulator. We directly reveal both edge and bulk properties of the system through the use of single-site and particle-resolved measurements, as well as non-local correlation functions. Continuously changing the Hubbard interaction strength of the system enables us to investigate the robustness of the phase to charge (density) fluctuations far from the regime of the Heisenberg model, using a novel correlator. A ladder-like arrangement of an ultracold gas of lithium atoms trapped in an optical lattice enables the observation of a symmetry-protected topological phase.

Differentiable programming is a fresh programming paradigm which composes parameterized algorithmic components and optimizes them using gradient search. The concept emerges from deep learning but is not limited to training neural networks. We present the theory and practice of programming tensor network algorithms in a fully differentiable way. By formulating the tensor network algorithm as a computation graph, one can compute higher-order derivatives of the program accurately and efficiently using automatic differentiation. We present essential techniques to differentiate through the tensor networks contraction algorithms, including numerical stable differentiation for tensor decompositions and efficient backpropagation through fixed-point iterations. As a demonstration, we compute the specific heat of the Ising model directly by taking the second-order derivative of the free energy obtained in the tensor renormalization group calculation. Next, we perform gradient-based variational optimization of infinite projected entangled pair states for the quantum antiferromagnetic Heisenberg model and obtain state-of-the-art variational energy and magnetization with moderate efforts. Differentiable programming removes laborious human efforts in deriving and implementing analytical gradients for tensor network programs, which opens the door to more innovations in tensor network algorithms and applications.

Synthetic quantum systems with interacting constituents play an important role in quantum information processing and in explaining fundamental phenomena in many-body physics. Following impressive advances in cooling and trapping techniques, ensembles of ultracold polar molecules have emerged as a promising platform that combines several advantageous properties 1 – 11 . These include a large set of internal states with long coherence times 12 – 17 and long-range, anisotropic interactions. These features could enable the exploration of intriguing phases of correlated quantum matter, such as topological superfluids 18 , quantum spin liquids 19 , fractional Chern insulators 20 and quantum magnets 21 , 22 . Probing correlations in these phases is crucial to understanding their properties, necessitating the development of new experimental techniques. Here we use quantum gas microscopy 23 to measure the site-resolved dynamics of quantum correlations of polar 23 Na 87 Rb molecules confined in a two-dimensional optical lattice. By using two rotational states of the molecules, we realize a spin-1/2 system with dipolar interactions between particles, producing a quantum spin-exchange model 21 , 22 , 24 , 25 . We study the evolution of correlations during the thermalization process of an out-of-equilibrium spin system for both spatially isotropic and anisotropic interactions. Furthermore, we examine the correlation dynamics of a spin-anisotropic Heisenberg model engineered from the native spin-exchange model by using periodic microwave pulses 26 – 28 . These experiments push the frontier of probing and controlling interacting systems of ultracold molecules, with prospects for exploring new regimes of quantum matter and characterizing entangled states that are useful for quantum computation 29 , 30 and metrology 31 . Experiments demonstrate the powerful capabilities of ultracold molecules to study dynamics in the context of quantum magnetism, and create new possibilities for studying quantum physics with ultracold molecules more broadly.

Pursuing fractionalized particles that do not bear properties of conventional measurable objects, exemplified by bare particles in the vacuum such as electrons and elementary excitations such as magnons, is a challenge in physics. Here we show that a machine-learning method for quantum many-body systems that has achieved state-of-the-art accuracy reveals the existence of a quantum spin liquid (QSL) phase in the region0.49≲J2/J1≲0.54convincingly in spin-1/2frustrated Heisenberg model with the nearest and next-nearest-neighbor exchanges,J1andJ2, respectively, on the square lattice. This is achieved by combining with the cutting-edge computational schemes known as the correlation ratio and level spectroscopy methods to mitigate the finite-size effects. The quantitative one-to-one correspondence between the correlations in the ground state and the excitation spectra found in the present analyses enables the reliable identification and estimation of the QSL and its nature. The spin excitation spectra containing both singlet and triplet gapless Dirac-like dispersions signal the emergence of gapless fractionalized spin-1/2Dirac-type spinons in the distinctive QSL phase. Unexplored critical behavior with coexisting and dual power-law decays of Néel antiferromagnetic and dimer correlations is revealed. The power-law decay exponents of the two correlations differently vary withJ2/J1in the QSL phase and thus have different values except for a single point satisfying the symmetry of the two correlations. The isomorph of excitations with the cuprated-wave superconductors revealed here implies a tight connection between the present QSL and superconductivity. This achievement demonstrates that the quantum-state representation using machine-learning techniques, which had mostly been limited to benchmarks, is a promising tool for investigating grand challenges in quantum many-body physics.

We study the spin-excitation spectrum (dynamic structure factor) of the spin-1/2 square-lattice Heisenberg antiferromagnet and an extended model (the J−Q model) including four-spin interactions Q in addition to the Heisenberg exchange J . Using an improved method for stochastic analytic continuation of imaginary-time correlation functions computed with quantum Monte Carlo simulations, we can treat the sharp (δ -function) contribution to the structure factor expected from spin-wave (magnon) excitations, in addition to resolving a continuum above the magnon energy. Spectra for the Heisenberg model are in excellent agreement with recent neutron-scattering experiments on Cu(DCOO)2·4D2O , where a broad spectral-weight continuum at wave vector q=(π,0) was interpreted as deconfined spinons, i.e., fractional excitations carrying half of the spin of a magnon. Our results at (π,0) show a similar reduction of the magnon weight and a large continuum, while the continuum is much smaller at q=(π/2,π/2) (as also seen experimentally). We further investigate the reasons for the small magnon weight at (π,0) and the nature of the corresponding excitation by studying the evolution of the spectral functions in the J−Q model. Upon turning on the Q interaction, we observe a rapid reduction of the magnon weight to zero, well before the system undergoes a deconfined quantum phase transition into a nonmagnetic spontaneously dimerized state. Based on these results, we reinterpret the picture of deconfined spinons at (π,0) in the experiments as nearly deconfined spinons—a precursor to deconfined quantum criticality. To further elucidate the picture of a fragile (π,0) -magnon pole in the Heisenberg model and its depletion in the J−Q model, we introduce an effective model of the excitations in which a magnon can split into two spinons that do not separate but fluctuate in and out of the magnon space (in analogy to the resonance between a photon and a particle-hole pair in the exciton-polariton problem). The model can reproduce the reduction of magnon weight and lowered excitation energy at (π,0) in the Heisenberg model, as well as the energy maximum and smaller continuum at (π/2,π/2) . It can also account for the rapid loss of the (π,0) magnon with increasing Q and the remarkable persistence of a large magnon pole at q=(π/2,π/2) even at the deconfined critical point. The fragility of the magnons close to (π,0) in the Heisenberg model suggests that various interactions that likely are important in many materials—e.g., longer-range pair exchange, ring exchange, and spin-phonon interactions—may also destroy these magnons and lead to even stronger spinon signatures than in Cu(DCOO)2·4D2O .

A bstract A general analysis of line defect renormalisation group (RG) flows in the ε expansion below d = 4 dimensions is undertaken. The defect beta function for general scalar-fermion bulk theories is computed to next-to-leading order in the bulk couplings. Scalar models as well as scalar-fermion models with various global symmetries in the bulk are considered at leading non-trivial order. Different types of potential infrared (IR) defect conformal field theories (dCFTs) and their RG stability are discussed. The possibility of multiple IR stable dCFTs is realised in specific examples with hypertetrahedral symmetry in the bulk. The one-point function coefficient of the order parameter in the stable IR dCFT of the cubic model is computed at next-to-leading order and compared with that in the IR dCFT of the Heisenberg model.