24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1470  |  回复: 11

linyanglin

木虫 (著名写手)

[求助] 请问有人知道怎样计算配合物基态的氧化电势,急用!如果解决了再送50金币

The rigorous way to obtain the GSOP (also termed ΔGox) is to compute the free
energy difference between the neutral and the oxidized GS species, (G0 - G+)GS.16 The Gibbs free energy in solution of a species i (Gsolv) is defined as Gsolv=Gvac+ ΔGsolv, where Gvac is the Gibbs free energy in gas phase and ΔGsolv
is the free energy of solvation. Gvac is obtained by performing a singlepoint calculation at the optimized geometry in vacuo, followed by frequency calculations in order to include the vibrational contribution to the total partition function. ΔGsolv
i is obtained by a single-point calculation in solution and a reference calculation in gas phase at the geometry optimized in solution.Geometry optimizations of the neutral and oxidized sensitizers were carried out in gas phase and in acetonitrile, that is the solvent experimentally employed for the electrochemical measures,
by using the B3LYP exchange-correlation functional and a 6-31G* basis set. Solvation effects were taken into account by means of the C-PCM solvation model as implemented in Gaussian 03 (G03). We employ the UAHF radii, expressly
optimized to reproduce the experimental solvation energy of a
series of organic molecules. For the calculation of ΔGox we adopted a 6-31+G* basis set, which includes both diffuse and polarization functions, since
in preliminary calculations we found a sizable effect (+0.2-0.3
eV) of the diffuse functions on the computed ΔGox. We tested
the performance of different GGA, hybrid GGA and meta-hybrid
GGA functionals, also combining different exchange and
correlation parts, in order to elucidate the combined effect on
the computed ΔGox of the local and nonlocal HF exchange as
well as of the correlation functional.
每个吉布斯自由能是如何得到的,输入文件的关键词应该怎样设置呢,谢谢。

[ Last edited by linyanglin on 2011-6-10 at 20:11 ]
回复此楼

» 猜你喜欢

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

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

gmy1990

荣誉版主 (著名写手)

优秀版主优秀版主

【答案】应助回帖

linyanglin(金币+5): 谢谢! 2011-06-13 09:37:06
做个优化,执行路径加上freq=noraman(不算拉曼),最后搜索sum,就会有你想要的吉布斯自由能!
2楼2011-06-13 08:42:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

linyanglin

木虫 (著名写手)

引用回帖:
Originally posted by gmy1990 at 2011-06-13 08:42:16:
做个优化,执行路径加上freq=noraman(不算拉曼),最后搜索sum,就会有你想要的吉布斯自由能!

但是上面的文献上应该说的是分别计算改配合物基态和氧化态的吉布斯自由能,而其基态吉布斯自由能是先在真空中优化,然后再计算真空中的单点能=Gvac然后加上该基态在溶剂中优化再分别计算它在气相和溶剂中的单点能,即ΔGsolv,故算出基态的吉布斯自由能。Gsolv=Gvac+ ΔGsolv,但是请问此处的气相指的是该配合物吗?谢谢
3楼2011-06-13 09:16:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gmy1990

荣誉版主 (著名写手)

优秀版主优秀版主

【答案】应助回帖

linyanglin(金币+5): 谢谢! 2011-06-13 10:01:14
引用回帖:
Originally posted by linyanglin at 2011-06-13 09:16:42:
但是上面的文献上应该说的是分别计算改配合物基态和氧化态的吉布斯自由能,而其基态吉布斯自由能是先在真空中优化,然后再计算真空中的单点能=Gvac然后加上该基态在溶剂中优化再分别计算它在气相和溶剂中的单点 ...

高斯中计算所谓气相和溶剂下都是指体系所处环境,所谓的溶剂下,也就是你计算的体系是液相环境,气相也就是真空环境.
4楼2011-06-13 09:49:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

linyanglin

木虫 (著名写手)

引用回帖:
Originally posted by gmy1990 at 2011-06-13 09:49:21:
高斯中计算所谓气相和溶剂下都是指体系所处环境,所谓的溶剂下,也就是你计算的体系是液相环境,气相也就是真空环境.

