I would like to show the identity of the 2 Stage Least Squares estimator and the control function estimator.
Assume a linear regression model
$$y = X\beta + u$$ where $X =[X_1 \ X_2]$ is $n \times k$ and where $X_2$ is endogenous ($\mathbb E[x_{2i}u_i]\not =0$). Let $Z$ be a $n \times l$ matrix $l\geq k$ of instruments (that includes $X_1$ which is assumed to exogenous $\mathbb E[x_{1i}u_i]=0$) and
$X = Z\Gamma + V$
where $\Gamma:=\mathbb E[z_iz_i^\top]^{-1}\mathbb E[z_i y_i]$ such that by construction $\mathbb E[z_iv_i^\top]=0$
The 2SLS Estimator
The two stage least squares estimator is defined as $\hat \beta_{2SLS}:=(\hat X^\top X)^{-1}(\hat X^\top y)$ where $\hat X = Z(Z^\top Z)^{-1} Z X = P_Z X$.
The control function approach regress $X$ on $Z$ to get residuals $\hat V$
$X = Z\hat \Gamma + \hat V$ and then includes these residuals in a regresssion of $y$ on $X$ and $\hat V$ to get the estimator
$$ \begin{bmatrix}\hat b \\ \hat \rho \end{bmatrix} = \begin{bmatrix} X^\top X & X^\top \hat V \\ \hat V^\top X & \hat V^\top \hat V\end{bmatrix}^{-1} \begin{bmatrix} X^\top y \\ \hat V^\top y\end{bmatrix} $$
where the result I am looking for then is that $\hat b = \hat \beta_{2SLS}$.