Total particle number can not be conserved using ALPS TEBD simulation
Dear ALPS users, When I follow the ALPS TEBD tutorial about hard-core boson model(http://alps.comp-phys.org/mediawiki/index.php/ALPS_2_Tutorials:TEBD-01_bhque...) by Python, I measured the average particle number for each site, i.e. 'Avg Number' except the 'Loschmidt Echo' and 'V' by the following python code: data = pyalps.load.loadTimeEvolution(pyalps.getResultFiles(prefix='1D_HardcoreBosonModel'), measurements=['Loschmidt Echo', 'V', 'Avg Number']) Then I found the 'Avg Number' equals 0.3! This means that the total particle number is 0.3*10=3, not 5. In the tutorial, the conserved quantum number is 'N_total': 5. When I change the 'N_total' to 1 and 10, I got the same 'Avg Number'=0.3 Can any one give me some help? Thanks a lot! Best regards, Rongyang Sun ----------------------------------------------------------------- Rongyang Sun Ph.D. Candidate Theoretical Quantum Matter Physics Institute for Advanced Study, Tsinghua University, Beijing 100084, China E-mail: sunry15@mails.tsinghua.edu.cn; sun-rongyang@outlook.com
I don’t really know the tebd code, but the mps_evolve has the same features and it is definitely able to conserve quantum numbers as well as running other models. https://alps.comp-phys.org/mediawiki/index.php/ALPS_2_Tutorials:_MPS-05_bhqu... We never really rewrote the text of the tutorial, but in the installation directory you will find a script that runs the bhquench example, as well as others. Best, Michele
On 27 Aug 2017, at 16:52, Sun Rongyang <sun-rongyang@outlook.com> wrote:
Dear ALPS users,
When I follow the ALPS TEBD tutorial about hard-core boson model(http://alps.comp-phys.org/mediawiki/index.php/ALPS_2_Tutorials:TEBD-01_bhque... <http://alps.comp-phys.org/mediawiki/index.php/ALPS_2_Tutorials:TEBD-01_bhquench>) by Python, I measured the average particle number for each site, i.e. 'Avg Number' except the 'Loschmidt Echo' and 'V' by the following python code:
data = pyalps.load.loadTimeEvolution(pyalps.getResultFiles(prefix='1D_HardcoreBosonModel'), measurements=['Loschmidt Echo', 'V', 'Avg Number'])
Then I found the 'Avg Number' equals 0.3! This means that the total particle number is 0.3*10=3, not 5. In the tutorial, the conserved quantum number is 'N_total': 5. When I change the 'N_total' to 1 and 10, I got the same 'Avg Number'=0.3
Can any one give me some help? Thanks a lot!
Best regards, Rongyang Sun
----------------------------------------------------------------- Rongyang Sun Ph.D. Candidate Theoretical Quantum Matter Physics Institute for Advanced Study, Tsinghua University, Beijing 100084, China E-mail: sunry15@mails.tsinghua.edu.cn <mailto:sunry15@mails.tsinghua.edu.cn>; sun-rongyang@outlook.com <mailto:sun-rongyang@outlook.com>
---- Comp-phys-alps-users Mailing List for the ALPS Project http://alps.comp-phys.org/ <http://alps.comp-phys.org/>
List info: https://lists.phys.ethz.ch//listinfo/comp-phys-alps-users <https://lists.phys.ethz.ch//listinfo/comp-phys-alps-users> Archive: https://lists.phys.ethz.ch//pipermail/comp-phys-alps-users <https://lists.phys.ethz.ch//pipermail/comp-phys-alps-users>
Unsubscribe by writing a mail to comp-phys-alps-users-leave@lists.phys.ethz.ch <mailto:comp-phys-alps-users-leave@lists.phys.ethz.ch>.
Thanks so much for your reply! I will try it. By the way, there are may be some missing INTEGER(HID_T) kind definitions of int data type at the Fortran source code about the tebd application, For example, at tebd/core/Hdf5Interface.f90 ------------------------------------------- INTEGER, INTENT(INOUT) :: groupId --> INTEGER(HID_T), INTENT(INOUT) :: groupId ------------------------------------------- When I compile the source code, I met some int kind mismatch errors. These errors disappeared after I change the source code. It may be a bug about ALPS. Best, Rongyang ----------------------------------------------------------------- Rongyang Sun Ph.D. Candidate Theoretical Quantum Matter Physics Institute for Advanced Study, Tsinghua University, Beijing 100084, China E-mail: sunry15@mails.tsinghua.edu.cn; sun-rongyang@outlook.com ________________________________ 发件人: Comp-phys-alps-users <comp-phys-alps-users-bounces@lists.phys.ethz.ch> 代表 Michele Dolfi <michele.dolfi@gmail.com> 发送时间: 2017年9月1日 16:42:03 收件人: comp-phys-alps-users@lists.phys.ethz.ch 主题: Re: [ALPS-users] Total particle number can not be conserved using ALPS TEBD simulation I don’t really know the tebd code, but the mps_evolve has the same features and it is definitely able to conserve quantum numbers as well as running other models. https://alps.comp-phys.org/mediawiki/index.php/ALPS_2_Tutorials:_MPS-05_bhqu... We never really rewrote the text of the tutorial, but in the installation directory you will find a script that runs the bhquench example, as well as others. Best, Michele On 27 Aug 2017, at 16:52, Sun Rongyang <sun-rongyang@outlook.com<mailto:sun-rongyang@outlook.com>> wrote: Dear ALPS users, When I follow the ALPS TEBD tutorial about hard-core boson model(http://alps.comp-phys.org/mediawiki/index.php/ALPS_2_Tutorials:TEBD-01_bhque...) by Python, I measured the average particle number for each site, i.e. 'Avg Number' except the 'Loschmidt Echo' and 'V' by the following python code: data = pyalps.load.loadTimeEvolution(pyalps.getResultFiles(prefix='1D_HardcoreBosonModel'), measurements=['Loschmidt Echo', 'V', 'Avg Number']) Then I found the 'Avg Number' equals 0.3! This means that the total particle number is 0.3*10=3, not 5. In the tutorial, the conserved quantum number is 'N_total': 5. When I change the 'N_total' to 1 and 10, I got the same 'Avg Number'=0.3 Can any one give me some help? Thanks a lot! Best regards, Rongyang Sun ----------------------------------------------------------------- Rongyang Sun Ph.D. Candidate Theoretical Quantum Matter Physics Institute for Advanced Study, Tsinghua University, Beijing 100084, China E-mail: sunry15@mails.tsinghua.edu.cn<mailto:sunry15@mails.tsinghua.edu.cn>; sun-rongyang@outlook.com<mailto:sun-rongyang@outlook.com> ---- Comp-phys-alps-users Mailing List for the ALPS Project http://alps.comp-phys.org/ List info: https://lists.phys.ethz.ch//listinfo/comp-phys-alps-users Archive: https://lists.phys.ethz.ch//pipermail/comp-phys-alps-users Unsubscribe by writing a mail to comp-phys-alps-users-leave@lists.phys.ethz.ch<mailto:comp-phys-alps-users-leave@lists.phys.ethz.ch>.
participants (2)
-
Michele Dolfi
-
Sun Rongyang