以晶体结构文件cif为基准计算基态激发态的相关问题

seantan521

大家好，我现在有问题如下：
我有XRD得到的晶体结构文件cif，导入chem3D转换成gaussian input文件gif格式。

然后， 如果我想知道基态的最小能量S0（错了，应该是此晶体结构的能量）？是不是应该用这个gif的input文件做一下energy计算还是freq计算？
energy计算方法：# b3lyp/6-311g(d) charge:0 Mult:1
freq计算方法：# freq b3lyp/6-311g(d) charge:0 Mult:1

然后， 如果我想继续优化S1和T1呢？是用cif转换后的gif文件直接做S1和T1的opt td，还是用energy或者freq计算后的chk文件做S1和T1的opt td？

S1优化方法：# opt td=(singlets) b3lyp/6-311g(d) charge:0 Mult:1
T1优化方法：# opt td=(triplets) b3lyp/6-311g(d) charge:0 Mult:1

最后，如果算出的晶体能量和opt基态得到的最小能量不一致的话，应该用哪一个呀？

谢谢大家解惑答疑。

h840473807

如果是分子晶体的话，建议还是以晶体中的一个分子为初始结构全优化，毕竟气态分子和晶体中的分子结构是不同的，后者有大量的分子间相互作用。
想对比基态和激发态的能量用单点能就可以了，不用freq计算，计算非常耗时，恐怕也得不到什么更有用的信息。

seantan521

h840473807 发表于 2018-4-5 22:39
如果是分子晶体的话，建议还是以晶体中的一个分子为初始结构全优化，毕竟气态分子和晶体中的分子结构是不同 ...

你的意思是，用这个晶体结构数据，重新进行一次基态优化？
这样的话，与直接用chemdraw结构式优化得到的结果不一样么？

h840473807

seantan521 发表于 2018-4-5 22:56
你的意思是，用这个晶体结构数据，重新进行一次基态优化？
这样的话，与直接用chemdraw结构式优化得到的 ...

也可以

望月凝香

seantan521 发表于 2018-4-5 22:56
你的意思是，用这个晶体结构数据，重新进行一次基态优化？
这样的话，与直接用chemdraw结构式优化得到的 ...

晶体结构更可能接近真实构象。从chemdraw出发容易陷在局部最小点。