我尝试复现文献中用deMon2k模拟环氧树脂吸附铝（铜）原子的XPS谱图的两篇论文（https://doi.org/10.1021/acsami.0c19616；https://doi.org/10.1021/acs.jpcc.9b07772）。这两篇论文模拟的效果挺好的（如图1），与实验数据比较吻合。其中模拟环氧树脂吸附铝原子的论文同时用了ΔKS和unrestricted Generalized Transition State theory (uGTS)两种方法；模拟环氧树脂吸附铜原子的，只用了uGTS的方法。我碰到了设置alpha轨道电荷的问题：在第一篇论文中，作者的原话是：2. Similarly, SP calculations for the ionized system with no electron in the α molecular orbital (MO) of the core−shell of C*。这句话是什么意思？是指设置碳原子核壳层的α分子轨道设置为0吗？该如何设置呢？完整的计算如图2（2.2.3）。
在第二篇文章中，作者的原话是：and (iii) finally, we perform a single-point calculation of the energy on the ionized molecule, with 0.33 electron on the 1s orbital for the atom of interest. We obtain the energy −εk(1/3) of the 1s orbital。如何将感兴趣的C原子和氧原子的1s轨道的电荷实现文中所描述的设置？完整的计算方法如图3（2.2.2）。 (, 下载次数 Times of downloads: 19)