Local Hilbert Space Research Articles

Comparative study of state-of-the-art matrix-product-state methods for lattice models with large local Hilbert spaces without U(1) symmetry

Lattice models consisting of high-dimensional local degrees of freedom without global particle-number conservation constitute an important problem class in the field of strongly correlated quantum many-body systems. For instance, they are realized in electron-phonon models, cavities, atom-molecule resonance models, or superconductors. In general, these systems elude a complete analytical treatment and need to be studied using numerical methods where matrix-product states (MPSs) provide a flexible and generic ansatz class. Typically, MPS algorithms scale at least quadratic in the dimension of the local Hilbert spaces. Hence, tailored methods, which truncate this dimension, are required to allow for efficient simulations. Here, we describe and compare three state-of-the-art MPS methods each of which exploits a different approach to tackle the computational complexity. We analyze the properties of these methods for the example of the Holstein model, performing high-precision calculations as well as a finite-size-scaling analysis of relevant ground-state observables. The calculations are performed at different points in the phase diagram yielding a comprehensive picture of the different approaches.

Nonergodic dynamics in disorder-free potentials

We review the dynamics of interacting particles in disorder-free potentials concentrating on a combination of a harmonic binding with a constant tilt. We show that a simple picture of an effective local tilt describes a variety of cases. Our examples include spinless fermions (as modeled by Heisenberg spin chain in a magnetic field), spinful fermions as well as bosons that enjoy a larger local on-site Hilbert space. We also discuss the domain-wall dynamics that reveals nonergodic features even for a relatively weak tilt as suggested by Doggen et al. (2020). By adding a harmonic potential on top of the static field we confirm that the surprising regular dynamics is not entirely due to Hilbert space shuttering. It seems better explained by the inhibited transport within the domains of identically oriented spins. Once the spin-1/2 restrictions are lifted as, e.g., for bosons, the dynamics involve stronger entanglement generation. Again for domain wall melting, the effect of the harmonic potential is shown to lead mainly to an effective local tilt.

Spectral statistics in constrained many-body quantum chaotic systems

We study the spectral statistics of spatially-extended many-body quantum systems with on-site Abelian symmetries or local constraints, focusing primarily on those with conserved dipole and higher moments. In the limit of large local Hilbert space dimension, we find that the spectral form factor $K(t)$ of Floquet random circuits can be mapped exactly to a classical Markov circuit, and, at late times, is related to the partition function of a frustration-free Rokhsar-Kivelson (RK) type Hamiltonian. Through this mapping, we show that the inverse of the spectral gap of the RK-Hamiltonian lower bounds the Thouless time $t_{\mathrm{Th}}$ of the underlying circuit. For systems with conserved higher moments, we derive a field theory for the corresponding RK-Hamiltonian by proposing a generalized height field representation for the Hilbert space of the effective spin chain. Using the field theory formulation, we obtain the dispersion of the low-lying excitations of the RK-Hamiltonian in the continuum limit, which allows us to extract $t_{\mathrm{Th}}$. In particular, we analytically argue that in a system of length $L$ that conserves the $m^{th}$ multipole moment, $t_{\mathrm{Th}}$ scales subdiffusively as $L^{2(m+1)}$. We also show that our formalism directly generalizes to higher dimensional circuits, and that in systems that conserve any component of the $m^{th}$ multipole moment, $t_{\mathrm{Th}}$ has the same scaling with the linear size of the system. Our work therefore provides a general approach for studying spectral statistics in constrained many-body chaotic systems.

Higher-order and fractional discrete time crystals in clean long-range interacting systems

Discrete time crystals are periodically driven systems characterized by a response with periodicity nT, with T the period of the drive and n > 1. Typically, n is an integer and bounded from above by the dimension of the local (or single particle) Hilbert space, the most prominent example being spin-1/2 systems with n restricted to 2. Here, we show that a clean spin-1/2 system in the presence of long-range interactions and transverse field can sustain a huge variety of different ‘higher-order’ discrete time crystals with integer and, surprisingly, even fractional n > 2. We characterize these (arguably prethermal) non-equilibrium phases of matter thoroughly using a combination of exact diagonalization, semiclassical methods, and spin-wave approximations, which enable us to establish their stability in the presence of competing long- and short-range interactions. Remarkably, these phases emerge in a model with continous driving and time-independent interactions, convenient for experimental implementations with ultracold atoms or trapped ions.

