Showing 1–50 of 115 results for author: McCulloch, I P

Fractals and spontaneous symmetry breaking with type-B Goldstone modes: a perspective from entanglement

Authors: Huan-Qiang Zhou, Qian-Qian Shi, John O. Fjærestad, Ian P. McCulloch

Abstract: The one-dimensional spin-$s$ ${\rm SU}(2)$ ferromagnetic Heisenberg model, as a paradigmatic example for spontaneous symmetry breaking (SSB) with type-B Goldstone modes (GMs), is expected to exhibit an abstract fractal underlying the ground state subspace. This intrinsic abstract fractal is here revealed from a systematic investigation into the entanglement entropy for a linear combination of fact… ▽ More The one-dimensional spin-$s$ ${\rm SU}(2)$ ferromagnetic Heisenberg model, as a paradigmatic example for spontaneous symmetry breaking (SSB) with type-B Goldstone modes (GMs), is expected to exhibit an abstract fractal underlying the ground state subspace. This intrinsic abstract fractal is here revealed from a systematic investigation into the entanglement entropy for a linear combination of factorized (unentangled) ground states on a fractal decomposable into a set of the Cantor sets. The entanglement entropy scales logarithmically with the block size, with the prefactor being half the fractal dimension of a fractal, as long as the norm for the linear combination scales as the square root of the number of the self-similar building blocks kept at each step $k$ for a fractal, under an assumption that the maximum absolute value of the coefficients in the linear combination is chosen to be around one, and the coefficients in the linear combination are almost constants within the building blocks. Actually, the set of the fractal dimensions for all the Cantor sets forms a {\it dense} subset in the interval $[0,1]$. As a consequence, the ground state subspace is separated into a disjoint union of countably infinitely many regions, each of which is labeled by a decomposable fractal. Hence, the interpretation of the prefactor as half the fractal dimension is valid for any support beyond a fractal, which in turn leads to the identification of the fractal dimension with the number of type-B GMs for the orthonormal basis states. Our argument may be extended to any quantum many-body systems undergoing SSB with type-B GMs. △ Less

Submitted 25 July, 2024; originally announced July 2024.

Comments: 23 pages, 8 figures

Meson Mass Sets Onset Time of Anomalous Dynamical Quantum Phase Transitions

Authors: Jesse J. Osborne, Johannes Knaute, Ian P. McCulloch, Jad C. Halimeh

Abstract: Dynamical quantum phase transitions (DQPTs) have been established as a rigorous framework for investigating far-from-equilibrium quantum many-body criticality. Although initially thought to be trivially connected to an order parameter flipping sign, a certain kind of \textit{anomalous} DQPTs have been discovered that exhibit no direct connection to the order parameter and have been shown to arise… ▽ More Dynamical quantum phase transitions (DQPTs) have been established as a rigorous framework for investigating far-from-equilibrium quantum many-body criticality. Although initially thought to be trivially connected to an order parameter flipping sign, a certain kind of \textit{anomalous} DQPTs have been discovered that exhibit no direct connection to the order parameter and have been shown to arise in the presence of confinement. Here, we show in two paradigmatic models how the onset time of anomalous DQPTs is directly connected, through a power law, to the meson mass in the confined regime of a global symmetry-broken phase. This relation becomes more prominent the closer the initial parameters are to the equilibrium quantum critical point, where a relativistic quantum field theory emerges. Our findings draw a direct connection between mesons and anomalous DQPTs, highlighting the power of the latter to classify exotic far-from-equilibrium criticality. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Comments: $7+2$ pages, $4+4$ figures

Exact matrix product state representations for a type of scale-invariant states

Authors: Huan-Qiang Zhou, Qian-Qian Shi, Ian P. McCulloch

Abstract: Exact matrix product state representations for a type of scale-invariant states are presented, which describe highly degenerate ground states arising from spontaneous symmetry breaking with type-B Goldstone modes in one-dimensional quantum many-body systems. As a possible application, such a representation offers a convenient but powerful means for evaluating the norms of highly degenerate ground… ▽ More Exact matrix product state representations for a type of scale-invariant states are presented, which describe highly degenerate ground states arising from spontaneous symmetry breaking with type-B Goldstone modes in one-dimensional quantum many-body systems. As a possible application, such a representation offers a convenient but powerful means for evaluating the norms of highly degenerate ground states. This in turn allows us to perform a universal finite system-size scaling analysis of the entanglement entropy. Moreover, this approach vividly explains why the entanglement entropy does not depend on what types of the boundary conditions are adopted, either periodic boundary conditions or open boundary conditions. Illustrative examples include the ${\rm SU}(2)$ spin-$s$ Heisenberg ferromagnetic model, the ${\rm SU}(2s+1)$ ferromagnetic model, and the staggered ${\rm SU}(3)$ spin-1 ferromagnetic biquadratic model. △ Less

Submitted 14 March, 2024; originally announced March 2024.

Comments: 12 pages, 5 figures, 3 tables

Quantum Many-Body Scarring in $2+1$D Gauge Theories with Dynamical Matter

Authors: Jesse Osborne, Ian P. McCulloch, Jad C. Halimeh

Abstract: Quantum many-body scarring (QMBS) has emerged as an intriguing paradigm of weak ergodicity breaking in nonintegrable quantum many-body models, particularly lattice gauge theories (LGTs) in $1+1$ spacetime dimensions. However, an open question is whether QMBS exists in higher-dimensional LGTs with dynamical matter. Given that nonergodic dynamics in $d{=}1$ spatial dimension tend to vanish in… ▽ More Quantum many-body scarring (QMBS) has emerged as an intriguing paradigm of weak ergodicity breaking in nonintegrable quantum many-body models, particularly lattice gauge theories (LGTs) in $1+1$ spacetime dimensions. However, an open question is whether QMBS exists in higher-dimensional LGTs with dynamical matter. Given that nonergodic dynamics in $d{=}1$ spatial dimension tend to vanish in $d{>}1$, it is important to probe this question. Using matrix product state techniques for both finite and infinite systems, we show that QMBS occurs in the $2{+}1$D $\mathrm{U}(1)$ quantum link model (QLM), as evidenced in persistent coherent oscillations in local observables, a marked slowdown in the growth of the bipartite entanglement entropy, and revivals in the fidelity. Interestingly, we see that QMBS is more robust when the matter degrees of freedom are bosonic rather than fermionic. Our results further shed light on the intimate connection between gauge invariance and QMBS, and highlight the persistence of scarring in higher spatial dimensions. Our findings can be tested in near-term analog and digital quantum simulators, and we demonstrate their accessibility on a recently proposed cold-atom analog quantum simulator. △ Less

Submitted 28 March, 2024; v1 submitted 13 March, 2024; originally announced March 2024.

Comments: $7+2$ pages, $4+4$ figures. v2 includes longer evolution times and updated references

Entanglement and logarithmic spirals in a quantum spin-1 many-body system with competing dimer and trimer interactions

Authors: Huan-Qiang Zhou, Qian-Qian Shi, Ian P. McCulloch, Murray T. Batchelor

Abstract: Spontaneous symmetry breaking (SSB) with type-B Goldstone modes is investigated in the macroscopically degenerate phase for a quantum spin-1 many-body system with competing dimer and trimer interactions. The SSB involves three distinct patterns. The first occurs at the dimer point, with the pattern from staggered ${\rm SU}(3)$ to ${\rm U}(1)\times{\rm U}(1)$. The second occurs at the trimer point,… ▽ More Spontaneous symmetry breaking (SSB) with type-B Goldstone modes is investigated in the macroscopically degenerate phase for a quantum spin-1 many-body system with competing dimer and trimer interactions. The SSB involves three distinct patterns. The first occurs at the dimer point, with the pattern from staggered ${\rm SU}(3)$ to ${\rm U}(1)\times{\rm U}(1)$. The second occurs at the trimer point, with the pattern from uniform ${\rm SU}(3)$ to ${\rm U}(1)\times{\rm U}(1)$. The third occurs in the dimer-trimer regime, with the pattern from uniform ${\rm SU}(2)$ to ${\rm U}(1)$. The number of type-B Goldstone modes is thus two, two and one for the three patterns, respectively. The ground state degeneracies arising from the three patterns are exponential with the system size, which may be recognized as sequences of integers relevant to self-similar logarithmic spirals. This in turn is attributed to the presence of an emergent symmetry operation tailored to a specific degenerate ground state. As a consequence, the residual entropy is non-zero, which measures the disorder present in a unit cell of highly degenerate ground state generated from a generalized highest weight state. An exact Schmidt decomposition exists for the highly degenerate ground states, thus exposing the self-similarities underlying an abstract fractal, described by the fractal dimension. The latter is extracted from performing a universal finite system-size scaling analysis of the entanglement entropy, which is identical to the number of type-B Goldstone modes. The model under investigation thus accommodates an exotic scale invariant quantum state of matter. △ Less

