gaussian固定坐标

诸葛壹次心

我想要在优化的时候用z-matrix的方法固定坐标，输入文件如下，希望固定键长B14，但是优化出来B14被改变了，我记得之前使用ccsd(t)的时候是成功了，但是为什么换成hf就不行了呢？

1. #p uhf/3-21g opt freq
   
2. 
   
3. test of uhf
   
4. 
   
5. 0 2
   
6. C
   
7. H,1,B1
   
8. H,1,B2,2,A1
   
9. H,1,B3,3,A2,2,D1 0
   
10. C,1,B4,3,A3,2,D2 0
    
11. H,5,B5,1,A4,3,D3 0
    
12. C,5,B6,1,A5,3,D4 0
    
13. H,7,B7,5,A6,1,D5 0
    
14. H,7,B8,5,A7,1,D6 0
    
15. C,5,B9,1,A8,7,D7 0
    
16. H,10,B10,5,A9,1,D8 0
    
17. H,10,B11,5,A10,1,D9 0
    
18. O,7,B12,5,A11,1,D10 0
    
19. H,13,B13,7,A12,5,D11 0
    
20. O,10,B14,5,A13,1,D12 0
    
21. O,15,B15,10,A14,5,D13 0
    
22. Variables:
    
23. B1=1.09331724
    
24. B2=1.09103427
    
25. B3=1.09270253
    
26. B4=1.53980793
    
27. B5=1.11826342
    
28. B6=1.52925655
    
29. B7=1.09913708
    
30. B8=1.09856152
    
31. B9=1.47120666
    
32. B10=1.08273318
    
33. B11=1.08216781
    
34. B12=1.42672228
    
35. B13=0.96129605
    
36. B15=1.23755496
    
37. A1=108.21528842
    
38. A2=108.77218725
    
39. A3=110.66888452
    
40. A4=107.21925057
    
41. A5=112.66540001
    
42. A6=109.24334854
    
43. A7=109.22114223
    
44. A8=113.3742421
    
45. A9=120.59855293
    
46. A10=121.04870225
    
47. A11=108.2641961
    
48. A12=108.34602049
    
49. A13=97.91817973
    
50. A14=108.49530908
    
51. D1=-116.86692823
    
52. D2=121.17143969
    
53. D3=-176.06339647
    
54. D4=-60.50298586
    
55. D5=-174.73599909
    
56. D6=-56.82232519
    
57. D7=128.87846502
    
58. D8=22.40409801
    
59. D9=-142.47354454
    
60. D10=64.28766706
    
61. D11=-176.65152912
    
62. D12=119.65744006
    
63. D13=-1.35903404
    
64. Constants:
    
65. B14=2.3
    

复制代码

sobereva

xinren 发表于 2018-12-16 23:28
优化激发态的结构，分子在激发态比较容易进行TICT,优化不出来LE的结构。试着尝试了固定二面角的方式得到了 ...

和主贴没直接关系的问题请单独发帖，或者在简单量化答疑转贴里跟帖提问

xinren

优化激发态的结构，分子在激发态比较容易进行TICT,优化不出来LE的结构。试着尝试了固定二面角的方式得到了一个相近的结构。请问有什么比较好的方法优化出这个结构，有没有扫描一下势能面的必要。也有调小过步长优化。请问老师一般这种情况怎么处理@sobereva 。另外有一个分子会出现能量震荡，有按着论坛的方法尝试，还是震荡，体系主要是偶氮苯的衍生物。谢谢老师

诸葛壹次心

sobereva 发表于 2015-2-14 11:26
CCSD(T)这种没有解析梯度的方法直接那么写就是基于内坐标优化的

原来如此！谢谢sob~

sobereva

诸葛壹次心 发表于 2015-2-14 10:25
好的，谢谢sob！

PS: 可是好诡异啊，我之前用CCSD(t)的时候opt后面什么也没加居然也可以有效果。。。

CCSD(T)这种没有解析梯度的方法直接那么写就是基于内坐标优化的

诸葛壹次心

sobereva 发表于 2015-2-14 00:16
由于是内坐标方式固定的，要写成opt=z-matrix。否则会在冗余内坐标下优化，内坐标的冻结设定不起作用。

好的，谢谢sob！

PS: 可是好诡异啊，我之前用CCSD(t)的时候opt后面什么也没加居然也可以有效果。。。

sobereva

由于是内坐标方式固定的，要写成opt=z-matrix。否则会在冗余内坐标下优化，内坐标的冻结设定不起作用。

诸葛壹次心

gaussian如何固定键长