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

arXiv:2610.03133

Infrared Memory of a Majorana Shutter

Ali Vahedi

quant-phcond-mat.mes-hallSubmitted Oct 2, 2026• Announced Oct 5, 2026• 0 citations
arXivHTMLPDF

Abstract

We study a finite-width Majorana shutter as a local, time-dependent scattering defect and determine which features of its infrared memory are insensitive to the microscopic spatial profile. For a smooth shutter, the zero-frequency scattering matrix is fixed by the integrated coupling, while finite-width corrections vanish in the long-wavelength limit. A time-dependent switch then produces an infrared Bogoliubov kernel whose singular part is controlled by the change in this scattering matrix. We examine the many-body response through three observables. The overlap of the pre- and post-shutter Gaussian vacua exhibits Anderson orthogonality, with an exponent set by the scattering angle. A lattice calculation of the post-quench quasiparticle number resolves its logarithmic infrared growth, reaching the predicted coefficient in the slow-switch test and reproducing the expected dependence on the scattering-angle change. When the entanglement cut passes through the shutter, the entropy follows the effective-central-charge description of a conformal defect, with parameter-free agreement at the percent level, improving for the wider shutter. When the shutter lies inside the interval, its contribution remains a subleading finite-size correction. A local bilinear correlator approaches the expected Majorana power law and shows no stable logarithmic amplitude. These results show that the same low-frequency scattering data control distinct global and spatial observables, while their infrared coefficients are not interchangeable. The remaining limits are also made explicit: the finite-memory correction to the continuum kernel and the controlled continuum limit of large-interval spatial entropy.

Comments
39 pages and 5 figures

Related papers

A Truncated Majorana

Ali Vahedi

A Majorana field cannot be split into two spatial pieces, but a propagating Majorana can be truncated by a time-dependent shutter. We study this truncation for the two chiral components of a massless Majorana field in one spatial dimension, modeling the shutter as a smooth time-dependent $SO(2)$ rotation of the two components. Although orthogonal at the level of the real field algebra, this rotation does not preserve the positive- and negative-frequency splitting that defines the Fock vacuum, so the in-vacuum becomes a fermionic squeezed state in outgoing modes. At any finite infrared regulator, a single incident Majorana is then a single outgoing quasiparticle on top of this squeezed vacuum. The squeezed background carries a universal soft response: its anomalous Bogoliubov kernel falls as the inverse sum of the two frequencies at low energy, so the anomalous Hilbert--Schmidt norm grows logarithmically with the infrared cutoff, with a coefficient fixed by the asymptotic shutter mismatch. The vacuum overlap exhibits a Fermi-edge type orthogonality exponent which, for the rotation-angle lift $Δθ$ with $|Δθ|<π$, is exactly $(Δθ/2π)^2$, in analogy with the Fermi-edge singularity and the Kondo and Luttinger-liquid problems. We prove local equivalence non-perturbatively for the squeezed-vacuum factor under a compact-support hypothesis, and give an exact local pullback identity for the full one-quasiparticle state on the even CAR algebra; combined with wavepacket-localization assumptions, these yield the expected local vacuum and one-Majorana descriptions. We verify profile independence using a smooth switching family and contrast the temporal benchmark with a number-conserving beam-splitter benchmark whose entropy is bounded by $\log 2$.

hep-thCross-listedarXiv:2608.30723•Sep 1, 2026•0 citations

The Junction That Remembers: Infrared Response of Two Coupled SYK Majorana Systems

Ali Vahedi

Two strongly interacting Majorana systems of the Sachdev--Ye--Kitaev type, each with its own independent disorder, are joined by a link that is switched on in time. We ask how the infrared structure of the link is imprinted on the response and whether the system remembers how it was opened. At large $N$ the leading cross-correlation between the two systems is an exact product of two equilibrium propagators and the switching profile; because the disorder is independent, no ladder dresses it. The link operator has scaling dimension one half, which gives its susceptibility a universal logarithm whose coefficient and phase we fix and verify numerically, while the gap of the coupled equilibrium problem grows more slowly than the naive crossover estimate. Exact dynamics for up to twelve Majoranas per side show strong parity transfer and entanglement several times the free-fermion ceiling. The work done by the switch has a first-order part that does not depend on the switching time and a second-order part that grows in magnitude with it. Once the work is held fixed, the protocol dependence of the entanglement, the link expectation value and a scrambling diagnostic nearly vanishes, and the small residuals do not depend on whether the protocols differ at low or at high frequency. These small systems have no infrared window, so they neither confirm nor exclude memory at large $N$.