Submitted 11 March, 2024; originally announced March 2024.

Comments: 14 pages, 3 figures, 2 tables

Comment on "Controlled Bond Expansion for Density Matrix Renormalization Group Ground State Search at Single-Site Costs" (Extended Version)

Authors: Ian P McCulloch, Jesse J Osborne

Abstract: In a recent Letter [Phys. Rev. Lett. 130, 246402 (2023)], Gleis, Li, and von Delft present an algorithm for expanding the bond dimension of a Matrix Product State wave function, giving accuracy similar to 2-site DMRG, but computationally more efficient, closer to the performance of 1-site DMRG. The Controlled Bond Expansion (CBE) algorithm uses the Hamiltonian projected onto two sites, and then fu… ▽ More In a recent Letter [Phys. Rev. Lett. 130, 246402 (2023)], Gleis, Li, and von Delft present an algorithm for expanding the bond dimension of a Matrix Product State wave function, giving accuracy similar to 2-site DMRG, but computationally more efficient, closer to the performance of 1-site DMRG. The Controlled Bond Expansion (CBE) algorithm uses the Hamiltonian projected onto two sites, and then further projected onto the two-site tangent space, to extract a set of $k$ vectors that are used to expand the basis between the two sites. CBE achieves this with a complicated sequence of five singular value decompositions (SVDs), in order to project onto the 2-site tangent space and reduce the bond dimension of the tensor network such that the contraction can be done in time $O(dwD^3)$. In this Comment, we show that (1) the projection onto the 2-site tangent space is unnecessary, and is generally not helpful; (2) the sequence of 5 SVDs can be replaced by a single $QR$ decomposition (optionally with one SVD as well), making use of the randomized SVD (RSVD) with high accuracy and significantly improved efficiency, scaling as $O(dwkD^2)$ i.e., the most expensive operations are only quadratic in the bond dimension $D$ and linear in the number of expansion vectors $k$; (3) several statements about the variational properties of the CBE algorithm are incorrect, and the variational properties are essentially identical to existing algorithms including 2-site DMRG and single-site subspace expansion (3S); (4) a similar RSVD approach can be applied to the 3S algorithm, which leads to many advantages over CBE, especially in systems with long range interactions. We also make some comments on the benchmarking MPS algorithms, and the overall computational efficiency with respect to the accuracy of the calculation. △ Less

Submitted 15 April, 2024; v1 submitted 1 March, 2024; originally announced March 2024.

Comments: Added some references; added appendix C

Probing Confinement Through Dynamical Quantum Phase Transitions: From Quantum Spin Models to Lattice Gauge Theories

Authors: Jesse Osborne, Ian P. McCulloch, Jad C. Halimeh

Abstract: Confinement is an intriguing phenomenon prevalent in condensed matter and high-energy physics. Exploring its effect on the far-from-equilibrium criticality of quantum many-body systems is of great interest both from a fundamental and technological point of view. Here, we employ large-scale uniform matrix product state calculations to show that a qualitative change in the type of dynamical quantum… ▽ More Confinement is an intriguing phenomenon prevalent in condensed matter and high-energy physics. Exploring its effect on the far-from-equilibrium criticality of quantum many-body systems is of great interest both from a fundamental and technological point of view. Here, we employ large-scale uniform matrix product state calculations to show that a qualitative change in the type of dynamical quantum phase transitions (DQPTs) accompanies the confinement-deconfinement transition in three paradigmatic models -- the power-law interacting quantum Ising chain, the two-dimensional quantum Ising model, and the spin-$S$ $\mathrm{U}(1)$ quantum link model. By tuning a confining parameter in these models, it is found that \textit{branch} (\textit{manifold}) DQPTs arise as a signature of (de)confinement. Whereas manifold DQPTs are associated with a sign change of the order parameter, their branch counterparts are not, but rather occur even when the order parameter exhibits considerably constrained dynamics. Our conclusions can be tested in modern quantum-simulation platforms, such as ion-trap setups and cold-atom experiments of gauge theories. △ Less

Submitted 18 October, 2023; originally announced October 2023.

Comments: 8 pages, 3 figures

Type-B Goldstone modes and a logarithmic spiral in the staggered $\rm SU(4)$ ferromagnetic spin-orbital model

Authors: Qian-Qian Shi, Huan-Qiang Zhou, Ian P. McCulloch, Murray T. Batchelor

Abstract: It is found that the staggered $\rm SU(4)$ ferromagnetic spin-orbital model accommodates highly degenerate ground states arising from spontaneous symmetry breaking with type-B Goldstone modes. The spontaneous symmetry breaking pattern is ${\rm SU(4)} \rightarrow {\rm U(1)} \times {\rm U(1)} \times {\rm U(1)}$, with three type-B Goldstone modes. An abstract fractal underlies the ground state subspa… ▽ More It is found that the staggered $\rm SU(4)$ ferromagnetic spin-orbital model accommodates highly degenerate ground states arising from spontaneous symmetry breaking with type-B Goldstone modes. The spontaneous symmetry breaking pattern is ${\rm SU(4)} \rightarrow {\rm U(1)} \times {\rm U(1)} \times {\rm U(1)}$, with three type-B Goldstone modes. An abstract fractal underlies the ground state subspace with the fractal dimension identified with the number of type-B Goldstone modes. This connection is established by evaluating the entanglement entropy, which exhibits a universal finite system-size scaling behaviour. The latter in turn implies that the entanglement entropy scales logarithmically with the block size in the thermodynamic limit. In addition, the ground state degeneracies, depending on the boundary conditions adopted, constitute the two Fibonacci-Lucas sequences. In the limit of large system size their asymptotic forms become a self-similar logarithmic spiral. As a result, the model has a non-zero residual entropy $S_{\!r} = -2 \ln R $, where $R=(\! \sqrt{6}-\!\sqrt{2})/2$. △ Less

Submitted 10 September, 2023; originally announced September 2023.

Comments: 16 pages, 5 figures, 1 table. arXiv admin note: substantial text overlap with arXiv:2302.13126

Extracting the number of type-B Goldstone modes and the dynamical critical exponent for a type of scale-invariant states

Authors: Huan-Qiang Zhou, Yan-Wei Dai, Qian-Qian Shi, Ian P. McCulloch, Murray T. Batchelor

Abstract: A generic scheme is proposed to perform a finite-entanglement scaling analysis for scale-invariant states, which appear to be highly degenerate ground states arising from spontaneous symmetry breaking with type-B Goldstone modes. This allows us to extract the number of type-B Goldstone modes and the dynamical critical exponent, in combination with a finite block-size scaling analysis, from numeric… ▽ More A generic scheme is proposed to perform a finite-entanglement scaling analysis for scale-invariant states, which appear to be highly degenerate ground states arising from spontaneous symmetry breaking with type-B Goldstone modes. This allows us to extract the number of type-B Goldstone modes and the dynamical critical exponent, in combination with a finite block-size scaling analysis, from numerical simulations of quantum many-body systems in the context of tensor network representations. The number of type-B Goldstone modes is identical to the fractal dimension, thus reflecting an abstract fractal underlying the ground state subspace. As illustrative examples, we investigate the spin-$s$ Heisenberg ferromagnetic model, the $\rm{SU}(3)$ ferromagnetic model and the $\rm{SO}(4)$ spin-orbital model. △ Less

Submitted 30 November, 2023; v1 submitted 10 September, 2023; originally announced September 2023.

Comments: 14 pages, 24 figures, 11 tables. Minor changes

Wavefunction branching: when you can't tell pure states from mixed states

Authors: Jordan K. Taylor, Ian P. McCulloch

Abstract: We propose a definition of wavefunction "branchings": quantum superpositions which can't be feasibly distinguished from the corresponding mixed state, even under time evolution. Our definition is largely independent of interpretations, requiring only that it takes many more local gates to swap branches than to distinguish them. We give several examples of states admitting such branch decomposition… ▽ More We propose a definition of wavefunction "branchings": quantum superpositions which can't be feasibly distinguished from the corresponding mixed state, even under time evolution. Our definition is largely independent of interpretations, requiring only that it takes many more local gates to swap branches than to distinguish them. We give several examples of states admitting such branch decompositions. Under our definition, we show that attempts to get relative-phase information between branches will fail without frequent active error correction, that branches are effectively the opposite of good error-correcting codes, that branches effectively only grow further apart in time under natural evolution, that branches tend to absorb spatial entanglement, that branching is stronger in the presence of conserved quantities, and that branching implies effective irreversibility. Identifying these branch decompositions in many-body quantum states could shed light on the emergence of classicality, provide a metric for experimental tests at the quantum/ classical boundary, and allow for longer numerical time evolution simulations. We see this work as a generalization of the basic ideas of environmentally-induced decoherence to situations with no clear system/ environment split. △ Less

