Arxiv:1806.09269V2 [Quant-Ph] 13 May 2019 Coordinate Momentum Transformation

Observation of a dynamical quantum phase transition by a superconducting qubit simulation

Xue-Yi Guo,1, 2 Chao Yang,1, 2 Yu Zeng,1, 2 Yi Peng,1, 2 He-Kang Li,1, 2 Hui Deng,3 Yi-Rong Jin,1, 4 Shu Chen,1, 2, ∗ Dongning Zheng,1, 2, 4, 5, † and Heng Fan1, 2, 4, 5, ‡ 1Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China 2School of Physical Sciences, UCAS, Beijing 100190, China 3Synergetic Innovation Centre in Quantum Information and Quantum Physics, USTC, Hefei 230026, China 4CAS Center for Excellence in Topological Quantum Computation, UCAS, Beijing 100190, China 5Songshan Lake Materials Laboratory, Dongguan 523808, China (Dated: May 14, 2019) A dynamical quantum phase transition can occur during time evolution of sudden quenched quantum systems across a phase transition. It corresponds to the nonanalytic behavior at a critical time of the rate function of the quantum state return amplitude, analogous to nonanalyticity of the free energy density at the critical temperature in macroscopic systems. A variety of many-body systems can be represented in momentum space as a spin-1/2 state evolving on the Bloch sphere, where each momentum mode is decoupled and thus can be simulated independently by a single qubit. Here, we report the observation of a dynamical quantum phase transition in a superconducting qubit simulation of the quantum quench dynamics of many-body systems. We take the Ising model with a transverse field as an example for demonstration. In our experiment, the spin state, which is initially polarized longitudinally, evolves based on a Hamiltonian with adjustable parameters depending on the momentum and strength of the transverse magnetic field. The time evolving quantum state is read out by state tomography. Evidence of dynamical quantum phase transitions, such as paths of time evolution states on the Bloch sphere, non-analytic behavior of the dynamical free energy and the emergence of Skyrmion lattice in momentum-time space, is observed. The experimental data agrees well with theoretical and numerical calculations. The experiment demonstrates for the first time explicitly the topological invariant, both topologically trivial and non-trivial, for dynamical quantum phase transitions. Our results show that the quantum phase transitions of this class of many-body systems can be simulated successfully with a single qubit by varying certain control parameters over the corresponding momentum range.

I. INTRODUCTION curring at critical temperature. The DQPT is intimately re- lated to quantum phase transitions in many-body systems [20– Quantum simulation can provide insight into quantum and 22, 24–36]. topological phases of matter, the role of entanglement, and Recently, experimental explorations of DQPT have been quantum dynamics [1–19]. It also constitutes one of the basic performed in ion-trap systems [3, 5] and cold atom sys- building blocks of quantum information processors. Systems tems [4, 6] with dozens of individual addressable qubits or a of high dimension or many-body systems can be simulated by cloud of fermionic atoms. Our experiment follows the DQPT the quantum processors with many coherently coupled qubits simulation approach by emulating a corresponding two-band [3, 4]. By increasing the number of qubits, the simulation model separately for each momentum mode with a single may outperform classical machines, and demonstrate quan- qubit. By ranging over the Brillouin zone of momentum tum advantage. On the other hand, a variety of many-body space, the results are equivalent to that of simulating many- systems with a large number of spin-1/2 states can be stud- body systems in space. The finite size effect can be ob- ied in momentum space by a two-band model with decoupled served for a finite number of momenta implemented experi- momentum modes, which is equivalent to a single spin-1/2 mentally. Our experimental system consists of superconduct- state evolving on the Bloch sphere for each mode. This fact ing Xmon qubits, which is one of the most promising plat- also provides a route of quantum simulation with one qubit forms for quantum simulation and quantum computation [37– and variables sweeping over momentum space, by means of 41]. We provide concrete evidence that the DQPT is success- arXiv:1806.09269v2 [quant-ph] 13 May 2019 coordinate momentum transformation. fully simulated. In particular, we demonstrate experimentally In this Letter, we emulate the dynamical quantum phase the topological invariant in DQPT, which was studied recently transition (DQPT) of the many-body systems by a single su- in Refs.[24, 33], and have obtained quantitatively the dynam- perconducting qubit. The DQPT is a phenomenon occurring ical free energy and Skyrmion lattice. in evolving quantum states [20–22] for isolated quantum systems far from equilibrium [23]. It is characterized by the nonanalyticity in dynamical free energy density at a critical time II. THE MODEL AND SCHEME FOR SIMULATION t = tc, which is analogous to traditional phase transitions oc- We begin with a two-band model with Hamiltonian written in momentum space as ∗ [email protected] † X [email protected] H = Ψ†h(k)Ψ , (1) ‡ [email protected] k k k 2 where Ψk denotes a spinor, which is a 2-dimensional column tomography. By ranging over the Brillouin zone of momen- vector formed by the fermion operators. The “first quantized” tum space for each mode k, we can obtain the rate function in Hamiltonian h(k) for k momentum mode takes the form Eq. (5). The occurrence of DQPT can be observed when the rotation path of |φ(k, t)i is a great circle on the Bloch sphere h(k) = d0(k) + d(k) · σ, (2) for mode k = k∗. In this case, |φ(k∗, t)i is orthogonal to the initial state at time t , hφ (k∗)|φ(k∗, t )i = 0, resulting in a non- where σ = (σx, σy, σz) is a vector of Pauli matrices and k is in c i c the Brillouin zone. This model can describe a variety of phys- analytic point of the rate function. For the full regime of k in ically different many-body systems [42–45], see appendix for the Brillouin zone, the time evolutions of states |φ(k, t)i will details. cover the full Bloch sphere when there exists DQPT, other- To study the quench dynamics, we first prepare the sys- wise only less than one half of the Bloch sphere is covered, tem in the ground state of the initial Hamiltonian h (k), i.e. as recently pointed out by our co-authors [24, 33]. This phe- h i i ρ (k) = |φ (k)ihφ (k)| = 1 1 − dˆ (k) · σ . Then with a sudden nomenon is observed, for the first time in experiment, as one i i i 2 i of the signatures in identifying the occurrence of DQPT. It is quench to the final Hamiltonian h (k), which determines d (k) f f actually a direct observation of the topological invariant. by Eq.(2), the state evolves as 1 h i ρ(k, t) = |φ(k, t)ihφ(k, t)| = 1 − dˆ(k, t) · σ , (3) 2 where −itd f (k)·σ itd f (k)·σ dˆ (k, t) · σ = e dˆ i (k) · σ e . (4) This is simply the spin precession on the Bloch sphere, that is, π dˆ i(k) rotating around dˆ f (k) with period . |d f (k)| Now we introduce the rate function of the dynamical free energy, 1 X f (t) = − log |hφ (k)|e−ith f (k)|φ (k)i|2. (5) N i i k The nonanalytic behavior of f (t) corresponds to DQPT, which −ith (k) is associated with zeros of hφi(k)|e f |φi(k)i for at least one ∗ critical momentum k at critical time tc. From the spin preces- ∗ sion picture, it is clear that spin vector di (k ) is perpendicular ∗ ∗ ∗ to the rotation axis d f (k ), where di (k ) , 0, d f (k ) , 0, and π tc repeats with period ∗ . |d f (k )| In this Letter, without loss of generality, we will investigate experimentally the DQPT of the Ising model with a transverse FIG. 1. Photography of qubit chip and external circuitries. (a) is the field, but the approach is applicable to other similar phenom- microscopic photography of our Xmon qubit chip. The red part is the ena of many-body systems. The Hamiltonian of the transverse Xmon qubit. Its frequency can be adjusted by applying DC current field Ising model is through its Z control line. Transmission line coupled to the readout N cavity is to measure the qubit state. Basic information of qubit are X x x z listed in the table. The experiment data of the energy relaxation time HIsing = − σi σi+1 + gσi , (6) T1, dephasing time T2∗ and spin echo dephasing time T2se are also i=1 shown. (b) is the sketch of our experiment circuit setup, the blue part where g is the strength of the field in the Z direction, and the is for Z bias, the green part is for XY control, and the brown part is periodic boundary condition is assumed. There are two phases for readout. (c) The qubit parameters are presented in the table. for this model, the ferromagnetic phase for g < 1, and the paramagnetic phase for g > 1; the phase transition critical point is gc = 1. It is proved that DQPT occurs if and only if the initial Hamiltonian with a gi field and the quenching III. EXPERIMENTAL SETUP Hamiltonian with a g f field belong to different phases [20]. The scheme for simulating DQPT in experiment is as fol- In the experiment we use a single qubit to simulate the dy- lows. We first prepare the initial qubit state |φi(k)i determined namics of the model. Figure 1 is the microscopic photograph by parameter gi for each mode k. By the sudden quench, of the superconducting Xmon qubit chip [37], the external state |φi(k)i evolves as |φ(k, t)i according to the Hamiltonian, circuitries and qubit parameters. In experiment, the Xmon y x h f (k) = (g f − cos k)σ + sin kσ , which depends on param- qubit is biased at its maximum frequency of 6.203 GHz, – eter g f , i.e., the spin vector dˆ i(k) rotates around axis dˆ f (k) also known as the sweet-spot. The measured anharmonicity on the Bloch sphere, see appendix for details. The time evo- is about −236 MHz, the measured energy relaxation time T1 lution state |φ(k, t)i will be read out experimentally by state about 8.3 µs, dephasing time T2∗ about 6.8 µs and spin echo 3

