## Abstract

Nonequilibrium quantum many-body systems, which are difficult to study via classical computation, have attracted wide interest. Quantum simulation can provide insights into these problems. Here, using a programmable quantum simulator with 16 all-to-all connected superconducting qubits, we investigate the dynamical phase transition in the Lipkin-Meshkov-Glick model with a quenched transverse field. Clear signatures of dynamical phase transitions, merging different concepts of dynamical criticality, are observed by measuring the nonequilibrium order parameter, nonlocal correlations, and the Loschmidt echo. Moreover, near the dynamical critical point, we obtain a spin squeezing of −7.0 ± 0.8 dB, showing multipartite entanglement, useful for measurements with precision fivefold beyond the standard quantum limit. On the basis of the capability of entangling qubits simultaneously and the accurate single-shot readout of multiqubit states, this superconducting quantum simulator can be used to study other problems in nonequilibrium quantum many-body systems, such as thermalization, many-body localization, and emergent phenomena in periodically driven systems.

## INTRODUCTION

Quantum simulation uses a controllable quantum system to mimic complex systems or solve intractable problems (*1*, *2*). Emergent phenomena in out-of-equilibrium quantum many-body systems (*3*), e.g., thermalization (*4*) versus localization (*5*), and time crystals (*6*), have all been recently studied using quantum simulation. Recently, the dynamical phase transition (DPT) and the nonequilibrium phase transition in transient time scales have been theoretically studied in the transverse-field Ising model with all-to-all interactions (*7*–*9*). These two transitions can be characterized by the nonequilibrium order parameter (*7*–*10*) and the Loschmidt echo associated with the Lee-Yang-Fisher zeros in statistical mechanics (*11*), respectively. Moreover, recent experimental progress has allowed the controllable simulation of these exotic phenomena with cold atoms (*12*, *13*) and trapped ions (*14*, *15*). Yet, experimental explorations for the dynamics of entanglement, as a valuable resource in quantum information processing, remain limited in the presence of a DPT.

In our experiments, applying a sudden change of the transverse field with a controllable strength, we drive the system, initially in its ground state, out of equilibrium. Accurate single-shot readout techniques enable us to synchronously record the dynamics of all qubits and to observe essential signatures of DPTs and spin squeezing from the dynamical criticality in the Lipkin-Meshkov-Glick (LMG) model. It is worth mentioning that our experimental system is a 16-qubit device featuring all-to-all connectivity, which complements the type of superconducting circuits used in other simulations of many-body physics (*16*–*20*), where neighboring couplings dominate. The presence of long-range interactions is essential for realizing the LMG model.

This work presents a systematic quantum simulation of DPTs with two different concepts, providing evidence of the relation between the nonequilibrium order parameter and the Loschmidt echo. We verify entanglement in spin-squeezed states generated from dynamical criticality, directly observing squeezing of −7.0 ± 0.8 dB for 16 qubits.

## RESULTS

Our quantum simulator is a superconducting circuit with 20 fully controllable transmon qubits capacitively coupled to a resonator bus ℛ (Fig. 1A). Sixteen qubits (*Q*_{1} to *Q*_{16}), with *XY*-control lines, are selected to perform experiments (see Materials and Methods). The resonant frequency of ℛ is fixed at about 5.51 GHz, while the qubit frequencies are individually tunable via their *Z*-control lines, enabling us to engineer the qubit-qubit interactions induced by ℛ. We detune all 16 qubits from ℛ by, e.g., Δ/2π≃ −450 MHz to switch on the resonator-mediated interactions between two arbitrary qubits (*21*). Simultaneously, identical resonant microwave drives, with a magnitude of *h ^{x}*, are imposed on all qubits to generate the local transverse fields for the control of a DPT (Fig. 1C). To ensure the uniformity of the local fields, the cross-talk effects of microwave pulses have been precisely corrected, and the microwave phase has been calibrated (see the Supplementary Materials). The effective Hamiltonian of the quenched system is

*N*= 16,

*g*represents the coupling strength between ℛ and

_{j}*Q*, and

_{j}*g*/Δ is the resonator-induced virtual coupling strength, which acts as a dominant part (parameters are shown in the Supplementary Materials). Because the values of λ

_{i}g_{j}*are nearly the same for most pairs of qubits and do not decay over distance ∣*

_{ij}*i*−

*j*∣ (Fig. 1B), the quenched system can be reasonably approximated by the LMG model, whose Hamiltonian is

*h*(see Materials and Methods). The LMG model was first introduced in nuclear physics (

^{x}*22*) and then used to describe two-mode Bose-Einstein condensates (

*23*). Recent studies (

*7*–

*10*,

*14*) have shown that

