24小时热门版块排行榜    

查看: 6589  |  回复: 10

zxwasd

金虫 (小有名气)

[交流] 【求助】Fukui指数计算问题? 已有6人参与

大家好!我想知道如何计算一个分子中相应原子的Fukui指数。例如我想计算甘氨酸分子(H2NH2CCOOH)中每个原子的fukui指数,我该怎么做:
(1)先优化和计算H2NH2CCOOH中每个原子的mulliken电荷分布;然后分别优化和计算H3N(+)H2CCOOH和H2NH2CCOO-的mulliken电荷分布,最后根据公式:
f+(r)=q(N+1)-q(N)和f-(r)=q(N)-q(N-1)计算的fukui指数。
(2)优化和计算H2NH2CCOOH,但在电荷输入时,分别输入0,+1和-1。
希望知道的大侠不惜赐教,先深表感谢,谢谢了!!!
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

端木青

金虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by sobereva at 2010-04-14 23:47:09:
计算福井函数(包括用布居方法计算的紧缩福井函数),只需对中性情况优化几何结构即可,带电情况下不应优化,应保持中性时的结构。因为福井函数的rho对N求导带着外势不变的条件,外势即原子核位置的函数。所以用差 ...

关于fukui指数的求算,我的理解是这样的两种情况
1.优化好中性分子后得到电子q(N),q(N+1)在这个中性分子基础上直接把电荷由零加1.
2.优化好中性分子后得到电子q(N),q(N+1)是在将中性分子修改成离子构型后优化计算得到
请问这两种哪种是合理正确的呢?
我做的是带有两个羧基的分子。。。
请大虫不吝赐教啊!
10楼2010-12-31 09:38:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

gnli

铜虫 (小有名气)

啥是fuku指数啊
争取早日成为合格的菜鸟
2楼2010-01-26 20:38:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

G03W

银虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
你自己不是很懂得啊,还问啊
3楼2010-01-27 09:02:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

CKX

