In a recent paper [Phys. Rev. Lett. 129, 120601 (2022)], we have shown that the dynamics of interfaces, in the symmetry-broken phase of the two-dimensional ferromagnetic quantum Ising model, displays a robust form of ergodicity breaking. In this paper, we elaborate more on the issue. First, we discuss two classes of initial states on the square lattice, the dynamics of which is driven by complementary terms in the effective Hamiltonian and may be solved exactly: (a) Strips of consecutive neighboring spins aligned in the opposite direction of the surrounding spins and (b) a large class of initial states, characterized by the presence of a well-defined ``smooth'' interface separating two infinitely extended regions with oppositely aligned spins. The evolution of the latter states can be mapped onto that of an effective one-dimensional fermionic chain, which is integrable in the infinite-coupling limit. In this case, deep connections with noteworthy results in mathematics emerge, as well as with similar problems in classical statistical physics. We present a detailed analysis of the evolution of these interfaces both on the lattice and in a suitable continuum limit, including the interface fluctuations and the dynamics of entanglement entropy. Second, we provide analytical and numerical evidence supporting the conclusion that the observed nonergodicity---arising from Stark localization of the effective fermionic excitations---persists away from the infinite-Ising-coupling limit, and we highlight the presence of a timescale $T\ensuremath{\sim}{e}^{cLlnL}$ for the decay of a region of large linear size $L$. The implications of our work for the classic problem of the decay of a false vacuum are also discussed.
In a recent paper [Phys. Rev. Lett. 129, 120601 (2022)], we have shown that the dynamics of interfaces, in the symmetry-broken phase of the two-dimensional ferromagnetic quantum Ising model, displays a robust form of ergodicity breaking. In this paper, we elaborate more on the issue. First, we discuss two classes of initial states on the square lattice, the dynamics of which is driven by complementary terms in the effective Hamiltonian and may be solved exactly: (a) Strips of consecutive neighboring spins aligned in the opposite direction of the surrounding spins and (b) a large class of initial states, characterized by the presence of a well-defined ``smooth'' interface separating two infinitely extended regions with oppositely aligned spins. The evolution of the latter states can be mapped onto that of an effective one-dimensional fermionic chain, which is integrable in the infinite-coupling limit. In this case, deep connections with noteworthy results in mathematics emerge, as well as with similar problems in classical statistical physics. We present a detailed analysis of the evolution of these interfaces both on the lattice and in a suitable continuum limit, including the interface fluctuations and the dynamics of entanglement entropy. Second, we provide analytical and numerical evidence supporting the conclusion that the observed nonergodicity---arising from Stark localization of the effective fermionic excitations---persists away from the infinite-Ising-coupling limit, and we highlight the presence of a timescale T\ensuremath∼ecLlnL for the decay of a region of large linear size L. The implications of our work for the classic problem of the decay of a false vacuum are also discussed.
We study the nonequilibrium evolution of coexisting ferromagnetic domains in the two-dimensional quantum Ising model-a setup relevant in several contexts, from quantum nucleation dynamics and false-vacuum decay scenarios to recent experiments with Rydberg-atom arrays. We demonstrate that the quantum-fluctuating interface delimiting a large bubble can be studied as an effective one-dimensional system through a "holographic" mapping. For the considered model, the emergent interface excitations map to an integrable chain of fermionic particles. We discuss how this integrability is broken by geometric features of the bubbles and by corrections in inverse powers of the ferromagnetic coupling, and provide a lower bound to the timescale after which the bubble is ultimately expected to melt. Remarkably, we demonstrate that a symmetry-breaking longitudinal field gives rise to a robust ergodicity breaking in two dimensions, a phenomenon underpinned by Stark many-body localization of the emergent fermionic excitations of the interface.
Understanding the effects of disorder in quantum many-body systems is one of the central challenges of contemporary statistical physics. Its interplay with interactions and quantum fluctuations profoundly modifies the nature of quantum states and their dynamics, giving rise to phenomena such as localization, ergodicity breaking, and the emergence of complex collective behavior. This thesis presents new results on two classes of disordered quantum systems. The first part is devoted to the study of the many-body localization transition, a paradigmatic example of disorder-induced ergodicity breaking. The critical properties of the transition are investigated using a numerical renormalization-group approach, establishing a connection between its scaling behavior and that of the Anderson localization transition on tree-like graphs. The second part focuses on the zero-temperature quantum spin-glass phase of the two-dimensional bond-disordered Heisenberg model. Combining large-scale numerical simulations with a semiclassical expansion, evidence is provided for the existence of the spin-glass phase. The effective model governing the semiclassical fluctuations is further analyzed to characterize the low-energy excitation spectrum and its localization properties. Finally, the weak-disorder regime is investigated and its implications for the ferromagnetic phase are discussed. Overall, the results presented in this thesis contribute to the understanding of the effects of disorder in interacting quantum many-body systems, encompassing both localization phenomena and the collective behavior of disordered quantum magnets.
The theory of statistical mechanics provides a powerful conceptual framework within which the relevant (macroscopic) features of systems at equilibrium can be described. As there is currently no equivalent capable of encompassing the much richer class of non-equilibrium phenomena, research in this direction proceeds mainly on an instance-by-instance basis. The aim of this Thesis is to describe in some detail three such attempts, which involve different dynamical aspects of classical and quantum systems. As summarised below, each of the last three Chapters of this document delves into one of these different topics, while Chapter 2 provides a brief introduction on the study of non-equilibrium dynamics. In Chapter 3 we investigate the purely relaxational dynamics of classical critical ferromagnetic systems in the proximity of surfaces, paying particular attention to the effects that the latter induce on the early stages of the evolution following an abrupt change in the temperature of the sample. When the latter ends close enough to the critical value which separates the paramagnetic from the ferromagnetic phase, it effectively introduces a temporal boundary which can be treated as if it were a surface. Within this picture, we highlight the emergence of novel effects near the effective edge formed by the intersection of the two spatial and temporal boundaries. Our findings are apparently in disagreement with previous predictions which were based on the assumption that the presence of such an edge would not affect the scaling behaviour of observables; in order to explain this discrepancy, we propose an alternative for the original power-counting argument which, at least, correctly predicts the emergence of novel field-theoretical divergences in our one-loop calculations. We show that said singularities are associated with the scaling at the edge. Moreover, by encoding our findings in a boundary renormalisation group framework, we argue that the new predicted behaviour represents a universal feature associated to the short-distance expansion of the order parameter of the transition near the edge; we also calculate explicitly its anomalous dimension at the first-order in a dimensional expansion. As a qualitative feature, this anomalous dimension depends on the type of phase transition occurring at the surface. We exploit this fact in order to provide numerical support to our predictions via Monte Carlo simulations of the dynamical behaviour of a three-dimensional Ising model. The main results reported in Chap. 3 have appeared in Ref. [1]. In Chapter 4 we revisit the Euclidean mapping to imaginary times which has been recently proposed [2, 3] as an alternative for approaching the problem of quantum dynamics following a quench. This is expected to allow one to reformulate the original problem as a static one confined in a film geometry. We show that this interpretation actually holds only if the initial state of the dynamics is pure. Statistical mixtures, instead, intertwine the effects due to the two boundaries, which therefore cannot be regarded as being independent. We emphasize that, although the aforementioned reinterpretation as a confined static problem fails, one is still able, in principle, to write down and solve the corresponding equations. We also discuss in some detail the relation between this approach and the real-time field-theoretical one which makes use of the two-time Keldysh contour. For this purpose, we study the analytical structure of relevant observables — such as correlation functions — in the complex plane of times, identifying a subdivision of this domain into several sectors which depend on the ordering of the imaginary parts of the involved time coordinates. Within each of these subdomains, the analytic continuation to the real axis provides in principle a different result. This feature allows one to reconstruct from the Euclidean formalism all possible non-time-ordered functions, which in particular include all those which can be calculated via the Keldysh two-time formalism. Moreover, we give a prescription on how to retrieve response functions, discussing some simple examples and rationalising some recent numerical data obtained for one of these observables in a one-dimensional quantum Ising chain [4]. We also highlight the emergence of a light-cone effect fairly similar to the one previously found for correlation functions [2], which therefore provides further confirmation to the fact that information travels across the system in the form of the entanglement of quasi-particles produced by the quenching procedure. We have reported part of this analysis in Ref. [5]. Chapter 5 presents part of our recent work on effective relaxation in quantum systems following a quench and on the observed prethermalisation. We analyse the effects caused by the introduction of a long-range integrability-breaking interaction in the early stages of the dynamics of an otherwise integrable quantum spin chain following a quench in the magnetic field. By employing a suitable transformation, we redefine the theory in terms of a fully-connected model of hard-core bosons, which allows us to exploit the (generically) low density of excitations for rendering our model exactly solvable (in a numerical sense, i.e., by numerically diagonalising an exact matrix). We verify that, indeed, as long as the parameters of the quench are not too close to the critical point, the low-density approximation captures the dynamical features of the elementary operators, highlighting the appearance of marked plateaux in their dynamics, which we reinterpret as the emergence of a prethermal regime in the original model. As expected, the latter behaviour is reflected also on extensive observables which can be constructed as appropriate combinations of the mode populations. For these quantities, the typical approach to the quasi-stationary value is algebraic with exponent a ≈ 3, independently of the size of the system, the strength of the interaction and the amplitude of the magnetic field (as long as it is kept far from the critical point). The plateaux mentioned above last until a recurrence time — which can be approximately identified with tR ≈ N/2 for single modes and t′R ≈ N/4 for extensive quantities — after which quantum oscillations due to the finite size of the chain reappear. Our procedure allows us to shed some light over prethermal features without having to considerably limit the size of the system, which we can choose to be quite large, as we discuss in Ref. [6].
We study the out-of-equilibrium dynamics of one-dimensional quantum Ising-like systems, arising from sudden quenches of the Hamiltonian parameter g driving quantum transitions between disordered and ordered phases. In particular, we consider quenches to values of g around the critical value gc, and mainly address the question whether, and how, the quantum transition leaves traces in the evolution of the transverse and longitudinal magnetizations during such a deep out-of-equilibrium dynamics. We shed light on the emergence of singularities in the thermodynamic infinite-size limit, likely related to the integrability of the model. Finite systems in periodic and open boundary conditions develop peculiar power-law finite-size scaling laws related to revival phenomena, but apparently unrelated to the quantum transition, because their main features are generally observed in quenches to generic values of g. We also investigate the effects of dissipative interactions with an environment, modeled by a Lindblad equation with local decay and pumping dissipation operators within the quadratic fermionic model obtainable by a Jordan-Wigner mapping. Dissipation tends to suppress the main features of the unitary dynamics of closed systems. We finally address the effects of integrability breaking, due to further lattice interactions, such as in anisotropic next-to-nearest-neighbor Ising (ANNNI) models. We show that some qualitative features of the post-quench dynamics persist, in particular, the different behaviors when quenching to quantum ferromagnetic and paramagnetic phases, and the revival phenomena due to the finite size of the system.
We analyze the quantum phases, correlation functions and edge modes for a class of spin-1/2 and fermionic models related to the one-dimensional Ising chain in the presence of a transverse field. These models are the Ising chain with anti-ferromagnetic long-range interactions that decay with distance r as 1 / r α , as well as a related class of fermionic Hamiltonians that generalize the Kitaev chain, where both the hopping and pairing terms are long-range and their relative strength can be varied. For these models, we provide the phase diagram for all exponents α , based on an analysis of the entanglement entropy, the decay of correlation functions, and the edge modes in the case of open chains. We demonstrate that violations of the area law can occur for α ≲ 1 , while connected correlation functions can decay with a hybrid exponential and power-law behavior, with a power that is α -dependent. Interestingly, for the fermionic models we provide an exact analytical derivation for the decay of the correlation functions at every α . Along the critical lines, for all models breaking of conformal symmetry is argued at low enough α . For the fermionic models we show that the edge modes, massless for α ≳ 1 , can acquire a mass for α < 1 . The mass of these modes can be tuned by varying the relative strength of the kinetic and pairing terms in the Hamiltonian. Interestingly, for the Ising chain a similar edge localization appears for the first and second excited states on the paramagnetic side of the phase diagram, where edge modes are not expected. We argue that, at least for the fermionic chains, these massive states correspond to the appearance of new phases, notably approached via quantum phase transitions without mass gap closure. Finally, we discuss the possibility to detect some of these effects in experiments with cold trapped ions.
Motivated by recent advances in quantum technologies, we investigate entanglement spreading and excitation transport within quantum complex networks. For this purpose, we adopt a description in which network nodes represent sites of the physical system (e.g. qubits), and links represent pairs of sites coupled by the Hamiltonian [Bentsen et al., Phys Rev Lett (2019)]. We consider a random tunable-range model [Gori et al., Phys Rev E (2017)] which, in terms of only two parameters, encompasses a large variety of graphs and reduces, in limiting cases, to well-studied benchmarks. Then, it is equipped with a paradigmatic Hamiltonian consisting of hopping terms on the network links. Building on known results [Gori et al., Phys Rev E (2017)], we first study the classical case, providing new semi-analytical estimates. We then address quantum dynamics: For the single-particle case, we analyse the asymptotic entanglement entropy between two halves of the system, computed by Exact Diagonalization. Surprisingly, we observe a new localization region for intermediate ranges of the network links. To understand this phenomenon, we introduce an approximate one-dimensional model and try to relate dynamics to network measures. The analysis is extended to half-filling using fermionic Gaussian states, obtaining results close to the single-particle ones that raise new interesting questions for future work. It is concluded that this model exhibits non-trivial behaviour, especially in the region where unexpected localization occurs. We infer that this is determined by the interplay between network features and quantum effects, possibly induced by structural disorder.