24小时热门版块排行榜    

查看: 583  |  回复: 0

appletree

铜虫 (小有名气)

[求助] 溶剂化能的计算

比如我想算3-溴苯甲酸(3-BBA)的电化学脱溴电位,通过以下的循环
3-BBA- + 2e- + H2O =BA- + OH- + Br-

3-BBA-gas(# opt freq b3lyp/6-311++g(d,p) pop=nbo)

Zero-point correction=                           0.090990 (Hartree/Particle)
Thermal correction to Energy=                    0.099484
Thermal correction to Enthalpy=                  0.100428
Thermal correction to Gibbs Free Energy=         0.055531
Sum of electronic and zero-point Energies=          -2993.856342
Sum of electronic and thermal Energies=             -2993.847848
Sum of electronic and thermal Enthalpies=           -2993.846903
Sum of electronic and thermal Free Energies=        -2993.891801

3-BBA-solution (# opt freq b3lyp/6-311++g(d,p) scrf=(cpcm,solvent=water,read) pop=nbo)

Zero-point correction=                           0.090932 (Hartree/Particle)
Thermal correction to Energy=                    0.099342
Thermal correction to Enthalpy=                  0.100287
Thermal correction to Gibbs Free Energy=         0.055291
Sum of electronic and zero-point Energies=          -2993.955046
Sum of electronic and thermal Energies=             -2993.946635
Sum of electronic and thermal Enthalpies=           -2993.945691
Sum of electronic and thermal Free Energies=        -2993.990687

溶剂化能=-2993.990687+-2993.891801=-0.098886 au

H20-gas(# opt freq b3lyp/6-311++g(d,p))
Zero-point correction=                           0.021289 (Hartree/Particle)
Thermal correction to Energy=                    0.024125
Thermal correction to Enthalpy=                  0.025069
Thermal correction to Gibbs Free Energy=         0.003647
Sum of electronic and zero-point Energies=            -76.437242
Sum of electronic and thermal Energies=               -76.434406
Sum of electronic and thermal Enthalpies=             -76.433462
Sum of electronic and thermal Free Energies=          -76.454884

For the salvation Gibbs free energy of water (∆G_solv^0(H2O)), we used the standard Gibbs free energy of vaporisation at 298.15 K and 1 atm with a concentration of 1m: ∆G_solv^0(1M)=-6.3 kcal mol-1.


BA-gas(# opt freq b3lyp/6-31++g(d,p))

Zero-point correction=                           0.101650 (Hartree/Particle)
Thermal correction to Energy=                    0.108621
Thermal correction to Enthalpy=                  0.109565
Thermal correction to Gibbs Free Energy=         0.069618
Sum of electronic and zero-point Energies=           -420.201997
Sum of electronic and thermal Energies=              -420.195027
Sum of electronic and thermal Enthalpies=            -420.194083
Sum of electronic and thermal Free Energies=         -420.234030

BA-solution (# opt freq b3lyp/6-311++g(d,p) scrf=(cpcm,solvent=water,read) )
Zero-point correction=                           0.101129 (Hartree/Particle)
Thermal correction to Energy=                    0.108074
Thermal correction to Enthalpy=                  0.109018
Thermal correction to Gibbs Free Energy=         0.068733
Sum of electronic and zero-point Energies=           -420.400986
Sum of electronic and thermal Energies=              -420.394041
Sum of electronic and thermal Enthalpies=            -420.393097
Sum of electronic and thermal Free Energies=         -420.433382

G_solv^0=-420.433382+420.234030 =-0.199352 au

Br-gas(# opt freq b3lyp/6-311++g(d,p))
Zero-point correction=                           0.000000 (Hartree/Particle)
Thermal correction to Energy=                    0.001416
Thermal correction to Enthalpy=                  0.002360
Thermal correction to Gibbs Free Energy=        -0.016176
Sum of electronic and zero-point Energies=          -2574.237796
Sum of electronic and thermal Energies=             -2574.236379
Sum of electronic and thermal Enthalpies=           -2574.235435
Sum of electronic and thermal Free Energies=        -2574.253971

Br-solution(# opt freq b3lyp/6-311++g(d,p) scrf=(cpcm,solvent=water,read))
Zero-point correction=                           0.000000 (Hartree/Particle)
Thermal correction to Energy=                    0.001416
Thermal correction to Enthalpy=                  0.002360
Thermal correction to Gibbs Free Energy=        -0.016176
Sum of electronic and zero-point Energies=          -2574.344613
Sum of electronic and thermal Energies=             -2574.343197
Sum of electronic and thermal Enthalpies=           -2574.342253
Sum of electronic and thermal Free Energies=        -2574.360789

G_solv^0= -2574.360789 +2574.253971=-0.106818

OH-gas(# opt freq b3lyp/6-311++g(d,p))
Zero-point correction=                           0.008535 (Hartree/Particle)
Thermal correction to Energy=                    0.010896
Thermal correction to Enthalpy=                  0.011840
Thermal correction to Gibbs Free Energy=        -0.007721
Sum of electronic and zero-point Energies=            -75.818913
Sum of electronic and thermal Energies=               -75.816552
Sum of electronic and thermal Enthalpies=             -75.815608
Sum of electronic and thermal Free Energies=          -75.835169

OH-solution (# opt freq b3lyp/6-311++g(d,p) scrf=(cpcm,solvent=water,read))
Zero-point correction=                           0.008337 (Hartree/Particle)
Thermal correction to Energy=                    0.010697
Thermal correction to Enthalpy=                  0.011642
Thermal correction to Gibbs Free Energy=        -0.007924
Sum of electronic and zero-point Energies=            -75.956417
Sum of electronic and thermal Energies=               -75.954056
Sum of electronic and thermal Enthalpies=             -75.953112
Sum of electronic and thermal Free Energies=          -75.972678

G_solv^0=-75.972678+75.835169= -0.137509 au

∆G_gas^0=-2574.253971--75.835169--420.234030+76.454884+2993.891801
      =0.023515 au=14.7554 kcal•mol-1

∆G_gas^0 (1M)= ∆G_gas^0 (1 atm)+RTln24.46=16.5924 kcal•mol-1

∆G_solv^0= -0.137509-0.106818-0.199352+0.098886 au+6.3 kcal•mol-1
      =-216.3576+6.3=-210.0576 kcal•mol-1

∆G_aq^0=16.5924-210.0576=-193.4652 kcal•mol-1


E=-(-193.4652)/2/23.061=4.1946 V


请教大神,这样算错误出在哪里,跪求
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子,希望有用哦~

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 appletree 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 Ocx3OKOzGU4q 2026-08-04 5/250 2026-08-04 13:31 by jzdxtiUmJ6PL
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 Ocx3OKOzGU4q 2026-08-03 7/350 2026-08-04 13:22 by jzdxtiUmJ6PL
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 Ocx3OKOzGU4q 2026-08-03 8/400 2026-08-04 13:15 by jzdxtiUmJ6PL
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Ocx3OKOzGU4q 2026-08-03 6/300 2026-08-04 11:34 by MWRYsWbJw4DE
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 Ocx3OKOzGU4q 2026-08-03 7/350 2026-08-04 11:31 by MWRYsWbJw4DE
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 Ocx3OKOzGU4q 2026-08-03 4/200 2026-08-04 11:26 by MWRYsWbJw4DE
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 17/850 2026-08-04 09:46 by archvillain
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 f1dmPx051StA 2026-08-03 7/350 2026-08-04 09:31 by v1TJC4qJauAB
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 09:17 by v1TJC4qJauAB
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 cu9Nq1xK233Z 2026-08-03 10/500 2026-08-04 09:15 by v1TJC4qJauAB
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +12 winnerche 2026-07-29 16/800 2026-08-04 08:18 by Equinoxhua
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 10/500 2026-08-04 08:05 by 88817753
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 07:37 by k7OM8YghWbkC
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 6/300 2026-08-04 06:10 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 04:47 by k7OM8YghWbkC
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
信息提示
请填处理意见