计算化学公社

标题: 求助计算化学位移张量的输出涉及超晶胞和晶体、价电子和总电子，该如何使用 [打印本页]

作者
Author: quilla 时间: 2023-11-22 22:31
标题: 求助计算化学位移张量的输出涉及超晶胞和晶体、价电子和总电子，该如何使用
参照官网的例子LCHIMAG - Vaspwiki，计算化学位移张量，打开LCHIMAG，OUTCAR关于这部分的输出如下图：
第一列和第二列分别为除去 G=0 贡献的列对于分子的超晶胞计算是有用的。包括 G=0 贡献的列适用于晶体。
(absolute, valence only):这一部分表示只考虑了价电子的化学位移。对于每个原子，给出了在排除G=0贡献时计算的化学位移（ISO_SHIFT）、SPAN和SKEW的数值。
(absolute, valence and core):这一部分表示考虑了价电子和内层电子（核心电子）的化学位移。对于每个原子，给出了在包含G=0贡献时计算的化学位移（ISO_SHIFT）、SPAN和SKEW的数值。

请教大家应该如何选用这些信息来计算化学位移？
如果应用在电解液体系，1M LiFSi/DME（所建立的盒子是周期性的），应该选取那一模块的信息来计算化学位移呢？

作者
Author: zjxitcc 时间: 2023-11-23 17:03
用右下方-6547.9848那一列数据，这是每个原子的各向同性屏蔽值的相反数。要得到核磁化学位移，还要计算该元素在核磁标准物中的各向同性屏蔽值。核磁化学位移=标准物中目标元素的各向同性屏蔽值 - 目标物种目标元素的各向同性屏蔽值。每种元素的核磁标准物不一样，也可以不只有一种，这个得看实验上用的什么对比。例如溶液核磁H^1-NMR的标准物是四甲基硅烷（TMS）。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)