\section{Introduction}
\label{sec:introduction}

In quantum field theories there exist operators that are defined by singular boundary conditions imposed on dynamical fields in the path integral.
Equivalently, such operators can be defined through the coupling to a singular background field, by which the dynamical field is shifted.
It is also meaningful to define an operator by coupling the theory to a singular background field even when there is no dynamical field to be shifted.
Much progress has been made in recent years in understanding operators of these kinds, the so-called disorder operators.
Duality often exchanges disorder operators with more conventional operators defined as functions of dynamical fields.
In supersymmetric field theories, some correlation functions of disorder operators preserving part of supersymmetry can be exactly computed by localization.

In this note we study the $\mathcal{N}=(2,2)$ theory on Euclidean plane consisting of a single chiral multiplet without a superpotential.
A ``local'' disorder operator, the flavor vortex operator, can be defined by coupling the theory to a background gauge field for the flavor $U(1)$ symmetry of the form
\begin{equation} \label{vortex-singularity}
v_\mu dx^\mu \sim \alpha d\varphi\,,
\end{equation}
where $\varphi$ is the angular coordinate around the insertion point and $\alpha$ a real parameter.
By also turning on the background auxiliary field $D$ in the non-dynamical vector multiplet, half of the supercharges can be preserved.
The analysis of broken supercharges will reveal the connection to mirror symmetry \cite{Hori:2000kt}.

A closely related paper \cite{Hosomichi:2015pia} studying gauge vortex operators appeared recently.

\section{Flavor vortex operator}
\label{sec:flav-vort-oper}

A chiral multiplet of two-dimensional $\mathcal{N}=(2,2)$ supersymmetry consists of a complex scalar $\phi$, fermions $\psi_\mp$, and a complex auxiliary field $F$.
We put the theory on Euclidean plane with metric $ds^2=|dz|^2$, where $z=x^1+ix^2$.
Their transformations under the supercharges $(Q_\pm, \overline{Q}_\pm)$ are given in (\ref{SUSY-chiral}).
The free theory of a single chiral multiplet without a superpotential, possibly with a twisted mass, has a $U(1)$ flavor symmetry.
We let the fields in the chiral multiplet have charge $+1$.
Fields $(\sigma, \overline{\sigma}, v_\mu, \lambda_\pm, \overline{\lambda}_\pm, D)$ in a vector multiplet transform according to (\ref{SUSY-vector}).
We are interested in configurations of a non-dynamical vector multiplet for which the transformations under some supercharges vanish.
The twisted mass is simply the constant vev of the complex scalar $\sigma$.

It is cleaner to start with a smooth disorder operator of finite size, as in references \cite{Kapustin:2012iw,Hook:2013yda} that studied the three-dimensional case.
It is defined by coupling the chiral multiplet to a smooth background gauge field $v=v_\mu dx^\mu$.
The configuration preserves those supercharges for which the variation of the fermions in the flavor gauge multiplet vanishes.
From (\ref{SUSY-vector})  we see that  $Q_-$ and $\overline{Q}_+$ are preserved when the condition~$D =2i v_{z\overline{z}}$, or equivalently $D dx^1\wedge dx^2+dv=0$, is satisfied.
(Similarly $Q_+$ and $\overline{Q}_-$ are preserved when  $D = -2i v_{z\overline{z}}$.)
We characterize the field configuration by a real function $\rho=\rho(x)$, the ``vorticity density'', with support in a compact region:
\begin{equation}
D=-2\pi \rho\,,\qquad v_{z\overline{z}} = \pi i \rho \,.
\end{equation}
Other fields in the vector multiplet are set to zero.
This defines a disorder operator; let us call it the smeared vortex operator~$V[\rho]$.

For the action $S=\int d^2x \mathcal L,$ the Lagrangian
\begin{equation}\label{Lag}
  \begin{aligned}
  \mathcal L&=  D^\mu\overline\phi D_\mu\phi +\overline\sigma  \sigma  \overline\phi\phi - \overline{\phi} D\phi  -\overline F F
\\
&\qquad \qquad
+2i \overline\psi{}_+ D_{ z} \psi_+ -2 i \overline\psi{}_- D_{ {\overline{z}}} \psi_-  +\overline\sigma \overline\psi{}_+\psi_-  +\sigma \overline\psi_- \psi_+ \,,
\end{aligned}
\end{equation}
and the transformations (\ref{SUSY-chiral}), the supersymmetry current $S^\mu$ corresponding to the linear combination $\delta= i\epsilon_+ Q_- - i\epsilon_- Q_+  -i \overline\epsilon{}_+ \overline Q{}_- + i\overline \epsilon{}_- \overline Q{}_+$ of supercharges is given as
\begin{equation} \label{SUSY-current-chiral}
  \begin{aligned}
  S^z =   -4 \epsilon _- (D_{\overline{z}} \bar{\phi })\psi _++4 \bar{\epsilon }_-   (D_{\overline{z}}\phi )\bar{\psi }_++2 i \sigma  \epsilon _+ \bar{\phi} \psi _++2 i \bar{\sigma } \bar{\epsilon }_+ \phi \bar{\psi }_+  \,,
\\
  S^{\overline{z}} = + 4 \epsilon _+ (D_z \bar{\phi })\psi _--4 \bar{\epsilon }_+ (D_z \phi) \bar{\psi }_-+2 i \epsilon _-   \bar{\sigma } \bar{\phi }\psi _-+2 i \sigma  \bar{\epsilon }_- \phi \bar{\psi }_-  \,.
  \end{aligned}
\end{equation}
Its derivation is reviewed in Appendix \ref{ST}.
Let $n_\mu$ be the outward normal vector to the boundary of the support of $\rho$, and $ds$ the line element.
By the standard argument \cite{Pol-Ch2} the integral
\begin{equation} \label{contour-integral-normal-vector}
-  \oint ds\, n_\mu S^\mu = -  \int d^2x\, \partial_\mu S^\mu
\end{equation}
in the defining background for $V[\rho]$ equals the variation $\delta V[\rho]$.
A straightforward computation gives
\begin{equation}
\delta  V[\rho] =   4\pi  \int  d^2x \, \rho(x) ( \epsilon_- \overline{\phi}\,\psi_+ + \overline\epsilon_+ \overline{\psi}_- \phi)  \cdot V[\rho]\,.
\end{equation}
Terms proportional to $\epsilon_+$ and $\overline{\epsilon}_-$ have dropped out because $Q_-$ and $\overline{Q}_+$ are preserved.