*H*

_{LMG}has a dynamical critical point separating the dynamical paramagnetic phase (DPP) and the dynamical ferromagnetic phase (DFP) with and without a global ℤ

_{2}symmetry, respectively.

First, we show that our programmable superconducting qubits can simulate and verify the DPT by measuring the magnetization and the spin correlation. The system is initialized at the eigenstate ∣00…0〉 of *H*_{1} with *h ^{x}* = 0, where ∣0〉 denotes the ground state of a qubit. Then, we quench the system by suddenly adding a transverse field and monitor its dynamics from

*t*= 0 to 600 ns. With the precise full control and the high-fidelity single-shot readout of each qubit, we are able to omnidirectionally track the evolutions of the average magnetization

*x*,

*y*,

*z*axes for different strengths of the quenched transverse fields, with α ∈ {

*x*,

*y*,

*z*}. By depicting the trajectory of the Bloch vector

*h*/2π ⋍ 2 MHz, 〈σ

^{x}*(*

^{z}*t*)〉 exhibits a slow relaxation (Fig. 2B). However, given a strong transverse field, e.g.,

*h*/2π ⋍ 8 MHz, 〈σ

^{x}*(*

^{z}*t*)〉 exhibits a large oscillation at an early time and approaches zero in the long-time limit (Fig. 2B). In Fig. 2C, we show the behavior of the time-averaged magnetization

*h*. For large

^{x}*h*,

^{x}Figure 2E shows the averaged spin correlation function*h ^{x}* with a final time

*t*= 600 ns, where the DPT is characterized by the local minimum of two-spin correlations. We can observe the critical behaviors of

_{f}Another perspective on dynamical criticality is based on the Loschmidt echo, defined as L(*t*)=∣⟨00…0∣^{2}, where the time *t*, satisfying L(*t*) = 0, is a Lee-Yang-Fisher zero. The zero will cause the non-analytical behavior of the rate function *r*(*t*) = − *N*^{−1} log [L(*t*)], regarded as the complex-plane generalization of the free-energy density (*11*, *24*). Recent numerical studies (*7*, *9*, *24*, *25*) have revealed that the existence of Lee-Yang-Fisher zeros closely relates to the DPT between the DFP and the DPP in long-range interacting systems. In Fig. 3A, we show the distinct behaviors of the Loschmidt echo in different dynamical phases. In the DPP (*h ^{x}*/2π ⋍ 8 MHz), the Loschmidt echo decays rapidly to a near-zero value, related to the occurrence of the non-analyticity of the rate function

*r*(

*t*). A clearer signature can be seen from the first minimum of the Loschmidt echo

*h*(Fig. 3B). In the Supplementary Materials, we clarify the reason of choosing

^{x}*r*(

*t*) using the real parameters of our quantum simulator are also presented as reference. The direct observation of the non-analytical points of the rate function

*r*(

*t*), as a diagnostic signature of the dynamical criticality, deserves further experimental investigations.

In addition to demonstrating a DPT, the LMG model is also useful for generating spin-squeezed states with twist-and-turn dynamics (*26*, *27*). Near the equilibrium critical point, spin squeezing can be achieved, originating from quantum fluctuations, according to the Heisenberg uncertainty principle (*28*). Similarly, we show that spin-squeezed states can also be generated from dynamical criticality. During the dynamics of the quenched Hamiltonian Eq. 1, we can visualize the spin-squeezed state by measuring the quasidistribution *Q*-function (*29*) *Q*(θ, ϕ) ∝ ⟨θ, ϕ∣ρ(*t*)∣θ, ϕ⟩, where *z* axis for each qubit before the joint readout. The experimental and numerical data of *Q*(θ, ϕ) are compared in Fig. 4A, which show spin squeezing with a large strength of the external field, due to stronger quantum fluctuations in the DPP (see also Fig. 2C).

We also measured the time-evolved spin-squeezing parameter (*26*) (see the Supplementary Materials)^{2} < 1, as a sufficient condition for particle entanglement (*30*, *31*), occurs in the time interval *t* ≲ 46 ns when *h ^{x}*/2π ⋍ 3 MHz, and for

*t*≲ 38 ns when

*h*/2π ⋍ 6 MHz. The minimum spin-squeezing parameter over time,

^{x}*h*is shown in Fig. 4C, where the minimum value

^{x}*N*

^{−2/3}, of the squeezing parameter for an

*N*-body one-axis twisting Hamiltonian (

*30*), our 16-qubit system achieves a spin-squeezing parameter satisfying

## DISCUSSION

