Search found 32 matches
- 02 Apr 2020, 05:39
- Forum: HowTos and FAQ for TeNPy
- Topic: Question about generating random MPS
- Replies: 1
- Views: 7271
Question about generating random MPS
Hi everyone, I'm trying to generate a random MPS by using the example in test module of tenpy (see the end of this post). I found that when L > 140, it will raise RuntimeError: SVD found no singluar values. I'm not sure where I make wrong, does anyone know something about this? thanks, Qicheng def r...
- 12 Mar 2020, 08:07
- Forum: HowTos and FAQ for TeNPy
- Topic: leading eigenvector of transfer matrix
- Replies: 6
- Views: 12863
Re: leading eigenvector of transfer matrix
You cannot/should not take the transfer matrix of a single site, if you have a two-site unit cell in your MPS. Sometimes, the charge values will not even allow you, and even if they do, it there is a different gauge on the two bonds - which can and will be the case if there are degenerate singular ...
- 11 Mar 2020, 04:32
- Forum: HowTos and FAQ for TeNPy
- Topic: leading eigenvector of transfer matrix
- Replies: 6
- Views: 12863
Re: leading eigenvector of transfer matrix
BTW, in this sense, the normalization factors of the left (to square singular value) and right eigenvector (identity) of right canonical form should cancel each other out, right? However, in a real case it seems will change the largest amplitude of the transfer matrix (the largest eigenvalue will b...
- 11 Mar 2020, 03:48
- Forum: HowTos and FAQ for TeNPy
- Topic: leading eigenvector of transfer matrix
- Replies: 6
- Views: 12863
Re: leading eigenvector of transfer matrix
It's just a matter of normalization. Eigenvalues get normalized to np.linalg.norm(v) = 1 , but in this case you want np.sum(v) = 1 . At the end of your code, add vdiag = v.diagonal() print(vdiag/np.sum(vdiag)) print(psi.get_SL(0)**2) and compare the results, the prints should give the same Thanks f...
- 10 Mar 2020, 13:48
- Forum: HowTos and FAQ for TeNPy
- Topic: leading eigenvector of transfer matrix
- Replies: 6
- Views: 12863
leading eigenvector of transfer matrix
Hi everyone, it is well known that, in left(right) canonical form infinite MPS, the leading left(right) eigenvector (corresponding to eigenvalue 1) of the transfer matrix is a unit matrix \mathbb{1} , and the leading right(left) eigenvector is \lambda^2 \mathbb{1} . However, when I trying to test th...
- 10 Mar 2020, 06:47
- Forum: HowTos and FAQ for TeNPy
- Topic: Real-time evolution of the GS of the U=0 Fermi-Hubbard model
- Replies: 5
- Views: 9166
Re: Real-time evolution of the GS of the U=0 Fermi-Hubbard model
Yes, that's for sure a good idea. Indeed, I thought about it and started to implement it already a few days ago, but I'm quite busy these days (moving to UC Berkeley :)), so I didn't have time to finish it. Hopefully, I can find a bit time next week - unless you volunteer to time to implement it ;-...
- 07 Dec 2019, 09:25
- Forum: HowTos and FAQ for TeNPy
- Topic: Real-time evolution of the GS of the U=0 Fermi-Hubbard model
- Replies: 5
- Views: 9166
Re: Real-time evolution of the GS of the U=0 Fermi-Hubbard model
Sorry for the late reply, I just see your message.
- 14 Nov 2019, 03:02
- Forum: HowTos and FAQ for TeNPy
- Topic: Real-time evolution of the GS of the U=0 Fermi-Hubbard model
- Replies: 5
- Views: 9166
Re: Real-time evolution of the GS of the U=0 Fermi-Hubbard model
it is a good idea to write a function to evaluate the variance of H, and maybe something like |<\psi|H|\psi>|^2 - <\psi|H^2|\psi> can be used in dmrg to be a control parameter? What order (second or fourth) and timestep do you use for TEBD? Does this go away if you simply use a smaller timestep/high...
- 01 Oct 2019, 03:14
- Forum: HowTos and FAQ for TeNPy
- Topic: print msg into file
- Replies: 6
- Views: 10280
Re: print msg into file
Thanks for your help Johannes, I actually do not know how to deal with this problem for a long time (not just for tenpy, also other code written by myself or others).
- 26 Sep 2019, 12:05
- Forum: HowTos and FAQ for TeNPy
- Topic: print msg into file
- Replies: 6
- Views: 10280
Re: print msg into file
maybe i didnot explain my point clearly enough. The point is to save the printing msg in terminal at time . If simply use python3 abc.py > log the printing msg will be saved after the task is done, i.e. during running of the task, you can not read the msg in this log file. In the case of using a pyt...
- 18 Sep 2019, 03:32
- Forum: HowTos and FAQ for TeNPy
- Topic: correlation length in idmrg
- Replies: 1
- Views: 3731
correlation length in idmrg
Hi, everyone. I'm trying to calculate the correlation length for TFI chain near or away from the critical point. I have test before, the calculated correlation length for critical point g=1 is agree with the pervious numeric presented by Frank Pollmann https://journals.aps.org/prl/abstract/10.1103/P...
- 18 Sep 2019, 03:22
- Forum: HowTos and FAQ for TeNPy
- Topic: how to write a three-site interaction?
- Replies: 4
- Views: 33811
Re: how to write a three-site interaction?
Thanks for your explanation, Umberto. It's very clear.
- 05 Sep 2019, 00:46
- Forum: HowTos and FAQ for TeNPy
- Topic: how to write a three-site interaction?
- Replies: 4
- Views: 33811
how to write a three-site interaction?
Hi everyone, I'm confused how to write a three-site interaction? I have read the example of toric code, but I donnot understand. If I want to write a O’Brien-Fendley model H=\sum_{j=1}^N [-\sigma^z_{j}\sigma^z_{j+1}-\sigma^x_{j}+\lambda(\sigma^x_{j}\sigma^z_{j+1}\sigma^z_{j+2}+\sigma^z_{j}\sigma^z_{...
- 04 Sep 2019, 07:14
- Forum: HowTos and FAQ for TeNPy
- Topic: print msg into file
- Replies: 6
- Views: 10280
print msg into file
I'm wondering if the msg can be printed into a log file?
It's not hard to modify the code, and I think this is necessary.
It's not hard to modify the code, and I think this is necessary.
- 02 Jul 2019, 13:51
- Forum: HowTos and FAQ for TeNPy
- Topic: Some questions about the svd_theta
- Replies: 7
- Views: 9990
Re: Some questions about the svd_theta
You can only split legs that were combined previously. In TeNPy, you can see that those legs store the charge data as tenpy.linalg.charges.LegPipe instead of tenpy.linalg.charges.LegCharge . Those pipes have the additional data needed to split the legs. The matrices generated by SVD, say A_ij = U_i...
- 02 Jul 2019, 13:29
- Forum: HowTos and FAQ for TeNPy
- Topic: Some questions about the svd_theta
- Replies: 7
- Views: 9990
Re: Some questions about the svd_theta
Thank you, Johannes, now I'm clear with the svd process, but here appears another question. As in the note said, the process of combination of two legs is something like mapping the indices for "vL" and "p" into something like "2*vL+p", so you can split out the physical...
- 01 Jul 2019, 04:16
- Forum: HowTos and FAQ for TeNPy
- Topic: Some questions about the svd_theta
- Replies: 7
- Views: 9990
Re: Some questions about the svd_theta
There's another thing related, and I donnot understand how to deal with it. For example, consider a spin-half chain with conserved Sz, we can write down the 3 sub-blocks with Sz "-2, 0, 2", In this case, we can do svd for each charge sector. But in the case of, for example, TFI chain where...
- 28 Jun 2019, 17:27
- Forum: HowTos and FAQ for TeNPy
- Topic: Some questions about the svd_theta
- Replies: 7
- Views: 9990
Re: Some questions about the svd_theta
Johannes, Thanks for your reply, I know that the transfer matrix method is only well-defined in a case with translational invariance, and for the contracting of the transfer matrix we need right matrice B or left matrice A. My question is about the implementation of the charge conservation in the ca...
- 28 Jun 2019, 06:18
- Forum: HowTos and FAQ for TeNPy
- Topic: Some questions about the svd_theta
- Replies: 7
- Views: 9990
Some questions about the svd_theta
Hi everyone, I want to write a little program to calculate the correlation length with a given "theta" (which is located in the center of the chain), so I need to do svd to theta and get the right matrice B. But there is a problem, I get the theta by using function "get_theta" fo...
- 20 Apr 2019, 17:13
- Forum: Algorithms
- Topic: question about tdvp
- Replies: 2
- Views: 15958
- 08 Apr 2019, 09:18
- Forum: Algorithms
- Topic: question about tdvp
- Replies: 2
- Views: 15958
question about tdvp
Hi, everyone. I'm trying to do tdvp, but there are something in the original paper which I cannot understand. Firstly, why the projector in this form? I know the projector is something like 1 - \sum_N |\psi_N><\psi_N| , but what's the target space here, it seems like we project it onto a basis orort...
- 08 Mar 2019, 10:50
- Forum: HowTos and FAQ for TeNPy
- Topic: Product state results
- Replies: 3
- Views: 5548
Re: Product state results
I believe it's easy to generate a random product state with a given total charge.
for example, for qubit, it's just a combination of 0 and 1 with a given sum of all elements in the list, you just need to give the number of 0 and 1, and do a random permutation.
for example, for qubit, it's just a combination of 0 and 1 with a given sum of all elements in the list, you just need to give the number of 0 and 1, and do a random permutation.
- 27 Feb 2019, 10:07
- Forum: HowTos and FAQ for TeNPy
- Topic: Product state results
- Replies: 3
- Views: 5548
Re: Product state results
I will ask you the question what Johannes asked me: What the "random" state means here? It should be noticed that the charge (quantum number) should be conserved in dmrg, so if your random state has a different total charge, the result will be wrong. Perhaps you can use "RandomUnitary...
- 22 Feb 2019, 09:34
- Forum: HowTos and FAQ for TeNPy
- Topic: Scaling of entanglement entropy for iDMRG
- Replies: 7
- Views: 9485
Re: Scaling of entanglement entropy for iDMRG
hi, everyone. I have to admit i was wrong, and what Umberto said is correct: The original formula is correct, and in infinite dmrg the boundary condition is open. BTW, bart, your scaling method may be wrong. Why you choose the data from a state which is not converged in dmrg? What I remember is: you...
- 09 Feb 2019, 13:30
- Forum: HowTos and FAQ for TeNPy
- Topic: How to generate a random MPS in tenpy
- Replies: 5
- Views: 9174
Re: How to generate a random MPS in tenpy
Thanks for your reply! In fact, I want to get some states with random entanglement entropy. I know the random unitary time evolution will be a good method. I was thinking about that is there a method can directly obtain these states? The first coming problem is about the charge conserve as you said,...