Submitted 29 August, 2023; v1 submitted 8 August, 2023; originally announced August 2023.

Comments: 18 pages, 8 figures; Added more references on quantum error correction with symmetries

Spin-$S$ $\mathrm{U}(1)$ Quantum Link Models with Dynamical Matter on a Quantum Simulator

Authors: Jesse Osborne, Bing Yang, Ian P. McCulloch, Philipp Hauke, Jad C. Halimeh

Abstract: 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… ▽ More 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. △ Less

Submitted 10 May, 2023; originally announced May 2023.

Comments: $15$ pages, $8$ figures

Entanglement entropy for scale-invariant states: universal finite-size scaling

Authors: Huan-Qiang Zhou, Qian-Qian Shi, Ian P. McCulloch, Murray T. Batchelor

Abstract: 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 const… ▽ More 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. △ Less

Submitted 22 April, 2023; originally announced April 2023.

Comments: 8 pages, 3 figures

Goldstone modes and the golden spiral in the ferromagnetic spin-1 biquadratic model

Authors: Huan-Qiang Zhou, Qian-Qian Shi, Ian P. McCulloch, Murray T. Batchelor

Abstract: 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 sponta… ▽ More 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. △ Less

Submitted 30 November, 2023; v1 submitted 25 February, 2023; originally announced February 2023.

Comments: 20 pages, 3 figures,1 table, additional references to earlier related work added, a minor change

Disorder-Free Localization in $2+1$D Lattice Gauge Theories with Dynamical Matter

Authors: Jesse Osborne, Ian P. McCulloch, Jad C. Halimeh

Abstract: 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… ▽ More 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. △ Less

Submitted 18 January, 2023; originally announced January 2023.

Comments: $7$ pages, $4$ figures

Ergodicity Breaking Under Confinement in Cold-Atom Quantum Simulators

Authors: Jean-Yves Desaules, Guo-Xian Su, Ian P. McCulloch, Bing Yang, Zlatko Papić, Jad C. Halimeh

Abstract: 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 $θ$-angle, which can be used to tune a confinement-deconfinement transiti… ▽ More 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 $θ$-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 $μ$ and confining potential $χ$, an emergent integrable phase for large $μ$, 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. △ Less

Submitted 24 February, 2024; v1 submitted 18 January, 2023; originally announced January 2023.

Comments: 28+12 pages, 16+8 figures

Journal ref: Quantum 8, 1274 (2024)

Deep Learning of Phase Transitions for Quantum Spin Chains from Correlation Aspects

Authors: Ming-Chiang Chung, Guang-Yu Huang, Ian P. McCulloch, Yuan-Hong Tsai

Abstract: 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… ▽ More 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. △ Less

Submitted 9 May, 2023; v1 submitted 16 January, 2023; originally announced January 2023.

Comments: 18 pages, 21 figures

Large-Scale $2+1$D $\mathrm{U}(1)$ Gauge Theory with Dynamical Matter in a Cold-Atom Quantum Simulator

Authors: Jesse Osborne, Ian P. McCulloch, Bing Yang, Philipp Hauke, Jad C. Halimeh

Abstract: 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 map… ▽ More 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. △ Less

Submitted 10 November, 2022; v1 submitted 2 November, 2022; originally announced November 2022.

Comments: 8+5 pages, 3+2 figures

Tuning the Topological $θ$-Angle in Cold-Atom Quantum Simulators of Gauge Theories

Authors: Jad C. Halimeh, Ian P. McCulloch, Bing Yang, Philipp Hauke

Abstract: The topological $θ$-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 o… ▽ More The topological $θ$-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 $θ$-angle. Here, we demonstrate how a tunable topological $θ$-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 $θ$-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. △ Less

Submitted 15 April, 2022; v1 submitted 13 April, 2022; originally announced April 2022.

Comments: 20 pages, 11 figures, journal article

Journal ref: PRX Quantum 3, 040316 (2022)

Absence of a critical nematic phase in the vicinity of the $\rm {SU}(3)$ ferromagnetic point for the one-dimensional spin-1 bilinear-biquadratic model

Authors: Yan-Wei Dai, Qian-Qian Shi, Huan-Qiang Zhou, Ian P. McCulloch

Abstract: The absence of a critical nematic phase in the vicinity of the $\rm {SU}(3)$ ferromagnetic point for the one-dimensional spin-1 bilinear-biquadratic model is demonstrated by means of the tensor network algorithms. As it turns out, the phase transition from the ferromagnetic phase to the dimerized phase at the $\rm {SU}(3)$ ferromagnetic point is direct, but not of the first-order. The transition p… ▽ More The absence of a critical nematic phase in the vicinity of the $\rm {SU}(3)$ ferromagnetic point for the one-dimensional spin-1 bilinear-biquadratic model is demonstrated by means of the tensor network algorithms. As it turns out, the phase transition from the ferromagnetic phase to the dimerized phase at the $\rm {SU}(3)$ ferromagnetic point is direct, but not of the first-order. The transition point features highly degenerate ground states, which are scale but not conformally invariant, with the fractal dimension being equal to 2. The conceptual developments in effective field theories - the fractal dimension and the counting rule of the Goldstone modes - play a pivotal role in clarifying the numerical artifacts arising from the finiteness of the bond dimension in the tensor network simulations, which are attributed to a proximity effect to a highly entangled scale or conformally invariant ground state. △ Less

Submitted 30 November, 2023; v1 submitted 4 January, 2022; originally announced January 2022.

Comments: 15 pages, 21 figures, 7 tables. References updated

Detecting quantum phase transitions in the quasi-stationary regime of Ising chains

Authors: Ceren B. Dağ, Philipp Uhrich, Yidan Wang, Ian P. McCulloch, Jad C. Halimeh

Abstract: 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 analytical… ▽ More 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. △ Less

Submitted 20 September, 2022; v1 submitted 6 October, 2021; originally announced October 2021.

Comments: 22 pages, 23 figures. Revised with new analytical calculations, new content and improved presentation

Journal ref: Phys. Rev. B 107, 094432 (2023)

What is a quantum shock wave?

Authors: S. A. Simmons, F. A. Bayocboc, Jr., J. C. Pillay, D. Colas, I. P. McCulloch, K. V. Kheruntsyan

Abstract: Shock waves are examples of the far-from-equilibrium behaviour of matter; they are ubiquitous in nature, yet the underlying microscopic mechanisms behind their formation are not well understood. Here, we study the dynamics of dispersive quantum shock waves in a one-dimensional Bose gas, and show that the oscillatory train forming from a local density bump expanding into a uniform background is a r… ▽ More Shock waves are examples of the far-from-equilibrium behaviour of matter; they are ubiquitous in nature, yet the underlying microscopic mechanisms behind their formation are not well understood. Here, we study the dynamics of dispersive quantum shock waves in a one-dimensional Bose gas, and show that the oscillatory train forming from a local density bump expanding into a uniform background is a result of quantum mechanical self-interference. The amplitude of oscillations, i.e., the interference contrast, decreases with the increase of both the temperature of the gas and the interaction strength due to the reduced phase coherence length. Furthermore, we show that vacuum and thermal fluctuations can significantly wash out the interference contrast, seen in the mean-field approaches, due to shot-to-shot fluctuations in the position of interference fringes around the mean. △ Less

Submitted 3 November, 2020; v1 submitted 27 June, 2020; originally announced June 2020.

Comments: Final published version, 6 pages, 3 figures, plus Supplementary Material

Journal ref: Phys. Rev. Lett. 125, 180401 (2020)

Robustness of gauge-invariant dynamics against defects in ultracold-atom gauge theories

Authors: Jad C. Halimeh, Robert Ott, Ian P. McCulloch, Bing Yang, Philipp Hauke

Abstract: Recent years have seen strong progress in quantum simulation of gauge-theory dynamics using ultracold-atom experiments. A principal challenge in these efforts is the certification of gauge invariance, which has recently been realized in [B.~Yang et al., arXiv:2003.08945]. One major but poorly investigated experimental source of gauge-invariance violation is an imperfect preparation of the initial… ▽ More Recent years have seen strong progress in quantum simulation of gauge-theory dynamics using ultracold-atom experiments. A principal challenge in these efforts is the certification of gauge invariance, which has recently been realized in [B.~Yang et al., arXiv:2003.08945]. One major but poorly investigated experimental source of gauge-invariance violation is an imperfect preparation of the initial state. Using the time-dependent density-matrix renormalization group, we analyze the robustness of gauge-invariant dynamics against potential preparation defects in the above ultracold-atom implementation of a $\mathrm{U}(1)$ gauge theory. We find defects related to an erroneous initialization of matter fields to be innocuous, as the associated gauge-invariance violation remains strongly localized throughout the time evolution. A defect due to faulty initialization of the gauge field leads to a mild proliferation of the associated violation. Furthermore, we characterize the influence of immobile and mobile defects by monitoring the spread of entanglement entropy. Overall, our results indicate that the aforementioned experimental realization exhibits a high level of fidelity in the gauge invariance of its dynamics at all evolution times. Our work provides strong evidence that ultracold-atom setups can serve as an extremely reliable framework for the quantum simulation of gauge-theory dynamics. △ Less

