24小时热门版块排行榜    

查看: 2532  |  回复: 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一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 j5mmowWXNWxv 2026-09-05 6/300 2026-09-06 18:41 by eObJPa4gWmp6
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:21 by eObJPa4gWmp6
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:01 by eObJPa4gWmp6
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 17:49 by eObJPa4gWmp6
[文学芳草园] 两块石头 +7 阿美_Lml888 2026-09-03 9/450 2026-09-06 17:48 by myrtle
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 17:21 by eObJPa4gWmp6
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 17:09 by eObJPa4gWmp6
[博后之家] 广西大学-广州大学招聘博士后 欢迎广大优秀人才!!! +5 SCSIOyxguo 2026-09-03 8/400 2026-09-06 15:06 by 可爱小花朵
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 13:09 by eObJPa4gWmp6
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 13:09 by eObJPa4gWmp6
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 12:49 by eObJPa4gWmp6
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:29 by eObJPa4gWmp6
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 7/350 2026-09-06 11:20 by eObJPa4gWmp6
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 11:09 by eObJPa4gWmp6
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 11:00 by eObJPa4gWmp6
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 4/200 2026-09-06 07:20 by E1iiBLMU2XnD
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 CRJ0Aq2FtVH0 2026-09-03 15/750 2026-09-06 05:48 by E1iiBLMU2XnD
[文学芳草园] 初秋的晨风 +4 阿美_Lml888 2026-09-04 6/300 2026-09-05 13:51 by 阿美_Lml888
[基金申请] 学科评审组评审是指会评吗? +6 瞬息宇宙 2026-08-31 9/450 2026-09-05 13:11 by haizai99
[基金申请] 要骂人了,新模版改版就是要淡化问题凝练这种虚的东西,结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
信息提示
请填处理意见