多类型原子时，Tc模拟出现错误

[1638927697]

老师您好，我在尝试两层不同类型的2D材料磁性Tc模拟的时候，计算得到的热容是发散的，老师您可以帮我看看是什么设置出现了问题吗，谢谢老师。
参数设置如下：
This is mcsolver's save file, version: 3.0
Lattice:
1.000000000 0.000000000 0.000000000
-0.500000000 0.866000000 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 0.435 pos [0.000000000 0.000000000 0.400000000] Dx 0.000000000 Dy 0.000000000 Dz 5.6 h 0.000000000
orb 1: type 0 spin 0.37 pos [0.333333333 0.666666667 0.500000000] Dx 0.000000000 Dy 0.000000000 Dz 5.6 h 0.000000000
Bonds:
17
id, source, target, overLat, exchange matrix elements of each bond:
bond 0: Jx -1494.300000000 Jy -1494.300000000 Jz -1494.300000000 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 1: Jx -1494.300000000 Jy -1494.300000000 Jz -1494.300000000 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 1 0]
bond 2: Jx -1494.300000000 Jy -1494.300000000 Jz -1494.300000000 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 3: Jx -1494.300000000 Jy -1494.300000000 Jz -1494.300000000 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 4: Jx -1494.300000000 Jy -1494.300000000 Jz -1494.300000000 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 -1 0]
bond 5: Jx -1494.300000000 Jy -1494.300000000 Jz -1494.300000000 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 6: Jx -1302.800000000 Jy -1302.800000000 Jz -1302.800000000 Jxy 0.000000000 Jxz 0.000000000 Jyz 0.000000000 Jyx 0.000000000 Jzx 0.000000000 Jzy 0.000000000 orb 1 to orb 1 over [0 1 0]
bond 7: Jx -1302.800000000 Jy -1302.800000000 Jz -1302.8300000000 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 1 0]
bond 8: Jx -1302.800000000 Jy -1302.800000000 Jz -1302.800000000 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 9: Jx -1302.800000000 Jy -1302.800000000 Jz -1302.800000000 Jxy 0.000000000 Jxz 0.000000000 Jyz 0.000000000 Jyx 0.000000000 Jzx 0.000000000 Jzy 0.000000000 orb 1 to orb 1 over [0 -1 0]
bond 10: Jx -1302.800000000 Jy -1302.800000000 Jz -1302.800000000 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 -1 0]
bond 11: Jx -1302.800000000 Jy -1302.800000000 Jz -1302.800000000 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 12: Jx 732.80 Jy 732.8 Jz 732.8 Jxy 0 Jxz 0 Jyz 0 Jyx 0 Jzx 0 Jzy 0 orb 0 to orb 1 over [0 0 0]
bond 13: Jx 732.80 Jy 732.8 Jz 732.8 Jxy 0 Jxz 0 Jyz 0 Jyx 0 Jzx 0 Jzy 0 orb 0 to orb 1 over [0 -1 0]
bond 14: Jx 732.80 Jy 732.8 Jz 732.8 Jxy 0 Jxz 0 Jyz 0 Jyx 0 Jzx 0 Jzy 0 orb 0 to orb 1 over [-1 -1 0]
bond 15: Jx 732.80 Jy 732.8 Jz 732.8 Jxy 0 Jxz 0 Jyz 0 Jyx 0 Jzx 0 Jzy 0 orb 1 to orb 0 over [0 1 0]
bond 16: Jx 732.80 Jy 732.8 Jz 732.8 Jxy 0 Jxz 0 Jyz 0 Jyx 0 Jzx 0 Jzy 0 orb 1 to orb 0 over [1 1 0]
Temperature scanning region:
Tmin 0.900000000 Tmax 800.000000000 nT 160
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
#group1 orb1-orb1

  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:
50000 100000 1
XAxis type:
T
Model type:
Heisenberg
Algorithm:
Wolff
Ncores:
10

[golddoushi]

温度太低，建议设置最低温度为1000

[1638927697]

温度太低，建议设置最低温度为1000

非常感谢老师您的回复。老师您是认为这里设置的J非常大，Tc应该会很高吗？我这里的spin是比较小的，Tc应该并不高，如果我只模拟单独一层，是可以正常模拟出Tc的，大约在300K左右，但两层一起模拟，就会导致计算的热容是10的20次方左右的数值，所以我比较困惑是什么地方设置错了导致这样的情况。再次感谢老师您的回复。