| 查看: 1765 | 回复: 8 | ||||
yaoqian0503木虫 (正式写手)
|
[求助]
molpro初学者，遇到问题实在不晓得怎么解决，求大神指点
|
|
初学者，遇到问题实在不晓得怎么解决，求大神指点 输入文件 ***,c4h10 memory,1800,m B1=1.09615979 ang B2=1.09615979 ang B3=1.09504935 ang B4=1.53155062 ang B5=1.09878447 ang B6=1.09878447 ang B7=1.53000083 ang B8=1.09878447 ang B9=1.09878447 ang B10=1.53155062 ang B11=1.09615979 ang B12=1.09504935 ang B13=1.09615979 ang A1=107.44288903 degree A2=107.65250504 degree A3=111.57441291 degree A4=109.51359724 degree A5=109.51359724 degree A6=113.32856139 degree A7=109.14466006 degree A8=109.14466006 degree A9=113.32856139 degree A10=111.16149499 degree A11=111.57441291 degree A12=111.16149499 degree D1=115.69334126 degree D2=-122.22242661 degree D3=-57.87636554 degree D4=57.87636554 degree D5=180.00000000 degree D6=57.67107344 degree D7=-57.67107344 degree D8=180.00000000 degree D9=59.81905623 degree D10=180.00000000 degree D11=-59.81905623 degree geometry={ C1; H1,C1,B1, H2,C1,B2,H1,A1, H3,C1,B3,H1,A2,H2,D1 C2,C1,B4,H3,A3,H1,D2; H4,C2,B5,C1,A4,H3,D3; H5,C2,B6,C1,A5,H3,D4; C3,C2,B7,C1,A6,H3,D5; H6,C3,B8,C2,A7,C1,D6; H7,C3,B9,C2,A8,C1,D7; C4,C3,B10,C2,A9,C1,D8; H8,C4,B11,C3,A10,C2,D9; H9,C4,B12,C3,A11,C2,D10; H10,C4,B13,C3,A12,C2,D11; } basis=AVDZ {hf;wf,34,1,0;} {casscf;occ,18;closed,16;} {optg,savexyz=c4h10-1.xyz inactive,B7} frequencies ecas=energy {rs2c;occ,18;closed,16;orbital,IGNORE_ERROR;maxit,500} emp2=energy table,B1,ecas,emp2 错误信息 对称数选为1的时候 报错信息 The request for symmetry elements XY, cannot be honoured ? Error ? Symmetry problem ? The problem occurs in zmatrix.f:zmat_orient 对称数选为2的时候，报错信息 SYYMETRY= 2 MS2= 0 ? Error ? Open-shell singlet not implemented ? The problem occurs in scfpro |
回复此楼» 猜你喜欢
江西理工大学南昌校区国际创新研究院无序物质科学研究中心李建荣课题组招调剂
已经有0人回复
-
一志愿中科大化材院物理化学专业，求调剂
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有154人回复
-
纳米生物医用材料招研究生
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【原创】盘点一些的理论化学书籍与给初学者的建议
已经有36人回复
- Multiwfn支持的弱相互作用的分析方法概览
已经有10人回复
- 大体系弱相互作用计算的解决之道
已经有42人回复
- Multiwfn 3.0波函数分析程序的意义、功能与用途
已经有38人回复
- 【拉登】RCS 期刊资源搜集 （请务必不要水帖）
已经有6人回复
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +453 - QC强帖: 12
- 应助: 85 (初中生)
- 金币: 10750.7
- 散金: 27307
- 红花: 295
- 帖子: 4586
- 在线: 916.7小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
yaoqian0503(zhou2009代发): 金币+10 2016-06-06 20:42:40
yaoqian0503(zhou2009代发): 金币+10 2016-06-06 20:42:40
|
变量设置有问题，不符合对称性要求：那些初值相同的参数，在优化过程中到底想等还是不等，程序不知道。 要么就降低对称性再算，要么就减少坐标变量 |
» 本帖已获得的红花（最新10朵）
2楼2016-06-05 22:17:47
lastzealot
新虫 (著名写手)
【乱世的奸雄】
- 应助: 87 (初中生)
- 金币: 2668.1
- 散金: 876
- 红花: 24
- 帖子: 2914
- 在线: 212.3小时
- 虫号: 730090
- 注册: 2009-03-24
- 性别: GG
- 专业: 高分子材料结构与性能
3楼2016-06-05 22:26:28
yaoqian0503
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 4351.2
- 红花: 8
- 帖子: 479
- 在线: 91.6小时
- 虫号: 1462930
- 注册: 2011-10-27
- 性别: MM
- 专业: 理论和计算化学
送红花一朵|
减少变量的那个输入文件我已经做过了，我的结构是下面这种的，但是还是一样的报错，减少对称性那个我也做了，就是改了些键长键角参数，之前也是看到你以前回的帖子，不晓得是不是我哪里做错了 geometry={angstrom; C1; H1,C1,R1, H2,C1,R1,H1,A1, H3,C1,R2,H1,A2,H2,D1; C2,C1,R3,H3,A3,H1,D2; H4,C2,R4,C1,A4,H3,D3; H5,C2,R4,C1,A4,H3,-D3; C3,C2,1.53,C1,A5,H3,D4; H6,C3,R4,C2,A6,C1,D5; H7,C3,R4,C2,A6,C1,-D5; C4,C3,R3,C2,A5,C1,D4; H8,C4,R1,C3,A7,C2,D6; H9,C4,R2,C3,A3,C2,D4; H10,C4,R1,C3,A7,C2,-D6; } R1=1.09615979 ang R2=1.09504935 ang R3=1.53155062 ang R4=1.09878447 ang A1=107.44288903 degree A2=107.65250504 degree A3=111.57441291 degree A4=109.51359724 degree A5=113.32856139 degree A6=109.14466006 degree A7=111.16149499 degree D1=115.69334126 degree D2=-122.22242661 degree D3=-57.87636554 degree D4=180.00000000 degree D5=57.67107344 degree D6=59.81905623 degree |
4楼2016-06-06 08:20:32
yaoqian0503
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 4351.2
- 红花: 8
- 帖子: 479
- 在线: 91.6小时
- 虫号: 1462930
- 注册: 2011-10-27
- 性别: MM
- 专业: 理论和计算化学
5楼2016-06-06 08:22:35
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +453 - QC强帖: 12
- 应助: 85 (初中生)
- 金币: 10750.7
- 散金: 27307
- 红花: 295
- 帖子: 4586
- 在线: 916.7小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
zhou2009: 金币+10 2016-06-06 20:43:08
zhou2009: 金币+10 2016-06-06 20:43:08
|
用同样的坐标，gaussian下做一个单点计算，然后找到“Deg. of freedom”，这就是独立变量的个数。 如果molpro优化的变量个数（包括没有特殊意义的常数1.53在内）大于此数，坐标肯定是错的 |
» 本帖已获得的红花（最新10朵）
6楼2016-06-06 08:38:28
yaoqian0503
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 4351.2
- 红花: 8
- 帖子: 479
- 在线: 91.6小时
- 虫号: 1462930
- 注册: 2011-10-27
- 性别: MM
- 专业: 理论和计算化学
7楼2016-06-06 08:49:30
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +453 - QC强帖: 12
- 应助: 85 (初中生)
- 金币: 10750.7
- 散金: 27307
- 红花: 295
- 帖子: 4586
- 在线: 916.7小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
yaoqian0503(zhou2009代发): 金币+10 2016-06-06 20:43:38
yaoqian0503(zhou2009代发): 金币+10 2016-06-06 20:43:38
|
先把HC4H骨架建立好，检查坐标个数是否符合当前的“Deg. of freedom”。 然后再加8个H，重新检查。 |
» 本帖已获得的红花（最新10朵）
8楼2016-06-06 10:29:40
yaoqian0503
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 4351.2
- 红花: 8
- 帖子: 479
- 在线: 91.6小时
- 虫号: 1462930
- 注册: 2011-10-27
- 性别: MM
- 专业: 理论和计算化学
9楼2016-06-06 10:59:54
10