(a)k = 0 (b)k = 0.2p (c)k = 0.4p (d)k = 0.6p (e)k = 0.8p (f) k = p (g)

(h)k = 0 (i)k = 0.2 p (j)k = 0.4 p (k)k = 0.6p (l)k = 0.8 p (m) k = p (n)

FIG. 2. Evolution of states on the Bloch sphere. The state evolves depending on a fixed momentum. The data of evolving path are presented on the Bloch sphere. Here gi = 0.2 is fixed, cases with g f = 1.5 are presented in (a)-(g) in upper panel, cases of g f = 0.5 are presented in (h)-(n) in lower panel. The momenta are chosen to be k = 0, 0.2π, 0.4π, 0.6π, 0.8π, π, presented respectively on up-down pairs of sub-figures, (a,h), (b,i), (c,j),(d,k), (e,l), (f,m). Data for two different cases are summarized together in (g) for g f = 1.5, and (n) for g f = 1.5, respectively. We can find that the whole Bloch sphere in (g) is covered, in contrast in (n), only partial region of the Bloch sphere is covered. We emphasize that the initial state is always prepared on the equator of the Bloch sphere in X-axis.

dephasing time T2se about 11.7 µs. The readout cavity fre- magnetic phase regime. The system is suddenly quenched to quency is about 6.793 GHz, which falls in the dispersive cou- the final Hamiltonian h f (k). Here two different strengths of pling regime. The frequency dispersive shift of the readout the field are chosen, g f = 0.5 and g f = 1.5, corresponding to cavity is κ/2π = −0.697 MHz. the ferromagnetic and paramagnetic phases, respectively. The energy gap of the qubit can be adjusted by an exter- The qubit is first rotated about the Y-axis by a π/2√ mi- nal flux bias. The Xmon qubit is capacitively coupled to a crowave pulse to the superposed state |φii = (|0i + |1i)/ 2. λ/4 coplanar waveguide (CPW) resonator that is coupled to a For a fixed k, a unitary operation based on the final Hamilto- CPW transmission line. In this device, the qubit state is read nian is applied to the initial state |φii as the quantum quench out by the dispersive method via the λ/4 resonator. The op- procedure. We then sweep mode k in the Brillouin zone from tical micrograph of this sample is shown in Fig. 1(a). The 0 to 2π with step length 2π/30. details of chip fabrication and the circuitry are presented in appendix. For each value of k, the state will be rotated for two cy- cles on the Bloch sphere, representing time evolution for two periods. The rotation axis is determined by the quench Hamil- tonian h f (k). The path of the state time evolution is presented IV. TIME EVOLUTION PATHS ON THE BLOCH SPHERE in Fig. 2, where only one cycle of data is presented. In the FOR DQPT figure, each dot represents the evolving state at a fixed time point read out experimentally by state tomography. For exam- Following our experimental scheme, we first prepare the ple, Figs. 2(a) and (h) represent k = 0 for different quenched initial state as the ground state of the Hamiltonian hi(k) for Hamiltonians. We can find that the initial state always stays a fixed mode k, then suddenly quench the system to the final at its original position, because the rotation axis is the X- Hamiltonian h f (k). For convenience, we actually√ always pre- direction determined by the corresponding Hamiltonian. pare the initial state as |φii = (|0i + |1i)/ 2, consequently Figure 2(a-g), 7 sub-figures in the upper panel, represent the the quenched Hamiltonian is changed accordingly. This is be- system is suddenly quenched to g f = 1.5, and the state evo- cause we can perform a rotation to both Hamiltonians, hi(k) lutions on the Bloch sphere for k = 0, 0.2π, 0.4π, 0.6π, 0.8π, π and h f (k), without changing the DQPT results. are presented in the first 6 sub-figures, respectively. All data The quenched quantum state will be read out at a sequence for this case are presented together in Fig. 2(g), where those of time points to obtain the time dependent density matrix k modes are for k ∈ [0, π] constituting a half region in the ρ(k, t). For a full rotation period, we can obtain a circular Brillouin zone. For each mode k, the state starts from |φii in evolution path of the state on the Bloch sphere. The same the original position and evolves like a circle on the Bloch procedure repeats by changing momentum k in the Brillouin sphere. In each cycle of time period, we take 70 time points zone. for state tomography readout. The experimental data are pre- In the experiment, we let gi = 0.2, which is in the ferro- sented as dots on the Bloch sphere, where each dot represents 4

