24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 6168  |  回复: 22
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

loovfnd

超级版主

优秀!!有木有!!!优秀!!有木有!!!优秀!!有木有!!!优秀!!有木有!!!

[交流] 【求助】Molecular electrostatic potential (MEP) maps 分子静电势怎样计算? 已有5人参与

J. AM. CHEM. SOC. 2006, 128, 1402-1403。
如下图所示:

[ Last edited by yjcmwgk on 2009-11-14 at 13:41 ]
回复此楼
莫强求!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

专家顾问

密度泛函·小卒

优秀版主

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
lei0736(金币+2,VIP+0):谢谢 7-18 19:26
引用回帖:
Originally posted by zhou2009 at 2009-7-17 09:00:
谢谢!yjcmwgk

我上面的跟帖,是从手册学习cube和实际操作做的笔记,对各参数作了说明,照贴出来了。真正运行时自然是你推荐的方便、简明。

zhou2009的方法其实对楼主更好。呵呵。
我这样一次写一长串命令,输入快,但是不利于学习
zhou2009大侠,本着诲人不倦的精神,依次列出各个选项的意义,极耐心,比我强多了
9楼2009-07-18 16:19:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 23 个回答

CKX

版主

优秀!!有木有!!!优秀!!有木有!!!优秀!!有木有!!!优秀!!有木有!!!

★ ★ ★ ★ ★ ★ ★ ★
lei0736(金币+3,VIP+0):谢谢 7-10 16:36
loovfnd(金币+5,VIP+0):谢谢你的热心帮助。 7-10 22:11
可以参阅《化学学报》 1992年05期 加入收藏 获取最新 静电势方法研究分子间的电子转移 Ⅰ.理论背景与方法曹轩  廖沐真  顾宪章  李林蜂  吴国是  
【摘要】:光电转换过程中起主导作用的分子间电子转移可视为电子在分子间势垒穿透的过程。我们在从头算SCF-MO波函数基础上将静电势方法用于确定分子间的势垒,并由此出发从理论上得到了各种动力学参数如电子穿透几率和弛豫时间等。文中给出了详细的数学公式和理论模型。
【作者单位】: 清华大学化学系 清华大学化学系 清华大学化学系 清华大学化学系 清华大学化学系
【关键词】: 分子静电势 分子间电子转移 势垒高度 弛豫时间 电子给体 穿透几率 电子穿透 波函数 光电转换 电荷分离
【基金】:863基金资助的课题
【DOI】:CNKI:SUN:HXXB.0.1992-05-004
2楼2009-07-10 14:48:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

新虫

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
lei0736(金币+8,VIP+0):谢谢 7-10 16:37
loovfnd(金币+5,VIP+0):那计算命令呢?输入行要怎么写呢。 7-10 22:11
你可能是问如何作这样的静电势图吧?
我没有作过这样的图。如果要我作这样的图,我会这样去作:

1、        Windows版cubegen的提示信息与回答:
进行g03对图示分子进行计算,得到.chk文件,改写成.fch文件
运行cubegen,
Property [Density]? Potential
……
这样就得到图示分子的分子静电势空间格点文件name.cube

2、        运行gsgrid程序,选择图示分子有静电势等值线图的那个截面:如果图面是xy,即选其为截面,z则选生成该静电势等值线图的那些原子的,它们都在一个平面,只有一个z值,且这个截面之外的原子都不会截取到,没有值。形成gsgrid.txt文件,即时改为自己需要的名字name.txt

3、        进入sigmaplot中作图。
图板的左边小方框图标展示了作图的十多种类型选项,选“Contour Plot”, 出现两个选择:“Contour”,即线条式等值线; “Filled Contour ”, 即色彩填充式等值线。选择:“Contour”,经过几步选择后即形成一个等值线表示的静电势等值线图。进一步选择等值线线型的表达形式,可以找到上面分子图的等值线的色带样式。得色彩等值线图。(这个图看似Filled Contour 的色彩填充式等值线,实应是等值线的彩色线型。)

4、        就在得到的图上面另行拷入一幅分子骨架图,调整好比例,覆盖在上面,二图叠加,并令这上面图层为“透明”,以便能看见下面的等值线图。拖动这骨架图,使之如上面分子图的位置。
这时会又调出等值线图,选择等值线的值,使之符合给定的范畴,这点还要推敲、实践,或者需要剪裁使之适合范畴,上面这个图有明显的剪裁痕迹。因为作者只想看这一圈原子对内的清晰静电势,这也不为过。图即完成了,输出作为插图即可。

当然,我这还是纸上谈兵。关键是这样的图,是作了一个静电势截面等值线图,镶在分子骨架之中。

[ Last edited by zhou2009 on 2009-7-10 at 19:22 ]

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

3楼2009-07-10 15:55:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

禁虫

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
lei0736(金币+6,VIP+0):谢谢 7-11 10:13
loovfnd(金币+6,VIP+0):太感谢了! 7-11 12:37
用G03作单点能计算,此时不需要用作CUBE的命令,用%chk=jobname得到CHK文件即成。再调用cubegen命令作分子的静电势空间格点文件。

