Search found 8 matches

by QS Wang
02 Apr 2021, 21:51
Forum: HowTos and FAQ for TeNPy
Topic: Calculate Energy Expectation on Ladder
Replies: 3
Views: 2453

Re: Calculate Energy Expectation on Ladder

Thanks! Another thing, now I am trying to calculate the spectrum function, for which I need to do tebd on the ladder model. However, it raises error when I try to calculate H_bond. When I try to get H_bond with calc_H_bond(), it says ------------------------------------------------------------------...
by QS Wang
28 Mar 2021, 19:45
Forum: HowTos and FAQ for TeNPy
Topic: Calculate Energy Expectation on Ladder
Replies: 3
Views: 2453

Re: Calculate Energy Expectation on Ladder

The Hamiltonian is

\(H_tot=-t \sum_i (c_i^\dagger c_{i+1} +h.c.) - h \sum_i (\chi_i^\dagger \chi_{i+1}+h.c.) +V\sum_i n_{c,i} n_{c,i+1} + g \sum_i n_{c,i} n_{x,i}\)
by QS Wang
27 Mar 2021, 22:54
Forum: HowTos and FAQ for TeNPy
Topic: Calculate Energy Expectation on Ladder
Replies: 3
Views: 2453

Calculate Energy Expectation on Ladder

Hello. I am working with an impurity-electron model on a ladder. The Hamiltonian is H_tot=-t \sum_i (c_i^\dagger c_{i+1} +h.c.) - h \sum_i (\chi_i^\dagger \chi_{i+1}+h.c.) +V\sum_i n_{c,i} n_{c,i+1} + g \sum_i n_{c,i} n_{x,i} where C is annihilator of a host Fermion and \chi that of an impority. To ...
by QS Wang
23 Mar 2021, 19:04
Forum: HowTos and FAQ for TeNPy
Topic: Particle Number Conservation
Replies: 7
Views: 3999

Re: Particle Number Conservation

Hi, it seems that this problem only appear at the first site (or the first two sites for ladder). That is, for example, if I am working with a ladder with length of 10, then the problem will arise at site 0 and up to now for i>1 the function MPS.apply_local_op(i,op) works well.
by QS Wang
23 Mar 2021, 18:55
Forum: HowTos and FAQ for TeNPy
Topic: Particle Number Conservation
Replies: 7
Views: 3999

Re: Particle Number Conservation

What I want to calculate is \(\langle GS|C_p H C_p^\dagger |GS \rangle \). I planed to apply Cd on each site and sum them up with some phase factors to construct \(C_p^\dagger GS \rangle \) and calculate the expectation of \(H\) for this state.
by QS Wang
23 Mar 2021, 15:22
Forum: HowTos and FAQ for TeNPy
Topic: Particle Number Conservation
Replies: 7
Views: 3999

Re: Particle Number Conservation

After attaining the ground state I tried to apply C^\dagger_p on it to calculate the energy spectrum. To do so I need to apply Cd to every sites and times the result with some phase factor. However, I was met with a error saying "s has wrong shape: (1,) instead of 0" (shown below). The cod...
by QS Wang
22 Mar 2021, 20:58
Forum: HowTos and FAQ for TeNPy
Topic: Particle Number Conservation
Replies: 7
Views: 3999

Re: Particle Number Conservation

All of plus_hc are set to be True since tenpy always raises TypeError once it is False. TypeError Traceback (most recent call last) <ipython-input-5-73d291e71b76> in <module>() 18 } 19 ---> 20 TestModel=GroupedBosonImpModel(test_params) 3 frames /usr/local/lib/python3.7/dist-packages/tenpy/models/mo...
by QS Wang
22 Mar 2021, 10:47
Forum: HowTos and FAQ for TeNPy
Topic: Particle Number Conservation
Replies: 7
Views: 3999

Particle Number Conservation

Hello, I am trying to find the ground state and calculate the ground state energy of a spinless fermi system using GroupedSite and DMRG. However, the problem is that I need the number of impurity to be fixed, while actually it always change and finally DMRG always give me the global minimum with ano...