()a average value of 5000 single-shot measurement results. Each 1.5 step of time evolution lasts 15 nanoseconds. Then one circle of period takes 1.05 µs, two circles are also performed exper- 1 imentally, they are within the coherence time. We have also 2 f(t) taken a normalization, | < |φ(k, t)|φ(k, t) > | = 1, at each time 0.5 point, implying pure states are assumed for time evolution, 2 Trρ (k, t) = 1. We take total 30 different momenta k in the 0 Brillouin zone in experiment, the evolution paths are given in 0 0.2 0.4 0.6 0.8 1 t/t Fig. 2(g). 0 ()b Figure 2(h-n), 7 sub-figures in the lower panel, represent 1.2 the case that the system is suddenly quenched to g f = 0.5. 1 Similar conventions are used as those in upper panel. 0.8 The occurrence of DQPT can be directly observed in Fig. 2. f(t) It is obvious that in upper panel of the figure, Fig. 2(a-g), the 0.6 full Bloch sphere is covered by states time evolution paths t c 0.4 shown explicitly in Fig. 2(g). This case is that g f = 1.5 and 0.45 0.5 0.55 0.6 0.65 t/t gi = 0.2 are located in two different phases, so DQPT hap- 0 pens. Since the full Bloch sphere is covered, it is apparent ∗ that there exists a k , the path of the evolving state√ is a great FIG. 3. The dynamical free energy. (a) Dynamical free energies ⊥ for different g f s are presented, while gi = 0.2 is fixed. We take 70 circle resulting in that state |φi i = (|0i − |1i)/ 2, located in the opposite direction of X-axis on the Bloch sphere, can be time points for each period. Error bar represents the deviation of the average value of 5000 single-shot measurements from the fitting reached at a critical time tc, shown in Fig. 2(b). The orthogo- nality leads to zero for overlap between the evolving state |φ⊥i value for the evolution path of 70 points on the Bloch sphere, see i appendix for details. (b) The dynamical free energies near the critical with the initial state |φii, leading to non-analyticity for loga- time τc for different number of momenta implemented in experiment rithm in the rate function (5). These results demonstrate the are given, corresponding to different sizes. The exact results for N → occurrence of DQPT at a critical time tc. In contrast, when ∞ are presented as solid line. Here we take gi = 0.2 and g f = 1.5. gi = 0.2 is quenched to g f = 0.5 but without going across the critical point gc = 1, we can observe in Fig. 2(h-n) that only less than one half of the Bloch sphere is covered for k in the proaches to the theoretical value for N → ∞ and demonstrates Brillouin zone, as summarized in Fig. 2(n). Then no DQPT non-analytic behavior. can happen. Figure 4 shows the emergence of Skyrmion lattice in momentum-time space for DQPT. We define the expectation value as, hdˆ(k, t)i = hφ(k, t)|dˆ i|φ(k, t)i, see also appendix. We V. THE RATE FUNCTION, FINITE SIZE EFFECT AND consider two different cases, g f = 1.5 in Fig. 4(c) and THE SKYRMION LATTICE g f = 0.5 in Fig. 4(d), both start from the initial condition gi = 0. The rate functions are presented respectively on Quantitatively, we can obtain the evolution of dynamical Fig. 4(a,b) for comparison. We find that when gi = 0 and free energy defined in Eq.(5). Figure 3(a) presents the time g f = 1.5 lie in different phases, the emergence of Skyrmion dependent rate functions for different g f , all with initial pa- lattice in momentum-time space can be seen obviously in Fig. rameter gi = 0.2. The experimental data are shown as dots, 4(c), which indicates the nontrivial dynamical Chern number the theoretical results are presented as lines. We can find implying the occurrence of DQPT. The time coordinates of that the rate functions f (t) have sharp peaks at the critical the center of Skyrmion is just the critical time tc. While if time tc for g f = 3.5, 2.5, 1.5, which lead to discontinuity for gi = 0 and g f = 0.5 lie in the same phase, the configuration derivative of f (t) at tc. This phenomenon corresponds to the of Skyrmion lattice does not appear as shown in Fig. 4(d), the DQPT. In comparison, it is obvious that the rate functions for corresponding dynamical Chern number is trivial. There is no g f = 0.9, 0.8, 0.7 are different from cases when g f > 1. The DQPT as shown Fig. 4(b). curves are much more smooth and no sharp peak appears, so no discontinuity is expected for derivative of the rate functions. Thus no DQPT will happen. The results agree well VI. CONCLUSION AND DISCUSSION with theoretical calculations. Figure 3(b) shows the results of different number N of In summary, we simulate successfully the two-band model modes for k, corresponding to size N of the Ising model. So of DQPT for the transverse field Ising model by a single su- experiments are performed for N equally separated momenta perconducting qubit. The DQPT is shown by state evolution for k ∈ [0, 2π]. Here gi = 0.2 and g f = 1.5 are fixed. It can be paths on the Bloch sphere, the dynamical free energy and the found that if t is far away from tc, the dynamical free energy Skyrmion lattice. The critical time of DQPT is quantitatively f (t) is quite close to theoretic value (pink curve) of N → ∞. identified. This approach is applicable in investigating various Near the critical time tc, f (t) is nearly smooth if the size is physical phenomena of the class of free fermionic many-body small, demonstrating finite size effect. As N increases, it ap- systems. The similar experimental scheme can be applied to 5

(a) (b) of China (Grant Nos. 11425419, 11404386, 11674376, 2 0.4 11774406), and Strategic Priority Research Program of Chi- 1.5 0.3 nese Academy of Sciences (Grant No. XDB28000000). 1 0.2 f(t) f(t)

