gaussian基本概念和用法(6)

3。无法收敛、优化和频率计算出错，这就很多了，比如L502，L716，L9999。这就只好具体问题具体分析了。 Range of M.O.s used for correlation: 20 167

NBasis= 167 NAE= 57 NBE= 57 NFC= 19 NFV= 0 NROrb= 148 NOA= 38 NOB= 38 NVA= 110 NVB= 110 Semi-Direct transformation.

ModeAB= 4 MOrb= 38 LenV= 38809382 LASXX= 49931981 LTotXX= 49931981 LenRXX= 101093547 LTotAB= 51161566 MaxLAS= 65294640 LenRXY= 0 NonZer= 151025528 LenScr= 265347674 LnRSAI= 65294640 LnScr1= 136751342 MaxDsk= -1 Total= 568487203 SrtSym= T ITran= 4

File extend in NtrExt1 failed; probably out of disk space. ************************* 呵呵，老大， LINK DIE容易找

但是，比如我上面的出错信息是什么意思啊？

硬盘空间不够了，哥们。 不过你说的对，像NtrExt1之类的，我的确看不懂：（ 我在用OPT的一个参数STEPSIZE的时候，从来没有成功过，不知道是为什么？ # UHF/6-31 OPT(Stepsize=30) ---------------------------

QPERR --- A SYNTAX ERROR WAS DETECTED IN THE INPUT LINE. # UHF/6-31 OPT(STEPSIZE=30) ' Last state=\TCursr= 9023 LCursr= 15

Error termination via Lnk1e in D:\\G98W\\l1.exe. Job cpu time: 0 days 0 hours 0 minutes 0.0 seconds.

File lengths (MBytes): RWF= 6 Int= 0 D2E= 0 Chk= 5 Scr= 1 这个贴子最后由helpme在 2004/03/29 04:58pm 第 1 次编辑]

我试了一下，OPT(STEPSIZE=30) 确实不可以用。（网上有人说过）可以用IOp(1/8)来实现。默认的IOP(1/8=30) IOp(1/8)

L103, L109, L112: MAXIMUM STEP SIZE ALLOWED DURING OPT.

0 DXMAXT = 0.1 BOHR OR RADIAN (L103, Estm or UnitFC). = 0.3 Bohr or Radian (L103, Read or CalcFC). = 0.2 Bohr or Radian (L105). = 0.3 Bohr or Radian (L113, L114). N DXMAXT = 0.01 * N

Davidson Disk Diagonalization: ConvIn= 1.00D-03 SkipCon=T Conv= 1.00D-03. Max sub-space: 200 roots to seek: 100 dimension of matrix: ****** Iteration 1 Dimension 100 NMult 100

PrmDim-Cnvrt allocation failure: iend,mxcore= 249095914 104856200

16

Error termination via Lnk1e in /usr/local/g03/l914.exe at Fri Feb 27 12:07:12 2004.

这个错误提示是什么意思，我是在SGI工作站上算的，应该没有硬盘问题，我给的内存是1G，NBasis是：1600。哪位大侠能说一下，我用的是zindo来算的。 半经验的方法不一定是内存越大越好

2，你用ZINDO算什么体系啊，很大吗？我也在用，是用来SCAN PES的，没有任何经验，如果您有什么感想，我们分享吧，呵呵

Optimization stopped.

-- Wrong number of Negative eigenvalues: Desired= 1 Actual= 0 -- Flag reset to prevent archiving.

负本征值数目不对，我做过渡态优化时也常遇到。我加上了noeigen，就解决了这个问题，结果也还可以。可是不知道根据是什么，请各位大虾指点

Gaussian98高级注释

Gaussian98高级注释

Mark A. Zottola，David C. Young（阿拉巴马超级计算中心） ****************** 目录

******************

