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:2403.18906 · Physical review. B./Physical review. B

Einstein relation for subdiffusive relaxation in Stark chains

Peter Prelovšek, Sourav Nandy, Marcin Mierzejewski

Published Aug 7, 2024• 1 citation
arXivHTMLPDF

Abstract

We investigate chains of interacting spinless fermions subject to a finite external field $F$ (also called Stark chains) and focus on the regime where the charge thermalization follows the subdiffusive hydrodynamics. First, we study reduced models conserving the dipole moment and derive an explicit Einstein relation which links the subdiffusive transport coefficient with the correlations of the dipolar current. This relation explains why the decay rate ${\mathrm{\ensuremath{\Gamma}}}_{q}$ of the density modulation with wave vector $q$ shows ${q}^{4}$ dependence. In the case of the Stark model, a similar Einstein relation is also derived and tested using various numerical methods. They confirm an exponential reduction of the transport coefficient with increasing $F$. On the other hand, our study of the Stark model indicates that upon increasing $q$ there is a crossover from subdiffusive behavior, ${\mathrm{\ensuremath{\Gamma}}}_{q}\ensuremath{\propto}{q}^{4}$, to the normal diffusive relaxation, ${\mathrm{\ensuremath{\Gamma}}}_{q}\ensuremath{\propto}{q}^{2}$, at the wave vector ${q}^{*}$ which vanishes for $F\ensuremath{\rightarrow}0$.

DOI
10.1103/physrevb.110.l081105

Related papers

Emergent dipole moment conservation and subdiffusion in tilted chains

Sourav Nandy, J. Herbrych, Zala Lenarčič +3 more

We study transport in an interacting tilted (Stark) chain. We show that the crossover between diffusive and subdiffusive relaxation is governed by $F\sqrt{L}$, where $F$ is the strength of the field, and $L$ is the wavelength of the excitation. While the subdiffusive behavior persists for large fields, the corresponding transport coefficient is exponentially suppressed with $F$ so that the finite-time dynamics appears almost frozen. We (i) explain the crossover scale between the diffusive and subdiffusive transport by bounding the dynamics of the dipole moment for arbitrary initial state, (ii) prove its emergent conservation at infinite temperature for $F\ensuremath{\gg}1/L$, and (iii) argue that the numerical data for the tilted chain are consistent with the hydrodynamics of fractons.

arXiv:2310.01862•Physical review. B./Physical review. B•Mar 12, 2024•12 citations

Spin subdiffusion in perturbed infinite-U Hubbard chain

Jakub Rękas, Marcin Mierzejewski, Zala Lenarčič +1 more

The $t$-model represents the Hubbard model in the limit $U \to \infty$ and is one of the basic models of strongly correlated electrons. On a one-dimensional chain, the model is integrable, and the charge dynamics corresponds to that of free spinless fermions. However, the sequence of spins is frozen, leading to the Hilbert space fragmentation and nontrivial spin dynamics. We consider integrable and perturbed models with perturbations that break integrability while preserving fragmentation, and show that they exhibit various types of spin dynamics, from ballistic transport to anomalous diffusion in the integrable case, and from diffusion to subdiffusion in the perturbed case. Due to fragmentation, in all cases considered, spin transport is mediated by charge transport, with a particular magnetization dependence, most notably leading to subdiffusion in the grandcanonical average of the perturbed model, with a mechanism distinct from subdiffusion in disordered or dipole-conserving models.

cond-mat.str-elCross-listedarXiv:2603.19743•Physical review. B./Physical review. B•Mar 23, 2026•0 citations

Density propagator for many-body localization: finite size effects, transient subdiffusion, (stretched-)exponentials

Soumya Bera, Giuseppe De Tomasi, Felix Weiner +1 more

