24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 681  |  回复: 2
当前主题已经存档。

penghcp

木虫 (小有名气)

可爱小猫咪

[交流] 【求助】计算溶剂效应(pcm模型)能量问题,高手请进

各位大虾,小弟在用pcm计算时发现输出文件少了一项能量,具体情况如下:
关键字:#P B3LYP/6-31G(d) scf=tight int=finegrid SCRF=(PCM,Read,Solvent=Water)

输出文件关于溶剂化能量如下


Variational PCM results
=======================
                     (a.u.) =    -155.029051
                     (a.u.) =    -155.039159
Total free energy in solution:
  with all non electrostatic terms            (a.u.) =    -155.030980
--------------------------------------------------------------------
(Polarized solute)-Solvent               (kcal/mol) =      -6.34
--------------------------------------------------------------------
Cavitation energy                        (kcal/mol) =      11.83
Dispersion energy                        (kcal/mol) =      -7.15
Repulsion energy                         (kcal/mol) =       0.45
Total non electrostatic                  (kcal/mol) =       5.13
--------------------------------------------------------------------

当我查看高斯使用指南时发现他例子中的能量输出比我多了一项(Total Electrostatic 静电项),他书中的例子
-

------------- VARIATIONAL PCM RESULTS -------------
------------------------------------------------------
(a.u.) = -100.024068
(a.u.) = -100.028947
(a.u.) = -100.029186
(a.u.) = -100.023819
(a.u.) = -100.029187
Total free energy in sol.
(with non electrost.terms) (a.u.) = -100.028994
------------------------------------------------------
(Unpol.Solute)-Solvent (kcal/mol) = -3.06 静电项(包括极性
(Polar.Solute)-Solvent (kcal/mol) = -3.37 和非极性作用)
Solute Polarization (kcal/mol) = 0.16
Total Electrostatic (kcal/mol) = -3.21
------------------------------------------------------
Cavitation energy (kcal/mol) = 3.20 非静电项(包括
Dispersion energy (kcal/mol) = -4.12 空腔能量、分散
Repulsion energy (kcal/mol) = 1.04 能和排斥能)
Total non electr. (kcal/mol) = 0.12
------------------------------------------------------

但是我又查了中文手册发现没这一项 输出如下
SCF Done: E(RHF) = -100.029187240 A.U. after 5 cycles
Convg = 0.4249D-05 -V/T = 2.0033
S**2 = 0.0000
--------------------------------------------------------------------
Variational PCM results
=======================
(a.u.) = -98.568013
(a.u.) = -98.573228
Total free energy in solution:
with all non electrostatic terms (a.u.) = -98.569083
--------------------------------------------------------------------
(Polarized solute)-Solvent (kcal/mol) = -3.27
--------------------------------------------------------------------
Cavitation energy (kcal/mol) = 5.34
Dispersion energy (kcal/mol) = -3.08
Repulsion energy (kcal/mol) = 0.34
Total non electrostatic (kcal/mol) = 2.60
--------------------------------------------------------------------

这是怎么回事啊,计算溶剂效应时应该包含静电项能量么,是不是我的关键字输入有误,还是其他,敬请高手解答
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

S07111072

铁杆木虫 (著名写手)

★ ★
penghcp(金币+1,VIP+0):这个我知道,但是问题还是没解决,问什么我的输出文件没有这一项 9-24 12:32
yjcmwgk(金币+1,VIP+0): 9-24 13:38
Variational PCM results              
=======================              
            
            
            
            
            
Total free energy in solution:      
with all non electrostatic terms   
-------------------------------------
(Unpolarized solute)-Solvent         
(Polarized solute)-Solvent           
Solute polarization                  
Total electrostatic                  
-------------------------------------
Cavitation energy                    
Dispersion energy                    
Repulsion energy                     
Total non electrostatic              
DeltaG (solv)                        
包含静电项能量,一般PCM计算输出包含上面的能量项
2楼2009-09-24 09:28:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

S07111072

铁杆木虫 (著名写手)


yjcmwgk(金币+1,VIP+0): 9-25 09:16
引用回帖:
Originally posted by penghcp at 2009-9-23 21:27:
各位大虾,小弟在用pcm计算时发现输出文件少了一项能量,具体情况如下:
关键字:#P B3LYP/6-31G(d) scf=tight int=finegrid SCRF=(PCM,Read,Solvent=Water)

输出文件关于溶剂化能量如下


Variational PC ...

在你输入文件最后加上

PCMDOC      
RADII=UAHF  
SCFVAC
3楼2009-09-25 09:15:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 penghcp 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 304求调剂(085602,过四级,一志愿985) +13 化工人999 2026-04-04 13/650 2026-04-05 13:42 by ujn_zhuj
[考研] 一志愿南航,数一英一学硕317求调剂!! +5 Acaciad 2026-04-04 5/250 2026-04-05 12:31 by 搏击518
[考研] 081700学硕,323分,一志愿中国海洋大学求调剂学校 +16 披星河 2026-04-04 16/800 2026-04-05 11:27 by 猪会飞
[有机交流] 甲醇/二氯 1:15过柱子 5+3 a哎y呦w喂 2026-03-31 3/150 2026-04-05 10:42 by 88817753
[考研] 一志愿江南大学085501机械工程专硕326分,本科佳木斯大学 +5 顾若浮生 2026-04-03 9/450 2026-04-05 09:57 by 1753564080
[考研] 11408,335分,本科211,求调剂,可转专业 +5 鳄梨大鳄鱼 2026-04-03 5/250 2026-04-04 22:49 by chongya
[考研] 086000生物与医药298调剂求助 +6 元元青青 2026-03-31 8/400 2026-04-04 22:15 by lys0704
[考研] 333求调剂 +12 wfh030413@ 2026-04-03 13/650 2026-04-04 21:02 by jj987
[考研] 22408求调剂 354分 可跨专业 +3 hannnnnnn 2026-04-04 3/150 2026-04-04 14:35 by 土木硕士招生
[考研] 350一志愿北京航空航天大学08500材料科学与工程求调剂 +5 kjnasfss 2026-04-03 5/250 2026-04-03 22:29 by 无际的草原
[考研] 322求调剂 +4 FZAC123 2026-04-03 4/200 2026-04-03 20:55 by zhq0425
[考研] 329求调剂,一志愿西北工业大学,材料工程(085601) +8 小小机灵虫 2026-03-29 14/700 2026-04-03 19:38 by lijunpoly
[考研] 数一英一285求调剂 +7 AZMK 2026-04-03 9/450 2026-04-03 13:03 by ms629
[考研] 316求调剂 +14 舟自梗 2026-04-01 18/900 2026-04-03 10:28 by linyelide
[考研] 274求调剂 +10 薛定谔的虎。 2026-04-01 10/500 2026-04-03 10:13 by tianyyysss
[考研] 309求调剂 +14 呆菇不是戴夫 2026-04-02 14/700 2026-04-03 09:42 by 蓝云思雨
[考研] 求调剂!生物与医药专硕 +4 逆转陆先生 2026-04-01 5/250 2026-04-03 08:33 by Jaylen.
[考研] 化学070300-总分378-求调剂 +5 挪椅子的泡泡糖 2026-04-02 5/250 2026-04-02 22:20 by ZXlzxl0425
[考研] 353求调剂 +4 拉钩不许变 2026-04-01 4/200 2026-04-01 18:10 by 记事本2026
[考研] 080200学硕,机械工程专业277分,求带走! +4 瓶子PZ 2026-03-31 4/200 2026-03-31 20:16 by vgtyfty
信息提示
请填处理意见