上个教程我们对 bs_loband 和 nband 两个参数进行了收敛分析，本教程我们针对 ecuteps 进行收敛分析，对现在的我们来说，单参数的收敛分析应该已经不具有难度了。
输入文件我们还是做如下对应修改
ndtset 4
ecuteps: 1 ecuteps+ 1
bs_loband 2 nband 7
bs_coulomb_term 11 # Full Coulomb term#irdwfk 1
#irdscr 1getwfk 99 # Trick to read the same file tbs_o3_DS20_WFK in each dataset
getscr 99 # Same trick for the SCR file
并且在外边创建长链接，
ln -s 444_shifted_WFK tbs_4o_DS99_WFK
ln -s 444_SCR tbs_4o_DS99_SCR
然后运行（这里需要注意，abinit 环境的 libnetcdf.so.15 会与 abipy 环境的冲突，所以运行的时候记得关闭 abipy_env ）
mpirun -n 8 abinit tbs_4.abi >& log_4
结束得到结果，然后画图
(abipy_env) (abinit_env)[hanjinzhuo@mu012 Work_bs$] gnuplotG N U P L O TVersion 6.0 patchlevel 3 last modified 2025-06-07Copyright (C) 1986-1993, 1998, 2004, 2007-2025Thomas Williams, Colin Kelley and many othersgnuplot home: http://www.gnuplot.infofaq, bugs, etc: type "help FAQ"immediate help: type "help" (plot window: hit 'h')Terminal type is now qt
gnuplot> p "tbs_4o_DS1_EXC_MDF" u 1:3 w l
gnuplot> rep "tbs_4o_DS2_EXC_MDF" u 1:3 w l
gnuplot> rep "tbs_4o_DS3_EXC_MDF" u 1:3 w l
gnuplot> rep "tbs_4o_DS4_EXC_MDF" u 1:3 w l
gnuplot>
看图发现， ecuteps=3和4基本无法区分，因此最终我们决定采用 ecuteps=3。
还要注意两个峰如何以不同的方式受到ecuteps变化的影响，第一个峰受到的影响最大。此行为是 与我们对硅的第一个峰具有很强的激子特性的肯定一致。