| 查看: 962 | 回复: 7 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】 优化结构后计算红外的问题
|
|||
|
我把优化后结构用GV导出后计算频率，在最后的结果中怎么会出现这句话？高手帮忙解释下，我刚开始学习，不太懂，谢谢！ Item Value Threshold Converged? Maximum Force 0.000541 0.000450 NO RMS Force 0.000088 0.000300 YES Maximum Displacement 0.197210 0.001800 NO RMS Displacement 0.043308 0.001200 NO Predicted change in Energy=-7.253577D-05 GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad [ Last edited by lei0736 on 2009-3-29 at 09:15 ] |
回复此楼» 猜你喜欢
085404（计算机技术）293分 希望调剂到材料相关专业
已经有0人回复
-
ELISA试验中不容忽视的细节盘点(二)
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有57人回复
-
细胞培养，这22个细节一定要注意！（一）
已经有0人回复
-
已知面上挂了，小学校实验方向无平台，感觉太难了
已经有26人回复
-
CSC访学博后项目获批，外方学校暂停合作该怎么办？
已经有51人回复
-
求助各位大佬，球球了，在这科萨头
已经有0人回复
-
分子动力学模拟合作需求
已经有3人回复
-
招生材料与化工、环境科学调剂研究生
已经有5人回复
iamddba
铁杆木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 7655.2
- 帖子: 343
- 在线: 111.1小时
- 虫号: 601558
- 注册: 2008-09-13
- 性别: GG
- 专业: 理论和计算化学
2楼2009-03-28 16:00:47
3楼2009-03-28 16:03:45
nyzhaoyin
新虫 (正式写手)
- 应助: 27 (小学生)
- 金币: 1044.9
- 散金: 112
- 红花: 9
- 帖子: 681
- 在线: 241.5小时
- 虫号: 563016
- 注册: 2008-05-24
- 性别: GG
4楼2009-03-28 17:30:33
|
我没有动过图像是直接导出的，而且导出的跟最后优化的不一样啊， 导出的文件： B1 5.84968305 B2 1.68199593 B3 1.84555824 B4 1.71934003 B5 1.72004643 B6 1.93328494 B7 1.70705285 B8 1.68199587 B9 1.48156639 B10 1.48091548 B11 1.48206373 B12 1.48091770 B13 1.48156085 B14 1.48206557 B15 0.95727458 B16 0.95594492 A1 5.41312099 A2 134.30399312 A3 110.73677025 A4 115.13606586 A5 128.08156939 A6 139.68553755 A7 5.41083569 A8 107.39383297 A9 109.68774667 A10 108.94025217 A11 104.32809223 A12 109.35295796 A13 112.24966032 A14 121.39099509 A15 126.24013892 D1 -105.33358002 D2 62.49084447 D3 -71.49841179 D4 169.19022883 D5 7.39914168 D6 0.00905043 D7 6.72768525 D8 -113.57763628 D9 126.23221987 D10 -171.94574084 D11 69.66113059 D12 -53.13419615 D13 -133.54395050 D14 -48.37114041 而优化的结果是： Variables: B1=5.84968288 B2=1.68199543 B3=1.84555929 B4=1.71934004 B5=1.72004623 B6=1.93328534 B7=1.70705296 B8=1.68199541 B9=1.48156589 B10=1.48091543 B11=1.48206339 B12=1.48091803 B13=1.48156102 B14=1.48206581 B15=0.95727505 B16=0.95594455 A1=5.41311053 A2=134.30399324 A3=110.73671073 A4=115.13605892 A5=128.08160563 A6=139.68550775 A7=5.41082467 A8=107.393842 A9=109.68772791 A10=108.94024537 A11=104.32810347 A12=109.35294536 A13=112.24966951 A14=121.39097284 A15=126.24009325 D1=-105.3335538 D2=62.49087201 D3=-71.49835638 D4=169.19027961 D5=7.3991061 D6=0.00907922 D7=175.9179047 D8=55.61260817 D9=-64.57751007 D10=-171.94569022 D11=69.66117654 D12=-53.13412297 D13=-133.54399572 D14=-48.37117686 为什么会这样呢？？？ |
5楼2009-03-28 18:00:57
6楼2009-03-31 00:01:47
alwens
铁杆木虫 (正式写手)
老木虫
- 应助: 0 (幼儿园)
- 贵宾: 0.45
- 金币: 5208.8
- 散金: 50
- 帖子: 994
- 在线: 806.8小时
- 虫号: 8486
- 注册: 2003-04-27
- 性别: GG
- 专业: 药物设计与药物信息
7楼2009-03-31 00:08:41
Jasminer
铁杆木虫 (著名写手)
- 应助: 52 (初中生)
- 金币: 7407.8
- 散金: 1845
- 红花: 27
- 帖子: 1002
- 在线: 728.8小时
- 虫号: 620156
- 注册: 2008-10-08
- 性别: GG
- 专业: 应用有机化学
★ ★ ★ ★
longyun80(金币+1):谢谢参与
yuhuobuku(金币+3):欢迎参加讨论
longyun80(金币+1):谢谢参与
yuhuobuku(金币+3):欢迎参加讨论
|
一方面，用gv导出时，由于对有效数字的截断，结构有细微的变化，这个比较优化和频率计算时算出的能量就知道，后几位是不一样的，解决方法是用geom=checkpoint从chk文件读入结构。 另一方面，优化时对力常数采用的是估算的方法，在一些情况下与用频率计算得到的力常数有较明显的差异。尽管如此，大多数情况下仍可以完成优化。但这时就会发现，原本4个yes，频率计算完再评价收敛结果，yes个数变少了。通常如果前两项标准还是yes的话，我们认为优化结构仍是可以接受的。但像你这样差别很大的，应重新优化，可以将所得结构用geom=checkpoint继续opt，同时加readfc参数从频率计算的chk中读取力常数。 [ Last edited by Jasminer on 2009-3-31 at 13:38 ] |
8楼2009-03-31 13:37:17