Submitted 11 August, 2020; v1 submitted 20 May, 2020; originally announced May 2020.

Comments: Accepted version. 20 pages, 11 figures. Videos of dynamics at https://www.youtube.com/watch?v=_FCrzldBtz8&list=PLoUsb3eaKix6xYrnNk68Bh3imd3II54gO

Journal ref: Phys. Rev. Research 2, 033361 (2020)

Hybrid infinite time-evolving block decimation algorithm for long-range multi-dimensional quantum many-body systems

Authors: Tomohiro Hashizume, Jad C. Halimeh, Ian P. McCulloch

Abstract: In recent years, the infinite time-evolution block decimation (iTEBD) method has been demonstrated to be one of the most efficient and powerful numerical schemes for time-evolution in one-dimensional quantum many-body systems. However, a major shortcoming of the method, along with other state-of-the-art algorithms for many-body dynamics, has been their restriction to one spatial dimension. We pres… ▽ More In recent years, the infinite time-evolution block decimation (iTEBD) method has been demonstrated to be one of the most efficient and powerful numerical schemes for time-evolution in one-dimensional quantum many-body systems. However, a major shortcoming of the method, along with other state-of-the-art algorithms for many-body dynamics, has been their restriction to one spatial dimension. We present an algorithm based on a \textit{hybrid} extension of iTEBD where finite blocks of a chain are first locally time-evolved before an iTEBD-like method combines these processes globally. This in turn permits simulating the dynamics of many-body systems in the thermodynamic limit in $d\geq1$ dimensions including in the presence of long-range interactions. Our work paves the way for simulating the dynamics of many-body phenomena that occur exclusively in higher dimensions, and whose numerical treatments have hitherto been limited to exact diagonalization of small systems, which fundamentally limits a proper investigation of dynamical criticality. We expect the algorithm presented here to be of significant importance to validating and guiding investigations in state-of-the-art ion-trap and ultracold-atom experiments. △ Less

Submitted 23 October, 2019; originally announced October 2019.

Comments: 10 pages, 6 figures, journal article

Journal ref: Phys. Rev. B 102, 035115 (2020)

Cumulants and Scaling Functions of Infinite Matrix Product States

Authors: Jason C. Pillay, Ian P. McCulloch

Abstract: The order parameter cumulants of infinite matrix product ground states are evaluated across a quantum phase transition. A scheme using the Binder cumulant, finite-entanglement scaling and scaling functions to obtain the critical point and exponents of the correlation length and cumulants is presented. Analogous to the scaling relations that relate the exponents of various thermodynamic quantities,… ▽ More The order parameter cumulants of infinite matrix product ground states are evaluated across a quantum phase transition. A scheme using the Binder cumulant, finite-entanglement scaling and scaling functions to obtain the critical point and exponents of the correlation length and cumulants is presented. Analogous to the scaling relations that relate the exponents of various thermodynamic quantities, a cumulant exponent relation is derived and used to check the consistency and relationship between the cumulant exponents. This scheme gives a numerically economical way of accurately obtaining the critical exponents. Examples of this scheme are shown for four one-dimensional models - the transverse field Ising model, the topological Kondo insulator, the S = 1 Heisenberg chain with single-ion anisotropy and the Bose-Hubbard model. A two-dimensional model is also exemplified in the square lattice transverse field Ising model on an infinite cylinder. These exemplary systems portray a variety of local and string order parameters as well as phase transition classes that can be studied with the scaling functions and infinite matrix product states. △ Less

Submitted 29 December, 2019; v1 submitted 10 June, 2019; originally announced June 2019.

Journal ref: Phys. Rev. B 100, 235140 (2019)

Dynamical phase transitions in the two-dimensional transverse-field Ising model

Authors: Tomohiro Hashizume, Ian P. McCulloch, Jad C. Halimeh

Abstract: We investigate two separate notions of dynamical phase transitions in the two-dimensional nearest-neighbor transverse-field Ising model on a square lattice using matrix product states and a new \textit{hybrid} infinite time-evolving block decimation algorithm, where the model is implemented on an infinitely long cylinder with a finite diameter along which periodic boundary conditions are employed.… ▽ More We investigate two separate notions of dynamical phase transitions in the two-dimensional nearest-neighbor transverse-field Ising model on a square lattice using matrix product states and a new \textit{hybrid} infinite time-evolving block decimation algorithm, where the model is implemented on an infinitely long cylinder with a finite diameter along which periodic boundary conditions are employed. Starting in an ordered initial state, our numerical results suggest that quenches below the dynamical critical point give rise to a ferromagnetic long-time steady state with the Loschmidt return rate exhibiting \textit{anomalous} cusps even when the order parameter never crosses zero. Within the accessible timescales of our numerics, quenches above the dynamical critical point suggest a paramagnetic long-time steady state with the return rate exhibiting \textit{regular} cusps connected to zero crossings of the order parameter. Additionally, our simulations indicate that quenching slightly above the dynamical critical point leads to a coexistence region where both anomalous and regular cusps appear in the return rate. Quenches from the disordered phase further confirm our main conclusions. Our work supports the recent finding that anomalous cusps arise only when local spin excitations are the energetically dominant quasiparticles. Our results are accessible in modern Rydberg experiments. △ Less

Submitted 16 January, 2022; v1 submitted 22 November, 2018; originally announced November 2018.

Comments: 9 pages, 6 figures, journal article

Journal ref: Phys. Rev. Research 4, 013250 (2022)

Topological nature of spinons and holons: Elementary excitations from matrix product states with conserved symmetries

Authors: V. Zauner-Stauber, L. Vanderstraeten, J. Haegeman, I. P. McCulloch, F. Verstraete

Abstract: We develop variational matrix product state (MPS) methods with symmetries to determine dispersion relations of one dimensional quantum lattices as a function of momentum and preset quantum number. We test our methods on the XXZ spin chain, the Hubbard model and a non-integrable extended Hubbard model, and determine the excitation spectra with a precision similar to the one of the ground state. The… ▽ More We develop variational matrix product state (MPS) methods with symmetries to determine dispersion relations of one dimensional quantum lattices as a function of momentum and preset quantum number. We test our methods on the XXZ spin chain, the Hubbard model and a non-integrable extended Hubbard model, and determine the excitation spectra with a precision similar to the one of the ground state. The formulation in terms of quantum numbers makes the topological nature of spinons and holons very explicit. In addition, the method also enables an easy and efficient direct calculation of the necessary magnetic field or chemical potential required for a certain ground state magnetization or particle density. △ Less

Submitted 19 April, 2019; v1 submitted 20 February, 2018; originally announced February 2018.

Comments: 13 pages, 4 pages appendix, 8 figures

Journal ref: Phys. Rev. B 97, 235155 (2018)

Phase diagram of the quantum Ising model with long-range interactions on an infinite-cylinder triangular lattice

Authors: S. N. Saadatmand, S. D. Bartlett, I. P. McCulloch

Abstract: Obtaining quantitative ground-state behavior for geometrically-frustrated quantum magnets with long-range interactions is challenging for numerical methods. Here, we demonstrate that the ground states of these systems on two-dimensional lattices can be efficiently obtained using state-of-the-art translation-invariant variants of matrix product states and density-matrix renormalization-group algori… ▽ More Obtaining quantitative ground-state behavior for geometrically-frustrated quantum magnets with long-range interactions is challenging for numerical methods. Here, we demonstrate that the ground states of these systems on two-dimensional lattices can be efficiently obtained using state-of-the-art translation-invariant variants of matrix product states and density-matrix renormalization-group algorithms. We use these methods to calculate the fully-quantitative ground-state phase diagram of the long-range interacting triangular Ising model with a transverse field on 6-leg infinite-length cylinders, and scrutinize the properties of the detected phases. We compare these results with those of the corresponding nearest neighbor model. Our results suggest that, for such long-range Hamiltonians, the long-range quantum fluctuations always lead to long-range correlations, where correlators exhibit power-law decays instead of the conventional exponential drops observed for short-range correlated gapped phases. Our results are relevant for comparisons with recent ion-trap quantum simulator experiments that demonstrate highly-controllable long-range spin couplings for several hundred ions. △ Less

Submitted 2 April, 2019; v1 submitted 1 February, 2018; originally announced February 2018.

