This is the most general way to express the nuclear gradient when using a real orthogonal basis that explicitly dependents upon R. Note that although we introduced (defined) the short-hand notation [∇ α V ] in Eq. 26, it can be justified as the matrix of [∇ α V ] using the chain rule and the completeness relation of basis as follows