用户工具

站点工具

本页面的其他翻译:
  • zh

adf:inductiondipole

如何计算诱导偶极矩

本例以HCl-H2O二聚体为例,演示诱导偶极的计算。三个及其以上的分子(片段),计算过程完全是相似的。计算过程可以考虑使用Excel表格的计算功能。

计算

首先应对二聚体进行结构优化,不过这种弱相互作用体系的优化,泛函的选择应该使用色散修正泛函,基组也应该使用高极化基组,比如TZ2P。本例中,省略该步骤,假定结构已经优化好了。

做三个计算:

  1. 二聚体的单点计算(需要设置details > Symmetry > Symbol:NoSym);
  2. 二聚体删掉HCl后的单点计算(需要设置details > Symmetry > Symbol:NoSym);
  3. 二聚体删掉H2O后的单点计算(需要设置details > Symmetry > Symbol:NoSym)

因为这里不关心键能,所以不使用色散修正泛函,而且使用色散修正泛函还可能导致各种范德华弱作用被重复考虑,因此这里使用普通泛函即可:

另外设置details > Symmetry > Symbol:NoSym,保存名为H2O-HCl.adf。

二聚体删掉HCl后的单点计算:(为了保证原子的坐标不变,必须)在第一个计算的ADFinput窗口,选中HCl分子,del键删除,之后另存任务,例如H2O.adf,并运行即可;

二聚体删掉H2O后的单点计算:(为了保证原子的坐标不变,必须)在第一个计算的ADFinput窗口,选中H2O分子,del键删除,之后另存任务,例如HCl.adf,并运行即可;

分析

打开H2O-HCl.out,按照分子的电偶极矩的计算中所示的方法,分别计算出二聚体中HCl和H2O的偶极矩。该计算得到的两个分子的总偶极矩,分别与H2O.out和HCl.out中所示的H2O和HCl的固有偶极矩相减(前者减去后者) 就分别得到两个分子诱导偶极矩

adf/inductiondipole.txt · 最后更改: 2020/12/01 14:03 由 liu.jun

© 2014-2022 费米科技(京ICP备14023855号