| 查看: 810 | 回复: 13 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】新手求助高斯几何优化的问题
|
|||
|
本人是新手，对磷酸镧进行几何优化，建了很多结构都不成功，主要是不收敛，能量没有出现震荡，希望高手给予指点．输入文件如下： %chk=24conventional cell.chk %mem=300MW %nproc=1 # B3LYP/LanL2DZ opt test Title Card Required ０１ La 0.00000000 1.49648495 2.12686250 La 0.00000000 6.46932505 6.38058750 La 3.93761000 5.47938995 2.12686250 La 3.93761000 2.48642005 6.38058750 P 5.28300573 0.00000000 4.25372500 P 5.28300573 7.96581000 4.25372500 P 6.52982427 3.98290500 8.50745000 P 6.52982427 3.98290500 0.00000000 O 0.47592026 2.75509384 4.03856895 O 0.47592026 5.21071616 4.46888105 O 2.21588310 4.18761771 3.01245189 O 2.21588310 3.77819229 5.49499811 O 5.65933690 4.18761771 1.24127311 O 5.65933690 3.77819229 7.26617689 O 3.46168974 1.22781116 8.29229395 O 3.46168974 6.73799884 0.21515605 O 1.72172690 0.20471271 1.24127311 O 1.72172690 7.76109729 7.26617689 O 4.41353026 1.22781116 4.46888105 O 4.41353026 6.73799884 4.03856895 O 7.39929974 5.21071616 8.29229395 O 7.39929974 2.75509384 0.21515605 O 6.15349310 7.76109729 5.49499811 O 6.15349310 0.20471271 3.01245189 |
回复此楼» 猜你喜欢
中国科学院成都山地灾害与环境研究所2026年特别研究助理（博士后）招聘启事
已经有0人回复
-
全奖博士 英国利物浦大学 × 台湾清华大学 联合培养
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有287人回复
-
寻博导
已经有13人回复
-
基金正文30页指的是报告正文还是整个申请书
已经有9人回复
-
哈尔滨工业大学（深圳）绿色化学团队招收2026年秋季入学博士生
已经有11人回复
-
澳大利亚昆士兰科技大学PhD博士奖学金
已经有0人回复
S07111072
铁杆木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 8405.7
- 散金: 9
- 红花: 3
- 沙发: 2
- 帖子: 2048
- 在线: 345.5小时
- 虫号: 491192
- 注册: 2008-01-05
- 性别: MM
- 专业: 理论和计算化学
2楼2009-06-17 16:20:39
hairan
木虫 (著名写手)
瞌睡木虫
- QC强帖: 3
- 应助: 61 (初中生)
- 贵宾: 1.06
- 金币: 4984.3
- 散金: 8072
- 红花: 44
- 帖子: 2900
- 在线: 733.1小时
- 虫号: 432069
- 注册: 2007-08-12
- 性别: GG
- 专业: 理论和计算化学
★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yuhuobuku(金币+4,VIP+0):谢谢交流 6-17 18:03
小木虫(金币+0.5):给个红包，谢谢回帖交流
yuhuobuku(金币+4,VIP+0):谢谢交流 6-17 18:03
|
首先你给出的输入文件格式有问题，不知道你在计算中是否就是用这样的文件 其次你搭的结构也有问题，哪有那么多氧原子堆在一块儿的 另外你给的信息太少了， 计算出错了吗？ 如果出错了，贴出输出文件最后的出错信息 如果没有出错，需要贴出不收敛的信息， 在你的输出文件中查找“Threshhold” 字段，找到最后一个 将其下面的四行贴出来， 顺便看看最后几个“Threshhold” 字段下面的情况 看你的chk文件的名字，好像还用了cell |
3楼2009-06-17 16:59:39
erylingjet
铁杆木虫 (著名写手)
- 应助: 0 (幼儿园)
- 贵宾: 5.372
- 金币: 9933.4
- 散金: 752
- 红花: 2
- 帖子: 2248
- 在线: 241.5小时
- 虫号: 415302
- 注册: 2007-06-28
- 性别: GG
- 专业: 计算化学
4楼2009-06-17 21:02:54
orangedp6597
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 32.9
- 散金: 1
- 帖子: 157
- 在线: 122.9小时
- 虫号: 777169
- 注册: 2009-05-22
- 专业: 理论和计算化学
5楼2009-06-18 01:49:35
yjcmwgk
禁虫 (文坛精英)
密度泛函·小卒
- QC强帖: 38
- 应助: 165 (高中生)
- 贵宾: 16.999
- 金币: 554.6
- 散金: 39880
- 红花: 395
- 沙发: 2
- 帖子: 10406
- 在线: 1528.1小时
- 虫号: 448368
- 注册: 2007-11-01
- 专业: 理论和计算化学
- 管辖: 科研工具资源
6楼2009-06-18 13:07:04
yjcmwgk
禁虫 (文坛精英)
密度泛函·小卒
- QC强帖: 38
- 应助: 165 (高中生)
- 贵宾: 16.999
- 金币: 554.6
- 散金: 39880
- 红花: 395
- 沙发: 2
- 帖子: 10406
- 在线: 1528.1小时
- 虫号: 448368
- 注册: 2007-11-01
- 专业: 理论和计算化学
- 管辖: 科研工具资源
★ ★ ★
yuhuobuku(金币+3,VIP+0):欢迎交流 6-18 21:26
yuhuobuku(金币+3,VIP+0):欢迎交流 6-18 21:26
|
Guess=indo的意思是：使用Gaussian98的默认初始猜测：对第一行元素是INDO，第二行元素是CNDO，第三行及以后的元素是Huckel。 如果不写Guess=indo，则使用Gaussian03的默认初始猜测：全部使用Huckel初始猜测。 Guess=indo的用处：使比较差的初始构型，能比较顺利地通过初始猜测。 想不到楼主的输入文件写的一塌糊涂，却居然知道Guess=indo！ [ Last edited by yjcmwgk on 2009-6-18 at 13:11 ] |
7楼2009-06-18 13:10:36
yjcmwgk
禁虫 (文坛精英)
密度泛函·小卒
- QC强帖: 38
- 应助: 165 (高中生)
- 贵宾: 16.999
- 金币: 554.6
- 散金: 39880
- 红花: 395
- 沙发: 2
- 帖子: 10406
- 在线: 1528.1小时
- 虫号: 448368
- 注册: 2007-11-01
- 专业: 理论和计算化学
- 管辖: 科研工具资源
8楼2009-06-18 13:30:21
9楼2009-06-18 13:59:49
|
非常感谢大家的回帖,我也知道问题,所以才拿出来让大家帮着看看,太感谢你们了 下面是最后三个Threshold的值,还请大家帮忙分析: Item Value Threshold Converged? Maximum Force 0.093072 0.000450 NO RMS Force 0.022360 0.000300 NO Maximum Displacement 0.074187 0.001800 NO RMS Displacement 0.016424 0.001200 NO Predicted change in Energy=-5.518068D-02 Item Value Threshold Converged? Maximum Force 0.082378 0.000450 NO RMS Force 0.020155 0.000300 NO Maximum Displacement 0.010794 0.001800 NO RMS Displacement 0.002636 0.001200 NO Predicted change in Energy=-2.249731D-05 Item Value Threshold Converged? Maximum Force 0.084392 0.000450 NO RMS Force 0.020264 0.000300 NO Maximum Displacement 0.199430 0.001800 NO RMS Displacement 0.062802 0.001200 NO Predicted change in Energy=-3.685631D-03 [ Last edited by huabangou on 2009-6-18 at 15:42 ] |
10楼2009-06-18 15:19:30
