24小时热门版块排行榜    

查看: 2496  |  回复: 5

Illusionist

银虫 (正式写手)

[求助] 如何用dmol3计算两个分子相互作用。

想要计算两个含有苯环的平面结构的有机分子相互作用
结果发现两个分子靠近后总能量并没有减少。
不知道是哪里设置有问题。
难道和对称性有关么?还是什么
急需求解!
不胜感激。

距离3.3A左右
使用的是MS5.5  DMOL3模块。
首先进行了几何优化:
# Task parameters
Calculate                     optimize
Opt_energy_convergence        1.0000e-005
Opt_gradient_convergence      2.0000e-003 A
Opt_displacement_convergence  5.0000e-003 A
Opt_iterations                50
Opt_max_displacement          0.3000 A
Initial_hessian               improved
Symmetry                      on
Max_memory                    2048

# Electronic parameters
Spin_polarization             restricted
Charge                        0
Basis                         dnp
Pseudopotential               none
Functional                    gga(p91)
Aux_density                   octupole
Integration_grid              fine
Occupation                    fermi
Cutoff_Global                 3.7000 angstrom
Scf_density_convergence       1.0000e-006
Scf_charge_mixing             0.2000
Scf_iterations                50
Scf_diis                      6 pulay

# Print options
Print                         eigval_last_it

# Calculated properties

然后计算能量

# Task parameters
Calculate                     energy
Symmetry                      on
Max_memory                    2048

# Electronic parameters
Spin_polarization             restricted
Charge                        0
Basis                         dnp
Pseudopotential               none
Functional                    pwc
Harris                        off
Aux_density                   octupole
Integration_grid              fine
Occupation                    fermi
Cutoff_Global                 3.7000 angstrom
Scf_density_convergence       1.0000e-006
Scf_charge_mixing             0.2000
Scf_iterations                50
Scf_diis                      6 pulay

# Print options
Print                         eigval_last_it

# Calculated properties

结果显示,
单个分子的总能量:
Energy components:
   Sum of atomic energies =    -1195.7413022Ha
                  Kinetic =      -12.9964458Ha
            Electrostatic =       -0.8916723Ha
     Exchange-correlation =        3.0765993Ha
        Spin polarization =        2.0531628Ha
两个分子的总能量:
Energy components:
   Sum of atomic energies =    -2391.4826044Ha
                  Kinetic =      -25.1271925Ha
            Electrostatic =       -2.6622324Ha
     Exchange-correlation =        6.1452670Ha
        Spin polarization =        4.1063257Ha
这样一来,两个分子的总能量并没有减少,就好象这两个分子并没有作用一样。
不知道怎么回事?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

redskywei

木虫 (著名写手)

【答案】应助回帖


uuv2010(金币+1): 多谢提示 2011-09-29 10:29:55
Illusionist(金币+1): 2011-09-29 14:46:47
3.3A的距离是范德华力作用范围吧,这种力第一性原理的计算效果本来就不好。dmol3估计更加不行了。要不你试试castep吧。
坠入无边的网~~
2楼2011-09-29 09:56:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Illusionist

银虫 (正式写手)

引用回帖:
2楼: Originally posted by redskywei at 2011-09-29 09:56:00:
3.3A的距离是范德华力作用范围吧,这种力第一性原理的计算效果本来就不好。dmol3估计更加不行了。要不你试试castep吧。

单算这两个分子的作用·
没有周期结构
castep好像算不了啊······
我看文献这个E(总的)-2*E(单个)=15Kj/mol左右
3楼2011-09-29 10:02:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

redskywei

木虫 (著名写手)

【答案】应助回帖


uuv2010(金币+1): 多谢提示 2011-09-29 10:30:04
Illusionist(金币+1): 2011-09-29 14:46:53
引用回帖:
3楼: Originally posted by Illusionist at 2011-09-29 10:02:58:
单算这两个分子的作用·
没有周期结构
castep好像算不了啊······
我看文献这个E(总的)-2*E(单个)=15Kj/mol左右

castep 里你建立超胞就可以了啊。先把单分子放在一个大超胞里算完后,再把两个放进去。
坠入无边的网~~
4楼2011-09-29 10:16:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Illusionist

银虫 (正式写手)

引用回帖:
4楼: Originally posted by redskywei at 2011-09-29 10:16:50:
castep 里你建立超胞就可以了啊。先把单分子放在一个大超胞里算完后,再把两个放进去。

有道理啊····
我试一试
一般而言,计算分子间的作用力也是这样算的么?
我说的是,有几个苯环这么大的有机物

DMOL3不使用周期结构计算和使用周期结构
会不会有什么不一样啊,好像
只要超胞足够大就结果就会一样,但是实际上是不是呢·····
5楼2011-09-29 10:32:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

redskywei

木虫 (著名写手)

【答案】应助回帖

★ ★
youzhizhe(金币+2): 谢谢交流。 2011-09-29 13:22:57
Illusionist(金币+1): 2011-09-29 14:46:58
Illusionist(金币+3): 2011-09-30 09:07:35
你问的我不知道怎么答了。建议你多看看相关文献吧。
DMOL3的使用周期与否的区别,自己算算不就知道了。我没算过。其实就算castep对分子间范德华力的描述也是不够的。最好描述都是自己写的程序,加入特别考虑足够。还是建议你看看文献。
坠入无边的网~~
6楼2011-09-29 10:57:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Illusionist 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 0XLacIJUOj8D 2026-08-24 7/350 2026-08-25 03:24 by BZKMTicpDhFj
[基金申请] 明天应该可查了!? +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 没有任何消息-是不是就凉了 +7 图啦图啦 2026-08-24 8/400 2026-08-24 22:00 by maomao_da
[基金申请] 人气不行了 +10 fansofjerry 2026-08-21 10/500 2026-08-24 21:03 by zhanghaozhu
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 19/950 2026-08-24 16:52 by iaeyuan
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +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
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见