Many-body localization in the interpolating Aubry-André-Fibonacci model

We investigate the localization properties of a spin chain with an antiferromagnetic nearest-neighbor coupling, subject to an external quasiperiodic on-site magnetic field. The quasiperiodic modulation interpolates between two paradigmatic models, namely the Aubry-Andr\'e and the Fibonacci models. We find that stronger many-body interactions extend the ergodic phase in the former, whereas they shrink it in the latter. Furthermore, the many-body localization transition points at the two limits of the interpolation appear to be continuously connected along the deformation of the quasiperiodic modulation. As a result, the position of the many-body localization transition depends on the interaction strength for an intermediate degree of deformation. Moreover, in the region of parameter space where the single-particle spectrum contains both localized and extended states, many-body interactions induce an anomalous effect: weak interactions localize the system, whereas stronger interactions enhance ergodicity. We map the model's localization phase diagram using the decay of the quenched spin imbalance in relatively long chains. This is accomplished employing a time-dependent variational approach applied to a matrix product state decomposition of the many-body state. Our model serves as a rich playground for testing many-body localization under tunable potentials.We investigate the localization properties of a spin chain with an antiferromagnetic nearest-neighbor coupling, subject to an external quasiperiodic on-site magnetic field. The quasiperiodic modulation interpolates between two paradigmatic models, namely the Aubry-Andr\'e and the Fibonacci models. We find that stronger many-body interactions extend the ergodic phase in the former, whereas they shrink it in the latter. Furthermore, the many-body localization transition points at the two limits of the interpolation appear to be continuously connected along the deformation of the quasiperiodic modulation. As a result, the position of the many-body localization transition depends on the interaction strength for an intermediate degree of deformation. Moreover, in the region of parameter space where the single-particle spectrum contains both localized and extended states, many-body interactions induce an anomalous effect: weak interactions localize the system, whereas stronger interactions enhance ergodicity. We map the model's localization phase diagram using the decay of the quenched spin imbalance in relatively long chains. This is accomplished employing a time-dependent variational approach applied to a matrix product state decomposition of the many-body state. Our model serves as a rich playground for testing many-body localization under tunable potentials.

arXiv:2106.13841Physical Review ResearchSep 17, 202132 citations

Many-body delocalization dynamics in long Aubry-André quasiperiodic chains

Elmer V. H. Doggen, A. D. Mirlin

We study quench dynamics in an interacting spin chain with a quasiperiodic on-site field, known as the interacting Aubry-Andr\'e model of many-body localization. Using the time-dependent variational principle, we assess the late-time behavior for chains up to $L=50$ . We find that the choice of periodicity $\mathrm{\ensuremath{\Phi}}$ of the quasiperiodic field influences the dynamics. For $\mathrm{\ensuremath{\Phi}}=(\sqrt{5}\ensuremath{-}1)/2$ (the inverse golden ratio) and interaction $\mathrm{\ensuremath{\Delta}}=1$, the model most frequently considered in the literature, we obtain the critical disorder ${W}_{c}=4.8\ifmmode\pm\else\textpm\fi{}0.5$ in units where the noninteracting transition is at $W=2$ . At the same time, for periodicity $\mathrm{\ensuremath{\Phi}}=\sqrt{2}/2$ we obtain a considerably higher critical value, ${W}_{c}=7.8\ifmmode\pm\else\textpm\fi{}0.5$. Finite-size effects on the critical disorder ${W}_{c}$ are much weaker than in the purely random case. This supports the enhancement of ${W}_{c}$ in the case of a purely random potential by rare ``ergodic spots,'' which do not occur in the quasiperiodic case. Further, the data suggest that the decay of the antiferromagnetic order in the delocalized phase is faster than a power law.We study quench dynamics in an interacting spin chain with a quasiperiodic on-site field, known as the interacting Aubry-Andr\'e model of many-body localization. Using the time-dependent variational principle, we assess the late-time behavior for chains up to $L=50$. We find that the choice of periodicity $\mathrm{\ensuremath{\Phi}}$ of the quasiperiodic field influences the dynamics. For $\mathrm{\ensuremath{\Phi}}=(\sqrt{5}\ensuremath{-}1)/2$ (the inverse golden ratio) and interaction $\mathrm{\ensuremath{\Delta}}=1$, the model most frequently considered in the literature, we obtain the critical disorder ${W}_{c}=4.8\ifmmode\pm\else\textpm\fi{}0.5$ in units where the noninteracting transition is at $W=2$. At the same time, for periodicity $\mathrm{\ensuremath{\Phi}}=\sqrt{2}/2$ we obtain a considerably higher critical value, ${W}_{c}=7.8\ifmmode\pm\else\textpm\fi{}0.5$. Finite-size effects on the critical disorder ${W}_{c}$ are much weaker than in the purely random case. This supports the enhancement of ${W}_{c}$ in the case of a purely random potential by rare ``ergodic spots,'' which do not occur in the quasiperiodic case. Further, the data suggest that the decay of the antiferromagnetic order in the delocalized phase is faster than a power law.

