高斯使用calcall关键词优化过渡态能否输出每一帧的频率？

星斗如盘

我们知道优化过渡态时使用calcall关键字，每一帧结构都会精确计算力常数，有没有办法让输出文件，输出每一帧的频率呢？

liyuanhe211

可以用非标准路径，比如下列输入文件相当于在 HF/STO-3G 级别下用 opt=(TS, CalcAll, NoEigenTest) 优化 H 自由基对 H2 取代的过渡态，但是每一步多输出一次频率和热力学量。

1. 
   
2. #p nonstd
   
3. 1/5=1,10=4,11=1,18=20,38=1/1,3;
   
4. 2/9=110,12=2,17=6,18=5,40=1/2;
   
5. 3/6=3,11=9,16=1,25=1,30=1,71=2,140=1/1,2,3;
   
6. 4//1;
   
7. 5/5=2,38=5/2;
   
8. 8/6=4,10=90,11=11/1;
   
9. 11/6=1,8=1,9=11,15=111,16=1/1,2,10;
   
10. 10/6=1/2;
    
11. 6/7=2,8=2,9=2,10=2,18=1,28=1/1;
    
12. 7/10=1,25=1/1,2,3,16;
    
13. 1/5=1,10=4,11=1,18=20/3(3);
    
14. 2/9=110/2;
    
15. 7/8=1,9=1,25=1,44=-1/16;
    
16. 99//99;
    
17. 2/9=110/2;
    
18. 3/6=3,11=9,16=1,25=1,30=1,71=2,140=1/1,2,3;
    
19. 4/5=5,16=3,69=1/1;
    
20. 5/5=2,38=5/2;
    
21. 8/6=4,10=90,11=11/1;
    
22. 11/6=1,8=1,9=11,15=111,16=1/1,2,10;
    
23. 10/6=1/2;
    
24. 7/10=1,25=1/1,2,3,16;
    
25. 7/8=1,9=1,25=1,44=-1/16;
    
26. 1/5=1,10=4,11=1,18=20/3(-9);
    
27. 2/9=110/2;
    
28. 6/7=2,8=2,9=2,10=2,18=1,19=2,28=1/1;
    
29. 7/8=1,9=1,25=1,44=-1/16;
    
30. 99//99;
    
31. 
    
32. Title Card Required
    
33. 
    
34. 0 2
    
35. H                  0.68980350   -1.59990556   -0.04189500
    
36. H                 -0.11729650   -1.59990556   -0.04189500
    
37. H                 -0.91999650   -1.59990556   -0.04189500
    
38. 
    

复制代码

这是一种简单的修改非标准路径的办法，见下图，其实只是在每个Overlay 7后面单独让它又算了一遍L716：

星斗如盘

liyuanhe211 发表于 2020-5-23 13:29
可以用非标准路径，比如下列输入文件相当于在 HF/STO-3G 级别下用 opt=(TS, CalcAll, NoEigenTest) 优化 H  ...

这是自己设置计算的路径对吧？

sobereva

如置顶的新社员必读贴和论坛首页的公告栏所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并反映出帖子具体内容，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题“高斯计算使用calcall关键字优化过渡态”改了，以后务必注意

另外，keyword是关键词不是关键字

sobereva

星斗如盘 发表于 2020-5-23 13:57
这是自己设置计算的路径对吧？

仔细看
Gaussian的Link、IOp与非标准计算路径
http://sobereva.com/57

星斗如盘

sobereva 发表于 2020-5-24 01:22
仔细看
Gaussian的Link、IOp与非标准计算路径
http://sobereva.com/57

sob老师，我在使用calcall关键词优化过渡态的时候，手动的把每一帧的chk文件保存下来，然后使用freq=readfc计算每一帧的频率，发现第一帧只有一个虚频，对应我想要的振动，但是后面的帧，我想要的虚频逐渐变小，逐渐多出来另一个虚频，最终优化到了其他的结构。

calcall不是精确计算力常数吗？为什么直接就走偏了？

sobereva

星斗如盘 发表于 2020-5-24 13:39
sob老师，我在使用calcall关键词优化过渡态的时候，手动的把每一帧的chk文件保存下来，然后使用freq=read ...

用这个工具监控比你一帧帧看频率方便得多得多
http://bbs.keinsci.com/thread-14081-1-1.html

sobereva

sobereva 发表于 2020-5-25 04:06
用这个工具监控比你一帧帧看频率方便得多得多
http://bbs.keinsci.com/thread-14081-1-1.html

初始结构Hessian有且只有一个负本征值，虽然是良好的初猜，但并不保证在实际优化过渡态过程中就一定能维持住

星斗如盘

sobereva 发表于 2020-5-25 04:06
用这个工具监控比你一帧帧看频率方便得多得多
http://bbs.keinsci.com/thread-14081-1-1.html

好的，谢谢社长，我那个样子确实很麻烦

星斗如盘

sobereva 发表于 2020-5-25 04:08
初始结构Hessian有且只有一个负本征值，虽然是良好的初猜，但并不保证在实际优化过渡态过程中就一定能维 ...

嗯嗯，我也感觉这是计算方法的问题，而不是初猜的问题，和计算的精度好像也没关系，我试试别的优化方法吧。

lanthanum

在正常的含有calcall的路径里添上 IOP(7/8=1) 即可，不需要自己编辑 nonstd。