交换关联能量、势与应力的计算

我们描述了一种在广义梯度近似下计算总能量、有效势和应力张量中交换关联 (XC) 贡献的方法。我们避免使用 $E_{xc}[\rho]$ 泛函导数的解析表达式，因为这些表达式依赖于电子密度 $\rho$ 的不连续二阶导数。相反，我们首先在实空间网格上用积分近似 $E_{xc}$，然后计算其在网格点上对密度的偏导数。这确保了计算出的总能量、势和应力之间的精确一致性，并避免了对二阶导数的需求。我们展示了该方法的一些应用，该方法仅需网格（可能非均匀）上的（自旋）电子密度值，并返回常规的（局域）XC 势。