\section{Introduction}

Let $\overbar{\mathcal{M}}_{g,n}$ be the moduli space of stable algebraic curves of genus~$g$ with~$n$ distinct marked points,
where $g$ and~$n$ are non-negative integers satisfying the stability condition $2g-2+n>0$.
For $1\le k\le n$ and $0\le j\le g$,
denote by~$\psi_k$ the first Chern class of the $k$-th tautological line bundle $\mathbb{L}_k$ on $\overbar{\mathcal{M}}_{g,n}$, and by~$\lambda_j$ the $j$-th Chern class of the Hodge bundle $\mathbb{E}_{g,n}$ on~$\overbar{\mathcal{M}}_{g,n}$.
The rational numbers defined by the formula
$$
\int_{\overbar{\mathcal{M}}_{g,n}}
\psi_1^{i_1}\cdots\psi_n^{i_n}\lambda_{j_1}\cdots\lambda_{j_m}
$$
are called the Hodge integrals. These numbers take zero value unless the degree-dimension counting matches: \begin{equation} \label{dd} i_1+\dots +i_n+ j_1+\dots +j_m = 3g-3+n. \end{equation}
Denote by $\mathcal{C}_g(z):= \sum_{j=0}^g \lambda_j z^j$ the Chern polynomial of~$\mathbb{E}_{g,n}$.
We will be particularly interested in the following class of Hodge integrals
defined via cubic products of Chern polynomials, called the cubic Hodge integrals:
\begin{equation}\label{cubHo}
\int_{\overbar{\mathcal{M}}_{g,n}}
\psi_1^{i_1}\cdots\psi_n^{i_n} \mathcal{C}_g(-p)\mathcal{C}_g(-q)\mathcal{C}_g(-r),
\end{equation}
where $p,q,r$ are complex parameters.
These integrals are called \emph{special} if $p,q,r$ satisfy the following
{\it local Calabi-Yau condition}:
\begin{equation}\label{calabi-yau}
pq+qr+rp=0.
\end{equation}
These Hodge integrals are important in the localization technique of computing Gromov-Witten invariants for toric three-folds \cite{GP, LLZ, OP}. Their significance was also manifested by the
Gopakumar-Mari\~no-Vafa conjecture regarding the Chern-Simons/string duality \cite{GV,MV}.

Let $\H=\mathcal{H}({\bf t}; p,q,r; \epsilon)$ be the cubic Hodge free energy defined by
\begin{align}
& \mathcal{H}({\bf t}; p,q,r; \epsilon) := \sum_{g\geq0}\epsilon^{2g-2}\mathcal{H}_g({\bf t}; p,q,r)  ,  \nonumber \\
& \mathcal{H}_g(\mathbf{t}; p,q,r) := \sum_{n\geq 0} \sum_{i_1,\dots,i_n\geq 0} \frac{t_{i_1}\cdots t_{i_n}}{n!}
 \int_{\overbar{\mathcal{M}}_{g,n}}  \psi_1^{i_1}\cdots\psi_n^{i_n} \mathcal{C}_g(-p)\mathcal{C}_g(-q)\mathcal{C}_g(-r) . \nonumber
\end{align}
Here $\mathcal{H}_g(\mathbf{t}; p,q,r)$ is called the genus~$g$ part of the free energy~$\H$. Then the exponential
\[e^{\mathcal{H}({\bf t};p,q,r;\epsilon)} =:Z_{\rm cubic}({\bf t};p,q,r;\epsilon)=:Z_{\rm cubic}\] is called the cubic Hodge partition function.
Clearly, $\H_g({\bf t}; p,q,r)\in \mathbb{C}[p,q,r][[{\bf t}]]$.
The genus zero free energy~$\H_0({\bf t})$ is actually independent of~$p,q,r$ and has the explicit expression
$$
\H_0({\bf t}) = \sum_{n\geq 3} \frac{1}{n(n-1)(n-2)} \sum_{i_1+\dots+i_n=n-3} \frac{t_{i_1}}{i_1!} \cdots \frac{t_{i_n}}{i_n!} .
$$

Define
\begin{equation}
v({\bf t}):=\partial_{t_0}^2 \H_0({\bf t})=\sum_{n\geq 1} \frac1n  \sum_{i_1+\dots+i_n=n-1} \frac{t_{i_1}}{i_1!} \cdots \frac{t_{i_n}}{i_n!},\label{vexpression}
\end{equation}
then it satisfies the following Riemann hierarchy:
\begin{equation}\label{kdv-genus-zero}
\frac{\partial v}{\partial t_i} = \frac{v^i}{i!} \frac{\partial v}{\partial t_0} , \quad i\geq 0.
\end{equation}
More generally, if one defines $w=\epsilon^2 \partial_{t_0}^2 \H({\bf t};p,q,r;\epsilon)$, then $w$ satisfies an integrable
hierarchy of Hamiltonian evolutionary PDEs \cite{DLYZ-1}, called the (cubic) Hodge hierarchy,
which is a deformation of the Riemann hierarchy. The first member of this integrable hierarchy reads
$$
w_{t_1} = w w_{t_0} + \frac{\epsilon^2}{12} \left(w_{t_0t_0t_0}  - (p+q+r) w_{t_0} w_{t_0t_0} \right) + {\rm O}\left(\epsilon^4\right).
$$

The Hodge-FVH correspondence is given by the following conjecture~\cite{LZZ}:
\begin{cnj}
The Hodge hierarchy for the special cubic Hodge integrals is equivalent, under a certain Miura type transformation,
to the {\it fractional Volterra hierarchy} (FVH). Furthermore, the corresponding cubic Hodge partition function
gives a tau function of the FVH.
\end{cnj}

For the case with $p=q$, the
 validity of the Hodge-FVH correspondence is
 implied by the Hodge-GUE correspondence established in~\cite{DLYZ-2, DY}.
The goal of this and the subsequent papers~\cite{LYZZ}
is to prove the Hodge-FVH correspondence.
In the present paper, we derive the Dubrovin-Zhang loop equation by studying the
Virasoro constraints for the special cubic Hodge partition function.
We show that this loop equation together with the genus zero free energy
uniquely determines the partition function.
In the subsequent paper, we will show that the fractional Volterra hierarchy admits the same
Virasoro constraints and that there exists a particular tau function of the integrable hierarchy
which is uniquely determined by the same loop equation, and we prove in this way the
Hodge-FVH correspondence.

From now on, we assume that $p,q,r$ satisfy the local Calabi-Yau condition~\eqref{calabi-yau}.
The case with $p,q,r\in \mathbb{Q}$ is called {\it rational}. Note that
the Virasoro constraints for the special cubic Hodge partition function in general case are quite complicated.
However, in the rational case, we find explicit expressions
of the Virasoro constraints which lead to the Dubrovin-Zhang loop
equation. It turns out that the loop equation for the general case can be deduced from
the one for the rational case.

Let us first consider the rational case. Due to the symmetry property and the homogeneity property
of the cubic Hodge integrals with respect to $p,q,r$,
we can assume that
\begin{equation}\label{rational-condition}
p=\frac{1}{K_1},\quad q=\frac{1}{K_2},
\quad r=-\frac{1}{h},
\end{equation}
where $K_1, K_2 \in \mathbb{N}$, $(K_1,K_2)=1$ and $h:=K_1+K_2$.
We denote, for $\ell\geq 0$,
\begin{align}
& b_{\alpha+h\ell} := \frac{\alpha}{K_1}+\ell, \quad c_{\alpha+h\ell} := \binom{b_{\alpha+h\ell} h}{b_{\alpha+h\ell} K_1} , \qquad \alpha=0,\dots,K_1-1,\label{bcdef1}\\
& b_{\alpha+h\ell} := \frac{-\alpha}{K_2}+\ell, \quad c_{\alpha+h\ell} := \binom{b_{\alpha+h\ell} h}{b_{\alpha+h\ell} K_2}, \qquad \alpha=-(K_2-1),\dots,-1,\label{bcdef2}
\end{align}
and
\[ \mathbb{N}_*=(\mathbb{N}-K_2) \backslash \left(\{0\} \cup (h\mathbb{N}-K_2)\right), \quad \mbox{where } a\mathbb{N}-K_2:=\{ak-K_2|k\in \mathbb{N}\}.\]
Define
\begin{equation}\label{Zxs}
Z(x,{\bf s};\epsilon):=\exp\left(\frac{A(x,\tilde {\bf s})}{\epsilon^2}\right)
Z_{\rm cubic}\left(\mathbf{t}(x,{\bf s}); \tfrac{1}{K_1}, \tfrac{1}{K_2}, -\tfrac{1}{h};\epsilon \right),
\end{equation}
where ${\bf s}:=(s_k)_{k\in \mathbb{N}_*}$ is an infinite vector of indeterminates,
$\tilde{s}_k=s_k-c_h^{-1}\delta_{k,h}$~($k\in \mathbb{N}_*$),
\begin{equation}\label{t-s}
t_i=t_i(x,{\bf s})=\sum_{k\in \mathbb{N}_*}  b_k^{i+1} c_k \tilde s_k+\delta_{i,1}+x \delta_{i,0},\quad i\geq 0,
\end{equation}
and~$A$ is the quadratic series
\begin{equation}
A:=A(x,{\bf s})=\frac{1}{2}\sum_{k_1,k_2\in \mathbb{N}_*} \frac{b_{k_1} b_{k_2}}{b_{k_1} + b_{k_2}} c_{k_1} c_{k_2} s_{k_1} s_{k_2} + x \sum_{k\in \mathbb{N}_*} c_k s_k.
\end{equation}
Note that for $g\geq0$,
$\mathcal{H}_g\Bigl(\mathbf{t}(x,\mathbf{s});\frac{1}{K_1}, \frac{1}{K_2}, -\frac{1}{h}\Bigr)$
is a well-defined formal power
series in $\mathbb{C}[[x-1,\mathbf{s}]]$.
Indeed, for each monomial $(x-1)^{k_0}s_{k_1}\cdots s_{k_m}$ in
$\mathcal{H}_g\Bigl(\mathbf{t}(x,\mathbf{s});\frac{1}{K_1}, \frac{1}{K_2}, -\frac{1}{h}\Bigr)$ only finitely many
 monomials $t_{i_1}\cdots t_{i_n}$ contribute to its coefficient for the dimension reason~\eqref{dd}.

Denote $I=\{-(K_2-1),\dots,K_1-1\}$ and $I_*=I\backslash \{0\}$, and define a family of linear operators $L_m=L_m\left(\epsilon^{-1}x,\epsilon^{-1}\mathbf{s},\epsilon\partial/\partial \mathbf{s}\right)$, $m\geq 0$ by
\begin{align}
& L_0=\sum_{k\in \mathbb{N}_*} b_k s_k \frac{\partial}{\partial s_k}+\frac{x^2}{2\epsilon^2}+\frac{1}{24}\left(\frac1h-\frac1{K_1}-\frac1{K_2}\right), \label{virasoro-target0} \\
& L_m=\sum_{k\in \mathbb{N}_*} b_k s_k \frac{\partial}{\partial s_{k+h m}}+x\frac{\partial}{\partial s_{hm}} \nonumber\\
& \qquad ~ +\frac{\epsilon^2}{2} \sum_{\ell=1}^{m-1}\frac{\partial^2}{\partial s_{h \ell} \partial s_{h(m-\ell)}} + \frac{\epsilon^2}{2}  \sum_{\alpha,\beta \in I_*} \sum_{\ell=0}^{m-1}  G^{\alpha\beta}  \frac{\partial^2}{\partial s_{\alpha+h\ell} \partial s_{\beta+h(m-1-\ell)}},\label{virasoro-target}
\end{align}
where $\left(G^{\alpha\beta}\right)_{\alpha,\beta\in I}$ is a symmetric non-degenerate constant matrix defined by
\begin{align}\label{Galphabeta}
G^{\alpha\beta}=
\begin{cases}
\frac{K_1}{h}\delta^{\alpha+\beta,-K_2} &\quad \alpha,\beta<0;\\
\quad 1                                       &\quad \alpha=\beta=0;\\
\frac{K_2}{h}\delta^{\alpha+\beta, K_1} &\quad \alpha,\beta>0;\\
\quad 0                                       &\quad \text{elsewhere}.
\end{cases}
\end{align}
It is easy to check that the operators $L_m$ satisfy the following Virasoro commutation relations:
$$
\bigl[ L_m,L_n \bigr]=(m-n)L_{m+n}, \quad \forall\, m,n\geq 0.
$$
\begin{thm}\label{thm-virasoro}
The series $Z(x,\mathbf{s};\epsilon)$ defined by~\eqref{Zxs} satisfies the Virasoro constraints
\begin{equation}\label{VirasoroZ}
L_m\left(\epsilon^{-1}x,\epsilon^{-1}\tilde{\mathbf{s}},\epsilon\partial/\partial \mathbf{s}\right)
Z(x,\mathbf{s};\epsilon)=0, \quad m\geq 0.
\end{equation}
\end{thm}
Using Theorem~\ref{thm-virasoro} and the technique developed in~\cite{DZ},
we derive in Section~\ref{sec-loop} the Dubrovin-Zhang loop equation for the special cubic Hodge
free energies in the rational case; see Theorem~\ref{main}.

We proceed  to the general case. Denote
\begin{equation}\label{sigma12}
\sigma_1 = -(p+q+r), \quad \sigma_3 = -2 \left(p^3+q^3+r^3\right).
\end{equation}
From Mumford's relations~\cite{M}, the local Calabi-Yau condition~\eqref{calabi-yau},
the fact that the integral~\eqref{cubHo} is symmetric in $p,q,r$,
and the well-known relationship between the Chern polynomial and the Chern character of~$\mathbb{E}_{g,n}$, we know that
\begin{equation}\label{polynomiality}
\H_g:=\H_g({\bf t};p,q,r)\in \mathbb{C}[\sigma_1,\sigma_3][[{\bf t}]], \quad g\ge 0.
\end{equation}
The following theorem is the main result of the present paper.
\begin{thm} \label{refineloop}
The equation
\begin{align}
&\sum_{i\geq0 }\left( \partial^i \Theta+\sum_{j=1}^i P_{j-1,i-j+1}\right) \frac{\partial \Delta H }{\partial z_i} \nonumber\\
= & \frac{\Theta^2}{16}-\left(\frac{1}{16}-\frac{\sigma_1}{24}\right)\Theta+
\epsilon^2\sum_{i\geq0}\partial^{i+2}\left(\frac{\Theta^2}{16}-\left(\frac{1}{16}-\frac{\sigma_1}{24}\right)\Theta\right)
\frac{\partial \Delta H}{\partial z_i}\nonumber\\
& + \frac{\epsilon^2}{2}\sum_{i,j\geq0} P_{i+1,j+1} \left(\frac{\partial^2\Delta H}{\partial z_i \partial z_j}
+\frac{\partial \Delta H}{\partial z_i}\frac{\partial \Delta H}{\partial z_j}\right),\label{loop-eq-0}
\end{align}
which is called the Dubrovin-Zhang loop equation for the special cubic Hodge integrals,
has a unique solution of the form
\[\Delta H:=\sum_{g\geq 1} \epsilon^{2g-2} H_g,\quad H_g:=H_g(z_0,\dots,z_{3g-2};\sigma_1,\sigma_3)\]
up to the addition of a constant to each $H_g$, $g\geq 1$.
These constants can be uniquely determined by the following conditions:
\begin{align}
& H_1= \frac1{24} \log z_1 + \frac{\sigma_1}{24} z_0, \quad
\sum_{j=1}^{3g-2} j z_j \frac{\partial H_g}{\partial z_j} = (2g-2)H_g, \quad g\geq 2.  \label{eq2233}
\end{align}
Here
\[\partial:=\sum_{k\geq 0} z_{k+1} \frac{\partial}{\partial z_k},\quad \Theta:= \frac1{1-e^{z_0}/\mu},\]
$P_{i,j}$ are certain polynomials in
$\Theta, z_1,z_2,\dots, \sigma_1,\sigma_3$
which are given in Section~\ref{section4}, and $\mu$ is an arbitrary parameter. Moreover,
let $v({\bf t})$ be defined in~\eqref{vexpression}, then the genus~$g$~($g\geq 1$) special cubic
Hodge free energy has the expression
\begin{equation}\label{HgHg}
\mathcal{H}_g
=H_g\left(v({\bf t}), \frac{\partial v({\bf t})}{\partial t_0},\cdots,\frac{\partial^{3g-2} v({\bf t})}{\partial t_0^{3g-2}};\sigma_1,\sigma_3\right).
\end{equation}
\end{thm}

We can recursively solve the loop equation to obtain the free energies $H_g$, $g\geq 1$.
For example, the first two $H_g$ are given by
\begin{align}
H_1=&\frac{1}{24}\log z_1+\frac{\sigma_1}{24} z_0, \label{explicitH1} \\
H_2
=&\frac{1}{1152 }\frac{z_4}{z_1^2} -\frac{7 }{1920 }\frac{z_2 z_3}{z_1^3}+\frac{1}{360}\frac{z_2^3}{z_1^4}
+\frac{\sigma_1 }{480 }\frac{z_3}{z_1}-\frac{11\sigma_1 }{5760}\frac{z_2^2}{z_1^2}
+\frac{7\sigma_1^2 }{5760}z_2\nonumber\\
&+\left(\frac{\sigma_1^3}{17280}-\frac{\sigma_2}{34560}\right)z_1^2. \label{explicitH2}
\end{align}
These expressions agree with the results of~\cite{DLYZ-1}.