We investigate charge relaxation in the spin-less disordered fermionic Hubbard chain ($t{-}V$-model). Our observable is the time-dependent density propagator, $\Pi_{\varepsilon}(x,t)$, calculated in windows of different energy density, $\varepsilon$, of the many-body Hamiltonian and at different disorder strengths, $W$, not exceeding the critical value $W_\text{c}$. The width $\Delta x_\varepsilon(t)$ of $\Pi_\varepsilon(x,t)$ exhibits a behavior $d\ln \Delta x_\varepsilon(t) / d\ln t {=} \beta_\varepsilon(t)$, where the exponent function $\beta_\varepsilon(t){\lesssim}1/2$ is seen to depend strongly on $L$ at all investigated parameter combinations. (i) We cannot confirm the existence of a region in phase space that exhibits (genuine) subdiffusive dynamics in the sense that $\beta_\varepsilon{<}1/2$ is numerically fixed in the limit of large $L$. Instead, subdiffusion might possibly be transient, only, finally giving way to conventional diffusive behavior with $\beta_\varepsilon{=}1/2$. (ii) Similarly, we cannot confirm the existence of many-body mobility edges deep in the delocalized phase. (iii) (Transient) subdiffusion $0<\beta_\varepsilon(t)\lesssim 1/2$, coexists with an enhanced probability for returning to the origin, $\Pi_\varepsilon(0,t)$, decaying much slower than $1/\Delta x_\varepsilon (t)$. Correspondingly, the spatial decay of $\Pi_\varepsilon(x,t)$ is far from Gaussian being exponential or even slower. On a phenomenological level, our findings are broadly consistent with effects of strong disorder and (fractal) Griffiths regions.

•Oct 10, 2016•4 citations

Second-order hydrodynamics for fermionic cold atoms: Detailed analysis of transport coefficients and relaxation times

Yuta Kikuchi, Kyosuke Tsumura, Teiji Kunihiro

We give a detailed derivation of the second-order (local) hydrodynamics for Boltzmann equation with an external force by using the renormalization group method. In this method, we solve the Boltzmann equation faithfully to extract the hydrodynamics without recourse to any ansatz. Our method leads to microscopic expressions of not only all the transport coefficients that are of the same form as those in Chapman-Enskog method but also those of the viscous relaxation times $τ_i$ that admit physically natural interpretations. As an example, we apply our microscopic expressions to calculate the transport coefficients and the relaxation times of the cold fermionic atoms in a quantitative way, where the transition probability in the collision term is given explicitly in terms of the $s$-wave scattering length $a_s$. We thereby discuss the quantum statistical effects, temperature dependence, and scattering-length dependence of the first-order transport coefficients and the viscous relaxation times: It is shown that as the temperature is lowered, the transport coefficients and the relaxation times increase rapidly because Pauli principle acts effectively. On the other hand, as $a_s$ is increased, these quantities decrease and become vanishingly small at unitarity because of the strong coupling. The numerical calculation shows that the relation $τ_π=η/P$, which is derived in the relaxation-time approximation and used in most of literature without almost any foundation, turns out to be satisfied quite well, while the similar relation for the relaxation time $τ_J$ of the heat conductivity is satisfied only approximately with a considerable error.

arXiv:1604.07458•Apr 25, 2016•1 citation

Density Propagator for Many-Body Localization: Finite-Size Effects, Transient Subdiffusion, and Exponential Decay

Soumya Bera, Giuseppe De Tomasi, Felix Weiner +1 more