Clock-dependent spacetime

Einstein’s theory of general relativity is based on the premise that the physical laws take the same form in all coordinate systems. However, it still presumes a preferred decomposition of the total kinematic Hilbert space into local kinematic Hilbert spaces. In this paper, we consider a theory of quantum gravity that does not come with a preferred partitioning of the kinematic Hilbert space. It is pointed out that, in such a theory, dimension, signature, topology and geometry of spacetime depend on how a collection of local clocks is chosen within the kinematic Hilbert space.

Disorder-Free Localization in an Interacting 2D Lattice Gauge Theory.

Disorder-free localization has been recently introduced as a mechanism for ergodicity breaking in low-dimensional homogeneous lattice gauge theories caused by local constraints imposed by gauge invariance. We show that also genuinely interacting systems in two spatial dimensions can become nonergodic as a consequence of this mechanism. This result is all the more surprising since the conventional many-body localization is conjectured to be unstable in two dimensions; hence the gauge invariance represents an alternative robust localization mechanism surviving in higher dimensions in the presence of interactions. Specifically, we demonstrate nonergodic behavior in the quantum link model by obtaining a bound on the localization-delocalization transition through a classical correlated percolation problem implying a fragmentation of Hilbert space on the nonergodic side of the transition. We study the quantum dynamics in this system by introducing the method of "variational classical networks," an efficient and perturbatively controlled representation of the wave function in terms of a network of classical spins akin to artificial neural networks. We identify a distinguishing dynamical signature by studying the propagation of line defects, yielding different light cone structures in the localized and ergodic phases, respectively. The methods we introduce in this work can be applied to any lattice gauge theory with finite-dimensional local Hilbert spaces irrespective of spatial dimensionality.

Scrambling and Lyapunov exponent in spatially extended systems

Scrambling of information in a quantum many-body system, quantified by the out-of-time-ordered correlator (OTOC), is a key manifestation of quantum chaos. A regime of exponential growth in the OTOC, characterized by a Lyapunov exponent, has so far mostly been observed in systems with a high-dimensional local Hilbert space and in weakly-coupled systems. Here, we propose a general criterion for the existence of a well-defined regime of exponential growth of the OTOC in spatially extended systems with local interactions. In such systems, we show that a parametrically long period of exponential growth requires the butterfly velocity to be much larger than the Lyapunov exponent times a microscopic length scale, such as the lattice spacing. As an explicit example, we study a random unitary circuit with tunable interactions. In this model, we show that in the weakly interacting limit the above criterion is satisfied, and there is a prolonged window of exponential growth. Our results are based on numerical simulations of both Clifford and universal random circuits supported by an analytical treatment.

Ergodic and Nonergodic Dual-Unitary Quantum Circuits with Arbitrary Local Hilbert Space Dimension.

Dual-unitary quantum circuits can be used to construct 1+1 dimensional lattice models for which dynamical correlations of local observables can be explicitly calculated. We show how to analytically construct classes of dual-unitary circuits with any desired level of (non-)ergodicity for any dimension of the local Hilbert space, and present analytical results for thermalization to an infinite-temperature Gibbs state (ergodic) and a generalized Gibbs ensemble (nonergodic). It is shown how a tunable ergodicity-inducing perturbation can be added to a nonergodic circuit without breaking dual unitarity, leading to the appearance of prethermalization plateaux for local observables.

Efficient and flexible approach to simulate low-dimensional quantum lattice models with large local Hilbert spaces

Quantum lattice models with large local Hilbert spaces emerge across various fields in quantum many-body physics. Problems such as the interplay between fermions and phonons, the BCS-BEC crossover of interacting bosons, or decoherence in quantum simulators have been extensively studied both theoretically and experimentally. In recent years, tensor network methods have become one of the most successful tools to treat such lattice systems numerically. Nevertheless, systems with large local Hilbert spaces remain challenging. Here, we introduce a mapping that allows to construct artificial U(1)U(1) symmetries for any type of lattice model. Exploiting the generated symmetries, numerical expenses that are related to the local degrees of freedom decrease significantly. This allows for an efficient treatment of systems with large local dimensions. Further exploring this mapping, we reveal an intimate connection between the Schmidt values of the corresponding matrixproductstate representation and the singlesite reduced density matrix. Our findings motivate an intuitive physical picture of the truncations occurring in typical algorithms and we give bounds on the numerical complexity in comparison to standard methods that do not exploit such artificial symmetries. We demonstrate this new mapping, provide an implementation recipe for an existing code, and perform example calculations for the Holstein model at half filling. We studied systems with a very large number of lattice sites up to L=501L=501 while accounting for N_{\rm ph}=63 phonons per site with high precision in the CDW phase.

