Emergence and Control of Complex Behaviours in Driven Systems of Interacting Qubits with Dissipation A. V. Andreev1, A. G. Balanov2,+, T. M. Fromhold3, M. T. Greenaway2, A. E. Hramov1, W. Li3, V. V. Makarov1, A. M. Zagoskin2,4,∗ 1Innopolis University, Universitetskaya Str. 1, Innopolis, 420500, Russia 2Department of Physics, Loughborough University, Loughborough LE11 3TU, United Kingdom 3School of Physics and Astronomy, University of Nottingham, Nottingham NG7 2RD, United Kingdom 4Department of Theoretical Physics and Quantum Technologies, National University of Science and Technology MISIS, 4 Leninsky Ave., Moscow 119049, Russia ∗ Corresponding author: [email protected] and +Corresponding author: [email protected] (Dated: November 13, 2020) Progress in the creation of large scale, artificial quantum coherent structures demands the inves- tigation of their nonequilibrium dynamics when strong interactions, even between remote parts, are non-perturbative. Analysis of multiparticle quantum correlations in a large system in the presence of decoherence and external driving is especially topical. Still, the scaling behaviour of dynamics and related emergent phenomena are not yet well understood. We investigate how the dynamics of a driven system of several quantum elements (e.g., qubits or Rydberg atoms) changes with increasing number of elements. Surprisingly, a two-element system exhibits chaotic behaviours. For larger system sizes a highly stochastic, far from equilibrium, hyperchaotic regime emerges. Its complexity systematically scales with the size of the system, proportionally to the number of elements. Finally, we demonstrate that these chaotic dynamics can be efficiently controlled by a periodic driving field. The insights provided by our results indicate the possibility of a reduced description for the behaviour of a large quantum system in terms of the transitions between its qualitatively different dynamical regimes. These transitions are controlled by a relatively small number of parameters, which may prove useful in the design, characterization and control of large artificial quantum structures. INTRODUCTION bly due to the absence of a convenient theoretical frame- work for studying these novel systems. One important Recent progress in experimental techniques for the fab- open question to be tackled is how the dynamics scales rication, control and measurement of quantum coherent with the increase of the system size (e.g., number of systems has allowed the routine creation of moderate-to- qubits), and what role in its change is played by emergent large arrays of controllable quantum units (e.g., super- phenomena. conducting qubits, trapped atoms and ions) expected to have revolutionary applications, e.g., in sensing, quantum communication and quantum computing [1–4]. They also allow the investigation of fundamental physics, such as the measurement problem [5–7], the link between the We address these questions by studying theoretically chaotic behaviour in classical systems and the corre- the far-from-equilibrium dynamics of a driven chain of sponding dynamics in their quantum counterparts [8–12] interacting two-level quantum systems with dissipation and the relative roles of entanglement and thermalization and dephasing [Fig. 1a]. We find that even a short chain arXiv:1907.03602v3 [quant-ph] 12 Nov 2020 [13–17]. comprising between two and four elements can demon- The deployment of modern quantum technologies gen- strate chaotic behaviour. Remarkably, in systems with erally requires both the improvement of unit quantum five or more elements a phenomenon known as hyper- elements (e.g., increasing their decoherence times) and chaos emerges, a dynamical regime characterized by two the scaling up of quantum structures [18]. or more positive Lyapunov exponents. Hyperchaos has It is known that large systems comprising many func- features that resemble thermalisation even though the tional elements tend to demonstrate emergent phenom- system is far from equilibrium. We show that, in con- ena in their dynamics, which can not be observed in trast to thermalisation, the deterministic origin of hyper- smaller systems [19, 20]. Examples include dynamical chaos allows it to be controlled either by tuning of the phases in driven Bose-Einstein condensates [21], synchro- system parameters, or by applying an external driving nization in micromechanical oscillator arrays [22] and co- force. Our results give new insights into the dynamics of herent THz emission from layered high-temperature su- large artificial quantum coherent structures, which will perconductors [23–25]. Nevertheless, the emergent phe- be important for the design and control of quantum sys- nomena are commonly neglected when designing and in- tems such as Rydberg molecules, quantum information vestigating large-scale coherent structures. This is possi- processors, simulators and detectors. 2 RESULTS partially quantum coherent system, for which it has been shown accurately to describe experimental measurements Mathematical model on qubits, up to the corrections due to two-point correla- tions [38, 42–45]. In such a case, the quantum coherence Here, we investigate the dynamics of a 1D chain of between spatially separated elements is disrupted, e.g. N qubits driven by an external electromagnetic field in by local ambient noise [46]. the presence of decoherence. The system under con- Rewriting the factorized density matrix ρ in terms of sideration is equivalent to a system of Rydberg atoms, the population inversion of the jth qubit, wj = (ρj )11 − ∗ by rescaling and reidentification of parameters [26–31]. (ρj )00, and its coherence, qj = (ρj )10 = (ρj )01, yields: Along with superconducting structures, Rydberg systems w˙ j = −2Ωℑqj − (wj + 1); provide a controllable testbed for studying fundamen- tal quantum dynamics [32–34] and for developing new 1 Ω q˙j = i ∆ − c (wk + 1) qj − qj + i wj . (4) quantum chemistry [35–37]. In particular, arrays of Ry- 2 2 Xk6=j dberg atoms can be controlled experimentally with very high precision [38–40]. Each atom has a ground state In (4) the dimensionless time τ = γt (i.e.x ˙ ≡ dx/dτ = −1 |gi and an electronically excited high-energy Rydberg γ dx/dt), Ω=ΩR/γ,∆= δω/γ, and c = V a/γ, where state |ei [Fig. 1b]. Rydberg states have large polariz- a is the lattice constant. In more detail the derivation of ability ∼ n7 (where n is the principal quantum number), (4) is discussed in Section A in Supplementary Materials. which generates strong and long-range interactions be- tween atoms [41]. To model an electromagnetically driven array of qubits Emergence of chaos and hyperchaos (Rydberg atoms), we use the Hamiltonian Here we study the dynamical regimes of chains of inter- N ΩR acting qubits. To characterise the complexity (random- H = −δω|eihe|j + (|eihg|j + |gihe|j ) (1) 2 ness) of the emergent chaos and hyperchaos we calculate Xj=1 the full spectrum of the Lyapunov exponents [47]. + Vi,j |eihe|j ⊗ |eihe|k, First we analyze two coupled qubits (N = 2). Pre- Xj<k viously, it has been shown that interplay between the energy pumping and dissipation can eventually trigger where δω = ωl −ω0 is the detuning between the laser and self-sustained state population oscillations in this system transition frequencies, ΩR is the Rabi frequency (tuned and even lead to the emergence of bistability, when homo- by the laser field amplitude), and Vi,j characterizes the geneous and antiferromagnetic states coexist [27]. Our interaction between ith and jth qubits. investigation reveals another interesting phenomenon as- The dynamics are described by the Liouville-von Neu- sociated with deterministic chaos, which emerges via a mann equation for the density matrix ρ, cascade of period-doubling bifurcations for a periodic os- cillations [47]. In Fig. 1c, we show a color map illustrat- ρ˙ = −i [H,ρ]+ L [ρ] , (2) ing the dependence of the largest Lyapunov exponent Λ1 where relaxation and dephasing processes are taken into on ∆ and Ω. A transition from yellow to black reflects the account through the appropriate Lindblad operator change from smaller to larger values of Λ1. The diagram clearly demonstrates the presence of chaotic dynamics, 1 for which Λ1 > 0, in considerable areas of the param- L[ρ]= γ |gihe|j ρ|eihg|j − {|eihe|j,ρ} , (3) 2 eter plane shaded by black color. This shows that the Xj discovered chaos is a robust phenomenon existing in our with γ to be the decay rate from the excited state to the system over wide parameter ranges. Notably, for certain ground state. Rydberg atom interactions can range over parameter values it coexists with the antiferromagnetic a few tens of micrometers, larger than the physical size of steady state. The detailed descriptions of the bifurca- the system. (The same is true for superconducting qubits tion transitions and analysis of the Lyapunov exponents interacting through a resonator field mode). Therefore, that characterize the stability of long-term dynamics are to simplify the analysis, we can initially assume that in- presented in Section B in Supplementary Materials. teractions are identical for any pair of qubits, Vi,j = V For N = 3 and 4, the same type of chaotic dynamics [see Fig. 1a,b]. For a system of N qubits, ρ has the di- exists in the system. However, surprisingly, for N ≥ 5 mensionality of its Hilbert space, 2N . a new type of chaotic dynamics appears, whose stability This problem can be greatly simplified by substitut- is characterized by more than one positive Lyapunov ex- ing a fully factorized density matrix approximation, ρ ≈ ponent. This type of chaos is known as hyperchaos [48]. j ⊗ρj , in (2). Here ρj is the (one-particle) density ma- A map of different dynamical regimes in the (∆, Ω) pa- Qtrix of the jth qubit. This approximation is justified for a rameter plane is shown in Fig. 2a for N = 5 and c = 3 5. This diagram was built by calculating the spectrum attractor. Accumulation of the corresponding bifurca- of Lyapunov exponents for the various limit sets that tions causes the second Lyapunov exponent to become exist in the (∆, Ω)-plane.