24小时热门版块排行榜

返回列表

迷茫的小飞侠

金虫 (正式写手)

应助: 5 (幼儿园)
金币: 758.4
帖子: 359
在线: 157.9小时
虫号: 1375147
注册: 2011-08-21
性别: MM
专业: 材料物理化学

[求助] 高斯l502.exe报错和关于中多种基组的命名

1.现在在做硅晶胞晶体中掺杂稀土金属的测试（主要是f区金属元素），在优化中一直出现l502.exe的报错而且并没有出现正常的优化结构（收敛部分基本是没有出现的），为此使用了各种解决方案，依旧没有解决.
使用的关键词是p B3lyp/Genecp NoSym Int=UltraFine SCF(dsymm,MaxCyc=400) punch(Coord) Opt(Cartesian,gdiis)，基组是根据近期的文献使用了ECP28MWB for La，6-31G* for Si H
下面内容是出现的502报错
Cycle 400  Pass 1  IDiag  1:
RMSU=  3.13D-09 CP:  1.00D+00 -2.61D+00 -2.82D+00 -9.22D-01 -2.97D+00
                  CP: -8.62D-02
E= -11109.0540824687    Delta-E=       0.000000003791 Rises=F Damp=F
DIIS: error= 1.56D-06 at cycle 400 NSaved= 7.
NSaved= 7 IEnMin= 2 EnMin= -11109.0540824750    IErMin= 7 ErrMin= 1.56D-06
ErrMax= 1.56D-06  0.00D+00 EMaxC= 1.00D+00 BMatC= 1.20D-10 BMatP= 1.20D-10
IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00
Large coefficients: NSaved=  7 BigCof= 29.34 CofMax= 10.00 Det=-2.94D-16
Inversion failed.  Reducing to  6 matrices.
Coeff-Com:  0.140D+01 0.324D+00-0.123D+01 0.499D+00-0.129D-01 0.178D-01
Coeff:    0.140D+01 0.324D+00-0.123D+01 0.499D+00-0.129D-01 0.178D-01
Gap=    0.017 Goal= None Shift= 0.000
RMSDP=1.53D-07 MaxDP=3.50D-05 DE= 3.79D-09 OVMax= 8.29D-05

>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(RB3LYP) =  -11109.0540825    A.U. after  401 cycles
         NFock=400  Conv=0.15D-06    -V/T= 2.0417
KE= 1.066459942440D+04 PE=-5.523222530108D+04 EE= 1.908083676521D+04
Convergence failure -- run terminated.
Error termination via Lnk1e in /home/programmes/gaussian/g09.D.01/g09/l502.exe at Sat Mar 12 18:25:01 2016.
Job cpu time:    0 days 23 hours 47 minutes 57.9 seconds.
File lengths (MBytes):  RWF= 232 Int=    0 D2E=    0 Chk=    29 Scr=    1

由于相关的解决措施已经尝试，现在依旧不能得到优化结果，渴求各位前辈指导！

2. 给予以上优化我想采用混合基组（两种基组的混合定义是会的），我能否对La和Si同时进行赝势基组的定义？因为适合f区元素的基组较少且文献中采用stuttgart RSC 1997 ECP，而Si则可以使用常规的SDD进行定义. 在测试中是采用同时对La和Si都适用的ECP基组（CRENBL ECP）还是分别定义两种赝势基组：ECP28MWB for La+SDD for Si. 如果是分别定义基组，在输入结构中就需要对三种混合基组进行定义（包括用于饱和的H原子），三种基组的和混合定义和两种基组的混合定义是一样的吗？先表示价层电子基组再表示赝势基组。（求组论坛的各位前辈）

非常谢谢

回复此楼