小木虫 --- 600万学术达人喜爱的学术科研平台
&&查看话题
求助高人指点，高斯错误代码L301是什么意思
这是我计算的结果的部分内容：
--------------------------------------------------------------------
Center& &&&Atomic& &&&Atomic& && && && &&&Coordinates (Angstroms)
Number& &&&Number& && &Type& && && && &&&X& && && &&&Y& && && &&&Z
---------------------------------------------------------------------
& & 1& && && & 6& && && && & 0& && & -0.143084& &-0.052209& &-0.040829
& & 2& && && & 7& && && && & 0& && & -0.980902& & 0.975376& & 0.033661
& & 3& && && & 6& && && && & 0& && & -0.452993& & 2.229130& & 0.053048
& & 4& && && & 6& && && && & 0& && &&&0.921885& & 2.429049& &-0.000254
& & 5& && && & 6& && && && & 0& && &&&1.730932& & 1.296900& &-0.072600
& & 6& && && & 7& && && && & 0& && &&&1.202940& & 0.044135& &-0.096391
& & 7& && && & 1& && && && & 0& && & -1.154130& & 3.052138& & 0.113047
& & 8& && && & 1& && && && & 0& && &&&1.346311& & 3.424889& & 0.015108
& & 9& && && & 1& && && && & 0& && &&&2.812981& & 1.334025& &-0.106934
& &10& && && & 1& && && && & 0& && &&&2.501331& &-1.172726& &-0.123540
& &11& && && & 8& && && && & 0& && &&&3.293581& &-1.798621& & 0.037992
& &12& && && & 1& && && && & 0& && &&&3.161934& &-2.011106& & 1.032014
& &13& && && & 1& && && && & 0& && & -2.748449& & 0.438515& & 0.103419
& &14& && && & 8& && && && & 0& && & -3.297237& &-0.421438& & 0.030811
& &15& && && & 1& && && && & 0& && & -4.242737& &-0.218085& & 0.348285
& &16& && && &16& && && && & 0& && & -0.834085& &-1.773827& &-0.070734
& &17& && && & 1& && && && & 0& && & -2.197332& &-1.259749& & 0.002834
---------------------------------------------------------------------
Rotational constants (GHZ):& && &1.4395154& && &0.9868850& && &0.5877872
Standard basis: 3-21G (6D, 7F)
There are& &101 symmetry adapted basis functions of A& &symmetry.
Integral buffers will be& & 262144 words long.
Raffenetti 2 integral format.
Two-electron integral symmetry is turned on.
& &101 basis functions,& &171 primitive gaussians,& &101 cartesian basis functions
& & 39 alpha electrons& && & 39 beta electrons
& && & nuclear repulsion energy& && & 513. Hartrees.
NAtoms=& &17 NActive=& &17 NUniq=& &17 SFac= 1.00D+00 NAtFMM=& &60 Big=F
------------------------------------------------------------------------------
United Atom Topological Model (UA0&&parameters set).
UA0: Hydrogen& &17 is unbound. Keep it explicit at all point on the
UA0: potential energy surface to get meaningful results.
Error termination via Lnk1e in C:\G03W\l301.exe at Tue Jul 12 18:09:29 2011.
Job cpu time:&&0 days&&0 hours&&6 minutes 26.0 seconds.
File lengths (MBytes):&&RWF=& & 203 Int=& && &0 D2E=& && &0 Chk=& && &8 Scr=& && &1
下面是我的输入文件：
%chk=0712-2watersolvent-2.chk
# opt mp2/3-21+g scrf=(iefpcm,solvent=water,admp)
Title Card Required
N& && && && && && &1& && && && &&&B1
C& && && && && && &2& && && && &&&B2& & 1& && && && &&&A1
C& && && && && && &3& && && && &&&B3& & 2& && && && &&&A2& & 1& && && && &&&D1
C& && && && && && &4& && && && &&&B4& & 3& && && && &&&A3& & 2& && && && &&&D2
N& && && && && && &5& && && && &&&B5& & 4& && && && &&&A4& & 3& && && && &&&D3
H& && && && && && &3& && && && &&&B6& & 2& && && && &&&A5& & 1& && && && &&&D4
H& && && && && && &4& && && && &&&B7& & 3& && && && &&&A6& & 2& && && && &&&D5
H& && && && && && &5& && && && &&&B8& & 4& && && && &&&A7& & 3& && && && &&&D6
H& && && && && && &6& && && && &&&B9& & 5& && && && &&&A8& & 4& && && && &&&D7
O& && && && && && &6& && && && & B10& & 5& && && && &&&A9& & 4& && && && &&&D8
H& && && && && &&&11& && && && & B11& & 6& && && && & A10& & 5& && && && &&&D9
H& && && && && && &2& && && && & B12& & 1& && && && & A11& & 6& && && && & D10
O& && && && && && &2& && && && & B13& & 1& && && && & A12& & 6& && && && & D11
H& && && && && &&&14& && && && & B14& & 2& && && && & A13& & 1& && && && & D12
S& && && && && && &1& && && && & B15& & 2& && && && & A14& & 3& && && && & D13
H& && && && && &&&16& && && && & B16& & 1& && && && & A15& & 2& && && && & D14
& &B1& && && && & 1.
& &B2& && && && & 1.
& &B3& && && && & 1.
& &B4& && && && & 1.
& &B5& && && && & 1.
& &B6& && && && & 1.
& &B7& && && && & 1.
& &B8& && && && & 1.
& &B9& && && && & 2.
& &B10& && && && &3.
& &B11& && && && &0.
& &B12& && && && &2.
& &B13& && && && &2.
& &B14& && && && &0.
& &B15& && && && &1.
& &B16& && && && &1.
& &A1& && && &&&115.
& &A2& && && &&&122.
& &A3& && && &&&116.
& &A4& && && &&&122.
& &A5& && && &&&116.
& &A6& && && &&&121.
& &A7& && && &&&121.
& &A8& && && && &97.
& &A9& && && && &96.
& &A10& && && & 102.
& &A11& && && & 122.
& &A12& && && & 108.
& &A13& && && & 140.
& &A14& && && & 116.
& &A15& && && & 109.
& &D1& && && && & 0.
& &D2& && && && & 0.
& &D3& && && && & 0.
& &D4& && && & -180.
& &D5& && && &&&180.
& &D6& && && &&&180.
& &D7& && && & -177.
& &D8& && && & -173.
& &D9& && && &&&119.
& &D10& && && &-178.
& &D11& && && &-179.
& &D12& && && & 172.
& &D13& && && & 180.
& &D14& && && && &0.
很可能是溶剂模型的问题，换下模型，或者修改下参数
我原来算过一个在苯里面的，同样问题，如上处理后，就正常收敛了
北京学而思教育科技有限公司 地址：北京市海淀区北三环甲18号中鼎大厦A座1层102室 电话：010-
浏览器进程
打开微信扫一扫
随时随地聊科研%chk=input.chk；%mem=600MB(使用600MB内存空间）；%nprocl=4（使用4个计算节点）；%nprocs=2（每个计算节点使用2颗CPU）；下面是一个输入文件；#TRHF/6-31G(d)Test；MyfirstGaussianjob:water；01；O-0.；H-0.；H0.441-
%chk=input.chk
%mem=600MB
(使用600MB内存空间）
（使用4个计算节点）
（每个计算节点使用2颗CPU）
下面是一个输入文件
#T RHF/6-31G(d) Test
My first Gaussian job: water single point energy
O -0.464 0.177 0.0
H -0.464 1.137 0.0
H 0.441 -0.143 0.0
第一行以#开头,是运行的说明行,#T表示指打印重要的输出部分,
#P表示打印更多的信息.后面的RHF表示限制性Hartree-Fock方法,这里要输入计算所选用的理论方法
6-31G(d)是计算所采用的基组,就是用什么样的函数组合来描述轨道
Test是指不记入Gaussian工作档案,对于单机版没有用处.
计算的理论,如HF(默认关键词,可以不写),B3PW91;
计算采用的基组,如6-31G, Lanl2DZ;
布局分析方法,如Pop=R
波函数自恰方法,如SCF=Tight.
Pop=Reg只在输出文件中打印出最高的5条HOMO轨道和最低的5条LOMU轨道,而采用 Pop=Full则打印出全部的分子轨道.
SCF设置是指波函数的收敛计算时的设定,一般不用写,SCF=Tight设置表示采用比一般方法较严格的收敛计算.
Gaussian 包含多种化学模型，比如
Gaussian 关键词
Hartree-Fock自恰场模型
Becke型3 参数密度泛函模型，采用Lee-Yang-Parr 泛函
二级Moller-Plesset
四级Moller-Plesset
一般采用开壳层的可能性是
存在奇数个电子，如自由基，一些离子
有多个单电子的体系
描述键的分裂过程
计算采用的基组，如6-31G，LANL2DZ
三种最常用的价层劈裂基组3-21G、6-31G 和6-311G
半经验方法
包括AM1，MINDO/3，PM3，常见的软件包有MOPAC，AMPAC，HyperChem，以及Gaussian。半经验方法采用了一些实验得来的参数，来帮助对薛定鄂方程的求解。
从头算，在解薛定鄂方程的过程中，只采用了几个物理常数，包括光速，电子和核的质量，普朗克常数，在求解薛定鄂方程的过程中采用一系列的数学近似，不同的近似也就导致了不同的方法。最经典的是Hartree-Fock 方法，缩写为HF。
几何优化的输入
#R RHF/6-31G(d) Opt Test
表明采用RHF 方法，6-31G(d)基组进行优化
几何优化和单点能计算都将原子理想化了，实际上原子一直处于振动状态。在平衡态，这些振动是规则的和可以预测的。
频率分析的计算要采用能量对原子位置的二阶导数。HF 方法，密度泛函方法(如B3LYP)，二阶Moller-Plesset 方法(MP2)和CASSCF 方法(CASSCF)都可以提供解析二阶导数
例中采用的是已经优化好的几何构型，输入格式
# RHF/6-31G(d) Freq Test
其他方法得到的频率同样存在系统误差，下面列出对于频率和零点能的矫正因子 方法
频率矫正因子
零点能矫正因子
HF/6-31G(d)
MP2(Full)/6-31G(d)
MP2(FC)/6-31G(d)
SVWN/6-31G(d)
BLYP/6-31G(d)
B3LYP/6-31G(d)
零点能(Zero Point Energy)和内能(Thermal Energy)
设置Freq=ReadIsotopes 关键词
下面是一个计算实例
Temperature 298.15 Kelvin。 Pressure 1.0000 Atm。
Zero-point correction= 0.029201
Thermal corection to Energy= 0.032054
Thermal correction to Enthalpy= 0.032999
Thermal correction to Gibbs Free Energy= 0.008244
Sum of electronic and zero-point Energies= -113.837130
Sum of electronic and thermal Energies= -113.834277
Sum of electronic and thermal Enthalpies= -113.833333
Sum of electronic and thermal Free Energies= -113.858087
后面四行的四个能量分别为E0，E，H，G。计算为
E0=E(elec) + ZPE
E=E0 + E(vib) + E(rol) + E(transl)
基组是体系内轨道的数学描述。大的基组由于对电子在空间上有小的限制而具有更大的精确性。
STO-3G 是最小基组
增大基组的第一个方法就是增加每个原子基函数的数量。分裂基组，比如3-21G 和6-31G，对于价键轨道都用两个函数来进行描述，
分裂基组允许轨道改变其大小，但不能改变形状。极化基组则取消了这样的限制，增加了角动量。比如在碳原子上增加d 轨道的成分，在过渡金属上增加f 轨道成分。有些在氢原子上增加p 轨道成分。一般的，常用的极化基组是6-31G(d)，这个基组来源与6-31G 基组，并在其基础上，对于重原子增加了d 轨道的成分。由于这个基组是中等大小的基组，在计算中很常用
弥散函数(Diffuse Functions)
高角动量基组
第三周期以后的原子的基组
选择合适的理论模型
使用半经验方法
电子相关和后SCF 方法
HF 方法可以解决很多问题，但是HF 方法也有缺陷。HF 方法在理论上没有考虑电子相关，自然电子相关对能量的贡献就没有考虑
密度泛函方法
密度泛函方法(Density Functional Theory Methods，DFT DFT 方法通过泛函来计算电子相关，其将电子能量分为几个部分，动能，电子-核相互作
用，库仑排斥，以及其余部分的交换相关项，最后一部分又根据密度泛函方法分解为交换和
原子化能的计算方法是三个原子的能量之和，减去二氧化碳能量和零点能之和。结果如下
高精度能量模型
原子化能是分子与组成分子的原子的能量差， 如对于PH2 ， 其原子化能为 E(P)+2(EH())-R(PH2)
电子亲和势
电子亲和势指体系增加一个电子后能量的变化，计算方法为中性分子和其阴离子的能量差。
离子化能指体系减少一个电子的能量的变化，计算方法为中性分子和其阳离子的能量差距
质子亲和能
质子亲和能为体系增加一个质子后的能量变化，计算方法为分子与在其基础上增加一个质子的体系的能量差距
三亿文库包含各类专业文献、应用写作文书、高等教育、生活休闲娱乐、专业论文、行业资料、72高斯输入文件的说明等内容。　
　《清华高斯讲义》《G03 使用说明》中; 译本 (ZORK 译) 和英文原本;Explorer ...所给的分 子结构(不管是手绘的还是来自晶体结构的)在转换为高斯输入的文件时,... 　对于本软件,当使用者希望计算某一确定的线性方程组时,现在提示的位 置输入方程...GPIO软件说明文档 10页 免费
高斯消元法类头文件及其... 暂无评价 5页 2下载... 　高斯软件分子内坐标的输入_计算机软件及应用_IT/...下面以甲烷为例说明: 首先使用 chem3D 画一个甲烷...使用记事本打开保存的文件,如果是内坐标表示完成;... 　在线互动式文档分享平台,在这里,您可以和千万网友分享自己手中的文档,全文阅读其他用户的文档,同时,也可以利用分享文档获取的积分下载文档 　高斯常见问题_军事/政治_人文社科_专业资料。一、混合...(1/8)=2, 此语句说明在找过度态的时候, 以 2A...四、奥奥的输入文件 %mem=8mw %rwf=1,2000mb,... 　在线互动式文档分享平台,在这里,您可以和千万网友分享自己手中的文档,全文阅读其他用户的文档,同时,也可以利用分享文档获取的积分下载文档 　shared processors 用于指定多核处理器,本实验 室计算工作站可以选择 4,高斯默认...文件,点击 Results→vibration,弹出如下对话框,检查有无虚频: 无虚频则说明得到... 　在线互动式文档分享平台,在这里,您可以和千万网友分享自己手中的文档,全文阅读其他用户的文档,同时,也可以利用分享文档获取的积分下载文档 　8000系列高斯计中文使用说明书_自然科学_专业资料。8000...对输入电缆发送供电请求之后, standby 按钮上的 LED...的描述和操 作指令请咨询你的应用程序的帮助文件。...小木虫 --- 600万学术达人喜爱的学术科研平台
&&查看话题
请问如何用Gaussain计算出化合物的分子电子密度体积V
各位虫友们大家好。
& & 我学习量化计算没多久，碰到问题的时候回先自己搜索文献帖子等等看看有没有解决方法。但还是有很多问题目前是解决不了。特此发帖前来求助。
& & 目前我需要计算化合物的& &分子电子密度体积V 。
& & 想先用简单的化合物例如 甲烷分子 计算其分子电子密度体积，文献中是这样说的：
我在Gaussian中进行了生成了如下输入文件：
%chk=SCRFch4
# hf/3-21g SCRF=(Solvent=water)
C& && && && && &&&-0.& & 1.& & 0.
H& && && && && &&&-0.& & 0.& & 0.
H& && && && && &&&-0.& & 2.& & 0.
H& && && && && &&&-0.& & 2.& &-0.
H& && && && && &&&-1.& & 1.& & 0.
1 2 1.0 3 1.0 4 1.0 5 1.0
得到了如下的计算内容：
Entering Gaussian System, Link 0=g09
Output=AAAAAAA.log
Initial command:
/root/gaussian/g09/l1.exe &/tmp/Gau-9611.inp& -scrdir=&/tmp/&
Entering Link 1 = /root/gaussian/g09/l1.exe PID=& && &9612.
Copyright (c) 92,98,13,
& && && && &Gaussian, Inc.&&All Rights Reserved.
This is part of the Gaussian(R) 09 program.&&It is based on
the Gaussian(R) 03 system (copyright 2003, Gaussian, Inc.),
the Gaussian(R) 98 system (copyright 1998, Gaussian, Inc.),
the Gaussian(R) 94 system (copyright 1995, Gaussian, Inc.),
the Gaussian 92(TM) system (copyright 1992, Gaussian, Inc.),
the Gaussian 90(TM) system (copyright 1990, Gaussian, Inc.),
the Gaussian 88(TM) system (copyright 1988, Gaussian, Inc.),
the Gaussian 86(TM) system (copyright 1986, Carnegie Mellon
University), and the Gaussian 82(TM) system (copyright 1983,
Carnegie Mellon University). Gaussian is a federally registered
trademark of Gaussian, Inc.
This software contains proprietary and confidential information,
including trade secrets, belonging to Gaussian, Inc.
This software is provided under written license and may be
used, copied, transmitted, or stored only in accord with that
written license.
The following legend is applicable only to US Government
contracts under FAR:
& && && && && && &&&RESTRICTED RIGHTS LEGEND
Use, reproduction and disclosure by the US Government is
subject to restrictions as set forth in subparagraphs (a)
and (c) of the Commercial Computer Software - Restricted
Rights clause in FAR 52.227-19.
Gaussian, Inc.
340 Quinnipiac St., Bldg. 40, Wallingford CT 06492
---------------------------------------------------------------
Warning -- This program may not be used in any manner that
competes with the business of Gaussian, Inc. or will provide
assistance to any competitor of Gaussian, Inc.&&The licensee
of this program is prohibited from giving any competitor of
Gaussian, Inc. access to this program.&&By using this program,
the user acknowledges that Gaussian, Inc. is engaged in the
business of creating and licensing software in the field of
computational chemistry and represents and warrants to the
licensee that it is not a competitor of Gaussian, Inc. and that
it will not use this program in any manner prohibited above.
---------------------------------------------------------------
Cite this work as:
Gaussian 09, Revision D.01,
M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria,
M. A. Robb, J. R. Cheeseman, G. Scalmani, V. Barone, B. Mennucci,
G. A. Petersson, H. Nakatsuji, M. Caricato, X. Li, H. P. Hratchian,
A. F. Izmaylov, J. Bloino, G. Zheng, J. L. Sonnenberg, M. Hada,
M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima,
Y. Honda, O. Kitao, H. Nakai, T. Vreven, J. A. Montgomery, Jr.,
J. E. Peralta, F. Ogliaro, M. Bearpark, J. J. Heyd, E. Brothers,
K. N. Kudin, V. N. Staroverov, T. Keith, R. Kobayashi, J. Normand,
K. Raghavachari, A. Rendell, J. C. Burant, S. S. Iyengar, J. Tomasi,
M. Cossi, N. Rega, J. M. Millam, M. Klene, J. E. Knox, J. B. Cross,
V. Bakken, C. Adamo, J. Jaramillo, R. Gomperts, R. E. Stratmann,
O. Yazyev, A. J. Austin, R. Cammi, C. Pomelli, J. W. Ochterski,
R. L. Martin, K. Morokuma, V. G. Zakrzewski, G. A. Voth,
P. Salvador, J. J. Dannenberg, S. Dapprich, A. D. Daniels,
O. Farkas, J. B. Foresman, J. V. Ortiz, J. Cioslowski,
and D. J. Fox, Gaussian, Inc., Wallingford CT, 2013.
******************************************
Gaussian 09:&&AM64L-G09RevD.01 24-Apr-2013
& && && && && & 12-Jul-2015
******************************************
%chk=AAAAAAA.chk
-------------------------------
# hf/3-21g SCRF=(solvent=water)
-------------------------------
2/12=2,17=6,18=5,40=1/2;
3/5=5,11=9,16=1,25=1,30=1,70=/1,2,3;
5/5=2,38=5,53=1/2;
6/7=2,8=2,9=2,10=2,28=1/1;
99/5=1,9=1/99;
-------------------
Title Card Required
-------------------
Symbolic Z-matrix:
Charge =&&0 Multiplicity = 1
C& && && && && && && &0.2832& & 1.09375& &0.
H& && && && && && && &0.63986& &0.08494& &0.
H& && && && && && && &0.63988& &1.59815& &0.87365
H& && && && && && && &0.63988& &1.59815&&-0.87365
H& && && && && && &&&-0.7868& & 1.09376& &0.
& && && && && && && && &&&Input orientation:& && && && && && && && &&&
---------------------------------------------------------------------
Center& &&&Atomic& && &Atomic& && && && & Coordinates (Angstroms)
Number& &&&Number& && & Type& && && && & X& && && &&&Y& && && &&&Z
---------------------------------------------------------------------
& && &1& && && & 6& && && &&&0& && &&&0.283203& & 1.093750& & 0.000000
& && &2& && && & 1& && && &&&0& && &&&0.639858& & 0.084940& & 0.000000
& && &3& && && & 1& && && &&&0& && &&&0.639876& & 1.598148& & 0.873652
& && &4& && && & 1& && && &&&0& && &&&0.639876& & 1.598148& &-0.873652
& && &5& && && & 1& && && &&&0& && & -0.786797& & 1.093763& & 0.000000
---------------------------------------------------------------------
& && && && && && &&&Distance matrix (angstroms):
& && && && && && &&&1& && && & 2& && && & 3& && && & 4& && && & 5
& &&&1&&C& & 0.000000
& &&&2&&H& & 1.070000& &0.000000
& &&&3&&H& & 1.070000& &1.747302& &0.000000
& &&&4&&H& & 1.070000& &1.747302& &1.747303& &0.000000
& &&&5&&H& & 1.070000& &1.747303& &1.747303& &1.747303& &0.000000
Stoichiometry& & CH4
Framework group&&T
Deg. of freedom& &&&1
Full point group& && && && && &&&T& && & NOp&&12
Largest Abelian subgroup& && && &D2& && &NOp& &4
Largest concise Abelian subgroup D2& && &NOp& &4
& && && && && && && && & Standard orientation:& && && && && && && && &
---------------------------------------------------------------------
Center& &&&Atomic& && &Atomic& && && && & Coordinates (Angstroms)
Number& &&&Number& && & Type& && && && & X& && && &&&Y& && && &&&Z
---------------------------------------------------------------------
& && &1& && && & 6& && && &&&0& && &&&0.000000& & 0.000000& & 0.000000
& && &2& && && & 1& && && &&&0& && &&&0.617765& & 0.617765& & 0.617765
& && &3& && && & 1& && && &&&0& && & -0.617765& &-0.617765& & 0.617765
& && &4& && && & 1& && && &&&0& && & -0.617765& & 0.617765& &-0.617765
& && &5& && && & 1& && && &&&0& && &&&0.617765& &-0.617765& &-0.617765
---------------------------------------------------------------------
Rotational constants (GHZ):& & 164.2463791& & 164.2463791& & 164.2463791
Standard basis: 3-21G (6D, 7F)
There are& &&&5 symmetry adapted cartesian basis functions of A& &symmetry.
There are& &&&4 symmetry adapted cartesian basis functions of B1&&symmetry.
There are& &&&4 symmetry adapted cartesian basis functions of B2&&symmetry.
There are& &&&4 symmetry adapted cartesian basis functions of B3&&symmetry.
There are& &&&5 symmetry adapted basis functions of A& &symmetry.
There are& &&&4 symmetry adapted basis functions of B1&&symmetry.
There are& &&&4 symmetry adapted basis functions of B2&&symmetry.
There are& &&&4 symmetry adapted basis functions of B3&&symmetry.
& & 17 basis functions,& & 27 primitive gaussians,& & 17 cartesian basis functions
& &&&5 alpha electrons& && &&&5 beta electrons
& && & nuclear repulsion energy& && &&&13. Hartrees.
NAtoms=& & 5 NActive=& & 5 NUniq=& & 2 SFac= 4.00D+00 NAtFMM=& &60 NAOKFM=F Big=F
Integral buffers will be& & 131072 words long.
Raffenetti 1 integral format.
Two-electron integral symmetry is turned on.
------------------------------------------------------------------------------
Polarizable Continuum Model (PCM)
=================================
Model& && && && && & : PCM (using non-symmetric T matrix).
Atomic radii& && && &: UFF (Universal Force Field).
Polarization charges : Total charges.
Charge compensation&&: None.
Solution method& && &: On-the-fly selection.
Cavity type& && && & : Scaled VdW (van der Waals Surface) (Alpha=1.100).
Cavity algorithm& &&&: GePol (No added spheres)
& && && && && && && && &Default sphere list used, NSphG=& & 5.
& && && && && && && && &Lebedev-Laikov grids with approx.&&5.0 points / Ang**2.
& && && && && && && && &Smoothing algorithm: Karplus/York (Gamma=1.0000).
& && && && && && && && &Polarization charges: spherical gaussians, with
& && && && && && && && && && && && && && && & point-specific exponents (IZeta= 3).
& && && && && && && && &Self-potential: point-specific (ISelfS= 7).
& && && && && && && && &Self-field& & : sphere-specific E.n sum rule (ISelfD= 2).
Solvent& && && && &&&: Water, Eps=&&78.355300 Eps(inf)=& &1.777849
------------------------------------------------------------------------------
Spheres list:
ISph&&on& &Nord& &&&Re0& & Alpha& && &Xe& && && && &Ye& && && && &Ze
& & 1& &C& && &1& & 1.& && &0.000000& && &0.000000& && &0.000000
& & 2& &H& && &2& & 1.& && &0.617765& && &0.617765& && &0.617765
& & 3& &H& && &3& & 1.& &&&-0.617765& &&&-0.617765& && &0.617765
& & 4& &H& && &4& & 1.& &&&-0.617765& && &0.617765& &&&-0.617765
& & 5& &H& && &5& & 1.& && &0.617765& &&&-0.617765& &&&-0.617765
------------------------------------------------------------------------------
One-electron integrals computed using PRISM.
NBasis=& & 17 RedAO= T EigKep=&&4.03D-02&&NBF=& &&&5& &&&4& &&&4& &&&4
NBsUse=& & 17 1.00D-06 EigRej= -1.00D+00 NBFU=& &&&5& &&&4& &&&4& &&&4
ExpMin= 1.83D-01 ExpMax= 1.72D+02 ExpMxC= 1.72D+02 IAcc=1 IRadAn=& && && &1 AccDes= 0.00D+00
Harris functional with IExCor=&&205 and IRadAn=& && & 1 diagonalized for initial guess.
HarFok:&&IExCor=&&205 AccDes= 0.00D+00 IRadAn=& && && &1 IDoV= 1 UseB2=F ITyADJ=14
ICtDFT=&&3500011 ScaDFX=&&1....000000
FoFCou: FMM=F IPFlag=& && && &&&0 FMFlag=& && &100000 FMFlg1=& && && &&&0
& && && &NFxFlg=& && && &&&0 DoJE=T BraDBF=F KetDBF=T FulRan=T
& && && &wScrn=&&0.000000 ICntrl=& &&&500 IOpCl=&&0 I1Cent=& & NGrid=& && && &&&0
& && && &NMat0=& & 1 NMatS0=& && &1 NMatT0=& & 0 NMatD0=& & 1 NMtDS0=& & 0 NMtDT0=& & 0
Petite list used in FoFCou.
Initial guess orbital symmetries:
& && & Occupied&&(A) (A) (T) (T) (T)
& && & Virtual& &(A) (T) (T) (T) (T) (T) (T) (A) (T) (T) (T) (A)
The electronic state of the initial guess is 1-A.
Keep R1 ints in memory in symmetry-blocked form, NReq=3156156.
Requested convergence on RMS density matrix=1.00D-08 within 128 cycles.
Requested convergence on MAX density matrix=1.00D-06.
Requested convergence on& && && && & energy=1.00D-06.
No special actions if energy rises.
Inv3:&&Mode=1 IEnd=& && &696972.
Iteration& & 1 A*A^-1 deviation from unit magnitude is 1.22D-15 for& & 122.
Iteration& & 1 A*A^-1 deviation from orthogonality&&is 1.18D-15 for& & 405& & 245.
Iteration& & 1 A^-1*A deviation from unit magnitude is 1.22D-15 for& & 352.
Iteration& & 1 A^-1*A deviation from orthogonality&&is 6.70D-16 for& & 449& & 401.
Error on total polarization charges =&&0.00401
SCF Done:&&E(RHF) =&&-39.& &&&A.U. after& & 8 cycles
& && && && &NFock=&&8&&Conv=0.29D-09& &&&-V/T= 2.0003
**********************************************************************
& && && && &Population analysis using the SCF density.
**********************************************************************
Orbital symmetries:
& && & Occupied&&(A) (A) (T) (T) (T)
& && & Virtual& &(A) (T) (T) (T) (T) (T) (T) (A) (T) (T) (T) (A)
The electronic state is 1-A.
Alpha&&occ. eigenvalues --&&-11.13818&&-0.95163&&-0.54852&&-0.54852&&-0.54852
Alpha virt. eigenvalues --& & 0.30063& &0.35728& &0.35728& &0.35728& &0.91467
Alpha virt. eigenvalues --& & 0.91467& &0.91467& &1.35371& &1.37195& &1.37195
Alpha virt. eigenvalues --& & 1.37195& &1.96121
& && && & Condensed to atoms (all electrons):
& && && && && &1& && && & 2& && && & 3& && && & 4& && && & 5
& &&&1&&C& & 5.324921& &0.372070& &0.372070& &0.372070& &0.372070
& &&&2&&H& & 0.372070& &0.502499&&-0.025956&&-0.025956&&-0.025956
& &&&3&&H& & 0.372070&&-0.025956& &0.502499&&-0.025956&&-0.025956
& &&&4&&H& & 0.372070&&-0.025956&&-0.025956& &0.502499&&-0.025956
& &&&5&&H& & 0.372070&&-0.025956&&-0.025956&&-0.025956& &0.502499
Mulliken charges:
& && && && && &1
& &&&1&&C& &-0.813200
& &&&2&&H& & 0.203300
& &&&3&&H& & 0.203300
& &&&4&&H& & 0.203300
& &&&5&&H& & 0.203300
Sum of Mulliken charges =& &0.00000
Mulliken charges with hydrogens summed into heavy atoms:
& && && && && &1
& &&&1&&C& & 0.000000
Electronic spatial extent (au):&&&R**2&=& && && && & 34.9726
Charge=& && && && &&&0.0000 electrons
Dipole moment (field-independent basis, Debye):
& & X=& && && && &&&0.0000& & Y=& && && && &&&0.0000& & Z=& && && && &&&0.0000&&Tot=& && && && &&&0.0000
Quadrupole moment (field-independent basis, Debye-Ang):
& &XX=& && && && & -8.3475& &YY=& && && && & -8.3475& &ZZ=& && && && & -8.3475
& &XY=& && && && &&&0.0000& &XZ=& && && && &&&0.0000& &YZ=& && && && &&&0.0000
Traceless Quadrupole moment (field-independent basis, Debye-Ang):
& &XX=& && && && &&&0.0000& &YY=& && && && &&&0.0000& &ZZ=& && && && &&&0.0000
& &XY=& && && && &&&0.0000& &XZ=& && && && &&&0.0000& &YZ=& && && && &&&0.0000
Octapole moment (field-independent basis, Debye-Ang**2):
&&XXX=& && && && &&&0.0000&&YYY=& && && && &&&0.0000&&ZZZ=& && && && &&&0.0000&&XYY=& && && && &&&0.0000
&&XXY=& && && && &&&0.0000&&XXZ=& && && && &&&0.0000&&XZZ=& && && && &&&0.0000&&YZZ=& && && && &&&0.0000
&&YYZ=& && && && &&&0.0000&&XYZ=& && && && &&&0.6960
Hexadecapole moment (field-independent basis, Debye-Ang**3):
XXXX=& && && && &-15.0133 YYYY=& && && && &-15.0133 ZZZZ=& && && && &-15.0133 XXXY=& && && && &&&0.0000
XXXZ=& && && && &&&0.0000 YYYX=& && && && &&&0.0000 YYYZ=& && && && &&&0.0000 ZZZX=& && && && &&&0.0000
ZZZY=& && && && &&&0.0000 XXYY=& && && && & -4.5795 XXZZ=& && && && & -4.5795 YYZZ=& && && && & -4.5795
XXYZ=& && && && &&&0.0000 YYXZ=& && && && &&&0.0000 ZZXY=& && && && &&&0.0000
N-N= 1.D+01 E-N=-1.D+02&&KE= 3.D+01
Symmetry A& & KE= 3.D+01
Symmetry B1& &KE= 1.D+00
Symmetry B2& &KE= 1.D+00
Symmetry B3& &KE= 1.D+00
1\1\GINC-317051\SP\RHF\3-21G\C1H4\ROOT\12-Jul-2015\0\\# hf/3-21g SCRF=
(solvent=water)\\Title Card Required\\0,1\C,0,0...
\H,0,0...\H,0,0...8736515\H,0,
0...8736515\H,0,-0...\\Versio
n=AM64L-G09RevD.01\State=1-A\HF=-39.9765944\RMSD=2.930e-10\Dipole=0.,0
.,0.\Quadrupole=0.,0.,0.,0.,0.,0.\PG=T \\@
THE FIRST AND LAST THING REQUIRED OF GENIUS IS THE LOVE OF TRUTH.
& && && && && && && && &&&-- GOETHE
Job cpu time:& && & 0 days&&0 hours&&0 minutes&&2.0 seconds.
File lengths (MBytes):&&RWF=& &&&10 Int=& && &0 D2E=& && &0 Chk=& && &1 Scr=& && &1
Normal termination of Gaussian 09 at Sun Jul 12 09:46:34 2015.
然后我就不知道该看输出文件中的那些数据才能计算出我想要的甲烷的分子电子密度体积V。
谢谢不吝赐教。
恕我愚钝，看了好几遍这两个关键词也没有找到与计算分子电子密度体积的联系。请问 ，Volume是用于计算分子体积的，分子体积与分子的电子密度体积之间的关系是怎么样的呢？他俩是等同的么？
我也用Volume关键词算了一个乙烷分子，但是我不知道应该怎样看输出文件。
以下便是输出文件：
Entering Gaussian System, Link 0=g09
Output=c2h6volume.log
Initial command:
/root/gaussian/g09/l1.exe &/tmp/Gau-9690.inp& -scrdir=&/tmp/&
Entering Link 1 = /root/gaussian/g09/l1.exe PID=& && &9691.
Copyright (c) 92,98,13,
& && && && &Gaussian, Inc.&&All Rights Reserved.
This is part of the Gaussian(R) 09 program.&&It is based on
the Gaussian(R) 03 system (copyright 2003, Gaussian, Inc.),
the Gaussian(R) 98 system (copyright 1998, Gaussian, Inc.),
the Gaussian(R) 94 system (copyright 1995, Gaussian, Inc.),
the Gaussian 92(TM) system (copyright 1992, Gaussian, Inc.),
the Gaussian 90(TM) system (copyright 1990, Gaussian, Inc.),
the Gaussian 88(TM) system (copyright 1988, Gaussian, Inc.),
the Gaussian 86(TM) system (copyright 1986, Carnegie Mellon
University), and the Gaussian 82(TM) system (copyright 1983,
Carnegie Mellon University). Gaussian is a federally registered
trademark of Gaussian, Inc.
This software contains proprietary and confidential information,
including trade secrets, belonging to Gaussian, Inc.
This software is provided under written license and may be
used, copied, transmitted, or stored only in accord with that
written license.
The following legend is applicable only to US Government
contracts under FAR:
& && && && && && &&&RESTRICTED RIGHTS LEGEND
Use, reproduction and disclosure by the US Government is
subject to restrictions as set forth in subparagraphs (a)
and (c) of the Commercial Computer Software - Restricted
Rights clause in FAR 52.227-19.
Gaussian, Inc.
340 Quinnipiac St., Bldg. 40, Wallingford CT 06492
---------------------------------------------------------------
Warning -- This program may not be used in any manner that
competes with the business of Gaussian, Inc. or will provide
assistance to any competitor of Gaussian, Inc.&&The licensee
of this program is prohibited from giving any competitor of
Gaussian, Inc. access to this program.&&By using this program,
the user acknowledges that Gaussian, Inc. is engaged in the
business of creating and licensing software in the field of
computational chemistry and represents and warrants to the
licensee that it is not a competitor of Gaussian, Inc. and that
it will not use this program in any manner prohibited above.
---------------------------------------------------------------
Cite this work as:
Gaussian 09, Revision D.01,
M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria,
M. A. Robb, J. R. Cheeseman, G. Scalmani, V. Barone, B. Mennucci,
G. A. Petersson, H. Nakatsuji, M. Caricato, X. Li, H. P. Hratchian,
A. F. Izmaylov, J. Bloino, G. Zheng, J. L. Sonnenberg, M. Hada,
M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima,
Y. Honda, O. Kitao, H. Nakai, T. Vreven, J. A. Montgomery, Jr.,
J. E. Peralta, F. Ogliaro, M. Bearpark, J. J. Heyd, E. Brothers,
K. N. Kudin, V. N. Staroverov, T. Keith, R. Kobayashi, J. Normand,
K. Raghavachari, A. Rendell, J. C. Burant, S. S. Iyengar, J. Tomasi,
M. Cossi, N. Rega, J. M. Millam, M. Klene, J. E. Knox, J. B. Cross,
V. Bakken, C. Adamo, J. Jaramillo, R. Gomperts, R. E. Stratmann,
O. Yazyev, A. J. Austin, R. Cammi, C. Pomelli, J. W. Ochterski,
R. L. Martin, K. Morokuma, V. G. Zakrzewski, G. A. Voth,
P. Salvador, J. J. Dannenberg, S. Dapprich, A. D. Daniels,
O. Farkas, J. B. Foresman, J. V. Ortiz, J. Cioslowski,
and D. J. Fox, Gaussian, Inc., Wallingford CT, 2013.
******************************************
Gaussian 09:&&AM64L-G09RevD.01 24-Apr-2013
& && && && && & 13-Jul-2015
******************************************
%chk=c2h6volume.chk
----------------------------------------------
#P B3LYP/6-31G(d,p) geom=connectivity # Volume
----------------------------------------------
1/38=1,57=2/1;
2/12=2,17=6,18=5,40=1/2;
3/5=1,6=6,7=101,11=2,16=1,25=1,30=1,74=-5/1,2,3;
5/5=2,38=5/2;
6/7=2,8=2,9=2,10=2,28=1,44=1/1,4;
99/5=1,9=1/99;
Leave Link& & 1 at Mon Jul 13 10:21:16 2015, MaxMem=& && && &&&0 cpu:& && && &0.0
(Enter /root/gaussian/g09/l101.exe)
-------------------
Title Card Required
-------------------
Symbolic Z-matrix:
Charge =&&0 Multiplicity = 1
C& && && && && && &&&-0.32227&&-1.81641& &0.
C& && && && && && && &1.17798&&-1.81641& &0.
H& && && && && && &&&-0.71738&&-0.77148& &0.
H& && && && && && &&&-0.7176& &-2.33875&&-0.90492
H& && && && && && &&&-0.71766&&-2.33879& &0.90486
H& && && && && && && &1.57348&&-2.86109& &0.00006
H& && && && && && && &1.57301&&-1.29386&&-0.90498
H& && && && && && && &1.57302&&-1.29378& &0.90492
NAtoms=& && &8 NQM=& && &&&8 NQMF=& && & 0 NMMI=& && &0 NMMIF=& && &0
& && && && && & NMic=& && & 0 NMicF=& && &0.
& && && && && && &&&Isotopes and Nuclear Properties:
(Nuclear quadrupole moments (NQMom) in fm**2, nuclear magnetic moments (NMagM)
&&in nuclear magnetons)
&&Atom& && && &1& && && &&&2& && && &&&3& && && &&&4& && && &&&5& && && &&&6& && && &&&7& && && &&&8
IAtWgt=& && && & 12& && && & 12& && && &&&1& && && &&&1& && && &&&1& && && &&&1& && && &&&1& && && &&&1
AtmWgt=&&12..0000000& &1.0078250& &1.0078250& &1.0078250& &1.0078250& &1.0078250& &1.0078250
NucSpn=& && && &&&0& && && &&&0& && && &&&1& && && &&&1& && && &&&1& && && &&&1& && && &&&1& && && &&&1
AtZEff=& &0.0000000& &0.0000000& &0.0000000& &0.0000000& &0.0000000& &0.0000000& &0.0000000& &0.0000000
NQMom=& & 0.0000000& &0.0000000& &0.0000000& &0.0000000& &0.0000000& &0.0000000& &0.0000000& &0.0000000
NMagM=& & 0.0000000& &0.0000000& &2.7928460& &2.7928460& &2.7928460& &2.7928460& &2.7928460& &2.7928460
AtZNuc=& &6.0000000& &6.0000000& &1.0000000& &1.0000000& &1.0000000& &1.0000000& &1.0000000& &1.0000000
Leave Link&&101 at Mon Jul 13 10:21:16 2015, MaxMem=& &
cpu:& && && &0.1
(Enter /root/gaussian/g09/l202.exe)
& && && && && && && && &&&Input orientation:& && && && && && && && &&&
---------------------------------------------------------------------
Center& &&&Atomic& && &Atomic& && && && & Coordinates (Angstroms)
Number& &&&Number& && & Type& && && && & X& && && &&&Y& && && &&&Z
---------------------------------------------------------------------
& && &1& && && & 6& && && &&&0& && & -0.322266& &-1.816406& & 0.000000
& && &2& && && & 6& && && &&&0& && &&&1.177984& &-1.816406& & 0.000000
& && &3& && && & 1& && && &&&0& && & -0.717380& &-0.771475& & 0.000000
& && &4& && && & 1& && && &&&0& && & -0.717602& &-2.338748& &-0.904922
& && &5& && && & 1& && && &&&0& && & -0.717661& &-2.338793& & 0.904863
& && &6& && && & 1& && && &&&0& && &&&1.573484& &-2.861085& & 0.000060
& && &7& && && & 1& && && &&&0& && &&&1.573006& &-1.293864& &-0.904978
& && &8& && && & 1& && && &&&0& && &&&1.573024& &-1.293782& & 0.904921
---------------------------------------------------------------------
& && && && && && &&&Distance matrix (angstroms):
& && && && && && &&&1& && && & 2& && && & 3& && && & 4& && && & 5
& &&&1&&C& & 0.000000
& &&&2&&C& & 1.500250& &0.000000
& &&&3&&H& & 1.117137& &2.164321& &0.000000
& &&&4&&H& & 1.117146& &2.164479& &1.809759& &0.000000
& &&&5&&H& & 1.117140& &2.164517& &1.809769& &1.809785& &0.000000
& &&&6&&H& & 2.164537& &1.117038& &3.100730& &2.518115& &2.518095
& &&&7&&H& & 2.164276& &1.117174& &2.517488& &2.517671& &3.100735
& &&&8&&H& & 2.164288& &1.117173& &2.517466& &3.100718& &2.517794
& && && && && && &&&6& && && & 7& && && & 8
& &&&6&&H& & 0.000000
& &&&7&&H& & 1.809772& &0.000000
& &&&8&&H& & 1.809755& &1.809899& &0.000000
Stoichiometry& & C2H6
Framework group&&C1
Deg. of freedom& & 18
Full point group& && && && && &&&C1& && &NOp& &1
Largest Abelian subgroup& && && &C1& && &NOp& &1
Largest concise Abelian subgroup C1& && &NOp& &1
& && && && && && && && & Standard orientation:& && && && && && && && &
---------------------------------------------------------------------
Center& &&&Atomic& && &Atomic& && && && & Coordinates (Angstroms)
Number& &&&Number& && & Type& && && && & X& && && &&&Y& && && &&&Z
---------------------------------------------------------------------
& && &1& && && & 6& && && &&&0& && & -0.750109& &-0.000008& &-0.000031
& && &2& && && & 6& && && &&&0& && &&&0.750141& &-0.000004& &-0.000045
& && &3& && && & 1& && && &&&0& && & -1.145214& & 0.077247& & 1.042043
& && &4& && && & 1& && && &&&0& && & -1.145453& & 0.863818& &-0.587844
& && &5& && && & 1& && && &&&0& && & -1.145506& &-0.941017& &-0.454085
& && &6& && && & 1& && && &&&0& && &&&1.145632& &-0.077300& &-1.041864
& && &7& && && & 1& && && &&&0& && &&&1.145165& & 0.941131& & 0.454155
& && &8& && && & 1& && && &&&0& && &&&1.145188& &-0.863808& & 0.588050
---------------------------------------------------------------------
Rotational constants (GHZ):& &&&76.5501021& &&&20.4292928& &&&20.4289749
Leave Link&&202 at Mon Jul 13 10:21:16 2015, MaxMem=& &
cpu:& && && &0.0
(Enter /root/gaussian/g09/l301.exe)
Standard basis: 6-31G(d,p) (6D, 7F)
Ernie: Thresh=&&0.10000D-02 Tol=&&0.10000D-05 Strict=F.
There are& & 60 symmetry adapted cartesian basis functions of A& &symmetry.
There are& & 60 symmetry adapted basis functions of A& &symmetry.
& & 60 basis functions,& & 98 primitive gaussians,& & 60 cartesian basis functions
& &&&9 alpha electrons& && &&&9 beta electrons
& && & nuclear repulsion energy& && &&&42. Hartrees.
IExCor=&&402 DFT=T Ex+Corr=B3LYP ExCW=0 ScaHFX=&&0.200000
ScaDFX=&&0....810000 ScalE2=&&1..000000
IRadAn=& && &0 IRanWt=& &&&-1 IRanGd=& && && && &0 ICorTp=0 IEmpDi=&&4
NAtoms=& & 8 NActive=& & 8 NUniq=& & 8 SFac= 1.00D+00 NAtFMM=& &60 NAOKFM=F Big=F
Integral buffers will be& & 131072 words long.
Raffenetti 2 integral format.
Two-electron integral symmetry is turned on.
Leave Link&&301 at Mon Jul 13 10:21:17 2015, MaxMem=& &
cpu:& && && &0.0
(Enter /root/gaussian/g09/l302.exe)
NPDir=0 NMtPBC=& &&&1 NCelOv=& &&&1 NCel=& && & 1 NClECP=& &&&1 NCelD=& && &1
& && && &NCelK=& && &1 NCelE2=& &&&1 NClLst=& &&&1 CellRange=& &&&0.0.
One-electron integrals computed using PRISM.
NBasis=& & 60 RedAO= T EigKep=&&7.75D-03&&NBF=& & 60
NBsUse=& & 60 1.00D-06 EigRej= -1.00D+00 NBFU=& & 60
Precomputing XC quadrature grid using
IXCGrd= 4 IRadAn=& && && &&&0 IRanWt=& && && & -1 IRanGd=& && && &&&0 AccXCQ= 0.00D+00.
Generated NRdTot=& && & 0 NPtTot=& && && &&&0 NUsed=& && && &&&0 NTot=& && && & 32
NSgBfM=& & 60& & 60& & 60& & 60& & 60 MxSgAt=& &&&8 MxSgA2=& &&&8.
Leave Link&&302 at Mon Jul 13 10:21:17 2015, MaxMem=& &
cpu:& && && &0.2
(Enter /root/gaussian/g09/l303.exe)
DipDrv:&&MaxL=1.
Leave Link&&303 at Mon Jul 13 10:21:17 2015, MaxMem=& &
cpu:& && && &0.0
(Enter /root/gaussian/g09/l401.exe)
ExpMin= 1.61D-01 ExpMax= 3.05D+03 ExpMxC= 4.57D+02 IAcc=1 IRadAn=& && && &1 AccDes= 0.00D+00
Harris functional with IExCor=&&402 and IRadAn=& && & 1 diagonalized for initial guess.
HarFok:&&IExCor=&&402 AccDes= 0.00D+00 IRadAn=& && && &1 IDoV= 1 UseB2=F ITyADJ=14
ICtDFT=&&3500011 ScaDFX=&&1....000000
FoFCou: FMM=F IPFlag=& && && &&&0 FMFlag=& && &100000 FMFlg1=& && && &&&0
& && && &NFxFlg=& && && &&&0 DoJE=T BraDBF=F KetDBF=T FulRan=T
& && && &wScrn=&&0.000000 ICntrl=& &&&500 IOpCl=&&0 I1Cent=& & NGrid=& && && &&&0
& && && &NMat0=& & 1 NMatS0=& && &1 NMatT0=& & 0 NMatD0=& & 1 NMtDS0=& & 0 NMtDT0=& & 0
Petite list used in FoFCou.
Harris En= -79.5& &
JPrj=0 DoOrth=F DoCkMO=F.
Leave Link&&401 at Mon Jul 13 10:21:17 2015, MaxMem=& &
cpu:& && && &0.5
(Enter /root/gaussian/g09/l502.exe)
Closed shell SCF:
Using DIIS extrapolation, IDIIS=&&1040.
Integral symmetry usage will be decided dynamically.
Keep R1 ints in memory in canonical form, NReq=2587035.
IVT=& && & 30172 IEndB=& && & 30172 NGot=& &
Requested convergence on RMS density matrix=1.00D-08 within 128 cycles.
Requested convergence on MAX density matrix=1.00D-06.
Requested convergence on& && && && & energy=1.00D-06.
No special actions if energy rises.
FoFCou: FMM=F IPFlag=& && && &&&0 FMFlag=& && && &&&0 FMFlg1=& && && &&&0
& && && &NFxFlg=& && && &&&0 DoJE=F BraDBF=F KetDBF=F FulRan=T
& && && &wScrn=&&0.000000 ICntrl=& &&&600 IOpCl=&&0 I1Cent=& && && &&&0 NGrid=& && && &&&0
& && && &NMat0=& & 1 NMatS0=& &1830 NMatT0=& & 0 NMatD0=& & 1 NMtDS0=& & 0 NMtDT0=& & 0
Symmetry not used in FoFCou.
Integral accuracy reduced to 1.0D-05 until final iterations.
Cycle& &1&&Pass 0&&IDiag&&1:
E= -79.5& &
DIIS: error= 4.34D-02 at cycle& &1 NSaved=& &1.
NSaved= 1 IEnMin= 1 EnMin= -79.5& &&&IErMin= 1 ErrMin= 4.34D-02
ErrMax= 4.34D-02&&0.00D+00 EMaxC= 1.00D-01 BMatC= 1.09D-01 BMatP= 1.09D-01
IDIUse=3 WtCom= 5.66D-01 WtEn= 4.34D-01
Coeff-Com:&&0.100D+01
Coeff-En:& &0.100D+01
Coeff:& && &0.100D+01
Gap=& &&&0.373 Goal=& &None& & Shift=& & 0.000
GapD=& & 0.373 DampG=1.000 DampE=0.500 DampFc=0.5000 IDamp=-1.
Damping current iteration by 5.00D-01
RMSDP=1.14D-02 MaxDP=1.36D-01& && && && &&&OVMax= 1.88D-01
Cycle& &2&&Pass 0&&IDiag&&1:
E= -79.1& &&&Delta-E=& && & -0. Rises=F Damp=T
DIIS: error= 1.15D-02 at cycle& &2 NSaved=& &2.
NSaved= 2 IEnMin= 2 EnMin= -79.1& &&&IErMin= 2 ErrMin= 1.15D-02
ErrMax= 1.15D-02&&0.00D+00 EMaxC= 1.00D-01 BMatC= 4.30D-03 BMatP= 1.09D-01
IDIUse=3 WtCom= 8.85D-01 WtEn= 1.15D-01
Coeff-Com:&&0.101D+00 0.899D+00
Coeff-En:& &0.160D+00 0.840D+00
Coeff:& && &0.108D+00 0.892D+00
Gap=& &&&0.464 Goal=& &None& & Shift=& & 0.000
RMSDP=1.25D-03 MaxDP=1.29D-02 DE=-5.36D-02 OVMax= 9.84D-02
Cycle& &3&&Pass 0&&IDiag&&1:
E= -79.9& &&&Delta-E=& && & -0. Rises=F Damp=F
DIIS: error= 2.02D-03 at cycle& &3 NSaved=& &3.
NSaved= 3 IEnMin= 3 EnMin= -79.9& &&&IErMin= 3 ErrMin= 2.02D-03
ErrMax= 2.02D-03&&0.00D+00 EMaxC= 1.00D-01 BMatC= 1.10D-04 BMatP= 4.30D-03
IDIUse=3 WtCom= 9.80D-01 WtEn= 2.02D-02
Coeff-Com:&&0.109D-01 0.821D-01 0.907D+00
Coeff-En:& &0.000D+00 0.000D+00 0.100D+01
Coeff:& && &0.107D-01 0.805D-01 0.909D+00
Gap=& &&&0.433 Goal=& &None& & Shift=& & 0.000
RMSDP=1.91D-04 MaxDP=1.60D-03 DE=-6.99D-02 OVMax= 4.05D-03
Cycle& &4&&Pass 0&&IDiag&&1:
E= -79.4& &&&Delta-E=& && & -0. Rises=F Damp=F
DIIS: error= 6.61D-04 at cycle& &4 NSaved=& &4.
NSaved= 4 IEnMin= 4 EnMin= -79.4& &&&IErMin= 4 ErrMin= 6.61D-04
ErrMax= 6.61D-04&&0.00D+00 EMaxC= 1.00D-01 BMatC= 9.20D-06 BMatP= 1.10D-04
IDIUse=3 WtCom= 9.93D-01 WtEn= 6.61D-03
Coeff-Com:&&0.157D-02-0.490D-02 0.212D+00 0.791D+00
Coeff-En:& &0.000D+00 0.000D+00 0.000D+00 0.100D+01
Coeff:& && &0.156D-02-0.487D-02 0.211D+00 0.793D+00
Gap=& &&&0.433 Goal=& &None& & Shift=& & 0.000
RMSDP=4.37D-05 MaxDP=3.77D-04 DE=-9.39D-05 OVMax= 1.03D-03
Cycle& &5&&Pass 0&&IDiag&&1:
E= -79.5& &&&Delta-E=& && & -0. Rises=F Damp=F
DIIS: error= 1.43D-05 at cycle& &5 NSaved=& &5.
NSaved= 5 IEnMin= 5 EnMin= -79.5& &&&IErMin= 5 ErrMin= 1.43D-05
ErrMax= 1.43D-05&&0.00D+00 EMaxC= 1.00D-01 BMatC= 7.99D-09 BMatP= 9.20D-06
IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00
Coeff-Com:&&0.135D-03-0.452D-03-0.792D-02-0.139D-01 0.102D+01
Coeff:& && &0.135D-03-0.452D-03-0.792D-02-0.139D-01 0.102D+01
Gap=& &&&0.433 Goal=& &None& & Shift=& & 0.000
RMSDP=1.38D-06 MaxDP=2.08D-05 DE=-8.56D-06 OVMax= 3.36D-05
Initial convergence to 1.0D-05 achieved.&&Increase integral accuracy.
Cycle& &6&&Pass 1&&IDiag&&1:
E= -79.2& &&&Delta-E=& && & -0. Rises=F Damp=F
DIIS: error= 6.62D-06 at cycle& &1 NSaved=& &1.
NSaved= 1 IEnMin= 1 EnMin= -79.2& &&&IErMin= 1 ErrMin= 6.62D-06
ErrMax= 6.62D-06&&0.00D+00 EMaxC= 1.00D-01 BMatC= 1.29D-09 BMatP= 1.29D-09
IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00
Coeff-Com:&&0.100D+01
Coeff:& && &0.100D+01
Gap=& &&&0.433 Goal=& &None& & Shift=& & 0.000
RMSDP=1.38D-06 MaxDP=2.08D-05 DE=-2.31D-05 OVMax= 1.68D-05
Cycle& &7&&Pass 1&&IDiag&&1:
E= -79.8& &&&Delta-E=& && & -0. Rises=F Damp=F
DIIS: error= 5.13D-07 at cycle& &2 NSaved=& &2.
NSaved= 2 IEnMin= 2 EnMin= -79.8& &&&IErMin= 2 ErrMin= 5.13D-07
ErrMax= 5.13D-07&&0.00D+00 EMaxC= 1.00D-01 BMatC= 1.81D-11 BMatP= 1.29D-09
IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00
Coeff-Com: -0.159D-01 0.102D+01
Coeff:& &&&-0.159D-01 0.102D+01
Gap=& &&&0.433 Goal=& &None& & Shift=& & 0.000
RMSDP=1.41D-07 MaxDP=1.93D-06 DE=-1.72D-09 OVMax= 2.13D-06
Cycle& &8&&Pass 1&&IDiag&&1:
E= -79.3& &&&Delta-E=& && & -0. Rises=F Damp=F
DIIS: error= 4.81D-07 at cycle& &3 NSaved=& &3.
NSaved= 3 IEnMin= 3 EnMin= -79.3& &&&IErMin= 3 ErrMin= 4.81D-07
ErrMax= 4.81D-07&&0.00D+00 EMaxC= 1.00D-01 BMatC= 1.86D-11 BMatP= 1.81D-11
IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00
Coeff-Com: -0.156D-01 0.520D+00 0.496D+00
Coeff:& &&&-0.156D-01 0.520D+00 0.496D+00
Gap=& &&&0.433 Goal=& &None& & Shift=& & 0.000
RMSDP=6.36D-08 MaxDP=1.13D-06 DE=-1.51D-12 OVMax= 1.16D-06
Cycle& &9&&Pass 1&&IDiag&&1:
E= -79.4& &&&Delta-E=& && & -0. Rises=F Damp=F
DIIS: error= 5.26D-08 at cycle& &4 NSaved=& &4.
NSaved= 4 IEnMin= 4 EnMin= -79.4& &&&IErMin= 4 ErrMin= 5.26D-08
ErrMax= 5.26D-08&&0.00D+00 EMaxC= 1.00D-01 BMatC= 1.73D-13 BMatP= 1.81D-11
IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00
Coeff-Com: -0.211D-02 0.184D-01 0.940D-01 0.890D+00
Coeff:& &&&-0.211D-02 0.184D-01 0.940D-01 0.890D+00
Gap=& &&&0.433 Goal=& &None& & Shift=& & 0.000
RMSDP=6.92D-09 MaxDP=1.04D-07 DE=-1.91D-11 OVMax= 1.23D-07
SCF Done:&&E(RB3LYP) =&&-79.& &&&A.U. after& & 9 cycles
& && && && &NFock=&&9&&Conv=0.69D-08& &&&-V/T= 2.0122
KE= 7.D+01 PE=-2.D+02 EE= 6.D+01
Leave Link&&502 at Mon Jul 13 10:21:23 2015, MaxMem=& &
cpu:& && && &5.4
(Enter /root/gaussian/g09/l601.exe)
Copying SCF densities to generalized density rwf, IOpCl= 0 IROHF=0.
**********************************************************************
& && && && &Population analysis using the SCF density.
**********************************************************************
Orbital symmetries:
& && & Occupied&&(A) (A) (A) (A) (A) (A) (A) (A) (A)
& && & Virtual& &(A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A)
& && && && && &&&(A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A)
& && && && && &&&(A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A)
& && && && && &&&(A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A)
& && && && && &&&(A) (A) (A)
The electronic state is 1-A.
Alpha&&occ. eigenvalues --&&-10.155&&-0.75072&&-0.60130&&-0.43058
Alpha&&occ. eigenvalues --& &-0.43058&&-0.36446&&-0.33384&&-0.33383
Alpha virt. eigenvalues --& & 0.09886& &0.15436& &0.15474& &0.15475& &0.18448
Alpha virt. eigenvalues --& & 0.18450& &0.24898& &0.49854& &0.55185& &0.55189
Alpha virt. eigenvalues --& & 0.61661& &0.61663& &0.67898& &0.84930& &0.84933
Alpha virt. eigenvalues --& & 0.85719& &0.89012& &0.89013& &0.93343& &1.06084
Alpha virt. eigenvalues --& & 1.32695& &1.32699& &1.44217& &1.69787& &1.69792
Alpha virt. eigenvalues --& & 1.89326& &2.03021& &2.03026& &2.04099& &2.04109
Alpha virt. eigenvalues --& & 2.04593& &2.04651& &2.26744& &2.26744& &2.33459
Alpha virt. eigenvalues --& & 2.47314& &2.47322& &2.65351& &2.72486& &2.76015
Alpha virt. eigenvalues --& & 2.76018& &2.87399& &2.87407& &3.09011& &3.19201
Alpha virt. eigenvalues --& & 3.34763& &3.34775& &3.48841& &3.48847& &4.32227
Alpha virt. eigenvalues --& & 4.58207
& && && & Condensed to atoms (all electrons):
& && && && && &1& && && & 2& && && & 3& && && & 4& && && & 5& && && & 6
& &&&1&&C& & 4.879104& &0.373327& &0.388080& &0.388071& &0.388073&&-0.037405
& &&&2&&C& & 0.373327& &4.879066&&-0.037423&&-0.037414&&-0.037413& &0.388091
& &&&3&&H& & 0.388080&&-0.037423& &0.614298&&-0.031292&&-0.031292& &0.005241
& &&&4&&H& & 0.388071&&-0.037414&&-0.031292& &0.614292&&-0.031290&&-0.004530
& &&&5&&H& & 0.388073&&-0.037413&&-0.031292&&-0.031290& &0.614286&&-0.004531
& &&&6&&H& &-0.037405& &0.388091& &0.005241&&-0.004530&&-0.004531& &0.614230
& &&&7&&H& &-0.037431& &0.388074&&-0.004537&&-0.004535& &0.005242&&-0.031292
& &&&8&&H& &-0.037428& &0.388073&&-0.004538& &0.005242&&-0.004533&&-0.031292
& && && && && &7& && && & 8
& &&&1&&C& &-0.037431&&-0.037428
& &&&2&&C& & 0.388074& &0.388073
& &&&3&&H& &-0.004537&&-0.004538
& &&&4&&H& &-0.004535& &0.005242
& &&&5&&H& & 0.005242&&-0.004533
& &&&6&&H& &-0.031292&&-0.031292
& &&&7&&H& & 0.614298&&-0.031276
& &&&8&&H& &-0.031276& &0.614297
Mulliken charges:
& && && && && &1
& &&&1&&C& &-0.304391
& &&&2&&C& &-0.304383
& &&&3&&H& & 0.101463
& &&&4&&H& & 0.101455
& &&&5&&H& & 0.101457
& &&&6&&H& & 0.101488
& &&&7&&H& & 0.101456
& &&&8&&H& & 0.101455
Sum of Mulliken charges =& &0.00000
Mulliken charges with hydrogens summed into heavy atoms:
& && && && && &1
& &&&1&&C& &-0.000016
& &&&2&&C& & 0.000016
Electronic spatial extent (au):&&&R**2&=& && && && &109.0983
Charge=& && && && &&&0.0000 electrons
Dipole moment (field-independent basis, Debye):
& & X=& && && && &&&0.0000& & Y=& && && && &&&0.0000& & Z=& && && && &&&0.0000&&Tot=& && && && &&&0.0000
Quadrupole moment (field-independent basis, Debye-Ang):
& &XX=& && && && &-15.5294& &YY=& && && && &-14.7542& &ZZ=& && && && &-14.7542
& &XY=& && && && &&&0.0000& &XZ=& && && && & -0.0004& &YZ=& && && && &&&0.0000
Traceless Quadrupole moment (field-independent basis, Debye-Ang):
& &XX=& && && && & -0.5168& &YY=& && && && &&&0.2585& &ZZ=& && && && &&&0.2584
& &XY=& && && && &&&0.0000& &XZ=& && && && & -0.0004& &YZ=& && && && &&&0.0000
Octapole moment (field-independent basis, Debye-Ang**2):
&&XXX=& && && && & -0.0004&&YYY=& && && && &&&0.0002&&ZZZ=& && && && &&&0.0000&&XYY=& && && && & -0.0001
&&XXY=& && && && & -0.0002&&XXZ=& && && && & -0.0014&&XZZ=& && && && &&&0.0004&&YZZ=& && && && & -0.0001
&&YYZ=& && && && &&&0.0001&&XYZ=& && && && &&&0.0000
Hexadecapole moment (field-independent basis, Debye-Ang**3):
XXXX=& && && && &-91.7168 YYYY=& && && && &-29.8248 ZZZZ=& && && && &-29.8228 XXXY=& && && && &&&0.0000
XXXZ=& && && && & -0.0018 YYYX=& && && && &&&0.2830 YYYZ=& && && && &&&0.0001 ZZZX=& && && && & -1.2543
ZZZY=& && && && &&&0.0001 XXYY=& && && && &-19.1708 XXZZ=& && && && &-19.1701 YYZZ=& && && && & -9.9411
XXYZ=& && && && &&&0.0001 YYXZ=& && && && &&&1.2534 ZZXY=& && && && & -0.2831
N-N= 4.D+01 E-N=-2.D+02&&KE= 7.D+01
No NMR shielding tensors so no spin-rotation constants.
Leave Link&&601 at Mon Jul 13 10:21:23 2015, MaxMem=& &
cpu:& && && &0.2
(Enter /root/gaussian/g09/l604.exe)
Evaluate density.
Using the total density.
Monte-Carlo method of calculating molar volume:
based on&&0.001 e/bohr**3 density envelope.
Number of points per bohr**3 =& &20 CutOff= 1.00D-04
Using the SCF density.
There are& && & 278 points.&&Will hold& & 115294 in memory.
Box volume =
fraction occupied=0.129
Integrated density= 8.5422D+00 error=-9.4578D+00
Molar volume =&&348.471 bohr**3/mol ( 31.097 cm**3/mol)
Recommended a0 for SCRF calculation =&&3.04 angstrom (&&5.75 bohr)
Leave Link&&604 at Mon Jul 13 10:21:23 2015, MaxMem=& &
cpu:& && && &0.0
(Enter /root/gaussian/g09/l9999.exe)
1\1\GINC-317051\SP\RB3LYP\6-31G(d,p)\C2H6\ROOT\13-Jul-2015\0\\#P B3LYP
/6-31G(d,p) geom=connectivity # Volume\\Title Card Required\\0,1\C,0,-
0...\C,0,1...\H,0,-0.7173796
4,-0..\H,0,-0...904922\H,0,-0.
,-2..904863\H,0,1...00006\H,0,1.5730063
6,-1..904978\H,0,1...904921\\Version=A
M64L-G09RevD.01\State=1-A\HF=-79.8363234\RMSD=6.920e-09\Dipole=0.00000
33,0..0000021\Quadrupole=-0....0
..0000379\PG=C01 \\@
YOU SHOULDN'T LET PEOPLE DRIVE YOU CRAZY
WHEN YOU KNOW IT'S WITHIN WALKING DISTANCE.
Job cpu time:& && & 0 days&&0 hours&&0 minutes&&6.7 seconds.
File lengths (MBytes):&&RWF=& && &5 Int=& && &0 D2E=& && &0 Chk=& && &1 Scr=& && &1
Normal termination of Gaussian 09 at Mon Jul 13 10:21:23 2015.
还望您试点一下，您的指点将对我有莫大的帮助。谢谢。
我不清楚这个“分子电子密度体积”具体指什么，文献里有表格列举吗？最好能知道单位什么的。。。。。。
输出文件的这段
(Enter /root/gaussian/g09/l604.exe)
Evaluate density.
Using the total density.
Monte-Carlo method of calculating molar volume:
based on&&0.001 e/bohr**3 density envelope.
Number of points per bohr**3 =& &20 CutOff= 1.00D-04
Using the SCF density.
There are& && & 278 points.&&Will hold& & 115294 in memory.
Box volume =
fraction occupied=0.129
Integrated density= 8.5422D+00 error=-9.4578D+00
Molar volume =&&348.471 bohr**3/mol ( 31.097 cm**3/mol)
Recommended a0 for SCRF calculation =&&3.04 angstrom (&&5.75 bohr)
Leave Link&&604 at Mon Jul 13 10:21:23 2015, MaxMem=& &
cpu:& && && &0.0
应该是在算的分子体积
您好 谢谢您的回复。
我是看到了一篇有关于计算化学的博士论文。里面有计算分子电子密度体积V。
此篇博士论文中还有对于V的解释：
我不太清楚这论文中所计算的体积是是否就是分子的体积呢？
您好 谢谢您的回复。
我是看到了一篇有关于计算化学的博士论文。里面有计算分子电子密度体积V。
此篇博士论文中还有对于V的解释：
我不太清楚这论文中所计算的体积是是否就是分子的体积呢？
从描述看应该是的，请看高斯手册中volume关键词的第一句话：
This keyword requests that the molecular volume be computed, defined as the volume inside a contour of 0.001 electrons/Bohr3 density.
谢谢您的回复。
通过您的说明以及又好好看了一遍Gaussian说明书后，还有几个问题想要咨询一下：
（1） 就关键词Volume而言，他是用于计算分子的体积的。这个分子体积的计算环境应该是在非溶剂中的吧？非溶剂环境指的是气态环境还是真空呢？
（2） 这篇博士论文中说明“在上述几何优化和 DFT 的计算中都使用自洽场(SCRF)和积分方程形式极化连续介质模型(IEFPCM)来考虑溶剂(水)效应”。我看过高斯说明书，里面对SCRF关键词的解释是“在溶剂存在的环境下进行相关计算”，可否理解为只是为各种计算加了一个周围环境而已？那么要在水环境之下计算分子的体积V以考虑水的效应，则高斯输入文件的计算语句应该怎么写呢？
是这样么：
%Chk=jobname.chk
#P B3LYP/6-31G(d,p) SCRF=(Solvent=water)
#Volume=(Tight)
这样算是在计算分子体积V的时候考虑到水的效应了么？
为什么要在计算的时候考虑水的效应呢？
一般计算分子体积的时候应该不用加各种溶剂环境吧。
附上这篇博士论文链接：ki.net/kns/detail/detail.aspx?QueryID=0&CurRec=1&recid=&FileName=.nh&DbName=CDFD0911&DbCode=CDFD&pr=
谢谢您的回复。 万分感谢。
2、他的自洽场应该是自洽反应场（self consistent reaction field）的意思，也就是通常所说的溶剂化
输入文件一般这样写即可：
%Chk=jobname.chk
#P B3LYP/6-31G(d,p) SCRF=(Solvent=water) Volume=(Tight) Density
考虑水主要是算结构优化和频率以及热力学数据的。和体积是否有关他并没有描述清楚。当然你可以分别算下气相和水相中的体积来验证，volume我用得不多，无法提供给你经验上的帮助，抱歉
3、文献无法打开
谢谢您的回复。我再去仔细研究一下。谢谢。
北京学而思教育科技有限公司 地址：北京市海淀区北三环甲18号中鼎大厦A座1层102室 电话：010-
浏览器进程
打开微信扫一扫
随时随地聊科研