We now take the limit
\begin{equation}
\rho \rightarrow  \alpha\cdot  \delta^2(x)\,,
\end{equation}
in which the smeared operator becomes a local operator inserted at the origin: $V[\rho]\rightarrow V_\alpha$.
It is the two-dimensional analog of the flavor vortex loop operator in three dimensions studied in \cite{MS,Kapustin:2012iw,Drukker:2012sr}.
(See also \cite{Yoshida:2014ssa,Assel:2015oxa,Beasley:2014jta} for the study of the same and the related operators in supersymmetric and non-supersymmetric settings.)

When $\alpha$ takes a non-integer real value, the parallel transport of matter fields around $V_\alpha$ produces a non-trivial phase; the fields and $V_\alpha$ are not mutually local.
In the terminology of \cite{Kapustin:2014gua,Aharony:2013hda,Gaiotto:2014kfa}, $V_\alpha$ is not a genuine local operator, but is an end of a topological line operator.

For the supersymmetry variations we get
\begin{equation}
 \overline{Q}_-\cdot V_\alpha = 4 \pi i \alpha\,    \overline{\psi }_- \phi\cdot  V_\alpha \,,\qquad
 Q_+\cdot  V_\alpha = 4 \pi i \alpha\, \overline{\phi}\, \psi_+ \cdot V_\alpha  \,.
\end{equation}
We see that the variations are proportional to $\alpha$ and $V_\alpha$.
This makes us propose that $V_\alpha$ can be written as
\begin{equation} \label{V-y}
  V_\alpha = e^{-\alpha y}\,,
\end{equation}
where the local operator $y$ is annihilated by $Q_-$ and $\overline{Q}_+$ and hence is a twisted chiral operator.
If we set
\begin{equation}\label{Y-Phi-relation-one}
\overline{\chi}_+ := i Q_+ y =  4 \pi \overline{\phi} \,\psi_+\,,
\quad
\chi_- :=- i \overline Q_- y = - 4 \pi \overline{\psi}_- \phi\,,
\quad
E := - 4 \pi \overline\psi_- \psi_+ \,,
\end{equation}
the fields $(y, \chi_-, \overline\chi_+, E)$ form a twisted chiral multiplet transforming as in (\ref{SUSY-twisted-chiral}).

The procedure reviewed in Appendix \ref{ST} gives the energy-momentum tensor
\begin{equation}\label{energy-momentum}
  \begin{aligned}
  T_{zz} &=   -4\pi D_z \overline{\phi} D_z \phi+\pi i   \overline{\psi}_- D_z \psi_- -\pi i   (D_z \overline{\psi}_-)\psi_-  \,,
\\
  T_{\overline{z}\overline{z}} &=   -4\pi D_{\overline{z}} \overline{\phi} D_{\overline{z}} \phi -\pi i \overline{\psi}_+ D_{\overline{z}} \psi_+ + \pi i (D_{\overline{z}} \overline{\psi}_+)\psi_+  \,,
\\
  T_{z\overline{z}} &= \pi |\sigma|^2 |\phi|^2 - \pi D |\phi|^2 + \frac{\pi }{2} \overline{\psi}_+\overline{\sigma}\psi_- + \frac{\pi }{2}\overline{\psi}_- \sigma \psi_+ \,.
  \end{aligned}
\end{equation}
This time we can use the Ward identity \cite{Pol-Ch2}
\begin{equation}
\frac{\partial}{\partial x^\mu} \mathcal O(x) = \frac{1}{2\pi} \int d^2x\, \partial_\nu T^\nu{}_\mu \cdot \mathcal O(x)
\end{equation}
to compute the derivatives $\partial_\mu V_\alpha$.  In terms of $y$, the result is
\begin{equation}\label{Y-Phi-relation-two}
  \partial_z y =4\pi ( D_z  \overline{\phi} )\phi - 2\pi i \overline{\psi}_- \psi_- \,, \qquad
  \partial_{\overline{z}} y =  4\pi   \overline{\phi} D_{\overline{z}} \phi  - 2\pi i \overline{\psi}_+ \psi_+ \,.
\end{equation}

Using $y$ we can also write the smeared vortex operator as
\begin{equation}
V[\rho] = \exp\left( \int d^2x \rho(x) y(x) \right)\,.
\end{equation}

\section{Mirror symmetry}
\label{sec:conn-mirr-symm}

The Hori-Vafa mirror symmetry \cite{Hori:2000kt} specialized to a single chiral multiplet is an abelian duality applied to the phase of the chiral superfield $\Phi = \phi+\theta^-\psi_-+\theta^+\psi_+ +\theta^+ \theta^- F + \ldots$.
The Landau-Ginzburg model mirror to the original theory contains a twisted chiral superfield $Y$, and is characterized by the twisted superpotential
\begin{equation}
  \widetilde{W} = \sigma Y + \mu\, e^{-Y}\,,
\end{equation}
where $\mu$ is a mass scale.
The superfield $Y$ and its conjugate $\overline{Y}$ are related to $\Phi$, its conjugate $\overline{\Phi}$, and the background vector superfield $\mathcal{V}$ as
\begin{equation}
  Y + \overline{Y} = 4 \pi \overline{\Phi} e^{\mathcal{V}} \Phi \,,
\end{equation}
up to an additive renormalization of ${\rm Re}\,y$.
Expansion in $\theta^\mp$ and $\overline{\theta}^\mp$ gives relations among the fields in the dual theories~\cite{Hori:2000kt}.
These include precisely (\ref{Y-Phi-relation-one}) and (\ref{Y-Phi-relation-two}), if $Y$ is related to $(y,\chi_-,\overline{\chi}_+,E)$ in (\ref{V-y}) and (\ref{Y-Phi-relation-one}) as
\begin{equation} \label{Y-y}
Y=y+\overline{\theta}{}^-\chi_- + \theta^+ \overline{\chi}{}_++ \theta^+\overline{\theta}{}^- E + \ldots\,.
\end{equation}
We thus claim that under mirror symmetry, the flavor vortex operator $V_\alpha=e^{-\alpha y}$ and the twisted chiral superfield $Y$ in the mirror theory are related as~(\ref{Y-y}).