0.5 0.1 APPENDIX A: THE QUBIT DEVICE AND EXTERNAL 0 0 CIRCUITRIES 0 1 2 0 1 2 t/t t/t 0 0 The sample was fabricated using a process involving (c) (d) 2 2 1 electron-beam-lithography (EBL) and double-angle evaporation. In brief, a 100 nm thick Al layer was ﬁrstly deposited 1.5 1.5 0.5 on a 10 × 10 mm sapphire substrate by means of electron-

1 1 0 beam evaporation, followed by EBL and wet etching to pro- k/ k/ duce large structures such as microwave coplanar-waveguide 0.5 0.5 -0.5 resonators/transmission lines, capacitors of Xmon qubit and electric leads. The EPL resist used was ZEP520 and wet etch- 0 0 -1 0 1 2 0 1 2 ing process was carried out using Aluminum Etchant Type A. t/t t/t 0 0 In the next step, the Josephson junctions of qubits were fabricated using the double-angle evaporation process. In this step, FIG. 4. Skyrmion and DQPT. (a) and (b) are respectively dynamical the under cut structure was created using a PMMA-MMA free energies for g f = 1.5 and g f = 0.5, both with gi = 0. The double layer EBL resist following a process similar to that re- number of momenta is N = 30. For each momentum mode k, states ported in Ref. [37]. During the evaporation, the bottom elec- at 2 × 70 = 140 time points are read out. Error bar is deviation of trode was about 30 nm thick while the top electrode was about the average value from the fitting value, see appendix for details. (c) 100 nm thick with intermediate oxidation. and (d) are expectation values of the initial spin operator about the In the measurements, the sample was mounted in an alu- evolving state |φ(k, t) >. The Skyrmion are shown obviously in (c), while no Skyrmion appears in (d). minum alloy sample box which is fixed on the mixing chamber stage of a dilution refrigerator. The temperature of the mixing chamber was below 15 mK during measurements. The readout input microwave lines and qubit XY control lines are simulate temporal topological phenomena by demonstrating heavily attenuated. Lines for qubit dc bias control are filtered that a single qubit is driven by two elliptically polarized peri- using filters (RLC ELECTRONICS F-10-200-R) that func- odic waves [45]. tions as combination of low-pass filter and copper powder fil- In our scheme, phenomena of many-body system are sim- ter. The microwave output signal from the transmission line is ulated by a single qubit at the expenses of repeating exper- amplified (≈ 39 dB) by a cryogenic HEMT amplifier mounted iments by ranging over the momentum space. On the other at the 4 K stage and a room temperature amplifier (≈ 38 dB) hand, besides the two-level system, it is known that the super- before being measured by a home-built heterodyne acquisition conducting Josephson junction can have controllable multiple system shown in Fig. 1(b) in the main text. energy levels. Then, this platform is promising for more simulating applications, such as the PT-symmetric physics, geo- metric quantum logic gates for quantum computation. Also, APPENDIX B: PHYSICAL DESCRIPTION OF THE the superconducting qubit or multi-level system can be cou- SUPERCONDUCTING QUBIT pled to bosonic modes by a resonator or cavity, simulations such as spin-boson phenomena, quantum random walks and quantum statistical models are expected. So, our results pave In the past decades, there has been a great progress in the the way for more applications of the superconducting quan- field of superconducting qubits. The main aims are to achieve tum platform. better control and longer coherent time for qubit or qubits. As a result, many types of superconducting qubits have been developed, each of which has its own advantages and lim- its. There are three main categories of the quantum supercon- ACKNOWLEDGMENTS ducting qubit working in different regimes according the ratio of the Josephson energy EJ to the charging energy EC [46– The first two authors, X.Y.G. and C.Y., contributed equally 49]: 1) charge qubit with EJ/EC ∼ 0.1; 2) flux qubit with 6 to this work. We thank Wuxin Liu and Haohua Wang EJ/EC ∼ 50; and 3) phase qubit EJ/EC ∼ 10 . By adding of Zhejiang University for technical support. We thank a large capacitor CS parallel to the superconducting quantum Ling-An Wu for careful reading our manuscript to improve interference device (SQUID) and thus shunting the later (cf. our presentation. This work was supported by National Fig. 5(a)), the transmon qubit works in the parameter regime Key Research and Development Program of China (Grant of EJ/EC being the order of several tens or several hundreds. Nos. 2016YFA0302104, 2016YFA0300600, 2014CB921401, It gains the advantage of exponentially suppressing the sen- 2017YFA0304300), National Natural Science Foundation sitivity to the charge noise at the expense of polynomial re- 6 duction of the anharmonicity [38, 46–49]. Notice that anhar- By varying the amplitude A and phase φ0 of the driving volt- monicity describes the variation of the energy level spacing age Vg which is applied through the gate capacitor Cg, we which ensures the possibility of addressing the lowest energy would have full control of the rotations of the qubit along the levels of the platform. The Hamiltonian of an isolated trans- X as well as the Y direction. The Z-direction control is exerted mon qubit is via the change control current Iext which adjusts the external 2 magnetic flux Φext thrusting through the SQUID loop. Hˆ0 = 4ECnˆ − EJ cos φ.ˆ (7) wheren ˆ is the operator corresponding to the number of Cooper pair tunneled through the Josephson junctions and φˆ APPENDIX C: THE MANY-BODY SYSTEMS AND THE denotes the gauge-invariant phase difference operator across TWO-BAND MODEL the Josephson junctions. They are mutually conjugate and satisfy the commutation relation [φ,ˆ nˆ] = i. The charg- The Hamiltonian of a two-band model is written as, ing energy E = e2/2C depends on the total capacitance C Σ X † CΣ = CS + Cg + CJ of the shunt capacitor CS , gate capacitor H = Ψkh(k)Ψk (14) k Cg and the Josephson junction capacitance CJ. The Joseph- son energy EJ = ~IC/2e is determined by the critical current where Ψk denotes a spinor, h(k) takes the form of the DC-SQUID, which is modulated by the external magnetic flux. In the transmon regime, CS is very large such that h(k) = d0(k) + d(k) · σ, (15) 20.EJ/EC.100, its `-th eigenenergy level should be [38, 46– 49] where σ = (σx, σy, σz) is a vector of Pauli matrices, as al- ready presented in the main text. This model can describe a p 1 E` ' −EJ + 8EJ EC(` + 1/2) − EC[`(` + 1) + 1/2]. (8) variety of physically different many-body systems. For exam- 2 ples, the Su-Schrieffer-Heeger (SSH) model [42] describes the The anharmonicity is big enough for the addressability of the simplest one-dimensional topological insulator. We have that two lowest energy levels and thus constitutes a qubit † † † Ψk = (ckA, ckB) with A and B referring two sub-lattices, d(k) = 1 [(t + δt) + (t − δt) cos k, (t − δt) sin k, 0)], with (t ± δt) being the Hˆ = ~ωσˆ (9) hopping amplitudes in the unit cell and between the adjacent 0 2 z √ cells, respectively. Another example is p-wave Kitaev chain with ω = 2( 8EJ EC −EC)/~ andσ ˆ z = |ei he|−|gi hg|. Here |gi [43] which describes a one-dimensional topological supercon- is the ground state of the transmon while |ei is its first excited † † ductor. For this case, we have that Ψk = (ck, c−k) and d(k) = eigenstate. [0, ∆ sin k, (−µ/2−t cos k)], where ∆ denotes the pairing poten- The quantum platform we employed is a superconducting tial, µ the chemistry potential, t the hopping amplitude. When Xmon qubit, which is designed on the basis of coplanar trans- t = ∆, it corresponds to the transverse field Ising model [44] mon. Essentially it is equivalent to a grounded transmon, see after mapping to the free fermions by Jordan-Wigner transfor- Fig. 5. Embedded in an uninterrupted ground plane, the Xmon mation. Here, we have d(k) = [0, sin k, (g−cos k)]. The details qubit can prolong the coherent time by further employing are as follows. coplanar waveguide made with high-quality material. Better The SSH model is the simplest two-band model describing connectivity can be accomplished via a cross-shaped capaci- polyacetylene, which is a one-dimensional topological insulator [37, 46–49]. The driving microwave applied as shown in tor. The Hamiltonian reads Fig. 5(b) is the following X † † HSSH = (t + δt)cA,lcB,l + (t − δt)cA,l+1cB,l + h.c. Vg = V cos(ωt + φ0). (10) l The frequency of the microwave is chosen to match the reso- (16) nant frequency of the isolated qubit in Eq.(9). As a result, the A and B refer to two sublattices. The hopping amplitude in the Hamiltonian of the Xmon is unit cell is t + δt while that between adjacent unit cell is t − δt. VgCgneˆ 1 P −ikl 2 Performing the Fourier transformation ckA = √ e cA,l Hˆ = 4ECnˆ + − EJ cos φ.ˆ (11) N l C P − Σ and c = √1 e iklc , where N is the number of sites, we kB N l B,l The two-level qubit Hamiltonian via truncating all the higher obtain energy levels is thus H = P (t + δt)(c† c + h.c.) AC e SSH k kA kB ˆ 1 g H = ~ (ω + ωΦ) σˆ z + cos(ωt + φ0)σ ˆ x, (12) − ik † 2 CΣ + (t δt)(e ckAckB + h.c.). (17) where ωΦ is the energy level shift caused by the external mag- † † † Introducing the spinor Ψk = (ckA, ckB), the Hamiltonian can netic flux Φext controlled by varying Iext. Moving to the inter- be written in a compact form, action picture with respect to Hˆ0, we would have P † − AC e HSSH = k Ψk[((t + δt) + (t δt) cos k)σx ˆ g Hr = ~ωΦσˆ z + cos φ0σˆ x + sin φ0σˆ y . (13) − CΣ + (t δt) sin kσy]Ψk, (18) 7