那么如果以乙醇作为溶剂的话,对基态进行优化,的命令设置为#b3lyp/lanll2dz opt freq scrf=(solvent=ch3ch2oh)这样正确吗?请问如果想要优化该化合物的氧化态,输入文件应怎样设置了?万分感谢!
5楼2011-06-13 10:00:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kaegi

金虫 (小有名气)

【答案】应助回帖

★ ★ ★
gmy1990(金币+3): 非常感谢你对论坛的贡献,奖励一下! 2011-06-13 10:50:41
linyanglin(金币+10): 2011-06-13 11:22:24
引用回帖:
Originally posted by linyanglin at 2011-06-13 10:00:53:
那么如果以乙醇作为溶剂的话,对基态进行优化,的命令设置为#b3lyp/lanll2dz opt freq scrf=(solvent=ch3ch2oh)这样正确吗?请问如果想要优化该化合物的氧化态,输入文件应怎样设置了?万分感谢!

如果用g09的话,应该是#p b3lyp/lanl2dz opt freq scrf=(solvent=Ethanol)
这里有一点疑问,楼主的计算体系里都有些什么元素,如果全部用Lanl2dz不一定合理。
另外如果要计算氧化态的络合物的话,楼主只需要改变计算输入文件中的电荷和自旋多重度就可以了。这两个参数到底是多少就取决于氧化态络合物整体的带电量和未成对电子数的多少了。
6楼2011-06-13 10:47:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

linyanglin

木虫 (著名写手)

引用回帖:
Originally posted by kaegi at 2011-06-13 10:47:39:
如果用g09的话,应该是#p b3lyp/lanl2dz opt freq scrf=(solvent=Ethanol)
这里有一点疑问,楼主的计算体系里都有些什么元素,如果全部用Lanl2dz不一定合理。
另外如果要计算氧化态的络合物的话,楼主只需要 ...

你说的这些我都明白,但具体怎样设置就不太确定了。我计算的体系主要有C、N、O、H、Zn或F等,使用LANL2dz应该可以吧?另外溶剂写入为分子式或英文名称应该是一样的吧!你帮我看看我的分子结构,顺便帮我猜猜它的氧化态的电荷和多重态呗!谢谢啦!呵呵

7楼2011-06-13 11:42:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kaegi

金虫 (小有名气)

引用回帖:
Originally posted by linyanglin at 2011-06-13 11:42:30:
你说的这些我都明白,但具体怎样设置就不太确定了。我计算的体系主要有C、N、O、H、Zn或F等,使用LANL2dz应该可以吧?另外溶剂写入为分子式或英文名称应该是一样的吧!你帮我看看我的分子结构,顺便帮我猜猜它的 ...

建议对C,H,O,N,F用6-31G(d0)基组,对Zn用lanl2dz基组。不知道你这个是几价的Zn。带电荷数和自选多重度是从图上看不出来的啊!再给点提示吧,呵呵。

» 本帖已获得的红花(最新10朵)

8楼2011-06-13 19:08:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

linyanglin

木虫 (著名写手)

送鲜花一朵
引用回帖:
Originally posted by kaegi at 2011-06-13 19:08:23:
建议对C,H,O,N,F用6-31G(d0)基组,对Zn用lanl2dz基组。不知道你这个是几价的Zn。带电荷数和自选多重度是从图上看不出来的啊!再给点提示吧,呵呵。

呵呵,应该是2价的,问了问我导师,他说设置为1 2就行。不过你说的用的是混合基组吗,没有用过,倒是听我师姐提到过。6-31G对Zn应该也是用吧?
9楼2011-06-13 19:59:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kaegi

金虫 (小有名气)

引用回帖:
Originally posted by linyanglin at 2011-06-13 19:59:38:
呵呵,应该是2价的,问了问我导师,他说设置为1 2就行。不过你说的用的是混合基组吗,没有用过,倒是听我师姐提到过。6-31G对Zn应该也是用吧?

你导师说的对,呵呵。我觉得用混合基组好些,Zn用6-31G我觉得有点小了。个人建议而已,不一定就对。
10楼2011-06-13 22:03:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 linyanglin 的主题更新
信息提示
请填处理意见