用G03作单点能计算,也可用关键字cube=(potential,cards),自己输入格点数、起始点、格点间隔。即得分子的静电势空间格点文件。

用cubegen命令

1、在计算中,用%chk=jobname----------将计算存为:jobname.chk;默认存于:scratch子文件夹。

2、调用Formchk 将二进制的checkpoint 文件(chk)转换成ASCII 码形式(fchk);得到jobname.fch

3、调用cubegen。

    在一个三维网格(立方)点上求解:Orbitals分子轨道( MO=n),potential静电势,Density电子密度,Gradient密度梯度,NormGradient密度梯度的范数和Laplacian电子密度的拉普拉斯算符(∇2ρ)。
默认计算电子密度(对应选项Density)。使用的密度由Density 关键字控制;若使用Density=Current 求解来自相关的或CI-Singles 波函的立方区域内电子密度,而不是默认的Hartree-Fock 密度。

Windows版cubegen的提示信息与回答:

Property [Density]? MO=n ----------指定要选取的MO数。一旦回答了MO=n,就不是作Density而是作Orbitals。分子轨道数n、关键字Homo,Lumo,All,OccA(全部是alpha占据),OccB(全部是beta 占据),Valence(所有价轨道) 和  Virtuals(所有虚轨道)也可以用来代替特定的轨道数。n 值没有默认值,省略的话会造成错误。

Formatted Checkpoint file? jobname ----------输入的fchk文件的名字,写路径,但不写后缀,  后缀默认为.fch。

Cube file [ ]? Jobname             ----------输出的cube文件的名字,写路径,但不写后缀,  后缀默认为.cube。

Approximate points per side [0]? 回车 ------点阵每边的点数。回车取0为默认的80×80×80。

Header in cube file [H]?  回车     ---------输出档案的格式:H表示包含标题(这是预设选择); GsGrid要求H。回答n,表示不要包含标题。若要产生无格式方块文件,这选项会被忽略。

如果是:
Property [Density]? Potential
得到的即是分子的静电势空间格点文件。
4楼2009-07-11 10:02:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 考研生物与医药调剂 +7 铁憨憨123425 2026-03-31 7/350 2026-04-01 08:45 by JourneyLucky
[考研] 359求调剂 +8 邓邓邓书书 2026-03-25 8/400 2026-04-01 08:20 by herarysara
[考研] 0703一志愿南师大334求调剂 +3 seven7yu 2026-03-30 3/150 2026-04-01 08:08 by zjbkx
[考研] 070300化学专业279调剂 +10 哈哈哈^_^ 2026-03-31 10/500 2026-03-31 23:13 by liu823948201
[考研] 086000生物与医药298调剂求助 +3 元元青青 2026-03-31 5/250 2026-03-31 19:58 by vgtyfty
[考研] 材料工程085601数二英一335求调剂 +5 双马尾痞老板2 2026-03-31 5/250 2026-03-31 19:07 by Wang200018
[考研] 0856 335分 | 封装or激光加工 老师看看我!! +9 cccchenso 2026-03-29 9/450 2026-03-31 16:37 by lishahe
[考研] 085600 295分求调剂 +13 W55j 2026-03-30 15/750 2026-03-31 13:29 by 王亮_大连医科大
[考研] 一志愿西电085401数一英一299求调剂 六级521 +4 爱吃大鸭梨 2026-03-31 4/200 2026-03-31 11:51 by 搏击518
[考研] 南京大学化学调剂 +11 景随风 2026-03-29 16/800 2026-03-31 10:14 by herarysara
[考研] 293分求调剂,外语为俄语 +5 加一一九 2026-03-31 5/250 2026-03-31 09:39 by zhshch
[考研] 085601一志愿西北工业大学初试346 +4 085601初试346 2026-03-30 4/200 2026-03-31 07:47 by jp9609
[考研] 293求调剂 +3 末未mm 2026-03-30 5/250 2026-03-30 17:23 by 王保杰33
[考研] 329求调剂 +8 星野? 2026-03-26 8/400 2026-03-30 13:41 by chemdavid
[考研] 材料与化工304求B区调剂 +4 邱gl 2026-03-26 7/350 2026-03-30 08:39 by 探123
[考研] 343求调剂085601 +3 要努力学习x 2026-03-29 3/150 2026-03-29 18:35 by wxiongid
[考研] 283求调剂 +3 A child 2026-03-28 3/150 2026-03-28 15:41 by ms629
[考研] 一志愿211院校 344分 东北农业大学生物学学硕,求调剂 +5 丶风雪夜归人丶 2026-03-26 8/400 2026-03-27 19:22 by 丶风雪夜归人丶
[考研] 调剂求收留 +7 果然有我 2026-03-26 7/350 2026-03-27 00:26 by wxiongid
[考研] 调剂 +4 柚柚yoyo 2026-03-26 4/200 2026-03-26 20:43 by fmesaito
信息提示
请填处理意见