Besides the above theorems, we also prove the following gap phenomenon for
the special cubic Hodge integrals in the rational case. This type of phenomenon was used
by M.-X.~Huang, A.~Klemm and S.~Quackenbush to compute the Gromov-Witten invariants
of the quintic Calabi-Yau three-fold up to genus~$51$ \cite{HKQ}.
\begin{cor}\label{gapcor} For the rational case, the free energies
$\mathcal{F}_g(x,{\bf s}):=\H_g\Bigl({\bf t}(x,{\bf s});\frac{1}{K_1}, \frac{1}{K_2}, -\frac{1}{h}\Bigr)$, $g\geq 1$ with ${\bf t}(x,{\bf s})$ defined in~\eqref{t-s}
have the following gap phenomenon:
\begin{align}
& \mathcal{F}_1(x,{\bf s}) - \frac{\sigma_1-1}{24} \log x  \in \mathbb{C}[x][[{\bf s}]],  \label{gap1}\\
& \mathcal{F}_g (x,{\bf s}) - \frac{R_g(\sigma_1,\sigma_3)}{x^{2g-2}} \in
\mathbb{C}[x] [[{\bf s}]], \quad g\geq 2. \label{gap2}
\end{align}
Here $R_g(\sigma_1,\sigma_3)$, $g\geq 2$ are certain polynomials of~$\sigma_1,\sigma_3$
satisfying the condition
\begin{equation}
\deg R_g\leq 3g-3,\quad \textrm{with}\  \deg \sigma_1=1, \deg \sigma_3=3,\label{zh-10}
\end{equation}
whose degree $(3g-3)$-parts are given by
\begin{equation} \label{Faber}
   \frac{(-1)^g}{2(2g-2)!}  \frac{|B_{2g}||B_{2g-2}|}{2g (2g-2)} \left(\frac{\sigma_1^3}3-\frac{\sigma_3}{6}\right)^{g-1},
\end{equation}
where $B_j$ denote the Bernoulli numbers.
\end{cor}

Note that in~\eqref{gap1} and \eqref{gap2}, the indeterminates of $R_g(\sigma_1,\sigma_3)$ take the values
\[\sigma_1=\frac1{K_1+K_2}-\frac1{K_1}-\frac1{K_2},\qquad
\sigma_3=\frac2{(K_1+K_2)^3}- \frac2{K_1^3} - \frac{2}{K_2^3}.\]
The first two $R_g$, $g\geq 2$ have the expressions
\begin{align*}
& R_2=-\frac{1}{1440}+\frac{13\sigma_1}{5760}-\frac{7 \sigma_1^2}{5760}+\frac{\sigma_1^3}{17280}-\frac{\sigma_3}{34560},\\
& R_3=\frac{1}{181440}-\frac{107 \sigma_1}{362880}+\frac{145\sigma_1^2}{290304}-\frac{961 \sigma_1^3}{4354560}
+\frac{31 \sigma_3}{2177280}+\frac{113 \sigma_1^4}{4354560}-\frac{113\sigma_1\sigma_3}{8709120}\\
& \qquad -\frac{\sigma_1^6}{13063680}+\frac{\sigma_1^3\sigma_3}{13063680}-\frac{\sigma_3^2}{52254720}.
\end{align*}
We hope that the above results will be useful in the study of Gromov-Witten invariants for toric Calabi-Yau varieties.

\paragraph{Organization of the paper}
In Section~2, we derive the explicit expressions of Virasoro constraints for~$Z(x,{\bf s}; \epsilon)$.
In Section~3, we obtain the Dubrovin-Zhang loop equation for the special cubic Hodge free energies in the rational case.
In Section~4, we prove Theorem~\ref{refineloop}.

\paragraph{Acknowledgements}
We would like to thank Boris Dubrovin and Don Zagier for several very helpful discussions.
This work is partially supported by NSFC No.\,11771238.
The work of S.-Q.~Liu is also supported by the National Science Fund for Distinguished Young Scholars No.\,11725104. Part of the work of C.~Zhou was done while he was a Ph.D. student of Tsinghua University; he acknowledges Tsinghua University for excellent working conditions and financial supports.

\section{Virasoro constraints}
In this section we first give two versions of Virasoro constraints for~$Z_{\rm cubic}$,
and then we prove Theorem~\ref{thm-virasoro}.
Denote by~$Z_{\rm WK}$ the Witten-Kontsevich partition function
$$
Z_{\rm WK}(\mathbf{t};\epsilon)=\exp
\left(\sum_{g\geq0}\epsilon^{2g-2}\sum_n
\frac{t_{i_1}\cdots t_{i_n}}{n!}
\int_{\overbar{\mathcal{M}}_{g,n}}\psi_1^{i_1}\cdots\psi_n^{i_n}\right).
$$
It is well known that $Z_{\rm WK}(\mathbf{t};\epsilon)$ satisfies the following Virasoro constraints~\cite{DVV,Witten,Kontsevich}:
\begin{equation}\label{vir-kdv}
L_m^{\rm KdV}\left(\epsilon^{-1}\tilde{\mathbf{t}},\epsilon\partial/\partial \mathbf{t}\right) Z_{\rm WK}(\mathbf{t};\epsilon)=0,\quad
\end{equation}
where $\tilde{t}_i=t_i-\delta_{i,1}$, and $L_m^{\rm KdV}:=L_m^{\rm KdV}\left(\epsilon^{-1}\mathbf{t},\epsilon\partial/\partial \mathbf{t}\right)$ are linear operators given by
\begin{align}
& L_{-1}^{\rm KdV} :=\sum_{i\geq1}t_i\frac{\partial}{\partial t_{i-1}}+\frac{t_0^2}{2\epsilon^2}, \nonumber \\
& L_0^{\rm KdV}  :=\sum_{i\geq0}\frac{2i+1}{2} t_i\frac{\partial}{\partial t_i}+\frac{1}{16}, \nonumber \\
& L_m^{\rm KdV} :=\sum_{i\geq0}\frac{(2i+2m+1)!!}{2^{m+1}(2i-1)!!} t_i\frac{\partial}{\partial t_{i+m}} +\frac{\epsilon^2}{2}\sum_{i+j=m-1}\frac{(2i+1)!!(2j+1)!!}{2^{m+1}}\frac{\partial^2}{\partial t_i\partial t_j},  ~ m\ge 1. \nonumber
\end{align}
These operators satisfy the Virasoro commutation relations
$$
\left[L_m^{\rm KdV},L_n^{KdV}\right]=(m-n)L_{m+n}^{KdV},\quad \forall\,m,n\geq -1.
$$

\begin{lem}[\cite{FP}] \label{lemmaFP}
The partition function~$Z_{\rm cubic}(\mathbf{t};p,q,r;\epsilon)$ has the expression
\begin{equation}\label{Z-cubic}
Z_{\rm cubic}(\mathbf{t};p,q,r;\epsilon)
=e^{G(\epsilon^{-1}\tilde{\mathbf{t}},\epsilon\partial/\partial \mathbf{t})} Z_{\rm WK}(\mathbf{t};\epsilon),
\end{equation}
where the operator~$G=G\left(\epsilon^{-1}\mathbf{t},\epsilon\partial/\partial \mathbf{t}\right)$ is defined by
$$
G\left(\epsilon^{-1}\mathbf{t},\epsilon\partial/\partial \mathbf{t}\right) = - \sum_{i\geq1}\frac{B_{2i}}{2i(2i-1)}\left(p^{2i-1}+q^{2i-1}+r^{2i-1}\right)D_i,
$$
$B_{2i}$ are Bernoulli numbers,  and $D_i=D_i\left(\epsilon^{-1}\mathbf{t},\epsilon\partial/\partial \mathbf{t}\right)$ are given by
$$
D_i := - \sum_{j\geq0}t_j\frac{\partial}{\partial t_{j+2i-1}}
+\frac{\epsilon^2}{2}\sum_{j=0}^{2i-2}\frac{\partial^2}{\partial t_j \partial t_{2i-2-j}}, \quad i\geq 1.
$$
\end{lem}

The lemma below follows from equation~\eqref{vir-kdv} and Lemma~\ref{lemmaFP}.

\begin{lem}[\cite{Z}] \label{lemmazhou}
Define a sequence of operators $L_m^{\rm cubic}=L_m^{\rm cubic}\left(\epsilon^{-1}\mathbf{t},\epsilon \partial / \partial \mathbf{t}\right)$ by
$$
L_m^{\rm cubic}\left(\epsilon^{-1}\mathbf{t},\epsilon \partial / \partial \mathbf{t}\right) := e^G \circ L_m^{\rm KdV}\left(\epsilon^{-1}\mathbf{t},\epsilon \partial / \partial \mathbf{t}\right) \circ e^{-G}, \quad m\geq -1.
$$
Then $L_m^{\rm cubic}$ satisfy the Virasoro commutation relations
$$
\Bigl[L_m^{\rm cubic},L_n^{\rm cubic}\Bigr]=(m-n)L_{m+n}^{\rm cubic},\quad m,n\geq -1.
$$
Moreover, $Z_{\rm cubic}(\mathbf{t};p,q,r;\epsilon)$ satisfies the equations
\begin{equation}\label{Virasorocubic}
L_m^{\rm cubic}\left(\epsilon^{-1}\tilde{\mathbf{t}},\epsilon\partial / \partial \mathbf{t}\right) Z_{\rm cubic}(\mathbf{t};p,q,r;\epsilon)=0,
\quad m\geq -1.
\end{equation}
\end{lem}
\noindent We call~\eqref{Virasorocubic} the Virasoro constraints for~$Z_{\rm cubic}(\mathbf{t};p,q,r;\epsilon)$.
Note that these Virasoro operators $L_m^{\rm cubic}$ are in general quite complicated, and it is difficult
to use~\eqref{Virasorocubic} directly for the computation of the cubic Hodge integrals. For example,
the operator $L_0^{\rm cubic}$ has the explicit expression
$$
L_0^{\rm cubic} = L_0^{\rm KdV} - \sum_{k\geq 1} \frac{B_{2k}}{2k}  \frac{\sigma_{2k-1}}{(2k-2)!} D_k,
$$
where $\sigma_{2k-1}=-(2k-2)! \left(p^{2k-1}+q^{2k-1}+r^{2k-1}\right)$,
and the corresponding constraint given in~\eqref{Virasorocubic} has the expression
\begin{align}
& \left(\frac32\frac{\partial}{\partial t_1} + \sum_{k\geq 1} \frac{B_{2k}}{2k} \frac{\sigma_{2k-1}}{(2k-2)!}  \frac{\partial}{\partial t_{2k}} \right)
Z_{\rm cubic} \nonumber\\
& = \left(\sum_{i\geq0}\frac{2i+1}{2} t_i\frac{\partial}{\partial t_{i}} +\frac1{16}
+ \sum_{k\geq 1} \frac{B_{2k}}{2k}  \frac{\sigma_{2k-1}}{(2k-2)!}
\left(\sum_{j\geq0}t_j\frac{\partial}{\partial t_{j+2k-1}}
 -
\sum_{j=0}^{2k-2}\frac{\partial^2}{\partial t_j \partial t_{2k-2-j}}\right)\right)
Z_{\rm cubic}. \nonumber
\end{align}

Following~\cite{DLYZ-2} we consider the following linear combinations of Virasoro operators:
\begin{equation}\label{vir-til}
\widetilde{L}_m^{\rm cubic}:=\sum_{k\geq-1}\frac{m^{k+1}}{(k+1)!}L_k^{\rm cubic},\quad m\geq 0.
\end{equation}
As it is shown in~\cite{DLYZ-2}, the operators $\widetilde{L}_m^{\rm cubic}$ also
satisfy the Virasoro commutation relations
$$
\left[\widetilde{L}_m^{\rm cubic},\widetilde{L}_n^{\rm cubic}\right]=(m-n)\widetilde{L}_{m+n}^{\rm cubic},
\quad \forall\, m,n\geq 0.
$$
From~\eqref{Virasorocubic} it follows that
\begin{equation}\label{Virasorocubic2nd}
\widetilde{L}_m^{\rm cubic}\left(\epsilon^{-1}\tilde{\mathbf{t}},\epsilon \partial / \partial \mathbf{t}\right) Z_{\rm cubic}(\mathbf{t};p,q,r;\epsilon)=0,
\quad m\geq 0.
\end{equation}
\noindent We call~\eqref{Virasorocubic2nd}
the second version of Virasoro constraints for~$Z_{\rm cubic}(\mathbf{t};p,q,r;\epsilon)$.

Let us proceed to derive the explicit expressions of~$\widetilde{L}_m^{\rm cubic}$.
To this end, we are to use the Givental quantization technique~\cite{G} to simplify the computation.
Let $(\mathcal{V},\omega)$ denote the Givental symplectic space, where $\mathcal{V}$ is the space of Laurent series over~$\mathbb{C}$,
and $\omega$ is the bilinear form on~$\mathcal{V}$ defined by
\begin{equation}
\omega(f,g):=-\Res_{z=\infty}f(-z)g(z)\frac{dz}{z^2},\quad \forall\, f,g\in\mathcal{V}.
\end{equation}
Denote by $q_i, p_i$~$(i\geq 0)$ the Darboux coordinates associated to $\omega$, and
define a family of infinitesimal symplectic transformations $\l_k$, $k\geq -1$ on~$\mathcal{V}$ by
\begin{equation}
l_k:=(-1)^{k+1}z^{3/2}\partial_z^{k+1} z^{-1/2},\quad k\geq -1.
\end{equation}
Then we have the following formulae~\cite{G}:
\begin{equation}\label{GiventalKdVlm}
L_k^{\rm KdV}
=\widehat{l_k}\big|_{q_i\mapsto t_i,\partial_{q_i}\mapsto\partial_{t_i}, \, i\geq 0} +\frac{\delta_{k,0}}{16}.
\end{equation}
Similarly, define a family of infinitesimal symplectic transformations $d_i$, $i\geq 1$ on~$\mathcal{V}$ by
\[d_i:=z^{1-2i},\]
then we have
\begin{equation}\label{diDi}
D_i \left(\epsilon^{-1}\mathbf{t},\epsilon\partial/\partial\mathbf{t}\right) =\widehat{d_i}\big|_{q_j\mapsto t_j,\partial_{q_j}\mapsto\partial_{t_j}, \, j\ge0}.
\end{equation}
Note that the above results of this section hold true for general $p,q,r$.
The local Calabi-Yau condition~\eqref{calabi-yau} for $p, q, r$ will be needed in what follows.

Similarly as in~\cite{DLYZ-2} we have the following lemma.
\begin{lem}\label{phi-expan}
Assume that $p,q,r$ satisfy the local Calabi-Yau condition~\eqref{calabi-yau}.
Define
\begin{equation}
\Psi(z)=\left(p^{1/p}q^{1/q}r^{1/r}\right)^{-z}
\sqrt{\frac{\Gamma(1-z/p)\Gamma(1-z/q)\Gamma(1-z/r)}
{\Gamma(1+z/p)\Gamma(1+z/q)\Gamma(1+z/r)}} .
\end{equation}
Then the asymptotic expansion of~$\log\Psi(z)$ as $z\rightarrow\infty$ within a properly chosen sector
is given by (up to an inessential constant $\pm \pi i/4$ )
\begin{equation}
\log \Psi(z) \sim -\sum_{i=1}^{\infty}\frac{B_{2i}}{2i(2i-1)} \left(p^{2i-1}+q^{2i-1}+r^{2i-1}\right)z^{1-2i} =: \log\Phi(z).
\end{equation}
\end{lem}

