版本选择和参数设置

[e-xx666]

@golddoushi 尊敬的老师您好！非常感谢您提供这么好用的软件！我刚开始学习您的软件，有许多地方不太懂，真的非常抱歉打扰您啦！这是我在用海森堡模型模拟二维方形铁磁材料居里温度时遇到以下问题：
1.我使用的是您的第二个版本（mcsolver 20.05.21）我需要使用最新版本（mcsolver 21.03.08）吗？
2.我的材料的易磁轴是面外，请问用海森堡模型可以忽略磁各向异性吗？我平时只算了一个J，是否只输入Jz一个值，其他都设置为0就可以了？还是说海森堡模型的这些J的分量也都要有数值吗？我再去想一想Jx，Jz怎么求。如果我的材料易磁轴在面内，不考虑磁各项异性影响大吗？
3.我用VASP计算了不加SOC的超胞铁磁和反铁磁的能量，求出了ΔE，在求J的时，我发现有的文献公式说J前面有1/2，有的文献却说没有，文献中S的取值也各不相同，有的文献取单个原子磁矩，有的取3/2，有的直接取1。我比较笨，有点懵了，我再想想。有文献在CrI3做图时磁矩直接用单个Cr的磁矩3随温度变化，我们的mcslover提供例子是单个Cr用1.5，请问老师，这里S和J只是单纯的倍数关系吗？我的材料单个磁性原子是4，我应该怎么设置ID 0的spin更加准确呢？
4.我重复了一些其他文献最近邻J，我得到的居里温度和文献有一些差别，可能因为VASP参数设置不一样，或者大家用的蒙特卡洛模拟器不一样。我设置初始温度为0.1K，在低温范围的比热有时会有凸起，这是因为我步数设置少还是某些材料体系本身导致的，还是我用的不是最新版本的mcslover？
感谢老师百忙之中抽时间看我的留言！我才疏学浅，今后一定好好学习！感谢老师提供如此棒的软件！谢谢老师！

[golddoushi]

您好，

1. 建议您使用新版本。新版本操作方式与老版相同，但是修复了一些问题
2. 如果您的材料是三维的，并且各向异性不是太强，可以忽略各向异性。但对于二维一维系统，忽略各向异性的海森堡模型无法给出正确的答案。设置各向同性的海森堡模型，需要把Jx，Jy，Jz都填为同样的值，都需要设置。各个方向的J的计算方法可以参考我们的论文 ：dio: 10.1021/acs.jpclett.0c01911
3. MCSOLVER采用了无double counting的算法，所以J前面没有1/2 。另外，因为我们没有考虑自旋的量子涨落，所以S可以随意设置，比如3μB的Cr原子的S，可以设置S=1，也可以设置S=3/2，也可以设置为3，等等。但是注意你采用的S必须前后一致，换句话说，用VASP计算的能量差ΔE与S算J的时候，以及设置软件的S与J算Tc的时候，用到的S必须是相同的。
4. 蒙特卡洛算法不适用于计算近基态的性质，所以建议不要计算太低的温度。极低温的反常比热峰有可能是数值计算的误差导致的。

祝您学习生活一切顺利，科研成功大丰收！

[e-xx666]

非常感谢老师提供的宝贵指导！谢谢老师耐心细致的解答！我再去试一试，真的非常感谢！也祝老师生活工作一切顺利！ 😊

[e-xx666]

@golddoushi 尊敬的老师您好！我又来请教您问题了，感谢您上次的耐心解答！
1.二维材料的海森堡模型，我们听从了您的建议，分别求了（各向异性）海森堡模型的JxJyJz,用Metropolis算法（tau=0）算出CrI3的居里温度为41K，很接近实验值，谢谢您提供的宝贵建议！我们现在想用wolff算法，但是在tau设置上遇到一点问题。您提到过tau表示两个相邻采样点之间的MC更新步数，tau越大，采样点间隔越大，采样点之间的关联越小，最终的统计结果越好，请问老师，用wolff算法时tau设置为多少比较好？（tau=0运算的太慢了，tau=1感觉又太小了）
2.二维材料的伊辛模型：文献里没有说过各向异性，只有一个J。我的问题是：若易磁轴是面外，我们不加SOC求的J设为Jz，此时Jx,Jy设为0吗？若易磁轴在面内，比如LaCl，我们Jx,Jy,Jz怎么设呢?（伊辛模型这个我们可能会用，希望老师能提供一点建议，谢谢老师！）
谢谢老师提供如此好的软件！希望老师有时间可以帮帮我。万分感谢！