计算电子激发态; 收敛问题的调整 ;自旋污染; 使用赝势; 耦合簇; Gaussian和CBS理论; 精度和CPU时间比较; 练习 计算电子激发态 有很多方法可以计算电子激发态，其中有些就在Gaussian中，虽然它并不是完成这项任务最好的程序。它们在精度、易用性和计算资源需求方面各不相同。全部考虑，有下面这些： 1. 通过改变多重度，计算不同于基态自旋的第一激发态。

2. CIS计算，使用CIS关键字。不用选择，自动包含Density关键字，所以布居数分析将反映激发态。对于闭壳层基态，你可以定义计算单重激发态还是三重激发态或是都计算 (CIS=Singlets, CIS=Triplets, CIS=50-50)。使用CIS=(Root=N)定义哪一个激发态用于几何优化，其中N 为1是第一激发态。

3. 对角区域哈密顿量允许定义每个对称性有多少电子态。但不能用于Gaussian。

4. 关键字ZINDO定义了ZINDO-1半经验方法。它的选项和上面CIS列出的相同。在Gaussian 98中，ZINDO不能用于几何优化。

5. 多重行列式计算中更高的根。在Gaussian中用于CASSCF，使用关键字CASSCF(NRoot=N)。Opt=Conical与CASSCF结合，寻找避免交叉或圆锥交叉点。关键字StateAverage可以给出更精确的激发能。关键字Spin在计算中包含两个态之间近似的自旋-轨道耦合。

6. 含时计算考虑两个态间的谐振，因而给出了得到激发态信息的方法。关键字TD 使用和上面CIS相同的选项。它可以用于HF和DFT计算。

7. 选择初始猜测。使用SCF=NoDIIS选项。

17

在Gaussian中，激发态优化不再利用冗余内坐标。因此，激发态优化对Z-矩阵结构非常敏感。 研究分子激发态的计算通常要多步处理，分别为 1. 寻找存在哪些激发态，以及要研究哪些。 2. 几何优化，如果需要的话。 3. 进行频率计算，以校验结构。 4. 特性的计算。 参见：

J. B. Foresman, ?. Frisch, \Chemistry with Electronic Structure Methods Second Edition\Gaussian, Pittsburgh (1996).

Adv. Chem. Phys. J. W. McGowan, Ed. vol. 45 (1981).

S. D. Peyerimhoff, R. J. Buenker, Adv. Quantum Chem. 9, 69 (1975). \

收敛问题的调整 如果SCF计算收敛失败，你首先会采取哪些技巧呢？这里是我们强烈推荐的首选方法。 1. 在Guess关键字中使用Core，Huckel或Mix选项，试验不同的初始猜测。

2. 对开壳层体系，尝试收敛到同一分子的闭壳层离子，接下来用作开壳层计算的初始猜测。添加电子可以给出更合理的虚轨道，但是作为普遍的经验规则，阳离子比阴离子更容易收敛。选项Guess=Read定义初始猜测从Gaussian计算生成的checkpoint文件中读取。

3. 另一个初始猜测方法是首先用小基组进行计算，由前一个波函得到用于大基组计算的初始猜测(Guess=Read自动进行)。 4. 尝试能级移动(SCF=Vshift)。

5. 如果接近SCF但未达到，收敛标准就会放松或者忽略收敛标准。这通常用于不是在初始猜测而是在平衡结构收敛的几何优化。SCF=Sleazy 放松收敛标准，Conver选项给出更多的控制。

6. 一些程序通过减小积分精度加速SCF。对于使用弥散函数，长程作用或者低能量激发态的体系，必须使用高积分精度：SCF=NoVarAcc。

7. 尝试改变结构。首先略微减小键长，接下来略微增加键长，接下来再对结构作一点改变。 8. 考虑使用不同的基组。

9. 考虑使用不同理论级别的计算。这并不总是实用的，但除此之外，增加迭代数量总是使得计算时间和使用更高理论级别差 …… 此处隐藏：2070字，全部文档内容请下载后查看。喜欢就下载吧 ……