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 ororthogonal to the original wavefunction. I can understand this minus part in twosite tdvp, it's just minus a onsite interaction (since we counter twice onsite interaction in twosite update), is this true?
Secondly, what's the zerosite Hamiltonian? I didn't do the onesite dmrg (I always do the twosite dmrg) before, so I have no knowledge about this. In tenpy, it's just the singular value, but I do not know why we need to minus the singular value. Is this the same in the twosite tdvp? If so, there's a question: in onesite update, we actually only update one singular value in one time. From this reason, I don't know why we need to minus the singular value, or it's different from the twosite tdvp?
Thirdly, if the twosite tdvp is what I said "just minus a onsite interaction", why not simply do lanczos for \(H_{bond}^i  H_{onesite}^{i+1} \otimes Id^i\). In fact, this's what we did before in tebd, it will be much faster than do lanczos twice (one for \(H_{bond}\), another for \(H_{onesite}\)), or I'm wrong with something?
Thanks, Qicheng
New release: v0.7.1.
question about tdvp

 Posts: 1
 Joined: 11 Apr 2019, 13:57
Re: question about tdvp
Hi Qicheng,
In order to answer your questions, I think it is best to go back to the mathematical construction behind TDVP. In my opinion, this is the clearest in the paper entitled "Unifying time evolution and optimization with matrix product states" (arXiv:1408.5056) from Jutho Haegeman et al.
The general idea is to insert a MPS with a given bond dimension \(\chi\) as an ansatz in the Schroedinger equation and to solve it for this ansatz. In order to illustrate, let's do it for three site:
\(\partial_t \sum_{s_1,s_2,s_3} A^{s_1}A^{s_2}A^{s_3} s_1,s_2,s_3>=i H\sum_{s_1,s_2,s_3} A^{s_1}A^{s_2}A^{s_3} s_1,s_2,s_3>\) with \(\hbar=1\).
This gives us:
\( \sum_{s_1,s_2,s_3} \dot{A}^{s_1}A^{s_2}A^{s_3} s_1,s_2,s_3>+A^{s_1}\dot{A}^{s_2}A^{s_3} s_1,s_2,s_3>+A^{s_1}A^{s_2}\dot{A}^{s_3} s_1,s_2,s_3>=i H\sum_{s_1,s_2,s_3} A^{s_1}A^{s_2}A^{s_3} s_1,s_2,s_3> \space (eq\space1)\)
where \(\dot{A}^{s_i}\) is the derivative with respect to time of \(A^{s_i}\).
Two remarks about equation 1:
\( \sum_{s_1,s_2,s_3} \dot{A}^{s_1}A^{s_2}A^{s_3} s_1,s_2,s_3>+A^{s_1}\dot{A}^{s_2}A^{s_3} s_1,s_2,s_3>+A^{s_1}A^{s_2}\dot{A}^{s_3} s_1,s_2,s_3>=i P_T H\sum_{s_1,s_2,s_3} A^{s_1}A^{s_2}A^{s_3} s_1,s_2,s_3> \space (eq\space1*)\)
where \(P_T\) is called the projector on the tangent space.
\( \sum_{s_1,s_2,s_3} B^{s_1}A^{s_2}A^{s_3} s_1,s_2,s_3>+A^{s_1}B^{s_2}A^{s_3} s_1,s_2,s_3>+A^{s_1}A^{s_2}B^{s_3} s_1,s_2,s_3>\) where \(B^{s_i}\) can be any tensor of compatible shape. Note that the left hand side of equations 1 and 1* belong to this space. Adding the projection allow the time evolution to stay confined in the tangent space.
The derivation of the precise form of the projector is too long to be explained here, but I refer to the supplementary material of the paper mentioned above for the details. The general idea is to find the vector of the tangent space which is the closest to \(iH\sum_{s_1,s_2,s_3} A^{s_1}A^{s_2}A^{s_3} s_1,s_2,s_3>\) (using the usual norm as a measure of the distance between two vectors).
\(
\dot{x}=f(x,y)
\)
\(
\dot{y}=g(x,y)
\) ?
Answer: First integrate the first line with small dt for y constant \(\dot{x}=f(x(t),y(t))\). Then integrate the second one using but using the update \(x(t+dt)\): \(\dot{y}=g(x(t+dt),y(t))\). This is called a splitting method, which is a good approximation if \(dt\) is small. We do the same in TDVP, which answer your last question:
Lastly:
I hope this is helpful, but if you really want to understand everything 100%, the best is to read the paper.
You can also look at the following lectures notes arXiv:1810.07006 which are interesting.
In order to answer your questions, I think it is best to go back to the mathematical construction behind TDVP. In my opinion, this is the clearest in the paper entitled "Unifying time evolution and optimization with matrix product states" (arXiv:1408.5056) from Jutho Haegeman et al.
The general idea is to insert a MPS with a given bond dimension \(\chi\) as an ansatz in the Schroedinger equation and to solve it for this ansatz. In order to illustrate, let's do it for three site:
\(\partial_t \sum_{s_1,s_2,s_3} A^{s_1}A^{s_2}A^{s_3} s_1,s_2,s_3>=i H\sum_{s_1,s_2,s_3} A^{s_1}A^{s_2}A^{s_3} s_1,s_2,s_3>\) with \(\hbar=1\).
This gives us:
\( \sum_{s_1,s_2,s_3} \dot{A}^{s_1}A^{s_2}A^{s_3} s_1,s_2,s_3>+A^{s_1}\dot{A}^{s_2}A^{s_3} s_1,s_2,s_3>+A^{s_1}A^{s_2}\dot{A}^{s_3} s_1,s_2,s_3>=i H\sum_{s_1,s_2,s_3} A^{s_1}A^{s_2}A^{s_3} s_1,s_2,s_3> \space (eq\space1)\)
where \(\dot{A}^{s_i}\) is the derivative with respect to time of \(A^{s_i}\).
Two remarks about equation 1:
 1. We do not know the derivative of the tensor a priori. Therefore the algorithm has to find them before integrating them in order to obtain \(A^{s_i}(t)\)
 2. Equation 1 is not really what we want since the left hand side does not leave on the same space than the right hand side. Indeed, the Hamiltonian applied to the MPS will increase the bond dimension.