arXiv:2609.05336•Sep 4, 2026•0 citations

Orthogonality Edges in Strong-Coupling Quantum Work Statistics

Atta ur Rahman, Muhammad Noman, S. M. Zangi +1 more

Strong coupling to a reservoir can do more than shift, broaden, or dress the work peaks of a driven quantum system. When the reservoir is infrared singular, a sudden change of a local control parameter can alter the boundary condition seen by infinitely many low-energy modes, converting a quasiparticle-like threshold line into a many-body edge. We demonstrate this mechanism for the inclusive work distribution of the biased spin-boson model under a sudden bias inversion. In the independent-boson limit, the problem is exactly solvable and gives a sharp infrared classification: a super-Ohmic bath can retain a finite elastic threshold weight, whereas Ohmic and sub-Ohmic baths extinguish the elastic line through boundary orthogonality. At the Ohmic fixed point, the same exponent controls both the vanishing elastic residue and the low-work continuum. We then ask how this edge is resolved away from the static-boundary limit. Using displaced-basis exact diagonalization of logarithmically discretized baths, we find that finite tunnelling leaves an edge-like continuum over the accessible energy window, while separating two operational diagnostics of the threshold: the cumulative-continuum exponent extracted from $z$-interleaved spectra lies above the elastic-overlap exponent extracted from $z$-averaged overlaps, $θ_C>θ_Z$. We interpret this separation as a finite-energy crossover away from the static-boundary fixed point, not as evidence for a new asymptotic fixed point. The separation survives fitting-window variation, oscillator-cutoff checks, spectrum-size checks, and leave-one-$z$-out tests, while time-domain characteristic functions provide a compatible but non-decisive diagnostic. Finally, the same threshold edge controls the sampling cost of Jarzynski-type exponential averages, making rare low-work events increasingly important at low temperature.

quant-pharXiv:2607.03950•Jul 7, 2026•0 citations

Memory-Dominated Quantum Criticality as a Universal Route to High-Temperature Superconductivity

Byung Gyu Chae

Understanding the dynamical origin of high-temperature superconductivity remains a central challenge in strongly correlated quantum matter. Near quantum criticality, diverging correlation times reorganize the infrared dynamics into a scale-free continuum of collective relaxation processes. We show that the infrared behavior of interacting electrons is generically controlled by the relaxation-rate spectrum of the underlying many-body dynamics. Starting from a microscopic fermionic theory, we derive that the Cooper-channel kernel admits a universal spectral representation in terms of the time-scale density of states (TDOS) of collective decay modes, without invoking a specific bosonic mediator. The superconducting instability follows directly from the vanishing of the quadratic kernel via a standard ladder resummation and Thouless criterion, with the pairing interaction determined entirely by the infrared structure of the relaxation spectrum. A finite TDOS at vanishing relaxation rate produces a memory-dominated regime characterized by long-time kernels $K(t)\sim 1/t$ and logarithmic enhancement of the retarded pairing interaction, leading to a BCS-like exponential transition scale set by infrared spectral weight. More generally, infrared-singular spectra generate power-law response and algebraic enhancement of the transition scale. The same relaxation spectrum controls normal-state dynamics, giving rise to long-time correlations, non-Markovian response, and strange-metal behavior. These results identify the spectral organization of relaxation modes as a universal organizing principle of quantum critical matter and establish memory-dominated criticality as a natural mechanism for enhanced pairing.

arXiv:2602.22626•Feb 26, 2026•0 citations

Entanglement entropy across the dynamical phase transition in the quantum $\mathcal{O}(N)$ model

Frederick del Pozo, Tangi Morvan, Irénée Frérot +1 more

We demonstrate that the dynamical phase transition of the quantum $\mathcal{O}(N)$ model at large $N$ leaves universal fingerprints in the infrared structure of the entanglement spectrum. While the leading contribution to the entanglement entropy at long time follows the conventional volume law associated with ballistic entanglement spreading, its subleading behavior sharply distinguishes the different dynamical regimes. Specifically, quenches at and below the critical point generate gapless low-energy entanglement modes together with logarithmic corrections to the long-time entanglement entropy, whose scaling is governed by the scaling exponent of the transition. Using an infinite-slab bipartition geometry and exact numerical correlation functions in the large-$N$ limit, we characterize these scaling laws across the dynamical phase diagram and relate them to the emergence of long-range correlations during the post-quench dynamics. We further show that the entanglement eigenmodes themselves reveal characteristic signatures of the dynamical phase transition through their spatio-temporal structure and degeneracy properties.