Cg Cg IC

CS CJ Vg CS CJ Φext Iext

(a)Isolated transmon circuit (b)Xmon circuit with control

FIG. 5. Equivalent circuits of (a) the isolated transmon circuit and (b) its Xmon variety with control. where d(k) = [(t+δt)+(t−δt) cos k, (t−δt) sin k, 0)] referring to For Hamiltonian (14), see also (1) in the main text, one can Eq.(15). The system is topologically nontrivial when |t + δt| < find that each k mode is decoupled, so we can investigate each |t − δt|. Otherwise it is topologically trivial. mode separately. The eigenvalues of h(k) are given by The p-wave Kitaev chain is a one-dimensional topological ± | | superconductor introduced by Kitaev [43]. The Hamiltonian ±(k) = d0(k) d(k) . (25) reads, The corresponding eigenvectors are denoted by |φ±(k)i, or X † † † written as density matrices HK = − t(cl cl+1 + h.c.) − ∆(cl cl+1 + h.c.) l 1 h i ρ±(k) = |φ±(k)ihφ±(k)| = 1 ± dˆ(k) · σ , (26) † 1 2 − µ(cl cl − ) (19) 2 ˆ d(k) where d(k) = |d(k)| corresponding to a unique vector on the t is the hopping amplitude; ∆ is the p-wave superconductor Bloch sphere. pairing potential and µ is the chemical potential. Performing 1 P −ikl To study the quench dynamics, we first prepare the system the Fourier transformation ck = √ l e cl and introducing N in ground state of the initial Hamiltonian hi(k), i.e. ρi(k) = † † 1 h i the spinor Ψ = (c , c− ), we obtain, ˆ k k k |φi(k)ihφi(k)| = 2 1 − di(k) · σ , corresponding to the minus X † eigenvector in (26). Then taking a sudden quench to the final HK = Ψ [∆ sin kσy + (−µ/2 − t cos k)σz]Ψk, (20) k Hamiltonian h f (k), which determines d f (k). The state evolves k −ith (k) as |φ(k, t)i = e f |φi(k)i. A more enlightening picture can where d(k) = [0, ∆ sin k, (−µ/2 − t cos k)] referring to Eq.(15). be presented as density matrix form, This model is mathematically equivalent to the transverse field 1 h i Ising model when t = ∆, and it is topologically nontrivial ρ(k, t) = |φ(k, t)ihφ(k, t)| = 1 − dˆ(k, t) · σ , (27) 2|t| 2 when |µ| > 1. The transverse field Ising model is described as, where X −itd f (k)·σ itd f (k)·σ x x z dˆ (k, t) · σ = e dˆ i (k) · σ e . (28) HIsing = − σl σl+1 + gσl . (21) l It is simply the spin precession on the Bloch sphere, that is, π g is the transverse field strength. The spin model can be dˆ i(k) rotates around dˆ f (k) with period . |d f (k)| mapped to the free-fermion model by using Jordan-Wigner transformation z † x Y † † APPENDIX D: EXPERIMENTAL SCHEME σl = 1 − 2cl cl, σl = (1 − 2c j c j)(cl + cl ). (22) j 1. tial state |φii in Eq.(29) is an eigenvector of σ , σ = 8

FIG. 6. Experiment control sequence. The initial state is prepared at the state initialization period by control quantity A0 cos(ωt + ϕ0). The state is prepared as |φii. Then for a quantum quench, by controlling ϕk depending on momentum k, we adjust the direction of the rotation axis, shown in Eq.(13).

