24小时热门版块排行榜    

查看: 3370  |  回复: 5
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 showcat 的 2 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

showcat

金虫 (小有名气)


[交流] 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文件时确实不需要进行溶剂化计算?
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

showcat

金虫 (小有名气)


不要沉啊!各路大神,鼓励一下初学者吧!
2楼2012-04-18 19:48:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 showcat 的主题更新
普通表情 高级回复(可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 论文投稿 5+3 love319 2024-09-20 5/250 2024-09-21 07:30 by nono2009
[论文投稿] PCCP投稿 +4 wuxiaohuang1 2024-09-20 4/200 2024-09-21 01:09 by 褚师宣密
[硕博家园] 现博士在读,硕士学校的硕导又让干活 +7 也就这样 2024-09-18 7/350 2024-09-20 21:53 by 望惜时W
[论文投稿] 应用设计类论文投稿求助 10+3 chenweigang 2024-09-16 8/400 2024-09-20 21:45 by abm456
[基金申请] 有感于请专家指导申报书 +8 huitong441 2024-09-16 14/700 2024-09-20 21:32 by huitong441
[论文投稿] 文章改了几轮,该投稿了,选了几篇杂志发给导师 +5 evan0489 2024-09-19 6/300 2024-09-20 20:13 by lyd782
[考博] 2025申博 +3 yqrxxxl 2024-09-18 4/200 2024-09-20 19:49 by ZHI1009
[有机交流] 核磁解谱 50+4 luffyyyyy 2024-09-16 15/750 2024-09-20 18:20 by luffyyyyy
[硕博家园] 关于博士选择 +3 翰墨古风 2024-09-20 6/300 2024-09-20 17:25 by 大泡沫
[找工作] 博士毕业16,两中两民; +5 hocoltyaya 2024-09-17 8/400 2024-09-20 16:51 by 93479180
[考博] Aggregate +6 学员eIXAWl 2024-09-15 8/400 2024-09-20 12:46 by yifeng11
[基金申请] 结项4年了,论文还能挂基金号吗 +10 xiandaiqi 2024-09-15 12/600 2024-09-19 22:00 by yuch
[论文投稿] 顶级水刊,几乎不拒稿的SCI! +4 a433587770 2024-09-18 4/200 2024-09-19 21:50 by ndyjjcb
[论文投稿] 作者署名:nature大子刊共一还是共同通讯? 10+7 Eriknpeng 2024-09-16 18/900 2024-09-19 17:32 by Eriknpeng
[论文投稿] 有没有便宜一点的开源期刊 100+3 zhuzhuyaya 2024-09-19 4/200 2024-09-19 16:43 by 莱茵润色
[基金申请] 有教育部人文项目在研,可以申报国家社科吗? +6 headtop 2024-09-15 6/300 2024-09-19 15:30 by 求易
[基金申请] 立项公示后,已经有几个老师及同学过来要本子,该怎么应对? +27 ty28ty28 2024-09-15 33/1650 2024-09-19 14:36 by sophiemyn
[有机交流] 酰胺缩合反应 5+3 我努力三百 2024-09-14 3/150 2024-09-19 10:04 by czyzsu
[论文投稿] 期刊主题不符是真的主题不行吗还是内容不行 5+6 AlexaHu 2024-09-15 11/550 2024-09-19 09:40 by AlexaHu
[基金申请] 最近这些情况,大家怎么看 +11 wjynus 2024-09-14 13/650 2024-09-16 14:05 by chemhua
信息提示
请填处理意见