mp-itebd + mp-iexpectation-cross

以2021年PRL中的xxz模型为例

其哈密顿量为

初态为

末态为

代码为

    cd ~/Wen/xxz/itebd/

    spinchain -o lattice

    mp-idmrg-s3e -w psi.t0.0 -H "lattice:-H_xx-*H_x"  -m 10x100 --create -u 2
    mp-itebd -H "lattice:1.2*H_xx+0.6*H_x+1.2*H_yy+0*H_y+4*H_zz+0.2*H_z" -w psi.t0.0 -o psi -t 0.1 -n 40

array=(0.0    0.1    0.2    0.3    0.4    0.5    0.6    0.7    0.8    0.9    1.0    1.1    1.2    1.3    1.4    1.5    1.6    1.7    1.8    1.9    2.0    2.1    2.2    2.3    2.4    2.5    2.6    2.7    2.8    2.9    3.0    3.1    3.2    3.3    3.4    3.5    3.6    3.7    3.8    3.9    4.0)

for(( i=0;i<=40;i++ ))
do
    cd ~/Wen/xxz/itebd/

    mp-iexpectation-cross psi.t${array[&Si]} lattice:"Sx(0)" psi.t${array[&Si]} > Sx.t${array[&Si]} 

    t=${array[i]}
    echo "==================="  ${t}
    ##显示输出的delta数值
    mkdir  -p  data/${i}
    ##创建一个里面是26个子文件夹的data文件夹
    echo ${t}   > data/${i}/t.txt

    cp -r psi.t${array[&Si]} data/${i}/psi
    cp -r Sx.t${array[&Si]} data/${i}/Sx.txt

done