24小时热门版块排行榜    

查看: 253  |  回复: 2
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

qzhaosdu

金虫 (著名写手)

[交流] 【求助】十二烷基二甲基氧化胺--电子结构

最近想做十二烷基二甲基氧化胺的一些分子模拟,主要想看看N-O和金属离子的作用。可是一上来就遇到了麻烦,用MS画了分子构型,简单做了下优化,可是用MS自带的2种粗略算电荷的方法得出的结果确不一样,一个说是0e,一个说是1e.
可是我感觉应该是不带电的,所以按照不带电的用gauss又优化了下,就想拿这个作为初始结构用gromacs去做下模拟,放了20个分子,但是最后grompp以后程序说还带着一些负电荷,不大,0.00626e。---这个.itp是通过prodrg生成的,当时prodrg就在氧上给我加了个H上去,我给去掉了,但是所有的电荷我都用gauss出来的密丽根电荷重新加的。
这样去做em就出错了。。。然后我去gromacs的mailinglist去看了下,好多这个错误的,用了几种方法去修改,还是不对。。。

于是我就想用discover是不是能跑下,可是无论用什么力场,都没有N-O这个,程序不给做。。。

请教做过这个N-O的高人指点。。。

[ Last edited by lei0736 on 2009-11-24 at 19:57 ]
回复此楼

» 猜你喜欢

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

Anewday,anewhour,anewminute,anewpeople.
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qiaobf

★ ★
lei0736(金币+2,VIP+0):鼓励新人
没太看懂:你的.itp的电荷总和是0,然后grompp认为不是0?这不可能呀!还有即使grompp认为总电荷不是0,也就是一个warnning,但是还是能跑的。
另外用MS计算电荷的时候,应该是指定总电荷,然后再分配电荷的吧
3楼2008-12-15 11:26:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

yyx19840628

木虫 (著名写手)

鹰的眼睛,狼的耳朵,豹的速 ...

没有做过,要不你用pop=nbo算算看再?
人生的最大遗憾莫过于错误地坚持了不该坚持的,轻易地放弃了不该放弃的……
2楼2008-12-13 18:37:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见