SciPost Submission Page
Evaluation of time-dependent correlators after a local quench in iPEPS: hole motion in the t-J model
by C. Hubig, A. Bohrdt, M. Knap, F. Grusdt, J. I. Cirac
This is not the current version.
|As Contributors:||Claudius Hubig|
|Submitted by:||Hubig, Claudius|
|Submitted to:||SciPost Physics|
|Subject area:||Condensed Matter Physics - Computational|
Infinite projected entangled pair states (iPEPS) provide a convenient variational description of infinite, translationally-invariant two-dimensional quantum states. However, the simulation of local excitations is not directly possible due to the translationally-invariant ansatz. Furthermore, as iPEPS are either identical or orthogonal, expectation values between different states as required during the evaluation of non-equal-time correlators are ill-defined. Here, we show that by introducing auxiliary states on each site, it becomes possible to simulate both local excitations and evaluate non-equal-time correlators in an iPEPS setting under real-time evolution. We showcase the method by simulating the t-J model after a single hole has been placed in the half-filled antiferromagnetic background and evaluating both return probabilities and spin correlation functions, as accessible in quantum gas microscopes.
Ontology / TopicsSee full Ontology or Topics database.
Submission & Refereeing History
Reports on this Submission
Anonymous Report 3 on 2019-12-3 Invited Report
- Cite as: Anonymous, Report on arXiv:1911.01159v1, delivered 2019-12-03, doi: 10.21468/SciPost.Report.1360
In this work the authors present a scheme to compute time-dependent correlation functions based on infinite projected-entangled pair states (iPEPS). Rather than just applying a local operator to the ansatz and time-evolving the state (which would lead to a loss of translational invariance) the scheme keeps the translational invariance of the ansatz by introducing auxiliary states on each site and an operator acting on both the physical and auxiliary space. The accuracy of the operator representation is depends on a small parameter \epsilon (controlling the density of excitations) which is chosen sufficiently small. They test their approach on the t-J model by computing the return probability and diagonal spin correlators of a single hole doped in the antiferromagnetic background at half filling, and they present comparisons with results obtained with matrix product states on cylinders.
Studying local quenches in 2D is a very challenging and important problem. This work introduces a very interesting and useful trick in this context which I believe will be very useful in future iPEPS simulations. As the authors mention, there is still room for improvement in the time evolution scheme used (i.e. a stable full update evolution rather than the less accurate simple update scheme), but even the present results clearly demonstrate the usefulness and applicability of their approach.
For these reasons I can recommend publication of this article in SciPost. The authors may want to take into account the comments and questions listed below when revising their paper.
Comments and questions
(1) In the results section the authors use an excitation density \epsilon of 10^-2, and they mention that ideally one would want to consider the case \epsilon -> 0. While I can imagine that most of the error comes from the truncation due to the finite D, it would nevertheless be good to discuss the dependence of the results on epsilon (e.g. showing some example data as a function of epsilon).
(2) It would be interesting to discuss also the dependence on the trotter step used.
(3) Is there any criterion which can be used in order to judge when the approach starts breaking down? For example, can one observe that at after a certain time scale the singular value spectrum (on the bonds in the simple update scheme) becomes very flat, suggesting that the entanglement is too large / the bond dimension D too small?
Report 2 by Didier Poilblanc on 2019-11-20 Invited Report
- Cite as: Didier Poilblanc, Report on arXiv:1911.01159v1, delivered 2019-11-20, doi: 10.21468/SciPost.Report.1332
Cutting-edge developments by experts in the field of tensor networks
Application of the method to an interesting well-studied problem
Part of the old literature on the one hole t-J model is ignored and may be relevant for interpreting the iPEPS results.
The authors introduce a new method to study, in a two dimensional correlated system, the time evolution after a sudden quench by adapting iPEPS techniques. These are clearly cutting-edge developments by experts in the field of tensor networks.
The authors benchmark the method by investigated the dynamics of a single hole in a quantum antiferromagnet - the so-called t-J model - a well-studied problem in the late 80's - early 90's after the discovery of high-Tc superconductors. Reliable techniques like Lanczos exact diagonalizations on small torii were used to obtain accurately the hole dispersion relation and to establish the composite nature of the hole excitation and the corresponding string-like force. These points were developed in early key references like:
E. Dagotto, R. Joynt, A. Moreo, S. Bacci and E. Gagliano, Phys. Rev. B 41, 9049 (1990).
Didier Poilblanc, H. J. Schulz, and Timothy Ziman, Phys. Rev. B 46, 6435 (1992)
Didier Poilblanc, H. J. Schulz, and Timothy Ziman, Phys. Rev. B 47, 3268 (1993)
D. Poilblanc, T. Ziman, H.J. Schulz and E. Dagotto, Phys. Rev. B 47, 14267 (1993).
P.Béran, D.Poilblanc and R.B.Laughlin, Nuclear Physics B 473, 707-720 (1996).
The string picture developed in the above references seems to be an essential ingredient to interpret the return probability computed by the authors.
The authors may consider (and cite) the above reference. They may also comment on whether any evidence for the holon-spinon string-like force is seen in their calculation.