Let
Ω
\Omega
be a strictly pseudoconvex domain in
C
n
\mathbb {C}^n
with
C
k
+
2
C^{k+2}
boundary,
k
≥
1
k \geq 1
. We construct a
∂
¯
\overline \partial
solution operator (depending on
k
k
) that gains
1
2
\frac 12
derivative in the Sobolev space
H
s
,
p
(
Ω
)
H^{s,p} (\Omega )
for any
1
>
p
>
∞
1>p>\infty
and
s
>
1
p
−
k
s>\frac {1}{p} -k
. If the domain is
C
∞
C^{\infty }
, then there exists a
∂
¯
\overline \partial
solution operator that gains
1
2
\frac 12
derivative in
H
s
,
p
(
Ω
)
H^{s,p}(\Omega )
for all
s
∈
R
s \in \mathbb {R}
. We obtain our solution operators via the method of homotopy formula. A novel technique is the construction of “anti-derivative operators” for distributions defined on bounded Lipschitz domains.