24小时热门版块排行榜    

查看: 3855  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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的回帖

孙玉洁

新虫 (初入文坛)



小木虫: 金币+0.5, 给个红包,谢谢回帖
老哥,计算cosmo文件遇到错误,怎么办呢
4楼2018-07-02 17:16:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

showcat

金虫 (小有名气)


不要沉啊!各路大神,鼓励一下初学者吧!
2楼2012-04-18 19:48:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

胡延杰

铁虫 (小有名气)



小木虫: 金币+0.5, 给个红包,谢谢回帖
你好!请问一下解决了吗?
3楼2017-07-07 17:10:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

孙玉洁

新虫 (初入文坛)



小木虫: 金币+0.5, 给个红包,谢谢回帖
请问你计算COSMO文件时,有遇到郭错误的,结果没cosmo文件得情况吗?
5楼2018-07-14 11:18:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0805 316求调剂 +6 大雪深藏 2026-03-18 6/300 2026-03-24 22:13 by peike
[考研] 07化学280分求调剂 +6 722865 2026-03-23 6/300 2026-03-24 20:58 by allen-yin
[考研] 一志愿北化315 求调剂 +3 akrrain 2026-03-24 3/150 2026-03-24 19:35 by 了了了了。。
[考研] 一志愿河北工业大学0817化工278分求调剂 +7 jhybd 2026-03-23 12/600 2026-03-24 09:03 by jhybd
[考研] 一志愿山东大学药学学硕求调剂 +3 开开心心没烦恼 2026-03-23 4/200 2026-03-24 00:06 by 开开心心没烦恼
[考研] 材料/农业专业,07/08开头均可,过线就行 +3 呵唔哦豁 2026-03-23 4/200 2026-03-23 22:30 by 汪!?!
[考研] 接收2026硕士调剂(学硕+专硕) +4 allen-yin 2026-03-23 6/300 2026-03-23 15:04 by 汪!?!
[考研] 276求调剂 +3 YNRYG 2026-03-21 4/200 2026-03-23 08:31 by 醉在风里
[考研] 一志愿华中农业071010,总分320求调剂 +5 困困困困坤坤 2026-03-20 6/300 2026-03-22 17:41 by hxsm
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +3 晨昏线与星海 2026-03-20 3/150 2026-03-22 16:00 by ColorlessPI
[考研] 275求调剂 +6 shansx 2026-03-22 8/400 2026-03-22 15:27 by barlinike
[考研] 生物学调剂 +5 Surekei 2026-03-21 5/250 2026-03-22 14:39 by tcx007
[考研] 一志愿深大,0703化学,总分302,求调剂 +4 七月-七七 2026-03-21 4/200 2026-03-21 18:20 by 学员8dgXkO
[考研] 材料学学硕080502 337求调剂-一志愿华中科技大学 +4 顺顺顺mr 2026-03-18 5/250 2026-03-21 10:22 by luoyongfeng
[考研] 华东师范大学-071000生物学-293分-求调剂 +3 研究生何瑶明 2026-03-18 3/150 2026-03-21 01:30 by JourneyLucky
[考研] 一志愿华中科技大学,080502,354分求调剂 +5 守候夕阳CF 2026-03-18 5/250 2026-03-21 01:06 by JourneyLucky
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 308求调剂 +3 阿姐阿姐家啊 2026-03-18 3/150 2026-03-20 23:24 by JourneyLucky
[考研] 330求调剂 +4 小材化本科 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
信息提示
请填处理意见