本教程针对 k-mesh 进行收敛分析，
gnuplot> p "/public/home/hanjinzhuo/local/src/abinit-10.4.7/tests/tutorial/Input/Work_bs/Work_tbs_5_444/tbs_5o_DS5_EXC_MDF" u 1:3 w l
gnuplot> rep "/public/home/hanjinzhuo/local/src/abinit-10.4.7/tests/tutorial/Input/Work_bs/Work_tbs_5_555/tbs_5o_DS5_EXC_MDF" u 1:3 w l
gnuplot> rep "/public/home/hanjinzhuo/local/src/abinit-10.4.7/tests/tutorial/Input/Work_bs/Work_tbs_5_666/tbs_5o_DS5_EXC_MDF" u 1:3 w l
可以看到计算时间
(abinit_env)[hanjinzhuo@mu012 Work_tbs_5_666$] grep "dataset:" log | grep "<<< TIME"dataset: 1 , wall: 1.70 [s] , cpu: 0.39 [s] <<< TIMEdataset: 2 , wall: 5.71 [s] , cpu: 1.39 [s] <<< TIMEdataset: 3 , wall: 4.09 [s] , cpu: 0.97 [s] <<< TIMEdataset: 4 , wall: 27.22 [s] , cpu: 6.76 [s] <<< TIMEdataset: 5 , wall: 02:16 [minutes] , cpu: 34.09 [s] <<< TIME
在 k-mesh = 6 6 6 的时候已经需要数分钟来计算，根据官方的数据
4x4x4: +Overall time at end (sec) : cpu= 112.4 wall= 112.4
5x5x5: +Overall time at end (sec) : cpu= 362.8 wall= 362.8
6x6x6: +Overall time at end (sec) : cpu= 914.8 wall= 914.8
8x8x8: +Overall time at end (sec) : cpu= 5813.3 wall= 5813.3
10x10x10: +Overall time at end (sec) : cpu= 20907.1 wall= 20907.1
12x12x12: +Overall time at end (sec) : cpu= 62738.2 wall= 62738.2
有限CPU上计算到 6x6x6 就已经是极限。官方依旧给图
正如预期的那样，谱随k点数收敛得很慢，我们最初做的4×4×4网格计算严重未收敛。用12×12×12 k网格得到的结果最为精确，然而即使是这种采样，也仅能使频率低于4.5 eV的结果收敛。这是所有Bethe - Salpeter（BS）计算的共同问题，即要实现全谱收敛极为困难。这一分析表明，我们仅能信任12×12×12结果在 [0:4.5] eV的能量范围内是可靠的；而要正确地描述更高能量处的谱，则需要包含更多的k点，且可能还需包含更多的带，以便正确考虑带的色散（即使在使用12×12×12时，RPA谱在高频率下也不收敛）。
应该强调的是，在这些收敛性测试中，zcut（虚部截断）起着非常重要的作用。例如，通过增加洛伦兹展宽，可以将8×8×8或10×10×10 k网格得到的结果与12×12×12的结果拉得更接近。在将理论与实验对比时，通常把zcut当作事后参数来处理，并将其选定为能与实验取得最佳一致性的值。