If we assume we are fitting the following regression model using ordinary least squares or maximum likelihood: $$ Y_i=X_i\hat{\beta} +\varepsilon_i \\ \varepsilon_i \sim N(0, \hat\sigma^2)$$
Computing $\hat\beta$, $\text{Var}(\hat\beta)$, and $\hat\sigma^2$ is straightforward using the usual formulas, but what are the formulas for $\text{Var}(\hat\sigma^2)$ and $\text{Cov}(\hat\beta, \hat\sigma^2)$? That is, what is the variance of the residual variance and what is its covariance with the coefficient estimates? The residual variance is an estimated quantity just like the coefficients, so I expect it to have a sampling variance and a sampling covariance with the model coefficients.