The imaginary part of $y$ in the mirror theory is a compact boson with identification ${\rm Im}\, y\sim {\rm Im}\, y+2\pi $.
Thus the exponential $ e^{  - \alpha y}$ for a generic real value of $\alpha$ is ill-defined as a local operator, as we found in the original theory.
Such an exponential operator with non-integer $\alpha$ was considered in \cite{Cecotti:2010fi,Cecotti:2013mba}.

More quantitative results and cases with gauge symmetries will be discussed elsewhere~\cite{with:Kazuo}.

\section*{Acknowledgements}

The author thanks Kazuo Hosomichi for useful discussions.
He is also grateful to the organizers and the participants of the workshop ``Gauge theories, supergravity and superstrings'' at Centro de Ciencias de Benasque Pedro Pascual for providing a stimulating atmosphere.
 This research is supported in part by Grant-in-Aid for Scientific Research~(B)~No.~25287049.


\section{Introduction}

\setcounter{equation}{0}
Given a sequence of positive integers $(x_n)$, which is such that $x_n|x_{n+1}$ for all $n$, the sum of the reciprocals is
the Engel series
\beq \label{engel}
\sum_{j=1}^\infty \frac{1}{x_j} = \sum_{j=1}^\infty \frac{1}{y_1y_2\cdots y_j},
\eeq
where $y_1=x_1$ and $y_{n+1}=x_{n+1}/x_n$ for $n\geq 1$. (It should be assumed that  $(x_n)$ is eventually increasing, which guarantees the convergence of the sum. A brief introduction to Engel series can be found in \cite{duve}.) In recent work \cite{curious}, we considered some particular series of this kind that are generated by certain nonlinear recurrences of second order, and are such that the sequences $(x_n)$ and $(y_n)$ appear interlaced in the continued fraction expansion. Here we start from a sequence with the stronger property that $x_n^2|x_{n+1}$. Any initial 1s can be ignored, but in what follows it will be convenient to start with $x_1=1$ and take $x_2\geq 2$, which implies $x_n\geq 2^{2^{n-2}}$ for $n\geq 2$, and we may write
\beq\label{xdef}
x_1=1, \qquad x_n =\prod_{j=2}^n z_j^{2^{n-j}}
\eeq
for some sequence of positive integers
$(z_n)$ with $z_2\geq 2$. The corresponding Engel series is
\beq \label{engels}
S:=\sum_{j=1}^\infty \frac{1}{x_j} =1+ \sum_{j=2}^\infty \frac{1}{z_2^{2^{j-2}} z_3^{2^{j-3}} \cdots z_j}.
\eeq
Our main result will be to show that in the generic case, when $z_2\geq 3$ and $z_n \geq 2$ for $n\geq 3$, the continued fraction expansion of $S$ has a universal structure,  which we present explicitly.

There is a precedent for these results in the work of Shallit, who first found the continued fraction expansion of the sum
\beq \label{kempner}
\sum_{k=0}^\infty \frac{1}{u^{2^k}}
\eeq
in \cite{sh1}, for integer $u\geq 3$ (with $u=2$ being a degenerate case), and went on \cite{sh3} to obtain the continued fraction for the more general series
\beq \label{shallit}
\sum_{k=0}^\infty \frac{1}{u^{c_k}},
\eeq
where $(c_k)$ is a sequence of positive integers with some non-negative $N$ such that
$
d_n:=c_{n+1}-2c_n\geq 0
$
for all $n\geq N$. If we set $z_2=u^{c_0}$ and $z_j=u^{d_{j-3}}$
for $j\geq 3$ in (\ref{engels}), and assume $d_n\geq 0$ for all
$n\geq 0$, then $S-1$ coincides with (\ref{shallit}).

\subsection{Outline of the paper}

In the next section we prove the main result, namely the expression for the continued fraction expansion of a generic series of the form (\ref{engels}). Section 3 is devoted to an infinite family of examples of series of this type, which are generated by nonlinear recurrences with the Laurent property. For such nonlinear recurrence sequences, we show that the sum of the corresponding series (\ref{engels}) is a transcendental number. In the fourth section we consider the continued fractions obtained from certain degenerate cases, when either
$z_2=2$ or $z_n =1$ for $n\geq 3$, and some particular examples of these degenerate cases are examined in more detail. The final section contains some conclusions.

\section{Continued fractions}
\setcounter{equation}{0}

We use the notation
$$
[a_0; a_1,a_2,a_3,\ldots,a_j,\ldots ]=a_0+\cfrac{1}{a_1+\cfrac{1}{a_2+\cfrac{1}{a_3+\cdots\cfrac{1}{a_j+\cdots}}}}
$$
for continued fractions, and for the
$n$th convergent of the continued fraction $[a_0;a_1,a_2,\ldots]$ we have
$$
\frac{p_n}{q_n} = [a_0;a_1,a_2,\ldots,a_n],
$$
where the numerators $p_n$ and denominators $q_n$ are given in terms of the coefficients $a_j$ according to the matrix identity
\beq\label{matid}
\left(\begin{array}{cc} a_0 & 1 \\ 1 & 0 \end{array}\right)
\left(\begin{array}{cc} a_1 & 1 \\ 1 & 0 \end{array}\right)
\ldots
\left(\begin{array}{cc} a_n & 1 \\ 1 & 0 \end{array}\right) =
\left(\begin{array}{cc} p_n & p_{n-1} \\ q_n & q_{n-1} \end{array}\right).
\eeq
For what follows, it will also be convenient to note the identity obtained by taking the determinant of each side of (\ref{matid}), that is
\beq\label{detid}
p_nq_{n-1}-p_{n-1}q_n = (-1)^{n+1}.
\eeq
For other basic results on continued fractions, the reader is referred to \cite{cassels}.

