24小时热门版块排行榜    

查看: 4028  |  回复: 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一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 wsYGJnURYvm6 2026-08-31 3/150 2026-09-01 01:52 by o3EMxHUAQDtS
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 EZ67DG6CvcKL 2026-08-31 4/200 2026-09-01 01:30 by o3EMxHUAQDtS
[考博] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 EZ67DG6CvcKL 2026-08-31 4/200 2026-09-01 01:30 by o3EMxHUAQDtS
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 EZ67DG6CvcKL 2026-08-31 3/150 2026-08-31 22:23 by o3EMxHUAQDtS
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
信息提示
请填处理意见