| 查看: 2094 | 回复: 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 |
回复此楼» 猜你喜欢
Tc数据拟合
已经有1人回复
-
求相关合成条件
已经有12人回复
-
物理化学论文润色/翻译怎么收费?
已经有228人回复
-
韩国釜山大学储能电池与功能材料实验室诚招博后（电池/电化学/器件制备/MOF合成）
已经有187人回复
-
祈福！
已经有79人回复
-
2027春季申博自荐
已经有6人回复
-
散金求好运
已经有315人回复
-
博士生招生, 材料模拟方向，University of South Carolina
已经有0人回复
-
Postdoc position, materials simulations/machine learning, U. of South Carolina
已经有0人回复
-
祈福青基
已经有34人回复
-
青基的材料青椒们，可以把硕士联培生送我这里，要研二男生
已经有57人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 新手请教：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人回复
|
最后结果 Predicted change in Energy=-4.623138D-13 Optimization completed. -- Stationary point found. ---------------------------- ! Optimized Parameters ! ! (Angstroms and Degrees) ! -------------------------- -------------------------- ! Name Definition Value Derivative Info. ! -------------------------------------------------------------------------------- ! R1 R(1,2) 1.6102 -DE/DX = 0.0 ! -------------------------------------------------------------------------------- GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad 1|1|UNPC-CHJ-PC|Freq|RMP2-Full|GenECP|H1I1|CHJ|10-Jul-2012|0||#N Geom= AllCheck Guess=TCheck SCRF=Check GenChk RMP2-Full/GenECP Freq||Title C ard Required||0,1|I,-0.9501980107,2.25373131,0.|H,-2.5603990293,2.2537 3131,0.||Version=IA32W-G09RevA.02|State=1-SG|HF=-295.2438808|MP2=-295. 4938741|RMSD=3.419e-010|RMSF=2.582e-007|ZeroPoint=0.0054477|Thermal=0. 0078083|Dipole=-0.1931939,0.,0.|DipoleDeriv=-0.0256766,0.,0.,0.,-0.063 4913,0.,0.,0.,-0.0634913,0.0256766,0.,0.,0.,0.0634913,0.,0.,0.,0.06349 13|Polar=37.0356885,0.,34.7230088,0.,0.,34.7230088|PG=C*V [C*(H1I1)]|N Imag=0||0.21637266,0.,0.00000015,0.,0.,0.00000015,-0.21637266,0.,0.,0. 21637266,0.,-0.00000015,0.,0.,0.00000015,0.,0.,-0.00000015,0.,0.,0.000 00015||0.00000045,0.,0.,-0.00000045,0.,0.|||@ TO SEE A WORLD IN A GRAIN OF SAND AND A HEAVEN IN A WILD FLOWER HOLD INFINITY IN THE PALM OF YOUR HAND AND ETERNITY IN AN HOUR -- WILLIAM BLAKE Job cpu time: 0 days 0 hours 1 minutes 28.0 seconds. File lengths (MBytes): RWF= 286 Int= 0 D2E= 0 Chk= 1 Scr= 1 Normal termination of Gaussian 09 at Tue Jul 10 02:53:34 2012. |
9楼2012-07-10 14:55:02
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