To begin with, it is convenient to treat the factors $z_2,z_3,\ldots$ of the sequence $(x_n)$ as variables. For the first few partial sums we find the continued fraction expansions
$$
S_1=1,\quad S_2=\frac{z_2+1}{z_2}=[1;z_2],\quad
S_3=\frac{z_2^2z_3+z_2z_3+1}{z_2^2z_3}=[1;z_2-1,1,z_3-1,z_2],
$$
where the
$n$th partial sum of (\ref{engels}) is denoted $S_n$. In general it is straightforward to write $S_n$ as a fraction, that is
\beq \label{numden}
S_n=\frac{\sum_{j=1}^{n-1}\prod_{k=2}^j z_k^{2^{n-k}-2^{j-k}}
\prod_{\ell=j+1}^n z_\ell^{2^{n-\ell}}+1}
{z_2^{2^{n-2}}z_3^{2^{n-3}}\cdots z_n},
\eeq
where the denominator is $x_n$ as given in (\ref{xdef});
but the continued fraction expansion of the $n$th partial sum
is best described recursively.

The basic pattern can be seen by looking at the continued fraction for the fourth partial sum, which  is
$$
S_4=[1;z_2-1,1,z_3-1,z_2,z_4-1,1,z_2-1,z_3-1,1,z_2-1].
$$
Observe that the first five coefficients are the same as those of $S_3$, followed by $z_4-1$, 1, and then four more coefficients which almost coincide with the last four in $S_3$ in reverse order, except that there is $z_2-1$ in place of the final $z_2$ in $S_3$.  This pattern persists, as described by the following

\begin{propn}\label{maincf}
Given the initial set of coefficients
$$
[a_0^{(3)}; a_1^{(3)},a_2^{(3)},a_3^{(3)},a_4^{(3)}]=[1;z_2-1,1,z_3-1,z_2]
$$
of the continued fraction expansion of $S_3$, of length $\ell_3=4$, define a sequence of sets of coefficients $(\{a_j^{(n)}\}_{j=0}^{\ell_{n}-1})$ with
$\ell_n=3\cdot 2^{n-2}-1$ for $n=3,4,5,\ldots$ recursively according to
$$
a_j^{(n+1)}=a_j^{(n)} \quad \mathrm{for}\quad j=0,\ldots, \ell_{n}-1
$$
$$
a_{\ell_{n}}^{(n+1)}=z_{n+1}-1, \qquad
a_{\ell_{n}+1}^{(n+1)}=1, \qquad a_{\ell_{n}+2}^{(n+1)}= a_{\ell_{n}-1}^{(n)}-1,
$$
and
$$
a_j^{(n+1)}=a_{2\ell_{n}-j+1}^{(n)} \quad \mathrm{for}\quad j=\ell_{n}+3,\ldots, 2\ell_{n}.
$$
Then the $n$th partial sum of the series (\ref{engels}) has the  continued fraction expansion
\beq\label{sncf}
S_n=[a_0^{(n)}; a_1^{(n)},\ldots, a_{\ell_n -1}^{(n)}].
\eeq
\end{propn}
\begin{prf} This can be done similarly
to the proof in \cite{sh3}, but we
prefer to use matrix computations, in the same vein as \cite{wu}. The case $n=3$ is easily verified directly. Proceeding by induction, suppose that the continued fraction expansion of $S_n$ is given by (\ref{sncf}), with  coefficients
$a_j^{(n)}$ for $j=0,\ldots \ell_n -1$ with $\ell_n =3\cdot 2^{n-2}-1$ defined according the prescription above, and denote the numerators and denominators of the convergents by $p_j$ and $q_j$ respectively; so the final convergent gives
$$
S_n = \frac{ {p}_{\ell_{n}-1}}{ {q}_{\ell_{n}-1}} , \qquad   {q}_{\ell_{n}-1}=x_n.
$$
Then for the next finite continued fraction defined by this  recursive procedure there are a total of $\ell_{n+1}=2\ell_n +1=3\cdot 2^{n-1}-1$ coefficients, that is $a_j^{(n+1)}$ for $j=0,\ldots \ell_{n+1} -1$, and for the convergents
we use $\tilde{p}_j$,$\tilde{q}_j$ to denote numerators/denominators, respectively. So by  (\ref{matid}) we have
\beq\label{tildes}
\bear{l}
\left(\begin{array}{cc} \tilde{p}_{\ell_{n+1}-1} &  \tilde{p}_{\ell_{n+1}-2} \\  \tilde{q}_{\ell_{n+1}-1} &  \tilde{q}_{\ell_{n+1}-2} \end{array}\right)  =
\bM_n
\left(\begin{array}{cc} z_{n+1}-1 & 1 \\ 1 & 0 \end{array}\right)\left(\begin{array}{cc} 1 & 1 \\ 1 & 0 \end{array}\right)
\tilde{\bM}_n^T
,
\eear
\eeq
where
$$
\bM_n= \left(\begin{array}{cc} a_0^{(n)} & 1 \\ 1 & 0 \end{array}\right)\cdots \left(\begin{array}{cc} a_{\ell_n-1}^{(n)} & 1 \\ 1 & 0 \end{array}\right)=
\left(\begin{array}{cc} p_{\ell_n-1} & p_{\ell_n-2} \\ q_{\ell_n-1} & q_{\ell_n-2} \end{array}\right)
$$
and
$$
\bear{l}
\tilde{\bM}_n^T=\left(\begin{array}{cc} a_{\ell_n-1}^{(n)}-1 & 1 \\ 1 & 0 \end{array}\right)\left(\begin{array}{cc}  a_{\ell_n-2}^{(n)} & 1 \\ 1 & 0 \end{array}\right) \cdots  \left(\begin{array}{cc} a_1^{(n)} & 1 \\ 1 & 0 \end{array}\right) \\
    \qquad  =
\left(\begin{array}{cc} a_{\ell_n-1}^{(n)}-1 & 1 \\ 1 & 0 \end{array}\right) \left(\begin{array}{cc} a_{\ell_n-1}^{(n)} & 1 \\ 1 & 0 \end{array}\right)^{-1}
\bM_n^T  \left(\begin{array}{cc} a_0^{(n)} & 1 \\ 1 & 0 \end{array}\right)^{-1}\\
 \qquad =
  \left(\begin{array}{cc} q_{\ell_n-1}-q_{\ell_n-2} & p_{\ell_n-1}-q_{\ell_n-1}+q_{\ell_n-2}-p_{\ell_n-2} \\ q_{\ell_n-2} & p_{\ell_n-2}-q_{\ell_n-2} \end{array}\right),
