| 查看: 2498 | 回复: 14 | |||
Illusionist银虫 (正式写手)
|
[求助]
计算两个分子相互作用
|
|
想要计算两个含有苯环的平面结构的有机分子相互作用 距离3.3A左右 结果发现两个分子靠近后总能量并没有减少。 用MS的DMOL3算的，结果算得0 再用g09算，也是不理想。文献是15Kj/mol左右，我只算了3Kj/mol。 应该怎么算才对。 还有这种计算是g09精度高还是MS的DMOL精度高。 |
回复此楼» 收录本帖的淘帖专辑推荐
 高斯 |
» 猜你喜欢
研究生招生
已经有0人回复
-
一志愿南大理论与计算化学方向（070300）342分，求调剂
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有70人回复
-
欢迎化学、化工、物理、材料、或计算机背景的考研咨询
已经有0人回复
-
山西大学资源与环境工程研究所化工分离/环境功能材料小组招收2026年博士生
已经有3人回复
-
一志愿南师大0703化学 275求调剂
已经有6人回复
-
A 区过线就可申！重庆双非公办高校 0856调剂，接受跨专业！纳米材料 / 润滑热门方向
已经有0人回复
-
欢迎调剂长江师范学院汤老师课题组的硕士研究生
已经有0人回复
-
调剂招生
已经有7人回复
-
河南理工大学化工学院招收调剂生啦！
已经有1人回复
-
中国电机工程学报的投稿系统
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 怎么能用量子化学计算含氮杂环芳香化合物与水分子的相互作用
已经有17人回复
- 分子间的两个苯环之间有Π...Π相互作用两个苯环之间的距离一般是多少
已经有5人回复
- 请问在Gaussian计算中怎样考虑弱相互作用？
已经有23人回复
- 如何用dmol3计算两个分子相互作用。
已经有5人回复
- 【求助】优化：怎样找到一个两个分子之间的相对位置
已经有26人回复
- 【求助】高斯计算分子结合能
已经有8人回复
- 【求助】请问有拟合分子力场参数的FORTRAN程序吗？
已经有13人回复
- 【求助】用ms计算两个分子的相互作用势
已经有20人回复
- 【讨论】如何计算两分子间的有效相互作用力
已经有12人回复
2楼2011-09-29 13:10:13
lihb734
铁杆木虫 (职业作家)
站在计算化学入门的门槛上
- QC强帖: 2
- 应助: 132 (高中生)
- 金币: 7495.4
- 散金: 9378
- 红花: 74
- 帖子: 3107
- 在线: 1177.3小时
- 虫号: 448403
- 注册: 2007-11-01
- 性别: GG
- 专业: 理论和计算化学
3楼2011-09-29 14:02:15
Illusionist
银虫 (正式写手)
- 应助: 20 (小学生)
- 金币: 232.8
- 散金: 1001
- 红花: 15
- 帖子: 431
- 在线: 329.1小时
- 虫号: 1123329
- 注册: 2010-10-15
- 专业: 半导体晶体与薄膜材料
|
我使用的是MS5.5 DMOL3模块算的，先算单分子，在算两个分子，然后E（两个分子）-2*E（单分子）。我并没有建立超胞使用周期结构，这样速度会快很多。 G09是别人帮我算的。我觉得对于有机物来说，可能G09会准确些，所以在这里求助。 首先进行了几何优化： # 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 这样一来，两个分子的总能量并没有减少，就好象这两个分子并没有作用一样。 |
4楼2011-09-29 14:45:11
lihb734
铁杆木虫 (职业作家)
站在计算化学入门的门槛上
- QC强帖: 2
- 应助: 132 (高中生)
- 金币: 7495.4
- 散金: 9378
- 红花: 74
- 帖子: 3107
- 在线: 1177.3小时
- 虫号: 448403
- 注册: 2007-11-01
- 性别: GG
- 专业: 理论和计算化学
5楼2011-09-29 14:52:08
Illusionist
银虫 (正式写手)
- 应助: 20 (小学生)
- 金币: 232.8
- 散金: 1001
- 红花: 15
- 帖子: 431
- 在线: 329.1小时
- 虫号: 1123329
- 注册: 2010-10-15
- 专业: 半导体晶体与薄膜材料
6楼2011-09-29 15:23:08
lihb734
铁杆木虫 (职业作家)
站在计算化学入门的门槛上
- QC强帖: 2
- 应助: 132 (高中生)
- 金币: 7495.4
- 散金: 9378
- 红花: 74
- 帖子: 3107
- 在线: 1177.3小时
- 虫号: 448403
- 注册: 2007-11-01
- 性别: GG
- 专业: 理论和计算化学
7楼2011-09-29 15:42:24
Illusionist
银虫 (正式写手)
- 应助: 20 (小学生)
- 金币: 232.8
- 散金: 1001
- 红花: 15
- 帖子: 431
- 在线: 329.1小时
- 虫号: 1123329
- 注册: 2010-10-15
- 专业: 半导体晶体与薄膜材料
8楼2011-09-29 15:53:09
Illusionist
银虫 (正式写手)
- 应助: 20 (小学生)
- 金币: 232.8
- 散金: 1001
- 红花: 15
- 帖子: 431
- 在线: 329.1小时
- 虫号: 1123329
- 注册: 2010-10-15
- 专业: 半导体晶体与薄膜材料
9楼2011-09-29 16:03:12
雨晨田木
金虫 (小有名气)
- 应助: 17 (小学生)
- 金币: 1177.7
- 散金: 91
- 红花: 1
- 帖子: 219
- 在线: 49.6小时
- 虫号: 1012487
- 注册: 2010-05-06
- 性别: GG
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★
小红豆(金币+3): 感谢交流 2011-09-29 20:37:19
Illusionist(金币+2): 谢谢 2011-09-30 09:07:06
小红豆(金币+3): 感谢交流 2011-09-29 20:37:19
Illusionist(金币+2): 谢谢 2011-09-30 09:07:06
| 可以试试cc-pvtz基组,能量要用ccsd算,还要做零点能校正和基组超越校正。具体你可以看这个文献：Interactions of Boranes and Carboranes with Aromatic Systems: CCSD(T) Complete Basis Set Calculations and DFT-SAPT Analysis of Energy Components （JPCA）10.1021/jp104411x |
10楼2011-09-29 18:43:24
