We investigate the effects of fluctuations on the dynamics of an isolated quantum system represented by a ${\left({\mathbf{\ensuremath{\phi}}}^{2}\right)}^{2}$ field theory with
O(N) symmetry after a quench in
d>2 spatial dimensions. A perturbative renormalization-group approach involving a dimensional expansion in $\ensuremath{\epsilon}=4\ensuremath{-}d$ is employed in order to study the evolution within a prethermal regime controlled by elastic dephasing. In particular, we focus on a quench from a disordered initial state to the critical point, which introduces an effective temporal boundary in the evolution. At this boundary, the relevant fields acquire an anomalous scaling dimension, while the evolution of both the order parameter and its correlation and response functions display universal aging. Since the relevant excitations propagate ballistically, a light cone in real space emerges. At longer times, the onset of inelastic scattering appears as secularly growing self-generated dissipation in the effective Keldysh field theory, with the strength of the dissipative vertices providing an estimate for the time needed to leave the prethermal regime.
We study correlations in space and time after a quantum quench in the Bose Hubbard model. We obtain a two particle irreducible effective action within the Schwinger-Keldysh formalism that allows us to study dynamics in both the superfluid and Mott insulating phases. We derive equations of motion for both the superfluid order parameter and two-point correlation functions. We discuss the numerical solutions of these equations and calculate the quasimomentum distribution for quenches in the Mott phase. We demonstrate light-cone like spreading of correlations in the Mott phase in both one and two dimensions and calculate the propagation velocity in each dimension.
We study a quantum quench in the Bose-Hubbard model where the tunneling rate J is suddenly switched from zero to a finite value in the Mott regime. In order to solve the many-body quantum dynamics far from equilibrium, we consider the reduced density matrices for a finite number of lattice sites and split them up into on-site density operators, i.e., the mean field, plus two-point and three-point correlations etc. Neglecting three-point and higher correlations, we are able to numerically simulate the time-evolution of the on-site density matrices and the two-point quantum correlations (e.g., their effective light-cone structure) for a comparably large number of lattice sites.
The response of physical systems to external perturbations can be used to probe both their equilibrium and nonequilibrium dynamics. While response and correlation functions are related in equilibrium by fluctuation-dissipation theorems, out of equilibrium they provide complementary information on the dynamics. In the past years, a method has been devised to map the quantum dynamics of an isolated extended system after a quench onto a static theory with boundaries in imaginary time; up to now, however, the focus was entirely on symmetrized correlation functions. Here we provide a prescription which, in principle, allows one to retrieve the whole set of relevant dynamical quantities characterizing the evolution, including linear response functions. We illustrate this construction with some relevant examples, showing in the process the emergence of light-cone effects similar to those observed in correlation functions.
We investigate the time evolution of density, spin, and pairing correlation functions in one-dimensional $t\text{\ensuremath{-}}J$ models following a quantum quench using the time-dependent density matrix renormalization group. While density and spin correlation functions show the typical light-cone behavior over a wide range of parameters, in pairing correlation functions it is strongly suppressed. This is supported by time-dependent BCS theory, where the light cone in the pairing correlation functions is found to be at least two orders of magnitude weaker than in the density correlator. These findings indicate that in global quantum quenches not all observables are affected equally by the excitations induced by the quench.
In recent years the out-of-equilibrium dynamics of many-body quantum systems has attracted growing interest, both for its fundamental implications in the understanding of quantum thermalization and for potential applications in quantum information and quantum computation. In particular, quantum quenches provide a natural setting to investigate the mechanisms of correlation propagation and entanglement growth far from equilibrium. Interest in these phenomena stems not only from the possibility of testing theoretical predictions coming from emergent hydrodynamic descriptions, but also from understanding whether, and to what extent, such descriptions remain valid in the presence of spatial gradients or interfaces. Two paradigmatic frameworks in this context are the Generalized Hydrodynamics (GHD) and quasi-particle picture. The former is an effective theory devised to describe the long-time evolution after a quench in the thermodynamic limit. Within this perspective, inhomogeneous quenches are a crucial testbed for extending hydrodynamic theories to more realistic situations where, for example, scattering at an interface can occur. This is the main goal of the present thesis. The reference physical model is the one-dimensional Kitaev chain, an integrable spinless p-wave superconductor that exhibits both a topological and a trivial phase. The bilinear structure of the Hamiltonian in the fermionic creation and annihilation operators admits an exact rewriting in terms of Majorana fermions, in which the Hamiltonian remains quadratic. This representation allows a compact description of the ground-state properties via the Majorana correlation matrix (Gamma) at time t=0: the entries of Gamma completely encode the quantum information of the Gaussian ground state. The time evolution of Gamma is studied in the Heisenberg picture, where unitary dynamics correspond to an orthogonal rotation of Gamma implemented by an orthogonal matrix determined by the Kitaev Hamiltonian in the Majorana basis. Consequently, the entanglement entropy dynamics following a quench, both homogeneous and inhomogeneous, can be simulated through simple matrix products, making the approach numerically efficient even for sizable systems. The evolution code for Gamma was validated on a homogeneous quench with periodic boundary conditions: subsystem-size (l) curves exhibit the expected finite-size scaling, (Δ \propto 1/l), and numerical growth rates in the linear regime (t < v/2l), extracted by linear fits at fixed l, match the quasi-particle predictions for various values of the new chemical potential. The rate discrepancy behaves asymptotically as (Δ_{rate}} \propto 1/t_{\min}\) for large (t_{\min}\), up to oscillations from sub-extensive terms. The physical problem of interest is then an inhomogeneous quench of the chemical potential on the infinite Kitaev chain ground state with open boundaries, (μ_0 to μ_L, μ_R), producing a central interface. Initial and final Hamiltonians are free in the quasiparticle sense (ballistic propagation), but the presence of two distinct energy bands implies interface scattering. Simulations varying only μ_R (within the same trivial phase) show a nonmonotonic dependence of the entanglement growth rate, with a maximum at a particular μ_R. To explain this, an analytic ansatz for the stationary Majorana wavefunction was constructed: plane waves multiplied by a Majorana spinor, with bulk spinors obtained from the homogeneous PBC problem via the time-independent Schrödinger equation. Imposing matching at the interface yields four nontrivial equations for two unknowns (overdetermined in general). In the special cases (Δ ̸= t) two equations become trivial, the system is determined, and one obtains an explicit transmission coefficient. Combining this result with known quasi-particle formulas for two semi-infinite chains yields a hydrodynamic prediction for entanglement growth with interface scattering. Energy conservation across the interface selects the momenta that can be transmitted and thus contribute to inter-region correlations. Numerical tests for several total lengths (L) confirm the finite-size scaling (Δ_{\mathrm{e.e.}} \propto 1/L\), and show excellent agreement between numerical and theoretical growth rates for both quenches restricted to the right half and quenches applied to both halves. For the right-only case the rate discrepancy again follows (Δ_{rate} \propto 1/t_{min}\) as
tmin→∞. Possible future directions include: the study of the general inhomogeneous quench with (Δ ̸= t) via a direct extension of the analytic ansatz; the analysis of quenches starting from the topological ground state; the investigation of universal effects near the quantum critical point (μ=2t); the introduction of multiple interfaces or regions with spatially varying parameters to explore interference phenomena between quasiparticle fronts (e.g. periodic alternation (μ_1,μ_2,μ_1,μ_2,....) and, finally, the extension to interacting models where the Gaussian formalism is no longer applicable.