Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Data structure of synthetic quantum matter - in and out of equilibrium

Riccardo Andreoni

Published Sep 25, 2026• 0 citations
OpenAlex

Abstract

Analog quantum simulators natively return large collections of single-shot projective measurements, the wave-function snapshots. These are high-dimensional configurations drawn from the Born distribution of the state they prepare. These collections carry far more information than is retained when they are simply compressed into low-order observables fixed in advance. Moreover, this procedure requires knowing beforehand what one is looking for. This thesis chooses to take a different route, asking what structure a dataset of snapshots possesses as a set of points in configuration space, and what kind of physics that structure reflects. The tools used here are drawn from un supervised learning and network theory, namely principal component analysis (PCA), intrinsic-dimension estimators (TWO-NN and BID), and the wave-function networks obtained by promoting snapshots to nodes and activating links according to the met ric structure of configuration space. They are chosen under two constraints, namely interpretability at every step, and the smallest possible amount of physical input: no order parameter, no symmetry, and no knowledge of the Hamiltonian. At equilibrium, wave-function networks are used to build a classification of one dimensional quantum phases based on snapshots alone. The scheme first selects, in the spirit of Occam’s razor, the measurement basis carrying the minimal information content, and then analyses the networks constructed in that basis. Paramagnetic and symmetry-broken phases display low complexity and “probability” networks, symmetry protected topological phases high complexity and Erdős–Rényi networks, and critical points high complexity and scale-free networks. Interestingly, gapped and gapless mat ter are treated on the same footing, and real-space decimation enters the scheme to resolve massive phases rather than flowing towards criticality. Out of equilibrium, the dynamical exponent governing transport of a conserved quan tity is shown to be imprinted on the singular-value spectrum of the snapshot matrix. In particular, it is shown that the PCA entropy, and in fact the leading eigenvalue alone, grows with the same exponent as the underlying spin or energy transport, at modest system sizes, evolution times and sample numbers. Ergodicity breaking is then ad dressed through the binary intrinsic dimension (BID) of the snapshots, in the quantum sun model: the BID grows and saturates at the Thouless time, reaching its maximum in the ergodic phase, remaining system-size independent in the non-ergodic one, and growing sub-maximally at criticality, where the data occupy an extensive but lower dimensional manifold. The associated networks reproduce the equilibrium dichotomy, being random in the ergodic phase and scale-free at the transition, whose location can be identified by the BID within a few percent of spectral estimates. Preliminary experimental observations, not reported here, show analogous features. Taken together, these results indicate that datasets orders of magnitude smaller than the Hilbert-space dimension suffice to identify collective properties of quantum matter through the study of the structure of the data points in configuration space. Moreover, the same structural classes recur across ground states and quenched dynamics, with no order parameter supplied and under realistic experimental imperfections. The final chapter elaborates on another aspect of quantum simulator capabilities, that is, that to realize gauge theory dynamics. There, one-dimensional U(1) lattice gauge theories with a topological θ term are mapped exactly onto kinetically con strained fermions in an emergent Stark potential. The mapping accounts for the slow relaxation reported in recent Rydberg experiments: prethermalization near an inte grable point at θ = 0, and Stark localization with Hilbert-space fragmentation at θ ̸ =0.

Related papers

Network theory classification of quantum matter based on wave function snapshots

Riccardo Andreoni, Vittorio Vitale, Cristiano Muzzi +3 more

