如何在Gaussian连续计算的第二步取消之前的Modredundant限制

liyuanhe211

设定一计算，希望首先锁定一些参数进行优化，随后放开这些限制继续优化。

1. %chk=***
   
2. #p opt=modredundant
   
3. 
   
4. title
   
5. 
   
6. 结构
   
7. 
   
8. B 14 2 F
   
9. 
   
10. --Link1--
    
11. %chk=***
    
12. #p opt=(calcfc,ts) freq geom=allcheck
    
13. 
    

复制代码

这时发现第二步计算中，第一步对结构的锁定仍然存在

1. <div>                           ----------------------------</div><div>                           !    Initial Parameters    !</div><div>                           ! (Angstroms and Degrees)  !</div><div> --------------------------                            --------------------------</div><div> ! Name  Definition              Value          Derivative Info.                !</div><div> --------------------------------------------------------------------------------</div><div> ! R1    R(1,2)                  1.4726         calculate D2E/DX2 analytically  !</div><div> ! R2    R(1,8)                  1.5501         calculate D2E/DX2 analytically  !</div><div> ! R3    R(1,23)                 1.0977         calculate D2E/DX2 analytically  !</div><div> ! R4    R(1,24)                 1.0966         calculate D2E/DX2 analytically  !</div><div> ! R5    R(2,3)                  1.2472         calculate D2E/DX2 analytically  !</div><div><font color="#ff0000"> ! R6    R(2,14)                 2.2            frozen, calculate D2E/DX2 analyt!</font></div><div> ! R7    R(3,4)                  1.4684         calculate D2E/DX2 analytically  !</div>

复制代码

根据Gaussian说明书的解释：

Checkpoint may be combined with the ModRedundant option if you want to retrieve and alter the molecule specification in a checkpoint file using redundant internal coordinate-style modifications.

似乎应当只有加入geom=modredundant的时候才提取相应modredundant声明，但这里不是这样。并且尝试在第二步加入geom=modredundant没有变化。

如何解决？

pigfool

opt=(calcfc,ts,nofreeze)

liyuanhe211

pigfool 发表于 2015-10-9 10:01
opt=(calcfc,ts,nofreeze)

啊可以了~~
谢谢。