Search

Quantum link models (QLMs) offer the realistic prospect for the practical implementation of lattice quantum electrodynamics (QED) on modern quantum simulators, and they provide a venue for exploring various nonergodic phenomena relevant to quantum many-body physics. In these models, gauge and electric fields are represented by spin-$S$ operators. So far, large-scale realizations of QLMs have been restricted to $S=1/2$ representations, whereas the lattice-QED limit is approached at $S\to\infty$. Here, we present a bosonic mapping for the representation of gauge and electric fields with effective spin-$S$ operators for arbitrarily large values of $S$. Based on this mapping, we then propose an experimental scheme for the realization of a large-scale spin-$1$ $\mathrm{U}(1)$ QLM using spinless bosons in an optical superlattice. Using perturbation theory and infinite matrix product state calculations, which work directly in the thermodynamic limit, we demonstrate the faithfulness of the mapping and stability of gauge invariance throughout all accessible evolution times. We further demonstrate the potential of our proposed quantum simulator to address relevant high-energy physics by probing the (de)confinement of an electron--positron pair by tuning the gauge coupling. Our work provides an essential step towards gauge-theory quantum simulators in the quantum-field-theory limit., $15$ pages, $8$ figures

A universal finite system-size scaling analysis of the entanglement entropy is presented for highly degenerate ground states arising from spontaneous symmetry breaking with type-B Goldstone modes in exactly solvable one-dimensional quantum many-body systems. These states appear to be scale-invariant, but not conformally invariant. Our findings are based on a physical argument, imposing three constraints on the entanglement entropy, in addition to further confirmation from an asymptotic analysis of the entanglement entropy for the ${\rm SU}(2)$ spin-$1/2$ ferromagnetic states. The resulting universal scaling form is demonstrated for three fundamental models -- the ${\rm SU}(2)$ spin-$s$ Heisenberg ferromagnetic model, the ${\rm SU}(N+1)$ ferromagnetic model, and the staggered ${\rm SU}(3)$ spin-1 ferromagnetic biquadratic model. The results point towards a classification for distinct types of scale-invariant states, relevant to a complete classification of quantum states of matter., 8 pages, 3 figures

Ferromagnetic ground states have often been overlooked in comparison to seemingly more interesting antiferromagnetic ground states. However, both the physical and mathematical structure of ferromagnetic ground states are particularly rich. We show that the highly degenerate and highly entangled ground states of the ferromagnetic spin-1 biquadratic model are scale invariant, originating from spontaneous symmetry breaking from ${\rm SU}(3)$ to ${\rm U}(1)\times {\rm U}(1)$ with two type-B Goldstone modes. The ground state degeneracies are characterized as the Fibonacci-Lucas sequences -- an ancient mathematical gem, under open and periodic boundary conditions, with the residual entropy being non-zero. This implies that the ground state degeneracies for this model are asymptotically the golden spiral. In addition, sequences of degenerate ground states generated from highest and generalized highest weight states are constructed to establish that the entanglement entropy scales logarithmically with the block size in the thermodynamic limit, with the prefactor being half the number of type-B Goldstone modes. The latter in turn is identified to be the fractal dimension., 20 pages, 3 figures,1 table, additional references to earlier related work added

Disorder-free localization (DFL) has been established as a mechanism of strong ergodicity breaking in $1+1$D lattice gauge theories (LGTs) with dynamical matter for quenches starting in homogeneous initial states that span an extensive number of gauge superselection sectors. Nevertheless, the fate of DFL in $2+1$D in the presence of dynamical matter has hitherto remained an open question of great interest in light of the instability of quenched-disorder many-body localization in higher spatial dimensions. Using infinite matrix product state calculations, we show that DFL survives in $2+1$D LGTs, albeit it is generally less pronounced than in $1+1$D, and highly depends on the matter configuration of the initial state. Through suitable matter configurations, we are able to relate and compare the $1+1$D and $2+1$D cases, showing that the main ingredient for the strength of DFL in our setup is the \textit{propagation directionality} of matter. Our results suggest that, generically, DFL is weakened with increasing spatial dimension, although it can be made independent of the latter by minimizing the propagation directionality of matter in the initial state., $7$ pages, $4$ figures

Using machine learning (ML) to recognize different phases of matter and to infer the entire phase diagram has proven to be an effective tool given a large dataset. In our previous proposals, we have successfully explored phase transitions for topological phases of matter at low dimensions either in a supervised or an unsupervised learning protocol with the assistance of quantum information related quantities. In this work, we adopt our previous ML procedures to study quantum phase transitions of magnetism systems such as the XY and XXZ spin chains by using spin-spin correlation functions as the input data. We find that our proposed approach not only maps out the phase diagrams with accurate phase boundaries, but also indicates some new features that have not observed before. In particular, we define so-called relevant correlation functions to some corresponding phases that can always distinguish between those and their neighbors. Based on the unsupervised learning protocol we proposed [Phys. Rev. B 104, 165108 (2021)], the reduced latent representations of the inputs combined with the clustering algorithm show the connectedness or disconnectedness between neighboring clusters (phases), just corresponding to the continuous or disrupt quantum phase transition, respectively., 18 pages, 21 figures

