The laws of nature encompass the small, the large, the few, and the many. In this thesis we are concerned with classical (i.e. not quantum) many-body systems, which refers to any microscopic or macroscopic system that contains a large number of interacting entities. The nearest-neighbor Ising model, originally developed in 1920 by Wilhelm Lenz, forms a cornerstone in our theoretical understanding of collective effects in classical many-body systems, and is to date a paradigm for statistical physics. Despite its elegant and simplistic description, exact analytical results in dimensions equal and larger than two are difficult to obtain. Therefore, much work has been done to con- struct methods that allow for approximate, yet accurate, analytical solutions. One of these methods is the Bethe-Guggenheim approximation, originally developed independently by Hans Bethe and Edward Guggenheim in 1935. This approximation goes beyond the well- known mean field approximation, and explicitly accounts for pair correlations between the spins in the Ising model. In this thesis, we embark on a journey to exploit the full capacity of the Bethe-Guggenheim approximation, in non-uniform and non-equilibrium settings. After we formally introduce the original Bethe-Guggenheim approximation for uniform systems, we will extend its scope to non-uniform systems, and derive a Cahn-Hilliard free energy functional. Here we find that the one-dimensional equilibrium concentration profile undergoes a delocalization-induced broadening transition at interaction strengths near and above the thermal energy. The broadening transition arises from a decreasing amplitude of capillary wave fluctuations, and is not accounted for in the mean field approximation. Finally, going beyond equilibrium properties, we also study the kinetics of the Ising model in and out of equilibrium. First, based on firm theoretical and experimental evidence, we con- struct an Ising-like minimal model for interacting cellular adhesion molecules. Here we find a dynamical critical point where adhesion cluster dissolution and formation are the fastest, and undergo a qualitative change in dynamics. Second, we demonstrate the existence of a finite-time dynamical phase transition for disordering quenches in the nearest-neighbor Ising model. Starkly different from the mean field approximation, the time at which the dynamical phase transition occurs, is bounded from below by a speed limit. Altogether, in this thesis, we unveil the non-trivial and a priori non-intuitive effects of pair correlations in the classical nearest-neighbor Ising model, which are taken into account in the Bethe-Guggenheim approximation, and neglected in the mean field approximation.
The laws of nature encompass the small, the large, the few, and the many. In this thesis we are concerned with classical (i.e. not quantum) many-body systems, which refers to any microscopic or macroscopic system that contains a large number of interacting entities. The nearest-neighbor Ising model, originally developed in 1920 by Wilhelm Lenz, forms a cornerstone in our theoretical understanding of collective effects in classical many-body systems, and is to date a paradigm for statistical physics. Despite its elegant and simplistic description, exact analytical results in dimensions equal and larger than two are difficult to obtain. Therefore, much work has been done to con- struct methods that allow for approximate, yet accurate, analytical solutions. One of these methods is the Bethe-Guggenheim approximation, originally developed independently by Hans Bethe and Edward Guggenheim in 1935. This approximation goes beyond the well- known mean field approximation, and explicitly accounts for pair correlations between the spins in the Ising model. In this thesis, we embark on a journey to exploit the full capacity of the Bethe-Guggenheim approximation, in non-uniform and non-equilibrium settings. After we formally introduce the original Bethe-Guggenheim approximation for uniform systems, we will extend its scope to non-uniform systems, and derive a Cahn-Hilliard free energy functional. Here we find that the one-dimensional equilibrium concentration profile undergoes a delocalization-induced broadening transition at interaction strengths near and above the thermal energy. The broadening transition arises from a decreasing amplitude of capillary wave fluctuations, and is not accounted for in the mean field approximation. Finally, going beyond equilibrium properties, we also study the kinetics of the Ising model in and out of equilibrium. First, based on firm theoretical and experimental evidence, we con- struct an Ising-like minimal model for interacting cellular adhesion molecules. Here we find a dynamical critical point where adhesion cluster dissolution and formation are the fastest, and undergo a qualitative change in dynamics. Second, we demonstrate the existence of a finite-time dynamical phase transition for disordering quenches in the nearest-neighbor Ising model. Starkly different from the mean field approximation, the time at which the dynamical phase transition occurs, is bounded from below by a speed limit. Altogether, in this thesis, we unveil the non-trivial and a priori non-intuitive effects of pair correlations in the classical nearest-neighbor Ising model, which are taken into account in the Bethe-Guggenheim approximation, and neglected in the mean field approximation.
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].
This thesis is concerned with non-equilibrium phenomena in interacting quantum many-body systems. Specifically, we investigate the time-evolution and relaxation dynamics of Bose gases in a highly restricted geometry, which constrains the dynamics to one spatial dimension. This leads to a description of the system in terms of a simple model Hamiltonian which permits exact many-body quantum mechanical solutions due to its integrability. In the first part of this thesis, a computational method is developed to obtain experimentally relevant correlation functions in the framework of the coordinate Bethe ansatz. We employ this method to compute exact ground-state correlation functions of the Lieb-Liniger gas for up to seven particles covering the whole regime of repulsive interactions. We also investigate the dynamics of the system after an instantaneous change of the interaction strength. This quantum quench deposits large amounts of energy that cannot be dissipated due to the system being closed, and so the dynamics far from equilibrium are probed. We prepare the system in two different initial states, and quench to the same final interaction strength. The latter is determined in such a way that the added energy due to the quench is the same for both scenarios. Conventional statistical mechanics predicts the same relaxed state, but due to the integrability of the system all considered correlation functions of the relaxed states differ from each other and also from the thermal ones. We then investigate the dynamics and relaxed state for a quench from zero to repulsive interactions in more detail, focussing on the mechanism of relaxation and the involved time-scales. We find that local correlation functions relax on time-scales determined by the interaction strength, in contrast to non-local correlation functions, whose relaxation time-scale is proportional to the system size. Next, we employ the same methodology to study the one-dimensional Bose gas with attractive interactions. In this case many-body bound states are permissible solutions of the Lieb-Liniger model. We compare exact ground-state correlation functions of up to seven particles to their corresponding mean-field solution. The latter displays a quantum phase transition at a critical interaction strength, marking the transition from a uniform-density state to a localized bright soliton. Our exact results agree remarkably well with the corresponding mean-field solution past the critical point. We also investigate the dynamics following an interaction strength quench, starting again from the non-interacting ground state. Bound states strongly influence correlation functions for all post-quench interaction strengths, and local correlation functions are largely increased compared to their initial value. In the last part of this thesis, we investigate the behavior of the one-dimensional Bose gas under periodic driving of the interaction strength. To this end, we extend the coordinate Bethe-ansatz formalism by employing Floquet theory to obtain solutions for the full time- dependent Hamiltonian. This realizes an extension to non-integrable systems in a way that allows controlled breaking of integrability. We compare the dynamics of the system following an interaction quench to a fixed final interaction strength to that of a system with periodically modulated post-quench interaction strength. For fast driving, the system evolves according to the time-averaged Hamiltonian and correlation functions are nearly identical to those of the undriven system. The response of the system qualitatively changes at resonances of the time-averaged Hamiltonian, where we observe energy absorption and a marked change in correlation functions. However, the system does not absorb energy indefinitely.
Recent works unraveled an intriguing finite-time dynamical phase transition in the thermal relaxation of the mean field Curie-Weiss model. The phase transition reflects a sudden switch in the dynamics. Its existence in systems with a finite range of interaction, however, remained unclear. Here we demonstrate the dynamical phase transition for nearest-neighbor Ising systems on the square and Bethe lattices through extensive computer simulations and by analytical results. Combining large-deviation techniques and Bethe-Guggenheim theory we prove the existence of the dynamical phase transition for arbitrary quenches, including those within the two-phase region. Strikingly, for any given initial condition we prove and explain the existence of non-trivial speed limits for the dynamical phase transition and the relaxation of magnetization, which are fully corroborated by simulations of the microscopic Ising model but are absent in the mean field setting. Pair correlations, which are neglected in mean field theory and trivial in the Curie-Weiss model, account for kinetic constraints due to frustrated local configurations that give rise to a global speed limit.
Recent advances in theoretical physics have increasingly focused on the dynamics of quantum systems with many interacting degrees of freedom, while ultracold atom platforms now enable the realization of diverse theoretical models with precise control over system parameters and advanced monitoring techniques. A prominent example is cavity quantum electrodynamics (CQED), where large atomic or molecular ensembles are coupled to high-finesse optical cavities, enabling strong interaction with confined electromagnetic modes and leading to phenomena such as the superradiant phase transition, in which atoms and photons form a coherent state. Most theoretical studies of CQED have concentrated on regimes of spatially uniform atom–cavity coupling, where collective behavior allows a description in terms of the classical dynamics of a few macroscopic variables. However, experimental capabilities now permit exploration beyond this limit, though theoretical progress is hindered by the exponential computational complexity of simulating large quantum systems, particularly when environmental coupling is included. This thesis addresses such challenges by developing efficient approaches, within the Keldysh formalism of non-equilibrium quantum field theory, to study systems with strongly disordered interactions coupled to external reservoirs. Two representative models are considered: a disordered spin–boson model describing light–matter interactions with quenched disorder, relevant to recent CQED experiments in confocal cavities, and a fermion–phonon system governed by the Yukawa–Sachdev–Ye–Kitaev (Yukawa-SYK) model. Both systems are quantum critical, exhibiting a continuum of low-energy modes above their stationary states, yet disorder produces qualitatively distinct dynamical behaviors: in the spin–boson system it induces slow, glassy relaxation, whereas in the Yukawa-SYK model it leads to rapid thermalization. These contrasting outcomes reveal the non-trivial role of disorder in quantum many-body systems and motivate further investigation into its effects on non-equilibrium quantum dynamics.
The study of quantum dynamics is an increasingly important aspect of atomic, molecular and optical physics, as well as condensed-matter physics and other disciplines. It is essential in order to understand non-equilibrium phenomena which do not fit within the framework of equilibrium thermodynamics, such as many-body localization, dynamical quantum phase transitions, and non-equilibrium steady states in driven-dissipative systems. However, the exact calculation of dynamical observables is extremely difficult in general for several reasons, which makes it challenging for theoretical predictions and analyses to keep pace with rapidly advancing experiments. In this thesis I use and develop several computational methods to investigate the dynamics of quantum many-body systems, with the two goals of better understanding the evolution of local observables and correlation functions in spin and Fermion systems as well as creating and testing more effective numerical methods for solving quantum dynamics problems. I explore quantum dynamics by using Hamiltonian parameter quenches in several commonly-used models in AMO physics and elsewhere: the Fermi-Hubbard model, the XXZ model, and the transverse-Ising model. The first represents a system of short-range interacting Fermions, while the latter two represent spin systems with different spin-spin couplings and magnetic field terms. All of these models can be realized in current experiments. First, I derive an exact solution for an extreme interaction quench in the Fermi-Hubbard model, from a strongly-interacting initial state to a noninteracting final Hamiltonian. From this result I show that there are nontrivial transient connected correlations, despite the fact that the system has very high temperature and a noninteracting Hamiltonian. Those features would suppress such correlations in equilibrium, but the dynamical correlations are still significant. Then, I develop and examine two numerical methods for solving dynamics of the spin expectation values and correlation functions in the transverse-Ising and XXZ models. Both methods rely on approximating a thermodynamically large system with a much smaller and more computationally tractable cluster of lattice sites, or a series of clusters. The first method is a dynamical numerical linked cluster expansion (NLCE). NLCEs have already been used successfully on systems in thermodynamic equilibrium, and I extended the technique to directly compute the dynamics of an observable after a quantum quench. I review the general NLCE method, describe my dynamical variant of it, and compute the quench dynamics of spin and correlation observables in the transverse-Ising and XXZ models. The second method, which I term the adaptive-boundary cluster method, relies on recently refined bounds on information propagation in quantum systems to choose a cluster of lattice sites that minimizes the finite size error bound of a dynamical computation. I briefly review these bounds, and describe a new algorithm which uses them to construct a cluster whose boundary is chosen according to a related error estimate. Such an adaptive-boundary cluster has the property that there exists a computable upper bound on the finite size error of observables measured at the cluster center, and this error bound is nearly minimal when compared to that of other clusters of the same size. I then compute the quench dynamics of spin and correlation observables in the transverse-Ising model using the adaptive boundary clusters. I apply this method to both uniform and strongly disordered systems. To evaluate the effectiveness of both dynamical NLCE and adaptive-boundary clusters, I compare the results obtained using them to results from a standard and general purpose numerical method: exact solvers applied a rectangular cluster with periodic boundary conditions. Both of the new methods usually converge faster and are more accurate than the standard method for comparable cluster sizes. This thesis illustrates that choosing numerical methods which take advantage of the geometric and entanglement structure of the initial state and Hamiltonian can substantially enhance our ability to study quantum dynamics, by reducing the computational cost of numerically exact solutions. I conclude with a discussion of various avenues for future research to build on this work.
The study of non-equilibrium behaviors of quantum matter is one of the research directions at the frontiers of condensed matter physics and, in general, of theoretical and experimental physics. Understanding the mechanisms which prevent many quantum particles from reaching thermal equilibrium can ultimately shed light on the foundations of statistical physics and provide insight into the anomalous long-time collective behaviors observed in experiments. In this thesis, we present new theoretical developments in this direction. In the first part of the thesis we discuss the non-equilibrium dynamics of the two-dimensional quantum Ising model, the prototypical model of magnetism. In particular, we investigate the dynamics of interfaces separating different spin domains. In some regimes, we can provide an exact description of their unitary time evolution, but we also discuss the robustness of our predictions away from the exactly-solvable limit. Moreover, we address the same problem in the presence of disorder and provide a hydrodynamic interpretation of our findings. In the second part of the thesis, we investigate the localization properties of disordered quantum systems. In particular, we develop a renormalization-group picture of the Anderson localization transition in generic geometries. This provides a new framework to understand the phase diagram of such systems and their critical behavior beyond the usual finite-size scaling analysis. We also provide further insights on this topic both by studying resonances and their connection with observables and by employing unsupervised non-parametric learning techniques to probe disordered systems in an experimentally accessible form. The results presented in this thesis constitute a step forward in understanding the non-equilibrium phases and regimes of quantum matter and they pave the way to interesting future developments.