Abstract We investigate the possibility of characterizing the different thermalization pathways through a large-deviation approach. Specifically, we consider clean, disordered and quasi-periodic harmonic chains under energy and momentum-conserving noise. For their associated master equations, describing the dynamics of normal modes energies, we compute the fluctuations of activity and dynamical entropy in the corresponding biased ensembles. First-order dynamical phase transition are found that originates from different activity regions in action space. At the transitions, the steady-state in the biased ensembles changes from extended to localized, yielding a kind of condensation in normal-modes space. For the disordered and quasi-periodic models, we argue that the phase-diagram has a critical point at a finite value of the disorder or potential strength.
Abstract We investigate the possibility of characterizing the different thermalization pathways through a large-deviation approach. Specifically, we consider clean, disordered and quasi-periodic harmonic chains under energy and momentum-conserving noise. For their associated master equations, describing the dynamics of normal modes energies, we compute the fluctuations of activity and dynamical entropy in the corresponding biased ensembles. First-order dynamical phase transition are found that originates from different activity regions in action space. At the transitions, the steady-state in the biased ensembles changes from extended to localized, yielding a kind of condensation in normal-modes space. For the disordered and quasi-periodic models, we argue that the phase-diagram has a critical point at a finite value of the disorder or potential strength.
We investigate thermalization in a tight-binding chain with an on-site defect subject to local dephasing noise implemented as random phase kicks. For a single linear defect of strength ε, we obtain an exact analytical description of the system spectrum and formulate the dephasing-induced dynamics in the eigenstate basis. We derive an approximate kinetic equation for mode populations that describes a continuous-time random walk in action space. The walk transition rates are defined by the overlap matrix encoding the spatial structure of eigenstates that can be computed exactly. Analyzing the spectral properties of the equation, we show that defect-induced localized modes act as bottlenecks that strongly slow down relaxation, with rates scaling as ε^{-2} for strong defects. Using large-deviation theory, we characterize rare dynamical trajectories and identify distinct relaxation pathways associated with low- and high-activity regimes in action space. We provide numerical evidence that the large-deviation function exhibits a dynamical phase transition in the limit ε→∞. We then extend our analysis to the nonlinear case, considering a single nonlinear defect embedded in either a linear or a fully nonlinear discrete Schrödinger equation. Numerical simulations reveal a qualitatively faster approach to equilibrium than in the corresponding linear-defect model, driven by the amplitude-dependent weakening of the defect. Our results provide a unified framework for understanding thermalization, rare fluctuations, and relaxation pathways in stochastic tight-binding systems.
We consider a model of quantum dynamical semigroup on a finite dimensional fermionic space, obtained as the continuous-time limit of a repeated interactions model between a system and several thermal baths, with a dynamic driven by quadratic Hamiltonians. We assume that there is a globally conserved observable which can be expressed as a sum of energies on the system and on each baths, and we study the energy fluxes between the baths and the system. First, we consider only the mean energy fluxes, and prove that every thermal machines on quasi-free fermions in trivial, in the sense that it is not possible to extract energy from the coolest bath, even when we dispose of several other baths at different temperatures. Then, we consider an unraveling of the semigroup as a random process, and we study the large deviations of the energy fluxes, following Jaksic, Pillet and Westrich (2014). We reduce the computation of the cumulant generating functional to the resolution of a Riccati equation (which is formally similar to the study of large deviations in classical networks of harmonic oscillators). We apply it to the numerical computation of the rate function for energy exchanges in a fermionic chain, and show that larger fluctuations are observed on a longer chain.
We study the dynamics of the statistics of the energy transferred across a point along a quantum chain which is prepared in the inhomogeneous initial state obtained by joining two identical semi-infinite parts thermalized at two different temperatures. In particular, we consider the transverse field Ising and harmonic chains as prototypical models of noninteracting fermionic and bosonic excitations, respectively. Within the so-called hydrodynamic limit of large space-time scales we first discuss the mean values of the energy density and current, and then, aiming at the statistics of fluctuations, we calculate exactly the scaled cumulant generating function of the transferred energy. From the latter, the evolution of the associated large deviation function is obtained. A natural interpretation of our results is provided in terms of a semiclassical picture of quasiparticles moving ballistically along classical trajectories. Similarities and differences between the transferred energy scaled cumulant and the large deviation functions in the cases of noninteracting fermions and bosons are discussed.
Disordered elastic systems are studied in the quantum, classical, and dynamic regime in the present thesis. The thesis consists of two main parts: In the first part we study the low temperature phase diagram of one-dimensional weakly disordered quantum systems like charge or spin density waves and Luttinger liquids by a full finite temperature renormalization group (RG) calculation. In the classical region, for vanishing quantum fluctuations those results are supplemented by an exact solution of the model in the case of strong disorder, described by the ground state and the correlation function. Furthermore, by a mapping of the problem onto a Burgers equation with noise, in the case of weak disorder, we can derive an expression for the correlation length. At zero temperature we reproduce the (quantum) phase transition between a pinned (localized) and an unpinned (delocalized) phase for weak and strong quantum fluctuations, respectively, as found previously by Fukuyama or Giamarchi and Schulz. At finite temperatures the localization transition is suppressed: the random potential is wiped out by thermal fluctuations on length scales larger than the thermal de Broglie wave length of the phason excitations. The existence of a zero temperature transition is reflected in a rich cross-over phase diagram determined by the correlation functions. In particular we find four different scaling regions: a classical disordered, a quantum disordered, a quantum critical, and a thermal region. The results can be transferred directly to the discussion of the influence of disorder in superfluids. Finally we extend the RG calculation to the treatment of a commensurate lattice potential, and for the case of density waves we discuss the influence of quantum phase-slips, which might lead to a new scenario for the unpinning (delocalization) transition at zero temperature. Additionally, we analyze the current noise in a classical one-dimensional charge density wave system in the weak pinning regime by solving the overdamped equation of motion numerically. At low temperatures and just above the zero temperature depinning threshold, the power spectrum of the current noise reveals the existence of 1/f or flicker noise. Our results are in agreement with experimental measurements in quasi-one-dimensional charge density wave systems. In part two, we examine the viscous motion of an interface driven by a periodically oscillating external field in a random medium. The velocity exhibits a smeared depinning transition showing a pronounced hysteresis, which is absent in the adiabatic case. Its width is determined by a new length scale, introduced by the frequency of the external drive and - in the low frequency regime - by the critical exponents of the zero frequency depinning transition. Scaling arguments and an approximate renormalization group treatment predict a power law behavior of the velocity as a function of frequency and field amplitude at the dc depinning threshold, which is confirmed numerically. Thermal fluctuations lead to an additional smearing of the transition. If the amplitude of the external driving force is smaller than the depinning field, the motion of the interface by avalanches has to be taken into account in order to describe the velocity hysteresis. To check the validity of our model with respect to experimental results, we calculate the complex susceptibilities using an adiabatic and non-adiabatic approach numerically for finite systems and compare it to measurements of a superferromagnetic granular multilayer system and find that the model of an ac driven interface in a disordered environment can describe the main features of the experimental results. Additionally, the influence of a strong surface potential on the critical depinning of an elastic system, driven in a random medium is considered. If the surface potential prevents depinning completely, the curvature of the displacement profile at zero temperature exhibits a pronounced rhombic hysteresis curve with width of two times the (bulk) depinning field. The hysteresis disappears at non-zero temperatures if the driving force is changed adiabatically. If the surface depinns by the applied force or thermal creep, the curvature is reduced with increasing velocity. The results apply, e.g., to driven magnetic domain walls, flux-line lattices, and charge-density waves. As addendum we examine low dimensional interacting, but clean Bose systems at low temperatures. The interaction leads to scattering events of the over-condensate particles for which we calculate the scattering times using Fermi's golden rule. With that we derive the thermal conductivity and the related weak localization corrections. This is done for short and long-range interactions.
Most of the previous studies on classical thermalization focus on the wave-vector space, encountering limitations when extended beyond quasi-integrable regions. In this study, we propose a scheme to study the thermalization of the classical Hamiltonian chain of interacting oscillators in real space by developing a thermalization indicator proposed by Parisi [Europhys. Lett. 40, 357 (1997)0295-507510.1209/epl/i1997-00471-9], which approaches zero in the thermal state. Upon reaching the steady state characterized by the generalized Gibbs ensemble for a harmonic chain, a quench protocol is implemented to change the Hamiltonian to a nonintegrable form instantaneously, thereby preparing nonequilibrium initial states. This approach enables investigations of thermalization in real space, particularly valuable for exploring regions beyond quasi-integrability. For the FPUT-β lattice, we observe that the thermalization time as a function of the nonintegrable strength follows a -2 scaling law in the quasi-integrable region and -1/4 in the strongly integrable region. Moreover, numerical results reveal the thermalization time is proportional to the Lyapunov time, which bridges microscopic chaotic dynamics and the macroscopic thermalization process.
The work presented in this thesis is divided into two main parts: strong zero modes in quantum spin chains with open boundaries and large deviations in classical spin systems. The results concerning the strong zero modes are presented in chapters 2 and 3, while chapter 4 contains the results about large deviations in systems of many, but independent, classical degrees of freedom. These two distinct classes of problems in spin chains are both about aspects of slow relaxation to equilibrium in many-body systems. Nowadays slow dynamics and non-ergodicity in many-body quantum systems are receiving much attention due to their notable technological potential for storing and processing, for very long times, the quantum information encoded in boundary degrees of freedom. In the first part of the work presented in this thesis we mainly focus on quantum spin chains with open boundaries, which can display coherence times for edge spins, that diverge with the system size, as a consequence of the emergence, in the thermodynamic limit, of conserved local quantities, called strong zero modes. First, we discuss the fate of these coherence times when dynamics becomes dissipative. Understanding the behavior of these long correlation times in open quantum systems, undergoing Markovian (memoryless) dynamics, would allow for the possibility of exploiting the protected information encoded in the boundary spins for applications in quantum devices in realistic non-equilibrium settings. Second, we study one possible simple generalization of the transverse field Ising model to a class of ladder spin systems with plaquette interactions in the presence of a transverse field. Also for these systems one can explicitly construct the strong zero mode operators and understand the behavior of the consequential long coherence times for the edge degrees of freedom under isolated (unitary) dynamics. This is possible because the models we consider, while appearing complex superficially, can be brought to simple forms by exploiting an extensive number of conserved quantities they possess. The results presented in chapters 2 and 3 are closely related, even if they focus on different models and investigate different aspects of their dynamics. Nevertheless, the underlying concept they are based on is the same one, that is, an eigenstate phase transition for the strong zero modes that has interesting consequences for the full energy spectrum. In the second part of the thesis we focus on a different type of phase transitions. Here we address the issue of dynamical large deviation phase transitions from a disordered phase to an ordered one, in classical systems which consist in a collection of independent Ising spins evolving stochastically. Surprisingly, by reweighing trajectories accordingly to a certain type of time-integrated observable, even a system of independently evolving spins can display some critical behavior, at the fluctuation level, in trajectory space in the large size and long time limit. By making use of the tilted generator techniques we study the statistics of two particular time-integrated observables. Their associated tilted generators have an exact Kramers-Wannier duality and they look like, up to a sign and an additive constant, some specific well known quantum Hamiltonians. In this sense, the dynamical large deviation phase transitions we observe in these classical stochastic systems translate in ground state quantum phase transitions taking place at the self dual point of the tilted generators we study.