Note that the total nuclear forces sum contributions due to the ES energy functional and the additional constraints, as listed in Eq. ( ), as well as the GS energy contributions E ζ GS . Regarding the discussed kernel options, additional terms for the gradient have to be considered for the transformed ADMM matrix of Eq. ( ), with all contributions implying the chain rule,