\( \sum_{s_1,s_2,s_3} \dot{A}^{s_1}A^{s_2}A^{s_3} s_1,s_2,s_3>+A^{s_1}\dot{A}^{s_2}A^{s_3} s_1,s_2,s_3>+A^{s_1}A^{s_2}\dot{A}^{s_3} s_1,s_2,s_3>=i P_T H\sum_{s_1,s_2,s_3} A^{s_1}A^{s_2}A^{s_3} s_1,s_2,s_3> \space (eq\space1*)\)
where \(P_T\) is called the projector on the tangent space.
\(\rightarrow\) We project on the tangent space i.e. the space of the vector of the forms:Firstly, why the projector in this form? [...]it seems like we project it onto a basis ororthogonal to the original wavefunction.
\( \sum_{s_1,s_2,s_3} B^{s_1}A^{s_2}A^{s_3} s_1,s_2,s_3>+A^{s_1}B^{s_2}A^{s_3} s_1,s_2,s_3>+A^{s_1}A^{s_2}B^{s_3} s_1,s_2,s_3>\) where \(B^{s_i}\) can be any tensor of compatible shape. Note that the left hand side of equations 1 and 1* belong to this space. Adding the projection allow the time evolution to stay confined in the tangent space.
The derivation of the precise form of the projector is too long to be explained here, but I refer to the supplementary material of the paper mentioned above for the details. The general idea is to find the vector of the tangent space which is the closest to \(iH\sum_{s_1,s_2,s_3} A^{s_1}A^{s_2}A^{s_3} s_1,s_2,s_3>\) (using the usual norm as a measure of the distance between two vectors).
\(\rightarrow\) Well... things are a bit compicated here. Once we find the projector on the tangent space, we want to integrate the equation of motion (see remark 1 above). It turns out that the projetor can be written elegantly as the sum over all sites of two terms (this is just a mathematical fact, nothing deep about it, the only way to see it is to go through the proof). The problem is that you get 2 N coupled differential equations. The variables are:I can understand this minus part in twosite tdvp, it's just minus a onsite interaction (since we counter twice onsite interaction in twosite update), is this true?
 \(A_C(n,t)\): the tensor on site n in left canonical form
 \(C(n,t)\): the singular value at bound n
\(
\dot{x}=f(x,y)
\)
\(
\dot{y}=g(x,y)
\) ?
Answer: First integrate the first line with small dt for y constant \(\dot{x}=f(x(t),y(t))\). Then integrate the second one using but using the update \(x(t+dt)\): \(\dot{y}=g(x(t+dt),y(t))\). This is called a splitting method, which is a good approximation if \(dt\) is small. We do the same in TDVP, which answer your last question:
We need first to update the tensor, then the singular value, as we did with x and y in the example above.Thirdly, if the twosite tdvp is what I said "just minus a onsite interaction", why not simply do lanczos for \(H_{ibond}−H_{i+1}\)one−site⊗Idi
Lastly:
The principle behind the two site TDVP is the same, except you take the derivative of two tensors at a time in equation (1*). Instead of updating only a singular value, you will update a singular value plus the tensor on the left side. The minus sign just comes from the form of the projector.From this reason, I don't know why we need to minus the singular value, or it's different from the twosite tdvp?
I hope this is helpful, but if you really want to understand everything 100%, the best is to read the paper.
You can also look at the following lectures notes arXiv:1810.07006 which are interesting.

 Posts: 32
 Joined: 08 Jan 2019, 03:03