We have presented clear signatures and entanglement behaviors of the DPT in the LMG model with a superconducting quantum simulator featuring all-to-all connectivity, including the nonequilibrium order parameter, Loschmidt echo, and spin squeezing. On the basis of its high degree of controllability, precise measurement, and long decoherence time, our platform with all-to-all connectivity is powerful for generating multipartite entanglement (*29*, *32*) and investigating nontrivial properties of out-of-equilibrium quantum many-body systems, such as many-body localization (*33*, *34*), quantum chaos in Floquet systems (*35*), and quantum annealing (*36*).

## MATERIALS AND METHODS

### Device information and system Hamiltonian

The device used here consists of 20 frequency-tunable superconducting qubits capacitively coupled to a central resonator bus. It is the same circuit presented in (*29*), where more details about the device, the qubit manipulation, and the readout can be found. In table S1, we present the characteristics for the quantum simulator involving 16 of the 20 qubits, with *XY*-control lines, which have been relabeled in the experiments.

The unused four qubits in this device, without *XY*-control lines, are detuned far off resonance from the other 16 qubits to avoid interacting with them during the experiments. Thus, they will not be included in the following descriptions. The system Hamiltonian, without applying external transverse fields, can be written as_{ℛ} and ω* _{j}* represent the fixed resonant frequency and the tunable frequency of

*Q*, respectively, while

_{j}*g*is the coupling strength between the

_{j}*Q*and resonator bus. The magnitude of the cross-talk coupling between

_{j}*Q*and

_{i}*Q*beyond the resonator-induced virtual coupling is denoted as

_{j}*g*/Δ being the magnitude of the resonator-mediated coupling between

_{i}g_{j}*Q*and

_{i}*Q*. It acts as a dominant part of the qubit-qubit interaction terms, because the cross-talk coupling

_{j}*for all the combinations of pairs of qubits. The individually controllable amplitude and the phase of the microwave drive on each*

_{ij}*Q*are represented by

_{j}*, respectively. In our experiments, we set the uniform amplitude and phase for all qubits, leading to the Hamiltonian in Eq. 1. To ensure this uniformity, the calibration process for the microwave drives is described in the Supplementary Materials.*

_{j}### Relation between the quantum simulator and the LMG model

Our device can be described by the Hamiltonian in Eq. 1. With uniform couplings *J* ≡ *N*λ. The second term can be directly rewritten as *h ^{x}*. According to [

*S*

^{2},

*S*

^{α}] = 0 (α ∈ {

*x*,

*y*,

*z*}), and the fact that the initial state ∣00…0〉 is an eigenstate of

*S*

^{2}, we have

*H*

_{1}can be approximately expressed as the ones of the LMG model

The location of the DPT critical point of the LMG model is μ* _{c}* = ∣

*J*∣/2, leading to

## SUPPLEMENTARY MATERIALS

Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/6/25/eaba4935/DC1

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is **not** for commercial advantage and provided the original work is properly cited.

## REFERENCES AND NOTES

**Acknowledgments:**We thank C. Song, Q. Guo, Z. Wang, and X. Zhang for technical support. Devices were made at the Nanofabrication Facilities at the Institute of Physics in Beijing and National Centre for Nanoscience and Technology in Beijing. The experiment was performed on the quantum computing platform at Zhejiang University.

**Funding:**This work was supported by the National Basic Research Program of China (grant nos. 2016YFA0302104, 2016YFA0300600, and 2017YFA0304300), the National Natural Science Foundation of China (grants nos. 11934018, 11725419, and 11904393), the Strategic Priority Research Program of Chinese Academy of Sciences (grant no. XDB28000000), the China Postdoctoral Science Foundation (grant no. 2018 M640055), the Zhejiang Province Key Research and Development Program (grant no. 2020C01019), the ARO (grant no. W911NF-18-1-0358), the JST Q-LEAP program, the JST CREST (grant no. JPMJCR1676), the JSPS Kakenhi (grant no. JP20H00134), the JSPS Postdoctoral Fellowship (grant no. P19326), FQXi, and the NTT PHI Lab.

**Author contributions:**D.Z., H.F., and H.W. conceived the research. K.X., Z.-H.S., Y.-R.Z., H.F., and H.W. designed the experiment. H.L. and D.Z. fabricated the device. K.X. and W.L. performed the experiment. K.X. and Z.-H.S. did numerical simulations. K.X., Z.-H.S., Y.-R.Z., and H.F. analyzed the results. K.X., Z.-H.S., Y.-R.Z., F.N., H.F., and H.W. wrote the paper. All authors contributed to the experimental setup, discussions of the results, and development of the manuscript.

**Competing interests:**The authors declare that they have no competing interests.

**Data and materials availability:**All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.

- Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).