|
1、请问在使用ABCluster搜索Mg10团簇的时候,我用atom通过abcinp.exe Mg 1 Gupta 10 100 100 4 100 10 Mg得到了84个Mg10结构能量(如下1,Mg.out),我看【Mg-LM】中的gjf文件显示的是#B3LYP/6-31+G(d) OPT FREQ(如图2),那我是不是可以直接使用能量最低的0.gif结构在B3LYP/def2svp上进行优化之后做吸附氢分子计算?还是说我要将产生输出的84个结构全部用Guassian在B3LYP/def2svp下优化比较出最低能量后再进行吸附氢分子计算?
我计算吸附氢分子的吸附能是在B3LYP/def2svp下进行优化,在B3LYP/def2tzvp下计算单点能得到的能量。
2、想问一下,图一Mg.out中输出的能量是不是就是在B3LYP/6-31+G(d)下得到的能量?
3、我看手册中说Gupta势不能用于计算小团簇,不知道我使用这个abcinp.exe Mg 1 Gupta 10 100 100 4 100 10 Mg计算对不对?如果不对的话,应该如何得到能量较低的小尺寸的金属原子团簇?
Mg.out
-- Results Report --
* 84 LMs will be saved in [ Mg-LM ].
# Energy Match-RMSD
0 -3.11386947 0.00000000
1 -3.07135338 0.77107043
2 -3.06817892 1.50210133
3 -3.06706532 1.66447736
4 -3.06678332 1.31504054
5 -3.06672519 1.36534950
6 -3.06424618 1.56697178
7 -3.06413506 1.59781632
8 -3.06363650 1.80416881
9 -3.06320488 1.26842174
...
...
79 -1.95390931 3.08333411
80 -1.94419472 2.81798429
81 -1.90561068 2.16965585
82 -1.81428867 2.89766949
83 -1.71360209 3.34466601
* Final Global Minimal Energy: -3.11386947
* The Global Minimum is saved as: [ Mg.(gjf/xyz) ]
【Mg-LM】中gjf文件:
#B3LYP/6-31+G(d) OPT FREQ
10, Energy = -3.11386947 au
0 1
Mg 3.30202068 -1.33001402 1.02378980
Mg -0.00918140 -2.03397818 1.85342613
Mg -0.32616611 -0.41371678 -3.66646590
Mg 1.36579014 1.07912586 2.84190930
Mg -2.13037450 -1.49816703 -0.88865905
Mg -1.98284319 0.86213249 1.70222737
Mg 0.85108742 0.55123648 -0.55065842
Mg -2.12003449 1.95975897 -1.66463566
Mg 1.12343167 -2.79980288 -1.41607587
Mg -0.07373023 3.62342507 0.76514230
|
|