24小时热门版块排行榜    

查看: 4025  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

Swunlly

木虫 (小有名气)

[求助] 能用高斯优化指定某个原子带电荷的离子么?

我想优化1图中的结构,然后用VIEW画图2的结构,设置电荷数为1,自选多重度也为1,这种做法可行么?

算出来的结果也看了看mulliken电荷,但感觉不太对啊。上网搜了搜,有人说mulliken电荷不是很精确,那该怎么判断分析呢?

或者大神给些建议吧?


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

yongleli

木虫 (正式写手)

★ ★
zhou2009(金币+2): 2011-04-27 17:35:02
引用回帖:
Originally posted by Swunlly at 2011-04-27 13:47:31:
首先,谢谢回复!
实际上,我不是想研究电荷,主要是想对这种离子进行模拟,然后通过看电荷分布的方式来验证,我优化出来的离子确实是N带正电荷。。。
现在,我首先得解决如何进行优化的问题。。。有什么好的 ...

另,友情提醒一下,这个体系好像是大范围共轭结构,
建议基组加“+”。
而且用密度泛函算共轭体系的话,结果可能会有问题,
结构跟实验有差别或者性质不对都是有可能的。建议找实验数据看看。
7楼2011-04-27 14:23:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

qzhaosdu

金虫 (著名写手)

什么电荷不太对?是说N上的不是+1吗?
Anewday,anewhour,anewminute,anewpeople.
2楼2011-04-27 11:27:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
gmy1990(金币+3): 2011-04-27 13:22:29
Swunlly(金币+2): 谢谢,先给你2分吧。那我想用高斯算这种离子有没有什么办法呢? 2011-04-27 13:32:59
小红豆(金币+2): 感谢交流 2011-04-27 17:01:34
1、        不能用高斯优化指定某个原子带电荷的离子,某个原子带电荷,分子别处电荷就流过来了;
2、        电荷数为1,自选多重度因为失去了一个电子,应为2;
3、        mulliken电荷在这样的体系不准,要设法按MO的电子密度分别计算出电荷。

[ Last edited by zhou2009 on 2011-4-27 at 17:33 ]
3楼2011-04-27 11:32:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yongleli

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
gmy1990(金币+3): 2011-04-27 13:22:41
Swunlly(金币+1): Thank you! 2011-04-27 15:21:42
小红豆(金币+2): 感谢交流 2011-04-27 17:01:41
Swunlly(金币+1): 2011-05-05 10:42:40
引用回帖:
Originally posted by Swunlly at 2011-04-27 10:28:51:
我想优化1图中的结构,然后用VIEW画图2的结构,设置电荷数为1,自选多重度也为1,这种做法可行么?

算出来的结果也看了看mulliken电荷,但感觉不太对啊。上网搜了搜,有人说mulliken电荷不是很精确,那该怎么判 ...

关于电荷,请google:
partial charge

里边详细讲解了四种电荷方案。
一般来说,NBO和AIM电荷就够用了。
不过这个体系由于有稠环,可能Gaussian不能算出来AIM电荷
(Gaussian的AIM程序是80年代一个东欧人写的,
后来这人走了这部分程序就扔下了没再更新)

Truhlar组里开发的CM2及以后的电荷也很好用,不过可能要自己写小程序了。
4楼2011-04-27 12:36:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +7 gy1nBQXYQJqL 2026-08-29 9/450 2026-08-31 01:15 by hZiFeudZyoGR
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-31 01:07 by hZiFeudZyoGR
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 gy1nBQXYQJqL 2026-08-29 9/450 2026-08-31 01:05 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-31 00:45 by hZiFeudZyoGR
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 11/550 2026-08-31 00:32 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 15/750 2026-08-31 00:10 by hZiFeudZyoGR
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4FFAWE8HcgUD 2026-08-29 4/200 2026-08-30 20:30 by hZiFeudZyoGR
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 4FFAWE8HcgUD 2026-08-29 3/150 2026-08-30 20:19 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 4FFAWE8HcgUD 2026-08-29 3/150 2026-08-30 20:10 by hZiFeudZyoGR
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 19:57 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 5/250 2026-08-30 19:23 by hZiFeudZyoGR
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 29号明天会评吗 +3 笨笨唐 2026-08-28 3/150 2026-08-30 09:29 by cww8181
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 06:57 by ZPa0EcMwuECS
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 05:48 by ZPa0EcMwuECS
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
信息提示
请填处理意见