arXiv:1901.06971Physical review. B./Physical review. BSep 13, 201980 citations

Many-body localization in a quasiperiodic Fibonacci chain

Nicolas Macé, Nicolas Laflorencie, Fabien Alet

We study the many-body localization (MBL) properties of a chain of interacting fermions subject to a quasiperiodic potential such that the non-interacting chain is always delocalized and displays multifractality. Contrary to naive expectations, adding interactions in this systems does not enhance delocalization, and a MBL transition is observed. Due to the local properties of the quasiperiodic potential, the MBL phase presents specific features, such as additional peaks in the density distribution. We furthermore investigate the fate of multifractality in the ergodic phase for low potential values. Our analysis is based on exact numerical studies of eigenstates and dynamical properties after a quench.We study the many-body localization (MBL) properties of a chain of interacting fermions subject to a quasiperiodic potential such that the non-interacting chain is always delocalized and displays multifractality. Contrary to naive expectations, adding interactions in this systems does not enhance delocalization, and a MBL transition is observed. Due to the local properties of the quasiperiodic potential, the MBL phase presents specific features, such as additional peaks in the density distribution. We furthermore investigate the fate of multifractality in the ergodic phase for low potential values. Our analysis is based on exact numerical studies of eigenstates and dynamical properties after a quench.

arXiv:1811.01912SciPost PhysicsApr 29, 201964 citations

Many-body localization and delocalization in large quantum chains

Elmer V. H. Doggen, Frank Schindler, Konstantin S. Tikhonov +4 more

We theoretically study the quench dynamics for an isolated Heisenberg spin chain with a random on-site magnetic field, which is one of the paradigmatic models of a many-body localization transition. We use the time-dependent variational principle as applied to matrix product states, which allows us to controllably study chains of a length up to $L=100$ spins, i.e., much larger than $L\ensuremath{\simeq}20$ that can be treated via exact diagonalization. For the analysis of the data, three complementary approaches are used: (i) determination of the exponent $\ensuremath{\beta}$ which characterizes the power-law decay of the antiferromagnetic imbalance with time; (ii) similar determination of the exponent ${\ensuremath{\beta}}_{\mathrm{\ensuremath{\Lambda}}}$ which characterizes the decay of a Schmidt gap in the entanglement spectrum; and (iii) machine learning with the use, as an input, of the time dependence of the spin densities in the whole chain. We find that the consideration of the larger system sizes substantially increases the estimate for the critical disorder ${W}_{c}$ that separates the ergodic and many-body localized regimes, compared to the values of ${W}_{c}$ in the literature. On the ergodic side of the transition, there is a broad interval of the strength of disorder with slow subdiffusive transport. In this regime, the exponents $\ensuremath{\beta}$ and ${\ensuremath{\beta}}_{\mathrm{\ensuremath{\Lambda}}}$ increase, with increasing $L$ , for relatively small $L$ but saturate for $L\ensuremath{\simeq}50$, indicating that these slow power laws survive in the thermodynamic limit. From a technical perspective, we develop an adaptation of the ``learning by confusion'' machine-learning approach that can determine ${W}_{c}$ .We theoretically study the quench dynamics for an isolated Heisenberg spin chain with a random on-site magnetic field, which is one of the paradigmatic models of a many-body localization transition. We use the time-dependent variational principle as applied to matrix product states, which allows us to controllably study chains of a length up to $L=100$ spins, i.e., much larger than $L\ensuremath{\simeq}20$ that can be treated via exact diagonalization. For the analysis of the data, three complementary approaches are used: (i) determination of the exponent $\ensuremath{\beta}$ which characterizes the power-law decay of the antiferromagnetic imbalance with time; (ii) similar determination of the exponent ${\ensuremath{\beta}}_{\mathrm{\ensuremath{\Lambda}}}$ which characterizes the decay of a Schmidt gap in the entanglement spectrum; and (iii) machine learning with the use, as an input, of the time dependence of the spin densities in the whole chain. We find that the consideration of the larger system sizes substantially increases the estimate for the critical disorder ${W}_{c}$ that separates the ergodic and many-body localized regimes, compared to the values of ${W}_{c}$ in the literature. On the ergodic side of the transition, there is a broad interval of the strength of disorder with slow subdiffusive transport. In this regime, the exponents $\ensuremath{\beta}$ and ${\ensuremath{\beta}}_{\mathrm{\ensuremath{\Lambda}}}$ increase, with increasing $L$, for relatively small $L$ but saturate for $L\ensuremath{\simeq}50$, indicating that these slow power laws survive in the thermodynamic limit. From a technical perspective, we develop an adaptation of the ``learning by confusion'' machine-learning approach that can determine ${W}_{c}$.