\begin{lem}\label{quantization} The operators~$\widetilde{L}_m^{\rm cubic}$, $m\geq 0$ have the expressions
\begin{align}
\widetilde{L}_m^{\rm cubic}\left(\epsilon^{-1}\mathbf{t},\epsilon\partial/\partial\mathbf{t}\right)
= \left(zV_m(z)e^{-m\partial_z}\right)^{\wedge} \Big|_{q_i\mapsto t_i,\partial_{q_i}\mapsto\partial_{t_i},\,i\geq 0}+\frac{m}{16}-\frac{p+q+r}{24}, \label{quantizationgivental}
\end{align}
where $V_m$, $m\geq 0$ are given by
\begin{equation}\label{v}
V_m(z)=\frac{\Phi(z)}{\Phi(z-m)}\sqrt{\frac{z}{z-m}}\,.
\end{equation}
\end{lem}
\begin{prf}
Denote
\[\widehat{\Phi}:=\exp\left(\left(\log\Phi(z)\right)^{\wedge}\big|_{q_i\mapsto t_i,\partial_{q_i}\mapsto\partial_{t_i},\,i\geq 0}\right).\]
It follows from~\eqref{diDi} and
Lemma~\ref{lemmazhou} that
\[L_m^{\rm cubic}\left(\epsilon^{-1}\mathbf{t},\epsilon\partial/\partial \mathbf{t}\right)
=\widehat{\Phi}\circ L_m^{\rm KdV}\left(\epsilon^{-1}\mathbf{t},\epsilon \partial/\partial \mathbf{t}\right)\circ \widehat{\Phi}^{-1}, \quad m\geq 0.\]
Then by using~\eqref{GiventalKdVlm} we obtain
\begin{align*}
L_m^{\rm cubic}\left(\epsilon^{-1}\mathbf{t},\epsilon\partial/\partial \mathbf{t}\right)
=&\left(\Phi\circ l_m \circ \Phi^{-1}\right)^{\wedge} \Big|_{q_i\mapsto t_i,\partial_{q_i}\mapsto \partial_{t_i}}+\frac{1}{16}\delta_{m,0}-\frac{p+q+r}{24}\delta_{m,-1}.
\end{align*}
Note that the term~$-\frac{p+q+r}{24}\delta_{m,-1}$ comes from the cocycle~\cite{G} in the quantization of~$\Phi\circ l_m \circ \Phi^{-1}$.
The lemma follows from~\eqref{vir-til} as well as the identity
\[z^{3/2}\Phi(z)\circ e^{-m\partial_z}\circ \frac{1}{\sqrt{z}\Phi(z)}=z\frac{\Phi(z)}{\Phi(z-m)}\sqrt{\frac{z}{z-m}} \, e^{-m\partial_z}.\]
\end{prf}

We are now ready to prove Theorem~\ref{thm-virasoro}.

\begin{prfn}{Theorem~\ref{thm-virasoro}}
In order to prove the validity of~\eqref{VirasoroZ}, it suffices to show that
the following identities hold true:
\begin{equation}\label{Lmconjugate}
L_m\left(\epsilon^{-1}x,\epsilon^{-1}\tilde{\mathbf{s}},\epsilon \partial / \partial \mathbf{s}\right)
=K^m e^{\frac{A(x,\tilde{{\bf s}})}{\epsilon^2}}\circ\widetilde{L}_m^{\rm cubic}\left(\epsilon^{-1}\tilde{\mathbf{t}},\epsilon\partial/\partial\mathbf{t}\right)\circ e^{-\frac{A(x,\tilde {\bf s})}{\epsilon^2}}, \quad m\ge 0,
\end{equation}
where
\begin{equation}\label{zh-9}
K=h^h K_1^{-K_1} K_2^{-K_2},
\end{equation}
 and
\begin{equation}\label{ttttt}
\tilde t_i=\tilde  t_i(x,{\bf s})=\sum_{k\in \mathbb{N}_*}  b_k^{i+1} c_k \tilde s_k+x \delta_{i,0},
\quad i\ge 0.
\end{equation}
For simplicity we denote
\[\bar{s}_k:=c_k s_k,\quad k\in \mathbb{N}_*.\]
Then the above substitution of the independent variables~\eqref{ttttt} reads
\begin{equation} \label{t-s-0}
\tilde t_i (x,{\bf s}) =  \sum_{k\in \mathbb{N}_*}b_k^{i+1}\tilde{\bar{s}}_k+x \delta_{i,0}, \quad \tilde{\bar{s}}_k = \bar{s}_k - \delta_{k,h}.
\end{equation}
Since both the left and the right hand sides of~\eqref{Lmconjugate} satisfy the Virasoro commutation relations,
we only need to prove~\eqref{Lmconjugate} for $m=0,1,2$.

\begin{lem}
The operator $\widetilde{L}_0^{\rm cubic}$ satisfies the following relation:
\begin{equation}
e^{\frac{A(x,\tilde{{\bf s}})}{\epsilon^2}}\circ \widetilde{L}_0^{\rm cubic} \left(\epsilon^{-1}\tilde{\mathbf{t}},\epsilon\partial/\partial\mathbf{t}\right) \circ e^{-\frac{A(x,\tilde{{\bf s}})}{\epsilon^2}}
=\sum_{k\in \mathbb{N}_*}b_k \tilde{\bar{s}}_k \frac{\partial}{\partial \bar{s}_k}+\frac{x^2}{2\epsilon^2}+\frac{\sigma_1}{24}.
\end{equation}
\end{lem}
\begin{prf}
By using Lemma~\ref{quantization} we have
\begin{align*}
\widetilde{L}_0^{\rm cubic} \left(\epsilon^{-1}\mathbf{t},\epsilon\partial/\partial\mathbf{t}\right)
=\widehat{z} \mid_{q_i\mapsto t_i,\partial_{q_i}\mapsto\partial_{t_i}}+\frac{\sigma_1}{24}
=\sum_{i\geq1}t_i\frac{\partial}{\partial t_{i-1}}+\frac{t_0^2}{2\epsilon^2}+\frac{\sigma_1}{24}\,.
\end{align*}
Under the substitution~\eqref{t-s-0}, we arrive at
\begin{align*}
&\widetilde{L}_0^{\rm cubic}\left(\epsilon^{-1}\tilde{\mathbf{t}},\epsilon\partial/\partial\mathbf{t}\right)=\sum_{i\geq 1}\left(\sum_{k\in \mathbb{N}_*}b_k^{i+1}\tilde{\bar{s}}_k\right)\frac{\partial}{\partial t_{i-1}}
+\frac{1}{2\epsilon^2}\left(\sum_{k\in \mathbb{N}_*} b_k\tilde{\bar{s}}_k +x\right)^2
+\frac{\sigma_1}{24}\\
=&\sum_{k\in \mathbb{N}_*}b_k \tilde{\bar{s}}_k \frac{\partial}{\partial \bar{s}_k}
+\frac{1}{2\epsilon^2}\left(\sum_{k\in \mathbb{N}_*} b_k\tilde{\bar{s}}_k +x\right)^2
+\frac{\sigma_1}{24}.
\end{align*}
The lemma is proved by applying the conjugation by~$e^{\frac{A(x,\tilde{{\bf s}})}{\epsilon^2}}$ to this equality.
\end{prf}

\begin{lem} The operator $\widetilde{L}_1^{\rm cubic}$ satisfies the following relation:
\begin{align}
& e^{\frac{A(x,\tilde{{\bf s}})}{\epsilon^2}}\circ
\widetilde{L}_1^{\rm cubic}\left(\epsilon^{-1}\tilde{\mathbf{t}},\epsilon\partial/\partial\mathbf{t}\right)\circ e^{-\frac{A(x,\tilde{{\bf s}})}{\epsilon^2}}\nonumber\\
& ~ =\sum_{k\in \mathbb{N}_*}b_k V_1(-b_k) \tilde{\bar{s}}_k \frac{\partial}{\partial \bar{s}_{k+h}}
+x V_1(0) \frac{\partial}{\partial \bar{s}_h}
+\frac{\epsilon^2}{2}\sum_{\alpha=1}^{K_1-1}
\frac{\Res_{z=b_{\alpha}}V_1(z)}{b_{\alpha}}
\frac{\partial^2}{\partial \bar{s}_{\alpha} \partial \bar{s}_{K_1-\alpha}}\nonumber\\
&  \qquad +\frac{\epsilon^2}{2}\sum_{\alpha=-(K_2-1)}^{-1}
\frac{\Res_{z=b_{\alpha}}V_1(z)}{b_{\alpha}}
\frac{\partial^2}{\partial \bar{s}_{\alpha} \partial \bar{s}_{-\alpha-K_2}}, \label{compare1}
\end{align}
where $V_1(z)$ is defined in~\eqref{v}.
\end{lem}

\begin{prf}
Since $p=1/K_1,q=1/K_2,r=-1/h$, the function~$V_1(z)$ becomes rational. By employing
Lemma~\ref{phi-expan} we find the following explicit expression of this rational function:
$$
V_1(z)=\frac{\prod_{i=1}^{K_1+K_2}\left(z-\frac{i}{K_1+K_2}\right)}
{\prod_{i=1}^{K_1}\left(z-\frac{i}{K_1}\right)\prod_{i=1}^{K_2}\left(z-\frac{i}{K_2}\right)}\,.
$$
Note that $V_1(z)$ has simple poles at $b_{-(K_2-1)},\dots,b_{-1}$, $b_1,\dots,b_{K_1-1}$, and $1$,
therefore,
\[V_1(z)=1 +\frac{\Res_{y=1}V_1(y)}{z-1} +\sum_{\alpha\in I_*} \frac{\Res_{y=b_{\alpha}}V_1(y)}{z-b_{\alpha}}.\]
As $z$ goes to infinity, the function $V_1(z)$ has the full asymptotic expansion
$$
V_1(z) \sim
1+\sum_{n\geq1}z^{-n} \left(\Res_{y=1}V_1(y)  + \sum_{\alpha\in I_*} b_{\alpha}^{n-1}\Res_{y=b_{\alpha}}V_1(y)
\right).
$$
Denote
$a_{\alpha}:=b_{\alpha}^{-1} \Res_{z=b_{\alpha}}V_1(z)$, $\alpha\in I_*$.
By taking $m=1$ in~\eqref{quantizationgivental} we obtain
\begin{align*}
&\widetilde{L}_1^{\rm cubic}  = \left(zV_1(z)e^{-\partial_z}\right)^{\wedge}\Big|_{q_i\mapsto t_i,\,\partial_{q_i}\mapsto \partial_{t_i}}
+\frac{1}{16}+\frac{\sigma_1}{24}\\
&=\sum_{i\geq0}\left(\sum_{j=0}^{i+1} t_j\binom{i+1}{j}
+\sum_{j=0}^i t_j\sum_{n=1}^{i+1-j}(-1)^n\binom{i+1}{j+n} \left(\Res_{y=1}V_1(y) + \sum_{\alpha\in I_*} a_{\alpha}b_{\alpha}^n\right)\right)
\frac{\partial}{\partial t_i}\\
& \quad +\frac{\epsilon^2}{2}\sum_{i,j\geq0}\left(\sum_{n=0}^{i+1} (-1)^n \binom{i+1}{n} \left(\Res_{y=1}V_1(y)+ \sum_{\alpha\in I_*} a_{\alpha} b_{\alpha}^{n+j+1}\right)\right)
\frac{\partial^2}{\partial t_i \partial t_j}
+\frac{t_0^2}{2\epsilon^2}+\frac{1}{16}+\frac{\sigma_1}{24}\,.
\end{align*}
Now by performing the variables substitution~\eqref{t-s-0} we arrive at
\begin{align*}
&\widetilde{L}_1^{\rm cubic}\left(\epsilon^{-1}\tilde{\mathbf{t}},\epsilon\partial/\partial\mathbf{t}\right) = \sum_{k\in \mathbb{N}_*}b_k V_1(-b_k) \tilde{\bar{s}}_k \frac{\partial}{\partial \bar{s}_{k+h}}
+x V_1(0) \frac{\partial}{\partial \bar{s}_h}\\
&+\frac{\epsilon^2}{2}\sum_{\alpha=1}^{K_1-1}a_{\alpha} \frac{\partial^2}{\partial \bar{s}_{\alpha} \partial \bar{s}_{K_1-\alpha}}
+\frac{\epsilon^2}{2}\sum_{\alpha=-(K_2-1)}^{-1} a_{\alpha} \frac{\partial^2}{\partial \bar{s}_{\alpha} \partial \bar{s}_{-\alpha-K_2}}\\
&+\sum_{\alpha=1}^{K_1-1}a_{\alpha}\frac{\partial A(x,\tilde {\bf s})}{\partial \bar{s}_{\alpha}}\frac{\partial}{\partial \bar{s}_{K_1-\alpha}}
 +\sum_{\alpha=-(K_2-1)}^{-1}a_{\alpha}\frac{\partial A(x,\tilde {\bf s})}{\partial \bar{s}_{\alpha}}\frac{\partial}{\partial \bar{s}_{-\alpha-K_2}} \nonumber\\
& +\frac{1}{2\epsilon^2}\left(\sum_{k\in \mathbb{N}_*} b_k \tilde{\bar{s}}_k+x\right)^2 +\frac{1}{16}+\frac{\sigma_1}{24}.
\end{align*}
Here we have used the identity
\[\sum_{j=0}^i\sum_{k=1}^{i+1-j}\binom{i+1}{j+k}\lambda^j \mu^k
=\mu\frac{(\lambda+1)^{i+1}-(\mu+1)^{i+1}}{\lambda-\mu}.\]
Dressing the operator $\widetilde{L}_1^{\rm cubic}\left(\epsilon^{-1}\tilde{\mathbf{t}},\epsilon\partial/\partial\mathbf{t}\right)$
by~$e^{\frac{A(x,\tilde {\bf s})}{\epsilon^2}}$ we obtain~\eqref{compare1} after a long but straightforward computation.
The lemma is proved.
\end{prf}

We can prove in a similar way the following lemma.
\begin{lem}
The operator $\widetilde{L}_2^{\rm cubic}$ satisfies the relation
\begin{align}
& e^{\frac{A(x,\tilde{{\bf s}})}{\epsilon^2}}\circ
\widetilde{L}_2^{\rm cubic}\left(\epsilon^{-1}\tilde{\mathbf{t}},\epsilon\partial/\partial\mathbf{t}\right)\circ e^{-\frac{A(x,\tilde{{\bf s}})}{\epsilon^2}} \nonumber \\
& =\sum_{k\in \mathbb{N}_*}b_k V_2(-b_k) \tilde{\bar{s}}_k \frac{\partial}{\partial \bar{s}_{k+2h}}
+x V_2(0) \frac{\partial}{\partial \bar{s}_{2h}} +\frac{\epsilon^2}{2}b_h^{-1}\Res_{z=b_h}V_2(z)\frac{\partial^2}{\partial \bar{s}_h \partial \bar{s}_h} \nonumber \\
& \qquad +\frac{\epsilon^2}{2}\sum_{\alpha=1}^{K_1-1}\sum_{\ell=0}^1 b_{\alpha}^{-1} \Res_{z=b_{\alpha}}V_2(z)
\frac{\partial^2}{\partial \bar{s}_{\alpha+h \ell} \partial \bar{s}_{K_1-\alpha+h(1-\ell)}} \nonumber \\
& \qquad +\frac{\epsilon^2}{2}\sum_{\alpha=-(K_2-1)}^{-1}\sum_{\ell=0}^1 b_{\alpha}^{-1} \Res_{z=b_{\alpha}}V_2(z)
\frac{\partial^2}{\partial \bar{s}_{\alpha+h \ell} \partial \bar{s}_{-\alpha-K_2+h(1-\ell)}}. \label{compare2}
\end{align}
\end{lem}

\begin{lem}\label{const}
The numbers $c_k$, $k\in \mathbb{N}_*$ defined in \eqref{bcdef1}--\eqref{bcdef2} have the following properties:
\begin{enumerate}[(i)]
\item
For $k\in \mathbb{N}_*$ and $\ell\geq 1$, $c_{k+h \ell}/c_k=K^{\ell} V_{\ell}(-b_k)$.
\item For $\ell\geq 1$, $c_{h \ell}=K^{\ell} V_{\ell}(0)$.
\item   Let $m,n\geq0$ be integers. Then
\begin{align}
& c_{\alpha+h m} \, c_{K_1-\alpha+ h n} = \frac{h}{K_2} \frac{K^{m+n+1}}{b_{\alpha+h m}}  \Res_{z=b_{\alpha+h m}}V_{m+n+1}(z), \quad \alpha=1,\dots,K_1-1,  \nonumber\\
& c_{\alpha+h m} \, c_{-\alpha-K_2+ h n} = \frac{h}{K_1} \frac{K^{m+n+1}}{b_{\alpha+h m}}  \Res_{z=b_{\alpha+h m}}V_{m+n+1}(z), \quad \alpha=-(K_2-1),\dots,-1,\nonumber \\
& c_{h (m+1)} \, c_{h (n+1)}=\frac{K^{m+n+2}}{b_{h (m+1)}}\Res_{z=h (m+1)}V_{m+n+2}(z).   \nonumber
\end{align}
\end{enumerate}
\end{lem}
The proof of Lemma~\ref{const} is elementary with the help of Lemma~\ref{phi-expan}, so we omit it here.

By comparing \eqref{compare1}--\eqref{compare2} with the definition of the left hand side of~\eqref{Lmconjugate} and by using Lemma~\ref{const}
we find that~\eqref{Lmconjugate} is true for $m=1, 2$. Thus Theorem~\ref{thm-virasoro} is proved.
\end{prfn}

\section{Loop equation: the rational case} \label{sec-loop}
In this section, we derive the Dubrovin-Zhang loop equation for the special cubic Hodge integrals in the rational case, namely, we take
\[p=1/K_1,\quad q=1/K_2,\quad r=-1/h,\]
where $K_1$ and $K_2$ are coprime positive integers, and $h=K_1+K_2$.