Separation of quantum, spatial quantum, and approximate quantum correlations

Quantum nonlocal correlations are generated by implementation of local quantum measurements on spatially separated quantum subsystems. Depending on the underlying mathematical model, various notions of sets of quantum correlations can be defined. In this paper we prove separations of such sets of quantum correlations. In particular, we show that the set of bipartite quantum correlations with four binary measurements per party becomes strictly smaller once we restrict the local Hilbert spaces to be finite dimensional, i.e., Cq(4,4,2,2)≠Cqs(4,4,2,2). We also prove non-closure of the set of bipartite quantum correlations with four ternary measurements per party, i.e., Cqs(4,4,3,3)≠Cqa(4,4,3,3).

Finite-temperature density-matrix renormalization group method for electron-phonon systems: Thermodynamics and Holstein-polaron spectral functions

We investigate the thermodynamics and finite-temperature spectral functions of the Holstein polaron using a density-matrix renormalization group method. Our method combines purification and local basis optimization (LBO) as an efficient treatment of phonon modes. LBO is a scheme which relies on finding the optimal local basis by diagonalizing the local reduced density matrix. By transforming the state into this basis, one can truncate the local Hilbert space with a negligible loss of accuracy for a wide range of parameters. In this work, we focus on the crossover regime between large and small polarons of the Holstein model. Here, no analytical solution exists and we show that the thermal expectation values at low temperatures are independent of the phonon Hilbert space truncation provided the basis is chosen large enough. We then demonstrate that we can extract the electron spectral function and establish consistency with results from a finite-temperature Lanczos method. We additionally calculate the electron emission spectrum and the phonon spectral function and show that all the computations are significantly simplified by the local basis optimization. We observe that the electron emission spectrum shifts spectral weight to both lower frequencies and larger momenta as the temperature is increased. The phonon spectral function experiences a large broadening and the polaron peak at large momenta gets significantly flattened and merges almost completely into the free-phonon peak.

Real-time dynamics in 2+1D compact QED using complex periodic Gaussian states

We introduce a class of variational states to study ground state properties and real-time dynamics in (2+1)-dimensional compact QED. These are based on complex Gaussian states which are made periodic in order to account for the compact nature of the $U(1)$ gauge field. Since the evaluation of expectation values involves infinite sums, we present an approximation scheme for the whole variational manifold. We calculate the ground state energy density for lattice sizes up to $20 \times 20$ and extrapolate to the thermodynamic limit for the whole coupling region. Additionally, we study the string tension both by fitting the potential between two static charges and by fitting the exponential decay of spatial Wilson loops. As the ansatz does not require a truncation in the local Hilbert spaces, we analyze truncation effects which are present in other approaches. The variational states are benchmarked against exact solutions known for the one plaquette case and exact diagonalization results for a $\mathbb{Z}_3$ lattice gauge theory. Using the time-dependent variational principle, we study real-time dynamics after various global quenches, e.g. the time evolution of a strongly confined electric field between two charges after a quench to the weak-coupling regime. Up to the points where finite size effects start to play a role, we observe equilibrating behavior.

Tightening the Lieb-Robinson Bound in Locally Interacting Systems

The Lieb-Robinson (LR) bound rigorously shows that in quantum systems with short-range interactions, the maximum amount of information that travels beyond an effective "light cone" decays exponentially with distance from the light-cone front, which expands at finite velocity. Despite being a fundamental result, existing bounds are often extremely loose, limiting their applications. We introduce a method that dramatically and qualitatively improves LR bounds in models with finite-range interactions. Most prominently, in systems with a large local Hilbert space dimension $D$, our method gives an LR velocity that grows much slower than previous bounds with $D$ as $D\to \infty$. For example, in the Heisenberg model with spin $S$, we find $v\leq$ const. compared to the previous $v\propto S$ which diverges at large $S$, and in multiorbital Hubbard models with $N$ orbitals, we find $v\propto \sqrt{N}$ instead of previous $v\propto N$, and similarly in the $N$-state truncated Bose-Hubbard model and Wen's quantum rotor model. Our bounds also scale qualitatively better in some systems when the spatial dimension or certain model parameters become large, for example in the $d$-dimensional quantum Ising model and perturbed toric code models. Even in spin-1/2 Ising and Fermi-Hubbard models, our method improves the LR velocity by an order of magnitude with typical model parameters, and significantly improves the LR bound at large distance and early time.

