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.commailto: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.cnmailto:sunry15@mails.tsinghua.edu.cn; sun-rongyang@outlook.commailto: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.chmailto:comp-phys-alps-users-leave@lists.phys.ethz.ch.
comp-phys-alps-users@lists.phys.ethz.ch