We investigate charge relaxation in quantum wires of spinless disordered fermions ($t\text{\ensuremath{-}}V$ model). Our observable is the time-dependent density propagator ${\mathrm{\ensuremath{\Pi}}}_{ϵ}(x,t)$, calculated in windows of different energy density $ϵ$ of the many-body Hamiltonian and at different disorder strengths $W$, not exceeding the critical value ${W}_{c}$. The width ${\mathrm{\ensuremath{\Delta}}x}_{ϵ}(t)$ of ${\mathrm{\ensuremath{\Pi}}}_{ϵ}(x,t)$ exhibits a behavior $d\mathrm{ln}{\mathrm{\ensuremath{\Delta}}x}_{ϵ}(t)/d\mathrm{ln}t={\ensuremath{\beta}}_{ϵ}(t)$, where the exponent function ${\ensuremath{\beta}}_{ϵ}(t)\ensuremath{\lesssim}1/2$ is seen to depend strongly on $L$ at all investigated parameter combinations. (i) We confirm the existence of a region in phase space that exhibits subdiffusive dynamics in the sense that ${\ensuremath{\beta}}_{ϵ}(t)<1/2$ in a large window of times. However, subdiffusion might possibly be transient, only, finally giving way to a conventional diffusive behavior with ${\ensuremath{\beta}}_{ϵ}=1/2$. (ii) We cannot confirm the existence of many-body mobility edges even in regions of the phase diagram that have been reported to be deep in the delocalized phase. (iii) (Transient) subdiffusion $0<{\ensuremath{\beta}}_{ϵ}(t)\ensuremath{\lesssim}1/2$ coexists with an enhanced probability for returning to the origin ${\mathrm{\ensuremath{\Pi}}}_{ϵ}(0,t)$, decaying much slower than $1/{\mathrm{\ensuremath{\Delta}}x}_{ϵ}(t)$. Correspondingly, the spatial decay of ${\mathrm{\ensuremath{\Pi}}}_{ϵ}(x,t)$ is far from Gaussian, being exponential or even slower. On a phenomenological level, our findings are broadly consistent with the effects of strong disorder and (fractal) Griffiths regions.

arXiv:1610.03085•Physical Review Letters•May 8, 2017•121 citations

Spin subdiffusion in perturbed infinite-? Hubbard chain

J. Rękas, Marcin Mierzejewski, Zala Lenarčič +1 more

The t model represents the Hubbard model in the limit U →∞and is one of the basic models of stronglycorrelated electrons. On a one-dimensional chain, the model is integrable, and its charge dynamics is equivalentto that of free spinless fermions. However, within its dynamics the sequence of spins is frozen, implyingHilbert-space fragmentation and nontrivial spin dynamics. We focus primarily on perturbed versions of themodel in which the perturbations break integrability while preserving fragmentation. Considered examplesof such perturbations are spin-dependent hopping and additional spin interaction. We show that these modelsexhibit a variety of spin-transport regimes. In contrast to the integrable model, which displays either ballistictransport or dissipationless diffusion, the perturbed systems exhibit either normal diffusion or subdiffusion. Dueto Hilbert-space fragmentation, spin transport is mediated by charge transport in all cases considered. This leadsto a characteristic dependence on the magnetization, most notably resulting in subdiffusive spin transport in thegrand-canonical ensemble of the perturbed model. The underlying mechanism is distinct from the subdiffusionobserved in disordered systems or in models with dipole conservation.