The quantum simulation of gauge theories on synthetic quantum matter devices has gained a lot of traction in the last decade, making possible the observation of a range of exotic quantum many-body phenomena. In this work, we consider the spin-$1/2$ quantum link formulation of $1+1$D quantum electrodynamics with a topological $\theta$-angle, which can be used to tune a confinement-deconfinement transition. Exactly mapping this system onto a PXP model with mass and staggered magnetization terms, we show an intriguing interplay between confinement and the ergodicity-breaking paradigms of quantum many-body scarring and Hilbert-space fragmentation. We map out the rich dynamical phase diagram of this model, finding an ergodic phase at small values of the mass $\mu$ and confining potential $\chi$, an emergent integrable phase for large $\mu$, and a fragmented phase for large values of both parameters. We also show that the latter hosts resonances that lead to a vast array of effective models. We propose experimental probes of our findings, which can be directly accessed in current cold-atom setups., Comment: 20+11 pages, 16+8 figures. Changes relative to v1: The section on quantum many-body scarring with confinement is expanded. In particular, an analytic derivation for the optimal confining potential is added. The MPS simulations of the Bose-Hubbard model are redone with a larger system size. Several references are added

The topological $\theta$-angle in gauge theories engenders a series of fundamental phenomena, including violations of charge-parity (CP) symmetry, dynamical topological transitions, and confinement--deconfinement transitions. At the same time, it poses major challenges for theoretical studies, as it implies a sign problem in numerical simulations. Analog quantum simulators open the promising prospect of treating quantum many-body systems with such topological terms, but, contrary to their digital counterparts, they have not yet demonstrated the capacity to control the $\theta$-angle. Here, we demonstrate how a tunable topological $\theta$-term can be added to a prototype theory with $\mathrm{U}(1)$ gauge symmetry, a discretized version of quantum electrodynamics in one spatial dimension. As we show, the model can be realized experimentally in a single-species Bose--Hubbard model in an optical superlattice with three different spatial periods, thus requiring only standard experimental resources. Through numerical calculations obtained from the time-dependent density matrix renormalization group method and exact diagonalization, we benchmark the model system, and illustrate how salient effects due to the $\theta$-term can be observed. These include charge confinement, the weakening of quantum many-body scarring, as well as the disappearance of Coleman's phase transition due to explicit breaking of CP symmetry. This work opens the door towards studying the rich physics of topological gauge-theory terms in large-scale cold-atom quantum simulators., Comment: 20 pages, 11 figures, journal article

A major driver of quantum-simulator technology is the prospect of probing high-energy phenomena in synthetic quantum matter setups at a high level of control and tunability. Here, we propose an experimentally feasible realization of a large-scale $2+1$D $\mathrm{U}(1)$ gauge theory with dynamical matter and gauge fields in a cold-atom quantum simulator with spinless bosons. We present the full mapping of the corresponding Gauss's law onto the bosonic computational basis. We then show that the target gauge theory can be faithfully realized and stabilized by an emergent gauge protection term in a two-dimensional single-species Bose--Hubbard optical Lieb superlattice with two spatial periods along either direction, thereby requiring only moderate experimental resources already available in current cold-atom setups. Using infinite matrix product states, we calculate numerical benchmarks for adiabatic sweeps and global quench dynamics that further confirm the fidelity of the mapping. Our work brings quantum simulators of gauge theories a significant step forward in terms of investigating particle physics in higher spatial dimensions, and is readily implementable in existing cold-atom platforms., 8+5 pages, 3+2 figures

Recently, single-site observables have been shown to be useful for the detection of dynamical criticality due to an emergence of a universal critically-prethermal temporal regime in the magnetization [arXiv:2105.05986]. Here, we explore the potential of single-site observables as probes of quantum phase transitions in integrable and nonintegrable transverse-field Ising chains (TFIC). We analytically prove the requirement of zero modes for a quasi-stationary temporal regime to emerge at a bulk probe site, and show how this regime gives rise to a non-analytic behavior in the dynamical order profiles. Our $t$-DMRG calculations verify the results of the quench mean-field theory for near-integrable TFIC both with finite-size and finite-time scaling analyses. We find that both finite-size and finite-time analyses suggest a dynamical critical point for a strongly nonintegrable and locally connected TFIC. We finally demonstrate the presence of a quasi-stationary regime in the power-law interacting TFIC, and extract local dynamical order profiles for TFIC in the long-range Ising universality class with algebraic light cones., Comment: 22 pages, 23 figures. Revised with new analytical calculations, new content and improved presentation

