In this arxiv paper (p. 11, eq. (3.2)) the authors claim that equation (3.2) is
... a quartic equation [...] which can be solved explicitly.
The equation in question is \begin{equation} %\tag{3.2} \frac{1}{m} = -z+\frac{\vartheta y}{\vartheta m-\varphi} - \frac{\left(\vartheta+\varphi\right)\left(1+\vartheta\varphi\right)y}{\left(\vartheta m-\varphi\right)\sqrt{\left[\left(1-\varphi\right)^{2}+m\left(1+\vartheta\right)^{2}\right]\left[\left(1+\varphi\right)^{2}+m\left(1-\vartheta\right)^{2}\right]}} \end{equation} with $\vartheta,\varphi\in\mathbb{R}$, $\left|\varphi\right|>0$, $y>0$ and $z,m\in\mathbb{C}^{+}$.
After some algebra, I arrive at \begin{align} \tag{2} \left(\left(\vartheta+\varphi\right)\left(1+\vartheta\varphi\right)y\right)m ={}& \overbrace{\left(\vartheta y m-\left(1+z m\right)\left(\vartheta m-\varphi\right)\right)}^{A} \\ &\underbrace{\cdot\sqrt{\left[\left(1-\varphi\right)^{2}+m\left(1+\vartheta\right)^{2}\right]\left[\left(1+\varphi\right)^{2}+m\left(1-\vartheta\right)^{2}\right]}}_{B} %&\underbrace{\cdot\sqrt{\left(1-\varphi\right)^{2}+m\left(1+\vartheta\right)^{2}} \sqrt{\left(1+\varphi\right)^{2}+m\left(1-\vartheta\right)^{2}}}_{B} \text{.} \end{align}
The term on the LHS is linear in $m$. Heuristically, the term on the RHS should behave as follows: The first factor $A$ is quadratic in $m$, the discriminant of the second factor is also quadratic under the root. By taking the root one gets orders $m$ and $m^{\frac{1}{2}}$. Multiplying this with $m^{2}$ and $m$ from factor $A$, one should get terms of orders $m^{3},m^{\frac{5}{2}},m^{2},m^{\frac{3}{2}},m,m^{\frac{1}{2}}$ and a constant. If one substitutes $x^{2}=m$, one has accordingly terms of order $x^{6},x^{5},x^{4},x^{3},x^{2},x$ and a constant.
Mathematica
solves (2) as a sextic equation and gives 6 roots (very long terms). Mathematica
also gives 6 roots if I assume
$\vartheta,\varphi\in\mathbb{R}$, $\left|\varphi\right|>0$, $y>0$ and $z,m\in\mathbb{C}^{+}$ as in the main paper.
Maybe the sextic equation collapses to quartic one by the assumption that $m\in\mathbb{C}^{+}$. If one substitutes $x^{2}=m$, the discriminant of $B$ is given by \begin{equation} %\tag{3} \left(\left(1-\varphi\right)^{2}+x^{2}\left(1+\vartheta\right)^{2}\right)\left(\left(1+\varphi\right)^{2}+x^{2}\left(1-\vartheta\right)^{2}\right) = \left(x\pm\frac{\left(\varphi+1\right)\mathsf{i}}{\vartheta-1}\right)\left(x\pm\frac{\left(\vartheta-1\right)\mathsf{i}}{\vartheta+1}\right) \text{,} \end{equation} where $\pm$ means that each factor occurs ones with positive and ones with negative sign. So in total we have on the RHS a product of four factors, each in $x$.
Because $\left|\varphi\right|<1$, the nominator of every root never changes in sign regardless of the value of $\varphi$. By the requirement that $m\in\mathbb{C}^{+}$, the possible values of $\vartheta$ get restricted. My argument goes like this: Since the roots of the discriminant are somehow related to the solution of the sextic equation, the sextic equation collapses to a quartic one?!
Is this reasoning right and can somebody please help to make it rigorous? And, most important, what are the four (assuming the authors in the linked paper are right) solutions to the first equation?