CDA分析求助

zhifanfan

各位老师好，我在做CDA分析时，使用的是Gaussian分析的波函数文件，报错提示：Warning: The coordinate of the fragments deviates from the complex distinctly! The result may be fully meaningless. Ple
ase check input files of your quantum chemistry code to make the coordinate of the fragments fully consistent with the complex. If you are a Gaussian user, please do not forget to add "nosymm" to avoid automatically translating and rotating the system to standard orientation Press ENTER button to continue，分析过程参考716博文，166博文，体系按照顺序建模，组合体的输入文件拆分成两个个体作为输入文件，不知道问什么会报错，同时分析结果也存在问题，转移电子：
      ========== Extended Charge decomposition analysis (ECDA) ==========
   Contribution to all occupied complex orbital:
Occupied, virtual orbitals of fragment  1:    2104.9361%        57.5451%
Occupied, virtual orbitals of fragment  2:    1409.2038%        57.6828%
   Contribution to all virtual complex orbital:
Occupied, virtual orbitals of fragment  1:    3664.2146%     11683.2908%
Occupied, virtual orbitals of fragment  2:    2038.3870%      3645.4902%
PL( 1) + CT( 1-> 2) =   73.2843      PL( 1) + CT( 2-> 1) =    1.1509
PL( 2) + CT( 1-> 2) =    1.1537      PL( 2) + CT( 2-> 1) =   40.7677
The net electrons obtained by frag. 2 = CT( 1-> 2) - CT( 2-> 1) =   72.1334

求助各位老师，该怎么修正结构或者优化过程。

zhifanfan

Uus/pMeC6H4-/キ 发表于 2025-1-9 14:31
Ce2O3-SiO2的孤立体系建模怎么都不像是符合化学直觉的样子，结果怎么可能合理。另外我很怀疑你是否完全清楚 ...

收到，谢谢老师

sobereva

使用Multiwfn做电荷分解分析(CDA)、绘制轨道相互作用图
http://sobereva.com/166

博文说得清楚得没法更清楚


写那堆IOp、pop=full完全多余

Uus/pMeC6H4-/キ

Ce2O3-SiO2的孤立体系建模怎么都不像是符合化学直觉的样子，结果怎么可能合理。另外我很怀疑你是否完全清楚Gaussian输入文件里写的关键词和iop的意思，特别是iop(5/13=1)。


Ce2O3和SiO2都是晶体的话，应当按周期性体系建模并采用CP2K等计算。（编辑：既然已经知道社长博文716了就去掉链接了。看前两帖似乎主要困难在建模上，在论坛搜“异质结 建模”可以找到讨论。）