Let Ω be a strictly pseudoconvex domain in C n with C k+2 boundary, k ≥ 1. We construct a ∂ solution operator (depending on k) that gains 1 2 derivative in the Sobolev space H s,p (Ω) for any 1 < p < ∞ and s > 1 p − k. If the domain is C ∞ , then there exists a ∂ solution operator that gains 1 2 derivative in H s,p (Ω) for all s ∈ R. We obtain our solution operators through the method of homotopy formula; a new feature is the construction of "anti-derivative operators" on distributions defined on bounded Lipschitz domains.