† U hi(k)U, which corresponds to hi(k) by a unitary transforma- the existence of zeroes of G(t) requires that there are at least tion U. Then for a quantum quench, the applied Hamiltonian one critical momentum k∗ satisfying † takes the form U h f (k)U. ˆ ∗ ˆ ∗ Now, let us show how to realize a qubit rotation in experi- di(k ) · d f (k ) = 0, (32) ment. In the Hamiltonian of a qubit Eq.(13), the two terms in ˆ ˆ the parentheses represent a rotation with axis in the XY plane. i.e., the vector di(k) is perpendicular to d f (k) at the critical k∗ The direction can be adjusted by controlling the parameter φ0 momentum , and the DQPT occurs at in Eq.(13), which depends on momentum k. Experimentally π 1 as shown in Fig. 6, by controlling ϕk in Ak cos(ωt + ϕk), we tc = (n + ), n = 0, 1, 2,..., (33) can realize the control of axis direction in the XY plane for a |d f (k)| 2 rotation. Explicitly in Fig.2 in the main text, we can find that ˆ ∗ −ˆ ∗ the rotation axes are in the XY plane. and the Bloch vector satisfies d(k , tc) = di(k ) [24]. For clarity, we employ the Ising model to elucidate the condition Our simulation scheme can be applied to general two-band of DQPT. One has, models. For example as shown in Ref.[45], the temporal topological phenomena can be simulated by a qubit subjected to a 2 2 di(k) · d f (k) = sin k + cos k − (gi + g f ) cos k + gig f two-frequency drive. The Hamiltonian takes the form, = 0. (34) x y Htemp. = v1 cos(ω1t + ϕ1)σ + v2 sin(ω2t + ϕ2)σ The solution exists when z +[m − b1 cos(ω1t + ϕ1) − b2 cos(ω2t + ϕ2)]σ , 1 + gig f (30) | cos k| = < 1. (35) gi + g f where the notations and the implication of this model can be One can obtain sgn[(1 − |gi|)(1 − |g f |)] = −1, i.e. DQPT oc- found in Ref.[45]. This Hamiltonian corresponds to Eq.(13), curs if and only if the initial Hamiltonian and the final Hamil- and can be realized by a superconducting qubit. The rotation tonian belong to different phases for the Ising model, and axis should be in arbitrary direction. The topological phenom- 1+g g k∗ = ± arccos i f . ena are described by whether the whole Bloch sphere of the gi+g f corresponding states are covered or not, which is similar with We also know from Refs. [24] and [26] that a dynamical our experiment performed. Chern number can be defined in momentum-time space in a quench process. First we should find the fixed points km that satisfying dˆ i(km) is parallel and anti-parallel to dˆ f (km). Here we just focus on the transverse field Ising model, there are APPENDIX E: RELATION BETWEEN DYNAMICAL QUANTUM PHASE AND DYNAMICAL CHERN NUMBER only two fixed points k = 0 and k = π. Then the dynamical Chern number is defined as Z π Z π In Ref. [50], it is shown that the Loschmidt amplitude of a 1 0 Cdyn = dk dt (dˆ × ∂t0 dˆ) · ∂kdˆ (36) two-band system can be written as 4π 0 0 Y ˆ ˆ where t0 = t is the rescaled time. G(t) = [cos(|d f (k)|t) + idi(k) · d f (k) sin(|d f (k)|t)], (31) d f k For the fixed point k = 0, we have dˆ(0) = (0, sgn(g − 1), 0), and for the fixed point k = π, we have dˆ(π) = (0, sgn(g+1), 0). where d (k) and d (k) correspond to the pre- and post-quench i f The dynamical Chern number is calculated, Hamiltonians, respectively. The dynamical quantum phase transition (DQPT) occurs when the Loschmidt amplitude 1 C = (cos θ − cos θ ), (37) reaches zero at a critical time tc. As we see from Eq. (31), dyn 2 k=0 k=π 9 where θk is the induced angle between dˆ i(k) and dˆ f (k). In our APPENDIX F: ERROR BAR SHOWN IN THE FIGURES experiment, we first choose gi = 0.2 and g f = 1.5, hence dˆ i(0) · dˆ f (0) = −1 and dˆ i(π) · dˆ f (π) = 1, the dynamical Chern The dynamical free energy can be expressed in terms of number is Cdyn = −1. As a result the Bloch sphere is fully dˆ i(k) and dˆ(k, t) covered as shown in Fig. 2(g) in the main text. From the continuity of the function dˆ i(k) · dˆ f (k), there must be a critical 1 X 1 + dˆ i(k) · dˆ(k, t) ∗ ∗ ∗ f (t) = − log , (39) momentum k between 0 and π satisfying dˆ i(k )·dˆ f (k ) = 0, so N 2 k we can draw a conclusion that the nontrivial dynamical Chern number ensures the occurrence of DQPT. In our experimental setup, dˆ (k) is a fixed unit vector. To es- ˆ i We also choose gi = 0.2 and g f = 0.5, we have di(0) · timate the experimental error of f (t), we need only to esti- ˆ ˆ ˆ d f (0) = di(π) · d f (π) = 1, and hence the dynamical Chern mate the fluctuation of dˆ i(k). Given a specific k, we have ob- number Cdyn = 0. In this case the Bloch sphere is not fully tained 70 state tomography data dˆ(k, t) corresponding to dif- covered as shown in Fig. 2(n) in the main text, and the DQPT ferent time points on the evolution path on the Bloch sphere. would not occur. Each of these tomography data dˆ(k, t) is an average of 5000 The nontrivial dynamical Chern number indicates the emer- raw data. We estimate the fluctuation of dˆ(k, t) by estimating gence of Skyrmion lattice in the momentum-time space. If the fluctuation of the radius of the evolution path traced on the gi = 0.2 and g f = 1.5, the dynamical Chern number is non- Bloch sphere. For each path, we choose three equally sepa- trivial, we consider the expectation value rated state points dˆ(k, t) and calculate the radius determined. Thus for k, we obtain 22 estimation of the evolution path. The ˆ ˆ hdˆ(k, t)i ≡ hφ(k, t)|dˆ · σ|φ(k, t)i = −dˆ(k, t) · dˆ i(k). (38) magnitude ∆d(k) of the fluctuation of d(k, t) is evaluated by the standard deviation of the 22 estimation of the radius. The ∗ At k = k and t = tc, hdˆ(k, t)i reaches the minimum −1 and error of the dynamical free energy is hence ∗ (k, t) = (k , tc) is the center in the texture of pseudospin as shown in Fig. 4(c). It forms a lattice during the time evolution 1 X ∆dˆ(k) ∆ f (t) = . (40) with the lattice spacing is just the period of DQPT π . In N ˆ ˆ |d f (k)| k 1 + di(k) · d(k, t) the case gi = 0.2 and g f = 0.5, the DQPT would not occur, the dynamical Chern number is trivial and Skyrmion lattices Those error bars are indicated in the Fig. 3 and Fig. 4 in the would not appear as shown in Fig. 4(d) in the main text. main text.

