Source author record

David J. Luitz

David J. Luitz appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

22works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2022arXiv

Avalanches and many-body resonances in many-body localized systems

We numerically study both the avalanche instability and many-body resonances in strongly-disordered spin chains exhibiting many-body localization (MBL). We distinguish between a finite-size/time MBL regime, and the asymptotic MBL phase, and identify some "landmarks" within the MBL regime. Our first landmark is an estimate of where the MBL phase becomes unstable to avalanches, obtained by measuring the slowest relaxation rate of a finite chain coupled to an infinite bath at one end. Our estimates indicate that the actual MBL-to-thermal phase transition, in infinite-length systems, occurs much deeper in the MBL regime than has been suggested by most previous studies. Our other landmarks involve system-wide resonances. We find that the effective matrix elements producing eigenstates with system-wide resonances are enormously broadly distributed. This means that the onset of such resonances in typical samples occurs quite deep in the MBL regime, and the first such resonances typically involve rare pairs of eigenstates that are farther apart in energy than the minimum gap. Thus we find that the resonance properties define two landmarks that divide the MBL regime in to three subregimes: (i) at strongest disorder, typical samples do not have any eigenstates that are involved in system-wide many-body resonances; (ii) there is a substantial intermediate regime where typical samples do have such resonances, but the pair of eigenstates with the minimum spectral gap does not; and (iii) in the weaker randomness regime, the minimum gap is involved in a many-body resonance and thus subject to level repulsion. Nevertheless, even in this third subregime, all but a vanishing fraction of eigenstates remain non-resonant and the system thus still appears MBL in many respects. Based on our estimates of the location of the avalanche instability, it might be that the MBL phase is only part of subregime (i).

preprint2022arXiv

Dissipation-Induced Order: The $S=1/2$ Quantum Spin Chain Coupled to an Ohmic Bath

We consider an $S=1/2$ antiferromagnetic quantum Heisenberg chain where each site is coupled to an independent bosonic bath with ohmic dissipation. The coupling to the bath preserves the global SO(3) spin symmetry. Using large-scale, approximation-free quantum Monte Carlo simulations, we show that any finite coupling to the bath suffices to stabilize long-range antiferromagnetic order. This is in stark contrast to the isolated Heisenberg chain where spontaneous breaking of the SO(3) symmetry is forbidden by the Mermin-Wagner theorem. A linear spin-wave theory analysis confirms that the memory of the bath and the concomitant retarded interaction stabilize the order. For the Heisenberg chain, the ohmic bath is a marginal perturbation so that exponentially large system sizes are required to observe long-range order at small couplings. Below this length scale, our numerics is dominated by a crossover regime where spin correlations show different power-law behaviors in space and time. We discuss the experimental relevance of this crossover phenomena.

preprint2022arXiv

Magnetization process and ordering of the $S=1/2$ pyrochlore magnet in a field