\eear
$$
with $T$ denoting transpose. Thus the equation  (\ref{tildes}) simplifies to yield
$$
\bear{l}
\left(\begin{array}{cc} \tilde{p}_{\ell_{n+1}-1} &  \tilde{p}_{\ell_{n+1}-2} \\  \tilde{q}_{\ell_{n+1}-1} &  \tilde{q}_{\ell_{n+1}-2} \end{array}\right)  =
\left(\begin{array}{cc} z_{n+1}q_{\ell_n-1} p_{\ell_n-1}+ 1 &
 z_{n+1} p_{\ell_n-1}\Delta_n -1 \\   z_{n+1} q_{\ell_n-1}^2 & z_{n+1} q_{\ell_n-1}\Delta_n-1\end{array}\right),
\eear
$$
with $\Delta_n= p_{\ell_n-1}- q_{\ell_n-1}$, where we have used the fact that
$$
\det \bM_n = \left|\begin{array}{cc} p_{\ell_n-1} & p_{\ell_n-2} \\ q_{\ell_n-1} & q_{\ell_n-2} \end{array}\right|=-1
$$
by  (\ref{detid}), since $\ell_n$ is odd. Hence we have
$$
\tilde{p}_{\ell_{n+1}-1}= z_{n+1}q_{\ell_n-1} p_{\ell_n-1}+ 1, \qquad  \tilde{q}_{\ell_{n+1}-1}= z_{n+1} q_{\ell_n-1}^2=x_{n+1},
$$
so that
$$
S_{n+1}=S_n+\frac{1}{x_{n+1}}= \frac{ {p}_{\ell_{n}-1}}{ {q}_{\ell_{n}-1}}+\frac{1}{z_{n+1}{q}_{\ell_{n}-1}^2}=\frac{\tilde{p}_{\ell_{n+1}-1}}{\tilde{q}_{\ell_{n+1}-1}}
$$
which is the required result.
\end{prf}
\begin{rem} Note that, mutatis mutandis, both the recursive structure of the partial sums and the above inductive proof hold for the partial sums of an Engel series (\ref{engel}) if, for some positive integer $n_0$, the sequence
$(x_n)$ satisfies the weaker  condition that $x_n^2\vert x_{n+1}$ for
$n\geq n_0$ only. This is the analogue of the fact that for the series
(\ref{shallit}) in \cite{sh3},
$c_{n+1}-2c_n\geq 0$ need only hold
for  $n\geq N$, for some $N$.
\end{rem}
The finite continued fraction expansions of the partial sums immediately yield the continued fraction for the full series (\ref{engels}), at least for a
generic choice of factors $z_2,z_3,\ldots$ of the sequence $(x_n)$.
\begin{thm}\label{main}
For integer factors $z_2\geq 3$ and $z_n\geq 2$ for all $n\geq 3$, the Engel  series (\ref{engels}) has the continued fraction expansion
\beq\label{scf}
S=[a_0; a_1,a_2,\ldots, a_j, \ldots ] = [1; z_2-1,1,z_3-1,z_2, z_4-1,1,\ldots]
\eeq
where the coefficients are given by
$
a_j = \lim_{n\to\infty}a_j^{(n)}.
$
\end{thm}
\begin{prf} The result follows from taking the limit $n\to\infty$ in (\ref{sncf}), provided that none of the coefficients in the finite continued fractions vanish. To see that this is so, note that only $1,z_2-1,z_2$ and $z_3-1$ appear as coefficients in the continued fraction for $S_3$, and all of these are non-zero with the above conditions on the $z_j$. At each step of the recursion in Proposition \ref{maincf} only  $z_{n+1}-1\geq 1$
and $ a_{\ell_{n}+2}^{(n+1)}= a_{\ell_{n}-1}^{(n)}-1$ are potentially new coefficients, so we must check that $ a_{\ell_{n}-1}^{(n)}-1$ cannot vanish.
For $n=3$  the last coefficient is $a_{4}^{(3)}=z_2$, so in $S_4$ this gives  $a_{7}^{(4)}=z_2-1$, while for $n\geq4$ we have
$ a_{\ell_{n}-1}^{(n)}=a_{1}^{(n-1)}=z_2-1$,
so $a_{\ell_{n}+2}^{(n+1)}=z_2-2$, which is why we require  $z_2\geq 3$. Thus the only numbers that appear as
coefficients  in the continued fraction expansion of
$S$ are $1,z_2,z_2-2$ and $z_j-1$ for $j\geq 2$, and none of these are zero.
\end{prf}

\section{Nonlinear recurrence sequences}
\setcounter{equation}{0}

Among nonlinear recurrences of the form
\beq
\label{recf}
x_{n+N} \, x_n = f(x_{n+1}, \ldots ,  x_{n+N-1}),
\eeq
where $f$ is a polynomial in $N-1$ variables,
there is a multitude of examples which surprisingly generate integer sequences.  In a wide variety of cases,  the recurrence (\ref{recf}) has the Laurent
property: for certain special  choices of $f$, all of the iterates belong to the ring $\Z[x_0^{\pm 1},\ldots ,x_{N-1}^{\pm 1}]$;
as a consequence,
if all the initial values are 1 (or $\pm 1$), then each term of the sequence is an integer.
Such sequences were popularized by Gale \cite{gale,gale2}, and subsequently  Fomin and Zelevinsky found  a useful technique -
 the Caterpillar Lemma \cite{fz} - which can be used to prove the   Laurent property in many cases, i.e. for recurrences coming from cluster algebras \cite{fz1}
or in the more general setting of  Laurent Phenomenon (LP) algebras \cite{lp}.

In \cite{honepla,numberpoly} we classified recurrences of second order,
of the form
\beq\label{rec2}
x_{n+2} \, x_n = f(x_{n+1}).
\eeq
For the Laurent property to hold,
the recurrence
(\ref{rec2}) must belong to one of  three classes, depending on the form of $f$: (i) $f(0)\neq 0$, in which case one can apply the framework of cluster algebras (when $f$ is a binomial) or LP algebras (when it is not);
(ii)  $f(0)= 0$, $f'(0)\neq 0$; (iii) $f(0)=f'(0)= 0$.
In the first two classes there are additional requirements on $f$, but in the third class one can take $f(x)=x^2 F(x)$ with arbitrary $F\in\Z[x]$.

