Contributor info: Prof. David J. Luitz

Maurits S. J. Tepaske, Dominik Hahn, David J. Luitz

SciPost Phys. 14, 073 (2023) · published 17 April 2023 | Toggle abstract · pdf

Near term quantum computers suffer from a degree of decoherence which is prohibitive for high fidelity simulations with deep circuits. An economical use of circuit depth is therefore paramount. For digital quantum simulation of quantum many-body systems, real time evolution is typically achieved by a Trotter decomposition of the time evolution operator into circuits consisting only of two qubit gates. To match the geometry of the physical system and the CNOT connectivity of the quantum processor, additional SWAP gates are needed. We show that optimal fidelity, beyond what is achievable by simple Trotter decompositions for a fixed gate count, can be obtained by compiling the evolution operator into optimal brickwall circuits for the $S=1/2$ quantum Heisenberg model on chains and ladders, when mapped to one dimensional quantum processors without the need of additional SWAP gates.

Dominik Hahn, Paul A. McClarty, David J. Luitz

SciPost Phys. 11, 074 (2021) · published 6 October 2021 | Toggle abstract · pdf

The fully frustrated ladder - a quasi-1D geometrically frustrated spin one half Heisenberg model - is non-integrable with local conserved quantities on rungs of the ladder, inducing the fragmentation of the Hilbert space into sectors composed of singlets and triplets on rungs. We explore the far-from-equilibrium dynamics of this model through the entanglement entropy and out-of-time-ordered correlators (OTOC). The post-quench dynamics of the entanglement entropy is highly anomalous as it shows clear non-damped revivals that emerge from short connected chunks of triplets and whose persistence is therefore a consequence of fragmentation. We find that the maximum value of the entropy follows from a picture where coherences between different fragments co-exist with perfect thermalization within each fragment. This means that the eigenstate thermalization hypothesis holds within all sufficiently large Hilbert space fragments. The OTOC shows short distance oscillations arising from short coupled fragments, which become decoherent at longer distances, and a sub-ballistic spreading and long distance exponential decay stemming from an emergent length scale tied to fragmentation.

David J. Luitz

SciPost Phys. 11, 021 (2021) · published 3 August 2021 | Toggle abstract · pdf

Periodically driven quantum many-body systems play a central role for our understanding of nonequilibrium phenomena. For studies of quantum chaos, thermalization, many-body localization and time crystals, the properties of eigenvectors and eigenvalues of the unitary evolution operator, and their scaling with physical system size $L$ are of interest. While for static systems, powerful methods for the partial diagonalization of the Hamiltonian were developed, the unitary eigenproblem remains daunting. In this paper, we introduce a Krylov space diagonalization method to obtain exact eigenpairs of the unitary Floquet operator with eigenvalue closest to a target on the unit circle. Our method is based on a complex polynomial spectral transformation given by the geometric sum, leading to rapid convergence of the Arnoldi algorithm. We demonstrate that our method is much more efficient than the shift invert method in terms of both runtime and memory requirements, pushing the accessible system sizes to the realm of 20 qubits, with Hilbert space dimensions $\geq 10^6$.

David J. Luitz, Ivan M. Khaymovich, Yevgeny Bar Lev

SciPost Phys. Core 2, 006 (2020) · published 30 April 2020 | Toggle abstract · pdf

The disordered XXZ model is a prototype model of the many-body localization transition (MBL). Despite numerous studies of this model, the available numerical evidence of multifractality of its eigenstates is not very conclusive due severe finite size effects. Moreover it is not clear if similarly to the case of single-particle physics, multifractal properties of the many-body eigenstates are related to anomalous transport, which is observed in this model. In this work, using a state-of-the-art, massively parallel, numerically exact method, we study systems of up to 24 spins and show that a large fraction of the delocalized phase flows towards ergodicity in the thermodynamic limit, while a region immediately preceding the MBL transition appears to be multifractal in this limit. We discuss the implication of our finding on the mechanism of subdiffusive transport.