[1] J. Zhang, P. W. Hess, A. Kyprianidis, P. Becker, A. Lee, J. “Observation of many-body localization of interacting fermions Smith, G. Pagano, I. D. Potirniche, A. C. Potter, A. Vishwanath, in a quasirandom optical lattice,” Science 349, 842-845 (2015). N. Y. Yao, and C. Monroe, “Observation of a discrete time crys- [8] G. A. Alvarez, D. Suter, and R. M. Kaiser, “Localization- tal,” Nature 543, 217-220 (2017). delocalization transition in the dynamics of dipolar-coupled nu- [2] S. Choi, J. Choi, R. Landig, G. Kucsko, H. Zhou, J. Isoya, F. clear spins,” Science 349, 846-848 (2015). Jelezko, S. Onoda, H. Sumiya, V. Khemani, C. Keyserlingk, N. [9] J. Y. Choi, S. Hild, J. Zeiher, P. Schauß, A. Rubio-Abadal, Y. Yao, E. Demler, and M. D. Lukin, “Observation of discrete T. Yefsah, V. Khemani, D. A. Huse, I. Bloch, and C. Gross, time-crystalline order in a disordered dipolar many-body sys- “Exploring the many-body localization transition in two dimen- tem,” Nature 543, 221-225 (2017). sions,” Science 352, 1547-1551 (2016). [3] J. Zhang, G. Pagano, P. W. Hess, A. Kyprianidis, P. Becker, H. [10] J. Smith, P. Richerme, B. Neyenhuis, P. W. Hess, P. Hauke, M. Kaplan, A. V. Gorshkov, Z.-X. Gong, and C. Monroe, “Obser- Heyl, M. A. Huse, and C. Monroe, “Many-body localization vation of a many-body dynamical phase transition with a 53- in a quantum simulator with programmable random disorder,” qubit quantum simulator,” Nature 551, 601-604 (2017). Nat. Physics 12, 907-911 (2016). [4] H. Bernie, S. Schwartz, A. Keesling, H. Levine, A. Omran, [11] K. Xu, J. J. Chen, Y. Zeng, Y. R. Zhang, C. Song, W. Liu, Q. H. Pichler, S. Choi, A. S. Zibrov, M. Endres, M. Greiner, V. Guo, P. Zhang, D. Xu, H. Deng, K. Huang, H. Wang, X. B. Zhu, Vuletic,´ and M. D. Lukin, “Probing many-body dynamics on a D. N. Zheng, and H. Fan, “Emulating many-body localization 51-atom quantum simulator,” Nature 551, 579-584 (2017). with a superconducting quantum processor,” Phys. Rev. Lett. [5] P. Jurcevic, H. Shen, P. Hauke, C. Maier, T. Brydges, C. 120, 050507 (2018). Hempel, B. P. Lanyon, M. Heyl, R. Blatt, and C. F. Roos, “Di- [12] P. Roushan, C. Neill, J. Tangpanitanon, V. M. Bastidas, A. rect observation of dynamical quantum phase transitions in an Megrant, R. Barends, Y. Chen, Z. Chen, B. Chiaro, A. interacting many-body system,” Phys. Rev. Lett. 119, 080501 Dunsworth, A. Fowler, B. Foxen, M. Giustina, E. Jeﬀrey, J. (2017). Kelly, E. Lucero, J. Mutus, M. Neeley, C. Quintana, D. Sank, [6] N. Flaschner,¨ D. Vogel, M. Tarnowski, B. S. Rem, D. Luhmann,¨ A. Vainsencher, J. Wenner, T. White, H. Neven, D. G. Ange- M. Heyl, J. C. Budich, L. Mathey, K. Sengstock, and C. Weit- lakis, and J. Martinis, “Spectroscopic signatures of localization enberg, “Observation of dynamical vortices after quenches in a with interacting photons in superconducting qubits,” Science system with topology,” Nature Physics 14, 265-268 (2018). 358, 1175-1179 (2017). [7] M. Schreiber, S. S. Hodgman, P. Bordia, H. P. Luschen,¨ M. [13] X. S. Tan, D. W. Zhang, Q. Liu, G. M. Xue, H. F. Yu, Y. H. Fischer, R. Vosk, E. Altman, U. Schneider, and I. Bloch, Q. Zhu, H. Yan, S. L. Zhu, and Y. Yu, “Topological Maxwell 10