Quantum computers and simulators offer unparalleled capabilities of probing quantum many-body states, by obtaining snapshots of the many-body wave function via collective projective measurements. The probability distribution obtained by such snapshots (which are fundamentally limited to a negligible fraction of the Hilbert space) is of fundamental importance to determine the power of quantum computations. However, its relation to many-body collective properties is poorly understood. Here, we develop a theoretical framework to link quantum phases of matter to their snapshots, based on a combination of data complexity and network theory analyses. The first step in our scheme consists of applying Occam's razor principle to quantum sampling: given snapshots of a wave function, we identify a minimal-complexity measurement basis by analyzing the information compressibility of snapshots over different measurement bases. The second step consists of analyzing arbitrary correlations using network theory, building a wave-function network from the minimal-complexity basis data. This approach allows us to stochastically classify the output of quantum computers and simulations, with no assumptions on the underlying dynamics, and in a fully interpretable manner. We apply this method to quantum states of matter in one-dimensional translational invariant systems, where such classification is exhaustive, and where it reveals an interesting interplay between algorithmic and computational complexity for many-body states. Our framework is of immediate experimental relevance, and can be further extended both in terms of more advanced network mathematics, including discrete homology, as well as in terms of applications to physical phenomena, such as time-dependent dynamics and gauge theories.

quant-pharXiv:2512.02121•Dec 3, 2025•0 citations

Many-body ergodicity breaking from wavefunction snapshots

Riccardo Andreoni, Devendra Singh Bhakuni, Marcello Dalmonte +2 more

Modern quantum experiments can probe many-body wavefunctions through projective measurements, providing snapshots of individual many-body configurations. A fundamental question is whether the intrinsic structure of their measurement distributions can reveal quantum dynamics beyond predefined observables. Here, we address this question at the many-body ergodicity-breaking transition using two complementary characteristics of nonequilibrium wavefunction snapshots: their binary intrinsic dimension (BID) and the topology of wavefunction networks constructed from Hamming distances. We show that the critical point is characterized by an extensive but submaximal BID, and scale-free network connectivity with highly connected hubs. These signatures distinguish the ergodic, critical, and nonergodic regimes, and provide experimentally accessible probes of ergodicity breaking from projective measurement snapshots.

quant-pharXiv:2608.21524•Aug 25, 2026•0 citations

Diagnosing quantum transport from wave function snapshots

Devendra Singh Bhakuni, Roberto Verdel, Cristiano Muzzi +3 more

We study nonequilibrium quantum dynamics of spin chains by employing principal component analysis on data sets of wave function snapshots and examine how information propagates within these data sets. The quantities we employ are derived from the spectrum of the sample second moment matrix, built directly from data sets. Our investigations on several interacting spin chains featuring distinct spin or energy transport reveal that the growth of data information spreading follows the same dynamical exponents as that of the underlying quantum transport of spin or energy. Specifically, our approach enables an easy, data-driven, and, importantly, interpretable diagnostic to track energy transport with a limited number of samples, which is usually challenging without any assumption on the Hamiltonian form. These observations are obtained at a modest finite-size and evolution time, which aligns with experimental and numerical constraints. Our framework directly applies to experimental quantum simulator data sets of dynamics in higher-dimensional systems, where classical simulation methods usually face significant limitations and apply equally to both near- and far-from-equilibrium quenches.

arXiv:2407.09092•Physical review. B./Physical review. B•Oct 21, 2024•2 citations

Intrinsic Dimension of Path Integrals: Data-Mining Quantum Criticality and Emergent Simplicity

T. Mendes-Santos, Adriano Angelone, Álex Rodríguez +2 more

Quantum many-body systems are characterized by patterns of correlations defining highly nontrivial manifolds when interpreted as data structures.Physical properties of phases and phase transitions are typically retrieved via correlation functions, that are related to observable response functions.Recent experiments have demonstrated capabilities to fully characterize quantum many-body systems via wavefunction snapshots, opening new possibilities to analyze quantum phenomena.Here, we introduce a method to data mine the correlation structure of quantum partition functions via their path integral (or equivalently, stochastic series expansion) manifold.We characterize path-integral manifolds generated via state-of-the-art quantum Monte Carlo methods utilizing the intrinsic dimension (ID) and the variance of distances between nearest-neighbor (NN) configurations: the former is related to data-set complexity, while the latter is able to diagnose connectivity features of points in configuration space.We show how these properties feature universal patterns in the vicinity of quantum criticality, that reveal how data structures simplify systematically at quantum phase transitions.This is further reflected by the fact that both ID and variance of NN distances exhibit universal scaling behavior in the vicinity of second-order and Berezinskii-Kosterlitz-Thouless critical points.Finally, we show how non-Abelian symmetries dramatically influence quantum data sets, due to the nature of (noncommuting) conserved charges in the quantum case.Complementary to neural-network representations, our approach represents a first elementary step towards a systematic characterization of path-integral manifolds before any dimensional reduction is taken, that is informative about universal behavior and complexity, and can find immediate application to both experiments and Monte Carlo simulations.

