crystal计算过渡金属基组选择问题 - 计算化学公社

我最近在用crystal计算过渡金属Fe slab 的吸附问题。我在crystal里优化好的结构带入到vasp里发现受力很大,

这说明我没有选择好过渡金属Fe的高斯基组。因为我之前算过很多过渡金属氧化物, crystal优化的结构在vasp里面

一般都比较稳定。通常我也是用这种方法测试一下自己自己基组是否选择的合理。我看过一些资料，里面说用高斯基组

计算金属体系比较麻烦，需要选择比较弥散的基组才行，但是在周期性体系里exponent不要小于0.1。下面就是我用的基组

我感觉我加了弥散的exponent和 f g 轨道的极化函数，而且Zeta 数目也较多了，可是还有问题。所以想请各位专家帮忙看

一下，谢谢。（我个人认为我对于crystal code 使用没有任何问题，因为我之前算过一些比价复杂的过渡金属氧化物都没问

题，只是对于这种金属体系没有经验。）

下面就是我选择的基组 Fe_pob_TZVP_2012

26 14

0 0 8 2.0 1.0

300784.84637 0.00022806273096

45088.970557 0.00176817887610

10262.516317 0.00919270834900

2905.2897293 0.03735549580700

946.11487137 0.12151108426000

339.87832894 0.28818881468000

131.94425588 0.41126612677000

52.111494077 0.21518583573000

0 0 4 2.0 1.0

329.48839267 -0.02474521647700

101.92332739 -0.11683089050000

16.240462745 0.55293621136000

6.8840675801 0.53601640182000

0 0 2 2.0 1.0

10.470693782 -0.22912708577000

1.7360039648 0.71159319984000

0 0 1 2.0 1.0

1.7565166800 1.00000000000000

0 0 1 0.0 1.0

0.77548354 1.00000000000000

0 0 1 0.0 1.0

0.1058918100 1.00000000000000

0 2 6 6.0 1.0

1585.3959970 0.00237939601790

375.38006499 0.01925315475500

120.31816501 0.09002183653600

44.788749031 0.25798172356000

17.829278584 0.41492649744000

7.2247153786 0.24207474784000

0 2 3 6.0 1.0

28.143219756 -0.02904175515200

3.8743241412 0.55312260343000

1.5410752281 0.96771136842000

0 2 1 0.0 1.0

0.9336590800 1.00000000000000

0 2 1 0.0 1.0

0.1999999900 1.00000000000000

0 3 4 6.0 1.0

61.996675034 0.01197197225500

17.873732552 0.07321013541000

6.2744782934 0.23103094314000

2.3552337175 0.39910706494000

0 3 1 0.0 1.0

0.8417158200 1.00000000000000

0 3 1 0.0 1.0

0.2611100300 1.00000000000000

0 4 1 0.0 1.0

1.5980000000 1.00000000000000