Sparse Domination for Bi-Parameter Operators Using Square Functions

Let $S$ be the dyadic bi-parameter square function $$Sf(x)^{2} = \sum_{R \in \mathcal{D}} |\langle f, h_{R} \rangle|^{2} \frac{1_{R}(x)}{|R|}.$$ We prove that if $T$ is a bi-parameter martingale transform and $f,g$ are suitable test functions, then there exists a sparse collection of rectangles $\mathcal{S}$ such that $$|\langle Tf, g \rangle| \lesssim \sum_{R \in \mathcal{S}} |R|(Sf)_{R}(Sg)_{R}.$$ We also extend this estimate to the case where $T$ is a bi-parameter cancellative dyadic shift and when $T$ is a paraproduct-free singular integral of Journ\'{e} type. Weighted estimates follow from the domination.