Publications

Order by: publication date number of citations

• Universality in Anderson localization on random graphs with varying connectivity

  Piotr Sierant, Maciej Lewenstein, Antonello Scardicchio

  SciPost Phys. 15, 045 (2023) · published 3 August 2023 | Toggle abstract · pdf

  We perform a thorough and complete analysis of the Anderson localization transition on several models of random graphs with regular and random connectivity. The unprecedented precision and abundance of our exact diagonalization data (both spectra and eigenstates), together with new finite size scaling and statistical analysis of the graph ensembles, unveils a universal behavior which is described by two simple, integer, scaling exponents. A by-product of such analysis is a reconciliation of the tension between the results of perturbation theory coming from strong disorder and earlier numerical works, which seemed to suggest that there should be a non-ergodic region above a given value of disorder $W_{E}$ which is strictly less than the Anderson localization critical disorder $W_C$, and that of other works which suggest that there is no such region. We find that, although no separate $W_{E}$ exists from $W_C$, the length scale at which fully developed ergodicity is found diverges like $|W-W_C|^{-1}$, while the critical length over which delocalization develops is $\sim |W-W_C|^{-1/2}$. The separation of these two scales at the critical point allows for a true non-ergodic, delocalized region. In addition, by looking at eigenstates and studying leading and sub-leading terms in system size-dependence of participation entropies, we show that the former contain information about the non-ergodicity volume which becomes non-trivial already deep in the delocalized regime. We also discuss the quantitative similarities between the Anderson transition on random graphs and many-body localization transition.

  7 citations

• The classical two-dimensional Heisenberg model revisited: An $SU(2)$-symmetric tensor network study

  Philipp Schmoll, Augustine Kshetrimayum, Jens Eisert, Román Orús, Matteo Rizzi

  SciPost Phys. 11, 098 (2021) · published 29 November 2021 | Toggle abstract · pdf

  The classical Heisenberg model in two spatial dimensions constitutes one of the most paradigmatic spin models, taking an important role in statistical and condensed matter physics to understand magnetism. Still, despite its paradigmatic character and the widely accepted ban of a (continuous) spontaneous symmetry breaking, controversies remain whether the model exhibits a phase transition at finite temperature. Importantly, the model can be interpreted as a lattice discretization of the $O(3)$ non-linear sigma model in $1+1$ dimensions, one of the simplest quantum field theories encompassing crucial features of celebrated higher-dimensional ones (like quantum chromodynamics in $3+1$ dimensions), namely the phenomenon of asymptotic freedom. This should also exclude finite-temperature transitions, but lattice effects might play a significant role in correcting the mainstream picture. In this work, we make use of state-of-the-art tensor network approaches, representing the classical partition function in the thermodynamic limit over a large range of temperatures, to comprehensively explore the correlation structure for Gibbs states. By implementing an $SU(2)$ symmetry in our two-dimensional tensor network contraction scheme, we are able to handle very large effective bond dimensions of the environment up to $\chi_E^\text{eff} \sim 1500$, a feature that is crucial in detecting phase transitions. With decreasing temperatures, we find a rapidly diverging correlation length, whose behaviour is apparently compatible with the two main contradictory hypotheses known in the literature, namely a finite-$T$ transition and asymptotic freedom, though with a slight preference for the second.

  7 citations

• Automatic contraction of unstructured tensor networks

  Adam S. Jermyn

  SciPost Phys. 8, 005 (2020) · published 15 January 2020 | Toggle abstract · pdf

  The evaluation of partition functions is a central problem in statistical physics. For lattice systems and other discrete models the partition function may be expressed as the contraction of a tensor network. Unfortunately computing such contractions is difficult, and many methods to make this tractable require periodic or otherwise structured networks. Here I present a new algorithm for contracting unstructured tensor networks. This method makes no assumptions about the structure of the network and performs well in both structured and unstructured cases so long as the correlation structure is local.

  7 citations