Comments: 16 pages, 12 figures, and no supplemental materials. v4: equivalent to the published version, except a few minor typographical and notation corrections

Journal ref: Phys. Rev. B 97, 155116 (2018)

Dynamics of multiple atoms in one-dimensional fields

Authors: Carlo Cascio, Jad C. Halimeh, Ian P. McCulloch, Alessio Recati, Ines de Vega

Abstract: We analyze the dynamics of a set of two-level atoms coupled to the electromagnetic environment within a waveguide. This problem is often tackled by assuming a weak coupling between the atoms and the environment as well as the associated Markov approximation. We show that the accuracy of such an approximation may be more limited than in the single-atom case and also be strongly determined by the pr… ▽ More We analyze the dynamics of a set of two-level atoms coupled to the electromagnetic environment within a waveguide. This problem is often tackled by assuming a weak coupling between the atoms and the environment as well as the associated Markov approximation. We show that the accuracy of such an approximation may be more limited than in the single-atom case and also be strongly determined by the presence of collective effects produced by atom-atom interactions. To this aim, we solve the full problem with exact diagonalization and also the time-dependent density matrix renormalization group method, and compare the result to that obtained within a weak-coupling master equation and with the Dicke approximation. Finally, we study the dynamics of the entanglement within the system when considering several inter-atomic distances and atomic frequencies. △ Less

Submitted 23 January, 2019; v1 submitted 24 January, 2018; originally announced January 2018.

Journal ref: Phys. Rev. A 99, 013845 (2019)

Topological phase transition and the effect of Hubbard interaction on the one-dimensional topological Kondo insulator

Authors: Jason C. Pillay, Ian P. McCulloch

Abstract: The effect of a local Kondo coupling and Hubbard interaction on the topological phase of the one-dimensional topological Kondo insulator (TKI) is numerically investigated using the infinite matrix-product state density-matrix renormalization group algorithm. The groundstate of the TKI is a symmetry-protected topological (SPT) phase protected by inversion symmetry. It is found that on its own, the… ▽ More The effect of a local Kondo coupling and Hubbard interaction on the topological phase of the one-dimensional topological Kondo insulator (TKI) is numerically investigated using the infinite matrix-product state density-matrix renormalization group algorithm. The groundstate of the TKI is a symmetry-protected topological (SPT) phase protected by inversion symmetry. It is found that on its own, the Hubbard interaction that tends to force fermions into a one-charge per site order is insufficient to destroy the SPT phase. However when the local Kondo Hamiltonian term that favors a topologically trivial groundstate with a one-charge per site order is introduced, the Hubbard interaction assists in the destruction of the SPT phase. This topological phase transition occurs in the charge sector where the correlation length of the charge excitation diverges while the correlation length of the spin excitation remains finite. The critical exponents, central charge and the phase diagram separating the SPT phase from the topologically trivial phase are presented. △ Less

Submitted 27 May, 2018; v1 submitted 4 January, 2018; originally announced January 2018.

Comments: 15 pages, 22 figures

Journal ref: Phys. Rev. B 97, 205133 (2018)

Symmetry between repulsive and attractive interactions in driven-dissipative Bose-Hubbard systems

Authors: Adil A. Gangat, Ian P. McCulloch, Ying-Jer Kao

Abstract: The driven-dissipative Bose-Hubbard model can be experimentally realized with either negative or positive onsite detunings, inter-site hopping energies, and onsite interaction energies. Here we use one-dimensional matrix product density operators to perform a fully quantum investigation of the dependence of the non-equilibrium steady states of this model on the signs of these parameters. Due to a… ▽ More The driven-dissipative Bose-Hubbard model can be experimentally realized with either negative or positive onsite detunings, inter-site hopping energies, and onsite interaction energies. Here we use one-dimensional matrix product density operators to perform a fully quantum investigation of the dependence of the non-equilibrium steady states of this model on the signs of these parameters. Due to a symmetry in the Lindblad master equation, we find that simultaneously changing the sign of the interaction energies, hopping energies, and chemical potentials leaves the local boson number distribution and inter-site number correlations invariant, and the steady-state complex conjugated. This shows that all driven-dissipative phenomena of interacting bosons described by the Lindblad master equation, such as "fermionization" and "superbunching", can equivalently occur with attractive or repulsive interactions. △ Less

Submitted 20 October, 2017; v1 submitted 22 June, 2017; originally announced June 2017.

Comments: single column 12 pages, 4 figures, 1 table

Journal ref: Scientific Reports 8, 3698 (2018)

Detection and characterization of symmetry-broken long-range orders in the spin-$\frac{1}{2}$ triangular Heisenberg model

Authors: S. N. Saadatmand, I. P. McCulloch

Abstract: We present new numerical tools to analyze symmetry-broken phases in the context of $SU(2)$-symmetric translation-invariant matrix product states (MPS) and density-matrix renormalization-group (DMRG) methods for infinite cylinders, and determine the phase diagram of the geometrically-frustrated triangular Heisenberg model with nearest and next-nearest neighbor (NN and NNN) interactions. The appeara… ▽ More We present new numerical tools to analyze symmetry-broken phases in the context of $SU(2)$-symmetric translation-invariant matrix product states (MPS) and density-matrix renormalization-group (DMRG) methods for infinite cylinders, and determine the phase diagram of the geometrically-frustrated triangular Heisenberg model with nearest and next-nearest neighbor (NN and NNN) interactions. The appearance of Nambu-Goldstone modes in the excitation spectrum is characterized by "tower of states" levels in the momentum-resolved entanglement spectrum. Symmetry-breaking phase transitions are detected by a combination of the correlation lengths and second and fourth cumulants of the magnetic order parameters (which we call the Binder ratio), even though symmetry implies that the order parameter itself is strictly zero. Using this approach, we have identified $120^{\circ}$ order, a columnar order, and an algebraic spin liquid (specific to width-6 systems), alongside the previously studied topological spin liquid phase. For the latter, we also demonstrate robustness against chiral perturbations. △ Less

Submitted 2 April, 2019; v1 submitted 11 April, 2017; originally announced April 2017.

Comments: 25 pages, 19 figures, and no supplemental materials. v5: equivalent to the published version, except added a note that $J_2=0.125$ in Fig 19

Journal ref: Phys. Rev. B 96, 075117 (2017)

Precursor of Laughlin state of hard core bosons on a two leg ladder

Authors: Alexandru Petrescu, Marie Piraud, Guillaume Roux, I. P. McCulloch, Karyn Le Hur

Abstract: We study hard core bosons on a two leg ladder lattice under the orbital effect of a uniform magnetic field. At densities which are incommensurate with flux, the ground state is a Meissner state, or a vortex state, depending on the strength of the flux. When the density is commensurate with the flux, analytical arguments predict the existence of a ground state of central charge $c = 1$, which displ… ▽ More We study hard core bosons on a two leg ladder lattice under the orbital effect of a uniform magnetic field. At densities which are incommensurate with flux, the ground state is a Meissner state, or a vortex state, depending on the strength of the flux. When the density is commensurate with the flux, analytical arguments predict the existence of a ground state of central charge $c = 1$, which displays signatures compatible with the expected Laughlin state at $ν=1/2$. This differs from the coupled wire construction of the Laughlin state in that there exists a nonzero backscattering term in the edge Hamiltonian. We construct a phase diagram versus density and flux in order to delimit the region where this precursor to the Laughlin state is the ground state, by using a combination of bosonization and numerics based on the density matrix renormalization group (DMRG) and exact diagonalization. We obtain the phase diagram from local observables and central charge. We use bipartite charge fluctuations to deduce the Luttinger parameter for the edge Luttinger liquid corresponding to the Laughlin state. The properties studied with local observables are confirmed by the long distance behavior of correlation functions. Our findings are consistent with a calculation of the many body ground state transverse conductivity in a thin torus geometry for parameters corresponding to the Laughlin state. The model considered is simple enough such that the precursor to the Laughlin state could be realized in current ultracold atom, Josephson junction array, and quantum circuit experiments. △ Less

Submitted 20 February, 2017; v1 submitted 15 December, 2016; originally announced December 2016.

Comments: 22 pages, 11 figures; version submitted to Phys. Rev. B

Journal ref: Phys. Rev. B 96, 014524 (2017)

Generic Construction of Efficient Matrix Product Operators

Authors: C. Hubig, I. P. McCulloch, U. Schollwöck

