Sharp quantitative integral inequalities for harmonic extensions

We prove a quantitative version of a sharp integral inequality by Hang, Wang, and Yan for both the Poisson operator and its adjoint. Our result has the strongest possible norm and the optimal stability exponent. This stability exponent is not necessarily equal to 2, displaying the same phenomenon that Figalli and Zhang observed for the $p$-Sobolev inequality.