Introduce a generating series of the operators~$L_m$, $m\geq0$ defined in \eqref{virasoro-target0},\,\eqref{virasoro-target} as follows:
\begin{equation}\label{Ltwo}
L(\lambda)=\sum_{m\geq 0}\frac{L_m}{\lambda^{m+1}}.
\end{equation}
\begin{lem}\label{virasoro-total}
The generating series~$L(\lambda)$ can be represented as
\begin{equation}
L(\lambda)=\biggl[J_1(\lambda)J_2(\lambda)+\frac{1}{2}J_2(\lambda)^2\biggr]_{{\rm reg},-}
+\frac{1}{\lambda}\left(\frac{x}{2\epsilon^2}+\frac{\sigma_1}{24}\right),
\end{equation}
where the operators $J_1(\lambda)$ and $J_2(\lambda)$ are defined by
\begin{align}
& J_1(\lambda) := \epsilon^{-1}\frac{x}{\sqrt{\lambda}} + \epsilon^{-1} \sum_{k \in \mathbb{N}_*}b_k \lambda^{b_k-\frac{1}{2}}s_k, \nonumber\\
& J_2(\lambda):=  \sum_{\ell\ge 1}\frac{\epsilon}{\lambda^{b_{h\ell}+1/2}}\frac{\partial}{\partial s_{h \ell}}
+ \sqrt{\frac{K_2}{h}}\sum_{\alpha=1}^{K_1-1} \sum_{\ell\geq0}\frac{\epsilon}{\lambda^{b_{\alpha+h\ell}+\frac{1}{2}}}\frac{\partial}{\partial s_{\alpha+h\ell}} \nonumber\\
&\qquad \qquad \qquad + \sqrt{\frac{K_1}{h}}\sum_{\alpha=-(K_2-1)}^{-1}\sum_{\ell\geq0}\frac{\epsilon}{\lambda^{b_{\alpha+h\ell}+\frac{1}{2}}}\frac{\partial}{\partial s_{\alpha+h\ell}}, \nonumber
\end{align}
and $[\,\cdot\,]_{{\rm reg},-}$ means keeping terms of the series with negative integer powers of~$\lambda$.
\end{lem}

Theorem~\ref{thm-virasoro} now implies that $L(\lambda)Z(x,\mathbf{s};\epsilon)=0$ holds true {\it identically} in~$\lambda$.
From the definition of~$Z(x,\mathbf{s};\epsilon)$ given by~\eqref{Zxs}
we know that its logarithm admits the genus expansion
\[\log Z(x,\mathbf{s};\epsilon)=\sum_{g\geq0}\epsilon^{2g-2}\mathcal{F}_g(x,\mathbf{s}).\]
Here we denote, as we do in Corollary~\ref{gapcor},
\begin{equation}\label{zh-1}
\mathcal{F}_g(x, {\bf s}) :=
\mathcal{H}_g\Bigl({\bf t}(x,{\bf s});\frac{1}{K_1}, \frac{1}{K_2}, -\frac{1}{h}\Bigr)+A(x,\tilde{\bf s}) \, \delta_{g,0},\quad \textrm{with}\ \tilde s_k=s_k-c_h^{-1}\delta_{k,h}.
\end{equation}
In what follows we will also use the notations
\[ u(x,\mathbf{s}):=\partial_x^2\mathcal{F}_0(x,\mathbf{s})\]
and
\begin{equation}\label{zh-5}
\Delta \mathcal{F}=\sum_{g\geq1}\epsilon^{2g-2}\mathcal{F}_g(x,\mathbf{s}).
\end{equation}
\begin{lem}\label{prp-1}
The function $u(x, \mathbf{s})$ satisfies the following equations
\begin{align}
&\partial_x^i u(x,\mathbf{s})=\partial^i_{t_0}v(\mathbf{t}(x,\mathbf{s})), \quad i\ge 0.\\
&\frac{\partial u}{\partial s_k}=c_k \partial_x e^{b_k u},
\quad \frac{\partial^2\mathcal{F}_0}{\partial x \partial s_k}=c_k e^{b_k u},
\quad k \in \mathbb{N}_*.\label{dispersionlesshierarchy}
\end{align}
\end{lem}
\begin{prf}
By using~\eqref{t-s} we obtain
$$
\partial_x=\partial_{t_0},\quad
\frac{\partial}{\partial s_k}=c_k \sum_{i\geq0}b_k^{i+1}\frac{\partial}{\partial t_i},\quad k\in\mathbb{N}_*.
$$
The lemma then follows from the Riemann hierarchy \eqref{kdv-genus-zero}.
\end{prf}

\begin{lem}[\cite{DLYZ-1}]\label{jetvariablesrep}
For any given $g\ge 1$, there exists a polynomial $H_g(z_0,\dots,z_{3g-2};p,q,r)$ of $p, q, r, z_2, z_3,\dots$ with coefficients depending smoothly on $z_0$ and rationally on~$z_1$, such that
\begin{equation}
\mathcal{H}_g(\mathbf{t};p,q,r)=
H_g\left(v({\bf t}), \frac{\partial v({\bf t})}{\partial t_0},\dots,\frac{\partial^{3g-2} v({\bf t})}{\partial t_0^{3g-2}};p,q,r\right). \label{HHgHg}
\end{equation}
\end{lem}
It should be noted that Lemma~\ref{jetvariablesrep} does not require $p,q,r$ satisfy the local Calabi-Yau condition.

Let us denote
\[u^{(i)}=\partial_x^i u(x,{\bf s}),\quad i\geq 0.\]
Then from~Lemmas~\ref{prp-1}--\ref{jetvariablesrep} we know that
\begin{equation}\label{zh-6}
\mathcal{F}_g(x,{\bf s})=H_g\left(u^{(0)}, u^{(1)},\dots,u^{(3g-2)}; p,q,r\right),\quad
g\geq 1,
\end{equation}
where $p=1/K_1$, $q=1/K_2$, $r=-1/(K_1+K_2)$.
Introduce a derivation~$D(\lambda)$ on $\mathbb{C}[[x-1,\mathbf{s}]]$ by
$$
D(\lambda)=
-\left[J_1(\lambda)J_2(\lambda)\right]_{{\rm reg},-}- \epsilon^{-2} \left[J_2(\lambda)(\Delta\mathcal{F})J_2(\lambda)\right]_{{\rm reg},-},
$$
where $J_1(\lambda), J_2(\lambda)$ are defined in Lemma \ref{virasoro-total}.

\begin{lem}\label{lemma35}
The series $D(\lambda)(\Delta\mathcal{F})$ has the following expression:
\begin{align}
&D(\lambda)(\Delta\mathcal{F})=\sum_{i\geq0} D(\lambda) \left(u^{(i)}\right)\frac{\partial \Delta\mathcal{F}}{\partial u^{(i)}}
=\frac{\sigma_1}{24\lambda} + \biggl[\frac{1}{2}\sum_{i\geq 0}J_2(\lambda)^2\left(u^{(i)}\right)\frac{\partial \Delta\mathcal{F}}{\partial u^{(i)}} \nonumber \\
&\ +\frac{1}{2}\sum_{i,j\geq0}J_2(\lambda)\left(u^{(i)}\right) J_2(\lambda)\left(u^{(j)}\right)
\left(\frac{\partial^2 \Delta\mathcal{F}}{\partial u^{(i)}\partial u^{(j)}}
+\frac{\partial \Delta\mathcal{F}}{\partial u^{(i)}}\frac{\partial \Delta\mathcal{F}}{\partial u^{(j)}}\right)+ \frac{1}{2\epsilon^2} J_2(\lambda)^2(\mathcal{F}_0)\biggr]_{{\rm reg},-}.  \label{loop-initial}
\end{align}
\end{lem}
\begin{prf}
From Lemma~\ref{virasoro-total} we know that the Virasoro constraints
$L(\lambda)Z(x,\mathbf{s};\epsilon)=0$ can be represented as
\begin{align}
&\biggl[J_1(\lambda)J_2(\lambda)(\mathcal{F}_0)
+\frac{1}{2 \epsilon^2}J_2(\lambda)(\mathcal{F}_0)^2\biggr]_{{\rm reg},-}
+\frac{x^2}{2\lambda}=0,\label{symmetry-0}\\
&\biggl[J_1(\lambda)J_2(\lambda)(\Delta\mathcal{F}) + \frac1{\epsilon^2}J_2(\lambda)(\mathcal{F}_0)J_2(\lambda)(\Delta\mathcal{F})\nonumber\\
& \quad +\frac{1}{2}\left(J_2(\lambda)^2(\Delta\mathcal{F})
+J_2(\lambda)(\Delta\mathcal{F})^2\right)
+\frac{1}{2\epsilon^2}J_2(\lambda)^2(\mathcal{F}_0)\biggr]_{{\rm reg},-}
+\frac{\sigma_1}{24\lambda}=0\label{symmetry-1}.
\end{align}
It is easy to see that~\eqref{symmetry-1} can be rewritten in the form
\begin{equation}
D(\lambda)(\Delta\mathcal{F})
=\frac{\sigma_1}{24\lambda} + \left[\frac{1}{2} \left(J_2(\lambda)^2(\Delta\mathcal{F})
+J_2(\lambda)(\Delta\mathcal{F})^2\right)
+\frac{1}{2\epsilon^2}J_2(\lambda)^2(\mathcal{F}_0)\right]_{{\rm reg},-}.
\end{equation}
Since $J_2(\lambda)$ is a derivation,
we obtain~\eqref{loop-initial}.
\end{prf}

