The Helmholtz decomposition theorem states that any sufficiently smooth rapidly decaying vector field can be uniquely resolved into the sum of a divergence-free rotational solenoidal component F_curl obeying ∇ · F_curl = 0 and a curl-free irrotational gradient component F_div