gaussian溶剂化计算求教

xiaomi7777

  高斯新手求教各位大神： 最近在计算柠檬酸和硫酸根两个分子的相互作用，想看一下两个分子在溶剂化的条件下会有什么相互的作用，看帖子说计算溶剂化效应：

“溶剂化的计算是单点计算，能量计算的正确与否由分子构型的合理性程度来决定的，所以能量计算是在优化好的结构上进行的，没有结构的优化（opt）和频率（freq）的输出。”根据这个观点我在气相条件下优化了两个分子的构型以及两个分子在气相条件下相互作用的构型，然后再用溶剂化效应做了单点的处理。可是从得到的图片看来的话，气相条件下，硫酸根会与柠檬酸的氢离子相互作用，形成硫酸氢跟。那么想请教各位大神，高斯溶剂化效应能否直接优化两个分子然后看两个分子在溶剂中有什么样的相互作用呢？

同时小弟是做光谱的，如果不能这样做的话，溶剂化效应如果想考虑溶剂化条件下两个分子的相互作用该怎么做呢？溶剂化效应下计算频率是怎样做的呢？

附件为我得到的气相条件下的两个分子作用的输出文件。做溶剂化效应时所添加的关键词如下：
%mem=8GB
%nproc=6
# B3LYP/6-311++G* scrf=(solvent=water,pcm)  iop(5/13=1) guess=read geom=allcheck

烦请各位大神指点。。

kevin

我想请问一个问题：#气相条件下，硫酸根会与柠檬酸的氢离子相互作用，形成硫酸氢跟#，我比较好奇，气相条件下，质子酸还能给出质子么？谢谢啦

xiaomi7777

sobereva 发表于 2016-4-7 16:44
你提到的那个观点并不确切。
许多情况下，气相和溶剂下优化结果基本一致，这时候可以用气相优化和频率计算 ...

好的。。谢谢sob老师的讲解。。

sobereva

你提到的那个观点并不确切。
许多情况下，气相和溶剂下优化结果基本一致，这时候可以用气相优化和频率计算的结果作为近似。但是很多情况，特别是对于带电荷的情况，有可能气相和溶剂下优化出的结构是明显不同的，这时候应当在溶剂下优化、计算频率，并基于此结构计算光谱、相互作用能等属性。

绝对不要用iop(5/13=1)，原因见此文文末
解决SCF不收敛问题的方法
http://sobereva.com/61