Let us proceed to simplify equation~\eqref{loop-initial}.
Introduce the following Puiseux series:
\begin{align} \label{eeplus}
E_{\alpha}=
\Bigg\{
\begin{array}{ll}
\sum_{n\geq0} c_{\alpha+hn} \zeta^{n+\frac{\alpha}{K_1}},
&\quad \alpha=0,\dots,K_1-1,  \\
&\\
\sum_{n\geq0} c_{\alpha+hn} \zeta^{n-\frac{\alpha}{K_2}},
& \quad  \alpha=-(K_2-1),\dots,-1,
\end{array}
\end{align}
where $\zeta:=e^{u(x,{\bf s})}/\lambda$.
Define
\begin{equation}\label{zh-2}
B_{i,j} := \sum_{\alpha,\beta\in I} \partial_x^i \left(E_{\alpha}\right) G^{\alpha\beta}\partial_x^j \left(E_{\beta}\right),
\end{equation}
where $G^{\alpha\beta}$ are the constants defined in~\eqref{Galphabeta}.
Note that  $B_{i,j}\in \mathbb{C}\bigl[u^{(1)},u^{(2)},\dots\bigr] [[\zeta]]$ for all~$i,j\geq 0$.

\begin{lem}\label{coef}
The following formulae hold true:
\begin{align}
& D(\lambda)\left(u^{(i)}\right)
=\frac{1}{\lambda}\left(\partial_x^iB_{0,0}+\sum_{j=1}^i\binom{i}{j}B_{j-1,i-j+1}\right), \quad \forall\,i\geq 0, \label{c111}\\
& \epsilon^{-2}\left[J_2(\lambda)\left(u^{(i)}\right)J_2(\lambda)\left(u^{(j)}\right)\right]_{{\rm reg},-}=\frac{1}{\lambda}B_{i+1,j+1}, \quad \forall\, i,j\geq 0. \label{c222}
\end{align}
\end{lem}
\begin{prf}
Using Lemma~\ref{prp-1} and the definition of~$J_2(\lambda)$ given in Lemma~\ref{virasoro-total}, we have
$$
\epsilon^{-1}J_2(\lambda)\left(\partial_x\mathcal{F}_0\right)+\frac{1}{\sqrt{\lambda}} =
\frac{1}{\sqrt{\lambda}}\left(E_0+\sqrt{\frac{K_2}{h}}\sum_{\alpha=1}^{K_1-1}E_{\alpha}
+\sqrt{\frac{K_1}{h}}\sum_{\alpha=-(K_2-1)}^{-1}E_{\alpha}\right).
$$
Acting~$\partial_x^2$ on the both sides of equation~\eqref{symmetry-0} yields
\[D(\lambda)(u)=\Bigl[\big(\epsilon^{-1}J_2(\lambda)(\partial_x\mathcal{F}_0)+\frac{1}{\sqrt{\lambda}}\big)^2\Bigr]_{{\rm reg},-}
=\frac{1}{\lambda} B_{0,0}, \]
which gives the validity of~\eqref{c111} for $i=0$. Assume that formula~\eqref{c111} is true for $i=k$;
then for $i=k+1$ we have
\begin{align*}
 &D(\lambda)\left(u^{(k+1)}\right)
= \partial_x D(\lambda)\left(u^{(k)}\right)+[D(\lambda),\,\partial_x]\left(u^{(k)}\right)\\
=& \frac{1}{\lambda}\left(\partial_x^{k+1}B_{0,0}+\sum_{j=1}^k\binom{k}{j}\left(B_{j,k-j+1}+B_{j-1,k-j+2}\right)\right)+ \epsilon^{-1}\Bigl[ J_2(\lambda)(\partial_x\mathcal{F}_0)
+\tfrac{1}{\sqrt{\lambda}} J_2(\lambda)\left(u^{(k)}\right)\Bigr]_{{\rm reg},-}\\
=& \frac{1}{\lambda}\left(\partial_x^{k+1}B_{0,0}+\sum_{j=1}^{k+1}\binom{k+1}{j}B_{j-1,k-j+2}\right).
\end{align*}
Hence by using mathematical induction we arrive at the formula~\eqref{c111}. Formula~\eqref{c222} can be verified directly. The lemma is proved.
\end{prf}

Let us now introduce a family of polynomials $f_{i,j} \in \mathbb{C}\bigl[u^{(1)},u^{(2)},\dots\bigr]$, $i,j\geq 0$ by
requiring the validity of the identity
\begin{equation}\label{deff}
\partial_x^i h(\zeta)=\sum_{j=0}^i f_{i,j}\left(\zeta\partial_\zeta\right)^j h(\zeta)
\end{equation}
for any smooth function $h(\zeta)$ of~$\zeta=e^{u(x,{\bf s})}/\lambda$. Clearly, $f_{i,j}$ vanishes if $i<j$.
For $i\geq j$, the polynomials $f_{i,j}$ can be uniquely determined by the following recursive relations:
\begin{align*}
&f_{i,0}=\delta_{i,0},\\
&f_{i+1,j+1}=\partial_x f_{i,j+1}+u^{(1)} f_{i,j}.
\end{align*}
Explicit expressions for~$f_{i,j}$ will be given in Section~\ref{section4}.
The functions $B_{i, j}$ defined in~\eqref{zh-2} can now be written as
\begin{equation}\label{f-0}
B_{i,j}=\sum_{k=0}^i\sum_{l=0}^jf_{i,k}f_{j,l}\widetilde{B}_{k,l},
\end{equation}
where
\begin{equation}
\widetilde{B}_{k,l}:= \sum_{\alpha,\beta\in I} \left(\zeta\partial_{\zeta}\right)^k \left(E_{\alpha}\right)
G^{\alpha\beta} \left(\zeta\partial_{\zeta}\right)^l \left(E_{\beta}\right). \label{widetb}
\end{equation}

\begin{lem}\label{algo}
The functions $\widetilde{B}_{i,j}$, $i,j\geq 0$ defined in~\eqref{widetb} satisfy the relations
\begin{align}
& \sum_{k\geq0}z^{-k}\widetilde{B}_{0,k} =\sqrt{z}\Phi(z)\frac{1}{1-K \zeta e^{-\partial_z}}\left(\frac{1}{\sqrt{z}\Phi(z)}\right), \label{f-1} \\
& \widetilde{B}_{i,j}=\widetilde{B}_{j,i},\quad
\zeta\partial_\zeta \widetilde{B}_{i,j}=\widetilde{B}_{i+1,j}+\widetilde{B}_{i,j+1}. \label{recur-f}
\end{align}
Moreover, $\widetilde{B}_{i,j}$ are polynomials of~$(1-K\zeta)^{-1}$ with degrees less than or equal to~$i+j+1$.
\end{lem}
\begin{prf}
By substituting~\eqref{eeplus} into~\eqref{widetb} we obtain
\begin{equation}
\widetilde{B}_{0,k}=\sum_{n\geq0}A_{k,n}\zeta^n,
\end{equation}
where
\begin{align*}
&A_{k,n}:=\sum_{\ell=1}^{n-1}b_{h\ell}^k c_{h\ell} c_{h(n-\ell)} +c_{h n} \delta_{k,0}\\
&+ \frac{K_2}{h}\sum_{\alpha=1}^{K_1-1}\sum_{\ell=0}^{n-1}b_{\alpha+h\ell}^k c_{\alpha+h\ell}c_{K_1-\alpha+h(n-1-\ell)}
+ \frac{K_1}{h}\sum_{\alpha=-(K_2-1)}^{-1}\sum_{\ell=0}^{n-1}b_{\alpha+h\ell}^k c_{\alpha+h\ell}c_{-\alpha-K_2+h(n-1-\ell)}.\\
\end{align*}
Using Lemma~\ref{const} we find that
\begin{align}
K^{-n} A_{k,n}
=&\sum_{\alpha\in I}\sum_{\ell=0}^{n-1}b_{\alpha+h\ell}^{k-1}\Res_{z=b_{\alpha+h\ell}}V_n(z)
 +\frac{\Res_{z=n}V_n(z)}{n}\delta_{k,0}.
\end{align}
In particular, we have
\[A_{0,n}=K^n\left(1-V_n(0)+\frac1n\Res_{z=n}V_n(z)\right)=K^n.\]
Hence
\begin{align}
& \sum_{k\geq0}z^{-k}\widetilde{B}_{0,k} =  \sum_{k,n\geq0}z^{-k}A_{k,n}\zeta^n  = \sum_{n\geq0}V_n(z) \left(K\zeta\right)^n \nonumber\\
& = \sum_{n\geq 0} \frac{\Phi(z)}{\Phi(z-n)}\sqrt{\frac{z}{z-n}} \left(K\zeta\right)^n = \sqrt{z}\Phi(z)\frac{1}{1-K\zeta e^{-\partial_z}}\left(\frac{1}{\sqrt{z}\Phi(z)}\right). \label{324}
\end{align}
Note that
\begin{equation}
\frac1{1-K \zeta e^{-\partial_z}} = \sum_{n\geq 0} \sum_{k=1}^{n+1} \frac{Q(n,k)}{(1-K\zeta)^k} \partial_z^n,\quad Q(n,k)=\frac{1}{k}\sum_{i=1}^k(-1)^{i-1}\binom{k}{i}i^{n+1}. \label{QQQQ}
\end{equation}
So from~\eqref{recur-f}, \eqref{324} it follows that $\widetilde{B}_{i,j}\in\mathbb{C}\bigl[(1-K\zeta)^{-1}\bigr]$.
The lemma is proved.
\end{prf}

\begin{rmk}
Introduce a trinomial curve
\begin{equation} \label{curve}
X Y^{K_1+K_2} - Y^{K_2} + 1 = 0 .
\end{equation}
Near $X=0$ we have two Puiseux series solutions:
\begin{align}
& Y_- = \sum_{m\geq 0} c_{-,m} X^m , \quad c_{-,0}=1, \nonumber \\
& Y_+ = X^{-\frac1{K_1}} \sum_{m\geq 0} c_{+,m} X^{m\frac{K_2}{K_1}},  \quad c_{+,0}=1. \nonumber
\end{align}
By using the Lagrange inversion we obtain the explicit formulae
\[
(Y_-)^{K_2} = \sum_{m\geq 0} \binom{m \frac{K_1+K_2}{K_2}}{m} \frac{X^m}{1+m\frac{K_1}{K_2}} ,
\quad (Y_+)^{-K_1} = \sum_{m\geq 0} \binom{m \frac{K_1+K_2}{K_1}}{m} \frac{X^{1+m \frac{K_2}{K_1}}}{1+m\frac{K_2}{K_1}}.
\]
We also observe the validity of the following relations
\begin{align}
& \sum_{\alpha=0}^{K_1-1} E_{\alpha} = \frac{d}{dX} \left((Y_+)^{-K_1}\right),
\qquad  \sum_{\alpha=-(K_2-1)}^{0} E_{\alpha} = \frac{K_1}{K_2} X \frac{d}{dX} \left((Y_-)^{K_2}\right) +(Y_-)^{K_2}
\end{align}
with $X=\zeta^{1/K_2}$.
\end{rmk}
\begin{cor}\label{J2-F0}
The following formula holds true:
\begin{equation}
\epsilon^{-2} \Bigl[J_2(\lambda)^2 \left(\mathcal{F}_0\right)\Bigr]_{{\rm reg},-}
=\frac{1}{\lambda}
\left(\frac{1}{8(1-K\zeta)^2}-\left(\frac{1}{8}-\frac{\sigma_1}{12}\right)\frac{1}{1-K \zeta}-\frac{\sigma_1}{12}\right).\label{zh-3}
\end{equation}
\end{cor}

\begin{prf} By using the definition of~$J_2$ given in Lemma~\ref{virasoro-total} and by using Lemma~\ref{prp-1} we obtain
\begin{equation}\label{zh-4}
\epsilon^{-2} \left[J_2(\lambda)^2(\mathcal{F}_0)\right]_{{\rm reg},-}=\int \frac{1}{\lambda}\frac{B_{1,1}}{u^{(1)}} dx =
\frac{1}{\lambda}\int \frac{\widetilde{B}_{1,1}}{K\zeta}d(K\zeta).
\end{equation}
Here we recall that $B_{1,1}$ is defined in~\eqref{zh-2}, and
the integration constant is chosen such that the left hand side goes to zero as $\lambda\rightarrow \infty$.
From Lemma~\ref{algo} it follows that
\begin{align*}
&\widetilde{B}_{1,1}
=\frac{1}{4}\frac{1}{(1-K\zeta)^3}
-\left(\frac{3}{8}-\frac{\sigma_1}{12}\right)\frac{1}{(1-K \zeta)^2}
+\left(\frac{1}{8}-\frac{\sigma_1}{12}\right)\frac{1}{1-K \zeta},
\end{align*}
which, together with~\eqref{zh-4}, leads to the formula~\eqref{zh-3}. The corollary is proved.\
\end{prf}

By using Lemmas~\ref{lemma35}--\ref{coef} and Corollary~\ref{J2-F0}
we arrive at the following theorem.
\begin{thm}\label{main}
The series $\Delta \mathcal{F}$ defined in \eqref{zh-5}, \eqref{zh-6} satisfies the Dubrovin-Zhang loop equation
\begin{align}
& \sum_{i\geq0}\frac{\partial \Delta \mathcal{F}}{\partial u^{(i)}} \partial_x^i \Theta
+\sum_{i\geq1} \sum_{j=1}^i \sum_{\alpha,\beta\in I} \frac{\partial \Delta \mathcal{F}}{\partial u^{(i)}}\binom{i}{j}
\partial_x^{j-1}\left(E_{\alpha}\right)G^{\alpha\beta}\partial_x^{i-j+1}\left(E_{\beta}\right)
 \nonumber \\
& =  \frac{\Theta^2}{16}-\left(\frac{1}{16}-\frac{\sigma_1}{24}\right)\Theta + \epsilon^2\sum_{i\geq0}
\frac{\partial \Delta \mathcal{F}}{\partial u^{(i)}} \partial_x^{i+2}\left(\frac{\Theta^2}{16}- \left(\frac{1}{16}-\frac{\sigma_1}{24}\right)\Theta\right) \nonumber \\
& \quad +\frac{\epsilon^2}{2}\sum_{i,j\geq0} \sum_{\alpha,\beta \in I}  \left(\frac{\partial^2\Delta \mathcal{F}}{\partial u^{(i)} \partial u^{(j)}}
+\frac{\partial \Delta \mathcal{F}}{\partial u^{(i)}}\frac{\partial \Delta \mathcal{F}}{\partial u^{(j)}}\right)
\partial_x^{i+1}\left(E_{\alpha}\right)G^{\alpha\beta}\partial_x^{j+1}\left(E_{\beta}\right), \label{loop-eq}
\end{align}
where $\Theta= \frac1{1- K \zeta}$, $K=h^h K_1^{-K_1} K_2^{-K_2}$, and
$\sigma_1= -(p+q+r)=\frac1{h}-\frac1{K_1}-\frac1{K_2}$.
\end{thm}
Note that each side of the loop equation~\eqref{loop-eq} is a power series of~$\Theta$.
It is understood that this equation for~$\Delta \mathcal{F}$ holds identically in~$\Theta$.
It is easy to check that $\mathcal{F}_g$ also satisfy the equations
\begin{equation}
\mathcal{F}_1= \frac1{24} \log u^{(1)} + \frac{\sigma_1}{24} u^{(0)}, \quad
\sum_{j=1}^{3g-2} j u^{(j)} \frac{\partial \mathcal{F}_g}{\partial u^{(j)}} = (2g-2) \mathcal{F}_g, \quad g\geq 2.  \label{eqfix}
\end{equation}

\begin{prp}\label{uniquenessprop}
The solution to equations \eqref{loop-eq}--\eqref{eqfix} is unique.
\end{prp}
\begin{prf}
For $g\geq 1$, by comparing the coefficients of~$\epsilon^{2g-2}$ in the both sides of~\eqref{loop-eq}, we obtain
\begin{align*}
&\sum_{i\geq 0}\left(\partial_x^i B_{0,0}
+\sum_{j=1}^i\binom{i}{j}B_{j-1,i-j+1}\right)
\frac{\partial \mathcal{F}_1}{\partial u^{(i)}}=\frac{\Theta^2}{16}-\left(\frac{1}{16}-\frac{\sigma_1}{24}\right) \Theta , \nonumber\\
&\sum_{i\geq 0}\left(\partial_x^i B_{0,0}
+\sum_{j=1}^i\binom{i}{j}B_{j-1,i-j+1}\right)\frac{\partial \mathcal{F}_g}{\partial u^{(i)}}\\
&=\sum_{i\geq0}\partial_x^{i+2}
\left(\frac{\Theta^2}{16}
- \left(\frac{1}{16}-\frac{\sigma_1}{24}\right)\Theta\right)
\frac{\partial\mathcal{F}_{g-1}}{\partial u^{(i)}},\\
& \quad + \frac12\sum_{i,j\geq0}B_{i+1,j+1}
\left(\frac{\partial^2 \mathcal{F}_{g-1}}{\partial u^{(i)}\partial u^{(j)}}
+ \sum_{k=1}^{g-1}\frac{\partial \mathcal{F}_k}{\partial u^{(i)}}\frac{\partial \mathcal{F}_{g-k}}{\partial u^{(j)}}\right), \quad g\geq 2. \nonumber
\end{align*}
By using the fact that
$\partial \mathcal{F}_g/\partial u^{(i)}=0$ for $i\geq 3g-1$,
and that $B_{i,j}$
are polynomials in~$\Theta$ of degrees~$i+j+1$, we arrive at the following system of equations:
\begin{align*}
&\left(\Theta,\cdots,\Theta^{3g-1}\right) \, M_g \,
\left(\frac{\partial \mathcal{F}_g}{\partial u^{(0)}},\cdots,\frac{\partial \mathcal{F}_g}{\partial u^{(3g-2)}}\right)=\left(\Theta,\cdots,\Theta^{3g-1}\right) \, N_g, \qquad g\geq 1,
\end{align*}
where $M_g$ is an invertible upper triangular $(3g-1)$ by $(3g-1)$ matrix, and $N_g$
is a column vector. All the elements of the matrix and the vector are differential polynomials in~$u$.
So the gradient of $\mathcal{F}_g$ can be uniquely determined.
The proposition is proved.
\end{prf}

\section{Loop equation: the general case}\label{section4}
In this section we drop the rational condition~\eqref{rational-condition} and consider the general case when $p, q, r$ are arbitrary complex numbers satisfying the {\it local Calabi-Yau condition \eqref{calabi-yau}}. We will give and prove
a refined version of the loop equation for the
corresponding Hodge free energies.

As it is pointed out in Lemma~\ref{jetvariablesrep}, there are polynomials
$H_g(z_0, z_1, \dots, z_{3g-2};p,q,r)$ in $p, q, r$ which satisfy the relations~\eqref{HHgHg}
and the relations~\eqref{eq2233} (see \cite{DLYZ-1}), here $g\geq 1$. It is then clear that
$H_g(z_0, \dots, z_{3g-2};p,q,r)$ are also polynomials of $\sigma_1,\sigma_3$ for~$g\geq 1$. For simplicity, we will denote them by
$H_g(z_0,\dots,z_{3g-2}; \sigma_1,\sigma_3)$.

We will need the following lemma.
\begin{lem}\label{elem}
Let $G(\sigma_1,\sigma_3)$ be a polynomial in $\mathbb{C}[\sigma_1,\sigma_3]$.
If $G$ vanishes for the values $(\sigma_1,\sigma_3)=\left(\frac1{K_1+K_2}-\frac1{K_1}-\frac1{K_2},
\frac2{(K_1+K_2)^3}- \frac{2}{K_1^3} - \frac{2}{K_2^3}\right)$, where
$K_1,K_2 \in \mathbb{N}$, $(K_1,K_2)=1$, then $G\equiv 0$.
\end{lem}
\begin{prf}
Suppose $G\not\equiv 0$.
Fix $K_1$ to be any positive integer. We observe that the points $\left(\frac1{K_1+K_2}-\frac1{K_1}-\frac1{K_2},
\frac2{(K_1+K_2)^3}- \frac2{K_1^3} - \frac{2}{K_2^3}\right)$, $K_2\in \mathbb{N}$ belong to the irreducible algebraic curve
\[2 K_1^3 \sigma_1^3 - 6 K_1 \sigma_1 - 6 - K_1^3 \sigma_3 =0\]
on the $(\sigma_1,\sigma_3)$ plane. It is easy to see that there are infinitely
many of such points with $(K_2,K_1)=1$. Hence the polynomial
\[2 K_1^3 \sigma_1^3 - 6 K_1 \sigma_1 - 6 - K_1^3 \sigma_3\]
must divide $G(\sigma_1,\sigma_3)$.
This contradicts with the fact that a polynomial must have a finite degree. The lemma is proved.
\end{prf}

It follows from Lemma~\ref{elem} that if a polynomial~$P$
in $\mathbb{C}[\sigma_1,\sigma_3][z_2,z_3,\dots]\bigl[z_1,z_1^{-1}\bigr]$
is equal to~$H_g(z_0,\dots,z_{3g-2};\sigma_1,\sigma_3)$
for
\[(\sigma_1,\sigma_3)=\left(\frac1{K_1+K_2}-\frac1{K_1}-\frac1{K_2},
\frac2{(K_1+K_2)^3}- \frac2{K_1^3} - \frac{2}{K_2^3}\right)\]
with
$K_1,K_2$ being arbitrary coprime positive integers, then $P\equiv H_g(z_0,\dots,z_{3g-2};\sigma_1,\sigma_3)$.
Here $g\geq 2$.
For $g=1$, we have the explicit expression
\[H_1(z_0,z_1;\sigma_1,\sigma_3)= \frac1{24} \log z_1 + \frac{\sigma_1}{24}z_0.\]

In order to give an efficient algorithm for the direct computation of $H_g(v,\dots,v_{3g-2};\sigma_1,\sigma_3)$, $g\geq 2$,
we proceed to derive the loop equation for these functions, which is a refinement of~\eqref{loop-eq}. To this end, we first introduce some notations.
Let $\mathcal{B}_{n,k}(X_1,\dots,X_{n-k+1})$ be the exponential Bell polynomials. They can be defined via the generating function
$$
\sum_{n,k\geq0}\mathcal{B}_{n,k} \left(X_1,\dots,X_{n-k+1}\right)\frac{y^n}{n!}z^k=\exp\left(z\sum_{j\geq1} X_j\frac{y^j}{j!}\right).
$$
Let $\mathcal{B}_n$ be the complete Bell polynomials, i.e.
\[
\mathcal{B}_0:=1,\quad \mathcal{B}_n(X_1,\dots,X_n) := \sum_{k=1}^n \mathcal{B}_{n,k} \left(X_1,\dots,X_{n-k+1}\right), \quad n\geq 1.
\]
Now we define $\widetilde{P}_{0,k}(\xi;\sigma_1,\sigma_3) \in \mathbb{C}[\sigma_1,\sigma_3] \bigl[(1-\xi)^{-1}\bigr]$, $k\geq0$ by the equation
\begin{equation}\label{zh-7}
\sum_{k\geq0}z^{-k}\widetilde{P}_{0,k}(\xi; \sigma_1,\sigma_3)
=\sqrt{z}\Phi(z;\sigma_1,\sigma_3)\frac{1}{1-\xi e^{-\partial_z}}\left(\frac{1}{\sqrt{z}\Phi(z;\sigma_1,\sigma_3)}\right),
\end{equation}
where
\[\Phi : =  \Phi(z;\sigma_1,\sigma_3)=\exp\left(-\sum_{i\geq1}\frac{B_{2i}}{2i(2i-1)}(p^{2i-1}+q^{2i-1}+r^{2i-1})z^{1-2i}\right).\]
Equation~\eqref{zh-7} can be written more explicitly as follows:
\begin{align}
&\sum_{n\geq0}z^{-n}\widetilde{P}_{0,n}(\xi;\sigma_1,\sigma_3)\nonumber\\
&=\sum_{n\geq0}\sum_{k=1}^{n+1}\frac{Q(n,k)}{(1-\xi)^k}
\sum_{m=0}^n \frac{(-1)^{n-m}(2n-2m-1)!!}{2^{n-m}m! (n-m)!}z^{-n+m}
\Phi\, \partial_z^m\left(\frac{1}{\Phi}\right), \label{B0xi}
\end{align}
where $Q(n,k)$ are the numbers defined in~\eqref{QQQQ}, and
\begin{align}
&\Phi \, \partial_z^m\left(\frac{1}{\Phi}\right)
=\mathcal{B}_m \left(-\partial_z\log\Phi(z),\dots,-\partial_z^m\log\Phi(z)\right), \quad m\geq0.\label{phiphi}
\end{align}
We define, for $i,j\geq0$, $\widetilde{P}_{i,j}=\widetilde{P}_{i,j}(\xi;\sigma_1,\sigma_3)$ by the following recursion:
\begin{align*}
\widetilde{P}_{i,j} &= \widetilde{P}_{j,i}, \\
\xi\partial_{\xi}\widetilde{P}_{i,j} &= \widetilde{P}_{i+1,j}+\widetilde{P}_{i,j+1}.
\end{align*}

It is easy to see that the functions $f_{i,k}$ defined in \eqref{deff} are just the Bell polynomials
$\mathcal{B}_{i,k}(z_1,\dots,z_{i-k+1})$. By using these functions we
introduce, for $i,j\geq 0$, the following notations:
$$
P_{i,j}=\sum_{k=0}^i\sum_{l=0}^j f_{i,k}(z_1,\dots,z_{i-k+1}) f_{j,l}(z_1,\dots,z_{j-l+1})\widetilde{P}_{k,l}(\xi;\sigma_1,\sigma_3).
$$

\medskip

\begin{prfn}{Theorem~\ref{refineloop}}
It is easy to verify that the functions $P_{i,j}$ defined above
coincide with the functions $B_{i,j}$ that are given in~\eqref{f-0} when $\xi =K \zeta$ and
\begin{align*}
& (\sigma_1, \sigma_3) = \left(\frac1{h}-\frac1{K_1}-\frac1{K_2},
\frac2{h^3}- \frac{2}{K_1^3} - \frac{2}{K_2^3}\right),\\
& (z_1,z_2,\dots) = \left(u^{(1)},u^{(2)},\dots\right),
\end{align*}
where $K_1,K_2$ are arbitrary coprime positive integers.
It then follows from Theorem~\ref{main} that
$\Delta H:=\sum_{g\geq 1} \epsilon^{2g-2} H_g(z_0,\dots,z_{3g-2};\sigma_1,\sigma_3)$ satisfies
equation~\eqref{loop-eq-0} for the particular values of $(\sigma_1,\sigma_3)$.
By definition we also see that $P_{i,j}\in \mathbb{C}\left[\sigma_1,\sigma_3][(1-\xi)^{-1};z_1,z_2,\dots\right]$.
The existence of solution of the loop equation~\eqref{refineloop} is then further ensured by Lemma~\ref{elem}.
As we do in the proof of Proposition~\ref{uniquenessprop}, we deduce from~\eqref{refineloop}
the following equation for each $g\geq 1$:
\begin{align*}
&\left(\Theta,\cdots,\Theta^{3g-1}\right) \, \overbar{M}_g(\sigma_1,\sigma_3) \,
\left(\frac{\partial H_g}{\partial z_0},\cdots,\frac{\partial H_g}{\partial z_{3g-2}}\right)
=\left(\Theta,\cdots,\Theta^{3g-1}\right) \, \overbar{N}_g(\sigma_1,\sigma_3),
\end{align*}
where $\overbar{M}_g(\sigma_1,\sigma_3)$ is an invertible $(3g-1)$ by $(3g-1)$ matrix, $\overbar{N}_g(\sigma_1,\sigma_3)$ is
a column vector, and their entries are polynomials
of $\sigma_1,\sigma_3$. Therefore we have verified that the gradient of~$H_g$ is uniquely determined by~\eqref{refineloop}.
The theorem is proved.
\end{prfn}

\begin{emp}
By taking the coefficient of $\epsilon^0$ in \eqref{loop-eq-0}, we obtain the following equation:
$$
P_{0,0}\frac{\partial H_1}{\partial z_0}
+\left(\partial P_{0,0}+P_{0,1}\right)\frac{\partial H_1}{\partial z_1}
=\frac{\Theta^2}{16}
-\left(\frac{1}{16}-\frac{\sigma_1}{24}\right) \Theta
$$
for~$H_1$.
The coefficients read
\[P_{0,0}=\Theta,\quad P_{0,1}
=\frac{z_1}{2}\left(\Theta^2-\Theta\right).\]
So we have
\begin{align*}
&\frac{3 z_1}{2}\frac{\partial H_1}{\partial z_1} \Theta^2
+\left(\frac{\partial H_1}{\partial z_0}-\frac{3 z_1}{2}\frac{\partial H_1}{\partial z_1}\right)
\Theta=\frac{\Theta^2}{16}
-\left(\frac{1}{16}-\frac{\sigma_1}{24}\right)\Theta.
\end{align*}
which gives
\[\frac{\partial H_1}{\partial z_0}=\frac{\sigma_1}{24},\quad
\frac{\partial H_1}{\partial z_1}=\frac{1}{24 z_1}.\]
Hence
\[H_1=\frac{1}{24}\log z_1+\frac{\sigma_1}{24} z_0\]
which gives~\eqref{explicitH1}.
In a similar way, we obtain formula~\eqref{explicitH2} and the following expression of~$H_3$:
\begin{align}
H_3
=&\frac{1}{82944}\frac{z_7}{z_1^3}-\frac{7}{46080}\frac{z_2 z_6}{z_1^4}
 -\frac{53}{161280}\frac{z_3 z_5}{z_1^4}
 +\frac{353}{322560}\frac{z_2^2 z_5}{z_1^5}
 -\frac{103}{483840}\frac{z_4^2}{z_1^4}\nonumber\\
 &+\frac{1273}{322560 }\frac{ z_2 z_3 z_4}{z_1^5}
-\frac{83 }{15120 }\frac{z_2^3 z_4}{z_1^6}+\frac{59 }{64512 }\frac{z_3^3}{z_1^5}
 -\frac{83  }{7168 }\frac{z_2^2 z_3^2}{z_1^6}+\frac{59}{3024 }\frac{ z_2^4 z_3}{z_1^7}\nonumber\\
&-\frac{5 }{648 }\frac{z_2^6}{z_1^8} +\frac{7 \sigma_1 }{138240}\frac{z_6}{ z_1^2}
 -\frac{383 \sigma_1 }{967680 }\frac{z_2 z_5}{z_1^3}
+\frac{41 \sigma_1^2 }{580608 }\frac{z_5}{z_1}-\frac{689 \sigma_1 }{967680}\frac{z_3 z_4}{ z_1^3}
+\frac{185 \sigma_1 }{96768 }\frac{z_2^2 z_4}{z_1^4}\nonumber\\
&-\frac{373 \sigma_1^2 }{1451520 }\frac{z_2 z_4}{z_1^2}
 +\left(\frac{23 \sigma_1^3}{580608}-\frac{11 \sigma_3}{2903040}\right)z_4
+\frac{869 \sigma_1 }{322560 }\frac{z_2 z_3^2}{z_1^4}
-\frac{61 \sigma_1^2 }{322560 }\frac{z_3^2}{z_1^2}\nonumber\\
& -\frac{9343 }{1451520 }\frac{z_2^3 z_3}{z_1^5}+\frac{151 \sigma_1^2 }{207360 }\frac{z_2^2 z_3}{z_1^3}
 -\left(\frac{19 \sigma_1^3}{1451520}-\frac{19 \sigma_3}{2903040}\right)\frac{z_2 z_3}{z_1} +\frac{131 \sigma_1 }{45360 }\frac{z_2^5}{z_1^6}\nonumber\\
&-\frac{19 \sigma_1^2 }{53760 }\frac{z_2^4}{z_1^4} +\left(\frac{41 \sigma_1^4}{4354560}-\frac{41 \sigma_1 \sigma_3}{8709120}\right)z_1 z_3
 +\left(\frac{\sigma_1^3}{108864}-\frac{\sigma_3}{217728}\right)\frac{z_2^3}{z_1^2} \nonumber \\
& +\left(\frac{31 \sigma_1^4}{4354560}-\frac{31 \sigma_1 \sigma_3}{8709120}\right)z_2^2
 -\left(\frac{\sigma_1^6}{13063680}-\frac{\sigma_1^3 \sigma_3}{13063680}
 +\frac{\sigma_3^2}{52254720}\right)z_1^4 . \label{explicitH3}
 \end{align}
\end{emp}

\smallskip

\begin{prfn}{Corollary~\ref{gapcor}}
The properties \eqref{gap1}--\eqref{gap2} follow from \eqref{eq2233}--\eqref{HgHg} and the simple fact that
 $v({\bf t}(x,{\bf s}))\in \log x + \mathbb{C}[x][[{\bf s}]]$, where it should be noted that due to Lemma~\ref{elem}
 the polynomials $R_g(\sigma_1,\sigma_3)$ must be unique. The degree estimate $\deg R_g(\sigma_1,\sigma_3)\leq 3g-3$ with~$g\geq 2$
 follows from the theorem~1.3 of~\cite{DLYZ-1}.
To show~\eqref{Faber} we first note that for $g\geq 2$,
\begin{align}
& \H_g({\bf t};p,q,r) \nonumber\\
& \quad = (-1)^{g-1} \sum_{n\geq 0} \sum_{i_1,\dots,i_n\geq 0} \frac{t_{i_1}\cdots t_{i_n}}{n!}
 \int_{\overbar{\mathcal{M}}_{g,n}}  \psi_1^{i_1}\cdots\psi_n^{i_n} \lambda_g \lambda_{g-1} \lambda_{g-2}  \left( p^g q^{g-1}r^{g-2} + (p\leftrightarrow q \leftrightarrow r) \right) \nonumber\\
& \quad\quad  + (-1)^{g-1} \sum_{n\geq 0} \sum_{i_1,\dots,i_n\geq 0} \frac{t_{i_1}\cdots t_{i_n}}{n!}
 \int_{\overbar{\mathcal{M}}_{g,n}}  \psi_1^{i_1}\cdots\psi_n^{i_n} \lambda_{g-1}^3 \, p^{g-1}q^{g-1}r^{g-1} + ~ \cdots ,
\end{align}
where
\[(p\leftrightarrow q \leftrightarrow r):=p^g r^{g-1} q^{g-2}+q^g p^{g-1}r^{g-2}+q^g r^{g-1}p^{g-2}+r^g p^{g-1}q^{g-2}+r^g q^{g-1}p^{g-2},\]
and the dots ``$\cdots$" denotes the terms in $\H_g({\bf t};p,q,r)$, as polynomials of $p,q,r$, having degree strictly less than $3g-3$.
Recalling~\eqref{sigma12}, \eqref{polynomiality} and \eqref{zh-10}, it is clear that the degree defined for polynomials in~$p,q,r$ and the one for polynomials in $\sigma_1$, $\sigma_3$ are consistent.
We further note that the following formula holds true for $g\geq 2$:
\[
\sum_{n\geq 0} \sum_{i_1,\dots,i_n\geq 0} \frac{t_{i_1}\cdots t_{i_n}}{n!}
 \int_{\overbar{\mathcal{M}}_{g,n}}  \psi_1^{i_1}\cdots\psi_n^{i_n} \lambda_g \lambda_{g-1} \lambda_{g-2} =
 \frac{1}{2(2g-2)!} \frac{|B_{2g}|}{2g} \frac{|B_{2g-2}|}{2g-2} \frac{\partial^{2g-2} v({\bf t})}{\partial t_0^{2g-2}},
\]
as it was proved in~\cite{DLYZ-1}.
The corollary is then proved by using Mumford's relation
\[\lambda_{g-1}^3= 2\lambda_g\lambda_{g-1}\lambda_{g-2}\] (see~\cite{M}),
and again the fact that $v({\bf t}(x,{\bf s}))\in \log x + \mathbb{C}[x][[{\bf s}]]$.
 \end{prfn}


\section{Introduction}

The notion of Heisenberg algebra and of Heisenberg group appear since long time in the description of  quantum dynamical systems and sophisticated techniques of algebra, topology and functional analysis are often used, in order to understand the behaviour of elementary particles or, more in general, of some quantum systems. When it is possible to detect algebraic structures like nilpotent Lie algebras, we are in a very good position to recognize symmetries and conservation laws of significant interest in physics. This is exactly the content of the well known Noether's Theorem, \cite{gold}. The reference of Snobl and others  \cite{snobl} show some recent techniques in the context of Lie algebras and theoretical physics, but here we will focus on a specific classification, due to the so called ``Berkovich--Zhou invariant'' (or ``corank'') of nilpotent finite dimensional Lie algebras (over a field of characteristic $0$).

We will refer to the recent works \cite{nr1, nr2, nr3}, which generalize \cite{bat, es, es2}, but, at the same time, the present contribution continues the efforts we made in \cite{bagrus2018}, where we looked at dynamical systems involving pseudo-bosons in terms of finite dimensional nilpotent Lie algebras for values of the dimension smaller than five.

 During our discussions, we will list a first classification
 of finite dimensional nilpotent Lie algebras by the size of
 their dimension, and, specifically, when the dimension is at
 most $5$ (see \cite{degraaf, morozov, turk} and Theorem
 \ref{classification} below). Then we will use  the notion of
 \textit{Schur multiplier} of a  Lie algebra, referring to \cite{weibel} for its classical approach via homological algebra.  Schur \cite{schur} introduced this concept in 1907 in group theory, but his approach was subject to generalizations of various nature for different algebraic and geometric structures.

 The Schur multiplier does not seem to be related with the classification of finite dimensional nilpotent Lie algebras by dimension, but this is only apparently true. Indeed we
 may define a non-negative integer (which turns out to be a
 numerical invariant), called ``corank''  and due to
 Berkovich and Zhou \cite{berkovich, zhou}. This number
 measures the size of the Schur multiplier of a nilpotent Lie
 algebra and we will see that the corank allows us to get a
 second classification of nilpotent Lie algebras, but this
 time the classification gives information on the way in
 which our Lie algebra may be written as the sum of two
 smaller abelian Lie algebras.

 This second classification is more useful if we
 want to focus on the structure of a general Lie algebra.
 Moreover,  this new classification presents again nilpotent
 Lie algebras of dimension $\le 5$ when the value of the
 corank is smaller than $6$ (see Theorem \ref{tl} below
 for more details).

 While the parallel approach of classification via corank and
 via dimension was known from a pure algebraic--geometric perspective (see
 \cite{bat, es, es2, mon, nr1, nr2}), one of the
 interesting result of the present paper is that explicit
 physical examples of dynamical systems can be related to
 this classification. This point is not known in literature,
 to the best of our knowledge.

  In fact we will show that all finite dimensional nilpotent
 Lie algebras, over the field of complex numbers may be
 realized  in terms of (pseudo-) bosonic operators  when
 their corank is at most $6$ (see Theorem \ref{main2}
 below). Moreover we will provide examples of significant
 interest in physics when the value of the corank is
 arbitrary big (see Theorem \ref{main1}).

Section 2 recalls some well known results and definitions of algebraic nature. They can be found in \cite{beltita,snobl,weibel} and are integrated by some recent material in \cite{nr1, nr2, nr3}.
 We refer to \cite[Section 3]{bagrus2018} for the terminology and notations which will be involved in the proofs of our main results, where a framework of functional analysis and of theory of the operators is applied. Finally, Section 3 illustrates  the main results with corresponding proofs.

\section{Some results of low dimensional cohomology}

In order to focus on the main topics, without dispersion of energies in additional details, we do not report the classical notion of \textit{Chevalley-Eilenberg complex} for a Lie algebra, since it can be found in \cite[Definition 7.7.1]{weibel}. Analogously, we will assume that the reader is familiar with the second cohomology Lie algebra with complex coefficients. This Lie algebra is widely described in \cite[Corollary 7.7.3 and Exercise 7.7.5]{weibel} and  we will call it  \textit{ Schur multiplier} of the Lie algebra $\mathfrak{l}$, using the notation $M(\mathfrak{l})$. In fact we will work only with finite dimensional  Lie algebras over the field of complex numbers.

The size of the Schur multiplier describes how many extensions we can actually form, beginning from an abelian Lie algebra with another one acting on it. Furthemore, we will see  that it is possible to introduce a numerical invariant, which we will call  ``corank'', or  ``Berkovich--Zhou invariant'', and this measures the size of the Schur multiplier from some ideal cases for which we known essentially all. The main results of this paper will indeed show that it is possible to find concrete examples for any possible value of this corank and, to the best of our knowledge, it is the first time that these examples can be found in physical contexts.

 We begin with the following notion, available in \cite{weibel}.

\begin{defn}\label{extensions}
A Lie algebra $\mathfrak{l}$ is an \textit{extension} of an abelian Lie algebra $\mathfrak{a}$ (called \textit{abelian kernel}) by another  Lie algebra $\mathfrak{b}$, if there is a short exact sequence
$$0 \ \longrightarrow \  \mathfrak{a} {\overset{\mu}{\longrightarrow}} \ \mathfrak{l} \ {\overset{\varepsilon}{\longrightarrow}} \mathfrak{b} \ \longrightarrow \ 0$$
such that  $\mathfrak{a}$ is an abelian ideal of $\mathfrak{l}$ and $\mathfrak{l}/\mathfrak{a} \simeq \mathfrak{b}$ is the Lie algebra quotient. Here $\mu$ is a monomorphism and $\varepsilon$ an epimorphism of Lie algebras such that $\ker \varepsilon = \mathrm{Im} \ \mu $. If in addition, $\mathfrak{a}$ is contained in the center \[Z(\mathfrak{l})=\{ a \in \mathfrak{l}  \ | \ [a,b]=0  \ \forall b \in \mathfrak{l}\} \]
 of $\mathfrak{l}$, the Lie algebra $\mathfrak{l}$ is called \textit{central extension} of $\mathfrak{a}$ by $\mathfrak{b}$.
\end{defn}

Definition \ref{extensions} is reported from \cite[Definition 7.6.1]{weibel}, where it is noted that $\mathfrak{a}$ may be endowed with the structure of $\mathfrak{b}$-module  defining $ba$ to be the product $[\widetilde{b},a]$ in $\mathfrak{l}$, where $\varepsilon(\widetilde{b})=b$ for all $a \in \mathfrak{a}$ and $b \in \mathfrak{b}$. Following \cite[Extension Problem 7.6.2]{weibel}, two extensions  $0 \ \longrightarrow \  \mathfrak{a} {\overset{\mu_1}{\longrightarrow}} \ \mathfrak{l}_1 \ {\overset{\varepsilon_1}{\longrightarrow}} \mathfrak{b} \ \longrightarrow \ 0$ and $0 \ \longrightarrow \  \mathfrak{a} {\overset{\mu_2}{\longrightarrow}} \ \mathfrak{l}_2 \ {\overset{\varepsilon_2}{\longrightarrow}} \mathfrak{b} \ \longrightarrow \ 0$ are equivalent if there is an isomorphism $\varphi$ such that the following diagram
$$\begin{CD}
0  @>>>   \mathfrak{a} @>\mu_1>> \mathfrak{l}_1 @>\varepsilon_1>>  \mathfrak{b} @>>> 0\\
 @. @| @VV\varphi V @| @. \\
0  @>>>   \mathfrak{a} @>\mu_2>> \mathfrak{l}_2 @>\varepsilon_2>>  \mathfrak{b} @>>> 0
\end{CD}
$$
is commutative. Then we can define the set $\mathrm{Ext}(\mathfrak{b},\mathfrak{a})$ of equivalence classes of the extensions of  $\mathfrak{a} $ by  $\mathfrak{b}$ realizing $\mathfrak{l}$. The problem of the description of $\mathrm{Ext}(\mathfrak{b},\mathfrak{a})$ is known as Extension Problem and a  classical result \cite[Classification Theorem 7.6.3]{ weibel} shows an important relation with the Schur multiplier.

\begin{lem}[See \cite{weibel}, Chapter 7]\label{h2}Given an extension $0 \ \rightarrow \  \mathfrak{a} \rightarrow \ \mathfrak{l} \rightarrow \mathfrak{b} \ \rightarrow \ 0$ of a Lie algebra $\mathfrak{l}$ as in Definition \ref{extensions}:
\begin{itemize}
\item[(i)]There is a 1-1 correspondence between $M(\mathfrak{l})$ and $\mathrm{Ext}(\mathfrak{b},\mathfrak{a})$;
\item[(ii)] If $M(\mathfrak{l})=0$, then $\mathfrak{l}$ cannot be realized as  extension with abelian kernel.
\end{itemize}
\end{lem}

As indicated in Lemma \ref{h2} (ii),  the size of the Schur multiplier actually measures how many extensions are possible. Further important examples of extensions are constructed below.

\begin{defn}\label{semidirect}
A Lie algebra $\mathfrak{l}$ is the {\em semidirect sum} of two of its Lie subalgebras $\mathfrak{a}$ and $\mathfrak{b}$ if  the following conditions are satisfied:
\begin{itemize}
\item[(i)]$\mathfrak{a}$  is an ideal of $\mathfrak{l}$,
\item[(ii)] $\mathfrak{l}= \mathfrak{a} + \mathfrak{b}$,
\item[(iii)] $\mathfrak{a} \cap \mathfrak{b} = 0$.
\end{itemize}
Of course, if $\mathfrak{a} \subseteq Z(\mathfrak{l})$, then $\mathfrak{l}$ is a central extension of $\mathfrak{a}$ by $\mathfrak{b}$.
\end{defn}

Having in mind Lemma \ref{h2} (i),  Definition \ref{semidirect} and Definition \ref{extensions}, it is clear that semidirect sums are special cases of extensions with abelian kernel, and so the Schur multiplier gives information of structure of a Lie algebra. In this perspective, we recall another classical result. A semisimple Lie algebra is a sum of simple Lie algebras. A well known theorem of Levi illustrates the importance of the semidirect sum in the structure of the finite dimensional Lie algebras.

\begin{thm}[See \cite{weibel}, Chapter 7, Levi Decomposition]
Any finite dimensional Lie algebra $\mathfrak{l}$ over $\mathbb{C}$  is semidirect sum of a solvable Lie algebra $\mathfrak{a}$ by a semisimple Lie algebra $\mathfrak{b}$.
\end{thm}

 Note that a solvable Lie algebra is a Lie algebra which
 contains a finite series of ideals in which each factor is
 an abelian Lie algebra. While the role of the Schur
 multiplier of solvable Lie algebras is quite difficult to
 compute, the  result above is classical and may help for the
 case of semisimple Lie algebras. An important consequence is
 the following.

\begin{cor}[See \cite{weibel}, Chapter 7]\label{semisimple} A semisimple Lie algebra $\mathfrak{l}$ over $\mathbb{C}$ of finite dimension has $M(\mathfrak{l})=0$.
\end{cor}

For a series of motivations, which will be clear later on, we need to investigate Schur multipliers of Lie algebras in which we may iterate the notion of center.

\begin{defn}\label{center}If $\mathfrak{l}$ is a Lie algebra, the \textit{derived} Lie subalgebra of $\mathfrak{l}$ is defined as the smallest Lie ideal containing all the commutators  $[a,b]=ab-ba$  of $\mathfrak{l}$, that is,
\[[\mathfrak{l},\mathfrak{l}]=\langle [a,b]  \ | \ a,b \in \mathfrak{l} \rangle. \]
Note that the Lie algebra quotient $\mathfrak{l}/[\mathfrak{l},\mathfrak{l}]$ is always abelian,
 \end{defn}

Both  $Z(\mathfrak{l})$ and $[\mathfrak{l},\mathfrak{l}]$ are ideals of $\mathfrak{l}$ and so the size of their quotients may  measure  how far we are from having an abelian Lie algebra. Since we do not know whether $\mathfrak{l}/Z(\mathfrak{l})$ is abelian or not in general,  there is a way to generalize the above notions via appropriate series (see \cite{ weibel}).

\begin{defn}\label{uppercentralseries}The \textit{upper central series} of $\mathfrak{l}$ is the series
\[0=Z_0(\mathfrak{l}) \leq Z_1(\mathfrak{l})=Z(\mathfrak{l}) \leq Z_2(\mathfrak{l}) \leq \ldots \leq Z_i(\mathfrak{l}) \le Z_{i+1}(\mathfrak{l}) \leq \ldots, \]
where each $Z_i(\mathfrak{l})$ turns out to be an ideal of $\mathfrak{l}$ (called $i$th $center$ of $\mathfrak{l}$) and
\[ \frac{Z_1(\mathfrak{l})}{Z_0(\mathfrak{l})}=Z(\mathfrak{l}), \frac{Z_2(\mathfrak{l})}{Z_1(\mathfrak{l})}=Z\left(\frac{\mathfrak{l}}{Z(\mathfrak{l})}\right),   \ldots, \frac{Z_{i+1}(\mathfrak{l})}{Z_i(\mathfrak{l})}=Z\left(\frac{\mathfrak{l}}{Z_i(\mathfrak{l})}\right), \ldots \]
is an abelian quotient Lie algebra. We say that $\mathfrak{l}$ is  \textit{nilpotent of class $c$} if $\mathfrak{l}=Z_c(\mathfrak{l})$, that is, if the upper central series of $\mathfrak{l}$ ends after a finite number $c \ge 0$ of steps.

By analogy, the \textit{lower central series} of $\mathfrak{l}$ is defined  inductively by
\[\gamma_1(\mathfrak{l})=\mathfrak{l} \geq \gamma_2(\mathfrak{l})=[\mathfrak{l},\mathfrak{l}]  \geq  \gamma_3(\mathfrak{l})=[[\mathfrak{l},\mathfrak{l}],\mathfrak{l}]  \geq \ldots \geq \gamma_i(\mathfrak{l}) \geq \gamma_{i+1}(\mathfrak{l}) \geq \ldots, \]
where each $\gamma_i(\mathfrak{l})$ turns out to be an ideal of $\mathfrak{l}$ (called $i$th $derived$ of $\mathfrak{l}$).
\end{defn}

The presence of finite dimensional nilpotent Lie algebras plays a fundamental role in mathematical physics, because these structures have been observed in several quantum dynamical systems (see \cite{beltita}).  Probably one of the first nilpotent Lie algebras, which was used for the description of the behaviour of complex dynamical systems, in which bosons are involved, is due to Heisenberg and we report the formal definition here.

\begin{defn}\label{Heisenberg}  A finite dimensional Lie algebra  $\mathfrak{l}$ is called $Heisenberg$ provided that $[\mathfrak{l},\mathfrak{l}]=Z(\mathfrak{l})$ and $\mathrm{dim}([\mathfrak{l},\mathfrak{l}]) = 1$. Such algebras are odd dimensional with
basis $v_1, \ldots , v_{2m}, v$ and the only non--zero commutator between basis elements is $[v_{2i-1}, v_{2i}] = -
[v_{2i}, v_{2i-1}]= v$ for $i = 1,2, \ldots ,m$. The symbol $\mathfrak{h}(m)$
denotes the Heisenberg algebra of dimension $2m + 1$.
\end{defn}

 One can see easily that $\mathfrak{h}(m)$ is a nilpotent
 Lie algebra of class $2$, that is, such that
 $Z_2(\mathfrak{h}(m)))=\mathfrak{h}(m)$, of
 arbitrary dimension of the form $2m+1$,  and such that the
 Lie algebra quotient
 $\mathfrak{h}(m)/Z(\mathfrak{h}(m))$ is abelian,
 where $Z(\mathfrak{h}(m))$ a 1-dimensional Lie algebra.
 It turns out that the structure of finite dimensional
 nilpotent Lie algebras of low dimension involves strongly
 Heisenberg algebras and may be described in terms of
 generators and relations in the following way.

