24小时热门版块排行榜    

查看: 2489  |  回复: 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的回帖

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的回帖
查看全部 6 个回答

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

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 能否退出参与的面上项目解除限项 +12 koalala 2026-08-24 14/700 2026-08-24 12:41 by S09300097
[基金申请] 建议基金发布提前给出明确的时间点 +12 kulium 2026-08-21 15/750 2026-08-24 12:02 by iaeyuan
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 明天应该可查了!? +4 chengyan1220 2026-08-23 4/200 2026-08-24 08:47 by gloomy6159
[基金申请] 什么时候开奖? +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 看来今天不会放榜了? +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 紫杉醇
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 今天放榜没戏了吧 +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
[基金申请] 明天放榜? +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 家与远方
信息提示
请填处理意见