We study the all-to-all Ising model in the presence of dissipation and periodic external field. The corresponding Lindblad equation has a time-periodic Liouvillian. The dynamics of magnetizations is explored by using both the mean-field theory and numerical simulation at a finite number of spins. The mean-field result exhibits a transition from the periodic response at small field amplitude to the chaotic dynamics at large amplitude. The Lyapunov exponents are calculated, which supports the existence of a classical chaotic phase. However, in the numerical simulation at a finite number of spins, the response is periodic at both small and large amplitudes. The scaling analysis of the Floquet Liouvillian spectrum suggests that the periodic response persists even in the thermodynamic limit. Further analysis shows that, in the classical chaotic regime, the density-matrix elements have a wide dispersion, instead of being localized. Such a delocalization of wave packets explains the failure of mean-field approximation.
We study the all-to-all Ising model in the presence of dissipation and periodic external field. The corresponding Lindblad equation has a time-periodic Liouvillian. The dynamics of magnetizations is explored by using both the mean-field theory and numerical simulation at a finite number of spins. The mean-field result exhibits a transition from the periodic response at small field amplitude to the chaotic dynamics at large amplitude. The Lyapunov exponents are calculated, which supports the existence of a classical chaotic phase. However, in the numerical simulation at a finite number of spins, the response is periodic at both small and large amplitudes. The scaling analysis of the Floquet Liouvillian spectrum suggests that the periodic response persists even in the thermodynamic limit. Further analysis shows that, in the classical chaotic regime, the density-matrix elements have a wide dispersion, instead of being localized. Such a delocalization of wave packets explains the failure of mean-field approximation.
We study a periodically driven macrospin system with anisotropic long-range interactions and collective dissipation, described by a Lindblad master equation. In the thermodynamic limit (N→∞), a mean-field treatment yields classical equations of motion, whose dynamics are characterized via the maximal Lyapunov exponent (MLE). Focusing on the thermodynamic limit, we map out chaotic, quasiperiodic, and periodic phases via bifurcation diagrams, MLEs, and Fourier spectra of evolved observables, identifying classic period-doubling bifurcations and fractal boundaries in the regions of attractors. Finite-size quantum simulations in the Dicke basis reveal that while both quantum and classical systems exhibit diverse dynamical phases, finite-size effects suppress some behaviors present in the thermodynamic limit. The sign of λmax serves as a key indicator of convergence between quantum and classical dynamics, which agree over timescales up to the Lyapunov time. Analysis of the density matrix shows that convergence occurs only when its nonzero elements are sharply localized. However, the nonconvergence does not imply a fundamental difference between quantum and classical dynamics: in chaotic regimes, although the evolution orbits of quantum and classical systems show significant differences, quantum evolution becomes mixed and diffusively explores the Hilbert space, signaling quantum chaos, which can be confirmed by the delocalized nature of the density matrix.
Abstract The Lindblad equation describes the dissipative time evolution of a density matrix that characterizes an open quantum system in contact with its environment. The widespread ensemble interpretation of a density matrix requires its time evolution to be linear. However, when the dynamics of the density matrix is of a quantum system results not only from the interaction with an external environment, but also with other quantum systems of the same type, the ensemble interpretation is inappropriate and nonlinear dynamics arise naturally. We therefore study the dynamical behavior of nonlinear Lindblad equations using the example of a two-level system. By using techniques developed for classical dynamical systems we show that various types of bifurcations and even chaotic dynamics can occur. As specific examples that display the various types of dynamical behavior, we suggest explicit models based on systems of interacting spins at finite temperature and exposed to a magnetic field that can change in dependence of the magnetization. Due to the interaction between spins, which is treated at mean-field level, the Hamiltonian as well as the transition rates of the Lindblad equation become dependent on the density matrix.
We study the critical dynamics arising from a time-dependent periodic homogenous source coupled to the order-parameter field, which drives a classical ferromagnetic system across a continuous transition. For this purpose, we consider the paradigmatic two-dimensional (2D) Ising model in the presence of a periodic magnetic field h(t)=−Acos(2πt/P), evolving under a purely relaxational dynamics at the critical temperature. We show that the periodic driving gives rise to a peculiar dynamic scaling behavior in the thermodynamic limit, arising from a nontrivial interplay among the time t, the amplitude A and period P of h(t). The relevant scaling variables are τ=t/P and σ=APκ, with κ=yh/z, where yh=(d+2−η)/2 is the critical dimension of the magnetic field, and z is dynamic exponent for the critical relaxational dynamics (κ≈0.865 for the 2D Ising model). The dynamic scaling behaviors of the magnetization and bond-energy density show an oscillatory behavior around a smooth curve which approaches a large-τ stationary behavior. We also briefly discuss the dynamic behavior of an Ising system driven across the critical point by a periodic time-varying temperature at zero magnetic field.
We investigate the dynamical phases and phase transitions arising in a classical two-dimensional anisotropic XY model under the influence of a periodically driven temporal external magnetic field in the form of a symmetric square wave. We use a combination of finite temperature classical Monte Carlo simulation, implemented within a CPU+GPU paradigm, utilizing local dynamics provided by the Glauber algorithm and a phenomenological equation-of-motion approach based on relaxational dynamics governed by the time-dependent free energy within a mean-field approximation to study the model. We investigate several parameter regimes of the variables (magnetic field, anisotropy, and the external drive frequency) that influence the anisotropic XY system. We identify four possible dynamical phases: Ising-SBO, Ising-SRO, XY-SBO, and XY-SRO. Both techniques indicate that only three of them (Ising-SRO, Ising-SBO, and XY-SRO) are stable dynamical phases in the thermodynamic sense. Within the Monte Carlo framework, a finite-size scaling analysis, shows that XY-SBO does not survive in the thermodynamic limit giving way to either an Ising-SBO or a XY-SRO regime. The finite-size scaling analysis further shows that the transitions between the three remaining dynamical phases either belong to the two-dimensional Ising universality class or are first-order in nature. Within the mean-field calculations yield three stable dynamical phases, i.e., Ising-SRO, Ising-SBO and XY-SRO, where the final steady state is independent of the initial condition chosen to evolve the equations of motion, as well as a region of bistability where the system flows to either Ising-SBO or XY-SRO (Ising-SRO) depending on the initial condition. Unlike the stable dynamical phases, the XY-SBO represents a transient feature that is eventually lost to either Ising-SBO or XY-SRO. Our mean-field analysis highlights the importance of the competition between switching of the stationary point(s) of the free energy after each half cycle of the external field and the two-dimensional nature of the phase space for the equations of motion.
In this thesis we study the dynamics of the long-range interacting Ising chain using a cluster mean-field approach, a mean-field theory that accounts for short range correlations. The principal result is that, if the system is isolated from an external environment, the model exhibits two different behaviors depending on the range of the interactions. Whenever the system is truly long-range, it exhibits a sharp mean-field dynamical phase transition from a dynamical ferromagnetic to a dynamical paramagnetic phase. Reducing the range of the interactions a critical region, showing hypersensitivity to initial conditions, appears. Interestingly, this chaotic region shares the same physics of a classical tossed coin that is allowed to bounce on the floor. In a second part of the work we derive the cluster mean-field equations of motion describing the dynamics of a fully connected Ising chain connected to an external bath. In particular we focused on the dynamics in presence of a dissipation generated by string of Glauber operators acting on one or more sites of the chain. This model, in presence of global dissipative processes, exhibits persistent oscillations in time revealing the existence of a boundary time-crystal and we studied the stability of the boundary time-crystal showing that global dissipative processes are a key ingredient for their existence.
We study correlations, transport, and chaos in a Heisenberg magnet as a classical model many-body system. By varying temperature and dimensionality, we can tune between settings with and without symmetry breaking and accompanying collective modes or quasiparticles (spin waves) which in the limit of low temperatures become increasingly long-lived. Changing the sign of the exchange interaction from a ferro- to an antiferromagnetic one varies the spin-wave spectrum, and hence the low-energy spectral density. We analyze both conventional and out-of-time-ordered spin correlators (decorrelators) to track the spreading of a spatiotemporally localized perturbation---the wing beat of the butterfly---as well as transport coefficients and Lyapunov exponents. We identify a number of qualitatively different regimes. Trivially, at T=0, there is no dynamics at all. In the limit of low temperature, T=0+, integrability emerges, with infinitely long-lived magnons; here the wave packet created by the perturbation propagates ballistically, yielding a light cone at the spin-wave velocity which thus subsumes the butterfly velocity; inside the light cone, a pattern characteristic of the free spin-wave spectrum is visible at short times. On top of this, residual interactions (nonlinearities in the equations of motion) lead to spin-wave lifetimes which, while divergent in this limit, remain finite at any nonzero T. At the longest times, this leads to a standard chaotic regime; for this regime, we show that the Lyapunov exponent is simply proportional to the (inverse) spin-wave lifetime. Visibly strikingly, between this and the short-time integrable regimes, a scarred regime emerges: Here, the decorrelator is spatiotemporally highly nonuniform, being dominated by rare and random scattering events seeding secondary light cones. As the spin correlation length decreases with increasing T, the distinction between these regimes disappears and at high temperature the previously studied chaotic paramagnetic regime emerges. For this, we elucidate how, somewhat counterintuitively, the ballistic butterfly velocity arises from a diffusive spin dynamics.