| 查看: 4264 | 回复: 5 | |||
[交流]
gaussian 生成cosmo文件
|
|
我在用高斯计算cosmo文件时,看到很多程序都是分三个步骤:1,构型优化;2,PCM计算(radii=klamt);3,生成cosmo文件。我想单纯为了计算cosmo文件的话应该只有第1、3步就好了,为什么非要在生成cosmo之前进行溶剂化计算呢?我猜是不是用溶剂化之后的分子构型生成cosmo文件?所以就在两种情况下分别计算了CO2的cosmo文件,.gjf文件如下: 三步的:%chk=co2.chk %mem=640mw %nproc=8 #p b3lyp/6-31g* opt cosmo 0 1 C 0 0.000000 0.015446 0.000000 O 0 0.000000 -1.237849 0.000000 O 0 0.000000 1.226261 0.000000 --Link1-- %chk=co2.chk %mem=640MW %nproc=8 #P BP86/TZVP scf=(tight,novaracc) SCRF=(CPCM,Read) guess=read geom=checkpoint opt=(MAXCYCLE=99) BP86/TZVP CPCM OPTIMlZATlON WlTH COSMO-RADII 0 1 RADII=KLAMT --Link1-- %chk=co2.chk #P BP86/TZVP scf=(tight,novaracc) SCRF=COSMORS guess=read geom=checkpoint BP86/TZVP COSMO SINGLE POINT 0 1 co2.cosmo 两步的:%chk=co2.chk %mem=64mw %nproc=1 #p b3lyp/6-31g* opt cosmo 0 1 C 0 0.000000 0.015446 0.000000 O 0 0.000000 -1.237849 0.000000 O 0 0.000000 1.226261 0.000000 --Link1-- %chk=co2.chk #P BP86/TZVP scf=(tight,novaracc) SCRF=COSMORS guess=read geom=checkpoint BP86/TZVP COSMO SINGLE POINT 0 1 co2.cosmo 从两者的.out文件中我注意到,第一种情况下CO2键长第一步从1.2533优化到1.1692,第二步从1.1692收敛为1.1721(klamt半径),第三步计算cosmo文件时键长仍然为1.1692;第二种情况下CO2键长第一步从1.2533优化到1.1692,第二步计算cosmo文件时键长为1.1692。也就是说,在生成COSMO文件的时候,利用的分子构型是第一步优化后的构型,与第二步溶剂化无关!? 各路大神,是我想多了,没有弄清计算cosmo文件的过程?还是说单纯计算cosmo文件时确实不需要进行溶剂化计算? |
» 猜你喜欢
一志愿南大理论与计算化学方向(070300)342分,求调剂
已经有1人回复
欢迎化学、化工、物理、材料、或计算机背景的考研咨询
已经有0人回复
物理化学论文润色/翻译怎么收费?
已经有190人回复
黑龙江省原子能研究院化学课题组招收2026级化学专业或化学相关专业硕士生
已经有0人回复
26博士申请
已经有0人回复
河南理工大学化工学院招收调剂生啦!
已经有1人回复
招收博士一名-宁波大学
已经有0人回复
微纳hub,微纳加工的宝库
已经有0人回复
微纳加工里的 lift-off 工艺
已经有0人回复
【接测试】TPD-TPR、TPO、脉冲吸附、原位红外,电化学原位红外,站内联系
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐,对您同样有帮助:
Diamond软件处理CIF文件如何只生成一个分子结构图
已经有12人回复
gaussian计算
已经有4人回复
gaussian 09 做好输入文件后总是关联不上,总说找不到文件,请高手指点!
已经有16人回复
提取Gaussian的rwf文件信息的工具rwfdump简介
已经有13人回复
gaussian09输入文件
已经有8人回复
用gaussian计算丙烷的标准生成热,结果与实验值相差很大,盼高手指导。
已经有7人回复
【sobereva个人文集】自写Link生成Gaussian的IRC任务中每个点的波函数文件
已经有10人回复
【求助】gaussian view 不能打开chk文件
已经有6人回复
【求助】gaussian 无法导入hyperchem和chem3D 的.ent文件
已经有6人回复
【求助】关于生成COSMO时出现的2070错误
已经有6人回复
2楼2012-04-18 19:48:09
3楼2017-07-07 17:10:24
4楼2018-07-02 17:16:16
5楼2018-07-14 11:18:41
6楼2018-07-14 11:23:45










回复此楼