metal bands in a superconducting qutrit,” Phys. Rev. Lett. 120, [33] C. Yang, Y. Wang, P. Wang, X. Gao, and S. Chen, “Dynami- 130503 (2018). cal signature of localization-delocalization transition in a one- [14] X. S. Tan, M. M. Li, D. Y. Li, K. Z. Dai, H. F. Yu, and Y. Yu, dimensional incommensurate lattice,” Phys. Rev. B 95, 184201 “Demonstration of Hopf-link semimetal bands with supercon- (2017). ducting circuits,” App. Phys. Lett. 112, 172601 (2018). [34] L. Zhou, Q. Wang, H. Wang, and J. Gong, “Dynamical quantum [15] X. Y. Guo, Y. Peng, C. N. Peng, H. Deng H, Y. R. Jin, C. phase transitions in non-Hermitian lattices,” arXiv:1711.10741 Tang, X. B. Zhu, D. Zheng, and H. Fan, “Demonstration of (2017). irreversibility and dissipation relation of thermodynamics with [35] M. Heyl, “Dynamical quantum phase transitions: a review,” a superconducting qubit,” arXiv:1710.10234 (2017). Rep. Prog. Phys. 81, 054001 (2018). [16] X. H. Peng, H. Zhou, B. B. Wei, J. Cui, J. F. Du, and R. [36] A. A. Zvyagin, “Dynamical quantum phase transitions,” Low B. Liu, “Experimental Observation of Lee-Yang Zeros,” Phys. Temp. Phys. 42, 971-994 (2016). Rev. Lett. 114, 010601 (2015). [37] R. Barends, J. Kelly, A. Megrant, D. Sank, E. Jeffrey, Y. Chen, [17] C. Y. Lu, W. B. Gao, O. Guhne,¨ X. Q. Zhou, Z. B. Chen, and J. Y. Yin, B. Chiaro, J. Mutus, C. Neill, P. OMalley, P. Roushan, W. Pan, “Demonstrating anyonic fractional statistics with a six- J. Wenner, T. C. White, A. N. Cleland, and J. M. Martinis, “Co- qubit quantum simulator,” Phys. Rev. Lett. 102, 030502 (2009). herent Josephson qubit suitable for scalable quantum integrated [18] Y. P. Zhong, D. Xu, P. Wang, C. Song, Q. J. Guo, W. X. Liu, K. circuits,” Phys. Rev. Lett. 111, 080502 (2013). Xu, B. X. Xia, C. Y. Lu, S. Han, J. W. Pan, and H. Wang, “Emu- [38] J. Koch, T. M. Yu, J. Gambetta, A. A. Houck, D. I. Schus- lating anyonic fractional statistical behavior in a superconduct- ter, J. Majer, A. Blais, M. H. Devoret, S. M. Girvin, and R. J. ing quantum circuit,” Phys. Rev. Lett. 117, 110501 (2016). Schoelkopf, “Charge-insensitive qubit design derived from the [19] E. A. Martinez, C. A. Muschik, P. Schindler, D. Nigg, A. Er- Cooper pair box,” Phys. Rev. A 76, 042319 (2007). hard, M. Heyl, P. Hauke, M. Dalmonte, T. Monz, P. Zoller, and [39] R. Barends, J. Kelly, A. Megrant, A. Veitia, D. Sank, E. Jef- R. Blatt, “Real-time dynamics of lattice gauge theories with a frey, T. C. White, J. Mutus, A. G. Fowler, B. Campbell, Y. few-qubit quantum computer,” Nature 534, 516-519 (2016). Chen, Z. Chen, B. Chiaro, A. Dunsworth, C. Neill, P. OMal- [20] M. Heyl, A. Polkovnikov, and S. Kehrein, “Dynamical quantum ley, P. Roushan, A. Vainsencher, J. Wenner, A. N. Korotkov, A. phase transitions in the transverse-field Ising model,” Phys. Rev. N. Cleland, and J. M. Martinis, “Superconducting quantum cir- Lett. 110, 135704 (2013). cuits at the surface code threshold for fault tolerance,” Nature [21] M. Heyl, “Dynamical quantum phase transitions in systems 508, 500-503 (2014). with broken-symmetry phases,” Phys. Rev. Lett. 113, 205701 [40] N. Ofek, A. Petrenko, R. Heeres, P. Reinhold, Z. Leghtas, B. (2014). Vlastakis, Y. Liu, L. Frunzio, S. M. Girvin, L. Jiang, M. Mir- [22] M. Heyl, “Scaling and universality at dynamical quantum phase rahimi, M. H. Devoret, and R. J. Schoelkopf, “Extending the transitions,” Phys. Rev. Lett. 115, 140602 (2015). lifetime of a quantum bit with error correction in superconduct- [23] J. Eisert, M. Friesdorf, and C. Gogolin, “Quantum many-body ing circuits,” Nature 536, 441-445 (2016). systems out of equilibrium,” Nat. Phys. 11, 124-130 (2015). [41] A. Blais, R.-S. Huang, A. Wallraff, S. M. Girvin, and R. J. [24] C. Yang, L. Li, and S. Chen, “Dynamical topological invariant Schoelkopf, “Cavity quantum electrodynamics for supercon- after a quantum quench,” Phys. Rev. B 97, 060304(R) (2018). ducting electrical circuits: An architecture for quantum com- [25] H. F. Lang, Y. X. Chen, Q. T. Hong, and H. Fan, “Dynami- putation,” Phys. Rev. A 69, 062320 (2004). cal quantum phase transition for mixed states in open systems,” [42] W. Y. Su, J. R. Schrieffer, and A. J. Heeger, “Solitons in poly- Phys. Rev. B 98, 134310 (2018). acetylene,” Phys. Rev. Lett. 42, 1698 (1979). [26] Z. Gong and M. Ueda, “Topological entanglement-spectrum [43] A. Yu. Kitaev, “Unpaired Majorana fermions in quantum crossing in quench dynamics,” arXiv:1710.05289 (2017). wires,” Phys. Usp. 44, 131 (2001). [27] E. Canovi, P. Werner, and M. Eckstein, “First-order dynamical [44] E. Barouch and B. M. Mccoy, “Statistical mechanics of XY- phase transitions,” Phys. Rev. Lett. 113, 265702 (2014). model. II. Spin-correlation functions,” Phys. Rev. A 3, 786 [28] J. C. Budich and M. Heyl, “Dynamical topological order param- (1971). eters far from equilibrium,” Phys. Rev. B 93, 085416 (2016). [45] I. Martin, G. Refael, and B. Halperin, “Topological frequency [29] S. Sharma, U. Divakaran, A. Polkovnikov, and A. Dutta, “Slow conversion in strongly driven quantum systems,” Phys. Rev. X quenches in a quantum Ising chain: Dynamical phase transi- 7, 041008 (2017). tions and topology,” Phys. Rev. B 93, 144306 (2016). [46] J. M. Fink, “Quantum nonlinearities in strong coupling circuit [30] Z. Huang and A. V. Balatsky, “Dynamical quantum phase tran- QED,” Ph.D thesis, ETH Zurich (2010). sitions: Role of topological nodes in wave function overlaps,” [47] S. M. Girvin, “Circuit QED: Superconducting Qubits Coupled Phys. Rev. Lett. 117, 086802 (2016). to Microwave Photons,” Lecture Notes of the Les Houches [31] T. Nemoto, R. L. Jack, and V. Lecomte, “Finite-size scaling of Summer School 96, Chap. 3, 113-255 (Oxford University Press, a first-order dynamical phase transition: Adaptive population 2011). dynamics and an effective model,” Phys. Rev. Lett. 118, 115702 [48] G. Wendin, “Quantum information processing with supercon- (2017). ducting circuits: a review,” Rep. Prog. Phys. 80, 106001 (2017). [32] M. Heyl and J. C. Budich, “Dynamical topological quantum [49] X. Gu, A. F. Kockum, A. Miranowicz, Y. X. Liu, and F. phase transitions for mixed states,” Phys. Rev. B 96, 180304(R) Nori, “Microwave photonics with superconducting quantum (2017). circuits,” Phys. Rep. 718-719, 1-102 (2017). [50] S. Vajna and B. Dora,´ “Topological classification of dynamical phase transitions,” Phys. Rev. B 91, 155127 (2015).