arXiv:1807.05051Physical review. B./Physical review. BNov 8, 2018186 citations

Possible many-body localized phases in Aubry-André systems

Sara Sharifan

Anderson localization, where quenched disorder leads to the localization of single-particle wave functions, has been studied extensively. However, the effect of interaction and non-random potentials on many-body localization (MBL) has gained attention recently. We investigate the dynamics of the interacting Aubry-André model to assess the stability of the many-body localized (MBL) phase. Our motivation to consider the Aubry–André potential as the disorder potential is to eliminate the avalanche mechanism, which often destabilizes MBL in systems with quenched disorder. We analyze the time evolution and the system-size scaling across a range of disorder strengths for the entanglement entropy, number entropy, and the Hartley number entropy. In the non-interacting case, our results confirm an Anderson localization transition at a critical disorder strength Dc = 2, which is consistent with earlier studies. However, when interactions are introduced, the number entropy appears to grow unboundedly as SN ~ ln ln t in the thermodynamic limit, even for large disorder. Meanwhile, the entanglement entropy grows logarithmically, S ~ ln t, and both entropies saturate in a correlated manner at finite times td ~ e^(L·D). These findings suggest that using a quasiperiodic potential to remove the avalanche mechanism does not stabilize true MBL behavior. Instead, the long-time dynamics are qualitatively similar to those found in systems with random disorder. In this thesis, we first introduce the fundamental concepts of thermalization, quantum quench, and symmetry-resolved entropies, then provide a study of localization, specifically in the Aubry-André model. Finally, we place our numerical results in the broader context of the current debate on the stability of many-body localized phases.Anderson localization, where quenched disorder leads to the localization of single-particle wave functions, has been studied extensively. However, the effect of interaction and non-random potentials on many-body localization (MBL) has gained attention recently. We investigate the dynamics of the interacting Aubry-André model to assess the stability of the many-body localized (MBL) phase. Our motivation to consider the Aubry–André potential as the disorder potential is to eliminate the avalanche mechanism, which often destabilizes MBL in systems with quenched disorder. We analyze the time evolution and the system-size scaling across a range of disorder strengths for the entanglement entropy, number entropy, and the Hartley number entropy. In the non-interacting case, our results confirm an Anderson localization transition at a critical disorder strength Dc = 2, which is consistent with earlier studies. However, when interactions are introduced, the number entropy appears to grow unboundedly as SN ~ ln ln t in the thermodynamic limit, even for large disorder. Meanwhile, the entanglement entropy grows logarithmically, S ~ ln t, and both entropies saturate in a correlated manner at finite times td ~ e^(L·D). These findings suggest that using a quasiperiodic potential to remove the avalanche mechanism does not stabilize true MBL behavior. Instead, the long-time dynamics are qualitatively similar to those found in systems with random disorder. In this thesis, we first introduce the fundamental concepts of thermalization, quantum quench, and symmetry-resolved entropies, then provide a study of localization, specifically in the Aubry-André model. Finally, we place our numerical results in the broader context of the current debate on the stability of many-body localized phases.

Aug 27, 20250 citations

Incommensurate many-body localization in the presence of long-range hopping and single-particle mobility edge

Ke Huang, DinhDuy Vu, Xiao Li +1 more