In \cite{curious} we considered the case that $f(x)=x^2 F(x)$,  where
$F$ has positive integer coefficients with  $F(0)=1$, and obtained the continued fraction expansion of the sum $\sum_{j=1}^\infty\frac{1}{x_j}$. In order to obtain an Engel series of the form
(\ref{engels}), we should instead choose $F$ so that (\ref{rec2}) becomes
\beq\label{nrec}
x_{n+2} \, x_n = x_{n+1}^{d_1} G(x_{n+1}),
\eeq
where
\beq\label{reqs}
d_1\geq 3, \quad
G(x)\in\Z_{\geq 0}[x], \quad \deg G =d_2\geq 0, \quad G(0)\neq 0, \quad G(1)\geq 3.
\eeq
From
$x_{n+2}/x_{n+1}^2 = \frac{x_{n+1}}{x_n}\cdot x_{n+1}^{d_1-3} G(x_{n+1})$,
we see  by induction that, starting with the initial values $x_0=x_1=1$,
$(x_n)$ is a sequence of positive integers such that $x_n^2|x_{n+1}$
with $x_2=G(1)\geq 3$; hence also $z_n=x_n/x_{n-1}^2\geq 3$ for $n\geq 2$.
\begin{exa}
Taking $d_1=3$, and $G(x)=3$ for all $x$ (so $d_2=0$), the recurrence (\ref{nrec}) becomes
$x_{n+2}x_n = 3x_{n+1}^3$,
which generates the sequence beginning
$
1,1,3,81,531441,5559060566555523,
\ldots
$.
In this case the recurrence can be solved explicitly to yield
$$
x_n =3^{s_n}, \qquad s_n = \frac{5-\sqrt{5}}{10}\left(\frac{3+\sqrt{5}}{2}\right)^n +  \frac{5+\sqrt{5}}{10}\left(\frac{3-\sqrt{5}}{2}\right)^n-1.
$$
The sum of the reciprocals is
$$
S=1+\frac{1}{3}+ \frac{1}{3^4}+\frac{1}{3^{12}}+\frac{1}{3^{33}}+\cdots =[1; 2,1,8,3,80,1,2,8,1,2,19682,\ldots].
$$
\end{exa}
\begin{rem} In the above example, $S-1$ is a sum of the type (\ref{shallit})
considered in \cite{sh3},
with $u=3$ and $c_n =s_{n-2}$. More generally, any recurrence of the form
$x_{n+2}x_n = ux_{n+1}^{d_1}$ with $d_1\geq 3$ and $x_0=x_1=1$ generates a sum of this type.
\end{rem}
For most choices of $G$  it is not possible to give the general solution of the nonlinear recurrence (\ref{nrec}) in closed form.
Nevertheless, one can adapt the methods
of Aho and Sloane \cite{ahos} to write a formula
giving precise asymptotic information. By
rewriting (\ref{nrec}) in terms of logarithms we find that $\La_n=\log x_n$ satisfies
\beq\label{lineq}
\La_{n+1}-(d_1+d_2)\La_{n}+\La_{n-1} =\log c +\al_n, \quad \mathrm{with}\quad
\al_n =\log\left(\frac{G(x_n)}{cx_n^{d_2}}\right),
\eeq
where
$ G(x) = cx^{d_2}+\,\mathrm{lower}\,\, \mathrm{order}.
$
Since $\al_n =\log(1+O(x_n^{-1}))=O(x_n^{-1})$ as $n\to\infty$, the leading order behaviour of $\La_n$ is determined by the linear expression on the left-hand side of (\ref{lineq}),
which has the characteristic  equation $\la^2-(d_1+d_2)\la +1=0$, with
largest root
\beq\label{root}
\la=\frac{d_1+d_2+\sqrt{(d_1+d_2)^2-4}}{2}>2.
\eeq
The next two statements are equivalent to analogous formulae obtained for the sequences considered in \cite{curious}.
\begin{propn} \label{exactf}
For the initial conditions $x_0=x_1=1$, the logarithm $\Lambda_n=\log x_n$ of each term
of the sequence satisfying (\ref{nrec}) is given by the formula
\beq\label{exact}
\La_n = \left(\frac{(1-\la^{-1})\la^n -(1-\la)\la^{-n}}{\la -\la^{-1}}-1\right)
\log c^{-\frac{1}{d_1+d_2-2}}
+\sum_{k=1}^{n-1} \left(\frac{\la^{n-k} -\la^{k-n}}{\la -\la^{-1}}\right)\al_k,
\eeq
where $\al_k$ is defined as in (\ref{lineq}) and $\la$ as in (\ref{root}).
\end{propn}
\begin{cor} \label{cor}
To leading order, the asymptotic approximation of
the logarithm $\La_n$ is given by
\beq\label{asymp}
\La_n\sim C\la^n,
\eeq
where
$$
C=  \frac{1}{d_1+d_2-2}\left(\frac{1-\la^{-1}}{\la -\la^{-1}}\right) \log c +\frac{1}{\la -\la^{-1}}\sum_{k=1}^\infty \la^{-k}\al_k,
$$
and for the terms of the sequence
$
x_n \sim c^{-\frac{1}{d_1+d_2-2}}\exp (C\la^n ).
$
\end{cor}
The asymptotic behaviour of these nonlinear recurrence sequences is enough to show that the sum of the corresponding Engel series is transcendental.
\begin{thm}\label{trans}
Suppose that the sequence $(x_n)$ with initial values $x_0=x_1=1$ is generated by the recurrence (\ref{nrec}) for some $G$ satisfying the conditions (\ref{reqs}). Then the sum $S$ in  (\ref{engels}) is a transcendental number.
\end{thm}
\begin{prf} This is essentially identical to the proof of Theorem 4 in \cite{curious},  so here we only sketch the argument.  Recall that
Roth's theorem says that if $\al$ is an irrational algebraic number then for an arbitrary fixed $\delta>0$ there are only
finitely many rational approximations $p/q$ for which
\beq\label{roth}
\left|\al -\frac{p}{q}\right|<\frac{1}{q^{2+\delta}}.
\eeq
The number $S$ has an infinite continued fraction expansion, so it is irrational.
From the asymptotics (\ref{asymp})  it follows that
for any $\epsilon>0$ the growth condition
\beq\label{growth}
x_{n+1}>x_n^{ \la-\epsilon}
\eeq
 holds for all sufficiently large $n$. By making a comparison with a geometric sum, this  gives
