24小时热门版块排行榜    

查看: 4033  |  回复: 9

Swunlly

木虫 (小有名气)

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

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

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

或者大神给些建议吧?


回复此楼

» 猜你喜欢

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

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

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的回帖

Swunlly

木虫 (小有名气)

引用回帖:
Originally posted by yongleli at 2011-04-27 12:36:05:
关于电荷,请google:
partial charge

里边详细讲解了四种电荷方案。
一般来说,NBO和AIM电荷就够用了。
不过这个体系由于有稠环,可能Gaussian不能算出来AIM电荷
(Gaussian的AIM程序是80年代一个东欧 ...

首先,谢谢回复!
实际上,我不是想研究电荷,主要是想对这种离子进行模拟,然后通过看电荷分布的方式来验证,我优化出来的离子确实是N带正电荷。。。
现在,我首先得解决如何进行优化的问题。。。有什么好的建议么?
5楼2011-04-27 13:47:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yongleli

木虫 (正式写手)

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

对不起!
优化的问题没什么更好的办法。。。
实际上我觉得Gaussian的优化需要实验数据。
我以前做的工作都是有晶体结构的。不用怎么在优化上费事。。。
6楼2011-04-27 14:20:36
已阅   回复此楼   关注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的回帖

Swunlly

木虫 (小有名气)

引用回帖:
Originally posted by yongleli at 2011-04-27 14:23:41:
另,友情提醒一下,这个体系好像是大范围共轭结构,
建议基组加“+”。
而且用密度泛函算共轭体系的话,结果可能会有问题,
结构跟实验有差别或者性质不对都是有可能的。建议找实验数据看看。

恩,谢谢帮助!
分少,就给你加1分,表达谢意啦!
8楼2011-04-27 15:21:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yumu0411

金虫 (著名写手)

打酱油

谢谢你们啊!受益匪浅啊!
学习学习再学习
9楼2011-04-27 19:09:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yongleli

木虫 (正式写手)

★ ★
zhou2009(金币+2): 2011-06-25 14:24:24
引用回帖:
Originally posted by yongleli at 2011-04-27 12:36:05:
关于电荷,请google:
partial charge

里边详细讲解了四种电荷方案。
一般来说,NBO和AIM电荷就够用了。
不过这个体系由于有稠环,可能Gaussian不能算出来AIM电荷
(Gaussian的AIM程序是80年代一个东欧 ...

更正一下,CM2是一个半经验电荷,只是基于AM1或者PM3算出来的Mulliken电荷和键级(好像是这个吧?)加上经验参数重新做出来的。可以申请索取他们的MOPAC(这个程序是开源的,很多组都有自己的版本)做。CM2的来源虽然粗糙,但是由于用了参数拟合,做出来的电荷在分子模拟中效果很好。Discovery Studio里估计药物性质用的就是这套电荷。
10楼2011-06-25 11:11:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Swunlly 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-09-02 00:19 by xiangy672
[基金申请] 面上合作单位盖章 +6 ssyjh 2026-08-27 9/450 2026-09-01 20:01 by huagongfeihu
[基金申请] 为什么国自然不能直接公布 +5 bjdxyxy 2026-08-26 5/250 2026-09-01 17:54 by 小伟大博士
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 中青基了要发朋友圈吗? +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
[基金申请] 有没有仍没收到信息的 +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 生科新手
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
信息提示
请填处理意见