We study the $S=1/2$ pyrochlore Heisenberg antiferromagnet in a magnetic field. Using large scale density-matrix renormalization group (DMRG) calculations for clusters up to $128$ spins, we find indications for a finite triplet gap, causing a threshold field to nonzero magnetization in the magnetization curve. We obtain a robust saturation field consistent with a magnon crystal, although the corresponding $5/6$ magnetization plateau is very slim and possibly unstable. Most remarkably, there is a pronounced and apparently robust $\frac 1 2$ magnetization plateau where the groundstate breaks (real-space) rotational symmetry, exhibiting oppositely polarised spins on alternating kagome and triangular planes. Reminiscent of the kagome ice plateau of the pyrochlore Ising magnet known as spin ice, it arises via a much more subtle `quantum order by disorder' mechanism.

preprint2020arXiv

Is there slow particle transport in the MBL phase?

We analyze the saturation value of the bipartite entanglement and number entropy starting from a random product state deep in the MBL phase. By studying the probability distributions of these entropies we find that the growth of the saturation value of the entanglement entropy stems from a significant reshuffling of the weight in the probability distributions from the bulk to the exponential tails. In contrast, the probability distributions of the saturation value of the number entropy are converged with system size, and exhibit a sharp cut-off for values of the number entropy which correspond to one particle fluctuating across the boundary between the two halves of the system. Our results therefore rule out slow particle transport deep in the MBL phase and confirm that the slow entanglement entropy production stems uniquely from configurational entanglement.

preprint2020arXiv

Multifractality and its role in anomalous transport in the disordered XXZ spin-chain

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.

preprint2020arXiv

Statistics of correlation functions in the random Heisenberg chain

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.

preprint2020arXiv

The pyrochlore S=1/2 Heisenberg antiferromagnet at finite temperature

Frustrated three dimensional quantum magnets are notoriously impervious to theoretical analysis. Here we use a combination of three computational methods to investigate the three dimensional pyrochlore $S=1/2$ quantum antiferromagnet, an archetypical frustrated magnet, at finite temperature, $T$: canonical typicality for a finite cluster of $2\times 2 \times 2$ unit cells (i.e. $32$ sites), a finite-$T$ matrix product state method on a larger cluster with $48$ sites, and the numerical linked cluster expansion (NLCE) using clusters up to $25$ lattice sites, which include non-trivial hexagonal and octagonal loops. We focus on thermodynamic properties (energy, specific heat capacity, entropy, susceptibility, magnetisation) next to the static structure factor. We find a pronounced maximum in the specific heat at $T = 0.57 J$, which is stable across finite size clusters and converged in the series expansion. This is well-separated from a residual amount of spectral weight of $0.47 k_B \ln 2$ per spin which has not been released even at $T\approx0.25 J$, the limit of convergence of our results. This is a large value compared to a number of highly frustrated models and materials, such as spin ice or the kagome $S=1/2$ Heisenberg antiferromagnet. We also find a non-monotonic dependence on $T$ of the magnetisation at low magnetic fields, reflecting the dominantly non-magnetic character of the low-energy spectral weight. A detailed comparison of our results to measurements for the $S=1$ material NaCaNi$_2$F$_7$ yields rough agreement of the functional form of the specific heat maximum, which in turn differs from the sharper maximum of the heat capacity of the spin ice material Dy$_2$Ti$_2$O$_7$, all of which are yet qualitatively distinct from conventional, unfrustrated magnets.

preprint2019arXiv

Hierarchy of relaxation timescales in local random Liouvillians

To characterize the generic behavior of open quantum systems, we consider random, purely dissipative Liouvillians with a notion of locality. We find that the positivity of the map implies a sharp separation of the relaxation timescales according to the locality of observables. Specifically, we analyze a spin-1/2 system of size $\ell$ with up to $n$-body Lindblad operators, which are $n$-local in the complexity-theory sense. Without locality ($n=\ell$), the complex Liouvillian spectrum densely covers a "lemon"-shaped support, in agreement with recent findings [Phys. Rev. Lett. 123, 140403;arXiv:1905.02155]. However, for local Liouvillians ($n<\ell$), we find that the spectrum is composed of several dense clusters with random matrix spacing statistics, each featuring a lemon-shaped support wherein all eigenvectors correspond to $n$-body decay modes. This implies a hierarchy of relaxation timescales of $n$-body observables, which we verify to be robust in the thermodynamic limit.

preprint2016arXiv

Anomalous thermalization in ergodic systems

It is commonly believed that quantum isolated systems satisfying the eigenstate thermalization hypothesis (ETH) are diffusive. We show that this assumption is too restrictive, since there are systems that are asymptotically in a thermal state, yet exhibit anomalous, subdiffusive thermalization. We show that such systems satisfy a modified version of the ETH ansatz and derive a general connection between the scaling of the variance of the offdiagonal matrix elements of local operators, written in the eigenbasis of the Hamiltonian, and the dynamical exponent. We find that for subdiffusively thermalizing systems the variance scales more slowly with system size than expected for diffusive systems. We corroborate our findings by numerically studying the distribution of the coefficients of the eigenfunctions and the offdiagonal matrix elements of local operators of the random field Heisenberg chain, which has anomalous transport in its thermal phase. Surprisingly, this system also has non-Gaussian distributions of the eigenfunctions, thus directly violating Berry's conjecture.

preprint2016arXiv

Bimodal entanglement entropy distribution in the many-body localization transition

We introduce the cut averaged entanglement entropy in disordered periodic spin chains and prove it to be a concave function of subsystem size for individual eigenstates. This allows us to identify the entanglement scaling as a function of subsystem size for individual states in inhomogeneous systems. Using this quantity, we probe the critical region between the many-body localized (MBL) and ergodic phases in finite systems. In the middle of the spectrum, we show evidence for bimodality of the entanglement distribution in the MBL critical region, finding both volume law and area law eigenstates over disorder realizations as well as within \emph{single disorder realizations}. The disorder averaged entanglement entropy in this region then scales as a volume law with a coefficient below its thermal value. We discover in the critical region, as we approach the thermodynamic limit, that the cut averaged entanglement entropy density falls on a one-parameter family of curves. Finally, we also show that without averaging over cuts the slope of the entanglement entropy \vs subsystem size can be negative at intermediate and strong disorder, caused by rare localized regions in the system.

preprint2016arXiv

Long tail distributions near the many body localization transition

The random field S=1/2 Heisenberg chain exhibits a dynamical many body localization transition at a critical disorder strength, which depends on the energy density. At weak disorder, the eigenstate thermalization hypothesis (ETH) is fulfilled on average, making local observables smooth functions of energy, whose eigenstate-to-eigenstate fluctuations decrease exponentially with system size. We demonstrate the validity of ETH in the thermal phase as well as its breakdown in the localized phase and show that rare states exist which do not strictly follow ETH, becoming more frequent closer to the transition. Similarly, the probability distribution of the entanglement entropy at intermediate disorder develops long tails all the way down to zero entanglement. We propose that these low entanglement tails stem from localized regions at the subsystem boundaries which were recently discussed as a possible mechanism for subdiffusive transport in the ergodic phase.

preprint2015arXiv

Entanglement entropies of the $J_1 - J_2$ Heisenberg antiferromagnet on the square lattice

Using a modified spin-wave theory which artificially restores zero sublattice magnetization on finite lattices, we investigate the entanglement properties of the Néel ordered $J_1 - J_2$ Heisenberg antiferromagnet on the square lattice. Different kinds of subsystem geometries are studied, either corner-free (line, strip) or with sharp corners (square). Contributions from the $n_G=2$ Nambu-Goldstone modes give additive logarithmic corrections with a prefactor ${n_G}/{2}$ independent of the Rényi index. On the other hand, corners lead to additional (negative) logarithmic corrections with a prefactor $l^{c}_q$ which does depend on both $n_G$ and the Rényi index $q$, in good agreement with scalar field theory predictions. By varying the second neighbor coupling $J_2$ we also explore universality across the Néel ordered side of the phase diagram of the $J_1 - J_2$ antiferromagnet, from the frustrated side $0<J_2/J_1<1/2$ where the area law term is maximal, to the strongly ferromagnetic regime $-J_2/J_1\gg1$ with a purely logarithmic growth $S_q=\frac{n_G}{2}\ln N$, thus recovering the mean-field limit for a subsystem of $N$ sites. Finally, a universal subleading constant term $γ_q^{\rm ord}$ is extracted in the case of strip subsystems, and a direct relation is found (in the large-S limit) with the same constant extracted from free lattice systems. The singular limit of vanishing aspect ratios is also explored, where we identify for $γ_q^\text{ord}$ a regular part and a singular component, explaining the discrepancy of the linear scaling term for fixed width {\it{vs}} fixed aspect ratio subsystems.

preprint2015arXiv

Universal logarithmic corrections to entanglement entropies in two dimensions with spontaneously broken continuous symmetries

We explore the Rényi entanglement entropies of a one-dimensional (line) subsystem of length $L$ embedded in two-dimensional $L\times L$ square lattice for quantum spin models whose ground-state breaks a continuous symmetry in the thermodynamic limit. Using quantum Monte Carlo simulations, we first study the $J_1 - J_2$ Heisenberg model with antiferromagnetic nearest-neighbor $J_1>0$ and ferromagnetic second-neighbor couplings $J_2\le 0$. The signature of SU(2) symmetry breaking on finite size systems, ranging from $L=4$ up to $L=40$ clearly appears as a universal additive logarithmic correction to the Rényi entanglement entropies: $l_q \ln L$ with $l_q\simeq 1$, independent of the Rényi index and values of $J_2$. We confirm this result using a high precision spin-wave analysis (with restored spin rotational symmetry) on finite lattices up to $10^5\times 10^5$ sites, allowing to explore further non-universal finite size corrections and study in addition the case of U(1) symmetry breaking. Our results fully agree with the prediction $l_q=n_G/2$ where $n_G$ is the number of Goldstone modes, by Metlitski and Grover [arXiv:1112.5166].

preprint2014arXiv

Critical Properties of the Superfluid - Bose Glass Transition in Two Dimensions

We investigate the superfluid (SF) to Bose glass (BG) quantum phase transition using extensive quantum Monte Carlo simulations of two-dimensional hard-core bosons in a random box potential. $T=0$ critical properties are studied by thorough finite-size scaling of condensate and SF densities, both vanishing at the same critical disorder $W_c=4.80(5)$. Our results give the following estimates for the critical exponents: $z=1.85(15)$, $ν=1.20(12)$, $η=-0.40(15)$. Furthermore, the probability distribution of the SF response $P(\lnρ_{\rm sf})$ displays striking differences across the transition: while it narrows with increasing system sizes $L$ in the SF phase, it broadens in the BG regime, indicating an absence of self-averaging, and at the critical point $P(\lnρ_{\rm sf}+z \ln L)$ is scale invariant. Finally, high-precision measurements of the local density rule out a percolation picture for the SF-BG transition.

preprint2014arXiv

Improving entanglement and thermodynamic Rényi entropy measurements in quantum Monte Carlo

We present a method for improving measurements of the entanglement Rényi entropies in quantum Monte Carlo simulations by relating them with measurements of participation Rényi entropies. Exploiting the capability of building improved estimators for the latter allows to obtain very good estimates for entanglement Rényi entropies. When considering a full system instead of a bipartition, the method can be further ameliorated providing access to the thermodynamic Rényi entropies with high accuracy. We also explore a recently-proposed method for the reconstruction of the entanglement spectrum from entanglement Rényi entropies and finally show how potential entanglement Hamiltonians may be tested for their validity using a comparison with thermal Rényi entropies.

preprint2014arXiv

Many-body localization edge in the random-field Heisenberg chain

We present a large scale exact diagonalization study of the one dimensional spin $1/2$ Heisenberg model in a random magnetic field. In order to access properties at varying energy densities across the entire spectrum for system sizes up to $L=22$ spins, we use a spectral transformation which can be applied in a massively parallel fashion. Our results allow for an energy-resolved interpretation of the many body localization transition including the existence of an extensive many-body mobility edge. The ergodic phase is well characterized by Gaussian orthogonal ensemble statistics, volume-law entanglement, and a full delocalization in the Hilbert space. Conversely, the localized regime displays Poisson statistics, area-law entanglement and non ergodicity in the Hilbert space where a true localization never occurs. We perform finite size scaling to extract the critical edge and exponent of the localization length divergence.

preprint2014arXiv

Shannon-Rényi entropies and participation spectra across 3d $O(3)$ criticality

Universal features in the scalings of Shannon-Rényi entropies of many-body groundstates are studied for interacting spin-$\frac{1}{2}$ systems across (2+1) dimensional $O(3)$ critical points, using quantum Monte Carlo simulations on dimerized and plaquettized Heisenberg models on the square lattice. Considering both full systems and line shaped subsystems, $SU(2)$ symmetry breaking on the Néel ordered side of the transition is characterized by the presence of a logarithmic term in the scaling of Shannon-Rényi entropies, which is absent in the disordered gapped phase. Such a difference in the scalings allows to capture the quantum critical point using Shannon-Rényi entropies for line shaped subsystems of length $L$ embedded in $L\times L$ tori, as the smaller subsystem entropies are numerically accessible to much higher precision than for the full system. Most interestingly, at the quantum phase transition an additive subleading constant $b_\infty^{*\rm line}=0.41(1)$ emerges in the critical scaling of the line Shannon-Rényi entropy $S_\infty^\text{line}$. This number appears to be universal for 3d $O(3)$ criticality, as confirmed for the finite-temperature transition in the 3d antiferromagnetic spin-$\frac{1}{2}$ Heisenberg model. Additionally, the phases and phase transition can be detected in several features of the participation spectrum, consisting of the diagonal elements of the reduced density matrix of the line subsystem. In particular the Néel ordering transition can be simply understood in the $\{S^z\}$ basis by a confinement mechanism of ferromagnetic domain walls.

preprint2014arXiv

Universal behavior beyond multifractality in quantum many-body systems

How many states of a configuration space contribute to a wave-function? Attempts to answer this ubiquitous question have a long history in physics and chemistry, and are keys to understand e.g. localization phenomena. Quantifying this aspect has often been overlooked for interacting many-body quantum systems, mainly due to the exponential growth of the configuration (Hilbert) space. Here, we introduce two Monte Carlo schemes to calculate Shannon-Renyi entropies for ground-states of large quantum many-body systems that are out of reach for any other exact method. Our simulations reveal that the very nature of quantum phases of matter and associated transitions is captured by universal subleading terms in these entropies, on top of a generic dominant multifractal behavior.

preprint2013arXiv

Single magnetic impurities in the Kane-Mele model

The realization of the spin-Hall effect in quantum wells has led to a plethora of studies regarding the properties of the edge states of a 2D topological insulator. These edge states constitute a class of one-dimensional liquids, called the helical liquid, where an electron's spin quantization axis is tied to its momentum. In contrast to one dimensional conductors, magnetic impurities - below the Kondo temperature - cannot block transport and one expects the current to circumvent the impurity. To study this phenomenon, we consider the single impurity Anderson model embedded into an edge of a Kane-Mele ribbon with up to 512x80 sites and use the numerically exact continuous time QMC method to study the Kondo effect. We present results on the temperature dependence of the spectral properties of the impurity and the bulk system that show the behaviour of the system in the various regimes of the Anderson model. A view complementary to the single particle spectral functions can be obtained using the spatial behaviour of the spin-spin correlation functions. Here we show the characteristic, algebraic decay in the edge channel near the impurity.

preprint2012arXiv

Understanding the Josephson current through a Kondo-correlated quantum dot

We study the Josephson current 0-$π$ transition of a quantum dot tuned to the Kondo regime. The physics can be quantitatively captured by the numerically exact continuous time quantum Monte Carlo method applied to the single-impurity Anderson model with BCS superconducting leads. For a comparison to an experiment the tunnel couplings are determined by fitting the normal-state linear conductance. Excellent agreement for the dependence of the critical Josephson current on the level energy is achieved. For increased tunnel couplings the Kondo scale becomes comparable to the superconducting gap and the regime of the strongest competition between superconductivity and Kondo correlations is reached; we predict the gate voltage dependence of the critical current in this regime.

preprint2011arXiv

Exact diagonalization study of the tunable edge magnetism in graphene

The tunable magnetism at graphene edges with lengths of up to 48 unit cells is analyzed by an exact diagonalization technique. For this we use a generalized interacting one-dimensional model which can be tuned continuously from a limit describing graphene zigzag edge states with a ferromagnetic phase, to a limit equivalent to a Hubbard chain, which does not allow ferromagnetism. This analysis sheds light onto the question why the edge states have a ferromagnetic ground state, while a usual one-dimensional metal does not. Essentially we find that there are two important features of edge states: (a) umklapp processes are completely forbidden for edge states; this allows a spin-polarized ground state. (b) the strong momentum dependence of the effective interaction vertex for edge states gives rise to a regime of partial spin-polarization and a second order phase transition between a standard paramagnetic Luttinger liquid and ferromagnetic Luttinger liquid.

preprint2009arXiv

A weak coupling CTQMC study of the single impurity and periodic Anderson models with s-wave superconducting baths

We apply the unbiased weak-coupling continuous time quantum Monte Carlo (CTQMC) method to review the physics of a single magnetic impurity coupled to s-wave superconducting leads described by the BCS reduced Hamiltonian. As a function of the superconducting gap $Δ$, we study the signature of the first order transition between the singlet and doublet (local moment) states on various quantities. In particular we concentrate on the Josephson current with 0 to $π$ phase shift, the crossing of the Andreev bound states in the single particle spectral function, as well as the local dynamical spin structure factor. Within DMFT, this impurity problem provides a link to the periodic Anderson model with superconducting conduction electrons (BCS-PAM). The first order transition observed in the impurity model is reproduced in the BCS-PAM and is signalized by the crossing of the low energy excitations in the local density of states. The momentum resolved single particle spectral function in the singlet state reveals the coherent, Bloch-like, superposition of Andreev bound states. In the doublet or local moment phase the single particle spectral function is characterized by incoherent quasiparticle excitations.