We study many-body localization (MBL) in the quasiperiodic ${t}_{1}\text{\ensuremath{-}}{t}_{2}$ model, focusing on the role of next-nearest-neighbor (NNN) hopping ${t}_{2}$ , which introduces a single-particle mobility edge. The calculated phase diagram can be divided into three distinct regimes, depending on the strength of the short-range interaction $U$ . For weak interactions ($U\ensuremath{\ll}{t}_{1}$), this model is always nonthermal. For intermediate interactions ($U\ensuremath{\sim}{t}_{1}$), the thermal-MBL phase transition in this model is qualitatively the same as that of the Aubry-Andre (AA) model, which is consistent with existing experimental observations. For strong interactions $(U\ensuremath{\gg}{t}_{1})$, the NNN hopping produces qualitatively new physics because it breaks down the Hilbert space fragmentation present in the AA model. The NNN hopping is thus irrelevant when the interaction is intermediate but relevant for strong interactions.We study many-body localization (MBL) in the quasiperiodic ${t}_{1}\text{\ensuremath{-}}{t}_{2}$ model, focusing on the role of next-nearest-neighbor (NNN) hopping ${t}_{2}$, which introduces a single-particle mobility edge. The calculated phase diagram can be divided into three distinct regimes, depending on the strength of the short-range interaction $U$. For weak interactions ($U\ensuremath{\ll}{t}_{1}$), this model is always nonthermal. For intermediate interactions ($U\ensuremath{\sim}{t}_{1}$), the thermal-MBL phase transition in this model is qualitatively the same as that of the Aubry-Andre (AA) model, which is consistent with existing experimental observations. For strong interactions $(U\ensuremath{\gg}{t}_{1})$, the NNN hopping produces qualitatively new physics because it breaks down the Hilbert space fragmentation present in the AA model. The NNN hopping is thus irrelevant when the interaction is intermediate but relevant for strong interactions.

arXiv:2205.15339Physical review. B./Physical review. BJan 17, 202316 citations

L=50

. We find that the choice of periodicity $\mathrm{\ensuremath{\Phi}}$ of the quasiperiodic field influences the dynamics. For $\mathrm{\ensuremath{\Phi}}=(\sqrt{5}\ensuremath{-}1)/2$ (the inverse golden ratio) and interaction $\mathrm{\ensuremath{\Delta}}=1$, the model most frequently considered in the literature, we obtain the critical disorder ${W}_{c}=4.8\ifmmode\pm\else\textpm\fi{}0.5$ in units where the noninteracting transition is at

W=2

. At the same time, for periodicity $\mathrm{\ensuremath{\Phi}}=\sqrt{2}/2$ we obtain a considerably higher critical value, ${W}_{c}=7.8\ifmmode\pm\else\textpm\fi{}0.5$. Finite-size effects on the critical disorder

{W}_{c}

are much weaker than in the purely random case. This supports the enhancement of

{W}_{c}

in the case of a purely random potential by rare ``ergodic spots,'' which do not occur in the quasiperiodic case. Further, the data suggest that the decay of the antiferromagnetic order in the delocalized phase is faster than a power law.

L=100

spins, i.e., much larger than $L\ensuremath{\simeq}20$ that can be treated via exact diagonalization. For the analysis of the data, three complementary approaches are used: (i) determination of the exponent $\ensuremath{\beta}$ which characterizes the power-law decay of the antiferromagnetic imbalance with time; (ii) similar determination of the exponent ${\ensuremath{\beta}}_{\mathrm{\ensuremath{\Lambda}}}$ which characterizes the decay of a Schmidt gap in the entanglement spectrum; and (iii) machine learning with the use, as an input, of the time dependence of the spin densities in the whole chain. We find that the consideration of the larger system sizes substantially increases the estimate for the critical disorder

{W}_{c}

that separates the ergodic and many-body localized regimes, compared to the values of

{W}_{c}

in the literature. On the ergodic side of the transition, there is a broad interval of the strength of disorder with slow subdiffusive transport. In this regime, the exponents $\ensuremath{\beta}$ and ${\ensuremath{\beta}}_{\mathrm{\ensuremath{\Lambda}}}$ increase, with increasing

L

, for relatively small

L

but saturate for $L\ensuremath{\simeq}50$, indicating that these slow power laws survive in the thermodynamic limit. From a technical perspective, we develop an adaptation of the ``learning by confusion'' machine-learning approach that can determine

{W}_{c}

We study many-body localization (MBL) in the quasiperiodic ${t}_{1}\text{\ensuremath{-}}{t}_{2}$ model, focusing on the role of next-nearest-neighbor (NNN) hopping

{t}_{2}

, which introduces a single-particle mobility edge. The calculated phase diagram can be divided into three distinct regimes, depending on the strength of the short-range interaction

U

. For weak interactions ($U\ensuremath{\ll}{t}_{1}$), this model is always nonthermal. For intermediate interactions ($U\ensuremath{\sim}{t}_{1}$), the thermal-MBL phase transition in this model is qualitatively the same as that of the Aubry-Andre (AA) model, which is consistent with existing experimental observations. For strong interactions $(U\ensuremath{\gg}{t}_{1})$, the NNN hopping produces qualitatively new physics because it breaks down the Hilbert space fragmentation present in the AA model. The NNN hopping is thus irrelevant when the interaction is intermediate but relevant for strong interactions.