24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

Illusionist

银虫 (正式写手)

1ST强帖: 1
应助: 20 (小学生)
金币: 232.8
散金: 1001
红花: 15
帖子: 429
在线: 328.1小时
虫号: 1123329
注册: 2010-10-15
专业: 半导体晶体与薄膜材料

[求助] 如何用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
这样一来，两个分子的总能量并没有减少，就好象这两个分子并没有作用一样。
不知道怎么回事？

回复此楼