求助QE计算单水分子O-H键解离势能面：关于开壳层自旋与smearing设置的合理性探讨

小葱拌豆腐

各位老师、前辈好：

我目前正在训练一个用于液态水分子动力学模拟的 Deep Potential (DP) 机器学习势函数。为了修复原始力场在长程外推时的非物理现象（MD模拟中水分子异常解离），我需要人为补充一部分 单水分子 O-H 键拉伸（0.7 Å ~ 2.5 Å）的势能面扫描数据 作为训练集的物理边界约束。

因为我需要对标原文献，选用的泛函是 SCAN。由于计算软件必须统一使用 Quantum ESPRESSO (QE) 来提供能量和受力，我在拉伸至 1.5 Å 及更远距离时遇到了 SCF 不收敛的问题（HOMO-LUMO gap 缩小导致电子占据震荡）。

为了解决这个问题并描述解离产生的双自由基（OH• 和 H•），我采用了以下计算策略：

• 开启了自旋极化 nspin = 2，并设置了 starting_magnetization 打破对称性。
• 【核心疑问点】 放弃了分子常用的 occupations = 'fixed'，改用 occupations = 'smearing'，配合冷展宽 smearing = 'mv'，并给了一个极小的展宽值 degauss = 0.005 Ry（约 0.07 eV）。
  
  

以下是我在拉伸至 2.0 Å 时的完整 QE input 文件：

&CONTROL
  calculation = 'scf'
  restart_mode = 'from_scratch'
  prefix = 'h2o_scan_14'
  pseudo_dir = '../pseudo'
  outdir = './tmp'
  disk_io = 'none'
  tprnfor = .true.
/
&SYSTEM
  ibrav = 0
  nat = 3
  ntyp = 2
  
  ecutwfc = 110
  input_dft = 'SCAN'
  
  ! 自旋极化设置
  nspin = 2
  starting_magnetization(1) = 0.2
  starting_magnetization(2) = 0.1
  
  ! 占据态与展宽设置
  occupations = 'smearing'
  smearing = 'mv'      
  degauss = 0.005      
  
  nosym = .true.
/
&ELECTRONS
  conv_thr = 1.0d-8
  mixing_beta = 0.3
  electron_maxstep = 200
/
ATOMIC_SPECIES
O  15.999  O.upf
H  1.008   H.upf
​
CELL_PARAMETERS (angstrom)
12.0 0.0 0.0
0.0 12.0 0.0
0.0 0.0 12.0
​
ATOMIC_POSITIONS (angstrom)
O  0.00000000  0.00000000  0.00000000
H  2.00000000  0.00000000  0.00000000
H  -0.24036323  0.92942083  0.00000000
K_POINTS (gamma)

我的主要问题如下，恳请各位前辈指点：

• Smearing 用于孤立分子的合理性： 在拉伸断键这种接近简并状态（过渡态/双自由基）的情况下，使用微小的 smearing 来替代 fixed 以促进 SCF 收敛，这种做法在量化计算（或为拟合势函数提供参考数据时）是否在物理上是可接受的妥协？能量误差是否在可控范围内？
• 自旋态的指定： 我仅给了初始磁矩让其自发极化（Broken Symmetry），没有显式指定总自旋多重度（比如通过 tot_magnetization）。对于 H2O 解离成基态的 OH 和 H 自由基，这样的设置能否大概率收敛到正确的基态物理面？
• 有无更好的 QE 替代方案？ 如果这种做法在审稿时容易受到质疑，在 QE 的框架下，有没有更严谨的控制参数可以保证这种开壳层断键体系的 SCF 顺利收敛？