Checking the MLE: From your specification of the problem, your log-likelihood function is:
$$\begin{equation} \begin{aligned}
\mathcal{l}_{\boldsymbol{x},\boldsymbol{y}}(\theta, \lambda)
&= \sum_{i=1}^m \ln p (x_i | \lambda) + \sum_{i=1}^n \ln p (y_i | \theta, \lambda) \\[8pt]
&= \sum_{i=1}^m (\ln \lambda - \lambda x_i) + \sum_{i=1}^n (\ln \theta + \ln \lambda - \theta \lambda y_i) \\[8pt]
&= m ( \ln \lambda - \lambda \bar{x} ) + n ( \ln \theta + \ln \lambda - \theta \lambda \bar{y}).
\end{aligned} \end{equation}$$
This gives the score functions:
$$\begin{equation} \begin{aligned}
\frac{\partial \mathcal{l}_{\boldsymbol{x},\boldsymbol{y}}}{\partial \theta}(\theta, \lambda)
&= n \Big( \frac{1}{\theta} - \lambda \bar{y} \Big), \\[8pt]
\frac{\partial \mathcal{l}_{\boldsymbol{x},\boldsymbol{y}}}{\partial \lambda}(\theta, \lambda)
&= m \Big( \frac{1}{\lambda} - \bar{x} \Big) + n \Big( \frac{1}{\lambda} - \theta \bar{y} \Big).
\end{aligned} \end{equation}$$
Setting both partial derivatives to zero and solving the resulting score equations yields the MLEs:
$$\hat{\theta}_{m,n} = \frac{\bar{x}}{\bar{y}} \quad \quad \quad \hat{\lambda}_{m,n} = \frac{1}{\bar{x}}.$$
(Note that in the case where $\bar{y} = 0$ the first of the score equations is strictly positive and so the MLE for $\theta$ does not exist.) This confirms your calculations of the MLE.
Adjusting the MLE to remove bias: Treating the MLE as a random variable we have:
$$\hat{\theta}_{m,n} = \frac{n}{m} \cdot \frac{\dot{X}}{\dot{Y}},$$
where $\dot{X} \equiv m \bar{X} \sim \text{Gamma} (m, \lambda)$ and $\dot{Y} \equiv n \bar{Y} \sim \text{Gamma} (n, \theta \lambda)$ are independent random variables. From this equation, the MLE is a scaled beta-prime random variable:
$$\hat{\theta}_{m,n} \sim \theta \cdot \frac{n}{m} \cdot \text{Beta-Prime}(m, n).$$
This estimator has expected value $\mathbb{E} (\hat{\theta}_{m,n}) = \frac{n}{n-1} \cdot \theta$, which means that it has positive bias. We can correct this bias by using the bias-adjusted MLE:
$$\tilde{\theta}_{m,n} = \frac{n-1}{n} \cdot \frac{\bar{X}}{\bar{Y}} \sim \theta \cdot \frac{n-1}{m} \cdot \text{Beta-Prime}(m, n).$$
Standard-Error of the adjusted MLE: The adjusted MLE is unbiased. It has variance:
$$\begin{equation} \begin{aligned}
\mathbb{V}(\tilde{\theta}_{m,n})
&= \int \limits_0^\infty \Big( \theta \cdot \frac{n-1}{m} \cdot r - \theta \Big)^2 \text{Beta-Prime} ( r | m, n) dr \\[8pt]
&= \theta^2 \cdot \frac{\Gamma(m) \Gamma(n)}{\Gamma(m+n)} \int \limits_0^\infty \Big( 1 - \frac{n-1}{m} \cdot r \Big)^2 r^{m-1} ( 1 + r )^{-m-n} dr \\[8pt]
&= \theta^2 \cdot \frac{n+m-1}{m(n-2)}.
\end{aligned} \end{equation}$$
The corresponding standard error is:
$$\text{se}(\tilde{\theta}_{m,n}) = \tilde{\theta}_{m,n} \cdot \sqrt{\frac{n+m-1}{m(n-2)}}.$$
Letting $\phi \equiv m/n$ and taking the limit as $n \rightarrow \infty$ we obtain the asymptotic approximation:
$$\text{se}(\tilde{\theta}_{m,n}) \approx \frac{\tilde{\theta}_{m,n}}{\sqrt{n-2}} \cdot \sqrt{\frac{1+\phi}{\phi}}.$$
This gives you both exact and approximate expressions for the standard error. I hope that is helpful. (Please make sure to review my algebra to make sure I haven't made a mistake!)
Best Answer
The first part is correct. Ok. The last one is wrong. The variance of the estimator is always unknown. One can estimate is using $\hat{\theta}$ but then it is random and varies for evey sample. Moreover, the variance of the estimator should converge to 0 as $n\to \infty$ (if not, the estimator would not be consitent and that is really not desirable)
There is just a mistype: You got: $$\sigma^2 = E[X^2] - E[X]^2 = \theta (1-\theta)$$ which is correct. But then $$SE_{\hat{\theta}} = \sqrt{Var(\hat{\theta})} = \sqrt{Var(2-\overline{X})} = \sqrt{\frac{\theta (1-\theta)}{n}}$$ where we use that $X_1,\dots, X_n$ are independent identically distributed.
So here I denoted by $SE_{\hat{\theta}}$ the standard error for the estimator $\hat{\theta}$ of $\theta$. So as you can see $SE_{\hat{\theta}}$ is unknown since $\theta$ is unknown, but one can estimate $\theta$ by $\hat{\theta} = 2-\overline{X}$ (stochastic). With our sample we have $\hat{\theta} = 1/3$ so an estimate of the standard error would be $$\hat{SE}_{\hat{\theta}} = \sqrt{\frac{\frac{1}{3} (1-\frac{1}{3})}{3}} = \sqrt{\frac{\frac{2}{9}}{3}} = \sqrt{2/27}.$$