Luis Colmenarez, Paul A. McClarty, Masudul Haque, David J. Luitz

SciPost Phys. 7, 064 (2019) · published 21 November 2019 | Toggle abstract · pdf

Ergodic quantum many-body systems satisfy the eigenstate thermalization hypothesis (ETH). However, strong disorder can destroy ergodicity through many-body localization (MBL) -- at least in one dimensional systems -- leading to a clear signal of the MBL transition in the probability distributions of energy eigenstate expectation values of local operators. For a paradigmatic model of MBL, namely the random-field Heisenberg spin chain, we consider the full probability distribution of eigenstate correlation functions across the entire phase diagram. We find gaussian distributions at weak disorder, as predicted by pure ETH. At intermediate disorder -- in the thermal phase -- we find further evidence for anomalous thermalization in the form of heavy tails of the distributions. In the MBL phase, we observe peculiar features of the correlator distributions: a strong asymmetry in $S_i^z S_{i+r}^z$ correlators skewed towards negative values; and a multimodal distribution for spin-flip correlators. A quantitative quasi-degenerate perturbation theory calculation of these correlators yields a surprising agreement of the full distribution with the exact results, revealing, in particular, the origin of the multiple peaks in the spin-flip correlator distribution as arising from the resonant and off-resonant admixture of spin configurations. The distribution of the $S_i^zS_{i+r}^z$ correlator exhibits striking differences between the MBL and Anderson insulator cases.

Francesca Pietracaprina, Nicolas Macé, David J. Luitz, Fabien Alet

SciPost Phys. 5, 045 (2018) · published 6 November 2018 | Toggle abstract · pdf

We provide a pedagogical review on the calculation of highly excited eigenstates of disordered interacting quantum systems which can undergo a many-body localization (MBL) transition, using shift-invert exact diagonalization. We also provide an example code at https://bitbucket.org/dluitz/sinvert_mbl/. Through a detailed analysis of the simulational parameters of the random field Heisenberg spin chain, we provide a practical guide on how to perform efficient computations. We present data for mid-spectrum eigenstates of spin chains of sizes up to $L=26$. This work is also geared towards readers with interest in efficiency of parallel sparse linear algebra techniques that will find a challenging application in the MBL problem.

David J. Luitz, Yevgeny Bar Lev, Achilleas Lazarides

SciPost Phys. 3, 029 (2017) · published 25 October 2017 | Toggle abstract · pdf

Using a numerically exact method we study the stability of dynamical localization to the addition of interactions in a periodically driven isolated quantum system which conserves only the total number of particles. We find that while even infinitesimally small interactions destroy dynamical localization, for weak interactions density transport is significantly suppressed and is asymptotically diffusive, with a diffusion coefficient proportional to the interaction strength. For systems tuned away from the dynamical localization point, even slightly, transport is dramatically enhanced and within the largest accessible systems sizes a diffusive regime is only pronounced for sufficiently small detunings.

David J. Luitz, Nicolas Laflorencie

SciPost Phys. 2, 011 (2017) · published 24 March 2017 | Toggle abstract · pdf

Using quantum Monte Carlo simulations, we compute the participation (Shannon-R\'enyi) entropies for groundstate wave functions of Heisenberg antiferromagnets for one-dimensional (line) subsystems of length $L$ embedded in two-dimensional ($L\times L$) square lattices. We also study the line entropy at finite temperature, i.e. of the diagonal elements of the density matrix, for three-dimensional ($L\times L\times L$) cubic lattices. The breaking of SU(2) symmetry is clearly captured by a universal logarithmic scaling term $l_q\ln L$ in the R\'enyi entropies, in good agreement with the recent field-theory results of Misguish, Pasquier and Oshikawa [arXiv:1607.02465]. We also study the dependence of the log prefactor $l_q$ on the R\'enyi index $q$ for which a transition is detected at $q_c\simeq 1$.