Abstract: Matrix Product Operators (MPOs) are at the heart of the second-generation Density Matrix Renormalisation Group (DMRG) algorithm formulated in Matrix Product State language. We first summarise the widely known facts on MPO arithmetic and representations of single-site operators. Second, we introduce three compression methods (Rescaled SVD, Deparallelisation and Delinearisation) for MPOs and show th… ▽ More Matrix Product Operators (MPOs) are at the heart of the second-generation Density Matrix Renormalisation Group (DMRG) algorithm formulated in Matrix Product State language. We first summarise the widely known facts on MPO arithmetic and representations of single-site operators. Second, we introduce three compression methods (Rescaled SVD, Deparallelisation and Delinearisation) for MPOs and show that it is possible to construct efficient representations of arbitrary operators using MPO arithmetic and compression. As examples, we construct powers of a short-ranged spin-chain Hamiltonian, a complicated Hamiltonian of a two-dimensional system and, as proof of principle, the long-range four-body Hamiltonian from quantum chemistry. △ Less

Submitted 5 January, 2017; v1 submitted 8 November, 2016; originally announced November 2016.

Comments: 13 pages, 10 figures

Journal ref: Phys. Rev. B 95, 035129 (2017)

Symmetry-broken states in a system of interacting bosons on a two-leg ladder with a uniform Abelian gauge field

Authors: S. Greschner, M. Piraud, F. Heidrich-Meisner, I. P. McCulloch, U. Schollwöck, T. Vekua

Abstract: We study the quantum phases of bosons with repulsive contact interactions on a two-leg ladder in the presence of a uniform Abelian gauge field. The model realizes many interesting states, including Meissner phases, vortex-fluids, vortex-lattices, charge-density-waves and the biased-ladder phase. Our work focuses on the subset of these states that break a discrete symmetry. We use density matrix re… ▽ More We study the quantum phases of bosons with repulsive contact interactions on a two-leg ladder in the presence of a uniform Abelian gauge field. The model realizes many interesting states, including Meissner phases, vortex-fluids, vortex-lattices, charge-density-waves and the biased-ladder phase. Our work focuses on the subset of these states that break a discrete symmetry. We use density matrix renormalization group simulations to demonstrate the existence of three vortex-lattice states at different vortex densities and we characterize the phase transitions from these phases into neighboring states. Furthermore, we provide an intuitive explanation of the chiral-current reversal effect that is tied to some of these vortex lattices. We also study a charge-density-wave state that exists at 1/4 particle filling at large interaction strengths and flux values close to half a flux quantum. By changing the system parameters, this state can transition into a completely gapped vortex-lattice Mott-insulating state. We elucidate the stability of these phases against nearest-neighbor interactions on the rungs of the ladder relevant for experimental realizations with a synthetic lattice dimension. A charge-density-wave state at 1/3 particle filling can be stabilized for flux values close to half a flux-quantum and for very strong on-site interactions in the presence of strong repulsion on the rungs. Finally, we analytically describe the emergence of these phases in the low-density regime, and, in particular, we obtain the boundaries of the biased-ladder phase, i.e., the phase that features a density imbalance between the legs. We make contact to recent quantum-gas experiments that realized related models and discuss signatures of these quantum states in experimentally accessible observables. △ Less

Submitted 29 December, 2016; v1 submitted 7 October, 2016; originally announced October 2016.

Comments: 25 pages, 26 figures

Journal ref: Phys. Rev. A 94, 063628 (2016)

Prethermalization and Persistent Order in the Absence of a Thermal Phase Transition

Authors: Jad C. Halimeh, Valentin Zauner-Stauber, Ian P. McCulloch, Inés de Vega, Ulrich Schollwöck, Michael Kastner

Abstract: We numerically study the dynamics after a parameter quench in the one-dimensional transverse-field Ising model with long-range interactions ($\propto 1/r^α$ with distance $r$), for finite chains and also directly in the thermodynamic limit. In nonequilibrium, i.e., before the system settles into a thermal state, we find a long-lived regime that is characterized by a prethermal value of the magneti… ▽ More We numerically study the dynamics after a parameter quench in the one-dimensional transverse-field Ising model with long-range interactions ($\propto 1/r^α$ with distance $r$), for finite chains and also directly in the thermodynamic limit. In nonequilibrium, i.e., before the system settles into a thermal state, we find a long-lived regime that is characterized by a prethermal value of the magnetization, which in general differs from its thermal value. We find that the ferromagnetic phase is stabilized dynamically: as a function of the quench parameter, the prethermal magnetization shows a transition between a symmetry-broken and a symmetric phase, even for those values of $α$ for which no finite-temperature transition occurs in equilibrium. The dynamical critical point is shifted with respect to the equilibrium one, and the shift is found to depend on $α$ as well as on the quench parameters. △ Less

Submitted 5 October, 2016; originally announced October 2016.

Comments: 6 pages, 4 figures

Journal ref: Phys. Rev. B 95, 024302 (2017)

Haldane insulator protected by reflection symmetry in the doped Hubbard model on the three-legged ladder

Authors: H. L. Nourse, I. P. McCulloch, C. Janani, B. J. Powell

Abstract: We demonstrate the existence of an insulating phase in the three-legged Hubbard ladder at two-thirds filling. In this phase chargons are bound because the physics within a unit cell favors the formation of triplets. The resultant moments lead to a ground state in the Haldane phase, a symmetry protected topological state of matter. In this purely fermionic model, reflection is protecting but time r… ▽ More We demonstrate the existence of an insulating phase in the three-legged Hubbard ladder at two-thirds filling. In this phase chargons are bound because the physics within a unit cell favors the formation of triplets. The resultant moments lead to a ground state in the Haldane phase, a symmetry protected topological state of matter. In this purely fermionic model, reflection is protecting but time reversal and dihedral symmetries are not, in contrast to spin models. △ Less

Submitted 15 December, 2016; v1 submitted 14 June, 2016; originally announced June 2016.

Comments: 6 pages, 6 figures

Journal ref: Phys. Rev. B 94, 214418 (2016)

Symmetry fractionalization in the topological phase of the spin-$\frac{1}{2}$ $J_1$-$J_2$ triangular Heisenberg model

Authors: S. N. Saadatmand, I. P. McCulloch

Abstract: Using density-matrix renormalization-group calculations for infinite cylinders, we elucidate the properties of the spin-liquid phase of the spin-$\frac{1}{2}$ $J_1$-$J_2$ Heisenberg model on the triangular lattice. We find four distinct ground-states characteristic of a non-chiral, $Z_2$ topologically ordered state with vison and spinon excitations. We shed light on the interplay of topological or… ▽ More Using density-matrix renormalization-group calculations for infinite cylinders, we elucidate the properties of the spin-liquid phase of the spin-$\frac{1}{2}$ $J_1$-$J_2$ Heisenberg model on the triangular lattice. We find four distinct ground-states characteristic of a non-chiral, $Z_2$ topologically ordered state with vison and spinon excitations. We shed light on the interplay of topological ordering and global symmetries in the model by detecting fractionalization of time-reversal and space-group dihedral symmetries in the anyonic sectors, which leads to coexistence of symmetry protected and intrinsic topological order. The anyonic sectors, and information on the particle statistics, can be characterized by degeneracy patterns and symmetries of the entanglement spectrum. We demonstrate the ground-states on finite-width cylinders are short-range correlated and gapped; however some features in the entanglement spectrum suggest that the system develops gapless spinon-like edge excitations in the large-width limit. △ Less

Submitted 12 April, 2018; v1 submitted 1 June, 2016; originally announced June 2016.

Comments: 6 pages, 4 figures, 1 table, and 10 pages of supplemental materials. v3: equivalent to the published version

Journal ref: Phys. Rev. B 94, 121111 (2016)

An efficient perturbation theory of density matrix renormalization group

Authors: Emanuele Tirrito, Shi-Ju Ran, Andrew J. Ferris, Ian P. McCulloch, Maciej Lewenstein

Abstract: 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… ▽ More 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 | ψ_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 ψ_i | \hat{H} | ψ_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. △ Less

Submitted 3 August, 2016; v1 submitted 3 May, 2016; originally announced May 2016.

Journal ref: Phys. Rev. B 95, 064110 (2017)

Optical conductivity of the Hubbard chain away from half filling

Authors: Alexander C. Tiegel, Thomas Veness, Piet E. Dargel, Andreas Honecker, Thomas Pruschke, Ian P. McCulloch, Fabian H. L. Essler

Abstract: We consider the optical conductivity $σ_1(ω)$ in the metallic phase of the one-dimensional Hubbard model. Our results focus on the vicinity of half filling and the frequency regime around the optical gap in the Mott insulating phase. By means of a density-matrix renormalization group implementation of the correction-vector approach, $σ_1(ω)$ is computed for a range of interaction strengths and dop… ▽ More We consider the optical conductivity $σ_1(ω)$ in the metallic phase of the one-dimensional Hubbard model. Our results focus on the vicinity of half filling and the frequency regime around the optical gap in the Mott insulating phase. By means of a density-matrix renormalization group implementation of the correction-vector approach, $σ_1(ω)$ is computed for a range of interaction strengths and dopings. We identify an energy scale $E_{\rm opt}$ above which the optical conductivity shows a rapid increase. We then use a mobile impurity model in combination with exact results to determine the behavior of $σ_1(ω)$ for frequencies just above $E_{\rm opt}$ which is in agreement with our numerical data. As a main result, we find that this onset behavior is not described by a power law. △ Less