Density matrix renormalization group (DMRG) is one of the most powerful numerical methods available for many-body systems. In this work, we develop a perturbation theory of DMRG (PT-DMRG) to largely increase its accuracy in an extremely simple and efficient way. By using the canonical matrix product state (MPS) representation for the ground state of the considered system, a set of orthogonal basis functions $\left\lbrace | \psi_i \rangle \right\rbrace$ is introduced to describe the perturbations to the ground state obtained by the conventional DMRG. The Schmidt numbers of the MPS that are beyond the bond dimension cut-off are used to define such perturbation terms. The perturbed Hamiltonian is then defined as $\tilde{H}_{ij}= \langle \psi_i | \hat{H} | \psi_j \rangle$; its ground state permits to calculate physical observables with a considerably improved accuracy as compared to the original DMRG results. We benchmark the second-order perturbation theory with the help of one-dimensional Ising chain in a transverse field and the Heisenberg chain, where the precision of DMRG is shown to be improved $\rm O(10)$ times. Furthermore, for moderate length $L$ the errors of DMRG and PT-DMRG both scale linearly with $L^{-1}$. The linear relation between the dimension cut-off of DMRG and that of PT-DMRG with the same precision shows a considerable improvement of efficiency, especially for large dimension cut-off's. In thermodynamic limit we show that the errors of PT-DMRG scale with $\sqrt{L^{-1}}$. Our work suggests an effective way to define the tangent space of the ground state MPS, which may shed lights on the properties beyond the ground state. Such second-order PT-DMRG can be readily generalized to higher orders, as well as applied to the models in higher dimensions.

Using numerical techniques, we study the miscible-immiscible quantum phase transition in a linearly coupled binary Bose-Hubbard model Hamiltonian that can describe low-energy properties of a two-component Bose-Einstein condensate in optical lattices. With the quantum many-body ground state obtained from density matrix renormalization group algorithm, we calculate the characteristic physical quantities of the phase transition controlled by the linear coupling between two components. Furthermore we calculate the Binder cumulant to determine the critical point and draw the phase diagram. The strong-coupling expansion shows that in the Mott insulator regime the model Hamiltonian can be mapped to a spin 1/2 XXZ model with a transverse magnetic field., Comment: 10 pages, 10 figures, submitted to Phys. Rev. A

We perform a density-matrix renormalization group (DMRG) study of the S=1/2 Heisenberg antiferromagnet on the kagome lattice to identify the conjectured spin liquid ground state. Exploiting SU(2) spin symmetry, which allows us to keep up to 16,000 DMRG states, we consider cylinders with circumferences up to 17 lattice spacings and find a spin liquid ground state with an estimated per site energy of -0.4386(5), a spin gap of 0.13(1), very short-range decay in spin, dimer and chiral correlation functions and finite topological entanglement gamma consistent with gamma=ld 2, ruling out gapless, chiral or non-topological spin liquids. At the same time, DMRG results provide strong evidence for a gapped topological Z2 spin liquid.

The problem of how complex quantum systems eventually come to rest lies at the heart of statistical mechanics. The maximum entropy principle put forward in 1957 by E. T. Jaynes suggests what quantum states one should expect in equilibrium but does not hint as to how closed quantum many-body systems dynamically equilibrate. A number of theoretical and numerical studies accumulate evidence that under specific conditions quantum many-body models can relax to a situation that locally or with respect to certain observables appears as if the entire system had relaxed to a maximum entropy state. In this work, we report the experimental observation of the non-equilibrium dynamics of a density wave of ultracold bosonic atoms in an optical lattice in the regime of strong correlations. Using an optical superlattice, we are able to prepare the system in a well-known initial state with high fidelity. We then follow the dynamical evolution of the system in terms of quasi-local densities, currents, and coherences. Numerical studies based on the time-dependent density-matrix renormalization group method are in an excellent quantitative agreement with the experimental data. For very long times, all three local observables show a fast relaxation to equilibrium values compatible with those expected for a global maximum entropy state. We find this relaxation of the quasi-local densities and currents to initially follow a power-law with an exponent being significantly larger than for free or hardcore bosons. For intermediate times the system fulfills the promise of being a dynamical quantum simulator, in that the controlled dynamics runs for longer times than present classical algorithms based on matrix product states can efficiently keep track of., 8 pages, 6 figures

We study the real-time dynamics of the local energy density in the spin-1/2 XXZ chain starting from initial states with an inhomogeneous profile of bond energies. Numerical simulations of the dynamics of the initial states are carried out using the adaptive time-dependent density matrix renormalization group method. We analyze the time dependence of the spatial variance associated with the local energy density to classify the dynamics as either ballistic or diffusive. Our results are consistent with ballistic behavior both in the massless and the massive phase. We also study the same problem within Luttinger Liquid theory and obtain that energy wave-packets propagate with the sound velocity. We recover this behavior in our numerical simulations in the limit of very weakly perturbed initial states., Comment: 15 pages, 21 figures, version as published, minor stylistic changes

We study the effects of the Coulomb interaction in the one dimensional Kondo lattice model on the phase diagram, the static magnetic susceptibility and electron spin relaxation. We show that onsite Coulomb interaction supports ferromagnetic order and nearest neighbor Coulomb interaction drives, depending on the electron filling, either a paramagnetic or ferromagnetic order. Furthermore we calculate electron quasiparticle life times, which can be related to electron spin relaxation and decoherence times, and explain their dependence on the strength of interactions and the electron filling in order to find the sweet spot of parameters where the relaxation time is maximized. We find that effective exchange processes between the electrons dominate the spin relaxation and decoherence rate., Comment: 9 pages, 7 figures