用G16广义内坐标（GIC）功能扫描两个反应坐标失败

九月九

各位老师好！我尝试扫描两个反应坐标，按照这篇文章的内容在Gaussian16中同时扫描两个反应坐标 - 腾讯云开发者社区-腾讯云 (tencent.com)，最后只算了一步单点就结束了，请问是什么回事呢？
输入文件如下：

1. %chk=methanal_water.chk
   
2. %mem=24GB
   
3. %nprocshared=24
   
4. #p M062X/6-31G(d) nosymm geom=addgic
   
5. 
   
6. Title Card Required
   
7. 
   
8. 0 1
   
9. C    -1.13256115   -0.04504504    0.00000000
   
10. H    -0.59939740   -0.97274996    0.00000000
    
11. H    -2.20256115   -0.04504504    0.00000000
    
12. O    -0.50552072    1.04600530    0.00000000
    
13. O    -0.71081810   -0.28742611    1.83667749
    
14. H    -0.42043150   -1.07235298    2.30695978
    
15. H    -0.23246558    0.54490727    1.83667749
    
16. 
    
17. R15=R(1,5)
    
18. R47=R(4,7)
    
19. R15(value=1.8)
    
20. R47(value=1.7)
    
21. 
    
22. --Link1--
    
23. %chk=methanal_water.chk
    
24. %mem=24GB
    
25. %nprocshared=24
    
26. #p opt M062X/6-31G(d) nosymm geom(allcheck,addGIC) guess=read
    
27. 
    
28. R15(freeze)
    
29. R47(freeze)
    
30. 
    
31. --Link1--
    
32. %chk=methanal_water.chk
    
33. %mem=24GB
    
34. %nprocshared=24
    
35. #p M062X/6-31G(d) nosymm geom(allcheck,addGIC) guess=read
    
36. 
    
37. R15(value=1.7)
    
38. R47(value=1.5)
    
39. 
    
40. --Link1--
    
41. %chk=methanal_water.chk
    
42. %mem=24GB
    
43. %nprocshared=24
    
44. #p opt M062X/6-31G(d) nosymm geom(allcheck,addGIC) guess=read
    
45. 
    
46. R15(freeze)
    
47. R47(freeze)
    
48. 
    
49. --Link1--
    
50. %chk=methanal_water.chk
    
51. %mem=24GB
    
52. %nprocshared=24
    
53. #p M062X/6-31G(d) nosymm geom(allcheck,addGIC) guess=read
    
54. 
    
55. R15(value=1.6)
    
56. R47(value=1.3)
    
57. 
    
58. --Link1--
    
59. %chk=methanal_water.chk
    
60. %mem=24GB
    
61. %nprocshared=24
    
62. #p opt M062X/6-31G(d) nosymm geom(allcheck,addGIC) guess=read
    
63. 
    
64. R15(freeze)
    
65. R47(freeze)
    
66. 
    
67. --Link1--
    
68. %chk=methanal_water.chk
    
69. %mem=24GB
    
70. %nprocshared=24
    
71. #p M062X/6-31G(d) nosymm geom(allcheck,addGIC) guess=read
    
72. 
    
73. R15(value=1.5)
    
74. R47(value=1.1)
    
75. 
    
76. --Link1--
    
77. %chk=methanal_water.chk
    
78. %mem=24GB
    
79. %nprocshared=24
    
80. #p opt M062X/6-31G(d) nosymm geom(allcheck,addGIC) guess=read
    
81. 
    
82. R15(freeze)
    
83. R47(freeze)

复制代码

Kalinite

参考http://sobereva.com/474
你既没有给出报错信息，抄的输入文件也莫名其妙，明明可以直接扫描，却分成了这么多任务

zjxitcc

这是过时的做法，很繁琐，不推荐使用。后来我们又发了两篇，你可能没看到
《一种用Gaussian 16中的GIC功能实现同时扫描多个坐标的方法》https://mp.weixin.qq.com/s/xrqrYCAig6HarUlLTaNMqg
《用Gaissian16中的GIC功能实现翻转过程的势能面扫描》https://mp.weixin.qq.com/s/bs-oJ6EDnFh6y1KDKqwD7A

九月九

Kalinite 发表于 2022-10-14 16:57
参考http://sobereva.com/474
你既没有给出报错信息，抄的输入文件也莫名其妙，明明可以直接扫描，却分成 ...

用的不是一般的扫描方法，一般的扫描比较耗时

九月九

zjxitcc 发表于 2022-10-14 17:10
这是过时的做法，很繁琐，不推荐使用。后来我们又发了两篇，你可能没看到
《一种用Gaussian 16中的GIC功能 ...

这方法太好了，多谢！

sobereva

同时扫两个坐标，利用GIC的灵活性令两个坐标之间满足约束关系就完了，下文都说过，参考“6.2 GIC扫描实例2：令1,3-丁二烯两个亚甲基同步旋转的扫描”
详谈使用Gaussian做势能面扫描
http://sobereva.com/474（http://bbs.keinsci.com/thread-12660-1-1.html）

九月九

sobereva 发表于 2022-10-14 20:30
同时扫两个坐标，利用GIC的灵活性令两个坐标之间满足约束关系就完了，下文都说过，参考“6.2 GIC扫描实例2 ...

已学习，谢谢老师！