MPToolkit使用说明(1)

一.建模

1.首先在目录下创建文件名

 cd ~/Software/mptoolkit/models

2.根据结构建模,即保存的lattice

主要注意原胞个数，有几种相互关系，和最后求和

求和时 sum_unit() 有括号

二.准备文件

1.进行编译

在 ~/Software/build/ 目录下编译

需要一个 compile.sh 文件，从models里面拷过来，只需要修改路径和生成的二进制文件的名称即可

编译指令如下：

 cp compile.sh ~/Software/build/
 cd ~/Software/build/
 ./compile.sh

2.给文件权限

若文件不是可编译的绿色，则需要给文件一定的权限

 chmod +x compile.sh

3.将二进制文件存入~/bin 目录

 cd   ~/Software/build
 cp -r spin_Heisenberg_Wen ~/bin/

三.得到波函数，计算基态能量

1.创建 lattice 文件

 spin_Heisenberg_Wen -o lattice

( spin_Heisenberg_Wen 是自己创建的文件名）

2.获得基态波函数

mp-idmrg-s3e -w psi -H lattice:H_J1_J2 -m 10..50x20,50x100 --create -u 2

【（1）lattice：H_J1_J2 每次需要根据lattice 文件中的名字修改

  （2）从10个扫描20次到50个扫描20次，到50个扫描100次

  （3）第一次算得时候需要--creat，之后都不需要 

  （4）如果是对称结构SU2，则

    mp-idmrg-s3e -w psi -H lattice:H_J1_J2 -m 10..50x20,50x100 --create -u 2 -q 0】

更多次迭代，进一步降低能量

  mp-idmrg-s3e -w psi -H lattice:H_J1_J2 -m 50x200

  mp-idmrg-s3e -w psi -H lattice:H_J1_J2 -m 50x200 --mix-factor 1e-5

  mp-idmrg-s3e -w psi -H lattice:H_J1_J2 -m 50x200 --mix-factor 1e-7

  mp-idmrg-s3e -w psi -H lattice:H_J1_J2 -m 50x200 --mix-factor 0

3.能量的偏差

  mp-imoments psi
 
  mp-imoments psi --power 2

下图是最终能量偏差图

【每次迭代完都需要看能量的偏差，观察是否有更精确】