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:2605.01595 · Journal of Statistical Mechanics Theory and Experiment

Entanglement dynamics after quenches with inhomogeneous Hamiltonians

Andrea Di Pasquale, Federico Rottoli, Vincenzo Alba

cond-mat.stat-mechcond-mat.quant-gascond-mat.str-elhep-thquant-phSubmitted May 2, 2026• Announced May 5, 2026• 0 citations• Revised Oct 13, 2026 (v2)
v1 · May 2, 2026v2 · Oct 13, 2026
arXivHTMLPDF

Abstract

We investigate entanglement dynamics in bipartite systems governed by inhomogeneous Hamiltonians of the form $H = H_L + H_R$, where $H_{L/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.

AI Summary

The paper studies entanglement growth across an interface joining two spin-chain regions with different uniform magnetic fields. It considers quenches in the free-fermion XX and transverse-field Ising chains, then explores interactions numerically in the XXZ chain. A scattering-based quasiparticle picture conjectures the hydrodynamic entropy growth rate. When only the left region produces excitations, it gives $S(t)=t\int_0^\pi \frac{dk}{2\pi}|v_L(k)|s_{YY}(T(k)n_L(k))$, where $T$ is the interface transmission probability and $s_{YY}$ is binary entropy. Transmission coefficients are obtained analytically from stationary lattice scattering equations. Bounded dispersion restricts transmission to overlapping energy bands; in the XX chain the linear growth rate vanishes for $|h_L-h_R|\geq 2J$. Gaussian correlation-matrix simulations, including chains of 1,000 sites, support the one-sided formula and an additive two-sided extension for the tested Ising quenches. That extension fails for XX quenches with nontrivial occupations on both sides; correlations beyond occupation densities remain unresolved. TEBD simulations of XXZ chains at $\Delta=0.7$ and $2$ show intermediate-time entanglement growth beyond the free-chain transmission threshold. At $\Delta=2$ and field difference $3$, growth persists through the converged window $t\lesssim15$, while magnetization and energy transport appear suppressed. Larger field differences yield apparent saturation. Whether the growth survives asymptotically is unresolved: later-time calculations remain bond-dimension dependent even at $\chi=2500$. The entropy formulas remain conjectural, assume pair-producing initial states, and apply to a bipartition located at the interface.

Little sign of AI-generated textIncrementalStrong assumptionsOverclaims

Central claim may not hold: The abstract describes the entropy formulas as derived, whereas the body presents them as numerically tested conjectures. The paper's own two-sided XX-chain Néel-state benchmark contradicts its proposed two-sided entropy formula, limiting the scope of the abstract's analytical claim. The entropy construction assumes initial states generating pairwise-entangled quasiparticles and does not establish the abstract's apparent broader initial-state scope. The printed Ising group velocity is twice the derivative of the printed dispersion, creating a factor-of-two inconsistency in the stated entropy prediction. The Ising transmission formula prints the velocity ratio backwards, although its final closed-form transmission coefficient uses the correct flux ratio.

Written by gpt-6.1-sol from the paper's full text. Automated: it can be wrong.

Journal reference
J. Stat. Mech. (2026) 103101
DOI
10.1088/1742-5468/aea554
Comments
20 pages, 9 figures

Related papers

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

Entanglement dynamics in confining spin chains

Stefano Scopa, Pasquale Calabrese, Alvise Bastianello

The confinement of elementary excitations induces distinctive features in the non-equilibrium quench dynamics. One of the most remarkable is the suppression of entanglement entropy, which in several instances turns out to oscillate rather than grow indefinitely. While the qualitative physical origin of this behavior is clear, till now no quantitative understanding away from the field theory limit was available. Here we investigate this problem in the weak quench limit, when mesons are excited at rest, hindering entropy growth and exhibiting persistent oscillations. We provide analytical predictions of the entire entanglement dynamics based on a Gaussian approximation of the many-body state, which captures numerical data with great accuracy and is further simplified to a semiclassical quasiparticle picture in the regime of weak confinement. Our methods are valid in general and we apply explicitly to two prototypical models: the Ising chain in a tilted field and the experimentally relevant long-range Ising model.

arXiv:2111.11483•Physical review. B./Physical review. B•Mar 16, 2022•29 citations

Spreading of entanglement and correlations after a quench with intertwined quasiparticles

Alvise Bastianello, Pasquale Calabrese

We extend the semiclassical picture for the spreading of entanglement and correlations to quantum quenches with several species of quasiparticles that have non-trivial pair correlations in momentum space. These pair correlations are, for example, relevant in inhomogeneous lattice models with a periodically-modulated Hamiltonian parameter. We provide explicit predictions for the spreading of the entanglement entropy in the space-time scaling limit. We also predict the time evolution of one- and two-point functions of the order parameter for quenches within the ordered phase. We test all our predictions against exact numerical results for quenches in the Ising chain with a modulated transverse field and we find perfect agreement.

arXiv:1807.10176•SciPost Physics•Oct 15, 2018•74 citations

Entanglement Oscillations near a Quantum Critical Point

Olalla A. Castro-Alvaredo, Máté Lencsés, István M. Szécsényi +1 more

We study the dynamics of entanglement in the scaling limit of the Ising spin chain in the presence of both a longitudinal and a transverse field. We present analytical results for the quench of the longitudinal field in the critical transverse field which go beyond current lattice integrability techniques. We test these results against a numerical simulation on the corresponding lattice model finding extremely good agreement. We show that the presence of bound states in the spectrum of the field theory leads to oscillations in the entanglement entropy and suppresses its linear growth on the time scales accessible to numerical simulations. For small quenches, we exactly determine these oscillatory contributions and demonstrate that their presence follows from symmetry arguments. For the quench of the transverse field at zero longitudinal field, we prove that the Rényi entropies are exactly proportional to the logarithm of the exponential of a time-dependent function, whose leading large-time behavior is linear, hence, entanglement grows linearly. We conclude that, in the scaling limit, linear growth and oscillations in the entanglement entropies can not be simply seen as consequences of integrability and its breaking, respectively.

arXiv:2001.10007•Physical Review Letters•Jun 11, 2020•75 citations

Entanglement growth after inhomogenous quenches

Tibor Rakovszky, Curt W. von Keyserlingk, Frank Pollmann

We study the growth of entanglement in quantum systems with a conserved quantity exhibiting diffusive transport, focusing on how initial inhomogeneities are imprinted on the entropy. We propose a simple effective model, which generalizes the minimal cut picture of Jonay, Huse, and Nahum [arXiv:1803.00089] in such a way that the ``line tension'' of the cut depends on the local entropy density. In the case of noisy dynamics, this is described by the Kardar-Parisi-Zhang (KPZ) equation coupled to a diffusing field. We investigate the resulting dynamics and find that initial inhomogeneities of the conserved charge give rise to features in the entanglement profile, whose width and height both grow in time as $\ensuremath{\propto}\sqrt{t}$. In particular, for a domain wall quench, diffusion restricts entanglement growth to be ${S}_{\text{vN}}\ensuremath{\lesssim}\sqrt{t}$. We find that for charge density wave initial states, these features in the entanglement profile are present even after the charge density has equilibrated. Our conclusions are supported by numerical results on random circuits and deterministic spin chains.

arXiv:1907.00869•Physical review. B./Physical review. B•Sep 18, 2019•28 citations

Molecular dynamics simulation of entanglement spreading in generalized hydrodynamics

Márton Mestyán, Vincenzo Alba

We consider a molecular dynamics method, the so-called flea gas for computing the evolution of entanglement after inhomogeneous quantum quenches in an integrable quantum system. In such systems the evolution of local observables is described at large space-time scales by the Generalized Hydrodynamics approach, which is based on the presence of stable, ballistically propagating quasiparticles. Recently it was shown that the GHD approach can be joined with the quasiparticle picture of entanglement evolution, providing results for entanglement growth after inhomogeneous quenches. Here we apply the flea gas simulation of GHD to obtain numerical results for entanglement growth. We implement the flea gas dynamics for the gapped anisotropic Heisenberg XXZ spin chain, considering quenches from globally homogeneous and piecewise homogeneous initial states. While the flea gas method applied to the XXZ chain is not exact even in the scaling limit (in contrast to the Lieb-Liniger model), it yields a very good approximation of analytical results for entanglement growth in the cases considered. Furthermore, we obtain the full-time dynamics of the mutual information after quenches from inhomogeneous settings, for which no analytical results are available.

arXiv:1905.03206•SciPost Physics•Apr 9, 2020•31 citations
We investigate entanglement dynamics in bipartite systems governed by inhomogeneous Hamiltonians of the form H=HL+HRH = H_L + H_RH=HL​+HR​, where HL/RH_{L/R}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 paper studies entanglement growth across an interface joining two spin-chain regions with different uniform magnetic fields. It considers quenches in the free-fermion XX and transverse-field Ising chains, then explores interactions numerically in the XXZ chain.
A scattering-based quasiparticle picture conjectures the hydrodynamic entropy growth rate. When only the left region produces excitations, it gives S(t)=t∫0πdk2π∣vL(k)∣sYY(T(k)nL(k))S(t)=t\int_0^\pi \frac{dk}{2\pi}|v_L(k)|s_{YY}(T(k)n_L(k))S(t)=t∫0π​2πdk​∣vL​(k)∣sYY​(T(k)nL​(k)), where TTT is the interface transmission probability and sYYs_{YY}sYY​ is binary entropy. Transmission coefficients are obtained analytically from stationary lattice scattering equations. Bounded dispersion restricts transmission to overlapping energy bands; in the XX chain the linear growth rate vanishes for ∣hL−hR∣≥2J|h_L-h_R|\geq 2J∣hL​−hR​∣≥2J. Gaussian correlation-matrix simulations, including chains of 1,000 sites, support the one-sided formula and an additive two-sided extension for the tested Ising quenches. That extension fails for XX quenches with nontrivial occupations on both sides; correlations beyond occupation densities remain unresolved.
TEBD simulations of XXZ chains at Δ=0.7\Delta=0.7Δ=0.7 and 222 show intermediate-time entanglement growth beyond the free-chain transmission threshold. At Δ=2\Delta=2Δ=2 and field difference 333, growth persists through the converged window t≲15t\lesssim15t≲15, while magnetization and energy transport appear suppressed. Larger field differences yield apparent saturation. Whether the growth survives asymptotically is unresolved: later-time calculations remain bond-dimension dependent even at χ=2500\chi=2500χ=2500. The entropy formulas remain conjectural, assume pair-producing initial states, and apply to a bipartition located at the interface.
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.
The confinement of elementary excitations induces distinctive features in the non-equilibrium quench dynamics. One of the most remarkable is the suppression of entanglement entropy, which in several instances turns out to oscillate rather than grow indefinitely. While the qualitative physical origin of this behavior is clear, till now no quantitative understanding away from the field theory limit was available. Here we investigate this problem in the weak quench limit, when mesons are excited at rest, hindering entropy growth and exhibiting persistent oscillations. We provide analytical predictions of the entire entanglement dynamics based on a Gaussian approximation of the many-body state, which captures numerical data with great accuracy and is further simplified to a semiclassical quasiparticle picture in the regime of weak confinement. Our methods are valid in general and we apply explicitly to two prototypical models: the Ising chain in a tilted field and the experimentally relevant long-range Ising model.
We extend the semiclassical picture for the spreading of entanglement and correlations to quantum quenches with several species of quasiparticles that have non-trivial pair correlations in momentum space. These pair correlations are, for example, relevant in inhomogeneous lattice models with a periodically-modulated Hamiltonian parameter. We provide explicit predictions for the spreading of the entanglement entropy in the space-time scaling limit. We also predict the time evolution of one- and two-point functions of the order parameter for quenches within the ordered phase. We test all our predictions against exact numerical results for quenches in the Ising chain with a modulated transverse field and we find perfect agreement.
We study the dynamics of entanglement in the scaling limit of the Ising spin chain in the presence of both a longitudinal and a transverse field. We present analytical results for the quench of the longitudinal field in the critical transverse field which go beyond current lattice integrability techniques. We test these results against a numerical simulation on the corresponding lattice model finding extremely good agreement. We show that the presence of bound states in the spectrum of the field theory leads to oscillations in the entanglement entropy and suppresses its linear growth on the time scales accessible to numerical simulations. For small quenches, we exactly determine these oscillatory contributions and demonstrate that their presence follows from symmetry arguments. For the quench of the transverse field at zero longitudinal field, we prove that the Rényi entropies are exactly proportional to the logarithm of the exponential of a time-dependent function, whose leading large-time behavior is linear, hence, entanglement grows linearly. We conclude that, in the scaling limit, linear growth and oscillations in the entanglement entropies can not be simply seen as consequences of integrability and its breaking, respectively.
We study the growth of entanglement in quantum systems with a conserved quantity exhibiting diffusive transport, focusing on how initial inhomogeneities are imprinted on the entropy. We propose a simple effective model, which generalizes the minimal cut picture of Jonay, Huse, and Nahum [arXiv:1803.00089] in such a way that the ``line tension'' of the cut depends on the local entropy density. In the case of noisy dynamics, this is described by the Kardar-Parisi-Zhang (KPZ) equation coupled to a diffusing field. We investigate the resulting dynamics and find that initial inhomogeneities of the conserved charge give rise to features in the entanglement profile, whose width and height both grow in time as $\ensuremath{\propto}\sqrt{t}$. In particular, for a domain wall quench, diffusion restricts entanglement growth to be ${S}_{\text{vN}}\ensuremath{\lesssim}\sqrt{t}$. We find that for charge density wave initial states, these features in the entanglement profile are present even after the charge density has equilibrated. Our conclusions are supported by numerical results on random circuits and deterministic spin chains.
We consider a molecular dynamics method, the so-called flea gas for computing the evolution of entanglement after inhomogeneous quantum quenches in an integrable quantum system. In such systems the evolution of local observables is described at large space-time scales by the Generalized Hydrodynamics approach, which is based on the presence of stable, ballistically propagating quasiparticles. Recently it was shown that the GHD approach can be joined with the quasiparticle picture of entanglement evolution, providing results for entanglement growth after inhomogeneous quenches. Here we apply the flea gas simulation of GHD to obtain numerical results for entanglement growth. We implement the flea gas dynamics for the gapped anisotropic Heisenberg XXZ spin chain, considering quenches from globally homogeneous and piecewise homogeneous initial states. While the flea gas method applied to the XXZ chain is not exact even in the scaling limit (in contrast to the Lieb-Liniger model), it yields a very good approximation of analytical results for entanglement growth in the cases considered. Furthermore, we obtain the full-time dynamics of the mutual information after quenches from inhomogeneous settings, for which no analytical results are available.