The expressions that need to be evaluated in this step are quite general, as they depend on both the nature of the solute, and the specific differentiation variable λ. For this reason, in the most general case, the step should be carried out in the host code, e.g. in case of the Fock matrix.