低温时的热容问题

[PG13s]

您好，感谢您一直更新软件，最新的版本发现一个问题：在测试CrI3的海森堡模型时，当温度低于25K，capacity反而会突然增大而不是趋于常数 ，不知道是什么原因造成的？

[golddoushi]

@PG13s
您好，感谢您的关注，我刚才测试了最新版，我这边结果是低于25K capacity逐渐趋近于零：
10.000 1.385354 0.033079 -77.998282 0.000531 0.999841
42.000 0.165011 0.052923 -29.421785 0.000779 0.174717
14.000 1.333637 0.033672 -73.605556 0.000547 0.999621
20.000 1.246949 0.034756 -66.683255 0.000587 0.998945
18.000 1.277217 0.034373 -69.041022 0.000570 0.999243
12.000 1.359972 0.033373 -75.822451 0.000535 0.999751
46.000 0.104372 0.048676 -24.093014 0.000546 0.209656
48.000 0.092841 0.046695 -21.994315 0.000473 0.231930
36.000 0.857736 0.042064 -43.402796 0.000978 0.968894
30.000 1.056193 0.037815 -53.526571 0.000722 0.994219
16.000 1.306049 0.034015 -71.347898 0.000560 0.999463
44.000 0.124345 0.050785 -26.526954 0.000643 0.182767
22.000 1.214645 0.035211 -64.247982 0.000604 0.998531
24.000 1.180100 0.035723 -61.743416 0.000628 0.997959
40.000 0.311218 0.053829 -33.185776 0.001155 0.320344
32.000 1.004058 0.038808 -50.462439 0.000774 0.991348
38.000 0.715362 0.045744 -38.950197 0.001245 0.879254
26.000 1.143060 0.036285 -59.135180 0.000653 0.997207
34.000 0.940461 0.040163 -47.116859 0.000857 0.985096
28.000 1.102094 0.036978 -56.410670 0.000682 0.996067
50.000 0.083784 0.044860 -20.179580 0.000416 0.230852

[PG13s]

感谢您的回复，我用的网盘里05.02update的samples\CrI3With2NNCoupling做的计算，参数没改动，算的heisenberg model，从1K到5K，结果如下，可以看到capacity这列随温度下降，呈上升趋势，提高nthermal与nsweep依然有这个问题，麻烦您百忙之中抽空看一下：
#Temp #Spin #Susc #energy #capacity #Binder cumulante
1.000 1.492994 0.020968 -160.292231 0.018065 0.999993
2.000 1.486235 0.020553 -159.289221 0.005976 0.999987
3.000 1.479480 0.020380 -158.284449 0.003742 0.999977
4.000 1.472666 0.020314 -157.276340 0.002959 0.999963
5.000 1.465757 0.020311 -156.265249 0.002563 0.999944

另外还有个小疑惑，如果把orbital list里的ID 0的spin改成0.5，ID 1的spin保持1.5，输出的result一栏spin随温度直接从0.5开始下降，请问这里怎么还原成晶格单胞的总自旋，直接乘4是否可行

[golddoushi]

@PG13s
我这边1-5K的数据看起来也挺正常的。近两天我把exe程序更新一下，平均自旋的问题大概增加一个UI来选择具体输出的自旋吧
#Temp #Spin #Susc #energy #capacity #Binder cumulante
1.000 1.489337 0.031875 -87.392472 0.000492 0.999999
5.000 1.444979 0.032407 -83.294346 0.000509 0.999967
4.000 1.456320 0.032283 -84.329267 0.000504 0.999979
3.000 1.467509 0.032139 -85.358852 0.000492 0.999989
2.000 1.478506 0.032010 -86.380293 0.000500 0.999995
11.000 1.372795 0.033221 -76.913455 0.000532 0.999795
7.000 1.421712 0.032677 -81.199584 0.000518 0.999931
8.000 1.409746 0.032827 -80.138530 0.000517 0.999906
10.000 1.385354 0.033079 -77.998282 0.000531 0.999841
6.000 1.433460 0.032532 -82.250892 0.000511 0.999951
9.000 1.397694 0.032939 -79.074600 0.000523 0.999874
13.000 1.346950 0.033517 -74.718671 0.000539 0.999697
12.000 1.359972 0.033373 -75.822451 0.000535 0.999751
15.000 1.319873 0.033862 -72.478680 0.000556 0.999546
14.000 1.333637 0.033672 -73.605556 0.000547 0.999621
18.000 1.277217 0.034373 -69.041022 0.000570 0.999243
17.000 1.291942 0.034170 -70.205923 0.000561 0.999358
16.000 1.306049 0.034015 -71.347898 0.000560 0.999463

[PG13s]

万分感谢！

[PG13s]

您好，再次打扰您，我拿exe的522版本做了默认CrI3With2NNCoupling的测试，仍然有比热反向变大的问题，与您之前计算的结果有出入，结果如下：
#Temp # # #Susc #energy #capacity #Binder cumulante
1.000 1.489074 1.489074 0.032660 -87.364138 0.001738 0.999998
2.000 1.478298 1.478298 0.032318 -86.357445 0.000804 0.999994
3.000 1.467376 1.467376 0.032269 -85.344385 0.000630 0.999988
4.000 1.456247 1.456247 0.032337 -84.323126 0.000571 0.999978
5.000 1.444969 1.444969 0.032413 -83.294889 0.000545 0.999966
小木虫上也有人遇到类似的问题：http://muchong.com/t-4803656-1
麻烦您抽空再看一下，谢谢！

[golddoushi]

@PG13s
您好，我这里是完全正常的。对于CrI3这么大的J，在1K的时候按照默认设置autocorrelation可能会非常大，造成显著统计误差。所以不同机器上可能得出相悖的结果。
现在程序做了一些更改，降低了内存消耗，同时在result.txt中增加了一栏auto corr，您可以对比一下低温区与近相变区的auto-corr，看看是不是这个原因造成的。
目前增大tau值（nthermal，nsweep行的最后一个）可以降低一些auto-corr，但是对于近基态的效果依然不太好。所以不建议计算近基态。
以后可能会考虑改变算法增加基态的取样效率。

[PG13s]

您好，感谢您这么快的更新，想请教一下tau值是什么参数，另外新的524版本exe有点问题，目前一点StartMC就不运行了，麻烦您抽空看一下，谢谢！

[golddoushi]

@PG13s
您好，tau表示两个相邻采样点之间的MC更新步数，tau越大则采样点间隔越大，采样点之间的关联越小，最终的统计结果越好。特别的，tau=0表示采样间隔等于超胞的格点总数，也就是文献中常用的sweep lattice。
另外524在特定环境下可能是有找不到库的问题，近期会修复。

[PG13s]

明白了，谢谢您！