Subdiffusion in strongly tilted lattice systems

The quantum dynamics away from equilibrium is of fundamental interest for interacting many-body systems. In this letter, we study tilted many-body systems using the effective Hamiltonian derived from the microscopic description. We first give general arguments for the density relaxation rate satisfying $1/\tau\propto k^4$ for a large class of systems, including the Fermi Hubbard model case as observed in the the recent experiment [1]. Here $k$ is the wave vector of the density wave. The main ingredients are the emergence of the reflection symmetry and dipole moment conservation to the leading non-trivial order of the large tilted strength. To support our analysis, we then construct a solvable model with large local Hilbert space dimension by coupling sites discribed by the Sachdev-Ye-Kitaev models, where the density response can be computed explicitly. The the tilt strength and the temperature dependence of the subdiffusion constant are also discussed.

Reconstructing quantum molecular rotor ground states

Nanomolecular assemblies of C$_{60}$ can be synthesized to enclose dipolar molecules. The low-temperature states of such endofullerenes are described by quantum mechanical rotors, which are candidates for quantum information devices with higher-dimensional local Hilbert spaces. The experimental exploration of endofullerene arrays comes at a time when machine learning techniques are rapidly being adopted to characterize, verify, and reconstruct quantum states from measurement data. In this paper, we develop a strategy for reconstructing the ground state of chains of dipolar rotors using restricted Boltzmann machines (RBMs) adapted to train on data from higher-dimensional Hilbert spaces. We demonstrate accurate generation of energy expectation values from an RBM trained on data in the free-rotor eigenstate basis, and explore the learning resources required for various chain lengths and dipolar interaction strengths. Finally, we show evidence for fundamental limitations in the accuracy achievable by RBMs due to the difficulty in imposing symmetries in the sampling procedure. We discuss possible avenues to overcome this limitation in the future, including the further development of autoregressive models such as recurrent neural networks for the purposes of quantum state reconstruction.

A model of quantum gravity with emergent spacetime

We construct a model of quantum gravity in which dimension, topology and geometry of spacetime are dynamical. The microscopic degree of freedom is a real rectangular matrix whose rows label internal flavours, and columns label spatial sites. In the limit that the size of the matrix is large, the sites can collectively form a spatial manifold. The manifold is determined from the pattern of entanglement present across local Hilbert spaces associated with column vectors of the matrix. With no structure of manifold fixed in the background, the spacetime gauge symmetry is generalized to a group that includes diffeomorphism in arbitrary dimensions. The momentum and Hamiltonian that generate the generalized diffeomorphism obey a first-class constraint algebra at the quantum level. In the classical limit, the constraint algebra of the general relativity is reproduced as a special case. The first-class nature of the algebra allows one to express the projection of a quantum state of the matrix to a gauge invariant state as a path integration of dynamical variables that describe collective fluctuations of the matrix. The collective variables describe dynamics of emergent spacetime, where multi-fingered times arise as Lagrangian multipliers that enforce the gauge constraints. If the quantum state has a local structure of entanglement, a smooth spacetime with well-defined dimension, topology, signature and geometry emerges at the saddle-point, and the spin two mode that determines the geometry can be identified. We find a saddle-point solution that describes a series of (3 + 1)-dimensional de Sitter-like spacetimes with the Lorentzian signature bridged by Euclidean spaces in between. The phase transitions between spacetimes with different signatures are caused by Lifshitz transitions in which the pattern of entanglement is rearranged across the system. Fluctuations of the collective variables are described by bi-local fields that propagate in the spacetime set up by the saddle-point solution.

Entanglement growth in diffusive systems

