| 查看: 1252 | 回复: 7 | ||
[求助]
GAMESS中几何优化和频率计算同时做应如何输入
|
|
如题，想在GAMESS中进行几何优化后马上做频率计算，看是否有虚频，就像gaussian中一样。请问输入文件应该怎么写？用MacMolPlt可以方便地构建几何优化的input file，那么频率计算如何输入？是否是runtyp=optimize后面加上runtyp=hessian，METHOD=ANALYTIC VIBANL=.TRUE.这些？ 用gaussian做freq可以得到化合物的零点能，焓，吉布斯自由能等数据，gamess中的freq是否可以？ |
回复此楼» 猜你喜欢
N-甲基苯胺烷基化
已经有14人回复
-
江西理工大学招收调剂
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有52人回复
-
上海工程技术大学招收调剂
已经有1人回复
-
江西科技师范大学硕士招生
已经有0人回复
-
与dbu相似结构的碱有哪些？
已经有2人回复
-
车辆美容镀膜镀晶光亮保护含硅类原料介绍（二）
已经有0人回复
-
课题组新增一中科院联培博士指标，材料设计，水凝胶 组织工程
已经有0人回复
-
想问一下各位大佬，这种光波导的图是怎么测出来的呀？
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 热工设备的仿真与优化
已经有9人回复
- 如何取出单晶金属配合物？
已经有7人回复
- castep计算完能带之后如何挑选有用的数据画出清晰的能带图
已经有16人回复
- 非稳态计算时time step和max iteration number的影响
已经有10人回复
- 求助！哪位高人帮忙优化一下分子结构并画出前线轨道模型~~~（重谢~）
已经有5人回复
- seebeck系数计算的一些问题
已经有15人回复
- VASP结构优化时出现的POTIM问题
已经有4人回复
- 弹性常数计算出现诡异。。。。来看看喽
已经有5人回复
- 纯理论计算投稿的问题
已经有21人回复
- 结构优化 a=b ----> a不等于b
已经有3人回复
- 几何优化 收敛问题
已经有6人回复
- 急急急：GC/MS条件优化
已经有15人回复
- 表面分子吸附的几何优化
已经有5人回复
- CASTEP几何优化计算中断，大侠们帮帮忙啊！
已经有5人回复
- 【求助】castep做几何优化时其参数如何设置？
已经有19人回复
yflchx
金虫 (正式写手)
- 应助: 19 (小学生)
- 金币: 1029.6
- 散金: 3140
- 红花: 14
- 帖子: 568
- 在线: 3042.2小时
- 虫号: 766553
- 注册: 2009-05-09
- 性别: GG
- 专业: 理论和计算化学
【答案】应助回帖
jianqin(金币+10): THX 2011-08-31 19:47:16
|
$CONTRL SCFTYP=RHF RUNTYP=OPTIMIZE DFTTYP=wB97X-D MAXIT=200 ICHARG=-1 MULT=1 MOLPLT=.TRUE. $END $SYSTEM TIMLIM=525600 MEMORY=10000000 $END $BASIS GBASIS=N31 NGAUSS=6 NDFUNC=1 NPFUNC=1 DIFFSP=.TRUE. $END $SCF NCONV=9 $END $STATPT OPTTOL=0.0001 NSTEP=100 HSSEND=.t. $END $DATA 分子结构省略。 |
2楼2011-08-30 22:27:38
3楼2011-08-31 19:42:38
yflchx
金虫 (正式写手)
- 应助: 19 (小学生)
- 金币: 1029.6
- 散金: 3140
- 红花: 14
- 帖子: 568
- 在线: 3042.2小时
- 虫号: 766553
- 注册: 2009-05-09
- 性别: GG
- 专业: 理论和计算化学
4楼2011-08-31 21:54:53
5楼2011-08-31 22:18:57
yflchx
金虫 (正式写手)
- 应助: 19 (小学生)
- 金币: 1029.6
- 散金: 3140
- 红花: 14
- 帖子: 568
- 在线: 3042.2小时
- 虫号: 766553
- 注册: 2009-05-09
- 性别: GG
- 专业: 理论和计算化学
|
$STATPT OPTTOL=0.0001 NSTEP=100 HESS=CALC HSSEND=.t. $END $FORCE METHOD=SEMINUM VIBSIZ=0.010000 VIBANL=.TRUE. $END 使用HSSEND=.t.之后，在输出文件的末尾就会有下面的输出结果。 挑选你想要的就可以了。 ------------------------------- THERMOCHEMISTRY AT T= 298.15 K ------------------------------- USING IDEAL GAS, RIGID ROTOR, HARMONIC NORMAL MODE APPROXIMATIONS. P= 1.01325E+05 PASCAL. ALL FREQUENCIES ARE SCALED BY 1.00000 THE MOMENTS OF INERTIA ARE (IN AMU*BOHR**2) 8.48840 466.28629 466.98211 THE ROTATIONAL SYMMETRY NUMBER IS 2.0 THE ROTATIONAL CONSTANTS ARE (IN GHZ) 212.41782 3.86691 3.86115 1 IMAGINARY FREQUENCY VIBRATION(S) IS(ARE) EXCLUDED FROM THE FOLLOWING THERMOCHEMICAL ANALYSIS. THE HARMONIC ZERO POINT ENERGY IS (SCALED BY 1.000) 0.028683 HARTREE/MOLECULE 6295.280937 CM**-1/MOLECULE 17.999129 KCAL/MOL 75.308358 KJ/MOL Q LN Q ELEC. 1.00000E+00 0.000000 TRANS. 1.56038E+07 16.563028 ROT. 7.69486E+03 8.948308 VIB. 1.95841E+00 0.672134 TOT. 2.35145E+11 26.183470 E H G CV CP S KJ/MOL KJ/MOL KJ/MOL J/MOL-K J/MOL-K J/MOL-K ELEC. 0.000 0.000 0.000 0.000 0.000 0.000 TRANS. 3.718 6.197 -41.059 12.472 20.786 158.498 ROT. 3.718 3.718 -22.182 12.472 12.472 86.872 VIB. 78.383 78.383 73.642 22.331 22.331 15.902 TOTAL 85.820 88.299 10.401 47.275 55.589 261.271 VIB. THERMAL CORRECTION E(T)-E(0) = H(T)-H(0) = 3074.972 J/MOL E H G CV CP S KCAL/MOL KCAL/MOL KCAL/MOL CAL/MOL-K CAL/MOL-K CAL/MOL-K ELEC. 0.000 0.000 0.000 0.000 0.000 0.000 TRANS. 0.889 1.481 -9.813 2.981 4.968 37.882 ROT. 0.889 0.889 -5.302 2.981 2.981 20.763 VIB. 18.734 18.734 17.601 5.337 5.337 3.801 TOTAL 20.512 21.104 2.486 11.299 13.286 62.445 VIB. THERMAL CORRECTION E(T)-E(0) = H(T)-H(0) = 734.936 CAL/MOL ......END OF NORMAL COORDINATE ANALYSIS...... |
6楼2011-08-31 22:24:43
7楼2019-03-18 23:27:38
|
input 文件如下： $CONTRL SCFTYP=UHF RUNTYP=OPTMIZE EXETYP=routine MAXIT=200 $END $SYSTEM TIMLIM=525600 MEMORY=100000000 $END $ BASIS GBASIS=N21 NGAUSS=3 $END $SCF DIRSCF=.TRUR. CONV=1.0d-06 RESET=.TRUE. $END $STATPT METHOD=SCHLEGEL OPTTOL=0.000005 NSTEP=100 IFOLOW=-1 RMAX=0.8 $END $ DATA 1st row is reserved for a comment C1 O 8.0 0.086 0.095 0.00 H 1.0 0.740 1.242 0.00 H 1.0 -1.07 0.430 0.00 $END |
8楼2019-03-18 23:36:19
10