Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

相同的交换参数在新版本mcsolver3.0模拟得不到有规律的结果。 #33

Open
pengxb1228 opened this issue Apr 19, 2024 · 5 comments

Comments

@pengxb1228
Copy link

您好,作者;
我采用了我的体系获得的交换参数,在mcsolver1.2版本模拟可以得到合理的结果,但是新版本mcsolver3.0在linux系统上模拟获得的磁矩与温度的关系是离散化的,没有演变的规律。请问是什么原因造成的?以
下贴出了我的输入文件,麻烦您帮忙解答一下,谢谢。
This is mcsolver's save file, version: 3.0
Lattice:
1.000000000 0.000000000 0.000000000
0.000000000 2.000000000 0.000000000
0.000000000 0.000000000 1.000000000
Supercell used in MC simulations:
30 30 1
Orbitals in cell:
2
Positions, initial spin states and onsite-anisotropy of every orbital:
orb 0: type 0 spin 2.000000000 pos [0.500000000 0.166666000 0.000000000] Dx 0.000000000 Dy 0.000000000 Dz -0.012000000 h 0.000000000
orb 1: type 0 spin 2.000000000 pos [0.000000000 0.666666667 0.000000000] Dx 0.000000000 Dy 0.000000000 Dz -0.012000000 h 0.000000000
Bonds:
11
id, source, target, overLat, exchange matrix elements of each bond:
bond 0: Jx -159.210000000 Jy -159.210000000 Jz -160.620000000 Jxy 0.000000000 Jxz 0.000000000 Jyz 0.000000000 Jyx 0.000000000 Jzx 0.000000000 Jzy 0.000000000 orb 1 to orb 1 over [1 0 0]
bond 1: Jx -159.210000000 Jy -159.210000000 Jz -160.620000000 Jxy 0.000000000 Jxz 0.000000000 Jyz 0.000000000 Jyx 0.000000000 Jzx 0.000000000 Jzy 0.000000000 orb 0 to orb 0 over [1 0 0]
bond 2: Jx -159.210000000 Jy -159.210000000 Jz -160.620000000 Jxy 0.000000000 Jxz 0.000000000 Jyz 0.000000000 Jyx 0.000000000 Jzx 0.000000000 Jzy 0.000000000 orb 1 to orb 0 over [0 1 0]
bond 3: Jx -159.210000000 Jy -159.210000000 Jz -160.620000000 Jxy 0.000000000 Jxz 0.000000000 Jyz 0.000000000 Jyx 0.000000000 Jzx 0.000000000 Jzy 0.000000000 orb 0 to orb 1 over [1 0 0]
bond 4: Jx -159.210000000 Jy -159.210000000 Jz -160.620000000 Jxy 0.000000000 Jxz 0.000000000 Jyz 0.000000000 Jyx 0.000000000 Jzx 0.000000000 Jzy 0.000000000 orb 0 to orb 1 over [0 0 0]
bond 5: Jx -159.210000000 Jy -159.210000000 Jz -160.620000000 Jxy 0.000000000 Jxz 0.000000000 Jyz 0.000000000 Jyx 0.000000000 Jzx 0.000000000 Jzy 0.000000000 orb 1 to orb 0 over [-1 1 0]
bond 6: Jx -10.140000000 Jy -10.140000000 Jz -10.610000000 Jxy 0.000000000 Jxz 0.000000000 Jyz 0.000000000 Jyx 0.000000000 Jzx 0.000000000 Jzy 0.000000000 orb 0 to orb 0 over [0 1 0]
bond 7: Jx -10.140000000 Jy -10.140000000 Jz -10.610000000 Jxy 0.000000000 Jxz 0.000000000 Jyz 0.000000000 Jyx 0.000000000 Jzx 0.000000000 Jzy 0.000000000 orb 0 to orb 1 over [2 -1 0]
bond 8: Jx -10.140000000 Jy -10.140000000 Jz -10.610000000 Jxy 0.000000000 Jxz 0.000000000 Jyz 0.000000000 Jyx 0.000000000 Jzx 0.000000000 Jzy 0.000000000 orb 0 to orb 1 over [-1 0 0]
bond 9: Jx -10.140000000 Jy -10.140000000 Jz -10.610000000 Jxy 0.000000000 Jxz 0.000000000 Jyz 0.000000000 Jyx 0.000000000 Jzx 0.000000000 Jzy 0.000000000 orb 0 to orb 1 over [-1 -1 0]
bond 10: Jx -10.140000000 Jy -10.140000000 Jz -10.610000000 Jxy 0.000000000 Jxz 0.000000000 Jyz 0.000000000 Jyx 0.000000000 Jzx 0.000000000 Jzy 0.000000000 orb 0 to orb 1 over [2 0 0]
Temperature scanning region:
Tmin 400.000000000 Tmax 1000.000000000 nT 50
Field scanning region (in unit 1.48872 T, only if Kelvin and uB is used for energy and spin):
Hmin 0.000000000 Hmax 0.100000000 nH 1
Dipole long-range coupling:
alpha 0.000000
Measurement:
measure the correlation function between orb0 and orb0 over [0 0 0]
Supergroup
OrbGroup:1
Supergroup
group0 orb0-orb0

  Topological section      <<<

LocalCircuit per cell: 0 (set to 0 to skip the calc. for topo. Q)

End of Topological section <<<
Distribution output frame: 0
Sweeps for thermalization and statistics, and relaxiation step for each sweep:
80000 640000 1
XAxis type:
T
Model type:
Heisenberg
Algorithm:
Metropolis
Ncores:
40

@golddoushi
Copy link
Owner

您好,我用百度网盘里最新的23.12.06版本跑您这个参数可以得到比较平滑的比热和磁矩演化,居里温度大概在750。您上面的设置主要问题在于统计采样数目过少,您可以尝试把Algorithm换成Wolff算法。或者保持Metropolis算法,但是把tau设为0(等价于tau=orbital总数=1800)

@pengxb1228
Copy link
Author

pengxb1228 commented Apr 20, 2024 via email

@pengxb1228
Copy link
Author

您好,刘老师,再次打扰了。
我发现,现在mcsolver官网提供百度网盘的链接 https://pan.baidu.com/s/1EaDqOOdB7AP9WXrwEIEaxQ不可访问. 另外我把tau改为了0,模拟出来的结果仍然是离散的,我的模拟是在Linux系统上,直接通过pip install mcsolver安装的,是不是我直接pip安装版本出现的问题?麻烦您帮忙解答一下,谢谢。祝您生活愉快!

@pengxb1228
Copy link
Author

pengxb1228 commented Apr 20, 2024 via email

@golddoushi
Copy link
Owner

您好,感谢您的测试,网盘链接超过5年失效了,现在修复了链接:https://pan.baidu.com/s/1kRtvPFXDBANga1StRxG7Wg?pwd=1234

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants