It is well known that if $x \sim \mathcal{N}(0, \sigma^2)$ and we have a sample size of $n$ observations of $x$, the distribution of $\frac{x}{\hat{\sigma}}$, where $\hat{\sigma}$ is the sample variance, will follow the Studnent-t law with $n-1$ degrees of freedom.
But what if the ML estimator of $\sigma$ is functionally penalized, i.e. $\sigma = f(\theta)$, and the estimator of $\sigma$, $\hat{\sigma}(\theta)$, can be only obtained numerically by means of MLEs of $\theta$? In that case, is there anything we can say about the distribution of $\frac{x}{\hat{\sigma}(\theta)}$?