\begin{thm}[Classification of Finite Dimensional Nilpotent Lie Algebras of Dimension 3, 4 and  5]\label{classification}
Let $\mathfrak{l}$ be a finite dimensional nilpotent Lie algebra (over any field of characteristic $\neq 2$) and $\mathfrak{i}$ an abelian Lie algebra of dimension 1. Then
\begin{itemize}
\item[(1)] $\mathrm{dim} \ \mathfrak{l} =3$ if and only if $\mathfrak{l}$ is isomorphic to one of the following Lie algebras:
 \begin{itemize}
\item[-] $\mathfrak{l}_{3,1}$ abelian Lie algebra of dimension 3,
\item[-] $\mathfrak{l}_{3,2}  \simeq \mathfrak{h}(1)$.
\end{itemize}
\item[(2)]$\mathrm{dim} \ \mathfrak{l} =4$ if and only if    $\mathfrak{l}$ is isomorphic to one of the following Lie algebras:
 \begin{itemize}
\item[-] $\mathfrak{l}_{4,1} =  \mathfrak{l}_{3,1} \oplus \mathfrak{i}$,
\item[-] $\mathfrak{l}_{4,2} =  \mathfrak{l}_{3,2} \oplus \mathfrak{i}$
\item[-] $\mathfrak{l}_{4,3} = \langle v_1, v_2, v_3, v_4 \ | \ [v_1,v_2]=v_3, [v_1,v_3]=v_4 \rangle$
\end{itemize}
\item[(3)]$\mathrm{dim} \ \mathfrak{l} =5$ if and only if $\mathfrak{l}$ is isomorphic to one of the following Lie algebras:
\begin{itemize}
\item[-] $\mathfrak{l}_{5,1} = \mathfrak{l}_{4,1} \oplus \mathfrak{i} $,
\item[-] $\mathfrak{l}_{5,2} = \mathfrak{l}_{4,2} \oplus \mathfrak{i} \simeq \mathfrak{h}(1) \oplus \mathfrak{i} \oplus \mathfrak{i}$,
\item[-] $\mathfrak{l}_{5,3} = \mathfrak{l}_{4,3} \oplus \mathfrak{i} $,
\item[-] $\mathfrak{l}_{5,4} =  \langle v_1, v_2, v_3, v_4, v_5 \ | \ [v_1,v_2]=[v_3,v_4]=v_5 \rangle \simeq \mathfrak{h}(2) $,
\item[-] $\mathfrak{l}_{5,5} =  \langle v_1, v_2, v_3, v_4, v_5 \ | \ [v_1,v_2]=v_3, [v_1,v_3]=[v_2,v_4]=v_5 \rangle$,
\item[-] $\mathfrak{l}_{5,6} =  \langle v_1, v_2, v_3, v_4, v_5 \ | \ [v_1,v_2]=v_3, [v_1,v_3]=v_4,$ $$[v_1,v_4]=[v_2,,v_3]=v_5  \rangle,$$
\item[-] $\mathfrak{l}_{5,7} =  \langle v_1, v_2, v_3, v_4, v_5 \ | \ [v_1,v_2]=v_3, [v_1,v_3]=v_4, [v_1,v_4]=v_5 \rangle$,
\item[-] $\mathfrak{l}_{5,8} =  \langle v_1, v_2, v_3, v_4, v_5 \ | \ [v_1,v_2]=v_4, [v_1,v_3]=v_5 \rangle$,
\item[-] $\mathfrak{l}_{5,9} =  \langle v_1, v_2, v_3, v_4, v_5 \ | \ [v_1,v_2]=v_3, [v_1,v_3]=v_4, [v_2,v_3]=v_5 \rangle$.
\end{itemize}
\end{itemize}
\end{thm}

