Finite-Gap 2D-Schrödinger Operators with Elliptic Coefficient

In general case the potential of the finite-gap Schrödinger operator $-\frac{\partial^2}{\partial x^2}+u(x)$ is expressed in terms of theta function of the spectral curve [5]. At the same time there are examples of finite-gap operators with elliptic potentials, for example, the Lamé operators $-\frac{\partial^2}{\partial x^2}+g(g+1)\wp(x)$ or the Treibich-Verdier operator $-\frac{\partial^2}{\partial x^2}+\sum_{i=0}^3a_i(a_i+1)\wp(x+\omega_i)$, where $\omega_i$ are semi-periods. Theorems 1 and 2 show that the same phenomena are possible in two-dimensional case.
$ $
Theorem 1. The Schrödinger operator
$$ \tag{1} H=\frac{\partial^2}{\partial z\partial\bar{z}} +a\bigg(\frac{\sqrt{g_0}-\wp'(az+b\bar{z})}{2\wp(az+b\bar{z})}\bigg)\frac{\partial}{\partial\bar{z}}-\frac{bg(g+1)\wp(az+b\bar{z})}{2a} $$
is finite-gap, where $\wp$ is elliptic Weierstrass function satisfying the equation
\begin{equation*}(\wp'(z))^2=\frac{2g(g+1)}{a^2}\wp(z)^3+g_2\wp(z)^2+g_1\wp(z)+g_0.\end{equation*}
The spectral curve of the operator $H$ is a hyperelliptic curve with genus $g$.
$ $
Thus for the operator $H$ theta functional formulas for the coefficients is reduced to the simpler formulas (1). Note that $H$ satisfies the identity
\begin{equation*}\bigg[H, -\frac{\partial^2}{\partial z^2}+g(g+1)\wp(az+b\bar{z})\bigg]=-2a\bigg(\frac{\partial}{\partial z}\bigg(\frac{\sqrt{g_0}-\wp'(az+b\bar{z})}{2\wp(az+b\bar{z})}\bigg)\bigg)H. \end{equation*}

Theorem 2. The Schrödinger operator
\begin{equation*}H=\frac{\partial^2}{\partial z\partial\bar{z}} +\frac{7a\wp'(az+b\bar{z})}{20g_2a^2-14\wp(az+b\bar{z})}\frac{\partial}{\partial\bar{z}}+ \frac{b\wp(az+b\bar{z})}{2a}\end{equation*}
is finite-gap, where $\wp$ is elliptic Weierstrass function satisfying the equation
\begin{equation*}(\wp'(z))^2=-\frac{1}{2a^2}\wp(z)^3+g_2\wp(z)^2-\bigg(\frac{7g_0}{10g_2a^2}+\frac{20g_2^2a^2}{49}\bigg)\wp(z)+g_0. \end{equation*}