Submitted 22 March, 2016; v1 submitted 4 January, 2016; originally announced January 2016.

Comments: 6 pages, 5 figures; final version

Journal ref: Phys. Rev. B 93, 125108 (2016)

Imaginary-time matrix product state impurity solver for dynamical mean-field theory

Authors: F. Alexander Wolf, Ara Go, Ian P. McCulloch, Andrew J. Millis, Ulrich Schollwöck

Abstract: We present a new impurity solver for dynamical mean-field theory based on imaginary-time evolution of matrix product states. This converges the self-consistency loop on the imaginary-frequency axis and obtains real-frequency information in a final real-time evolution. Relative to computations on the real-frequency axis, required bath sizes are much smaller and less entanglement is generated, so mu… ▽ More We present a new impurity solver for dynamical mean-field theory based on imaginary-time evolution of matrix product states. This converges the self-consistency loop on the imaginary-frequency axis and obtains real-frequency information in a final real-time evolution. Relative to computations on the real-frequency axis, required bath sizes are much smaller and less entanglement is generated, so much larger systems can be studied. The power of the method is demonstrated by solutions of a three band model in the single and two-site dynamical mean-field approximation. Technical issues are discussed, including details of the method, efficiency as compared to other matrix product state based impurity solvers, bath construction and its relation to real-frequency computations and the analytic continuation problem of quantum Monte Carlo, the choice of basis in dynamical cluster approximation, and perspectives for off-diagonal hybridization functions. △ Less

Submitted 30 July, 2015; originally announced July 2015.

Comments: 8 pages + 4 pages appendix, 9 figures

Journal ref: Phys. Rev. X 5, 041032 (2015)

Chebyshev matrix product state approach for time evolution

Authors: Jad C. Halimeh, Fabian Kolley, Ian P. McCulloch

Abstract: We present and test a new algorithm for time-evolving quantum many-body systems initially proposed by Holzner et al. [Phys. Rev. B 83, 195115 (2011)]. The approach is based on merging the matrix product state (MPS) formalism with the method of expanding the time-evolution operator in Chebyshev polynomials. We calculate time-dependent observables of a system of hardcore bosons quenched under the Bo… ▽ More We present and test a new algorithm for time-evolving quantum many-body systems initially proposed by Holzner et al. [Phys. Rev. B 83, 195115 (2011)]. The approach is based on merging the matrix product state (MPS) formalism with the method of expanding the time-evolution operator in Chebyshev polynomials. We calculate time-dependent observables of a system of hardcore bosons quenched under the Bose-Hubbard Hamiltonian on a one-dimensional lattice. We compare the new algorithm to more standard methods using the MPS architecture. We find that the Chebyshev method gives numerically exact results for small times. However, the reachable times are smaller than the ones obtained with the other state-of-the-art methods. We further extend the new method using a spectral-decomposition-based projective scheme that utilizes an effective bandwidth significantly smaller than the full bandwidth, leading to longer evolution times than the non-projective method and more efficient information storage, data compression, and less computational effort. △ Less

Submitted 5 July, 2015; originally announced July 2015.

Comments: 14 pages, 8 figures

Journal ref: Phys. Rev. B 92, 115130 (2015)

Strongly interacting bosons on a three-leg ladder in the presence of a homogeneous flux

Authors: F. Kolley, M. Piraud, I. P. McCulloch, U. Schollwoeck, F. Heidrich-Meisner

Abstract: We perform a density-matrix renormalization-group study of strongly interacting bosons on a three-leg ladder in the presence of a homogeneous flux. Focusing on one-third filling, we explore the phase diagram in dependence of the magnetic flux and the inter-leg tunneling strength. We find several phases including a Meissner phase, vortex liquids, a vortex lattice, as well as a staggered-current pha… ▽ More We perform a density-matrix renormalization-group study of strongly interacting bosons on a three-leg ladder in the presence of a homogeneous flux. Focusing on one-third filling, we explore the phase diagram in dependence of the magnetic flux and the inter-leg tunneling strength. We find several phases including a Meissner phase, vortex liquids, a vortex lattice, as well as a staggered-current phase. Moreover, there are regions where the chiral current reverses its direction, both in the Meissner and in the staggered-current phase. While the reversal in the latter case can be ascribed to spontaneous breaking of translational invariance, in the first it stems from an effective flux increase in the rung direction. Interactions are a necessary ingredient to realize either type of chiral-current reversal. △ Less

Submitted 15 September, 2015; v1 submitted 24 June, 2015; originally announced June 2015.

Journal ref: New J. Phys. 17 (2015) 092001

Spontaneous increase of magnetic flux and chiral-current reversal in bosonic ladders: Swimming against the tide

Authors: S. Greschner, M. Piraud, F. Heidrich-Meisner, I. P. McCulloch, U. Schollwöck, T. Vekua

Abstract: The interplay between spontaneous symmetry breaking in many-body systems, the wavelike nature of quantum particles and lattice effects produces an extraordinary behavior of the chiral current of bosonic particles in the presence of a uniform magnetic flux defined on a two-leg ladder. While non-interacting as well as strongly interacting particles, stirred by the magnetic field, circulate along the… ▽ More The interplay between spontaneous symmetry breaking in many-body systems, the wavelike nature of quantum particles and lattice effects produces an extraordinary behavior of the chiral current of bosonic particles in the presence of a uniform magnetic flux defined on a two-leg ladder. While non-interacting as well as strongly interacting particles, stirred by the magnetic field, circulate along the system's boundary in the counterclockwise direction in the ground state, interactions stabilize vortex lattices. These states break translational symmetry, which can lead to a reversal of the circulation direction. Our predictions could readily be accessed in quantum gas experiments with existing setups or in arrays of Josephson junctions. △ Less

Submitted 19 December, 2016; v1 submitted 24 April, 2015; originally announced April 2015.

Comments: 5 pages + 5 pages of supplementary material

Journal ref: Phys. Rev. Lett. 115, 190402 (2015)

Non-thermal melting of Neel order in the Hubbard model

Authors: Karsten Balzer, F. Alexander Wolf, Ian P. McCulloch, Philipp Werner, Martin Eckstein

Abstract: We study the unitary time evolution of antiferromagnetic order in the Hubbard model after a quench starting from the perfect Néel state. In this setup, which is well suited for experiments with cold atoms, one can distinguish fundamentally different pathways for melting of long-range order at weak and strong interaction. In the Mott insulating regime, melting of long-range order occurs due to the… ▽ More We study the unitary time evolution of antiferromagnetic order in the Hubbard model after a quench starting from the perfect Néel state. In this setup, which is well suited for experiments with cold atoms, one can distinguish fundamentally different pathways for melting of long-range order at weak and strong interaction. In the Mott insulating regime, melting of long-range order occurs due to the ultra-fast transfer of energy from charge excitations to the spin background, while local magnetic moments and their exchange coupling persist during the process. The latter can be demonstrated by a local spin-precession experiment. At weak interaction, local moments decay along with the long-range order. The dynamics is governed by residual quasiparticles, which are reflected in oscillations of the off-diagonal components of the momentum distribution. Such oscillations provide an alternative route to study the prethermalization phenomenon and its influence on the dynamics away from the integrable (noninteracting) limit. The Hubbard model is solved within nonequilibrium dynamical mean-field theory, using the density matrix-renormalization group as an impurity solver. △ Less

Submitted 9 April, 2015; originally announced April 2015.

Comments: 11 pages

Journal ref: Phys. Rev. X 5, 031039 (2015)

Phase Diagram of the Spin-$1/2$ Triangular $J_1$-$J_2$ Heisenberg Model on a 3-leg Cylinder

Authors: S. N. Saadatmand, B. J. Powell, I. P. McCulloch

Abstract: We study the phase diagram of the frustrated Heisenberg model on the triangular lattice with nearest and next-nearest neighbor spin exchange coupling, on 3-leg ladders. Using the density-matrix renormalization-group method, we obtain the complete phase diagram of the model, which includes quasi-long-range $120^\circ$ and columnar order, and a Majumdar-Ghosh phase with short-ranged correlations. Al… ▽ More We study the phase diagram of the frustrated Heisenberg model on the triangular lattice with nearest and next-nearest neighbor spin exchange coupling, on 3-leg ladders. Using the density-matrix renormalization-group method, we obtain the complete phase diagram of the model, which includes quasi-long-range $120^\circ$ and columnar order, and a Majumdar-Ghosh phase with short-ranged correlations. All these phases are non-chiral and planar. We also identify the nature of phase transitions. △ Less