$$
\left|S-\frac{p_{\ell_n-1}}{q_{\ell_n-1}}\right|=\sum_{j=n+1}^\infty \frac{1}{x_j} <
\frac{1}{x_{n}^{\la-\epsilon-\epsilon'}}=\frac{1}{q_{\ell_n-1}^{\la-\epsilon-\epsilon'}}
$$
for any $\epsilon'>0$ and $n$  large enough.
 So if $\epsilon$ and $\epsilon'$ are chosen such that $\la-\epsilon-\epsilon'=2+\delta >2$, then $\al=S$ has infinitely many
rational  approximations satisfying (\ref{roth}), and hence must be transcendental.
\end{prf}
\begin{exa}
Taking $d_1=3$, and $G(x)=2x+1$ for all $x$ (so $d_2=1$), the recurrence (\ref{nrec}) becomes
\beq \label{nex}
x_{n+2}x_n = x_{n+1}^3(2x_{n+1}+1),
\eeq
which generates the sequence
$$
1,1,3,189,852910317,
5599917937724687764238078261637795,
\ldots
$$
with leading order asymptotics
$$
x_n \sim \frac{e^{ C (2+\sqrt{3})^n}}{\sqrt{2}}, \qquad C\approx 0.107812043.
$$
The sum of the reciprocals is the  transcendental number
$$
S =[1; 2,1,20,3,23876,1,2,20,1,2,7697947188058154,\ldots]\approx 1.3386243
$$
\end{exa}

Engel series of the form (\ref{engels}) can also be generated by nonlinear recurrences of higher order. For instance, one can take a recurrence of third order,
\beq\label{third}
X_{n+3}X_n =X_{n+1}^{e_1}X_{n+2}^{e_2}H(X_{n+1},X_{n+2}),\quad\mathrm{with}\,\, e_1\geq 1,\,e_2\geq 2,
\eeq
where the polynomial $H(X,Y)\in\Z_{\geq 0}[X,Y]$ is not divisible by either of its arguments. It is straightforward
to show that
the Laurent property holds for this recurrence, and from
$X_{n+3}/X_{n+2}^2 =\frac{X_{n+1}}{X_{n}}\cdot X_{n+1}^{e_1-1}X_{n+2}^{e_2-2}H(X_{n+1},X_{n+2})$,
it is easy to see by induction that the initial values $X_0=X_1=X_2=1$ generate an integer sequence with
$X_n^2|X_{n+1}$ for all $n\geq 0$. Thus the sum of reciprocals starting from the index 2, that is
\beq\label{sp}
S':=\sum_{j=2}^\infty \frac{1}{X_j},
\eeq
 is an Engel series of the form (\ref{engels}). Note that the condition $H(1,1)\geq 3$ should be imposed, in order for Theorem \ref{main} to apply to this series.

A particular class of recurrences of the form (\ref{third}) can be obtained by factorizing the terms of a sequence satisfying (\ref{nrec}) as
$x_n =X_n X_{n+1}$,
which lifts the second order recurrence to
\beq\label{thirds}
X_{n+3}X_n =(X_{n+1}X_{n+2})^{d_1-1}G(X_{n+1}X_{n+2}).
\eeq
For a generic polynomial $H(X_{n+1},X_{n+2})$ on the right-hand-side of (\ref{third}), it is not immediately
obvious which term
will be dominant as $n\to\infty$, but in the special case (\ref{thirds}) the same techniques as for the second order recurrence can be applied directly, to show that the  leading order asymptotics is
$
\log X_n \sim C' \la^n
$
for some $C'>0$, where $\la$ is given by (\ref{root}). This means that Theorem \ref{trans} applies to the series (\ref{sp})
as well.
\begin{exa}
Setting $x_n=X_nX_{n+1}$ in (\ref{nex}) gives the recurrence
\beq \label{nexlift}
X_{n+3}X_n = X_{n+1}^2X_{n+2}^2(2X_{n+1}X_{n+2}+1),
\eeq
which generates the sequence beginning
$$
1,1,1,3,63,13538259,
413636490314204194515563505,
\ldots .
$$
To leading order,
$
\log X_n \sim C' (2+\sqrt{3})^n$, with $C'\approx 0.0227833$.
The sum of the reciprocals in (\ref{sp}) is  the  transcendental number
$$
S' =[1; 2,1,6,3,3410,1,2,6,1,2,2256800700104,\ldots] \approx  1.3492064.
$$
\end{exa}
For other examples of transcendental numbers whose complete continued fraction expansion is known, see \cite{ds} and references.

\section{Degenerate cases}
\setcounter{equation}{0}

If either $z_2=2$ or  $z_n=1$ for some $n\geq3$, then one of the coefficients in the continued fraction
becomes zero, and Theorem \ref{main} is no longer valid. To
obtain a continued fraction
with non-zero coefficients, one can use the replacement rule
$
[\ldots, a,0,b, \ldots ] \longrightarrow [\ldots, a+b,\ldots]
$
(see Proposition 3 in \cite{vdpsh}) to remove the zero. Each such replacement, decreases the length of a finite continued
fraction  by two, so in degenerate cases the length of the continued fraction expansion of
$S_n$ is typically shorter than the generic value $\ell_n=3\cdot 2^{n-2}-1$. Here we present the expansion for two particular
degenerate cases, omitting details of the proof.

\subsection{The case $z_2=2$}

