Po原子衰变为Pb原子后的动力学过程该如何模拟

冬天一起堆雪人

11个Po原子的团簇，其中一个Po原子衰变为Pb原子后具备一个动能（初速度），想要模拟PbPo10的动力学过程该如何模拟？用什么软件？

之前用Gaussian中的BOMD进行了模拟，一直没有结果，不知道是不是方向不对，求助各位大佬。（以下是Gaussian模拟的输入文件）

%nprocs=32
%mem=20GB
#p PBE1PBE/def2svpp bomd(stepsize=2000,maxpoints=2000,ReadVelocity) nosymm

Title Card Required

-2 1
Po                -1.99457800   -0.17343000    1.96452800
Po                 1.48296500   -2.22757600    1.94291900
Po                 0.57396100   -3.26485100   -0.59563800
Po                -0.40388100   -0.69905900   -1.67202500
Po                -3.19989900   -1.23435000   -2.13794700
Po                 2.53384200    2.99542800   -0.61300300
Po                -0.35745700    3.14383400   -0.82261600
Po                -1.38479400    2.66933100    1.85894200
Po                 4.03248400   -1.11030700    1.01162100
Po                 3.14411900    0.24680600   -1.39089700
Pb                -4.42676000   -0.34582600    0.45411600

0
0.0    0.0     0.0
0.0    0.0     0.0
0.0    0.0     0.0
0.0    0.0     0.0
0.0    0.0     0.0
0.0    0.0     0.0
0.0    0.0     0.0
0.0    0.0     0.0
0.0    0.0     0.0
0.0    0.0     0.0
51400000000000    24700000000000     7890000000000

wzkchem5

建议用ORCA。一方面ORCA的BOMD做的比高斯好，另一方面ORCA对相对论方法的支持比较好，这么重的体系用赝势未必靠谱

冬天一起堆雪人

wzkchem5 发表于 2024-7-3 14:35
建议用ORCA。一方面ORCA的BOMD做的比高斯好，另一方面ORCA对相对论方法的支持比较好，这么重的体系用赝势未 ...

感谢您的回复，基组您有什么建议吗？

sobereva

冬天一起堆雪人 发表于 2024-7-3 14:51
感谢您的回复，基组您有什么建议吗？

可以先尝试def2-SVP赝势基组。如果遇到实际问题或者结果不理想再说

冬天一起堆雪人

sobereva 发表于 2024-7-3 19:49
可以先尝试def2-SVP赝势基组。如果遇到实际问题或者结果不理想再说

好的，谢谢。