• Static and dynamical signatures of Dzyaloshinskii-Moriya interactions in the Heisenberg model on the kagome lattice

  Francesco Ferrari, Sen Niu, Juraj Hasik, Yasir Iqbal, Didier Poilblanc, Federico Becca

  SciPost Phys. 14, 139 (2023) · published 1 June 2023 | Toggle abstract · pdf

  Motivated by recent experiments on Cs$_2$Cu$_3$SnF$_{12}$ and YCu$_{3}$(OH)$_{6}$Cl$_{3}$, we consider the ${S=1/2}$ Heisenberg model on the kagome lattice with nearest-neighbor super-exchange $J$ and (out-of-plane) Dzyaloshinskii-Moriya interaction $J_D$, which favors (in-plane) ${{\bf Q}=(0,0)}$ magnetic order. By using both variational Monte Carlo and tensor-network approaches, we show that the ground state develops a finite magnetization for $J_D/J \gtrsim 0.03 \mathrm{-} 0.04$; instead, for smaller values of the Dzyaloshinskii-Moriya interaction, the ground state has no magnetic order and, according to the fermionic wave function, develops a gap in the spinon spectrum, which vanishes for $J_D \to 0$. The small value of $J_D/J$ for the onset of magnetic order is particularly relevant for the interpretation of low-temperature behaviors of kagome antiferromagnets, including ZnCu$_{3}$(OH)$_{6}$Cl$_{2}$. For this reason, we assess the spin dynamical structure factor and the corresponding low-energy spectrum, by using the variational Monte Carlo technique. The existence of a continuum of excitations above the magnon modes is observed within the magnetically ordered phase, with a broad peak above the lowest-energy magnons, similarly to what has been detected by inelastic neutron scattering on Cs$_{2}$Cu$_{3}$SnF$_{12}$.

  7 citations

• Many versus one: The disorder operator and entanglement entropy in fermionic quantum matter

  Weilun Jiang, Bin-Bin Chen, Zi Hong Liu, Junchen Rong, Fakher F. Assaad, Meng Cheng, Kai Sun, Zi Yang Meng

  SciPost Phys. 15, 082 (2023) · published 7 September 2023 | Toggle abstract · pdf

  Motivated by recent development of the concept of the disorder operator and its relation with entanglement entropy in bosonic systems, here we show the disorder operator successfully probes many aspects of quantum entanglement in fermionic many-body systems. From both analytical and numerical computations in free and interacting fermion systems in 1D and 2D, we find the disorder operator and the entanglement entropy exhibit similar universal scaling behavior, as a function of the boundary length of the subsystem, but with subtle yet important differences. In 1D they both follow the $\log{L}$ scaling behavior with the coefficient determined by the Luttinger parameter for disorder operator, and the conformal central charge for entanglement entropy. In 2D they both show the universal $L\log L$ scaling behavior in free and interacting Fermi liquid states, with the coefficients depending on the geometry of the Fermi surfaces. However at a 2D quantum critical point with non-Fermi-liquid state, extra symmetry information is needed in the design of the disorder operator, so as to reveal the critical fluctuations as does the entanglement entropy. Our results demonstrate the fermion disorder operator can be used to probe quantum many-body entanglement related to global symmetry, and provide new tools to explore the still largely unknown territory of highly entangled fermion quantum matter in 2 or higher dimensions.

  7 citations

• Density resolved wave packet spreading in disordered Gross-Pitaevskii lattices

  Yagmur Kati, Xiaoquan Yu, Sergej Flach

  SciPost Phys. Core 3, 006 (2020) · published 8 October 2020 | Toggle abstract · pdf

  We perform novel energy and norm density resolved wave packet spreading studies in the disordered Gross-Pitaevskii (GP) lattice to confine energy density fluctuations. We map the locations of GP regimes of weak and strong chaos subdiffusive spreading in the 2D density control parameter space and observe strong chaos spreading over several decades. We obtain a renormalization of the ground state due to disorder, which allows for a new disorder-induced phase of disconnected insulating puddles of matter due to Lifshits tails. Inside this Lifshits phase, the wave packet spreading is substantially slowed down.

  6 citations

• Non-coplanar magnetism, topological density wave order and emergent symmetry at half-integer filling of moiré Chern bands

  Patrick H. Wilhelm, Thomas C. Lang, Mathias S. Scheurer, Andreas M. Läuchli

  SciPost Phys. 14, 040 (2023) · published 20 March 2023 | Toggle abstract · pdf

  Twisted double- and mono-bilayer graphene are graphene-based moiré materials hosting strongly correlated fermions in a gate-tunable conduction band with a topologically non-trivial character. Using unbiased exact diagonalization complemented by unrestricted Hartree-Fock calculations, we find that the strong electron-electron interactions lead to a non-coplanar magnetic state, which has the same symmetries as the tetrahedral antiferromagnet on the triangular lattice and can be thought of as a skyrmion lattice commensurate with the moiré scale, competing with a set of ferromagnetic, topological charge density waves featuring an approximate emergent O(3) symmetry, "rotating" the different charge density wave states into each other. Direct comparison with exact diagonalization reveals that the ordered phases are accurately described within the unrestricted Hartree-Fock approximation. Exhibiting a finite charge gap and Chern number $|C|=1$, the formation of charge density wave order which is intimately connected to a skyrmion lattice phase is consistent with recent experiments on these systems.

  6 citations

