24小时热门版块排行榜    

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

shangshui123

新虫 (职业作家)

[求助] DMOL3求吸附能遇到问题 已有2人参与

最近在计算4原子铜簇A与有机大分子B(约100个原子),大分子上含苯环、碳氢氧元素,再求铜簇与苯环上吸附时,发现吸附能过大约3ev,后来在高手指导下我把簇模型和大分子放开有几十个埃米左右再优化,发现能量值并不等于之前铜簇和大分子的能量和,前后相差1.6ev.从理论上讲应该是相等的啊,因为没有相互作用吗。后来我把他们簇模型和大分子整体固定住,优化后发现能量基本没有变化,和不固定时能量基本相等。想问问哪位大神遇到过此种问题,怎么解决的,谢谢。
另外,我算过两个同簇模型放在一块但是相距很远,优化后能量是原来单个簇的能量的两倍,且把同样的两个大分子放在一起但是相聚很远优化后,能量是单个分子的两倍。
我的参数如下:Spin_polarization             unrestricted
Charge                        0
Basis                         dnp
Pseudopotential               dspp
Functional                    gga(p91)
Aux_density                   octupole
Integration_grid              medium
Occupation                    thermal 0.0050
Cutoff_Global                 4.0000 angstrom
Scf_density_convergence       1.0000e-005
Scf_charge_mixing             0.2000
Scf_spin_mixing               0.5000
Scf_iterations                999
Scf_diis                      6 pulay
回复此楼

» 猜你喜欢

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

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

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

shangshui123

新虫 (职业作家)

引用回帖:
2楼: Originally posted by mywai520 at 2014-04-27 22:20:35
放开几十A那还是一样的结构,都不在一个体系,怎么能一样。可能是我理解错了,能不能讲清楚些!!!!!

放开几十个A是想让他们没有作用,在一个体系中的,DMOL是无限延伸的,放进的话,比如十个A,结果也是一样的,我没有建板模型,因此是在一个体系中的
3楼2014-04-27 22:30:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 16 个回答

mywai520

铁杆木虫 (著名写手)


【答案】应助回帖

感谢参与,应助指数 +1
放开几十A那还是一样的结构,都不在一个体系,怎么能一样。可能是我理解错了,能不能讲清楚些!!!!!
2楼2014-04-27 22:20:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mywai520

铁杆木虫 (著名写手)


引用回帖:
3楼: Originally posted by shangshui123 at 2014-04-27 22:30:39
放开几十个A是想让他们没有作用,在一个体系中的,DMOL是无限延伸的,放进的话,比如十个A,结果也是一样的,我没有建板模型,因此是在一个体系中的...

你放在一起的时候,他们都成键了。成键要放去能量,或者比较近的情况下有范德瓦尔斯力作用,这些都会影响体系。

你放那么远,就是两个没有干涉的体系。这怎么能说是同一体系?????
4楼2014-04-27 22:35:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mywai520

铁杆木虫 (著名写手)


不知道你的所谓理论上应该如此,这是什么理论?若是有根据,我想看看。
5楼2014-04-27 22:36:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 让我中一个面上吧! +12 大萍1987 2026-08-20 15/750 2026-08-23 20:44 by 大萍1987
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 19:22 by lXgU0RvNTDC9
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 19:22 by lXgU0RvNTDC9
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 2JOx3r2CYEgw 2026-08-21 9/450 2026-08-23 18:46 by lXgU0RvNTDC9
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 4/200 2026-08-23 15:30 by KM5EcsNQRBPn
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 6/300 2026-08-23 14:21 by KM5EcsNQRBPn
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 5/250 2026-08-23 04:07 by OEbVnUOu01ol
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 04:04 by OEbVnUOu01ol
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 03:52 by OEbVnUOu01ol
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:43 by OEbVnUOu01ol
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] 范进中举一文的中心思想 +3 炎黄贵胄 2026-08-22 4/200 2026-08-22 20:02 by mycsru
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见