quant-pharXiv:2605.23600•May 25, 2026•0 citations

Entanglement dynamics after spatially inhomogeneous quenches.

ANDREA DI PASQUALE

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 \(t_{\min} \to \infty\). 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.

•Dec 4, 2025•0 citations
We study a finite-width Majorana shutter as a local, time-dependent scattering defect and determine which features of its infrared memory are insensitive to the microscopic spatial profile. For a smooth shutter, the zero-frequency scattering matrix is fixed by the integrated coupling, while finite-width corrections vanish in the long-wavelength limit. A time-dependent switch then produces an infrared Bogoliubov kernel whose singular part is controlled by the change in this scattering matrix.
We examine the many-body response through three observables. The overlap of the pre- and post-shutter Gaussian vacua exhibits Anderson orthogonality, with an exponent set by the scattering angle. A lattice calculation of the post-quench quasiparticle number resolves its logarithmic infrared growth, reaching the predicted coefficient in the slow-switch test and reproducing the expected dependence on the scattering-angle change. When the entanglement cut passes through the shutter, the entropy follows the effective-central-charge description of a conformal defect, with parameter-free agreement at the percent level, improving for the wider shutter. When the shutter lies inside the interval, its contribution remains a subleading finite-size correction. A local bilinear correlator approaches the expected Majorana power law and shows no stable logarithmic amplitude.
These results show that the same low-frequency scattering data control distinct global and spatial observables, while their infrared coefficients are not interchangeable. The remaining limits are also made explicit: the finite-memory correction to the continuum kernel and the controlled continuum limit of large-interval spatial entropy.
A Majorana field cannot be split into two spatial pieces, but a propagating Majorana can be truncated by a time-dependent shutter. We study this truncation for the two chiral components of a massless Majorana field in one spatial dimension, modeling the shutter as a smooth time-dependent SO(2)SO(2)SO(2) rotation of the two components. Although orthogonal at the level of the real field algebra, this rotation does not preserve the positive- and negative-frequency splitting that defines the Fock vacuum, so the in-vacuum becomes a fermionic squeezed state in outgoing modes. At any finite infrared regulator, a single incident Majorana is then a single outgoing quasiparticle on top of this squeezed vacuum. The squeezed background carries a universal soft response: its anomalous Bogoliubov kernel falls as the inverse sum of the two frequencies at low energy, so the anomalous Hilbert--Schmidt norm grows logarithmically with the infrared cutoff, with a coefficient fixed by the asymptotic shutter mismatch. The vacuum overlap exhibits a Fermi-edge type orthogonality exponent which, for the rotation-angle lift ΔθΔθΔθ with ∣Δθ∣<π|Δθ|<π∣Δθ∣<π, is exactly (Δθ/2π)2(Δθ/2π)^2(Δθ/2π)2, in analogy with the Fermi-edge singularity and the Kondo and Luttinger-liquid problems. We prove local equivalence non-perturbatively for the squeezed-vacuum factor under a compact-support hypothesis, and give an exact local pullback identity for the full one-quasiparticle state on the even CAR algebra; combined with wavepacket-localization assumptions, these yield the expected local vacuum and one-Majorana descriptions. We verify profile independence using a smooth switching family and contrast the temporal benchmark with a number-conserving beam-splitter benchmark whose entropy is bounded by log⁡2\log 2log2.
Two strongly interacting Majorana systems of the Sachdev--Ye--Kitaev type, each with its own independent disorder, are joined by a link that is switched on in time. We ask how the infrared structure of the link is imprinted on the response and whether the system remembers how it was opened. At large NNN the leading cross-correlation between the two systems is an exact product of two equilibrium propagators and the switching profile; because the disorder is independent, no ladder dresses it. The link operator has scaling dimension one half, which gives its susceptibility a universal logarithm whose coefficient and phase we fix and verify numerically, while the gap of the coupled equilibrium problem grows more slowly than the naive crossover estimate. Exact dynamics for up to twelve Majoranas per side show strong parity transfer and entanglement several times the free-fermion ceiling. The work done by the switch has a first-order part that does not depend on the switching time and a second-order part that grows in magnitude with it. Once the work is held fixed, the protocol dependence of the entanglement, the link expectation value and a scrambling diagnostic nearly vanishes, and the small residuals do not depend on whether the protocols differ at low or at high frequency. These small systems have no infrared window, so they neither confirm nor exclude memory at large NNN.
Strong coupling to a reservoir can do more than shift, broaden, or dress the work peaks of a driven quantum system. When the reservoir is infrared singular, a sudden change of a local control parameter can alter the boundary condition seen by infinitely many low-energy modes, converting a quasiparticle-like threshold line into a many-body edge. We demonstrate this mechanism for the inclusive work distribution of the biased spin-boson model under a sudden bias inversion. In the independent-boson limit, the problem is exactly solvable and gives a sharp infrared classification: a super-Ohmic bath can retain a finite elastic threshold weight, whereas Ohmic and sub-Ohmic baths extinguish the elastic line through boundary orthogonality. At the Ohmic fixed point, the same exponent controls both the vanishing elastic residue and the low-work continuum. We then ask how this edge is resolved away from the static-boundary limit. Using displaced-basis exact diagonalization of logarithmically discretized baths, we find that finite tunnelling leaves an edge-like continuum over the accessible energy window, while separating two operational diagnostics of the threshold: the cumulative-continuum exponent extracted from zzz-interleaved spectra lies above the elastic-overlap exponent extracted from zzz-averaged overlaps, θC>θZθ_C>θ_ZθC​>θZ​. We interpret this separation as a finite-energy crossover away from the static-boundary fixed point, not as evidence for a new asymptotic fixed point. The separation survives fitting-window variation, oscillator-cutoff checks, spectrum-size checks, and leave-one-zzz-out tests, while time-domain characteristic functions provide a compatible but non-decisive diagnostic. Finally, the same threshold edge controls the sampling cost of Jarzynski-type exponential averages, making rare low-work events increasingly important at low temperature.
Understanding the dynamical origin of high-temperature superconductivity remains a central challenge in strongly correlated quantum matter. Near quantum criticality, diverging correlation times reorganize the infrared dynamics into a scale-free continuum of collective relaxation processes. We show that the infrared behavior of interacting electrons is generically controlled by the relaxation-rate spectrum of the underlying many-body dynamics. Starting from a microscopic fermionic theory, we derive that the Cooper-channel kernel admits a universal spectral representation in terms of the time-scale density of states (TDOS) of collective decay modes, without invoking a specific bosonic mediator. The superconducting instability follows directly from the vanishing of the quadratic kernel via a standard ladder resummation and Thouless criterion, with the pairing interaction determined entirely by the infrared structure of the relaxation spectrum. A finite TDOS at vanishing relaxation rate produces a memory-dominated regime characterized by long-time kernels K(t)∼1/tK(t)\sim 1/tK(t)∼1/t and logarithmic enhancement of the retarded pairing interaction, leading to a BCS-like exponential transition scale set by infrared spectral weight. More generally, infrared-singular spectra generate power-law response and algebraic enhancement of the transition scale. The same relaxation spectrum controls normal-state dynamics, giving rise to long-time correlations, non-Markovian response, and strange-metal behavior. These results identify the spectral organization of relaxation modes as a universal organizing principle of quantum critical matter and establish memory-dominated criticality as a natural mechanism for enhanced pairing.
Entanglement entropy across the dynamical phase transition in the quantum O(N)\mathcal{O}(N)O(N) model
We demonstrate that the dynamical phase transition of the quantum O(N)\mathcal{O}(N)O(N) model at large NNN leaves universal fingerprints in the infrared structure of the entanglement spectrum. While the leading contribution to the entanglement entropy at long time follows the conventional volume law associated with ballistic entanglement spreading, its subleading behavior sharply distinguishes the different dynamical regimes. Specifically, quenches at and below the critical point generate gapless low-energy entanglement modes together with logarithmic corrections to the long-time entanglement entropy, whose scaling is governed by the scaling exponent of the transition. Using an infinite-slab bipartition geometry and exact numerical correlation functions in the large-NNN limit, we characterize these scaling laws across the dynamical phase diagram and relate them to the emergence of long-range correlations during the post-quench dynamics. We further show that the entanglement eigenmodes themselves reveal characteristic signatures of the dynamical phase transition through their spatio-temporal structure and degeneracy properties.
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⁡→∞t_{\min} \to \inftytmin​→∞. 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.