An important observation can be done here: this is designed for our purposes.

\begin{lem}\label{technical}
The nilpotent Lie algebras $\mathfrak{l}_{4,3}$, $\mathfrak{l}_{5,4}$, $\mathfrak{l}_{5,5}$,  $\mathfrak{l}_{5,6}$, $\mathfrak{l}_{5,7}$, $\mathfrak{l}_{5,8}$ and $\mathfrak{l}_{5,9}$ may be realized as central extensions.
\end{lem}

\begin{proof} We begin to consider $\mathfrak{l}_{4,3}$, which is  described at \cite[Page 3530]{har1}. We note that  $\mathfrak{l}_{4,3}$ is nilpotent of class $2$ and is realized as central extension of a  1-dimensional ideal $\mathfrak{i}=\mathfrak{a} \subseteq Z(\mathfrak{l}_{4,3})$ by the  quotient  $\mathfrak{b}=\mathfrak{l}_{4,3}/\mathfrak{a}=\mathfrak{h}(1)$.  Definition \ref{extensions} is satisfied. On the other hand, it is clear that the Heisenberg algebra of dimension 5 is realized as central extension, just looking at the way in which it is defined. Then we look at \cite[Pages 3531--3532, Case 2]{har1} and see that $\mathfrak{l}_{5,5}$ is also nilpotent of class $2$ and  realized as central extension of a  1-dimensional ideal $\mathfrak{i}=\mathfrak{a} \subseteq Z(\mathfrak{l}_{5,5})$ by the Lie algebra quotient  $\mathfrak{b}=\mathfrak{l}_{5,5}/\mathfrak{a}=\mathfrak{h}(1)+\mathfrak{i}$. Again Definition \ref{extensions} is satisfied. On the other hand, we may look at \cite[Pages 3531--3532]{har1} and see that $\mathfrak{l}_{5,8}$ may be  realized as central extension of a  1-dimensional ideal $\mathfrak{i}=\mathfrak{a} \subseteq Z(\mathfrak{l}_{5,5})$ by the Lie algebra quotient  $\mathfrak{b}=\mathfrak{l}_{5,5}/\mathfrak{a}=\mathfrak{h}(1)+\mathfrak{i}$. Again Definition \ref{extensions} is satisfied.  Then we consider $\mathfrak{l}_{5,6}$. We  look  at \cite[Page 4206, Case 2.1]{har} and see that it is nilpotent of class $2$ and realized as central extension of a 1-dimensional ideal  $\mathfrak{i}=\mathfrak{a} \subseteq \mathfrak{l}_{5,6}=\mathfrak{l}_{4,3}$. Definition \ref{extensions} is satisfied and, in the same way, we can construct  $\mathfrak{l}_{5,7}$ and $\mathfrak{l}_{5,9}$ as central extensions with central ideal of dimension one. These three central extensions of a central ideal of dimension one by $\mathfrak{l}_{4,3}$ are not equivalent, as indicated by the computations at \cite[Page 4206, Case 2.1]{har}. The result follows.
\end{proof}

 We will quote now some recent results on Schur multipliers of Lie algebras, in order to understand numerically the size of $M(\mathfrak{l})$

\begin{thm}[See \cite{nr1}, Lemmas 2.4 and 2.5] \label{ab} \

