这里会显示出您选择的修订版和当前版本之间的差别。
两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版 | ||
adf:sfotombo [2022/12/07 19:27] – liu.jun | adf:sfotombo [2022/12/12 22:15] (当前版本) – liu.jun | ||
---|---|---|---|
行 6: | 行 6: | ||
其中 | 其中 | ||
* A、B为2个原子 | * A、B为2个原子 | ||
- | * D、S分别为密度矩阵、原子轨道重叠矩阵,在AMSinput → Details → Output Details,勾选Density matrix、SFO overlap matrix,同时AMSinput → Details → SCF → Use old SCF即可输出密度矩阵(基于SFO)、原子轨道(即SFO)重叠矩阵。当然还需要取消对称性(AMSinput → Details → Symmetry → NOSYM),否则重叠积分数值可能不正确。在输出文件中搜索“Pmat-Low”即可看到密度矩阵,搜索“SFO Overlap Matrix”即可看到SFO重叠矩阵 | + | * D、S分别为密度矩阵、原子轨道重叠矩阵。 |
* μ、λ是SFO的编号序号,在AMS计算输出结果中,SCM → Output → Properties → SFO construction可以看到SFO的编号列表 | * μ、λ是SFO的编号序号,在AMS计算输出结果中,SCM → Output → Properties → SFO construction可以看到SFO的编号列表 | ||
- | 如此,即可分析具体某个轨道重叠,对MBO的贡献大小。 | + | ADF中,可以使用**EXTENDEDPOPAN**关键词,给出每个角动量对Mayer键级(以及Mulliken atom-atom布居)的贡献: |
+ | < | ||
+ | Engine ADF | ||
+ | Basis | ||
+ | Type DZP | ||
+ | Core None | ||
+ | End | ||
+ | XC | ||
+ | GGA PBE | ||
+ | End | ||
+ | Relativity | ||
+ | Level None | ||
+ | End | ||
+ | BondOrders | ||
+ | TypeForAMS Mayer | ||
+ | End | ||
+ | ExtendedPopan Yes | ||
+ | EndEngine | ||
+ | </ | ||
+ | |||
+ | 另外,on the resulting adf.rkf(dmpkf adf.rkf > ASCII.txt命令转化为文本文件,TAPE15需要类似转化为文本文件) one can find the density matrix on AO (SCF%Pmat_A,此处指上一行为SCF,下一行为Pmat_A,A%SFO类似), | ||
+ | |||
+ | 从这个角度,应该是可以计算SFO对Mayer键级的贡献的。 | ||
+ | |||
+ | 密度矩阵与重叠矩阵的存储顺序: | ||
+ | |||
+ | Smat and Pmat are symmetric stored in packed storage mode, meaning consecutively: | ||
+ | |||
+ | Suppose one has nsfo SFOs and nbas basis functions (AOs,BAS) A%SFO is stored as coef(nbas, | ||
+ | |||
+ | SFO_i = sum_j coef_ij bas_j, and first all coefficient for SFO_1, etcetera. Note that one has to take into account the so called npart indices, which are stored in A%npart in case of symmetry NOSYM, the npart indices are indices that give for each of the nbas functions, the number of the basis function in the list of all basis functions. |