Submitted 12 April, 2018; v1 submitted 25 March, 2015; originally announced March 2015.

Comments: 12 pages, 15 figures, 1 table, and no supplemental materials. v3: equivalent to the published version

Journal ref: Phys. Rev. B 91, 245119 (2015)

Spectral functions and time evolution from the Chebyshev recursion

Authors: F. Alexander Wolf, Jorge A. Justiniano, Ian P. McCulloch, Ulrich Schollwöck

Abstract: We link linear prediction of Chebyshev and Fourier expansions to analytic continuation. We push the resolution in the Chebyshev-based computation of $T=0$ many-body spectral functions to a much higher precision by deriving a modified Chebyshev series expansion that allows to reduce the expansion order by a factor $\sim\frac{1}{6}$. We show that in a certain limit the Chebyshev technique becomes eq… ▽ More We link linear prediction of Chebyshev and Fourier expansions to analytic continuation. We push the resolution in the Chebyshev-based computation of $T=0$ many-body spectral functions to a much higher precision by deriving a modified Chebyshev series expansion that allows to reduce the expansion order by a factor $\sim\frac{1}{6}$. We show that in a certain limit the Chebyshev technique becomes equivalent to computing spectral functions via time evolution and subsequent Fourier transform. This introduces a novel recursive time evolution algorithm that instead of the group operator $e^{-iHt}$ only involves the action of the generator $H$. For quantum impurity problems, we introduce an adapted discretization scheme for the bath spectral function. We discuss the relevance of these results for matrix product state (MPS) based DMRG-type algorithms, and their use within dynamical mean-field theory (DMFT). We present strong evidence that the Chebyshev recursion extracts less spectral information from $H$ than time evolution algorithms when fixing a given amount of created entanglement. △ Less

Submitted 27 April, 2015; v1 submitted 28 January, 2015; originally announced January 2015.

Comments: 12 pages + 6 pages appendix, 11 figures

Journal ref: Phys. Rev. B 91, 115144 (2015)

A Strictly Single-Site DMRG Algorithm with Subspace Expansion

Authors: Claudius Hubig, Ian P. McCulloch, Ulrich Schollwöck, F. Alexander Wolf

Abstract: We introduce a strictly single-site DMRG algorithm based on the subspace expansion of the Alternating Minimal Energy (AMEn) method. The proposed new MPS basis enrichment method is sufficient to avoid local minima during the optimisation, similarly to the density matrix perturbation method, but computationally cheaper. Each application of $\hat H$ to $|Ψ\rangle$ in the central eigensolver is reduce… ▽ More We introduce a strictly single-site DMRG algorithm based on the subspace expansion of the Alternating Minimal Energy (AMEn) method. The proposed new MPS basis enrichment method is sufficient to avoid local minima during the optimisation, similarly to the density matrix perturbation method, but computationally cheaper. Each application of $\hat H$ to $|Ψ\rangle$ in the central eigensolver is reduced in cost for a speed-up of $\approx (d + 1)/2$, with $d$ the physical site dimension. Further speed-ups result from cheaper auxiliary calculations and an often greatly improved convergence behaviour. Runtime to convergence improves by up to a factor of 2.5 on the Fermi-Hubbard model compared to the previous single-site method and by up to a factor of 3.9 compared to two-site DMRG. The method is compatible with real-space parallelisation and non-abelian symmetries. △ Less

Submitted 12 April, 2015; v1 submitted 22 January, 2015; originally announced January 2015.

Comments: 9 pages, 6 figures; added comparison with two-site DMRG

Journal ref: Phys. Rev. B 91, 155115 (2015)

Fast convergence of imaginary time evolution tensor network algorithms by recycling the environment

Authors: Ho N. Phien, Ian P. McCulloch, Guifré Vidal

Abstract: We propose an environment recycling scheme to speed up a class of tensor network algorithms that produce an approximation to the ground state of a local Hamiltonian by simulating an evolution in imaginary time. Specifically, we consider the time-evolving block decimation (TEBD) algorithm applied to infinite systems in 1D and 2D, where the ground state is encoded, respectively, in a matrix product… ▽ More We propose an environment recycling scheme to speed up a class of tensor network algorithms that produce an approximation to the ground state of a local Hamiltonian by simulating an evolution in imaginary time. Specifically, we consider the time-evolving block decimation (TEBD) algorithm applied to infinite systems in 1D and 2D, where the ground state is encoded, respectively, in a matrix product state (MPS) and in a projected entangled-pair state (PEPS). An important ingredient of the TEBD algorithm (and a main computational bottleneck, especially with PEPS in 2D) is the computation of the so-called environment, which is used to determine how to optimally truncate the bond indices of the tensor network so that their dimension is kept constant. In current algorithms, the environment is computed at each step of the imaginary time evolution, to account for the changes that the time evolution introduces in the many-body state represented by the tensor network. Our key insight is that close to convergence, most of the changes in the environment are due to a change in the choice of gauge in the bond indices of the tensor network, and not in the many-body state. Indeed, a consistent choice of gauge in the bond indices confirms that the environment is essentially the same over many time steps and can thus be re-used, leading to very substantial computational savings. We demonstrate the resulting approach in 1D and 2D by computing the ground state of the quantum Ising model in a transverse magnetic field. △ Less

Submitted 27 March, 2015; v1 submitted 3 November, 2014; originally announced November 2014.

Comments: 17 pages, 28 figures

Journal ref: Phys. Rev. B 91, 115137 (2015)

Phase diagram of the J1-J2 Heisenberg model on the kagome lattice

Authors: Fabian Kolley, Stefan Depenbrock, Ian P. McCulloch, Ulrich Schollwöck, Vincenzo Alba

Abstract: We perform an extensive density matrix renormalization group (DMRG) study of the ground-state phase diagram of the spin-1/2 J_1-J_2 Heisenberg model on the kagome lattice. We focus on the region of the phase diagram around the kagome Heisenberg antiferromagnet, i.e., at J_2=0. We investigate the static spin structure factor, the magnetic correlation lengths, and the spin gaps. Our results are cons… ▽ More We perform an extensive density matrix renormalization group (DMRG) study of the ground-state phase diagram of the spin-1/2 J_1-J_2 Heisenberg model on the kagome lattice. We focus on the region of the phase diagram around the kagome Heisenberg antiferromagnet, i.e., at J_2=0. We investigate the static spin structure factor, the magnetic correlation lengths, and the spin gaps. Our results are consistent with the absence of magnetic order in a narrow region around J_2\approx 0, although strong finite-size effects do not allow us to accurately determine the phase boundaries. This result is in agreement with the presence of an extended spin-liquid region, as it has been proposed recently. Outside the disordered region, we find that for ferromagnetic and antiferromagnetic J_2 the ground state displays signatures of the magnetic order of the \sqrt{3}\times\sqrt{3} and the q=0 type, respectively. Finally, we focus on the structure of the entanglement spectrum (ES) in the q=0 ordered phase. We discuss the importance of the choice of the bipartition on the finite-size structure of the ES. △ Less

Submitted 24 March, 2015; v1 submitted 29 October, 2014; originally announced October 2014.

Comments: 9 pages, 10 figures. Figure added, minor modifications, as published

Journal ref: Phys. Rev. B 91, 104418 (2015)

Solving nonequilibrium dynamical mean-field theory using matrix product states

Authors: F. Alexander Wolf, Ian P. McCulloch, Ulrich Schollwöck

Abstract: We solve the nonequilibrium dynamical mean-field theory (DMFT) using matrix product states (MPS). This allows us to treat much larger bath sizes and by that reach substantially longer times (factor $\sim$ 2 -- 3) than with exact diagonalization. We show that the star geometry of the underlying impurity problem can have substantially better entanglement properties than the previously favoured chain… ▽ More We solve the nonequilibrium dynamical mean-field theory (DMFT) using matrix product states (MPS). This allows us to treat much larger bath sizes and by that reach substantially longer times (factor $\sim$ 2 -- 3) than with exact diagonalization. We show that the star geometry of the underlying impurity problem can have substantially better entanglement properties than the previously favoured chain geometry. This has immense consequences for the efficiency of an MPS-based description of general impurity problems: in the case of equilibrium DMFT, it leads to an orders-of-magnitude speedup. We introduce an approximation for the two-time hybridization function that uses time-translational invariance, which can be observed after a certain relaxation time after a quench to a time-independent Hamiltonian. △ Less

Submitted 19 December, 2014; v1 submitted 13 October, 2014; originally announced October 2014.

Comments: 11 pages + 3 pages appendix, 14 figures

Journal ref: Phys. Rev. B 90, 235131 (2014)