• Stripes in the extended $t-t^\prime$ Hubbard model: A Variational Monte Carlo analysis

  Vito Marino, Federico Becca, Luca F. Tocchio

  SciPost Phys. 12, 180 (2022) · published 1 June 2022 | Toggle abstract · pdf

  By using variational quantum Monte Carlo techniques, we investigate the instauration of stripes (i.e., charge and spin inhomogeneities) in the Hubbard model on the square lattice at hole doping $\delta=1/8$, with both nearest- ($t$) and next-nearest-neighbor hopping ($t^\prime$). Stripes with different wavelengths $\lambda$ (denoting the periodicity of the charge inhomogeneity) and character (bond- or site-centered) are stabilized for sufficiently large values of the electron-electron interaction $U/t$. The general trend is that $\lambda$ increases going from negative to positive values of $t^\prime/t$ and decreases by increasing $U/t$. In particular, the $\lambda=8$ stripe obtained for $t^\prime=0$ and $U/t=8$ [L.F. Tocchio, A. Montorsi, and F. Becca, SciPost Phys. 7, 21 (2019)] shrinks to $\lambda=6$ for $U/t\gtrsim 10$. For $t^\prime/t<0$, the stripe with $\lambda=5$ is found to be remarkably stable, while for $t^\prime/t>0$, stripes with wavelength $\lambda=12$ and $\lambda=16$ are also obtained. In all these cases, pair-pair correlations are highly suppressed with respect to the uniform state (obtained for large values of $|t^\prime/t|$), suggesting that striped states are not superconducting at $\delta=1/8$.

  6 citations

• jVMC: Versatile and performant variational Monte Carlo leveraging automated differentiation and GPU acceleration

  Markus Schmitt, Moritz Reh

  SciPost Phys. Codebases 2 (2022) · published 23 August 2022 | Toggle abstract · pdf

  The introduction of Neural Quantum States (NQS) has recently given a new twist to variational Monte Carlo (VMC). The ability to systematically reduce the bias of the wave function ansatz renders the approach widely applicable. However, performant implementations are crucial to reach the numerical state of the art. Here, we present a Python codebase that supports arbitrary NQS architectures and model Hamiltonians. Additionally leveraging automatic differentiation, just-in-time compilation to accelerators, and distributed computing, it is designed to facilitate the composition of efficient NQS algorithms.

  6 citations

• Anisotropy-mediated reentrant localization

  Xiaolong Deng, Alexander L. Burin, Ivan M. Khaymovich

  SciPost Phys. 13, 116 (2022) · published 24 November 2022 | Toggle abstract · pdf

  We consider a 2d dipolar system, $d=2$, with the generalized dipole-dipole interaction $\sim r^{-a}$, and the power $a$ controlled experimentally in trapped-ion or Rydberg-atom systems via their interaction with cavity modes. We focus on the dilute dipolar excitation case when the problem can be effectively considered as single-particle with the interaction providing long-range dipolar-like hopping. We show that the spatially homogeneous tilt $\beta$ of the dipoles giving rise to the anisotropic dipole exchange leads to the non-trivial reentrant localization beyond the locator expansion, $a<d$, unlike the models with random dipole orientation. The Anderson transitions are found to occur at the finite values of the tilt parameter $\beta = a$, $0<a<d$, and $\beta = a/(a-d/2)$, $d/2<a<d$, showing the robustness of the localization at small and large anisotropy values. Both exact analytical methods and extensive numerical calculations show power-law localized eigenstates in the bulk of the spectrum, obeying recently discovered duality $a\leftrightarrow 2d-a$ of their spatial decay rate, on the localized side of the transition, $a>a_{AT}$. This localization emerges due to the presence of the ergodic extended states at either spectral edge, which constitute a zero fraction of states in the thermodynamic limit, decaying though extremely slowly with the system size.

  6 citations