Derivative-order ambiguity in covariant Maxwell potential equations (${\nabla }_b{\nabla }^a{A}^b-{\nabla }^a{\nabla }_b{A}^b=-C^a{}_b{A}^b$)

Use the negative curvature convention defined in curvature-sign convention in Killing derivative identities, and $C_{ab}=C^c{}_{acb}$. Covariant derivatives of a vector obey $[{\nabla }_b,{\nabla }^a]A^b=-C^a{}_b{A}^b$. Hence the flat-space equivalent expressions $\square A^a-{\partial }_b{\partial }^a{A}^b$ and $\square A^a-{\partial }^a{\partial }_b{A}^b$ acquire a difference $+C^a{}_b{A}^b$ under the comma-to-semicolon prescription. Covariantizing $F_{ab}=2{\nabla }_{[a}{A}_{b]}$ and its Maxwell equations gives $\square A^a-{\nabla }^a({\nabla }_b{A}^b)+C^a{}_b{A}^b={\mu }_0{J}^a$. This form respects gauge invariance. Dropping the Ricci term is generally not gauge invariant. The two expressions coincide in a Ricci-flat vacuum spacetime; vanishing matter alone with nonzero cosmological constant does not imply Ricci flatness.