A list of error messages and possible solutions -
Gaussian calculations can fail with various error messages. Some error messages from .out and .log files - and possible solutions - have been compiled here to facilitate problem solving.-
These are divided into:-
Syntax and similar errors-语法类错误
Memory and similar errors-内存类错误
Convergence problems -不收敛错误
Errors in solvent calculations -溶剂中的计算错误
Errors in log files-错误文件
-
ERROR MESSAGES IN OUTPUT FILES-
Syntax and similar errors:End of file in ZSymb.-
Error termination via Lnk1e in /global/apps/gaussian/g03.e01/g03/l101.exe Solution: The blank line after the coordinate section in the .inp file is missing. (输入文件空行丢失)
Unrecognized layer "X".-（不识别层X）
Error termination via Lnk1e in /global/apps/gaussian/g03.e01/g03/l101.exeSolution: Error due to syntax error(s) in coordinate section (check carefully). If error is "^M", it is caused by DOS end-of-line characters (e.g. if coordinates were written under Windows). Remove ^M from line ends using e.g. emacs. To process .inp files from command line, use sed -i 's/^M//' File.inp (Important: command does not work if ^M is written as characters - generate ^M on command line using ctrl-V ctrl-M).-
QPERR --- A SYNTAX ERROR WAS DETECTED IN THE INPUT LINE.-
Solution: Check .inp carefully for syntax errors in keywords -
RdChkP: Unable to locate IRWF=0 Number= 522.-
Error termination via Lnk1e in /global/apps/gaussian/g03.e01/g03/l401.exe or-
FileIO operation on non-existent file.-
[...] Error termination in NtrErr:-
NtrErr Called from FileIO.Solution: Operation on .chk file was specified (e.g.
geom=check, opt=restart), but .chk was not found. Check that:-
%chk= was specifed in .inp-
.chk has the same name as .inp-
.chk is in the same directory as .inp -
run script transports .chk to temporary folder upon job start. Run scripts downloaded here should do this. -
The combination of multiplicity N and M electrons is impossible.-（多重性）Error termination via Lnk1e in /global/apps/gaussian/g03.e01/g03/l301.exeSolution: Either the charge or the multiplicity of the molecule was not specified correctly
in .inp.-（电荷和多重性指定错误）
Memory and similar errors: Out-of-memory error in routine RdGeom-1 (IEnd= 1200001 MxCore= 2500)-
Use %mem=N MW to provide the minimum amount of memory required to complete this step-
Error termination via Lnk1e in /global/apps/gaussian/g03.e01/g03/l101.exe or-
Not enough memory to run CalDSu, short by 1000000 words.-
Error termination via Lnk1e in /global/apps/gaussian/g03.e01/g03/l401.exe or-
[...] allocation failure: -（表示配分失败）
Error termination via Lnk1e in /global/apps/gaussian/g03.e01/g03/l1502.exe Solution: Specify more memory in .inp (%mem=Nmb). Possibly, also increase pvmem value in run script. Especially solvent calculations can exhibit allocation failures and explicit amounts of memory should be specified.-
galloc: could not allocate memory.-（无法分配内存）
Solution: The %mem value in .inp is higher than pvmem value in run script. Increase pvmem or decrease %mem. -
Probably out of disk space(磁盘空间). Write error in NtrExt1 Solution: /scratch space is most likely full. Delete old files in temporary folder. -
Convergence problems: Density matrix is not changing but DIIS error= 1.32D-06 CofLast= 1.18D-02.-（收敛问题）
The SCF is confused. Error termination via Lnk1e in
/global/apps/gaussian/g03.e01/g03/linda-exe/l502.exel Solution: Problem with DIIS. Turn it off completely, e.g. using SCF=qc, or partly by using
SCF=(maxconventionalcycles=N,xqc), where N is the number of steps DIIS should be used (see SCF keyword). -
Convergence criterion not met. SCF Done: E(RHF) = NNNNNNN A.U. after 129 cycles -
[...] Convergence failure -- run terminated. Error termination via Lnk1e in
/global/apps/gaussian/g03.e01/g03/linda-exe/l502.exe Solution: One SCF cycle has a default of maximum 128 steps, and this was exceeded without convergence achieved. Possible solution: In the route section of input file, specify SCF=(MaxCycle=N), where N is the number of steps per SCF cycles. Alternatively, turn of DIIS (e.g. by SCF=qc) (see SCF keyword).-
-
Problem with the distance matrix.-（距离矩阵）
Error termination via Lnk1e in /pkg/gaussian/g03/l202.exe Solution: Try to restart optimization from a different input geometry. -（重新不同几何异构体的输入优化）
New curvilinear step not converged（新曲线步骤不收敛）. Error imposing constraints-
Error termination via Lnk1e in /pkg/gaussian/g03/l103.exe-
Solution: Problem with constrained coordinates (e.g. in OPT=modredun calculation). Try to restart optimization from a slightly different input geometry. -（一种稍微不同的输入几何）
-
Optimization stopped. -- Number of steps exceeded, NStep= N-
[..] Error termination request processed by link 9999.-
Error termination via Lnk1e in /global/apps/gaussian/g03.e01/g03/l9999.exe Solution: Maximum number of optimization steps is twice the number of variables to be optimized. Try increasing the value by specifying OPT=(MaxCycle=N) in .inp file, where N is the number of optimization steps (see OPT keyword). Alternatively, try to start optimization from different geometry.-
-
Errors in solvent calculations: AdVTs1: ISph= 2543 is engulfed by JSph= 2544 but Ae( 2543) is not yet zero!-
Error termination via Lnk1e in /global/apps/gaussian/g03.e01/g03/l301.exe Solution: Problem is related to building of the cavity in solvent calculations（溶剂效应优化计算错误）. One possible solution is to change the cavity（腔） model (default in g03 is UAO, can be changed by adding RADII keyword in section below coordinates in the .inp file, e.g. RADII=UFF, see SCRF keyword).-
-
Hydrogen X has 2 bounds. Keep it explicit at all point on the-
potential energy surface to get meaningful results.Solution: In UAO cavity model, spheres are placed on groups of atoms, with hydrogens assigned to the heavy atom, they are bound to. If assignment fails (e.g. because heavy atom-H bond is elongated), cavity building fails. Possible solutions: a) use cavity model that also assigns spheres to hydrogens (e.g. RADII=UFF) or b) Assign a sphere explicity on problematic H atom (use SPHEREONH=N, see SCRF keyword)-
-
ERROR MESSAGES IN LOGFILES =>> PBS: job killed: wall time N exceeded limit M-
signal number 15 received. Solution: Job did not finish within specified wall time. Retrieve .out and .chk files from temporary folder /global/work/$USER/$JOB (or $PBS_JOBID) and restart calculation if possible (using e.g. opt=restart or scf=restart). -
cp: cannot stat $JOB.inp: No such file or directory Solution: The .inp file is not in the directory from where the job was submitted (or its name was misspelled during submission. If error reads: cp: cannot stat $JOB .inp .inp, the .inp file was submitted with extension).-
ntsnet: unable to schedule the minimum N workers Solution: The value of %N proc Linda=N in the .inp file is higher than the number of nodes asked for during submission. Make sure these values match.
Connection refused [...] died without ever signing in-
Sign in timed out after 0 worker connections. Did not reach minimum (N), shutting downSolution: Error appears if you run parallel calculations but did not add this file to your $HOME directory: .tsnet.config containing only the line: Tsnet.Node.lindarsharg: ssh (see also guidelines for submission). -
Density matrix is not changing but DIIS error - Suggested solutions
1/
- SCF=qc will probably solve the problem, albeit at a cost
- Change the SCF converger to either SD, Quadratic or Fermi
2/
- lower the symmetry of optimize with and optimize
with the "nosymm" keyword
I solved the problem using a variation on the first suggestion. Normally the scf took less than 80 cycles to converge. So i used scf=(Maxconventionalcycles=100,xqc) which resulted in a good compromise between using scf=qc and optimisation speed.
In the case of the DIIS error the scf always took more than 100 cycles before the error,
so by adding scf=(Maxconventionalcycles=100,xqc) the scf switched to qc after 100 cycles in the standard DIIS mode.
l9999错误是优化圈数不够，把out文件保存成gjf，修改后接着优化。这样比较省事。或者在原来输入文件opt中加入maxcyc=500,把优化圈数加大到500
L1002错误贋势基组不对，换基组，或者用guess=Indo
L9999
Optimization stopped.
-- Wrong number of Negative eigenvalues: Desired= 1 Actual= 3
-- Flag reset to prevent archiving.
这是你原来贴的出错信息。
显然，gaussian检测了虚频个数，发现是3个而非1个——过渡态有且只有一个虚频，从而终止，显示出错。
加上noeigentest的意思是在优化过程中不必每步都去检测虚频个数，而是按照最可能的方式继续寻找下去，到跑完为止，所以如果你对过渡态的初始构型不是非常有把握的情况下，加上noeigen是比较不错的选择。
但你现在继续出状况，却没有贴新的出错信息，我没法猜这次是什么问题。。。
同时你那个命令行# b3lyp/6-31++g(d,p) opt=(ts,calcfc,noeigen) freq test
改成#p b3lyp/6-31++g(d,p) opt=(ts,calcfc,noeigen) freq
也就是任何时候#后面都要加一个p，结尾那个test不要写，被人笑话的。Inaccurate quadrature in CalDSu 错误解决方法
已有 3175 次阅读 2012-7-4 16:27 |个人分类:Gaussian|系统分类:科研笔记
Inaccurate quadrature in CalDSu 错误解决方法
g09/g03 输出文件最后几行大致如下：
Spurious integrated density or basis function:
NE= 149 NElCor= 0 El error=2.83D-03 rel=1.89D-05 Tolerance=1.00D-03
Shell 67 absolute error=2.30D-02 Tolerance=1.20D-02
Shell 67 signed error=2.30D-02 Tolerance=1.00D-01
Inaccurate quadrature in CalDSu.
Error termination via Lnk1e in /home/soft/g03/l502.exe at Mon Feb 23 09:42:55 2009. Job cpu time: 0 days 0 hours 0 minutes 3.1 seconds.
File lengths (MBytes): RWF= 13 Int= 0 D2E= 0 Chk= 13 Scr= 1
解决方法
1. 对某些分子，可以设置guess=indo。(Guess=indo的意思是：使用Gaussian98的默认初始猜测：对第一行元素是INDO，第二行元素是CNDO，第三行及以后的元素是Huckel。如果不写Guess=indo，则使用Gaussian03的默认初始猜测：全部使用Huckel初始猜测。Guess=indo的用处：使比较差的初始构型，能比较顺利地通过初始猜测。)
2. 对于对称分子，设置scf(dsymm) 强制使用密度对称。
3. 对于对称分子，设置全局nosymm，降低对称性。
4. 使用scf(novaracc)也可能会产生或消除这个错误。
5. 加上int(untrafine)，不过这个改变了积分网格，这个计算出的能量不能与同类计算直接比较。
6. 最佳解决方法：加关键词int(NoXCTest)网格问题(G09默认的精度检测方法不适用,加上Int=NoXCTest), Skip tests of numerical accuracy of XC quadrature (g09手册)。需要检测结果波函数是否正确。可同时加guess=indo。
Error in internal coordinate system.
可在opt中加入cartesian
例如：
%chk=TS3-D-itm3-eoc-4.chk
%mem=45000MB
%NProcShared=8
#B3LYP/6-31G** iop(5/13=1) iop(2/11=1)
#scf=(maxcycle=200)
#opt=(cartesian,TS,calcfc,noeigen,maxcycle=200) freq
#geom=allcheck
Error in internal coordinates
c r e a t e d:2009-07-3114:38:48
T h i s t o p i c c o v e r s h o w t o s o l v e t h e"E r r o r i n i n t e r n a l c o o rd i n a t e s"p ro b l e m.
S o m e t i m e s G a u s s i a n q u i t s a t t h e B e r n y o p t i m i z a t i o n s t a g e wi t h t h e"Er r o r i n
i n t e r n a l c o o rd i n a t e s"m e s s a g e.T h i s h a p p e n s a s B e r n y o p t i m i z a t i o n b y
d e f a u l t u s e s r e d u n d a n t i n t e r n a l c o o rd i n a t e s.Ca r t e s i a n c o o rd s a r e e a s y t o
d e f i n e a n d u s e b u t c a n b e s t ro n g l y c o u p l e d t o o n e a n o t h e r.I n t e r a l c o o r d s
d e s c r i b e m o l e c u l a r p r o p e r t i e s(b o n d l e n g t h s,a n g l e...)n a t u r a l l y a n d t h u s
i n v o l v e l e s s c o u p l i n g.Th e r e a re a f e w t y p e s o f r e d u n d a n t i n t e r n a l c o o r d s
wh i c h c a n a c c e l e r a t e o p t p r o c e s s.Bu t wh e n a m o l e c u l e h a s f l a t g e o m e t r y,i t
c o u l
d b
e h a r d t o c a l c u l a t e t h e r e v e r s e m a t r i x o
f r e d u n d a n t i n t e r n a l c o o r d s.
T h e r e f o r e,i t i s b e t t e r t o s wi t c h b a c k t o C a r t e s i a n(o p t=(...,C a r t e s i a n,...))o r u s e"n o s y m m".
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad Berny optimization.
NTrRot= -1 NTRed= 29 NAtoms= 10 NSkip= 5 IsLin=F
Error in internal coordinate system.
Error termination via Lnk1e in /home/cast/program/g09/l103.exe at Tue Nov 23
13:44:47 2010.
本人计算时出现如上错误，哪位达人知道是怎么回事？我觉得我的结构应该没问题
加上关键词OPT=Cartesian！
高斯对不收敛问题的对策
首先，我们必须理解收敛是什么意思。在自洽场(SCF)计算中，自洽循环中，首先产生一个轨道占据的初始猜测，
1)然后根据此轨道占据构造电荷密度和哈密顿量。
2)对角化哈密顿量，得到新的轨道能级和占据。
3)产生新的电荷分布和哈密顿量，重复步骤2）
经过一定次数的循环后，某次循环前和循环后的电荷密度差别小于一定的标准，我们称之为收敛。
如果以上过程不能收敛，则gaussian给出convergence failure的警告。
如果SCF计算收敛失败，你首先会采取哪些技巧呢？这里是我们强烈推荐的首选方法。
1 考虑使用更小的基组
由于一定的基组对应于一定精度和速度，所以更换基组并不在所有的情况下都适用。方法是首先用小基组进行计算，由前一个波函得到用于大基组计算的初始猜测(Guess=Read自动进行)。
2 增加最大循环步数
Gaussian默认的最大循环步数为64 (SCF=DM或SCF=QC方法则为512)，如果循环次数超过这个数目则会汇报convergence failure。在一定的情况下，不收敛的原因仅仅是因为最大循环步数不够。可以通过设置maxcyc来增大最大循环步数。更多的SCF迭代
(SCF(MaxCycle=N)，其中N是迭代数)。这很少有帮助，但值得一试。
3
放宽收敛标准
如果接近SCF但未达到，收敛标准就会放松或者忽略收敛标准。这通常用于不是在初始猜测而是在平衡结构收敛的几何优化。SCF=Sleazy放松收敛标准，Conver选项给出更多的控制。在Gaussian中设置收敛标准的命令行是Conv=8 表示10-8,你可以将这个标准降低，比如改为10-5等等。微调结构是比较有效的方法，能量收敛标准不宜太低，最好不低于10-7.
4
尝试改变初始构型
首先略微减小键长，接下来略微增加键长，接下来再对结构作一点改变。
5
尝试能级移动Level shifting (SCF=Vshift)
如果不收敛的原因是波函数的震荡行为，通常是因为在相近的能量上的泰的混合。对于这种情况，我们可以采用level shifting的方法。Level shifting的含义是人工的升高非占据轨道的能级，以防止和最高占据轨道之间的混合，以达到收敛的目的。在Gaussian中此方法的关键词为SCF=Vshift
6
使用强制的收敛方法SCF=QC
SCF=QC通常最佳，但在极少数情况下SCF=DM更快。此关键字将大大增加计算时间，但是收敛的机会更大。不要忘记给计算额外增加一千个左右的迭代。应当测试这个方法获得的波函，保证它最小，并且正好不是稳定点(使用Stable关键字)。
7
对开壳层体系，尝试收敛到同一分子的闭壳层离子，接下来用作开壳层计算的初始猜测。添加电子可以给出更合理的虚轨道，但是作为普遍的经验规则，阳离子比阴离子更容易收敛。选项Guess=Read定义初始猜测从Gaussian计算生成的checkpoint文件中读取。
8
一些程序通过减小积分精度加速SCF。对于使用弥散函数，长程作用或者低能量激发态的体系，必须使用高积分精度：SCF=NoVarAcc。
9
改变模型或方法
可以考虑改变模型方法。比较常见的方法有HF,GVB,MCSCF,CASSCF,MPn等。改变模型方法通常也会收敛性质。通常，精度更高的方法更难收敛。精度比较低的方法产生的计算结果可以作为高精度计算的初始猜测。考虑使用不同理论级别的计算。这并不总是实用的，但除此之外，增加迭代数量总是使得计算时间和使用更高理论级别差不多。
10
关闭DIIS外推(SCF=NoDIIS)。同时进行更多的迭代(SCF=(MaxCycle=N))。
11
试着改用DIIS之外其它方法(SCF=SD或SCF=SSD)。
Gaussian不收敛的可能原因及对策：
1 由于体系有很多能量相近的能级，导致计算不收敛。
如果计算中采用的是Hartree-Fock方法或者其他的混合形式的交换相关势（如
B3LYP）,则可以尝试一下的方法进行改进。
scf=(noincfock,conver=11,maxcyc=1025) iop(5/22=20)
用以上的命令行将强制Gaussian采用EDIIS的算法。此算法计算量更大，但是更加稳定。在命令行里也增加了最大循环数，以增加收敛的到基态的可能性。
2．检查是否有初始文件错误
常见初级错误：
a. 自旋多重度错误
b. 变量赋值为整数
c. 变量没有赋值或多重赋值
d. 键角小于等于0度，大于等于180度
e. 分子描述后面没有空行
f. 二面角判断错误，造成两个原子距离过近
g. 分子描述一行内两次参考同一原子,或参考原子共线
3． SCF(自洽场)不收敛
则一般是L502错误，省却情况做64个cycle迭代(G03缺省128 cycles)
a. 修改坐标，使之合理
b. 改变初始猜
Guess=Huckel 或其他的，看Guess关键词。
c. 增加叠代次数 SCFCYC=N （对小分子作计算时最好不要增加，很可能结构不合理）
d. iop（5/13=1）这样忽略不收敛，继续往下做。
4．分子对称性改变
a. 修改坐标，强制高对称性或放松对称性
b. 给出精确的、对称性确定的角度和二面角。如CH4的角度给到109.47122
c. 放松对称性判据
Symm=loose
d. 不做对称性检查 iop(2/16=1) （最好加这个选项）
iop(2/16=2) 则保持新的对称性来计算
5． Opt时收敛的问题
a. 修改坐标，使之合理
b. 增加叠代次数optcyc=N
6．优化过渡态，若势能面太平缓，则不好找到。
iop(1/8=10) 默认30（下一个结构和该结构的差别0.3?），可改成10。如果每一步都要用到小的步长，应该加opt(notrustupdate)
7．在CI（组态）方法中
如QCISD(T),CCSD(T),CID方法中，省却最大循环50，若出错（L913错误）
解决方法：#P QCISD(maxcyc=N) 注：N≤512
8．优化过渡态
opt=TS (给出过渡态)
opt=qst2 (给出反应物和产物)
opt=qst3 (给出反应物和产物和过渡态)
a.
用G03时的出错opt=ts 必须加FC (force constant)
写法：opt=(TS, calcFc)
or opt=(TS,calchffc)计算HF力常数，对QCISD，CCSD等方法用；
or opt=(TS,modRedundant)
(最好写这个)
b. 如果计算采用QCISD计算（不好计算FC）
则写为QCISD opt=(TS, calcHFFC) (用HF计算FC)
9. 无法写大的Scratch文件RWF
a. 劈裂RWF文件 %rwf=loc1,size1,loc2,size2,……..,locN,-1
b. 改变计算方法 MP2=Direct可以少占硬盘空间
c. 限制最大硬盘 maxdisk=N GB，****MB,有些系统写2GB会出错，可以写2000MB
10. FOPT出错原因是变量数与分子自由度数不相等。可用POPT 或直接用OPT
11. 优化过渡态只能做一个STEP 原因是负本征数目不对添加 iop（1/11）=1或者noeigentest (eigentest 是表示优化过渡态检测分子振动的本征值，过渡态只有一个负值，但优化的时候往往出现两个或者更多的情况，默认的是如果出现多于一个的情况就停止优化计算，这往往时不必要的，noeigentest表示优化时不坐此检测，实际上优化过渡态的时候出现两个三个很正常的，只是第一个一般负的很大，接着的比较小，不会影响过渡态的搜寻，当然有时候需要自己判断是不是接近真实的过渡态。其实经验上也是过渡态的本征虚频越小，相对越难找。)
对于 (L502, L508, L9999)出错的对策
对于一个优化计算，它的过程是先做一个SCF计算，得到这个构型下的能量，然后优化构型，再做SCF，然后再优化构型。。。因此，会有两种不收敛的情况：一是在某一步的SCF 不收敛(L502错误)，或者构型优化没有找到最后结果（L9999错误）。
预备知识:计算时保存chk文件,可以在后续计算中使用guess=read读初始猜测.
对于SCF不收敛，通常有以下的解决方法：
1. 使用小基组，或低级算法计算，得到scf收敛的波函数，用guess=read读初始波函数。
2. 使用scf=qc，这个计算会慢，而且需要用stable关键字来测试结果是否波函数稳定。如果这个还不收敛，会提示L508错误。
3. 改变键长，一般是缩小一点，有时会有用。
4. 计算相同体系的其他电子态，比如相应的阴离子、阳离子体系或单重态体系，得到的收敛波函数作为初始猜测进行计算。
对于优化不收敛,即L9999错误，实际上是在规定的步数内没有完成优化，即还没有找到极小值点。（或者对于过渡态优化，还没有找到过渡态）
这有几种可能性：
1. 看一下能量的收敛的情况，可能正在单调减小，眼看有收敛的趋势，这样的情况下，只要加大循环的步数（opt(maxcycle=200)），可能就可以解决问题了。
2. 加大循环步数还不能解决的（循环步数有人说超过200再不收敛，再加也不会有用了，这虽然不一定绝对正确，但200步应该也差不多了），有两种可能。一是查看能量，发现能量在振荡了，且变化已经很小了，这时可能重新算一下，或者构型稍微变一下，继续优化，就可以得到收敛的结果（当然也有麻烦的，看运气和经验了）；二是构型变化太大，和你预计的差别过大，这很可能是你的初始构型太差了，优化不知道到哪里去了，这时最好检查一下初始构型，再从头优化。
3. 对于L9999快达到收敛时,考虑减小优化步长有时对于能量振荡的情况也是有用的，opt(maxstep=1).(flyingheart )
https://www.doczj.com/doc/272790861.html,/chemistry/resources/messages/2002/10/16.005-dir/index.html
Hello,
here is the summary of all the answers I got for my problem. Thanks (in alphabetically order) to Laurence Cuffe, Kerwin Dobbs, Paul Kiprof, Liviu Mirica, Cory Pye and Dave Shobe.
The problem was solved by increasing the number of SCF cycles, as
suggested by many. This seems very often necessary with transition metals. The two other suggestions (using first a small basis set before switching
to the larger, or turning off "vshift", did not work, unfortunately.
Thanks again,
R.Herrmann.
The original question:
---------------------------------------------------------------------------- Hello,
I recently tried to calculate the NMR chemical shifts of VOCl3
with Gaussian 98:
---------------------
#B3LYP/6-31G* SP NMR
VOCl3
0 1
cl
cl 1 clcl2
cl 2 clcl3 1 clclcl3
o 3 ocl4 2 oclcl4 1 dih4
v 3 vcl5 2 vclcl5 1 dih5
clcl3 3.506243
clclcl3 59.999
ocl4 2.938735
oclcl4 53.376
dih4 -64.586
vcl5 2.118170
vclcl5 34.144
dih5 -31.640
-----------------------
The result (even with several variations in the input) is always an
error termination when trying the "incremental Fock formation"
(convergence criterion not met). The structure was optimized with
GAMESS.
Any suggestion what I could do is wellcome - thank you!
------------------------------
108 basis functions 288 primitive gaussians
41 alpha electrons 41 beta electrons
nuclear repulsion energy 561.9817715949 Hartrees.
One-electron integrals computed using PRISM.
NBasis= 108 RedAO= T NBF= 108
NBsUse= 108 1.00D-04 NBFU= 108
Projected Huckel Guess.
Requested convergence on RMS density matrix=1.00D-08 within 64 cycles.
Requested convergence on MAX density matrix=1.00D-06.
Virtual orbitals will be shifted by 0.200 hartree.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
>>>>>>>>>> Convergence criterion not met.
SCF Done: E(RB+HF-LYP) = -2359.59549049 A.U. after 65 cycles
Convg = 0.5851D-
03 -V/T = 1.9657
S**2 = 0.0000
Convergence failure -- run terminated.
Error termination via Lnk1e in /usr/local/g98_a7/l502.exe.
Job cpu time: 0 days 0 hours 16 minutes 33.0 seconds.
File lengths (MBytes): RWF= 14 Int= 0 D2E= 0
Chk= 1 Scr= 1
-------------------------------------------------------------------------- The answers:
1. From: Kerwin D Dobbs
To: "Rudolf Herrmann
Subject: Re: CCL:G98: VOCl3 problem: convergence criterion not met
Dr. Herrmann,
Try turning off the "vshift" by changing the line,
#B3LYP/6-31G* SP NMR
to
# B3LYP/6-31G* SP NMR IOP(5/10=-1)
When there is a TM in the molecule, G98 automatically turns
on "VSHIFT" which usually makes SCF convergence worse. The
above addition to your line turns off "VSHIFT". Give it a try and
see what happens.
Kerwin
------------------------------------------------------------------
2. From: Cory Pye
To: Rudolf Herrmann
Subject: Re: CCL:G98: VOCl3 problem: convergence criterion not met
Hello,
your SCF is not converging. Try running, at first, HF/STO-3G. If it works, then
Then try B3LYP/6-31G* afterwards with guess=read in the input. It usually works like a charm for me, I always use something like this.
-Cory
----------------------------------------------------------------------
3. From: Liviu Mihail Mirica
To: Rudolf Herrmann
Subject: Re: CCL:G98: VOCl3 problem: convergence criterion not met
Try adding SCFcyc=1000in your command lines. This will increase the # of SCF cycles up to the desired number (1000 in this case). that should be enough for your system.
------------------------------------------------------------------------ 4. From: Laurence Cuffe
To: Rudolf Herrmann
Subject: Re: CCL:G98: VOCl3 problem: convergence criterion not met
And it should run these keywords change the default SCF method
and then increase the number of SCF cycles allowed. I got the job
to complete in about 50 minutes running it in core on a linux box
with a memory limit of 600MB but it this doesn't make sense to
you don't worry about it. Try using the modified line above and if
there are any problems don't hesitate to come back to me.
All the best
Laurence Cuffe
-------------------------------------------------------------------------
5. From: Paul Kiprof
To: herrmann-: at :-itqb.unl.pt
Subject: VOCL3
Dear Dr. Herrmann,
I would add to the input
SCF(direct,qc,maxcycle=1000)
This will increase the number of convergence cycles and hopefully converge
the calculation. I always do that when I have transition metal compounds.
Best wishes,
Paul Kiprof
---------------------------------------------------------------------------
6. From: "Shobe, Dave"
To: 'Rudolf Herrmann'
Subject: RE: G98: VOCl3 problem: convergence criterion not met
The standard 64 SCF "cycles"(or iterations of the variational
procedure) is rarely enough for transition metal compounds, even closed shell ones like
VOCl3. Read the section in the Gaussian manual describing the "scf" keyword, and you will find various options which may improve the
convergence. A lot of it is trial and error, with these different scf options. The simplest is scf(maxcyc=999)which increases the number of cycles/iterations to 999 before the program "gives up" and prints the
error message. Scf=qc and scf(vshift=###) have also been useful in my experience.
More advice: once you get a good wavefunction, make a backup of the
checkpoint file, so if the wave function becomes corrupted in subsequent calculations you can always go back to it.
----------------------------------------------------------------------------
--
Dr. Rudolf Herrmann
Instituto de Tecnologia Quimica e Biologica (ITQB)
Av. da Republica, EAN, Apt. 127
P-2781-901 Oeiras
e-mail: herrmann-: at :-itqb.unl.pt
－－－－－－－－－－－－－－用Gaussian计算不收敛的问题
问题一：
Problem detected with inexpensive integrals.
Switching to full accuracy and repeating last cycle.
Cycle 2 Pass 1 IDiag 1:
Warning! Spurious integrated density:
NE= 188 NElCor= 0 Integral= -76.06821 Tolerance=1.00D-03
Consistency failure #2 in CalDSu.
Error termination via Lnk1e in /home/gaussian/g98/l502.exe.
Job cpu time: 1 days 3 hours 21 minutes 30.2 seconds.
File lengths (MBytes): RWF= 118 Int= 0 D2E= 0 Chk= 17 Scr= 1
解决：(from ccl)
As far as I know it's a dft-related error. Gaussian checks the charge
resulting from the integration of the density to the number of electrons, and in case of a large difference prints this eror message. Try a denser grid with the "integral" keyword.
可是我采用integral（Grid=ultrafine ）网格仍是出现以下错误信息：
Raffenetti 2 integral format.
Two-electron integral symmetry is turned off.
306 basis functions 778 primitive gaussians
51 alpha electrons 51 beta electrons
nuclear repulsion energy 865.0704012775 Hartrees.
One-electron integrals computed using PRISM.
1 Symmetry operations used in ECPInt.
ECPInt: NShTT= 17391 NPrTT= 123061 LenC2= 14761 LenP2D= 57336.
LDataN: DoStor=F MaxTD1= 3 Len= 28
LDataN: DoStor=T MaxTD1= 3 Len= 28
NBasis= 306 RedAO= T NBF= 306
NBsUse= 306 1.00D-04 NBFU= 306
Projected Huckel Guess.
of initial guess= 0.0000
Requested convergence on RMS density matrix=1.00D-06 within 800 cycles.
Requested convergence on MAX density matrix=1.00D-04.
Integral accuracy reduced to 1.0D-05 until final iterations.
Warning! Spurious integrated density:
NE= 102 NElCor= 0 Integral= 102.20958 Tolerance=1.00D-03
Consistency failure #2 in CalDSu.
Error termination via Lnk1e in /usr/people/gaussian/g98/l502.exe.
Job cpu time: 0 days 0 hours 13 minutes 24.7 seconds.
File lengths (MBytes): RWF= 90 Int= 0 D2E= 0 Chk= 2 Scr= 1
请问各位大侠有没有什么解决方法？
问题二：
>>>>>>>>>> Convergence criterion not met.
SCF Done: E(UB-PW91) = -1306.37542885 A.U. after 801 cycles
Convg = 0.1221D+01 -V/T = 2.0190
S**2 = 0.0000
KE= 1.281978320748D+03 PE=-8.302694580633D+03 EE= 3.210822249139D+03
Annihilation of the first spin contaminant:
S**2 before annihilation 0.0000, after 0.0000
Convergence failure -- run terminated.
Error termination via Lnk1e in /usr/people/gaussian/g98/l502.exe.
Job cpu time: 2 days 2 hours 6 minutes 40.2 seconds.
File lengths (MBytes): RWF= 97 Int= 0 D2E= 0 Chk= 19 Scr= 1
一般而言，解决收敛的问题有以下几种，第一，波函数初猜，不稳定的波函数，积分难以收敛是肯定的，解决方案是作stable ＝
opt，做完以后，在此基础上从check文件中读入波函数。
－－link1－－
%chk
#p hf opt guess＝read geom＝check
第二，过渡金属原子会出现d轨道分裂问题，如果分裂能不够大，电子显然会在高能轨道
和低能轨道之间跳动，能量会上下波动，难以收敛，
解决方案是加入vshift关键词
第三，氢键，过渡金属，默认的129cycle所难以达到收敛，所以加大scf(cyc=500)可能可
以解决问题。
第四，初始结构一直容易被人所忽略，所以很多方法没有解决问题的时候，还是希望能
仔细检查一遍结构
第五， scf的时候，检查一下力常数也是好的，通常的解决方案作一下calcfc
对于第二个收敛的问题，我正在用SCF(vshift,QC) 尝试，得过几天才知道行不行。至少我已尝试了将内坐标改为直角坐标；加大循环次数（800仍不收敛）从check文件中读入，这些都不行。
对于第一个问题，我是用massage考虑分子间相互作用，可是好多都出现第一个问题，根本不让你算，用了integral(grid=ultrafine) 还是不行，输入文件如下：（这是一个不含金属原子的，有好几个含金属的输入，也是出现这个情况。
芝麻582：
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. 考虑使用不同理论级别的计算。这并不总是实用的，但除此之外，增加迭代数量总是使得计算时间和使用更高理论级别差不多。 10. 关闭DIIS外推(SCF=NoDIIS)。同时进行更多的迭代( SCF=(MaxCycle=N) )。
11. 更多的SCF迭代( SCF(MaxCycle=N)，其中N是迭代数)。这很少有帮助，但值得一试。
12. 使用强制的收敛方法。SCF=QC通常最佳，但在极少数情况下SCF=DM 更快。不要忘记给计算额外增加一千个左右的迭代。应当测试这个方法获得的波函，保证它最小，并且正好不是稳定点(使用Stable关键字)。
13. 试着改用DIIS之外其它方法(SCF=SD或SCF=SSD)。
L9999错误，实际上是在规定的步数内没有完成优化，即还没有找到极小值点。（或者对于过渡态优化，还没有找到过渡态）
这有几种可能性：
1。看一下能量的收敛的情况，可能正在单调减小，眼看有收敛的趋势，这样的情况下，只要加大循环的步数（opt(maxcycle=200)），可能就可以解决问题了。
2。加大循环步数还不能解决的（循环步数有人说超过200再不收敛，再加也不会有用了，这虽然不一定绝对正确，但200步应该也差不多了），有两种可能。一是查看能量，发现能量在振荡了，且变化已经很小了，这时可能重新算一下，或者构型稍微变一下，继续优化，就可以得到收敛的结果（当然也有麻烦的，看运气和经验了）；二是构型变化太大，和你预计的差别过大，这很可能是你的初始构型太差了，优化不知道到哪里去了，这时最好检查一下初始构型，再从头优化。
一个建议是，对于大体系，难收敛体系，先用小基组，低精度算法优化一下，以得到较高的初始构型，再用高精度的计算接着算。如果前面的方法保留了chk文件，重新计算时需要使用 geom=allcheck 读入构型（就不必麻烦地写构型了）， guess=read（读入初始波函数，可以加快第一步SCF收敛）。
确实是不收敛！有几种方法参考：1、添加关键词scf=qc；2、添加关键词scf=vshift或scf=vshift=n （n取1~200），有时候会成功收敛，但成功率不高，不常用；3、改变初始结构（楼上的已说，常用）；4、改用不同的方法（常用）5、增加计算步长（也不常用）
我就用过这些，仅供参考！希望能有帮助.
使用了SCF=QC后，计算量要远远超过同等级的非QC计算。
另外，QC计算得到的结果还需要进行波函数稳定性分析，这也是为什么不推荐
使用QC的原因之一，实在要用，可以参照我前面的帖子用
scf=(maxcycle=n,xqc), n是一个你认为合适的数值，推荐不要超过100，这样能节省绝大部分的计算量，因为这个参数只在必要时使用qc。
承诺书 我们仔细阅读了中国大学生数学建模竞赛的竞赛规则. 我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。 我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。 我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。 我们参赛选择的题号是（从A/B/C/D中选择一项填写）： 我们的参赛报名号为（如果赛区设置报名号的话）： 所属学校（请填写完整的全名）： 参赛队员(打印并签名) ：1. 2. 3. 指导教师或指导教师组负责人(打印并签名)： 日期：年月日赛区评阅编号（由赛区组委会评阅前进行编号）：
编号专用页 赛区评阅编号（由赛区组委会评阅前进行编号）： 全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：
基于系统综合评价的城市表层土壤重金属污染分析 摘要 本文针对城市表层土壤重金属污染问题，首先对各重金属元素进行分析，然后对各种重金属元素的基本数据进行统计分析及无量纲化处理，再对各金属元素进行相关性分析，最后针对各个问题建立模型并求解。 针对问题一，我们首先利用EXCEL 和 SPSS 统计软件对各金属元素的数据进行处理，再利用Matlab 软件绘制出该城区内8种重金属元素的空间分布图最后通过内梅罗污染 模型：2 /12 max 22?? ? ? ??+=P P P 平均综，其中平均P 为所有单项污染指数的平均值，max P 为土壤环境中 针对问题二,我们首先利用EXCELL 软件画出8种元素在各个区内相对含量的柱状图，由图可以明显地看出各个区内各种元素的污染情况，然后再根据重金属元素污染来源及传播特征进行分析，可以得出工业区及生活区重金属的堆积和迁移是造成污染的主要原因，Cu 、Hg 、Zn 主要在工业区和交通区如公路、铁路等交通设施的两侧富集，随时间的推移，工业区、交通区的土壤重金属具有很强的叠加性，受人类活动的影响较大。同时城市人口密度，土地利用率，机动车密度也是造成重金属污染的原因。 针对问题三，我们从两个方面考虑建模即以点为传染源和以线为传染源。针对以点为传染源我们建立了两个模型：无约束优化模型()[]()[]() 22y i y x i x m D -+-=，得到污染源的位置坐标()6782,5567；有衰减的扩散过程模型得位置坐标（8500,5500），模型为： u k z u c y u b x u a h u 222 2222222-??+??+??=??， 针对以线为传染源我们建立了l c be u Y ?-+=0模型，并通过线性拟合分析线性污染源的位置。 针对问题四，我们在已有信息的基础上，还应收集不同时间内的样点对应的浓度以及各污染源重金属的产生率。根据高斯浓度模型建立高斯修正模型，得到浓度关于时间和空间的表达式ut e C C -?=0。 在本题求解过程中，我们所建立的模型与实际紧密联系，有很好的通用性和推广性。但在求点污染源时，我们假设只有一个污染源，而实际上可能有多个点污染源，从而使得误差增大，或者使污染源的位置够不准确。 关键词 内梅罗污染模型 无量纲化 相关性 回归模型 高斯浓度模型
高斯收敛问题 首先，我们必须理解收敛是什么意思。在自洽场(SCF)计算中，自洽循环中，首先产生一个轨道占据的初始猜测， 1)然后根据此轨道占据构造电荷密度和哈密顿量。 2)对角化哈密顿量，得到新的轨道能级和占据。 3)产生新的电荷分布和哈密顿量，重复步骤2） 经过一定次数的循环后，某次循环前和循环后的电荷密度差别小于一定的标准，我们称之为收敛。 如果以上过程不能收敛，则gaussian给出convergence failure的警告。 如果SCF计算收敛失败，你首先会采取哪些技巧呢？这里是我们强烈推荐的首选方法。1 考虑使用更小的基组 由于一定的基组对应于一定精度和速度，所以更换基组并不在所有的情况下都适用。方法是首先用小基组进行计算，由前一个波函得到用于大基组计算的初始猜测(Guess=Read自动进行)。 2 增加最大循环步数 Gaussian默认的最大循环步数为64 (SCF=DM或SCF=QC方法则为512)，如果循环次数超过这个数目则会汇报convergence failure。在一定的情况下，不收敛的原因仅仅是因为最大循环步数不够。可以通过设置maxcyc来增大最大循环步数。更多的SCF迭代(SCF(MaxCycle=N)，其中N是迭代数)。这很少有帮助，但值得一试。 3 放宽收敛标准 如果接近SCF但未达到，收敛标准就会放松或者忽略收敛标准。这通常用于
不是在初始猜测而是在平衡结构收敛的几何优化。SCF=Sleazy放松收敛标准，Conver选项给出更多的控制。在Gaussian中设置收敛标准的命令行是Conv=8 表示10-8,你可以将这个标准降低，比如改为10-5等等。微调结构是比较有效的方法，能量收敛标准不宜太低，最好不低于10-7. 4 尝试改变初始构型 首先略微减小键长，接下来略微增加键长，接下来再对结构作一点改变。 5 尝试能级移动Level shifting (SCF=Vshift) 如果不收敛的原因是波函数的震荡行为，通常是因为在相近的能量上的泰的混合。对于这种情况，我们可以采用level shifting的方法。Level shifting的含义是人工的升高非占据轨道的能级，以防止和最高占据轨道之间的混合，以达到收敛的目的。在Gaussian中此方法的关键词为SCF=Vshift 6 使用强制的收敛方法SCF=QC SCF=QC通常最佳，但在极少数情况下SCF=DM更快。此关键字将大大增加计算时间，但是收敛的机会更大。不要忘记给计算额外增加一千个左右的迭代。应当测试这个方法获得的波函，保证它最小，并且正好不是稳定点(使用Stable 关键字)。 7 对开壳层体系，尝试收敛到同一分子的闭壳层离子，接下来用作开壳层计算的初始猜测。添加电子可以给出更合理的虚轨道，但是作为普遍的经验规则，阳离子比阴离子更容易收敛。选项Guess=Read定义初始猜测从Gaussian计算生成的checkpoint文件中读取。 8 一些程序通过减小积分精度加速SCF。对于使用弥散函数，长程作用或者低能量激发态的体
大气污染扩散 第一节大气结构与气象 有效地防止大气污染的途径，除了采用除尘及废气净化装置等各种工程技术手段外，还需充分利用大气的湍流混合作用对污染物的扩散稀释能力，即大气的自净能力。污染物从污染源排放到大气中的扩散过程及其危害程度，主要决定于气象因素，此外还与污染物的特征和排放特性，以及排放区的地形地貌状况有关。下面简要介绍大气结构以及气象条件的一些基本概念。 一、大气的结构 气象学中的大气是指地球引力作用下包围地球的空气层，其最外层的界限难以确定。通常把自地面至1200 km左右范围内的空气层称做大气圈或大气层，而空气总质量的98．2％集中在距离地球表面30 km以下。超过1200 km的范围，由于空气极其稀薄，一般视为宇宙空间。 自然状态的大气由多种气体的混合物、水蒸气和悬浮微粒组成。其中，纯净干空气中的氧气、氮气和氩气三种主要成分的总和占空气体积的99.97％，它们之间的比例从地面直到90km高空基本不变，为大气的恒定的组分；二氧化碳由于燃料燃烧和动物的呼吸，陆地的含量比海上多，臭氧主要集中在55～60km高空，水蒸气含量在4％以下，在极地或沙漠区的体积分数接近于零，这些为大气的可变的组分；而来源于人类社会生产和火山爆发、森林火灾、海啸、地震等暂时性的灾害排放的煤烟、粉尘、氯化氢、硫化氢、硫氧化物、氮氧化物、碳氧化物为大气的不定的组分。 大气的结构是指垂直（即竖直）方向上大气的 密度、温度及其组成的分布状况。根据大气温度在 垂直方向上的分布规律，可将大气划分为四层：对 流层、平流层、中间层和暖层，如图5－1所示。 1. 对流层 对流层是大气圈最靠近地面的一层，集中了大 气质量的75％和几乎全部的水蒸气、微尘杂质。受 太阳辐射与大气环流的影响，对流层中空气的湍流 运动和垂直方向混合比较强烈，主要的天气现象云 雨风雪等都发生在这一层，有可能形成污染物易于 扩散的气象条件，也可能生成对环境产生有危害的 逆温气象条件。因此，该层对大气污染物的扩散、输送和转化影响最大。 大气对流层的厚度不恒定，随地球纬度增高而降低，且与季节的变化有关，赤道附近约
关于金属汞扩散的问题 引言： 我们都知道，重金属丢弃到土地后会严重污染环境，同时对人体健康造成危害。著名的秦始皇陵墓，据专家在陵墓周围取数据观测，周围的汞含量呈现出外渗的趋势。也就是说，随着外围半径的扩大，汞含量浓度递减，并且随着时间的增加，汞渗透的半径越来越大。这就证明了汞金属在泥土中会发生扩散。因此，我们就提出，能否通过在外部取样的观察数据，建立一个数学模型，来判断陵墓中心处汞的浓度呢？ 模型的提出： 由于汞的扩散快慢跟本身的化学性质,物理性质有关。还有，由于在土堆里头，在各个方向上受到的力不相同和各种因素的影响，因此扩散的速度也会有差异。例如东西方向和南北方向会因为地球的自传而扩散速度会不一样。另一方面，汞在扩散的过程，由于泥土的吸收，化学反应等因数的影响，也会影响到汞的扩散。 为此我们引入一个函数u(x, y, z, t)，它表示t时刻在（x,y,z）处汞的浓度。我们的目标就是利用所观测到的数据，来推断出这个函数的表达式。 模型符号的引入： 为了表示汞在想x,y,z 方向上的扩散速度，我们在此引
入扩散系数： 2 a ：x 方向上的扩散系数 2 b ：y 方向上的扩散系数 2 c ：z 方向上的扩散系数 2 k ：由于泥土吸收，化学反应而引起的衰减系数 M ：扩散源汞的质量 模型假设： 1。假设有一汞扩散源，汞从扩散源沿 x ，y ，z 三个方向向四周扩散。 2。扩散前周围空间此物质的浓度为零。 3。扩散过程中没有人为因素的影响。 模型建立： u(x, y, z, t) 是 t 时刻点 (x, y , z) 处某物质的浓度。任取一个闭曲面 S ，它所围的区域是Ω，由于扩散，从 t 到 t t +? 时刻这段时间内，通过 S 流入Ω的质量为 1 M 2 2 2 1(cos cos cos )d d t t t S u u u M a b c S t x y z αβγ+????= ++???? ?? 其中 2 a ,2 b ,2 c 分别是沿 x ，y ，z 方向的扩散系数。 由高斯公式 ： ? ??? ?+Ω ??+??+??= t t t t z y x z u c y u b x u a M d d d d )(2 2 2 2 2 2 2 2 2 1
Abaqus遇到的问题小结 1.Abaqus量纲系统（SI）：m、N、Kg、s、Pa、J、Kg/m3。 2.有限元单元类型（Element Type）：实体单元、壳单元、杆件单元（梁Beam、桁架Truss） 怎样将部件的不同部分设定为不同的子集 tools >partition cell/plane（如果是要分层使用，则进行分割，分别赋予不同的属性;或者对其中的单元进行设置Set；如果要在运算中替换则需要在关键词里添加语句定义） 4.建模的过程中输入的尺寸被自动降低精度（四舍五入），请问如何设置尺寸精度（小数点位数） 在 sketch option里打开sketch模块(左下角最下面)打开sketchoption 对话框,其中decimal places 表示小数点位数,默认为2,最高能调到6 。 中出现，Dependent partinstances connot be edited ，怎么办 mesh on part , mesh on instance一个非独立实体只是原始部件的一个指针,可以对原始部件划分网格，但是不能对一个非独立实体划分网格，即mesh on part。 方法：左边模型树里面，Assembly->instances里面的子选项右击，点Make independent;或者点击model旁边的Part,对part进行网格划分。 6.材料属性中拉伸强度和压缩强度怎么定义 计算结果有每个计算节点的应力，计算中用不到强度值（弹性力学三组基本方程不涉及强度，求解过程用不到）。（强度值是用来校核的：断裂，屈服或超过线性变形） 7.如何选择创建独立实体还是非独立实体 如果集合中包含许多具有相同性质的部件，则创建非独立实体（Dependent）。反之，如果集合中包含许多不同性质的部件，则创建独立实体（Independent）要有优势些。 8.地基土的分层问题 方法1：定义为一个part，然后partition分层赋材料属性，不需要设置接触。（超级赞！！！）
A list of error messages and possible solutions - Gaussian calculations can fail with various error messages. Some error messages from .out and .log files - and possible solutions - have been compiled here to facilitate problem solving.- These are divided into:- Syntax and similar errors- 语法类错误 Memory and similar errors- 内存类错误 Convergence problems - 不收敛错误 Errors in solvent calculations - 溶剂中的计算错误 Errors in log files- 错误文件 - ERROR MESSAGES IN OUTPUT FILES- Syntax and similar errors:End of file in ZSymb.- Error termination via Lnk1e in /global/apps/gaussian/g03.e01/g03/l101.exe Solution: The blank line after the coordinate section in the .inp file is missing. (输入文件空行丢失) Unrecognized layer "X".-（不识别层X） Error termination via Lnk1e in /global/apps/gaussian/g03.e01/g03/l101.exeSolution: Error due to syntax error(s) in coordinate section (check carefully). If error is "^M", it is caused by DOS end-of-line characters (e.g. if coordinates were written under Windows). Remove ^M from line ends using e.g. emacs. To process .inp files from command line, use sed -i 's/^M//' File.inp (Important: command does not work if ^M is written as characters - generate ^M on command line using ctrl-V ctrl-M).-QPERR --- A SYNTAX ERROR WAS DETECTED IN THE INPUT LINE.- Solution: Check .inp carefully for syntax errors in keywords - RdChkP: Unable to locate IRWF=0 Number= 522.- Error termination via Lnk1e in /global/apps/gaussian/g03.e01/g03/l401.exe or-FileIO operation on non-existent file.- [...] Error termination in NtrErr:- NtrErr Called from FileIO.Solution: Operation on .chk file was specified (e.g. geom=check, opt=restart), but .chk was not found. Check that:-
Initialization Files The Gaussian system includes initialization files to set up the user environment for running the program. These files are $g03root/g03/bsd/g03.login C shell $g03root/g03/bsd/g03.profile Bourne shell Note that the g03root environment variable must be set up by the user. Thus, it is customary to include lines like the following within the .login or .profile file for Gaussian users: .login files: setenv g03root location source $g03root/g03/bsd/g03.login .profile files: g03root=location export g03root . $g03root/g03/bsd/g03.profile Once things are set up correctly, the g03 command is used to execute Gaussian 03 (see below). # .bash_profile # Get the aliases and functions if [ -f ~/.bashrc ]; then . ~/.bashrc fi # User specific environment and startup programs PATH=$PATH:$HOME/bin:/home/g03: export USERNAME BASH_ENV PATH g03root=/home/ export g03root . $g03root/g03/bsd/g03.profile GAUSS_SCRDIR=/home/tmp export GAUSS_SCRDIR LD_LIBRARY_PATH=/home/g03 GAUSS_EXEDIR=/home/g03 ----------------------------- 运行: g03 try ,结果如下: ----------------------------- Entering Gaussian System, Link 0=g03 Input=https://www.doczj.com/doc/272790861.html, Output=try.log Initial command:
初学Gaussian03常见出错分析 最初级错误 1. 自旋多重度错误 2. 变量赋值为整数 3. 变量没有赋值 4. 键角小于等于0度，大于等于180度 5. 分子描述后面没有空行 6. 二面角判断错误，造成两个原子距离过近 7. 分子描述一行内两次参考同一原子,或参考原子共线 运行出错 1. 自洽场不收敛 SCF a. 修改坐标，使之合理 b. 改变初始猜 Guess c. 增加叠代次数 SCFCYC=N d. iop（5/13=1） 2. 分子对称性改变 a. 修改坐标，强制高对称性或放松对称性 b. 给出精确的、对称性确定的角度和二面角 c. 放松对称性判据 Symm=loose d. 不做对称性检查 iop(2/16=1) 3. 无法写大的Scratch文件RWF a. 劈裂RWF文件 %rwf=loc1,size1,loc2,size2,……..,locN,-1 b. 改变计算方法 MP2=Direct可以少占硬盘空间 c. 限制最大硬盘 maxdisk=N GB 4. FOPT出错原因是变量数与分子自由度数不相等。可用POPT 或直接用OPT 5. 优化过渡态只能做一个STEP 原因是负本征数目不对添加 iop（1/11）=1 6. 组态相互作用计算中相关能叠代次数不够，增加叠代次数 QCISD（Maxcyc=N） Default.Rou设置 ?在Scratch文件夹中的Default.Rou文件中设置G03程序运行的省缺参数： ? -M- 200MW ? -P- 4 ? -#- MaxDisk=10GB ? -#- SCF=Conventional or Direct ? -#- MP2=NoDirect or Direct ? -#- OPTCYC=200 ? -#- SCFCYC=200 ? -#- IOPs 设置如iop（2/16=1） Default.Rou设置中的冲突 ? Default route: MaxDisk=2GB SCF=Direct MP2=Direct OPTCYC=200 SCFcyc=100 iop(2/16=1) iop(5/13=1) ? ------------------
1 云团扩散模型 根据物质泄漏后所形成的气云的物理性质的不同，可以将描述气云扩散的模型分为非重气云模型和重气云模型两种[5-13]。 1.1 非重气云模型 高斯模型是一种常用的非重气扩散模型，高斯烟羽(Plume model)模型又称高架点连续点源扩散模型，适用于连续源的扩散，即连续源或泄放时间大于或等于扩散时间的扩散。 高斯烟团(Puff model)模型适用于短时间泄漏的扩散，即泄放时间相对于扩散时间比较短的情形，如突发性泄放等。若假设气体云内空间上的分布为高斯分布，则地面地处风向的烟团浓度分布算式为 式中， c(x，y，H)——点(x，y，H)处浓度值，mg/m3； Q——源强，即单位时问的排放量，mg/s； u——环境平均风速，m/s； σx，σy，σz——扩散参数； H——源高(烟团高度)，m； x——下方向到泄漏原点的距离，m； y，z——侧风方向、垂直向上方向离泄漏原点的距离，m。 高斯模式的实际应用效果很大程度上依赖于如何给定模式中的一些参数，尤其要注意源强、扩散参数等的确定。 源强与污染物的物理化学属性、扩散方式、释放点的地理环境等有关。扩散参数表征大气边界层内
湍流扩散的强弱，是高斯模式的一项重要数据。高斯扩散模式所描述的扩散过程(实质上也包含了在实际应用中对高斯模式的一些限制)主要有： 1)下垫面平坦、开阔、性质均匀，平均流场稳定，不考虑风场的切变。 2)扩散过程中，污染物本身是被动、保守的，即污染物和空气无相对运动，且扩散过程中污染物无损失、无转化，污染物在地面被反射。 3)扩散在同一温度层结中发生，平均风速大于1.0 m/s。 4)适用范围一般小于10～20 km。 1.2 重气云模型 由于重气本身的特殊性，在重气扩散领域也有大量基于不同理论的模型。鉴于重气扩散与中性或浮性气体扩散有着明显的区别，目前国内外已开发大量的不同复杂程度的重气扩散模型，如箱模型、相似模型、LTA－HGDM模型、CFD模型等。 1.2.1 箱(BOX)模型 箱模型是指假定浓度、温度和其他场，在任何下风横截面处为矩形分布等简单形状，这里的矩形分布是指在某些空间范围内场是均匀的，而在其他地方为零。该类模型预报气云的总体特征，如平均半径、平均高度和平均气云温度，而不考虑其在空间上的细节特征。重气效应消失后其行为表现为被动气体扩散，所以该类模型还包括被动扩散的高斯模型及对它的修正。 1.2.2 层流及湍流大气环境中的重气扩散(LTA－HGDM)模型 LTA－HGDM模型(Heavy Gas Dispersion Model in Lsaminar and Turbulent Atmosphere层流及湍流大气环境中的重气扩散模型)以箱模型为基础，结合虚点源模型，能描述重气泄漏扩散整个过程。模型同三维有限元模型相比，具有形式简单、原始输入数据运算速度快等优点。 LTA－HGDM模型的建立基于以下几点假设： 1)危险性气体初时泄漏时，其外形呈正圆柱形(H=2R)。 2)初始时刻泄漏源即此核电站内部的浓度、温度呈均匀分布。 3)扩散过程不考虑泄漏源即此核电站内部温度的变化，忽略热传递、热对流及热辐射。
【分享】高斯对不收敛问题的对策 首先，我们必须理解收敛是什么意思。在自洽场(SCF)计算中，自洽循环中，首先产生一个轨道占据的初始猜测， 1)然后根据此轨道占据构造电荷密度和哈密顿量。 2)对角化哈密顿量，得到新的轨道能级和占据。 3)产生新的电荷分布和哈密顿量，重复步骤2） 经过一定次数的循环后，某次循环前和循环后的电荷密度差别小于一定的标准，我们称之为收敛。 如果以上过程不能收敛，则gaussian给出convergence failure的警告。 如果SCF计算收敛失败，你首先会采取哪些技巧呢？这里是我们强烈推荐的首选方法。 1 考虑使用更小的基组 由于一定的基组对应于一定精度和速度，所以更换基组并不在所有的情况下都适用。方法是首先用小基组进行计算，由前一个波函得到用于大基组计算的初始猜测(Guess=Read自动进行)。 2 增加最大循环步数 Gaussian默认的最大循环步数为64 (SCF=DM或SCF=QC方法则为512)，如果循环次数超过这个数目则会汇报convergence failure。在一定的情况下，不收敛的原因仅仅是因为最大循环步数不够。可以通过设置maxcyc来增大最大循环步数。更多的SCF迭代(SCF(MaxCycle=N)，其中N是迭代数)。这很少有帮助，但值得一试。 3 放宽收敛标准 如果接近SCF但未达到，收敛标准就会放松或者忽略收敛标准。这通常用于不是在初始猜测而是在平衡结构收敛的几何优化。SCF=Sleazy放松收敛标准，Conver选项给出更多的控制。 4 尝试改变初始构型 首先略微减小键长，接下来略微增加键长，接下来再对结构作一点改变。 5 尝试能级移动Level shifting (SCF=Vshift) 如果不收敛的原因是波函数的震荡行为，通常是因为在相近的能量上的泰的混合。对于这种情况，我们可以采用level shifting的方法。Level shifting的含义是人工的升高非占据轨道的能级，以防止和最高占据轨道之间的混合，以达到收敛的目的。在Gaussian中此方法的关键词为SCF=Vshift 6 使用强制的收敛方法SCF=QC SCF=QC通常最佳，但在极少数情况下SCF=DM更快。此关键字将大大增加计算时间，但是收敛的机会更大。不要忘记给计算额外增加一千个左右的迭代。应当测试这个方法获得的波函，保证它最小，并且正好不是稳定点(使用Stable关键字)。 7 对开壳层体系，尝试收敛到同一分子的闭壳层离子，接下来用作开壳层计算的初始猜测。添加电子可以给出更合理的虚轨道，但是作为普遍的经验规则，阳离子比阴离子更容易收敛。选项Guess=Read定义初始猜测从Gaussian计算生成的checkpoint文件中读取。 8 一些程序通过减小积分精度加速SCF。对于使用弥散函数，长程作用或者低能量激发态的体系，必须使用高积分精度：SCF=NoVarAcc。 9 改变模型或方法 可以考虑改变模型方法。比较常见的方法有HF,GVB,MCSCF,CASSCF,MPn等。改变模型方法通常也会收敛性质。通常，精度更高的方法更难收敛。精度比较低的方法产生的计算结果可以作为高精度计算的初始猜测。考虑使用不同理论级别的计算。这并不总是实用的，但除此之外，增加迭代数量总是使得计算时间和使用更高理论级别差不多。 10 关闭DIIS外推(SCF=NoDIIS)。同时进行更多的迭代(SCF=(MaxCycle=N))。 11 试着改用DIIS之外其它方法(SCF=SD或SCF=SSD)。
赛区评阅编号（由赛区组委会填写）： 2015高教社杯全国大学生数学建模竞赛 承诺书 我们仔细阅读了《全国大学生数学建模竞赛章程》和《全国大学生数学建模竞赛参赛规则》（以下简称为“竞赛章程和参赛规则”，可从全国大学生数学建模竞赛网站下载）。 我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。 我们知道，抄袭别人的成果是违反竞赛章程和参赛规则的，如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。 我们郑重承诺，严格遵守竞赛章程和参赛规则，以保证竞赛的公正、公平性。如有违反竞赛章程和参赛规则的行为，我们将受到严肃处理。 我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。 我们参赛选择的题号（从A/B/C/D中选择一项填写）： C 我们的报名参赛队号（12位数字全国统一编号）： 参赛学校（完整的学校全称，不含院系名）：温州医科大学 参赛队员 (打印并签名) ：1. 章成俊 2. 杨超 3. 谢锦 指导教师或指导教师组负责人 (打印并签名)： 日期：年月日
赛区评阅编号（由赛区组委会填写）： 2015高教社杯全国大学生数学建模竞赛 编号专用页 送全国评阅统一编号（由赛区组委会填写）： 全国评阅随机编号（由全国组委会填写）：
对垃圾处理厂污染的动态监控及居民补偿 摘要 城市垃圾处理问题是一个世界性难题。目前垃圾焚烧正逐步成为中国垃圾处理的主要手段之一。本论文构根据题目设置的垃圾处理厂规模，建立了环境动态监控体系，并根据潜在污染风险对周围居民进行了合理经济补偿的设计。 对于问题（1），为了实现对垃圾焚烧厂烟气排放及相关环境影响状况的动态监控，本论文在高斯烟羽模型的基础上进行改进，引入温度、降雨对污染物扩散的影响，建立了新的污染物扩散模型。本论文创新性的提出了风雨影响指数M,用来衡量风向、降雨对颗粒物扩散的影响。本论文将抽象的污染物含量形象化，利用空气污染指数API描述具体的污染程度及其给周围居民带来的影响。并且从不同角度给出了模型检验，验证了所建模型的准确性。 对于问题（1）具体赔偿方案的制定，在综合考虑了不同方位风向频率、受污染时间、受污染程度的基础上，本论文使用了层次分析法，并且进行了一致性检验，使得赔偿方案具有说服力。通过MATLAB编程，计算出当政府和垃圾处理厂共支付风险赔偿金为N时，得出居住地的每位居民应得的赔偿金额计算公式。对于监测点的设置，经计算共需21个，具体布置情况见后文。 对于问题（2），在题目所述的发生事故的情况下，对污染物的具体含量进行了合理的预测与假设。模拟出酸性物质与颗粒物的影响范围，并根据具体的污染程度设置不同的污染区。对每个污染区的不同情况设置更改监测点的设置，并且在问题（1）的基础上对居民的经济补偿进行合理修改。 关键词：高斯烟羽模型，层次分析法，空气污染指数，烟气抬升公式 一、问题重述 “垃圾围城”是世界性难题，在今天的中国显得尤为突出。数据显示，目前全国三分之二以上的城市面临“垃圾围城”问题，垃圾堆放累计侵占土地75万亩。因此，垃圾焚烧正逐步成为中国垃圾处理的主要手段之一。然而，由于政府监管不力、投资者目光短浅等多方面的原因，致使前些年各地建设的垃圾焚烧电厂在运营中出现了环境污染问题，给垃圾焚烧技术在我国的推广造成了很大阻力，许多城市的新建垃圾焚烧厂选址都出现因居民反对而难以落地的局面。在垃圾焚烧厂运行监管方面，目前主要是在垃圾焚烧厂内进行测量监控，缺少从周边环境视角出发的外围动态监控，因而难以形成为民众所信服的全方位垃圾焚烧厂环境监控体系。 深圳市某地点计划建立一个中型的垃圾焚烧厂，计划处理垃圾量1950吨/天（设置三台可处理垃圾650吨/天的焚烧炉，排烟口高度80米，每天24小时运转）。从构建环境动态监控体系、并根据潜在污染风险对周围居民进行合理经济补偿的需求出发，有关部门希望能综合考虑垃圾焚烧厂对周围带来环境污染以及其他危害的多种因素（例如，焚烧炉的污染物排放量、居住点离开垃圾焚烧厂的距离、风力和风向及降雨等气象条件、地形地貌以及建筑物的遮挡程度等等），在进行科学定量分析的基础
优化 1在优化时采用Scf(tight)的选项，增加收敛的标准。再去计算频率。如果还有虚频，参见下一步。 2.对称性的影响，很多情况下的虚频是由于分子本身的对称性造成的。这样， 在优化时，如果必要，要将对称性降低，还有，输入文件有时是用内坐标。 建议如果有虚频的话，将内坐标改成直角坐标优化。 3.如果上述方法还有虚频，看一下虚频，找到强度较大的，将在频率中 产生的原子的振动坐标加到相应的输入文件中。这样，重新计算。直到 虚频没有。 4实际上，如果分子柔性较大，很难找到最低点，这是电子结构计算的问题， 这种情况下，需要动力学的东西，用构象搜寻的办法解决。 如：模拟退火，最陡下降法，淬火法等。将得到的能量最低的构象做 一般的电子结构计算，这样，应当没有问题。 不要讲你还没有得到最稳定的结构，那么，是你的分子有问题， 要么计算错了，要么就是游离出现代计算的范畴 1。低频率振动模式时，力常数很小，必须使用opt＝tight以确保适当的收敛和随后任务步骤中频率计算的可靠度。 比如我曾经算过一个体系，第一个频率振动的力常数只有10～20，很不放心，加入tight后，优出来的结构却与前面大相径庭，而且还出现了一个虚频。
2。做IRC计算的时候，如果步长很小，必须要用verytight。就曾有人步长设为1，而不加verytight，算出来的反应路径是v字形倒海鸥状的。 3。使用int指定用于数值积分的积分网格来用于消虚频的时候，必须注意，在比较能量的时候，对所有的计算要使用相同的积分网格。 检查错误 1．检查是否有初始文件错误 在命令行中加入%kJob L301 or %kJob L302 如果通过则一般初始文件ok。 常见初级错误： a.自旋多重度错误 b.变量赋值为整数 c.变量没有赋值或多重赋值 d.键角小于等于0度，大于等于180度 e.分子描述后面没有空行 f.二面角判断错误，造成两个原子距离过近 g.分子描述一行内两次参考同一原子,或参考原子共线 2．SCF(自洽场)不收敛 则一般是L502错误 省却情况做64个cycle迭代(G03缺省128 cycles) a.修改坐标，使之合理 b.改变初始猜Guess=Huckel或其他的，看Guess关键词。
放射性气体扩散浓度预估模型 【摘要】本文是以日本地震引起的福岛核电站的核泄漏为背景，并以给出的数据为基础，研究某一假设核电站的核泄漏问题。我们通过收集相关的资料，并结合题目给出的数据，建立了高斯模型、连续点源高斯扩散模型解决了题目提出的四个问题。 针对问题一：考虑到泄漏源是连续、均匀和稳定的，我们运用散度、梯度、流量等数学概念，通过“泄漏放射性物质质量守恒”、“气体泄漏连续性定理”、 Guass 公式及积分中值定理得到了无界区域的抛物线型偏微分方程，然后再通过电源函数解出空间任意一点的放射性物质浓度的表达式，把此表达式定为模型一的前身。鉴于放射性物质的扩散受到诸多因素的影响，如：泄漏源的实际高度、地面反射等。我们以泄漏口为坐标原点建立三维坐标系，通过“像源法”处理地面反射对放射性物质浓度的影响，并由此对模型一的前身进行修正完善，得到模型一：高斯模型，即放射性物质浓度的预测模型。最后我们模拟了放射性物质无风扩散仿真图。 针对问题二：当风速为k m/s 时，我们根据放射性核素云团在大气中迁移和扩散的数值计算的基本方法和步骤，并以泄漏点源在地面的投影点为坐标原点，以风向方向为x 轴，铅直方向为z 轴，与x 轴水平面垂直方向为y 轴建立三维坐标系，地面的反射作用同样利用“像源法”进行处理，得到连续点源高斯扩散模型。考虑到地面反射、烟云抬升、放射性物质自身的沉降及雨水的吸附等对浓度的影响，我们对连续点源高斯扩散模型进行了修正，建立了修正的连续点源高斯扩散模型。最后利用大气稳定度确定了扩散参数，进而求解了模型。 针对问题三：经分析，问题三的提出是以问题二为基础的，模型三的建立只需要将模型二加以调整即可。我们以风速方向为x 轴正方向，将风速与放射性物质的扩散速度进行矢量运算，此问题则转化为求(,0,)L z 和(,0,)L z -两点处的放射性物质浓度，由此建立模型三，即上风和下风L 公里处放射性物质浓度浓度的预测模型。 针对问题四：首先，我们通过网络收集了相关数据，然后，我们结合模型二、模型三对数据进行整理代入，算出了日本福岛核电站泄漏的放射性物质扩散到中国东海岸和美国西海岸的浓度分别为334.242910/g m -?、432.385410/g m -?。 关键词：高斯模型 连续点源高斯扩散模型 核泄漏
如何从下面的Gaussian输出文件中找出轨道系数及轨道能!!(新手多谢),请帮忙 标出来!!! 求助]如何从下面的Gaussian输出文件中找出轨道系数及轨道能!!(新手多谢),请帮忙标出来!!! The electronic state is 1-A1. Alpha occ. eigenvalues -- -20.58265 -11.33946 -1.39265 -0.87259 -0.69715 Alpha occ. eigenvalues -- -0.63950 -0.52294 -0.44073 Alpha virt. eigenvalues -- 0.13573 0.24842 0.33338 0.37329 0.73660 Alpha virt. eigenvalues -- 0.80783 0.84685 0.94689 1.10445 1.10700 Alpha virt. eigenvalues -- 1.13937 1.27145 1.33529 1.62050 1.78192 Alpha virt. eigenvalues -- 1.79416 1.99239 2.18347 2.23684 2.45514 Alpha virt. eigenvalues -- 2.64513 2.87165 2.97616 3.27576 4.09792 Alpha virt. eigenvalues -- 4.47637 Molecular Orbital Coefficients 1 2 3 4 5 (A1)--O (A1)--O (A1)--O (A1)--O (B2)--O EIGENVALUES -- -20.58265 -11.33946 -1.39265 -0.87259 -0.69715 1 1 C 1S 0.00000 0.99566 -0.11060 -0.1626 2 0.00000 2 2S 0.00047 0.02675 0.20981 0.33995 0.00000 3 2PX 0.00000 0.00000 0.00000 0.00000 0.00000 4 2PY 0.00000 0.00000 0.00000 0.00000 0.42017 5 2PZ -0.00007 0.0006 6 0.17259 -0.18451 0.00000 6 3S -0.00024 -0.00743 0.08051 0.31309 0.00000 7 3PX 0.00000 0.00000 0.00000 0.00000 0.00000 8 3PY 0.00000 0.00000 0.00000 0.00000 0.15760 9 3PZ -0.00048 0.00135 -0.01160 -0.07970 0.00000 10 4XX -0.00002 -0.00272 -0.01628 -0.01333 0.00000 11 4YY -0.00006 -0.00202 -0.01365 0.03019 0.00000 12 4ZZ -0.00074 -0.00123 0.03302 -0.00166 0.00000 13 4XY 0.00000 0.00000 0.00000 0.00000 0.00000 14 4XZ 0.00000 0.00000 0.00000 0.00000 0.00000 15 4YZ 0.00000 0.00000 0.00000 0.00000 -0.01394 16 2 O 1S 0.99472 -0.00038 -0.19672 0.08889 0.00000 17 2S 0.02094 0.00025 0.44184 -0.20351 0.00000 18 2PX 0.00000 0.00000 0.00000 0.00000 0.00000 19 2PY 0.00000 0.00000 0.00000 0.00000 0.32122 20 2PZ -0.00153 -0.00029 -0.13537 -0.14216 0.00000 21 3S 0.00436 -0.00058 0.37895 -0.27048 0.00000 22 3PX 0.00000 0.00000 0.00000 0.00000 0.00000 23 3PY 0.00000 0.00000 0.00000 0.00000 0.17976
高斯扩散模型及其适用条件 （1）一般表达式 根据质量守恒原理和梯度输送理论，污染物在大气中一般运动规律为：（3分） 1N x y z p p c c c c c c c u v w k k k S t x y z x x y y z z =????????????????+++=+++ ? ? ?????????????????∑ C ：污染物质平均浓度； X ，y ，z ：三个方向坐标； u ，v ，w ：三个方向速度分量； k x ，k y ，k z ：三个方向扩散系数； t ：为污染物扩散时间； S P ：污染物源、汇强度。 （2）高斯模型的适用条件：①大气流动稳定，表明污染物浓度不随时间改变，即0t ?=?； ②有主导风向，表明u=常数，且v=w=0； ③污染物在大气中只有物理运动，物化学 和生物变化，且预测范围内无其他同类 污染的源和汇。表明S P =0（p=1,2,….n ） 此时三维的动态模型就可简化为三维的稳态模型，得： x y z c c c c u k k k x x x y y z z ?????????????=++ ? ? ???????????? ?? （3分） ④有主导风情况下，主导风对污染物输送 应远远大于湍流运动引起污染物在主导风方
向上扩散。即c u x ??（平流输送作用）远远大于x c k x x ???? ????? （湍流弥散作用）。 此时方程又可以简化为： y z c c c u k k x y y z z ?????????=+ ? ???????? ?? （2分） （3）由于y 和z 方向上污染物浓度不发生变化，故规定y k 与y 无关，z k 与z 无关，即： 22z 22z y c c c u k k x y ???=+??? （1分） （4）由质量守恒原，理运用连续点源源强计算方式，按照单元体积(3)简化得到的方程进行积分ucdydz=Q ∞∞ -∞-∞??，结合边界条件 {0c=x y z c=0x y z ===∞ →∞时，，，时，对方程进行求解。（2分） （5）设x=ut ，令22y y z z =2k t =2k t σσ；。化简求解得到高斯扩散模型的标准 形式： ()2222y z 1c ,,exp 22y z Q y z x y z u πσσσσ????=-+?? ? ??????? （1分）