For generic values of the factors $z_j$, the sequence of lengths $\ell_n$ of partial sums begins $1,2,5,11,23$ for $n=1,2,3,4,5$. When $z_2=2$ and  $z_j\geq2$ for all $j\geq3$, the first few continued fractions for the partial sums of $S$ are $S_1=1$,
$$
S_2=[1;2],\quad
S_3=[1;1,1,z_3-1,2], \quad
S_4=[1;1,1,z_3-1,2,z_4-1,1,1,z_3-1,2],
$$
which are of the same length as in the generic case, except for $S_4$ being of length 10, since at the end $[\ldots, 1,1]\to
[\ldots, 2]$. The first zero appears in $S_5$, which contains  a single coefficient $z_2-2$, so removing this
and making the final replacement $[\ldots, 1,1]\to [\ldots, 2]$,
results in the length being 20:
$$
S_5=[1;1,1,z_3-1,2,z_4-1,1,1,z_3-1,1,1,z_5-1,2,z_3-1,1,1,z_4-1,2,z_3-1,2].
$$
Thereafter the pattern continues with the continued fraction doubling in length at each step, as described by the
following
\begin{thm}\label{z2}
When $z_2=2$ and $z_j\neq 1$ for $j\geq 3$,
the  Engel series (\ref{engels}) has the continued fraction expansion
$$
S=[a_0; a_1,a_2,\ldots, a_j, \ldots ] = [1; 1,1,z_3-1,2, z_4-1,1,\ldots]
$$
with coefficients given by
$
a_j = \lim_{n\to\infty}a_j^{(n)},
$
where  $(\{a_j^{(n)}\}_{j=0}^{\ell_{n}-1})$, the  sequence of sets of coefficients of
the finite continued fractions for partial sums $S_n$, of length
$\ell_n=5\cdot 2^{n-3}$ for $n=4,5,\ldots$, is defined  by starting from
$$
[a_0^{(4)}; a_1^{(4)},\ldots,a_9^{(4)}]=[1;1,1,z_3-1,2,z_4-1,1,1,z_3-1,2],
$$
and obtaining subsequent coefficients according to
$$
a_j^{(n+1)}=a_j^{(n)} \quad \mathrm{for}\quad j=0,\ldots, \ell_{n}-2
$$
$$
a_{\ell_{n}-1}^{(n+1)}=1,\qquad a_{\ell_{n}}^{(n+1)}=1,\qquad
a_{\ell_{n}+1}^{(n+1)}=z_{n+1}-1,
$$
$$
a_j^{(n+1)}=a_{2\ell_{n}-j+1}^{(n)} \quad \mathrm{for}\quad j=\ell_{n}+2,\ldots, 2\ell_{n}-2,\qquad\mathrm{and}\quad a_{2\ell_{n}-1}^{(n+1)}=2.
$$
\end{thm}
In order to obtain a sequence $(x_n)$ of this degenerate type from a second order recurrence of the form (\ref{nrec}), the
conditions (\ref{reqs}) should be modified so that $G(1)=2$, which requires that
$G(x)=x^{d_2}+1$ for some non-negative integer $d_2$. So the recurrence becomes
\beq\label{mrec}
x_{n+2} \, x_n = x_{n+1}^{d_1} (x_{n+1}^{d_2}+1),
\eeq
with $d_1\geq 3$, $d_2\geq 0$. The results of Proposition \ref{exactf}, Corollary \ref{cor} and Theorem \ref{trans}
all apply without alteration to  sequences obtained from (\ref{mrec}).
\begin{exa}
Taking $d_1=3$, and $G(x)=x+1$ for all $x$, the recurrence (\ref{nrec}) becomes
$x_{n+2}x_n = x_{n+1}^3(x_{n+1}+1)$,
which generates the sequence
$$
1,1,2,24,172800,37150633525248000000,
\ldots
$$
with asymptotics
$x_n \sim e^{C(2+\sqrt{3})^n}$,
$C\approx 0.06224548$.
The sum of the reciprocals is the transcendental number
$$
S =[1;1,1,5, 2,299,1,1,5,1,1,1244167199,2,5,1,1,299,\ldots]\approx  1.54167245.
$$
\end{exa}

\subsection{The case $z_j=1$ for $j\geq 3$}

If we set $z_2=u$ and all other factors $z_j=1$ then $x_n=u^{2^{n-2}}$ for $n\geq 2$; the expansion (\ref{scf}) is no longer valid because each coefficient
$z_j-1$ becomes zero for $j\geq 3$. In that case, the sum of the reciprocals is
$S=1+\sum_{k=0}^\infty u^{-2^k}$, so that $S-1$
coincides with (\ref{kempner}).
The continued fractions for the partial sums were first obtained in \cite{sh1}, and a nonrecursive description was given in \cite{sh2}.
The sequence of lengths begins 1,2,3,5,9,17, with
$\ell_n=2^{n-2}+1$ for $n\geq 3$, and the full continued fraction is
$$
S=[1;u-1,u+2,u,u,u-2,u,u+2,u,u-2,u+2,u,u-2,u,u,u+2,u,\ldots]
$$
for $u\neq 2$.
The only numbers that appear as coefficients in this continued fraction
are $1,u-2,u-1,u,u+2$.

However, the case $u=2$ is special, since some of the coefficients in the
above expansion become zero.
The sequence of lengths of partial sums starts with 1,2,3,5,7,11, and
$\ell_n=2^{n-3}+3$ for $n\geq 4$. Only the numbers 1,2,4,6 appear as coefficients in
the continued fraction for the series, which is
$$
S=[1;1,4,2,4,4,6,4,2,4,6,2,4,6,4,4,2,4,6,\ldots].
$$

The argument used to prove Theorem \ref{trans},
based on Roth's theorem,
does not apply to the partial sums of the series  (\ref{kempner}),
since the sequence $(x_n)$ does not grow fast enough: in contrast to
(\ref{growth}), it satisfies the recurrence $x_{n+1}=x_n^2$.
However,
a direct proof of transcendence of (\ref{kempner}),
valid for all integers $u\geq 2$, was first given by Kempner \cite{kempner}; various alternative proofs are collected in \cite{adamc}.

\section{Conclusions}

We have found the continued fraction expansion for an Engel series of the special type (\ref{engels}). In some cases, coming from nonlinear recurrence sequences, it has been shown that this produces transcendental numbers. We expect that the sum (\ref{engels}) should be transcendental for any choice of the factors $z_2\geq 2$ and $z_j\geq 1$ for $j\geq 3$.
However, we do not know of a simple way to prove this  in general.

\vspace{.05in}
\noindent
{\bf Acknowledgments:}
This work is supported by Fellowship EP/M004333/1
from the EPSRC.
The original inspiration came from Paul Hanna's  observations concerning the nonlinear recurrence sequences
described in  \cite{curious}, which were communicated via the Seqfan mailing list.
The author is grateful to Jeffrey Shallit for pointing out some of his work, as well as Kempner's proof \cite{kempner}.

\small