计算化学公社

标题: g09wD.01如何用自定义非标准路径实现opt=recalc [打印本页]
作者
Author: ldatea    时间: 2020-4-23 22:25
标题: g09wD.01如何用自定义非标准路径实现opt=recalc
本帖最后由 ldatea 于 2020-6-28 12:05 编辑

请问Windows版G09 D.01能用自定义非标准路径实现opt=recalc吗？
如果能实现的话，有谁能给一个G16几何优化使用reclac的路径。  谢谢


我用testret.exe看了一下#p opt=calcfc 6-311g(d,p) nosymm geom=connectivity m062x em=gd3
  1. 1/10=4,14=-1,18=20,19=15,26=3,38=1,57=2/1,3;
  2. 2/9=110,12=2,15=1,17=6,18=5,40=1/2;
  3. 3/5=4,6=6,7=101,11=2,16=1,25=1,30=1,71=2,74=-55,124=31,140=1/1,2,3;
  4. 4//1;
  5. 5/5=2,38=5/2;
  6. 8/6=4,10=90,11=11/1;
  7. 11/6=1,8=1,9=11,15=111,16=1,31=1/1,2,10;
  8. 10/6=1,13=1,31=1/2;
  9. 6/7=2,8=2,9=2,10=2,28=1/1;
  10. 7/10=1,18=20,25=1,30=1/1,2,3,16;
  11. 1/10=4,14=-1,18=20,19=15,26=3/3(2);
  12. 2/9=110,15=1/2;
  13. 99//99;
  14. 2/9=110,15=1/2;
  15. 3/5=4,6=6,7=101,11=2,16=1,25=1,30=1,71=1,74=-55,124=31/1,2,3;
  16. 4/5=5,16=3,69=1/1;
  17. 5/5=2,38=5/2;
  18. 7/30=1/1,2,3,16;
  19. 1/14=-1,18=20,19=15,26=3/3(-5);
  20. 2/9=110,15=1/2;
  21. 6/7=2,8=2,9=2,10=2,19=2,28=1/1;
  22. 99/9=1/99;
复制代码
#p opt=calcall 6-311g(d,p) nosymm geom=connectivity m062x em=gd3
  1. 1/10=4,14=-1,18=20,19=15,26=3,38=1,57=2/1,3;      ！和calcfc第 1 行一致
  2. 2/9=110,12=2,15=1,17=6,18=5,40=1/2;                      ！一致
  3. 3/5=4,6=6,7=101,11=2,16=1,25=1,30=1,71=2,74=-55,124=31,140=1/1,2,3;     ！一致
  4. 4//1;         ！一致
  5. 5/5=2,38=5/2;       ！一致
  6. 8/6=4,10=90,11=11/1;        ！一致
  7. 11/6=1,8=1,9=11,15=111,16=1,31=1/1,2,10;                 ！一致
  8. 10/6=1,31=1/2;                                    ！与同样位置的calcfc相比少了13=1
  9. 6/7=2,8=2,9=2,10=2,18=1,28=1/1;        ！与同样位置的calcfc相比多了18=1
  10. 7/10=1,25=1,30=1/1,2,3,16;                 ！与同样位置的calcfc相比少了18=20
  11. 1/10=4,14=-1,18=20,19=15,26=3/3(3);   ！一致
  12. 2/9=110,15=1/2;       ！一致
  13. 7/8=1,9=1,25=1,30=1,44=-1/16;       ！比calcfc多出来的一行，大概是算频率用的
  14. 99//99;      ！一致
  15. 2/9=110,15=1/2;    ！一致
  16. 3/5=4,6=6,7=101,11=2,16=1,25=1,30=1,71=2,74=-55,124=31,140=1/1,2,3; ！和  callall  和  calcfc  第 3 行一致，calcfc 倒数第 8 行相比，这里是 71=1（calcfc是71=2），多了140=1
  17. 4/5=5,16=3,69=1/1;  ！一致
  18. 5/5=2,38=5/2;         ！一致
  19. 8/6=4,10=90,11=11/1;    ！calcfc没有这行
  20. 11/6=1,8=1,9=11,15=111,16=1,31=1/1,2,10;     ！calcfc没有这行
  21. 10/6=1,31=1/2;     ！calcfc没有这行
  22. 7/10=1,25=1,30=1/1,2,3,16;          ! 和calcfc倒数第 5 行相比，多了10=1,25=1,30=1
  23. 1/10=4,14=-1,18=20,19=15,26=3/3(-8);          ! 和calcfc倒数第 4 行相比，多了10=4
  24. 2/9=110,15=1/2;         ！和calcfc倒数第 3 行一致
  25. 6/7=2,8=2,9=2,10=2,18=1,19=2,28=1/1;      ！和calcfc倒数第 2 行一致
  26. 7/8=1,9=1,25=1,30=1,44=-1/16;                 ！比calcfc多出来的一行，大概是算频率用的
  27. 99//99;         ！和calcfc的倒数第1行相比，少了9=1
复制代码
sob老师给过.help文件，但是这些不太容易看懂，目前在尝试看懂。
作者
Author: sobereva    时间: 2020-4-23 23:57
显然能，但会牵扯到好多地方，自己写起来费劲。所以建议拿一个标准的带recalc的自动产生的非标准路径信息作为模板，在里面修改对应计算级别之类的IOp
作者
Author: ldatea    时间: 2020-4-24 00:34
sobereva 发表于 2020-4-23 23:57
显然能，但会牵扯到好多地方，自己写起来费劲。所以建议拿一个标准的带recalc的自动产生的非标准路径信息作 ...

谢谢sob老师。我再去研究一下，这个主要还是学习过程，用起来不一定效果好。




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3