arXiv:2103.02640•PRX Quantum•Aug 23, 2021•34 citations

A study on non-equilibrium dynamics in classical and quantum systems

Matteo Marcuzzi

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].

•Oct 30, 2013•0 citations

Probing Complexity, Universality and Dimensional Reduction in Nonequilibrium Statistical Physics

Cristiano Muzzi

This thesis investigates universality and complexity in classical and quantum statistical physics, with a particular emphasis on nonequilibrium phenomena. While the emergence of universal behavior in equilibrium systems stands as one of the cornerstones of statistical physics, the extension of these ideas to nonequilibrium settings presents fundamental and still open challenges. The inherently diverse and system-specific nature of nonequilibrium dynamics calls for a broad spectrum of conceptual and analytical tools. Furthermore, nonequilibrium dynamics is nowadays accessible to experiments in both classical and quantum regimes. These modern experiments produce a large amount of experimental data, which calls for analysis of physics through data-driven methods, complexity theory and high-dimensional statistics. The possibility of probing and accessing correlated states of matter out of equilibrium in modern experiments raises several questions. Can we understand the output of experiments in terms of complexity theory and dimensional reduction? Is it possible to extend concepts of universality to data structures arising in simulations? Can we describe phase transitions (either classical or quantum) concepts of modern data science? How to control entanglement in quantum dynamics? Which nonequilibrium phenomena can we hope to observe in quantum simulators? The research work presented in this thesis is an attempt to answer to these questions. In this thesis, we argue that order parameters and diagnostics of criticality in nonequilibrium classical systems can be found by unsupervised learning and dimensional reduction, which allow to reduce the complexity of datasets by projections on lower-dimensional spaces, while still preserving essential universal features. We then move to the study of models of quantum dynamics relevant for quantum simulators. In particular, we show that the same concepts of dimensional reduction apply also to quantum systems, allowing for the extraction of information from outputs of quantum simulators. We then study a protocol to prepare complex, entangled states via unitary evolution and enhance entanglement via measurements in a fermionic model. Finally, we study the emergence of anomalous transport in one-dimensional chiral systems with non-Abelian symmetries.

