基于簇模型计算Ni-Pt(1:1)合金表面吸附原子氢时自旋多重度的设置问题

SantoUranus

仔细阅读了http://sobereva.com/540，例子中Ag团簇很容易设置成闭壳层体系，我自己的体系是Ni-Pt(1:1)合金，我建立了合适的团簇模型作为吸附表面（S），吸附物为原子氢（A）。现在在设置自旋多重度时遇到了问题，Ni-Pt这边成键情况略复杂，我不清楚自旋多重度该怎么设，暂时认为应该把各种自旋多重度从低到高都试一遍（2,4,6,...,N）? 假设自旋多重度为n时得到 S+A 的最低的能量，那我随后进行 S 的单点能计算时需要重新再把所有可能的自旋多重度（1,3,5,...,N-1）都试一遍么？还是应该直接用n-1，哪怕对应得到的能量不是最低的?


P.S. 可能有人会问为什么要用Gaussian而不用VASP等第一性原理计算软件算周期性体系，我这里选择Gaussian是因为要考虑溶剂效应，而我使用的溶剂是离子液体，看了http://sobereva.com/431认为Gaussian中SMD对离子液体的隐形溶剂模型描述还是比较好的，相比较下用VASPsol的话只是改改静态介电常数感觉不太好。

SantoUranus

自己顶一下