计算化学公社

标题: gaussian固定坐标 [打印本页]
作者
Author: 诸葛壹次心    时间: 2015-2-13 23:58
标题: gaussian固定坐标
我想要在优化的时候用z-matrix的方法固定坐标，输入文件如下，希望固定键长B14，但是优化出来B14被改变了，我记得之前使用ccsd(t)的时候是成功了，但是为什么换成hf就不行了呢？

  1. #p uhf/3-21g opt freq

  3. test of uhf

  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
复制代码


作者
Author: 诸葛壹次心    时间: 2015-2-13 23:59
gaussian如何固定键长
作者
Author: sobereva    时间: 2015-2-14 00:16
由于是内坐标方式固定的，要写成opt=z-matrix。否则会在冗余内坐标下优化，内坐标的冻结设定不起作用。
作者
Author: 诸葛壹次心    时间: 2015-2-14 10:25
sobereva 发表于 2015-2-14 00:16
由于是内坐标方式固定的，要写成opt=z-matrix。否则会在冗余内坐标下优化，内坐标的冻结设定不起作用。

好的，谢谢sob！

PS: 可是好诡异啊，我之前用CCSD(t)的时候opt后面什么也没加居然也可以有效果。。。
作者
Author: sobereva    时间: 2015-2-14 11:26
诸葛壹次心 发表于 2015-2-14 10:25
好的，谢谢sob！

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

CCSD(T)这种没有解析梯度的方法直接那么写就是基于内坐标优化的
作者
Author: 诸葛壹次心    时间: 2015-2-15 16:40
sobereva 发表于 2015-2-14 11:26
CCSD(T)这种没有解析梯度的方法直接那么写就是基于内坐标优化的

原来如此！谢谢sob~
作者
Author: xinren    时间: 2018-12-16 23:28
优化激发态的结构，分子在激发态比较容易进行TICT,优化不出来LE的结构。试着尝试了固定二面角的方式得到了一个相近的结构。请问有什么比较好的方法优化出这个结构，有没有扫描一下势能面的必要。也有调小过步长优化。请问老师一般这种情况怎么处理@sobereva 。另外有一个分子会出现能量震荡，有按着论坛的方法尝试，还是震荡，体系主要是偶氮苯的衍生物。谢谢老师
作者
Author: sobereva    时间: 2018-12-17 00:27
xinren 发表于 2018-12-16 23:28
优化激发态的结构，分子在激发态比较容易进行TICT,优化不出来LE的结构。试着尝试了固定二面角的方式得到了 ...

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



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