| 查看: 2097 | 回复: 8 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
用Gaussian03或09，MP2(full)/aug-cc-pvdz如何优化HI？
|
|||
|
用Gaussian优化含I的分子总是出错，请各位高人帮帮忙，看看怎么能优化出来 gjf文件 # mp2(full)/gen opt(maxcycle=400) freq scf(maxcycle=400) Title Card Required 0 1 I -0.94029852 2.25373131 0.00000000 H -2.57029852 2.25373131 0.00000000 H 0 aug-cc-pvdz **** I 0 S 6 1.00 2.449790E+03 4.190000E-04 3.598080E+02 2.240000E-03 1.440580E+01 3.972230E-01 9.076320E+00 -9.322490E-01 2.088100E+00 9.371380E-01 1.034980E+00 3.920860E-01 S 6 1.00 2.449790E+03 1.750000E-04 3.598080E+02 1.057000E-03 1.440580E+01 1.690000E-01 9.076320E+00 -4.217930E-01 2.088100E+00 6.388640E-01 1.034980E+00 3.201150E-01 S 1 1.00 3.162840E-01 1.0000000 S 1 1.00 1.217190E-01 1.0000000 S 1 1.00 4.200000E-02 1.0000000 P 5 1.00 1.953010E+01 5.893400E-02 1.108820E+01 -2.309300E-01 2.715630E+00 6.648010E-01 1.204300E+00 4.506730E-01 3.399450E-01 2.898000E-02 P 5 1.00 1.953010E+01 -1.883600E-02 1.108820E+01 8.000600E-02 2.715630E+00 -3.066520E-01 1.204300E+00 -1.475940E-01 3.399450E-01 6.075060E-01 P 1 1.00 1.108810E-01 1.0000000 P 1 1.00 3.380000E-02 1.0000000 D 5 1.00 4.547650E+01 4.266000E-03 1.319280E+01 -1.362500E-02 4.227410E+00 3.097560E-01 1.942800E+00 5.097720E-01 8.397710E-01 2.974610E-01 D 1 1.00 3.000000E-01 1.0000000 D 1 1.00 1.191000E-01 1.0000000 **** out出错信息 >>>>>>>>>> Convergence criterion not met. SCF Done: E(RHF) = -1631.30068077 A.U. after 401 cycles Convg = 0.4417D+00 -V/T = 9.3532 S**2 = 0.0000 Convergence failure -- run terminated. Error termination via Lnk1e in e:\G03W\l502.exe at Sun Jul 08 16:59:17 2012. Job cpu time: 0 days 0 hours 0 minutes 55.0 seconds. File lengths (MBytes): RWF= 11 Int= 0 D2E= 0 Chk= 4 Scr= 1 |
回复此楼» 猜你喜欢
分子动力学 Plumed教程求助
已经有0人回复
-
International Journal of Nanoscience属于中科院4区吗？
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有292人回复
-
研究生学术新锐奖｜宁波微纳材料研讨会参会通知
已经有2人回复
-
申请2027年材料/化学类博士
已经有6人回复
-
澳大利亚南昆士兰大学量子点课题组招收博士生
已经有0人回复
-
求相关合成条件
已经有12人回复
-
2027春季申博自荐
已经有6人回复
-
香港中文大学（深圳）2026年招聘公告
已经有0人回复
-
青基的材料青椒们，可以把硕士联培生送我这里，要研二男生
已经有61人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 新手请教：turbomole的基组如何选择？
已经有5人回复
- 柔性扫描两个H2分子之间的势能曲线，该怎么写gauss配置文件？
已经有7人回复
- 用MP2或者MP4进行单点能以及opt报错
已经有4人回复
- Pb2H4优化选择什么基组合适
已经有10人回复
- 分子优化总是出错，请大家帮忙
已经有9人回复
- 高斯计算分子间相互作用
已经有11人回复
- gaussian不同半经验优化得到的能量相差较大
已经有4人回复
- cc-pvdz-pp基集的使用
已经有7人回复
- 急求Gaussian 03或者09软件，
已经有4人回复
- Gaussian 03优化结构
已经有12人回复
- Gaussian.03W-for.VC
已经有200人回复
- 利用DFT-D进行Gaussian优化，输入文件应该怎么设置呢
已经有4人回复
- 电子轨道能量用什么方法计算比较好？
已经有6人回复
- this version requires full pseudpotential generation information 这是什么情况？
已经有5人回复
- 农药增效剂有哪些，天然化合物和化学合成的都行
已经有6人回复
- 全双工通信 full duplex在未来有前途么？可以作为研究方向么？
已经有4人回复
- gaussian 03和gaussian 09 优化分子的区别
已经有5人回复
- vasp优化结构中出现的问题
已经有10人回复
- 求问一般优化小分子结构用什么方法和机组
已经有10人回复
- 【讨论】用gaussian进行SCS-mp2计算
已经有12人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
1maomao: 金币+10, ★★★★★最佳答案, 解决了我的燃眉之急 2012-07-11 11:29:14
1maomao: 金币+10, ★★★★★最佳答案, 解决了我的燃眉之急 2012-07-11 11:29:14
|
%mem=1000MB %nproc=2 # mp2(full)/genecp opt(maxcycle=400) freq scf(maxcycle=400) Title Card Required 0 1 I -0.94029852 2.25373131 0.00000000 H -2.57029852 2.25373131 0.00000000 H 0 aug-cc-pVDZ **** I 0 S 9 1.00 5.546500E+03 1.560000E-04 8.382140E+02 9.860000E-04 1.821870E+02 2.792000E-03 3.121230E+01 -4.325100E-02 1.953140E+01 2.341340E-01 8.240990E+00 -7.509430E-01 2.194550E+00 8.829680E-01 1.109110E+00 4.620610E-01 8.106100E-02 1.102000E-03 S 9 1.00 5.546500E+03 -7.300000E-05 8.382140E+02 -5.080000E-04 1.821870E+02 -1.158000E-03 3.121230E+01 1.219300E-02 1.953140E+01 -8.785400E-02 8.240990E+00 3.382000E-01 2.194550E+00 -5.765500E-01 1.109110E+00 -4.092980E-01 8.106100E-02 1.432310E-01 S 1 1.00 3.746410E-01 1.0000000 S 1 1.00 1.770800E-01 1.0000000 S 1 1.00 1.001000E-01 1.0000000 S 1 1.00 4.120000E-02 1.0000000 P 9 1.00 1.889880E+02 5.850000E-04 2.128680E+01 3.692300E-02 1.003960E+01 -2.353240E-01 3.451800E+00 3.414830E-01 1.974560E+00 5.347880E-01 1.024200E+00 2.651410E-01 4.494370E-01 2.578700E-02 1.866480E-01 5.220000E-04 7.348100E-02 6.060000E-04 P 9 1.00 1.889880E+02 -2.560000E-04 2.128680E+01 -1.168200E-02 1.003960E+01 8.319200E-02 3.451800E+00 -1.569700E-01 1.974560E+00 -2.245180E-01 1.024200E+00 -1.144510E-01 4.494370E-01 3.753560E-01 1.866480E-01 5.751360E-01 7.348100E-02 2.459170E-01 P 1 1.00 5.981000E-01 1.0000000 P 1 1.00 9.618000E-02 1.0000000 P 1 1.00 3.590000E-02 1.0000000 D 7 1.00 1.326620E+02 5.720000E-04 3.760540E+01 4.402000E-03 1.038910E+01 -4.092200E-02 6.490170E+00 9.966100E-02 3.454510E+00 3.226630E-01 1.844130E+00 4.003430E-01 9.624780E-01 2.683060E-01 D 1 1.00 4.728530E-01 1.0000000 D 1 1.00 1.932000E-01 1.0000000 D 1 1.00 7.820000E-02 1.0000000 F 1 1.00 4.064000E-01 1.0000000 F 1 1.00 1.848000E-01 1.0000000 **** I 0 I-ECP 4 28 g-ul potential 1 2 1.00000000 0.00000000 s-ul potential 3 2 40.03337600 49.98964900 2 17.30057600 281.00655600 2 8.85172000 61.41673900 p-ul potential 4 2 15.72014100 67.41623900 2 15.20822200 134.80769600 2 8.29418600 14.56654800 2 7.75394900 28.96842200 d-ul potential 4 2 13.81775100 35.53875600 2 13.58780500 53.33975900 2 6.94763000 9.71646600 2 6.96009900 14.97750000 f-ul potential 4 2 18.52295000 -20.17661800 2 18.25103500 -26.08807700 2 7.55790100 -0.22043400 2 7.59740400 -0.22164600 |
7楼2012-07-10 14:52:45
Siriusgoul
铜虫 (初入文坛)
- 应助: 10 (幼儿园)
- 金币: 266
- 散金: 200
- 帖子: 48
- 在线: 27.4小时
- 虫号: 1303115
- 注册: 2011-05-22
- 性别: GG
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
gmy1990: 金币+2 2012-07-08 21:55:47
1maomao: 金币+5, ★有帮助 2012-07-09 21:52:06
1maomao: 金币+2, ★有帮助 2012-07-11 08:10:45
感谢参与，应助指数 +1
gmy1990: 金币+2 2012-07-08 21:55:47
1maomao: 金币+5, ★有帮助 2012-07-09 21:52:06
1maomao: 金币+2, ★有帮助 2012-07-11 08:10:45
|
你这个是自洽场没有收敛的问题，可以在本论坛搜一搜解决方法。 但刚看过类似的问题在如下的帖子里 http://muchong.com/bbs/viewthread.php?tid=4696851&fpage=1 先可以在scf里加xqc等关键字试试，然后是再看看对碘原子的赝势处理对不对，一般碘这样的重原子容易出问题 |
2楼2012-07-08 20:55:10
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +459 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13336.5
- 散金: 27307
- 红花: 301
- 帖子: 4683
- 在线: 928.6小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
1maomao
3楼2012-07-08 21:32:47
【答案】应助回帖
★ ★ ★
感谢参与，应助指数 +1
gmy1990: 金币+3 2012-07-09 09:47:53
感谢参与，应助指数 +1
gmy1990: 金币+3 2012-07-09 09:47:53
|
SCF不能收敛，若你输入正确则可能存在不稳定态，用ZMAT坐标或者目前的CART坐标改变一下结构先跳到“阱”内试一下，若还不行，若还不行 ，换MP2到其它更多包含关联能方法尝试一下，还不行就有可能是排斥态。 但从你的输入看，似乎有问题，I元素你用的为aug-cc-pVDZ的ECP基组，你只输入了价电子基组，对核电子的ECP没有输入，这样不会收敛的，即使你做了上述的尝试，把这段按照格式加在后边应该就没问题了，注意空格和**** I 0 I-ECP 4 28 g-ul potential 1 2 1.00000000 0.00000000 s-ul potential 3 2 40.03337600 49.98964900 2 17.30057600 281.00655600 2 8.85172000 61.41673900 p-ul potential 4 2 15.72014100 67.41623900 2 15.20822200 134.80769600 2 8.29418600 14.56654800 2 7.75394900 28.96842200 d-ul potential 4 2 13.81775100 35.53875600 2 13.58780500 53.33975900 2 6.94763000 9.71646600 2 6.96009900 14.97750000 f-ul potential 4 2 18.52295000 -20.17661800 2 18.25103500 -26.08807700 2 7.55790100 -0.22043400 2 7.59740400 -0.22164600 |
4楼2012-07-09 09:22:52
