What XCint doesΒΆ
XCint integrates the exchange-correlation (XC) energy \(E_\text{xc}\) and the elements of the XC potential matrix \(V_\text{xc}\), as well as their derivatives with respect to electric field and/or geometric perturbations. The integration is performed on a standard numerical grid.
Before the XC energy can be computed, we require the densities. This is done in two steps:
Then the XC energy is computed using the XC energy density \(\epsilon_\text{xc}\) evaluated with the help of XCFun:
A similar strategy is used to compute \(V_\text{xc}\) matrix elements, again in two steps:
Note that XCFun is called only once and returns \(\epsilon_\text{xc}\) and \(v_\text{xc}\) in one go.
In the above scheme we work with batches of points in order to exploit vectorization and screening, making use of BLAS level 3 libraries to compute \(X_{kb}\) and \((V_\text{xc})_{kl}\).