【求助】geom=allcheck 到底读取的哪里的构型？
我在计算过程中首先第一个工作分两步：
第一步，做气相下的优化，用了％check=1.chk
第二步，做溶液下的优化，但是这一步优化过程非常慢，运行了一段时间，然后就停掉了。
接着做了第二个工作，在溶液中进行频率计算，用了％check=1.chk
分子说明部分用了geom=allcheck，那么这一步读取的构型是哪个？
是气相下优化好了的构象，还是溶液中优化到半路的构型啊？
请高手赐教！谢谢！万分感激！
[ Last edited by yuhuobuku on 2009-2-20 at 09:53 ]
返回小木虫查看更多
今日热帖
京公网安备 11010802022153号
你这个应该是把原来的1.chk文件覆盖了.
溶液下的没有优化完,做频率就没什么意思了吧?
我主要是想得到Variational C-PCM results
=========================
<psi(0)| H |psi(0)> (a.u.) = -733.133752
<psi(0)|H+V(0)/2|psi(0)> (a.u.) = -733.193523
<psi(0)|H+V(f)/2|psi(0)> (a.u.) = -733.196729
<psi(f)| H |psi(f)> (a.u.) = -733.129748
<psi(f)|H+V(f)/2|psi(f)> (a.u.) = -733.196756
Total free energy in solution:
with all non electrostatic terms (a.u.) = -733.192171
--------------------------------------------------------------------
(Unpolarized solute)-Solvent (kcal/mol) = -37.51
(Polarized solute)-Solvent (kcal/mol) = -42.05
Solute polarization (kcal/mol) = 2.51
Total electrostatic (kcal/mol) = -39.54
--------------------------------------------------------------------
Cavitation energy (kcal/mol) = 27.87
Dispersion energy (kcal/mol) = -29.66
Repulsion energy (kcal/mol) = 4.66
Total non electrostatic (kcal/mol) = 2.88
DeltaG (solv) (kcal/mol) = -36.66
--------------------------------------------------------------------
这个结果，用的SCF=TIHGT，
版主，第二步读取的哪个构象呢？ 是气相下优化好的构象还是液相下优化到半路的构象呢？
楼主，在进行第二步溶剂下的结构优化的时候，chk文件名是什么设置的呢？
要是还是采用％check=1.chk，那么你在进行第二步的时候，生成的check文件就把你第一步的check文件给覆盖了。在你后来读取check文件的时候，你读取的应该的在你第二步在液相中没有弄好的结构式的信息。
要是你在进行第二步的时候，给出了另外的check文件名，例如：％check=2.chk，那么，在进行第二步，溶剂条件下优化的时候，就自己会生成check文件，你以前气相的check文件保留着，那么你后来用％check=1.chk来读取结构式，应该是读取的第一步气相下优化好的结构信息。
明白一些了！非常感谢