\begin{itemize}
\item[(i)]$\mathrm{dim} \ M(\mathfrak{h}(1))=2$.
\item[(ii)]$\mathrm{dim} \ M(\mathfrak{h}(m)) =2m^2-m-1$ for all $m\geq 2$.
\item[(iii)]A nilpotent finite dimensional Lie algebra $\mathfrak{l}$ of dimension $n$ has
$$\mathrm{dim} \ M(\mathfrak{l}) \le \frac{1}{2}n(n-1)$$
and the bound is exactly $\frac{1}{2}n(n-1)
$ if and only if $\mathfrak{l}$ is abelian.
\end{itemize}
\end{thm}

The condition (iii) above is due to Moneyhun \cite{mon} and allows us to introduce the following notion, which has been studied in \cite{nr3}.

\begin{defn}\label{corank} Given a nilpotent Lie algebra $\mathfrak{l}$ of dimension $n$, there exists a (unique) integer  $t(\mathfrak{l})\geq 0$ such that
\begin{equation}
t(\mathfrak{l})=\frac{1}{2}n(n - 1) - \mathrm{dim}\ M(\mathfrak{l}).
\end{equation}
This number is called \textit{corank} of $\mathfrak{l}$ (in the sense of Berkovich--Zhou).
\end{defn}
The  invariant of Definition \ref{corank} is originally due to Berkovich \cite{berkovich} and Zhou \cite{zhou}, who  introduced it for finite $p$-groups. It clearly gives information on the gap between the size of $M(\mathfrak{l})$ and the critical value $\frac{1}{2}n(n - 1)$ in which we find abelian Lie algebras only (see Theorem \ref{ab}). The same information can be rephrased in the context of groups as well. Of course, the case of Lie algebra is more rich than the case of finite $p$-groups and we have more instruments than the group case; this has allowed to describe the structure of all nilpotent Lie algebras in \cite{har} and \cite{har1} with $t(\mathfrak{l}) = 3, 4, 5, 6, 7 ,8$.

 \begin{thm}[Classification of Finite Dimensional Nilpotent Lie algebras via $t(\mathfrak{l}) \le 6$]\label{tl} Let $\mathfrak{l}$ be an $n$-dimensional nilpotent Lie algebra. Then
\begin{itemize}
\item[(i)]$t(\mathfrak{l})=0$ if and only if $\mathfrak{l}$ is abelian;
\item[(ii)]$t(\mathfrak{l})=1$ if and only if $ \mathfrak{l} \simeq \mathfrak{h}(1)$;
\item[(iii)]$t(\mathfrak{l})=2$ if and only if $ \mathfrak{l}  \simeq \mathfrak{l}_{4,2} \simeq \mathfrak{h}(1) \oplus \mathfrak{i}$;
\item[(iv)]$t(\mathfrak{l})=3$ if and only if $ \mathfrak{l} \simeq \mathfrak{l}_{5,2} \simeq \mathfrak{l}_{4,2} \oplus \mathfrak{i}$;
\item[(v)]$t(\mathfrak{l})=4$ if and only if either $ \mathfrak{l} \simeq \mathfrak{l}_{5,2} \oplus \mathfrak{i}$, or $ \mathfrak{l} \simeq  \mathfrak{l}_{4,3}$ or $ \mathfrak{l} \simeq \mathfrak{l}_{5,8}$;
\item[(vi)] $t(\mathfrak{l})=5$ if and only if $ \mathfrak{l} \simeq \mathfrak{l}_{7,2} \simeq \mathfrak{l}_{5,2}\oplus \mathfrak{i} \oplus \mathfrak{i}$ or $ \mathfrak{l} \simeq \mathfrak{h}(2)$.
\item[(vii)] $t(\mathfrak{l})=6$ if and only if either  $ \mathfrak{l} \simeq \mathfrak{l}_{4,2}\oplus \mathfrak{i}$, or $\mathfrak{l} \simeq \mathfrak{l}_{5,5}$, or $ \mathfrak{l} \simeq \mathfrak{h}(2) \oplus \mathfrak{i}$, or $\mathfrak{l} \simeq \mathfrak{l}_{5,8} \oplus \mathfrak{i}$, or $\mathfrak{l} \simeq \mathfrak{l}_{8,2} \simeq \mathfrak{l}_{7,2} \oplus \mathfrak{i}$;
\end{itemize}
\end{thm}

There are a series of open problems which prevent us to have a more general approach for higher dimensions.

\begin{op}The classification for $t ( \mathfrak{l} ) \ge 9$ is not known in general. There are some partial classification, imposing restrictions of various nature on $\mathfrak{l}$.
\end{op}

On the other hand, \cite{har} shows that for  $t ( \mathfrak{l} ) \in \{7,8\}$ we have nilpotent Lie algebras of dimension higher than $5$ for which an argument as in  Lemma \ref{technical} is not available. Therefore the main results of this paper will show that it is possible to realize physically the above nilpotent Lie algebras using concepts of quantum physics. This is the first time that the classification above may be interpreted in terms of concrete dynamical systems, and we will see that this approach is so sophisticated that may cover in fact all the Lie algebras of Theorem \ref{tl}.

\section{Main results}
In the present section we mainly refer to pseudo-bosons, considering bosons as a particular case of those and assume the reader is familiar with the terminology which is used in \cite{baginbagbook, a2, bagBS, bagnew1, bagnew2, barnett, gold, schu}.

In particular, one needs to have clear the notions of

\begin{itemize}
\item[-] bounded operators in Hilbert spaces (see \cite[Definition 3.1]{bagrus2018})
\item[-]Riesz basis (see \cite{schu}),
\item[-] $\mathcal{G}$-quasi basis in \cite[Definition 3.3]{bagrus2018},
\item[-]canonical commuting relations (see \cite{baginbagbook}),
\item[-] distribution (see \cite[Eq. (3.2)]{bagrus2018}),
\item[-]Assumption $\mathcal{D}$-pb 1,  Assumption $\mathcal{D}$-pb 2, Assumption $\mathcal{D}$-pb 3, Assumption $\mathcal{D}$-pbw 3, Assumption $\mathcal{D}$-pbs 3, reported in \cite[Section 3]{bagrus2018}.
\end{itemize}

The above concepts are strictly related to the field of the functional analysis and allow us to construct the framework  in which pseudo-bosons may form the structure of a Lie algebra, which is described in the previous section in abstract.

\begin{thm}\label{main1}There is at least one family of  pseudo-bosonic operators and a nilpotent Lie algebra $\mathfrak{l}$ of finite dimension for arbitrary values of $t(\mathfrak{l})$.
\end{thm}

\begin{proof} About the notations and the terminology for pseudo-bosonic operators, we will use the same in \cite[Section 3]{bagrus2018}. We begin to consider the nilpotent Lie algebra \begin{equation}
\mathfrak{l} = \mathfrak{h}(m) \oplus \underbrace{\mathfrak{i} \oplus \ldots \oplus \mathfrak{i}}_{k-\mbox{times}}.
\end{equation}
This Lie algebra has dimension $(2m+1)+k$, where $m \ge 1$ and $k \ge 0$ are arbitrary integers.

If $k=0$, then $\mathfrak{l} \simeq \mathfrak{h}(m)$. In this case a possible example  can be constructed with $m$-dimensional bosonic operators, and with their pseudo-bosonic versions.

Let $x_j$, $j=1,2,\ldots,m$ be $m$ position operators, and let $p_j=-i\frac{d}{dx_j}$ the related momentum operator. Then the generators of $\mathfrak{h}(m)$ can be written as follows:
$$
v_{2j-1}=\frac{1}{\sqrt{2}}(x_j+ip_j), \qquad v_{2j}=\frac{1}{\sqrt{2}}(x_j-ip_j),
$$
$j=1,2,\ldots,m$, with $v_{2m+1}=\mathbb{I}$, the identity operator. There are just lowering and raising operators associated to an $m$-dimensional quantum harmonic oscillator, or (alternatively) to $m$ modes of bosons.

A different representation can also be constructed easily in terms of pseudo-bosonic operators, simply by considering a shifted version of the $v_j$, that is,

$$w_{2j-1}=v_{2j-1}+\alpha_j \mathbb{I}, \qquad  w_{2j}=v_{2j}+\beta_j \mathbb{I},$$
$j=1,2,\ldots,m,$, for any scalar $\alpha_j, \beta_j \in \mathbb{C}$ with $\alpha_j\neq\overline{\beta_j}.$

On the other hand,  we may consider the genuine abelian case $\mathfrak{l}$ of dimension $k$ putting $m=0$. Therefore
\begin{equation}\label{heisenberg+abelian}
\mathfrak{l}=\langle v_1, v_2, \ldots, v_{2m}, v_{2m+1}, \ldots, v_{2m+k+1} \ | \ [v_{2j-1},v_{2j}]=\mathbb{I} \  \ \forall j =1, 2, \ldots, m,
\end{equation}
$$ \mbox{and} \ v_{2m+1}= v_{2m+2}= \ldots = v_{2m+k+1}=\mathbb{I} \rangle. $$

This means that there is no loss of generality in assuming both $m \ge 0$ and $k\ge 0$.

From \cite[Theorem 4.2, Page 44]{nr3}, we can see that
\begin{equation}
t(\mathfrak{l}) = \mathrm{dim} \ \mathfrak{l} = 2m+k+1,
\end{equation}
where $m$ and $k$ can grow arbitrarily.
This allows us to conclude that $\mathfrak{l}$ is always realized by pseudo-bosons and its corank may be arbitrary big, so the result follows.

\end{proof}

Now we provide concrete examples for the Lie algebras with $t(\mathfrak{l}) \le 6$.

\begin{thm}\label{main2}All the nilpotent Lie algebras of Theorem \ref{tl} are realized by pseudo-bosonic operators.
\end{thm}

\begin{proof}
We begin to consider $t(\mathfrak{l})$ for the values $0,1,2,3,4$. Looking at Theorem \ref{tl} (i), (ii), (iii), (iv) and (v), we identify just two cases which do not fit into the form of the Lie algebra \eqref{heisenberg+abelian}. These are $\mathfrak{l}_{4,3}$ and $\mathfrak{l}_{5,8}$. We are going to prove that both of them may be described in terms of pseudo-bosonic operators. About the notations and the terminology for pseudo-bosonic operators, we will use again that  in \cite[Section 3]{bagrus2018}.

Let $c$ be a bosonic operator such that $[c,c^\dagger]=cc^\dagger-c^\dagger\,c=\mathbb{I}$.  We define
\begin{equation}\label{add1}
v_1=c,\quad v_2=\frac{1}{2}\,(c^\dagger)^2,\quad v_3=c^\dagger,\quad v_4=\mathbb{I}.
\end{equation}
It is easy to check that all the commutators are zero, except that
$$ [v_1,v_2]=v_3, \qquad [v_1,v_3]=v_4.$$
Hence we get $\mathfrak{l}_{4,3}$. These operators can be used in quantum optics to construct two-photons Hamiltonians which are relevant, for instance, in connection with squeezed states, \cite{barnett}. The Hamiltonian is $H_0=\omega c^\dagger c+i(g{c^\dagger}^2-\overline{g}c^2)$ which, in terms of $v_j$ above can be written as
$$
H=\omega v_3v_1+i\left(2gv_2-\overline{g}v_1^2\right).
$$
Here $\omega$ is a real quantity related to the free energy of the bosons, while the complex $g$ {\em measures} the degree of squeezing of the system, \cite{barnett}.

Another example of the same algebra can be constructed using pseudo-bosonic operators $a$ and $b$ satisfying $[a,b]=\mathbb{I}$. In this case we put
$$
v_1=a,\quad v_2=\frac{1}{2}\,b^2,\quad v_3=b,\quad v_4=\mathbb{I}.
$$ They also describe $\mathfrak{l}_{4,3}$.

Of course, the same Hamiltonian $H$ as above could be used simply replacing $v_j$ with their pseudo-bosonic counterparts. This could be relevant in the analysis of bi-squeezed states, which extend squeezed states to a non hermitian context, as bi-coherent states do for coherent states.

Now we pass to describe $\mathfrak{l}_{5,8}$.
In this case one single bosonic or pseudo-bosonic operator is not enough. For this reason, let us consider two bosonic operators $c_1$ and $c_2$ such that $$
[c_k,c_l^\dagger]=\delta_{k,l}\mathbb{I}, \qquad [c_k,c_l]=0,
$$
where $\delta_{k,l}$ denotes the well known symbol of Kronecker. Then, defining
$$
v_1=c_1,\quad v_2=c_1^\dagger,\quad v_3=c_1^\dagger\,c_2,\quad v_4=\mathbb{I}, \quad v_5=c_2,
$$
the only non zero commutators are the following:
$$
[v_1,v_2]=v_4, \qquad [v_1,v_3]=v_5.$$
Hence we recover $\mathfrak{l}_{5,8}$. The same result can be found if we use two pseudo-bosonic operators $[a_1,b_1]=[a_2,b_2]=\mathbb{I}$, all the other commutators being zero. In this case we have to define
$$ v_1=a_1,\quad v_2=b_1,\quad v_3=b_1\,a_2,\quad v_4=\mathbb{I}, \quad v_5=a_2.
$$
As for a physical appareance of these (bosonic or pseudo-bosonic) operators we could consider the operator
 $$
 H=\lambda(v_3+v_2v_5),
 $$
for real $\lambda$. $H$ is manifestly not self-adjoint either if we use bosonic or pseudo-bosonic operators. The meaning of $H$ is simple: it describes interactions between the two modes of particles. Notice that no kinetic term is present in $H$. This is not strange, in concrete models, since quite often the most relevant part of an Hamiltonian is that term which describes the interactions.

Then we pass to consider the case of $t(\mathfrak{l})=5$, looking at Theorem \ref{tl} (vi). Again here we may realize pseudo-bosonic operators via the Lie algebra \eqref{heisenberg+abelian}.

Passing to the case $t(\mathfrak{l})=6$,  Theorem \ref{tl} (vii) shows that $\mathfrak{l}_{5,5}$ and $\mathfrak{l}_{5,8} \oplus \mathfrak{i}$ are the cases not involved in \eqref{heisenberg+abelian}. While $\mathfrak{l}_{5,8} \oplus \mathfrak{i}$ may be written easily in terms of pseudo-bosonic operators adding an operator which commutes with all those of $\mathfrak{l}_{5,8}$, an explicit construction for $\mathfrak{l}_{5,5}$  must be given.

For that, we can invoke Lemma \ref{technical} and see that it is possible to construct this Lie algebra with a family of pseudo-bosonic operators realizing  the Lie algebra $\mathfrak{h}(1) \oplus \mathfrak{i}$ and a central ideal $\mathfrak{a}=\mathfrak{i}$. Therefore $\mathfrak{l}_{5,5}$ is the central extension of $\mathfrak{a}$ by $\mathfrak{b}=\mathfrak{l}_{5,5}/\mathfrak{a}$. This is a first possible abstract approach.

A different realization of $\mathfrak{l}_{5,5}$ can be constructed using again two modes of bosonic (or pseudo-bosonic) operators $c_1$ and $c_2$ as for $\mathfrak{l}_{5,8}$. Now we define
$$
v_1=c_1,\quad v_2=c_2+\frac{1}{2}(c_1^\dagger)^2,\quad v_3=c_1^\dagger,\quad v_4=c_2^\dagger,\quad v_5=\mathbb{I}.
$$
The only non zero commutators are $[v_1,v_2]=v_3$, and $[v_1,v_3]=[v_2,v_4]=v_5$. These commutators can also be recovered using pseudo-bosonic operators $(a_1,b_1)$ and $(a_2,b_2)$. In this case we define
$$
v_1=a_1,\quad v_2=a_2+\frac{1}{2}b_1^2,\quad v_3=b_1,\quad v_4=b_2,\quad v_5=\mathbb{I}.
$$
A possible Hamiltonian constructed out of these last operators is the following:
$$
H=\omega v_3v_1+\lambda v_4v_2,
$$
which is manifestly non self-adjoint, and describes a sort of "always increasing"
 interaction between two modes, without any decreasing.
\end{proof}

As we have seen in the proof of the previous theorem, the notion of central extension in Definition \ref{extensions} is appropriate to describe nilpotent Lie algebras in terms of pseudo-bosons. In fact we can describe all those of dimension at most five.

\begin{cor}\label{dimension}
All the nilpotent Lie algebras of Theorem \ref{classification} may be realized by central extensions of families of pseudo-bosonic operators.
\end{cor}

\begin{proof}
Using Lemma \ref{technical} and the fact that \eqref{heisenberg+abelian} is described by pseudo-bosons, the result follows.
\end{proof}

 The conclusions of the above corollary may be rephrased in a
 more significant way, showing the the theory of the
 pseudo-boson operators is really a powerful tool both in
 mathematics and physics.

 \begin{cor}\label{final}
 If there exists a finite dimensional nilpotent Lie algebras
 which cannot be realized by central extensions of families
 of pseudo-bosonic operators, then its dimension must be
 $\ge 6$.
 \end{cor}

\section*{Acknowledgements}
F.B. acknowledges partial support from Palermo University and from the Gruppo Nazionale di Fisica Matematica (GNFM) of the Istituto Nazionale di Alta Matematica (INdAM).