We investigate the non-equilibrium dynamics of the transverse field quantum Ising chain evolving from an inhomogeneous initial state given by joining two macroscopically different semi-infinite chains. We obtain integral expressions for all two-point correlation functions of the Jordan--Wigner Majorana fermions at any time and for any value of the transverse field. Using this result, we analytically compute the profiles of various physical observables in the space-time scaling limit and show that they can be obtained from a hydrodynamic picture based on ballistically propagating quasiparticles. Going beyond the hydrodynamic limit, we analyze the approach to the non-equilibrium steady state and find that the leading late time corrections display a lattice effect. We also study the fine structure of the propagating fronts which are found to be described by the Airy kernel and its derivatives. Near the front we observe the phenomenon of energy back-flow where the energy locally flows from the colder to the hotter region.
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.
Thermal rectification is the phenomenon by which the flux of heat depends on the direction of the flow. It has attracted much interest in recent years due to the possibility of devising thermal diodes. In this paper, we consider the rectification phenomenon in the quantum XXZ chain subject to an inhomogeneous field. The chain is driven out of equilibrium by the contact at its boundaries with two different reservoirs, leading to a constant flow of magnetization from one bath to the other. The nonunitary dynamics of this system, which is modeled by a Lindblad master equation, is treated exactly for small sizes and numerically for larger ones. The functional dependence of the rectification coefficient on the model parameters (anisotropy, field amplitude, and out of equilibrium driving strength) is investigated in full detail. Close to the XX point and at small inhomogeneity and low driving, we have found an explicit expression for the rectification coefficient that is valid at all system sizes. In particular, it shows that the phenomenon of rectification persists even in the thermodynamic limit. Finally, we prove that in the case of the XX chain, there is no rectification.
Results are presented for the nonequilibrium dynamics of a quantum XXZ -spin chain whose spins are initially arranged in a domain wall profile via the application of a magnetic field in the z direction, which is spatially varying along the chain. The system is driven out of equilibrium in two ways: a). by rapidly turning off the magnetic field, b). by rapidly quenching the interactions at the same time as the magnetic field is turned off. The time evolution of the domain wall profile as well as various two-point spin correlation functions are studied by the exact solution of the fermionic problem for the XX chain and via a bosonization approach and a mean-field approach for the XXZ chain. At long times the magnetization is found to equilibrate (reach the ground state value), while the two-point correlation functions in general do not. In particular, for quenches within the gapless XX phase, the spin correlation function transverse to the z direction acquires a spatially inhomogeneous structure at long times whose details depend on the initial domain wall profile. The spatial inhomogeneity is also recovered for the case of classical spins initially arranged in a domain wall profile and shows that the inhomogeneities arise due to the dephasing of transverse spin components as the domain wall broadens. A generalized Gibbs ensemble approach is found to be inadequate in capturing this spatially inhomogeneous state.
We investigate entanglement dynamics in bipartite systems governed by inhomogeneous Hamiltonians of the form
H=HL+HR, where
HL/R acts only on the left or right region and is homogeneous within each region. Focusing on the XX chain and the transverse-field Ising chain, we derive analytical formulas for the entanglement entropy between the two regions in the hydrodynamic limit of long times. In this regime, fermions incident on the interface undergo scattering, generating entanglement between reflected and transmitted modes. The resulting quasiparticle picture is controlled by the transmission coefficient, which we obtain analytically by solving the stationary lattice Schrödinger equation. Due to the bounded dispersion, strong inhomogeneity suppresses both transport and entanglement growth. We benchmark our analytical predictions against numerical simulations in paradigmatic setups. Finally, we extend the analysis to the interacting XXZ chain using tDMRG. The numerical data show qualitative agreement with the quadratic case: entanglement growth remains suppressed in the strongly inhomogeneous limit. Notably, however, entanglement continues to increase even when transport is suppressed, at least at intermediate times.
The equilibrium dynamics of the spin-
21 XX chain is reexamined within a recently developed formalism based on the quantum transfer matrix and a thermal form-factor expansion. The transversal correlation function is evaluated in real time and space. The high-accuracy calculation reproduces several exact results in limiting cases as well as the well-known asymptotic formulas obtained by the matrix Riemann-Hilbert approach. Furthermore, comparisons to numerical data based on a direct evaluation of the Pfaffian as well as to asymptotic formulas obtained within nonlinear Luttinger liquid theory are presented.