•Sep 16, 2026•0 citations
We investigate chains of interacting spinless fermions subject to a finite external field FFF (also called Stark chains) and focus on the regime where the charge thermalization follows the subdiffusive hydrodynamics. First, we study reduced models conserving the dipole moment and derive an explicit Einstein relation which links the subdiffusive transport coefficient with the correlations of the dipolar current. This relation explains why the decay rate \ensuremathΓq{\mathrm{\ensuremath{\Gamma}}}_{q}\ensuremathΓq​ of the density modulation with wave vector qqq shows q4{q}^{4}q4 dependence. In the case of the Stark model, a similar Einstein relation is also derived and tested using various numerical methods. They confirm an exponential reduction of the transport coefficient with increasing FFF. On the other hand, our study of the Stark model indicates that upon increasing qqq there is a crossover from subdiffusive behavior, \ensuremathΓq\ensuremath∝q4{\mathrm{\ensuremath{\Gamma}}}_{q}\ensuremath{\propto}{q}^{4}\ensuremathΓq​\ensuremath∝q4, to the normal diffusive relaxation, \ensuremathΓq\ensuremath∝q2{\mathrm{\ensuremath{\Gamma}}}_{q}\ensuremath{\propto}{q}^{2}\ensuremathΓq​\ensuremath∝q2, at the wave vector q∗{q}^{*}q∗ which vanishes for F\ensuremath→0F\ensuremath{\rightarrow}0F\ensuremath→0.
We study transport in an interacting tilted (Stark) chain. We show that the crossover between diffusive and subdiffusive relaxation is governed by FLF\sqrt{L}FL​, where FFF is the strength of the field, and LLL is the wavelength of the excitation. While the subdiffusive behavior persists for large fields, the corresponding transport coefficient is exponentially suppressed with FFF so that the finite-time dynamics appears almost frozen. We (i) explain the crossover scale between the diffusive and subdiffusive transport by bounding the dynamics of the dipole moment for arbitrary initial state, (ii) prove its emergent conservation at infinite temperature for $F\ensuremath{\gg}1/L$, and (iii) argue that the numerical data for the tilted chain are consistent with the hydrodynamics of fractons.
The ttt-model represents the Hubbard model in the limit U→∞U \to \inftyU→∞ and is one of the basic models of strongly correlated electrons. On a one-dimensional chain, the model is integrable, and the charge dynamics corresponds to that of free spinless fermions. However, the sequence of spins is frozen, leading to the Hilbert space fragmentation and nontrivial spin dynamics. We consider integrable and perturbed models with perturbations that break integrability while preserving fragmentation, and show that they exhibit various types of spin dynamics, from ballistic transport to anomalous diffusion in the integrable case, and from diffusion to subdiffusion in the perturbed case. Due to fragmentation, in all cases considered, spin transport is mediated by charge transport, with a particular magnetization dependence, most notably leading to subdiffusion in the grandcanonical average of the perturbed model, with a mechanism distinct from subdiffusion in disordered or dipole-conserving models.
We investigate charge relaxation in the spin-less disordered fermionic Hubbard chain (t−Vt{-}Vt−V-model). Our observable is the time-dependent density propagator, Πε(x,t)\Pi_{\varepsilon}(x,t)Πε​(x,t), calculated in windows of different energy density, ε\varepsilonε, of the many-body Hamiltonian and at different disorder strengths, WWW, not exceeding the critical value WcW_\text{c}Wc​. The width Δxε(t)\Delta x_\varepsilon(t)Δxε​(t) of Πε(x,t)\Pi_\varepsilon(x,t)Πε​(x,t) exhibits a behavior dln⁡Δxε(t)/dln⁡t=βε(t)d\ln \Delta x_\varepsilon(t) / d\ln t {=} \beta_\varepsilon(t)dlnΔxε​(t)/dlnt=βε​(t), where the exponent function βε(t)≲1/2\beta_\varepsilon(t){\lesssim}1/2βε​(t)≲1/2 is seen to depend strongly on LLL at all investigated parameter combinations. (i) We cannot confirm the existence of a region in phase space that exhibits (genuine) subdiffusive dynamics in the sense that βε<1/2\beta_\varepsilon{<}1/2βε​<1/2 is numerically fixed in the limit of large LLL. Instead, subdiffusion might possibly be transient, only, finally giving way to conventional diffusive behavior with βε=1/2\beta_\varepsilon{=}1/2βε​=1/2. (ii) Similarly, we cannot confirm the existence of many-body mobility edges deep in the delocalized phase. (iii) (Transient) subdiffusion 0<βε(t)≲1/20<\beta_\varepsilon(t)\lesssim 1/20<βε​(t)≲1/2, coexists with an enhanced probability for returning to the origin, Πε(0,t)\Pi_\varepsilon(0,t)Πε​(0,t), decaying much slower than 1/Δxε(t)1/\Delta x_\varepsilon (t)1/Δxε​(t). Correspondingly, the spatial decay of Πε(x,t)\Pi_\varepsilon(x,t)Πε​(x,t) is far from Gaussian being exponential or even slower. On a phenomenological level, our findings are broadly consistent with effects of strong disorder and (fractal) Griffiths regions.
We give a detailed derivation of the second-order (local) hydrodynamics for Boltzmann equation with an external force by using the renormalization group method. In this method, we solve the Boltzmann equation faithfully to extract the hydrodynamics without recourse to any ansatz. Our method leads to microscopic expressions of not only all the transport coefficients that are of the same form as those in Chapman-Enskog method but also those of the viscous relaxation times τiτ_iτi​ that admit physically natural interpretations. As an example, we apply our microscopic expressions to calculate the transport coefficients and the relaxation times of the cold fermionic atoms in a quantitative way, where the transition probability in the collision term is given explicitly in terms of the sss-wave scattering length asa_sas​. We thereby discuss the quantum statistical effects, temperature dependence, and scattering-length dependence of the first-order transport coefficients and the viscous relaxation times: It is shown that as the temperature is lowered, the transport coefficients and the relaxation times increase rapidly because Pauli principle acts effectively. On the other hand, as asa_sas​ is increased, these quantities decrease and become vanishingly small at unitarity because of the strong coupling. The numerical calculation shows that the relation τπ=η/Pτ_π=η/Pτπ​=η/P, which is derived in the relaxation-time approximation and used in most of literature without almost any foundation, turns out to be satisfied quite well, while the similar relation for the relaxation time τJτ_JτJ​ of the heat conductivity is satisfied only approximately with a considerable error.
We investigate charge relaxation in quantum wires of spinless disordered fermions ($t\text{\ensuremath{-}}V$ model). Our observable is the time-dependent density propagator ${\mathrm{\ensuremath{\Pi}}}_{ϵ}(x,t)$, calculated in windows of different energy density ϵϵϵ of the many-body Hamiltonian and at different disorder strengths WWW, not exceeding the critical value Wc{W}_{c}Wc​. The width ${\mathrm{\ensuremath{\Delta}}x}_{ϵ}(t)$ of ${\mathrm{\ensuremath{\Pi}}}_{ϵ}(x,t)$ exhibits a behavior $d\mathrm{ln}{\mathrm{\ensuremath{\Delta}}x}_{ϵ}(t)/d\mathrm{ln}t={\ensuremath{\beta}}_{ϵ}(t)$, where the exponent function ${\ensuremath{\beta}}_{ϵ}(t)\ensuremath{\lesssim}1/2$ is seen to depend strongly on LLL at all investigated parameter combinations. (i) We confirm the existence of a region in phase space that exhibits subdiffusive dynamics in the sense that ${\ensuremath{\beta}}_{ϵ}(t)<1/2$ in a large window of times. However, subdiffusion might possibly be transient, only, finally giving way to a conventional diffusive behavior with ${\ensuremath{\beta}}_{ϵ}=1/2$. (ii) We cannot confirm the existence of many-body mobility edges even in regions of the phase diagram that have been reported to be deep in the delocalized phase. (iii) (Transient) subdiffusion $0<{\ensuremath{\beta}}_{ϵ}(t)\ensuremath{\lesssim}1/2$ coexists with an enhanced probability for returning to the origin ${\mathrm{\ensuremath{\Pi}}}_{ϵ}(0,t)$, decaying much slower than $1/{\mathrm{\ensuremath{\Delta}}x}_{ϵ}(t)$. Correspondingly, the spatial decay of ${\mathrm{\ensuremath{\Pi}}}_{ϵ}(x,t)$ is far from Gaussian, being exponential or even slower. On a phenomenological level, our findings are broadly consistent with the effects of strong disorder and (fractal) Griffiths regions.
The t model represents the Hubbard model in the limit U →∞and is one of the basic models of stronglycorrelated electrons. On a one-dimensional chain, the model is integrable, and its charge dynamics is equivalentto that of free spinless fermions. However, within its dynamics the sequence of spins is frozen, implyingHilbert-space fragmentation and nontrivial spin dynamics. We focus primarily on perturbed versions of themodel in which the perturbations break integrability while preserving fragmentation. Considered examplesof such perturbations are spin-dependent hopping and additional spin interaction. We show that these modelsexhibit a variety of spin-transport regimes. In contrast to the integrable model, which displays either ballistictransport or dissipationless diffusion, the perturbed systems exhibit either normal diffusion or subdiffusion. Dueto Hilbert-space fragmentation, spin transport is mediated by charge transport in all cases considered. This leadsto a characteristic dependence on the magnetization, most notably resulting in subdiffusive spin transport in thegrand-canonical ensemble of the perturbed model. The underlying mechanism is distinct from the subdiffusionobserved in disordered systems or in models with dipole conservation.