The residual network outputs the vector potential A; the correction is ΔU = ∇×A. By the mathematical identity ∇·(∇×A) ≡ 0, the field is divergence-free by construction — the divergence visible in the right panel is at the floor of finite-difference truncation error, not a loss penalty. Pick a potential and see for yourself.