•Sep 22, 2025•0 citations
Analog quantum simulators natively return large collections of single-shot projective measurements, the wave-function snapshots. These are high-dimensional configurations drawn from the Born distribution of the state they prepare. These collections carry far more information than is retained when they are simply compressed into low-order observables fixed in advance. Moreover, this procedure requires knowing beforehand what one is looking for. This thesis chooses to take a different route, asking what structure a dataset of snapshots possesses as a set of points in configuration space, and what kind of physics that structure reflects. The tools used here are drawn from un supervised learning and network theory, namely principal component analysis (PCA), intrinsic-dimension estimators (TWO-NN and BID), and the wave-function networks obtained by promoting snapshots to nodes and activating links according to the met ric structure of configuration space. They are chosen under two constraints, namely interpretability at every step, and the smallest possible amount of physical input: no order parameter, no symmetry, and no knowledge of the Hamiltonian. At equilibrium, wave-function networks are used to build a classification of one dimensional quantum phases based on snapshots alone. The scheme first selects, in the spirit of Occam’s razor, the measurement basis carrying the minimal information content, and then analyses the networks constructed in that basis. Paramagnetic and symmetry-broken phases display low complexity and “probability” networks, symmetry protected topological phases high complexity and Erdős–Rényi networks, and critical points high complexity and scale-free networks. Interestingly, gapped and gapless mat ter are treated on the same footing, and real-space decimation enters the scheme to resolve massive phases rather than flowing towards criticality. Out of equilibrium, the dynamical exponent governing transport of a conserved quan tity is shown to be imprinted on the singular-value spectrum of the snapshot matrix. In particular, it is shown that the PCA entropy, and in fact the leading eigenvalue alone, grows with the same exponent as the underlying spin or energy transport, at modest system sizes, evolution times and sample numbers. Ergodicity breaking is then ad dressed through the binary intrinsic dimension (BID) of the snapshots, in the quantum sun model: the BID grows and saturates at the Thouless time, reaching its maximum in the ergodic phase, remaining system-size independent in the non-ergodic one, and growing sub-maximally at criticality, where the data occupy an extensive but lower dimensional manifold. The associated networks reproduce the equilibrium dichotomy, being random in the ergodic phase and scale-free at the transition, whose location can be identified by the BID within a few percent of spectral estimates. Preliminary experimental observations, not reported here, show analogous features. Taken together, these results indicate that datasets orders of magnitude smaller than the Hilbert-space dimension suffice to identify collective properties of quantum matter through the study of the structure of the data points in configuration space. Moreover, the same structural classes recur across ground states and quenched dynamics, with no order parameter supplied and under realistic experimental imperfections. The final chapter elaborates on another aspect of quantum simulator capabilities, that is, that to realize gauge theory dynamics. There, one-dimensional U(1) lattice gauge theories with a topological θ term are mapped exactly onto kinetically con strained fermions in an emergent Stark potential. The mapping accounts for the slow relaxation reported in recent Rydberg experiments: prethermalization near an inte grable point at θ = 0, and Stark localization with Hilbert-space fragmentation at θ ̸ =0.
Quantum computers and simulators offer unparalleled capabilities of probing quantum many-body states, by obtaining snapshots of the many-body wave function via collective projective measurements. The probability distribution obtained by such snapshots (which are fundamentally limited to a negligible fraction of the Hilbert space) is of fundamental importance to determine the power of quantum computations. However, its relation to many-body collective properties is poorly understood. Here, we develop a theoretical framework to link quantum phases of matter to their snapshots, based on a combination of data complexity and network theory analyses. The first step in our scheme consists of applying Occam's razor principle to quantum sampling: given snapshots of a wave function, we identify a minimal-complexity measurement basis by analyzing the information compressibility of snapshots over different measurement bases. The second step consists of analyzing arbitrary correlations using network theory, building a wave-function network from the minimal-complexity basis data. This approach allows us to stochastically classify the output of quantum computers and simulations, with no assumptions on the underlying dynamics, and in a fully interpretable manner. We apply this method to quantum states of matter in one-dimensional translational invariant systems, where such classification is exhaustive, and where it reveals an interesting interplay between algorithmic and computational complexity for many-body states. Our framework is of immediate experimental relevance, and can be further extended both in terms of more advanced network mathematics, including discrete homology, as well as in terms of applications to physical phenomena, such as time-dependent dynamics and gauge theories.
Modern quantum experiments can probe many-body wavefunctions through projective measurements, providing snapshots of individual many-body configurations. A fundamental question is whether the intrinsic structure of their measurement distributions can reveal quantum dynamics beyond predefined observables. Here, we address this question at the many-body ergodicity-breaking transition using two complementary characteristics of nonequilibrium wavefunction snapshots: their binary intrinsic dimension (BID) and the topology of wavefunction networks constructed from Hamming distances. We show that the critical point is characterized by an extensive but submaximal BID, and scale-free network connectivity with highly connected hubs. These signatures distinguish the ergodic, critical, and nonergodic regimes, and provide experimentally accessible probes of ergodicity breaking from projective measurement snapshots.
We study nonequilibrium quantum dynamics of spin chains by employing principal component analysis on data sets of wave function snapshots and examine how information propagates within these data sets. The quantities we employ are derived from the spectrum of the sample second moment matrix, built directly from data sets. Our investigations on several interacting spin chains featuring distinct spin or energy transport reveal that the growth of data information spreading follows the same dynamical exponents as that of the underlying quantum transport of spin or energy. Specifically, our approach enables an easy, data-driven, and, importantly, interpretable diagnostic to track energy transport with a limited number of samples, which is usually challenging without any assumption on the Hamiltonian form. These observations are obtained at a modest finite-size and evolution time, which aligns with experimental and numerical constraints. Our framework directly applies to experimental quantum simulator data sets of dynamics in higher-dimensional systems, where classical simulation methods usually face significant limitations and apply equally to both near- and far-from-equilibrium quenches.
Quantum many-body systems are characterized by patterns of correlations defining highly nontrivial manifolds when interpreted as data structures.Physical properties of phases and phase transitions are typically retrieved via correlation functions, that are related to observable response functions.Recent experiments have demonstrated capabilities to fully characterize quantum many-body systems via wavefunction snapshots, opening new possibilities to analyze quantum phenomena.Here, we introduce a method to data mine the correlation structure of quantum partition functions via their path integral (or equivalently, stochastic series expansion) manifold.We characterize path-integral manifolds generated via state-of-the-art quantum Monte Carlo methods utilizing the intrinsic dimension (ID) and the variance of distances between nearest-neighbor (NN) configurations: the former is related to data-set complexity, while the latter is able to diagnose connectivity features of points in configuration space.We show how these properties feature universal patterns in the vicinity of quantum criticality, that reveal how data structures simplify systematically at quantum phase transitions.This is further reflected by the fact that both ID and variance of NN distances exhibit universal scaling behavior in the vicinity of second-order and Berezinskii-Kosterlitz-Thouless critical points.Finally, we show how non-Abelian symmetries dramatically influence quantum data sets, due to the nature of (noncommuting) conserved charges in the quantum case.Complementary to neural-network representations, our approach represents a first elementary step towards a systematic characterization of path-integral manifolds before any dimensional reduction is taken, that is informative about universal behavior and complexity, and can find immediate application to both experiments and Monte Carlo simulations.
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 investigates universality and complexity in classical and quantum statistical physics, with a particular emphasis on nonequilibrium phenomena. While the emergence of universal behavior in equilibrium systems stands as one of the cornerstones of statistical physics, the extension of these ideas to nonequilibrium settings presents fundamental and still open challenges. The inherently diverse and system-specific nature of nonequilibrium dynamics calls for a broad spectrum of conceptual and analytical tools. Furthermore, nonequilibrium dynamics is nowadays accessible to experiments in both classical and quantum regimes. These modern experiments produce a large amount of experimental data, which calls for analysis of physics through data-driven methods, complexity theory and high-dimensional statistics. The possibility of probing and accessing correlated states of matter out of equilibrium in modern experiments raises several questions. Can we understand the output of experiments in terms of complexity theory and dimensional reduction? Is it possible to extend concepts of universality to data structures arising in simulations? Can we describe phase transitions (either classical or quantum) concepts of modern data science? How to control entanglement in quantum dynamics? Which nonequilibrium phenomena can we hope to observe in quantum simulators? The research work presented in this thesis is an attempt to answer to these questions. In this thesis, we argue that order parameters and diagnostics of criticality in nonequilibrium classical systems can be found by unsupervised learning and dimensional reduction, which allow to reduce the complexity of datasets by projections on lower-dimensional spaces, while still preserving essential universal features. We then move to the study of models of quantum dynamics relevant for quantum simulators. In particular, we show that the same concepts of dimensional reduction apply also to quantum systems, allowing for the extraction of information from outputs of quantum simulators. We then study a protocol to prepare complex, entangled states via unitary evolution and enhance entanglement via measurements in a fermionic model. Finally, we study the emergence of anomalous transport in one-dimensional chiral systems with non-Abelian symmetries.