木虫 (著名写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
fuku指数可以查文献的啊,文献中有说明的
4楼2010-01-27 09:07:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主


小木虫(金币+0.5):给个红包,谢谢回帖交流
你写的挺对的啊
你还想问啥
5楼2010-01-27 21:30:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuxianlv

金虫 (正式写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
aylayl08(金币+2):呵呵谢谢回帖交流 2010-04-20 18:14
我在练习《用量子化学密度泛函理论研究环状含氮化合物分子结构与缓蚀性能的关系》第24卷第4期,中国腐蚀与防护学报,这篇文献时有疑惑,按照文献“富克指数的计算是,对中性分子和分别带一个正,负电荷的体系在几何结构全优化的基础上,进行了Mullike布局分析计算”我对这话理解有困难,我是用Gaussview建构型先对中性分子用Gaussion03优化并算频率,若无虚频(得原子电荷M)再将电荷数改为1或-1优化带电分子的构型(得M+或M-),从而可以算富克指数。我的困惑是用Gaussview建构型先对中性分子用Gaussion03优化并算频率,若无虚频(得原子电荷M)后算带电分子的原子电荷时要不要优化构型?即是直接改电荷数和自旋度做sp计算呢还是要改电荷数和自旋度后做opt计算啊?
           以上两种情况我都试了,分子中性时M值与原文符合,用改电荷数和自旋度后做opt计算时M+值和M值与原文都符合,就是M-始终不符合;中性分子用Gaussion03优化并算频率,若无虚频(得原子电荷M)后直接改电荷数和自旋度做sp计算时M+和M-都对不上。
       所以我想学习计算富克指数的细节 设置。
    另外为什么初始结构略微变化后经相同基组优化后,键长,键角等略微不同啊,但两次的收敛都是4个Yes频率计算也无虚频啊,应该找到的就是体系能量最低点啊,最低点应该只有一个啊,别人的高斯计算文献数据为什么采用相同基组却很难重复呢?
      我翘首以盼各位的指点。
6楼2010-04-14 22:59:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫


★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
yjcmwgk(金币+3): 2010-04-15 13:26
计算福井函数(包括用布居方法计算的紧缩福井函数),只需对中性情况优化几何结构即可,带电情况下不应优化,应保持中性时的结构。因为福井函数的rho对N求导带着外势不变的条件,外势即原子核位置的函数。所以用差分近似来计算福井函数时前后结构应保持一致。

优化时为获得精确结果应设严收敛限。
7楼2010-04-14 23:47:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qzhaosdu

金虫 (著名写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by sobereva at 2010-04-14 23:47:09:
计算福井函数(包括用布居方法计算的紧缩福井函数),只需对中性情况优化几何结构即可,带电情况下不应优化,应保持中性时的结构。因为福井函数的rho对N求导带着外势不变的条件,外势即原子核位置的函数。所以用差 ...

膜拜S大,您简直就是爱因斯坦再世,薛定谔下凡~~
Anewday,anewhour,anewminute,anewpeople.
8楼2010-04-16 11:08:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zxwasd

金虫 (小有名气)

★ ★
aylayl08(金币+2):谢谢回帖交流 2010-04-20 18:14
引用回帖:
Originally posted by liuxianlv at 2010-04-14 22:59:09:
我在练习《用量子化学密度泛函理论研究环状含氮化合物分子结构与缓蚀性能的关系》第24卷第4期,中国腐蚀与防护学报,这篇文献时有疑惑,按照文献“富克指数的计算是,对中性分子和分别带一个正,负电荷的体系在几何 ...

呵呵,你好!我也遇到了类似的问题!
      我知道的情况是这样的:高斯优化的结果与输入的初始结构密切相关,虽然你优化得到了所谓的最低点,但只是相对于你输入的结构,而不一定是分子的最低点!由于输入的初始结构可能不同,所以,别人的高斯计算文献数据采用相同基组却很难重复。
      判断自己是否算正确,可以比较一下自己的能量是否比文献更低,而且在量子化学研究者眼中,中国腐蚀与防护学报肯定没有分量,所以正确与否,还需要自己判断!
9楼2010-04-18 17:38:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zxwasd 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 mibUvS8DDCwf 2026-09-15 3/150 2026-09-18 02:54 by cITaGg3p5edr
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 vZfe6xYu34yj 2026-09-14 4/200 2026-09-18 02:28 by cITaGg3p5edr
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-18 02:16 by cITaGg3p5edr
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 Aj1rhIDL5ixY 2026-09-14 6/300 2026-09-18 01:54 by cITaGg3p5edr
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-18 01:04 by cITaGg3p5edr
[考博] 申博发邮件 +6 Lmengk 2026-09-15 7/350 2026-09-17 23:36 by Lmengk
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-13 8/400 2026-09-17 23:28 by cuPZTDXS3VOd
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-17 22:39 by cuPZTDXS3VOd
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 Aj1rhIDL5ixY 2026-09-14 6/300 2026-09-17 14:17 by ObVzqQQh1rrE
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-17 05:38 by k06BKNblrGRH
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 Aj1rhIDL5ixY 2026-09-14 4/200 2026-09-17 04:52 by k06BKNblrGRH
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-17 03:16 by Tql5LQhh5rLK
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-17 02:38 by Tql5LQhh5rLK
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 s3fFTmArrBt6 2026-09-13 7/350 2026-09-17 02:14 by Tql5LQhh5rLK
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +7 3n8v2C8RimXI 2026-09-13 8/400 2026-09-17 01:38 by Tql5LQhh5rLK
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 QUjhNVAcOSff 2026-09-13 6/300 2026-09-16 14:27 by i3KBRD19ERkG
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-16 02:14 by DgNGHc3h5tPl
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 17:39 by CgyNCDVNhGVg
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:23 by CgyNCDVNhGVg
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 07:02 by 5BDX0d0WFp7t
信息提示
请填处理意见