adf:如何为对称性的分子指定电子的占据方式
差别
这里会显示出您选择的修订版和当前版本之间的差别。
| 两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版 | ||
| adf:如何为对称性的分子指定电子的占据方式 [2020/11/12 13:31] – [如何指定电子的占据方式] liu.jun | adf:如何为对称性的分子指定电子的占据方式 [2023/10/26 18:05] (当前版本) – 移除 liu.jun | ||
|---|---|---|---|
| 行 1: | 行 1: | ||
| - | ======如何指定电子的占据方式====== | ||
| - | 在ADFinput > details > Run Script中，或run文件中中添加IrrepOccupations关键字。格式示范如下： | ||
| - | < | ||
| - | ...... | ||
| - | Engine ADF | ||
| - | Basis | ||
| - | Type DZP | ||
| - | PerRegion Region=Region_1 Type=DZP Core=None | ||
| - | PerRegion Region=Region_2 Type=DZP Core=Large | ||
| - | PerRegion Region=Region_3 Type=TZP Core=None | ||
| - | End | ||
| - | ALLOW FRACTIONAL | ||
| - | XC | ||
| - | GGA PBE | ||
| - | End | ||
| - | IrrepOccupations | ||
| - | DELTA.g 0 0 | ||
| - | DELTA.u 0 | ||
| - | PI.g 4 0 | ||
| - | PI.u 4 2 0 | ||
| - | SIGMA.g 2 0 | ||
| - | SIGMA.u 4 0 | ||
| - | End | ||
| - | EndEngine | ||
| - | ...... | ||
| - | </ | ||
| - | |||
| - | ====Restricted方法==== | ||
| - | 再次强调：如果体系是高自旋态，例如三重态，那么有2个电子自旋朝同一个方向。这个时候，如果用Restricted方法，往往都很难收敛到正确的占据方式上去，经常出现LUMO比HOMO还低的情况，这要尤其注意，这也是导致[[adf: | ||
| - | <code bash> | ||
| - | OCCUPATIONS | ||
| - | A1 4.0 | ||
| - | B1 2.0 | ||
| - | B2 2.0 | ||
| - | END | ||
| - | </ | ||
| - | 其中A1、B1、B2是分子轨道的不可约表示名，后面的数字表示该不可约表示的电子个数。 | ||
| - | |||
| - | <code bash> | ||
| - | OCCUPATIONS | ||
| - | A 14.0 1.0 1.0 1.0 | ||
| - | END | ||
| - | </ | ||
| - | <color blue> | ||
| - | |||
| - | 对Restricted方法的计算，用户也可以分别指定alpha电子、beta电子的个数，例如上面的占据方式和下面是等价的： | ||
| - | <code bash> | ||
| - | OCCUPATIONS | ||
| - | A1 2.0//2.0 | ||
| - | B1 1.0//1.0 | ||
| - | B2 1.0//1.0 | ||
| - | END | ||
| - | </ | ||
| - | |||
| - | 如果是开壳层体系（有的电子的自旋没有配对），Restricted方法也允许用这种alpha、beta电子分开指定的方式，来实现高自旋态，这种情况，就只能手写类似如下的内容，到run文件或者User Input里面： | ||
| - | <code bash> | ||
| - | OCCUPATIONS | ||
| - | A 7.0//10.0 | ||
| - | END | ||
| - | </ | ||
| - | 表示有7个alpha电子、10个beta电子。但遗憾的是，这样指定，基本上绝大部分情况，都不能收敛到正确的占据方式上去。例如，最后可能电子的占据方式是： | ||
| - | |||
| - | <code bash> | ||
| - | OCCUPATIONS | ||
| - | A 7.0//7.0 0.0 1.0 1.0 1.0 | ||
| - | END | ||
| - | </ | ||
| - | 也就是：最低空轨道（本例中是第8个beta轨道）很可能比最高占据轨道（本例中为第11个beta轨道），能级还要低。在进行开壳层片段分析的时候，尤其要注意这一点。如果真实的占据方式变成这样，而这种占据就是你想要的，那么就需要按照这种方式设置Fragoccupations（关于该关键字的作用，参考[[adf: | ||
| - | ====Unrestricted方法==== | ||
| - | 用户只能分别指定alpha、beta电子的个数： | ||
| - | |||
| - | <code bash> | ||
| - | OCCUPATIONS | ||
| - | A1 2.0//2.0 | ||
| - | B1 1.0//1.0 | ||
| - | B2 1.0//1.0 | ||
| - | END | ||
| - | </ | ||
| - | 表示A1不可约表示，Spin up占据2个电子，Spin down2个电子；B1不可约表示，Spin up占据1个电子，Spin down1个电子；B2不可约表示，Spin up占据1个电子，Spin down1个电子。 | ||
| - | 或者： | ||
| - | <code bash> | ||
| - | OCCUPATIONS | ||
| - | A1 11.0//10.0 | ||
| - | B1 1.0//1.0 | ||
| - | B2 1.0//1.0 | ||
| - | END | ||
| - | </ | ||
| - | 表示A1不可约表示，Spin up的电子前11个轨道占据11个电子，Spin down前10个轨道占据10个电子；B1不可约表示，Spin up占据1个电子，Spin down1个电子；B2不可约表示，Spin up占据1个电子，Spin down1个电子。 | ||
| - | |||
| - | <code bash> | ||
| - | OCCUPATIONS | ||
| - | A1 10.0 0.33333333333333 0.33333333333333 0.33333333333333// | ||
| - | B1 1.0//1.0 | ||
| - | B2 1.0//1.0 | ||
| - | END | ||
| - | </ | ||
| - | 表示A1不可约表示，Spin up的电子前10个轨道占据10个电子，第11、12、13轨道分别占据1/ | ||
| - | |||
| - | 对于没有对称性的分子，只有一个不可约表示A： | ||
| - | <code bash> | ||
| - | OCCUPATIONS | ||
| - | A 2.0 0.0 1.0// 4.0 | ||
| - | END | ||
| - | </ | ||
| - | 表示A不可约表示Spin up的占据方式是：1 1 0 1，Spin Down的占据方式是1 1 1 1。 | ||
adf/如何为对称性的分子指定电子的占据方式.1605159078.txt.gz · 最后更改: 2020/11/12 13:31 由 liu.jun