Entanglement helps in understanding diverse phenomena, going from quantifying complexity to classifying phases of matter. Here we study the influence of conservation laws on entanglement growth. Focusing on systems with U(1) symmetry, i.e., conservation of charge or magnetization, that exhibits diffusive dynamics, we theoretically predict the growth of entanglement, as quantified by the Rényi entropy, in lattice systems in any spatial dimension d and for any local Hilbert space dimension q (qudits). We find that the growth depends both on d and q, and is in generic case first linear in time, similarly as for systems without any conservation laws. Exception to this rule are chains of 2-level systems where the dependence is a square-root of time at all times. Predictions are numerically verified by simulations of diffusive Clifford circuits with upto ~ 105 qubits. Such efficiently simulable circuits should be a useful tool for other many-body problems.

Correlation measures and the entanglement wedge cross-section after quantum quenches in two-dimensional conformal field theories

We consider the time evolution of mixed state correlation measures in two-dimensional conformal field theories, such as logarithmic negativity, odd entropy, and reflected entropy, after quantum quenches of various kinds. These correlation measures, in the holographic context, are all associated to the entanglement wedge cross section. We contrast various classes of conformal field theories, both rational and irrational (pure) conformal field theories. First, for rational conformal field theories, whose dynamics can be well described by the quasi-particle picture, we find all four quantities for disjoint intervals to be proportional, regardless of the specific quench protocol. Second, using the light cone boot- strap, we generalize our results to irrational conformal field theories where we find sharp distinctions from the quasi-particle results and striking differences between mutual information and the other measures. The large surplus of logarithmic negativity relative to mutual information forces us to reconsider what mutual information and logarithmic negativity really measure. We interpret these results as a signature of information scrambling and chaos in irrational theories. These CFT results perfectly agree with our gravitational (holographic) calculations. Furthermore, using holography, we are able to generalize the results to outside of the light cone limit. Finally, due to the breakdown of the quasi-particle picture for irrational theories, we appeal to the “line-tension picture,” motivated by random unitary circuits, as a phenomenological description. We observe that random unitary circuits, with local Hilbert space dimension determined by the Cardy formula, have precisely the same entanglement dynamics as irrational (including holographic) conformal field theories.

Theory of the phase transition in random unitary circuits with measurements

We present a theory of the entanglement transition tuned by measurement strength in qudit chains evolved by random unitary circuits and subject to either weak or random projective measurements. The transition can be understood as a nonanalytic change in the amount of information extracted by the measurements about the initial state of the system, quantified by the Fisher information. To compute the von~Neumann entanglement entropy $S$ and the Fisher information $\mathcal{F}$, we apply a replica method based on a sequence of quantities $\tilde{S}^{(n)}$ and $\mathcal{F}^{(n)}$ that depend on the $n$-th moments of density matrices and reduce to $S$ and $\mathcal{F}$ in the limit $n\to 1$. These quantities with $n\ge 2$ are mapped to free energies of a classical spin model with $n!$ internal states in two dimensions with specific boundary conditions. In particular, $\tilde{S}^{(n)}$ is the excess free energy of a domain wall terminating on the top boundary, and $\mathcal{F}^{(n)}$ is related to the magnetization on the bottom boundary. Phase transitions occur as the spin models undergo ordering transitions in the bulk. Taking the limit of large local Hilbert space dimension $q$ followed by the replica limit $n\to 1$, we obtain the critical measurement probability $p_c=1/2$ and identify the transition as a bond percolation in the 2D square lattice in this limit. Finally, we show there is no phase transition if the measurements are allowed in an arbitrary nonlocal basis, thereby highlighting the relation between the phase transition and information scrambling. We establish an explicit connection between the entanglement phase transition and the purification dynamics of a mixed state evolution and discuss implications of our results to experimental observations of the transition and simulability of quantum dynamics.

On Complex Gamma-Function Integrals

It was observed recently that relations between matrix elements of certain operators in the ${\rm SL}(2,\mathbb R)$ spin chain models take the form of multidimensional integrals derived by R.A. Gustafson. The spin magnets with ${\rm SL}(2,\mathbb C)$ symmetry group and ${\rm L}_2(\mathbb C)$ as a local Hilbert space give rise to a new type of $\Gamma$-function integrals. In this work we present a direct calculation of two such integrals. We also analyse properties of these integrals and show that they comprise the star-triangle relations recently discussed in the literature. It is also shown that in the quasi-classical limit these integral identities are reduced to the duality relations for Dotsenko-Fateev integrals.