Genus expansion of matrix models and (cid:126) expansion of B KP hierarchy

We continue the investigation of the connection between the genus expansion of matrix models and the (cid:126) expansion of integrable hierarchies started in [APSZ20]. In this paper, we focus on the B KP hierarchy, which corresponds to the inﬁnite-dimensional Lie algebra of type B . We consider the genus expansion of such important solutions as Br´ezin-Gross-Witten (BGW) model, Kontsevich model, and generating functions for spin Hurwitz numbers with completed cycles. We show that these partition functions with inserted parameter (cid:126) , which controls the genus expansion, are solutions of the (cid:126) - B KP hierarchy with good quasi-classical behavior. (cid:126) - B KP language implies the algorithmic prescription for (cid:126) -deformation of the mentioned models in terms of hypergeometric B KP τ - functions and gives insight into the similarities and diﬀerences between the models. Firstly, the insertion of (cid:126) into the Kontsevich model is similar to the one in the BGW model, though the Kontsevich model seems to be a very speciﬁc example of hypergeometric τ -function. Secondly, generating functions for spin Hurwitz numbers appear to possess a diﬀerent prescription for genus expansion. This property of spin Hurwitz numbers is not the unique feature of B KP: already in the KP hierarchy, one can observe that generating functions for ordinary Hur-witz numbers with completed cycles are deformed diﬀerently from the standard matrix model examples.


Introduction
The theory of matrix models has a long history and extensive applications in physics, mathematics, and many other fields of knowledge.Some motivations, history, and applications can be found, for example, in [Mor94;Mir94;DGZ95;EKR15].This theory has deep connections with integrable systems because matrix models form a large class of solutions (τ -functions) of integrable hierarchies of KP/Toda type [DJKM82;JM83].In this paper we continue to develop deeper connections between the two, i.e. we study the similarity between genus expansion, coming from the matrix model side, and expansion of integrable hierarchies on the other side.
In the literature, genus expansion is also known as large N expansion which goes back to 't Hooft [Hoo74].It comes from the perturbative expansion of matrix integrals as a sum over ribbon graphs.Each connected graph comes with a factor of N 2−2g , where g is a genus of the surface where it can be drawn and N is the size of matrices in the ensemble.That is, in the connected part of the partition function (free energy) one can distinguish contribution from surfaces of different genera (1) Large N factorization of correlators allows one to think of the large N limit of matrix models as a quasi-classical expansion.From different points of view, it can be also understood as a perturbative calculation of string amplitudes or as a WKB approximation.
In the presence of an external field in the action things get slightly more complicated [KMMM95].In general, the partition function depends both on the traces of the external matrix and on N .But quite often the models are defined to be independent on N , such that an explicit factor N 2−2g does not appear.Nevertheless, genus expansion for such partition function still exists.Usually, they can be viewed as generating functions of some geometric quantities (for example, integrals over moduli spaces), which provide the genus expansion of the initial matrix integral.One can distinguish contributions from surfaces of different genera with the help of the formal parameter .The free energy is then (2) Such an expansion generalizes (1) by taking = 1/N and the transition from the partition function to free energy given by F = 2 log Z .The higher genus part of the expansion can be calculated from genus zero one-and two-point functions with the help of spectral curve topological recursion procedure [AMM04; AMM05; CE06; EO08].
From the side of integrable hierarchies, one can independently insert a formal parameter into equations.At this moment has nothing to do with the genus expansion parameter yet.Such a "deformation" of the KP hierarchy was considered in [TT95] in order to study the dispersionless limit → 0 of the hierarchy.The parameter is inserted in such a way that the standard KP hierarchy is restored at = 1.At the first glance, such a deformation of the hierarchy seems trivial: if one looks only at the equations of the hierarchy, then the insertion of is given just by the rescaling of "times": t k → t k .To satisfy such -deformed equations one needs simply to rescale times in the τ -functions in the same way.However, one can notice that can also be inserted, for example, in the Plücker coefficients so that the -KP equations are still satisfied.Thus, the single KP τ -function can be -deformed to satisfy the -KP equations in many different ways.
Since was inserted to study the dispersionless limit, it is natural to consider only those -deformations of τ -functions for which the → 0 limit exists.That is, the free energy should contain only non-negative powers of .And the trivial rescaling of times is not enough to satisfy this condition: 2 log τ t would be very singular with respect to as → 0. That is, properly deformed τ -functions should contain non-trivially.In [TT95] it was shown which insertion of into the solutions is required in order to possess the good quasi-classical limit (see also [APSZ20]).Even so, -deformation of the τ -functions is still not unique under this additional condition.Nevertheless, we are interested only in the one particular -deformation of the τ -function which coincides with its genus expansion.
In [APSZ20] several particular examples of matrix model KP τ -functions were considered.Among the examples were the Gaussian Hermitian model, Brézin-Gross-Witten model, the Kontsevich model, and the generating function for simple Hurwitz numbers.All the models were deformed by the insertion of the parameter responsible for the genus expansion.By definition of the genus expansion, they have good quasi-classical behavior of the form (2). It was shown that, firstly, -deformed models are τ -functions of -deformed KP hierarchy.That is, the genus expansion parameter from the matrix model side precisely coincides with the quasi-classical parameter on the integrable system side.Secondly, since almost all the mentioned partition functions (except for the Kontsevich model and BGW model in the Kontsevich phase) belong to the one family of hypergeometric τ -functions, the general prescription for insertion into these models was obtained.
Let us consider this algorithmic prescription for deformation in more detail.The hypergeometric family of KP τ -functions has the following expansion in Schur polynomials [OS01]: where f λ = (i,j)∈λ f (i − j).Here i, j are the coordinates of all the boxes in the Young diagram λ, and f (n) is an arbitrary function.Choosing certain f (n) gives a particular representative of the family.Let us perform the following insertion of parameter : Such a prescription appears to perform the genus expansion for the Gaussian Hermitian model, simple Hurwitz numbers, and the BGW model in the character phase.Naively, the genus expansion is very specific for a certain model.Nevertheless, it was shown that there exists a common pattern for it.However, one can see that the Kontsevich model and BGW model in the Kontsevich phase do not fit into this elegant picture in the context of KP hierarchy because they do not belong to hypergeometric KP τ -functions.The recent progress in the investigation of superintegrability in these models [MM21; Ale21; Ale23] revealed that they belong to the hypergeometric family of BKP τ -functions.The BKP hierarchy, which is similar to the KP hierarchy, has the parametrization of solutions by an infinite-dimensional algebra of type B. That is why one may seek for the natural -deformation of these models in the context of BKP hierarchy.
In this paper, firstly, we generalize the -deformation approach of [APSZ20] to the BKP hierarchy in accordance with [Tak93] and obtain the simple algorithmic prescription of -deformation for BGW and Kontsevich models in terms of family of hypergeometric BKP τ -functions.That is, given that the family has the following expansion in Q-Schur polynomials [Orl03] (see also section 3.1): where r λ = (i,j)∈λ r(j), we obtain the following prescription for -deformation: We show that genus expansion of both the Kontsevich model and BGW model in the Kontsevich phase is governed by this prescription of -deformation, which implies that they are -BKP solutions.Also, we explain why a shift on 1/2 is necessary in the -BKP case.Secondly, we consider -deformation of another set of examples from BKP hypergeometric family, i.e. generating functions for spin Hurwitz numbers with completed cycles.These functions receive a lot of attention recently [GKL21; AS21; MMN20; MMNO21; MMZ21; MMZ22] and they require special consideration.First of all, it is important to mention that the genus expansion of these models is done in the following way.The partition function is rewritten as a generating function of integrals over moduli spaces of curves with the help of cohomological representation motivated by Gromov-Witten theory [GKL21].This implies that several Hurwitz numbers, possibly counting surfaces of different genera, contribute to one such integral.The contribution of a certain genus g to the partition function comes from the integrals over moduli spaces M g,n rather than from Hurwitz numbers themselves.Now, despite the fact that these τ -functions belong to the hypergeometric family, their genus expansion is not governed by (6).Instead, they have their own prescription for -deformation.While the non-deformed spin Hurwitz τ -functions look like where Casimirs (or, r-th completed cycles) are p r (λ) = (λ) i=1 λ r i , their -deformation prescription is given by To compare this to the deformation (6), one needs to rewrite e urpr(λ) = (i,j)∈λ e urpr(j) , which implies rewriting the -deformation (8) as: We show that p r (n) = (n r − (n − 1) r ), thus, prescription (9) is different from (6).Nevertheless, we show that -deformed spin Hurwitz partition functions are still solutions of -BKP with good quasi-classical limit (2).Such a discrepancy between the prescription of -deformation tends us to revisit the KP case (4) and consider generating functions for ordinary Hurwitz numbers with completed cycles, which were not mentioned in [APSZ20].The proper -deformation of such partition functions was considered, for example, in [BDKS20] and, indeed, we observe the same picture.While all the "simple" examples of [APSZ20] are deformed by (4), ordinary Hurwitz numbers with completed cycles have their own -deformation prescription.The only standing out example is the generating function for simple Hurwitz numbers for which both prescriptions coincide.That is, this discrepancy in -deformation of Hurwitz numbers appears both in KP and BKP cases.
The reason for the difference in the -deformation lies in the form of the function f (n) (or r(n)).Starting from the KP/BKP point of view the difference is not visible.However, the theory of [BDKS20] gives a recipe for deformation starting from the spectral curve.And for spectral curves corresponding to Hurwitz numbers, one can see the more complicated way of deformation than for the other simple examples.The function f (n) (or r(n)) appears to be tightly connected with the spectral curve data, which implies the difference in the genus expansion.Thus, the deformation recipe is sensitive to either function f (n) (or r(n)) is polynomial or exponential, as it is in BGW and spin Hurwitz examples respectively.From this point of view, the Kontsevich model stands aside because of the very special form of r(n), depending on n mod 3 (see section 5.2).Even so, we observe that the Kontsevich model indeed deforms as simply as the BGW model.And this fact suggests that there exists some generalization of [BDKS20] on the BKP case which should include the Kontsevich model as a special case.
The paper is organized as follows.In section 2 we introduce the notations which we use throughout the paper, such as Q-Schur polynomials, algebra go(∞) and neutral fermions.Section 3 is devoted to the BKP hierarchy and, in particular, to the important class of hypergeometric solutions which arises in examples.In section 4 we consider -formulation of BKP hierarchy.In section 5 we explicitly introduce such solutions of -BKP as BGW and Kontsevich models.For each model we explain the insertion of the parameter and prove that they are, indeed, solutions of -BKP.Moreover, we show the common pattern of insertion of into these models in terms of hypergeometric τ -functions.Finally, in section 6 we consider generating functions of both ordinary and spin Hurwitz numbers.Firstly, we revisit the KP case and discuss the difference in -deformation of these functions with the other examples.Then we switch to the BKP case and observe the same feature for the spin Hurwitz numbers.

BKP boson-fermion correspondence
In this section, we introduce the language that we are using throughout the paper: we briefly review the most important facts about Q-Schur polynomials, infinite-dimensional Lie algebra of type B, and boson-fermion correspondence in the context of BKP integrable hierarchy.The latter allows one to identify bosonic and fermionic descriptions and to use the more convenient one.Extensive information about Q-Schur polynomials can be found in [Mac98], here we summarize only what we do need.

Q-Schur polynomials
Firstly, let us introduce the Q-Schur polynomials.Let us consider an ordered set of nonnegative integers λ 1 ≥ λ 2 ≥ • • • ≥ λ ≥ 0. We denote this set by λ = [λ 1 , λ 2 , . . ., λ ] and call it a Young diagram.Each Young diagram corresponds to a partition of an integer |λ| = λ 1 + λ 2 + • • • + λ onto (λ) non-zero parts λ i .Graphical representation of Young diagrams is a finite collection of boxes, arranged in left-justified rows, with the length of each row equal to λ 1 , λ 2 , . . ., λ .For example, the diagram [5, 3, 2] has the following graphical representation: Consider an infinite set of variables with odd indices t = {t 1 , t 3 , . ..}. Q-Schur polynomials Q λ (t) are labeled by Young diagrams and are defined via the pfaffian formula: where antisymmetric matrix M λ (t) is defined as for even (λ); if (λ) is odd, add exactly one line P 0,λ i (t) (as if one adds zero length line to the Young diagram λ).In turn, polynomials P n,m (t) are evaluated with the help of the following generating function: This definition implies that the Q-Schur polynomials vanish if λ contains two lines of equal length.The diagrams that do not have lines of equal length are called strict partitions (SP).
To be precise, let us introduce the first few non-trivial examples of Q-Schur poly-nomials: 2.2 Fock space, neutral fermions and go(∞) algebra Again, here we summarize only the necessary results.The original approach to neutral fermion Fock space is described in great detail in [DKM81; DJKM82; You89; vdLeu95; Orl03].In our notations, we follow a more recent summary of [Ale23].There is a natural way to describe solutions of BKP hierarchy in terms of neutral fermions.Firstly, let us introduce an infinite-dimensional Clifford algebra with generators φ k , k ∈ Z and commutation relations: Note that φ 2 0 = 1/2.One can make neutral fermions from regular free fermions Let us introduce generating series for neutral fermions: Neutral fermion Fock space F (and its dual F * ) is defined by the action of Clifford algebra on the vacuum vector |0 (and respectively corresponding co-vacuum vector 0|): and the elements φ The linear space F splits into two subspaces where F 0 and F 1 denote the subspaces with even and odd number of generators respectively.If we consider only the space F 0 , the basis in it can be labelled by strict partitions λ ∈ SP in the following way: With the help of commutation relations (14) it is easy to see that where We denote the normal ordering of fermionic operators as :(. ..):, which means that all annihilation operators are moved to the right and all creation operators to the left, with respect to (−1) with each transposition of fermions.For example, :φ −2 φ 1 : = −φ 1 φ −2 .Note that it is not the same as the transposition of fermions with the help of commutation relations (14).
Let's consider a Lie algebra of matrices go(∞): each matrix A ∈ go(∞) is an infinite matrix, with two additional requirements: • only finitely many diagonals are non-zero, Let us introduce matrices E ij that have 1 on i, j's place and 0 everywhere else, i. e. (E ij ) kl = δ ik δ jl .Then the standard basis for the algebra go(∞) consists of matrices It is easy to show that F ij satisfy standard commutation relations Now let us consider go (∞)'s central extension, go(∞).As a linear space, it is go(∞) ⊕

Cc. The commutator of two arbitrary elements
where α(A, B) is linear in each variable and therefore can be defined on the basis elements F ij : where 1 is the identity element.Now the go(∞) representation on space F is, in terms of basis elements: It is straightforward to check that r is indeed a representation: commutation relations (25) for r(A), r(B) are preserved.The central charge c in this representation is equal to 1.

Now let us consider operators
which satisfy the following commutation relations and thus generate the Heisenberg subalgebra A ∈ go(∞).
It appears that one can construct an isomorphic representation of go(∞) in bosonic Fock space with the help of maps Φ i : F i → B i = C t 1 , t 3 , t 5 , . . .for i = 0, 1, which are homomorphisms of representations [You89].Under this homomorphism the operators H k map into multiplication and differentiation w. r. t. times Maps Φ i provide the BKP variant of boson-fermion correspondence.In what follows we consider BKP hierarchy, thus we are interested only in the action of GO(∞) on vacuum vector |0 .Group GO(∞) is a standard exponential map from the algebra go(∞).Maps Φ i , in this case, can be explicitly written as vacuum expectation value: where 1| = √ 2 0| φ 0 and The explicit isomorphism between spaces F 0 and B 0 is given by That is, the boson-fermion correspondence provides an explicit realization for the states (20) in terms of Q-Schur polynomials.

BKP hierarchy
In this section, we briefly review the main facts about BKP equations and their solutions.For a detailed explanation see the original papers [DJKM82; JM83; You89] or [Orl03].BKP hierarchy is an infinite set of non-linear differential equations with the first equation given by It is more common to work with τ -function τ (t) = exp F (t).We assume that τ (t) is at least a formal power series in times t k , and maybe it is even a convergent series.
The entire set of equations of the hierarchy can be written in terms of τ -function using Hirota bilinear identity where and Contour integration dk 2πi here means that we expand integrand at the point k = ∞ and take the coefficient of k −1 .
With the change of variables t j = t j + j , t j = t j − j one can rewrite bilinear identity in terms of Hirota derivatives Expanding integrand in powers of j and taking the coefficient of k −1 one obtains BKP equations.
On the one hand, all formal power series solutions of BKP hierarchy can be decomposed over the basis of Q-Schur polynomials Function written as a formal sum over Q-Schur polynomials is a BKP solution if and only if coefficients C λ satisfy the BKP Plücker relations: The first non-trivial relation is We call τ -functions of the form (39) as τ -functions in bosonic representation.
On the other hand, τ -function is an image under the boson-fermion correspondence of a point on the orbit of the vacuum |0 under the action of some element G of GO(∞): Such τ -functions are called τ -functions in fermionic representation.In what follows we use both representations for BKP τ -functions.

Hypergeometric τ -functions
In this paper, we are interested in a subset of BKP τ -functions of hypergeometric type, or simply hypergeometric τ -functions.This relatively simple set of BKP solutions contains surprisingly many physical examples.It was first introduced in [Orl03].In fermionic representation these τ -functions have the form: where function r(n) has to satisfy and β = {β 1 , β 2 , . ..} is an arbitrary set of parameters.Matrix B k has non-zero elements of a specific form on the k-and (−k)-th diagonals.It is easy to obtain another form of the matrix B k which is more convenient in some cases: where D = z∂/∂z, therefore r(D)z n = r(n)z n .Using the explicit form of the fermionic fields (17) one can obtain (45).Bosonic representation of hypergeometric τ -functions requires a notion of a BKPcontent c(w) of a box w of Young diagram λ: For example, boxes of the diagram [5, 3, 2] have the following contents: Hypergeometric τ -functions in bosonic representation are given by where

-formulation of BKP hierarchy
In this section, we introduce a formal parameter in the BKP hierarchy (we shortly call it -BKP).The idea to study -BKP was first formulated in [Tak93] in order to investigate the dispersionless limit of the hierarchy.In our work we follow the way analogous to -formulation of the KP hierarchy in [TT95; TT99; NZ16].Let us define -BKP hierarchy as an ordinary BKP hierarchy with rescaled variables t k → t k / , and redefined F = 2 log τ .The first equation of the hierarchy is then of the form As one can see, it is not a "deformation" in any sense, it is just simply a rescaled original BKP hierarchy that allows F -functions to depend on arbitrary powers of formal parameter .One can obtain τ -functions of -BKP from BKP ones with the change of variables t k → t k / .
Non-triviality appears when we restrict free energy F not to be singular in .That allows to perform dispersionless limit → 0. This requirement imposes restrictions on F : parameter should be inserted into the logarithm of τ -function and equations "properly", that is, F -function must not contain any negative powers of : The dispersionless free energy should satisfy the dispersionless hierarchy with the first equation given by The non-triviality in -dependence of F -function can be obtained by inserting in the Plücker coefficients as well as a rescaling of times: Proposition ( -BKP solution criterion).τ -function of the form (52) satisfies -BKP equations if and only if coefficients C λ satisfy classical Plücker identities.
Proof.If τ (t) solves -BKP then τ (t ) solves BKP, therefore C λ satisfy the classical Plücker relations.If C λ satisfy the classical Plücker relations then by the same logic τ (t) solves -BKP.
In the following sections, we show that for known solutions of the BKP hierarchy, their genus expansion (given by non-trivial insertion of in C λ ), however, solves -BKP.In this sense, we truly deform our original BKP τ -functions with the formal parameter and obtain -BKP τ -functions.

Examples of -BKP solutions
In this section, we discuss separately two solutions of BKP: Brézin-Gross-Witten model and Kontsevich model.For each of these τ -functions we show two things: • genus expanded τ -function satisfies the -BKP hierarchy, • genus expansion of both τ -functions can be obtained by following one simple prescription.
Even though this prescription works for the Kontsevich and BGW models, it changes for the generating functions for spin Hurwitz numbers.For now, we leave the discussion of the difference until the next section.

Brézin-Gross-Witten model
Firstly, let us discuss the Brézin-Gross-Witten model and its genus expansion.This model has two phases and both of them were already considered in [APSZ20] in the context of the KP hierarchy.It was shown that the -deformation of the model (which reveals the genus expansion) is a solution of the -KP hierarchy in both phases.While the genus expansion of the character phase is governed by the general prescription for -deformation of hypergeometric KP τ -functions (4), the Kontsevich phase does not fit into this family and therefore into this prescription of deformation.Here we treat the Kontsevich phase of the model as the simplest non-trivial hypergeometric solution of the BKP hierarchy and develop a prescription for -deformation of such τ -functions.In particular, we show that in this phase the deformed BGW model is a solution of the -BKP hierarchy.

BGW model was first introduced as a partition function of 2D lattice gauge theories [GW80; BG80]:
where the integration is over N × N unitary matrices with the Haar measure DU of the unitary group U (N ), normalized by DU = 1, and N → ∞.
Since the Haar measure is invariant with respect to the group action, Z BGW J, J † depends only on N parameters, the eigenvalues of the matrix JJ † .Depending on the choice of variable t k , there are two phases [MMS96]: We focus only on the Kontsevich phase since the BGW model in this phase is a solution to the BKP hierarchy and has a simple expansion in Q-Schur polynomials [Ale23] So we see that the partition function of the BGW model in the properly normalized times is a hypergeometric τ -function of the BKP hierarchy with β k = 2δ k,1 .
In fermionic formalism BGW τ -function looks like

BGW as a solution of -BKP
The BGW model is a generating function for intersection numbers of Θ-classes and ψ-classes on compactified moduli spaces M g,n of complex curves of genus g with n marked points [Nor17].These intersection numbers are rational numbers, which are not equal to zero only if Let us write the generating function with the parameter enumerating contributions of different genera: From (60) one can easily see that the genus expansion is obtained by rescaling of times: One can show from definition (10) that Q-Schur polynomials in variables t k k are homogeneous in and moreover We know that coefficients satisfy BKP Plücker relations because BGW τ -function solves BKP.Rescaling C λ by h |λ| does not change the relations, since they are homogeneous by the sum |λ satisfy the classical Plücker relations and, hence, expanded τ -function of spin Hurwitz numbers is a solution of -BKP.By the same argument as in (62) we achieve that Now we can write the -deformed τ -function as From ( 48) and (56) we have Let us define Now it seems reasonable to introduce and we are ready to see that rule allows us to get r (n) from r(n) defined in (56).That is, we see that -deformation of the BGW model is a solution of -BKP hierarchy and the prescription to the deformation is given by In fermionic formalism -BGW τ -function looks like

Kontsevich model
Next, we discuss the Kontsevich model and its genus expansion.Similarly to the BGW model, the Kontsevich phase of the model was considered in [APSZ20] in the context of the KP hierarchy.This phase does not fit into the family of hypergeometric KP solutions, hence, its deformation in terms of KP hierarchy seems odd.Natural language for the description of the Kontsevich model is the language of hypergeometric BKP τfunctions.Here we consider the deformed Konsevich model in these terms and show that it is a solution of the -BKP hierarchy.The recipe for the deformation in the context of hypergeometric BKP τ -functions is the same as for the BGW model.

Classical Kontsevich model
Let us consider the bosonic representation of the Kontsevich model.τ -function of the Kontsevich model [Kon92] in Kontsevich phase [MMS96] is defined by the matrix integral where integration is taken over hermitian matrices X.From the point of view of BKP, it depends only on odd times t k , which are just powers of Λ According to [MM21; LY22], the character expansion of the Kontsevich τ -function in the basis of the Schur Q-functions is Therefore, we have For our purposes, we have to find r(n) function for given r λ .Now let us derive it.The answer is formulated in (86).From [Ale23] we know that Therefore, we have 1 16 Now it is convenient for us to express any partition λ in the following form λ = (3k 1 , . . ., 3k p , 3m 1 + 1, . . ., 3m q + 1, 3n 1 + 2, . . ., 3n r + 2). (79) The remain part of r λ is given by the ratio of Q-Schur polynomials at the point δ k,3 (see [LY22] or [MMNO21]): Our goal is to rewrite this expression in terms of the product over the boxes of the Young diagram.It can be thought of in the following way.We can by default fill λ with 1 in each box and associate particular multipliers in the last equation with particular boxes.Each of k i , m i , n i , corresponds to the particular row of λ and it is natural for beginning to associate products up to k i or m i or n i with the corresponding row.In turn, for each picked row we can associate the particular multiplier with the particular box.We propose to fill three types of rows (associated with k i or m i or n i ) in the following way (81) All minuses are here to get the overall factor (−1) r .Because in fact the number of all minuses in the partition λ is The last equality holds since [LY22] Q-Schur polynomials Q λ (δ k,3 ) are nonzero only for diagrams with q = r.The first example of strict partition with 3n boxes, but which has q = r is [7, 4, 1] And for such diagram Schur Q-function at the point δ k,3 is zero.After we calculate the product of all minuses, we get After all Note that rows of λ (81) are filled in agreement with the function on the r.h.s.excluding power of 1/3, which is trivial to insert.
BKP-content is given by ( 46), and now we can restore After looking once again on proposed way of filling partition λ (81) (remember that q = r) the simplified function r(n) looks like In the fermionic formalism the Kontsevich τ -function has the following form

Kontsevich model as a solution of -BKP
The Kontsevich model is the generating function for intersection numbers of Chern classes on compactified moduli spaces M g,n of complex curves of genus g with n marked points.Intersection numbers of Chern classes are rational numbers, which are not equal to zero only if n i=1 Let us write the generating function with inserted parameter enumerating contributions of different genera [Ale14]: From (90) we can see that the genus expansion is obtained by rescaling times: After this rescaling one can obtain Now we have r λ = 2|λ|/3 r λ and using the same technique as in case of BGW model we claim that genus expanded Kontsevich model solves -BKP hierarchy.Next, by the same arguments as in the previous subsection we get Note that this insertion inside r(n) is true only for such λ that have q = r, and the other Young diagrams do not contribute to the partition function.That is, the recipe for deformation of hypergeometric BKP τ -functions (71) holds for the Kontsevich model as well.
In the fermionic formalism -deformed Kontsevich τ -function has the form 6 Spin Hurwitz numbers as -BKP solution and its -KP counterpart This section is devoted to the generating functions for Hurwitz numbers with completed cycles, and their spin counterpart.Ordinary Hurwitz numbers, counting ramified coverings of a Riemann surface with imposed conditions on the ramifications, were defined by Hurwitz in [Hur91; Hur01].In the more recent years, Hurwitz numbers again became an object of interest, due to strong ties with the integrable hierarchies [Oko00], Gromov-Witten theory [OP06], the intersection theory of the moduli spaces of curves via ELSV type formulae [ELSV01], topological recursion [EO08; BM08; BEMS11], and W -representation of partition functions [MMN11].
In this section, we in particular consider a type of Hurwitz numbers called spin Hurwitz numbers, introduced by Eskin-Okounkov-Pandharipande [EOP08].The defining feature of these numbers is the presence of a spin structure (or theta characteristic) on the surfaces, and the counting of coverings is weighted by the parity of this theta characteristic.We denote spin Hurwitz numbers with a superscript ϑ, to emphasize the role of the theta characteristic.
In this section we find f (n) and r(n) functions for both generating functions for ordinary Hurwitz numbers with completed cycles and their spin counterpart.We formulate the prescription of insertion into the functions f (n) and r(n).The prescription differs from (71).We explicitly check that the genus expansion of these models in terms of satisfies the -BKP hierarchy.

Hurwitz numbers with completed cycles
This part of the section is devoted to a topic that is directly connected with the KP hierarchy and its -deformation.As it was mentioned in the introduction, such a revisit of -KP examples is important in understanding the -BKP case.We refer the reader to [APSZ20] for a detailed introduction to the subject.Here we only fix some notation and immediately after that we switch to the discussion of Hurwitz numbers.

Classical ordinary Hurwitz numbers with completed cycles
Schur polynomials S λ (t) are defined with the help of generating function and determinant formula: Partition functions Z = Z(t, u) for ordinary Hurwitz numbers with completed cycles [MMN11] is given by where the symmetric sums (or Casimirs) are defined as Now let us find the f (n) function for To do this we can start with the natural ansatz Identity (KP content differs from BKP analogue (46) and equals c(w) = j − i) must hold for all λ and u k , therefore And it is easy to check that function solves this functional equation.One can see that the set of p r (n), r > 0 is the basis in space of infinitely differentiable functions and after all f (n) for Hurwitz numbers with completed cycles is an arbitrary infinitely differentiable function.Therefore, every partition function Z for Hurwitz numbers with completed cycles is hypergeometric KP τ -function and vice versa.

Ordinary Hurwitz numbers with completed cycles as a solution of -KP
Hurwitz numbers with (r + 1)-completed cycles defined by partition µ and genus g of the covering surface again can be expressed as the integrals over the moduli space of curves [SSZ15; DKPS19] where and C is the Chiodo class [Chi06].The generating function for Hurwitz numbers with (r + 1)-completed cycles and insertion is then given by where p µ = (µ) k=1 p µ k and p-variables (commonly used for Hurwitz generating functions) are related to the t-variables (commonly used for integrable hierarchies) via simple rescaling p k = kt k .The corresponding partition function with inserted is then Note that the partition function is the generating function for the disconned Hurwitz numbers h •,r g,µ , where g is defined for disconnected surfaces via (104).Now let us rearrange the terms in the -deformed partition function to bring it to the form (96).One of the most important ingredients in our calculations is the Frobenius formula [Fro96] where m k is the number of lines of length k and χ λ µ are the characters of the symmetric group.Combinatorial definition of disconnected Hurwitz numbers with (r + 1)completed cycles [OP06; SSZ12] can be used for the simplification of the partition function Next, similarly to the non-deformed case, in order to obtain the KP τ -function of the form (96) one has to go back to the t-variables For (r + 1)-completed cycles KP Plücker coefficients are So every term in Plücker relations will be of the form One can show that for the KP Plücker relations the first factor is the same for all terms.Therefore, since Schur polynomials satisfy Plücker relations, τ -function of Hurwitz numbers with (r + 1)-completed cycles solves -KP.This fact easily generalizes to the arbitrary linear combinations of completed cycles.
To summarize, Z solves -KP and, in particular, for each p r (n): We observe the following deformation rule for the ordinary Hurwitz numbers with completed cycles: and it differs from (4).Let us mention that for the simple Hurwitz numbers, considered in [APSZ20], u r = u 2 δ r,2 and the deformation rule (114) simplifies to but f (n) = f ( n) for any other choices of u with higher r.The reason for such a difference in -deformation prescriptions is in the form of function f (n).From [BDKS20] we know that f (n) is tightly connected with the spectral curve data.Since the spectral curve (in fact, the whole procedure of topological recursion) contains the information about correlators corresponding to different genera, the -deformation procedure can be encoded into this data as well.However, from this point of view, the form of the function f (n) (either it is rational or exponential) appears to be extremely important and it changes the way of -deformation.Finally, let us write the -deformed Hurwitz τ -function in the fermionic formalism: In the fermionic formalism p r (n) corresponds to the operator p r (D) and we can write the generating function of (r + 1)-completed Hurwitz numbers as where Such insertion implies that τ -function has good quasi-calssical limit according to [TT95].

Spin Hurwitz numbers with completed cycles
Finally, we consider the last family of BKP τ -functions, in full analogy with the KP case.

Classical spin Hurwitz numbers with completed cycles
Partition function Z = Z(t, u) for spin Hurwitz numbers with completed cycles [MMN20] is given by where the BKP symmetric sum is given by As for the non-spin case, we are going to find the r(n) function for Let us make an analogous ansatz The identity must hold for all λ and u r , therefore And we obtain Note that the function p r (n) differs from the function p r (n), given by (102), by the simple shift n → n − 1 2 .Now, the set of p r (n), r ∈ Z + odd is the basis in space of infinitely differentiable functions which are symmetric with respect to 1/2 and r(n) is an arbitrary function in this space.Therefore, every partition function Z for spin Hurwitz numbers with completed cycles is hypergeometric τ -function of BKP and vice versa.

Spin Hurwitz numbers with completed cycles as a solution of -BKP
Spin Hurwitz numbers with (r +1)-completed cycles also have a representation in terms of integrals over the moduli spaces of curves [GKL21; AS21] where we use definitions from (104) and C ϑ is the Chiodo class twisted by the 2-spin Witten class.Free energy and partition function for spin Hurwitz numbers with (r + 1)completed cycles are defined in full analogy with the non-spin case.The difference is that one should replace h r g,µ → h r,ϑ g,µ , µ should run over the set of odd partitions (OP) and λ over the set of strict partitions (SP).Here we need to introduce the characters of the Sergeev group ζ λ µ in a similar to the Frobenius formula way [Ser85] where δ(λ) is equal to 0 or 1 if (λ) is even or odd respectively.Note that we have the factor of 2 − 1 2 (λ) in the definiton of Q-Schur polynomials (10).Here, again, p k = kt k .Combinatorial definition for disconnected spin Hurwitz numbers with (r+1)-completed cycles [Gun16; Lee20] allows us to rewrite the partition function in the t-variables (and standard for BKP rescaling t k → t k /2) as Thus, BKP Plücker coefficients for (r + 1)-completed cycles are and from the general form of BKP Plücker relations (40) it is easy to see that every term is of the form The first factor is the same for all the terms.Using that Q-Schur polinomials solve the BKP Plücker relations, we obtain that τ -function of spin Hurwitz numbers with (r +1)completed cycles solves -BKP.The generalization to arbitrary linear combination of completed cycles is straightforward.
To summarize, Z solves -BKP and, in particular, for each p r (n): p r (n) = r−1 (n r − (n − 1) r ). (133) We observe the following deformation rule for the spin Hurwitz numbers with completed cycles: The deformation rule (134) in terms of function p r (n) is identical to the one of the ordinary Hurwitz numbers (114).However, similarly to the KP case, it differs from BGW-and Kontsevich-type deformation prescription (71) because for any choice of u but the trivial case u r = uδ r,1 .So, for the spin Hurwitz numbers we observe the same discrepancy in -deformation as it was for the ordinary Hurwitz numbers.The reasonings in this case should be similar: either rational or exponential form of the function r(n) implies different insertion into the spectral curve data, which is in agreement with [AS21].
Finally, let us write the -deformed spin Hurwitz τ -function in the fermionic formalism:

Conclusion
To summarize, the main results of the paper are: • We have considered the genus expansion of several BKP τ -functions governed by parameter .Among the examples were Kontsevich and BGW models and generating functions for spin Hurwitz numbers with completed cycles.We have shown, that all these τ -functions with inserted parameter are solutions of -BKP with the correct quasi-classical behavior.
• We have considered all the mentioned examples as members of the hypergeometric BKP family, each member of which is parametrized by a single function r(n).We have performed the -deformation in terms of insertion of into the function r(n).However, the prescription for deformation depends on the form of function.
• We have revisited the case of KP hierarchy and observed that for some members of the hypergeometric family, in particular, Hurwitz numbers with completed cycles, there is a unique deformation prescription (114) that does not coincide with (4).The reason for that difference is in the form (polynomial or exponential) of function f (n), which carries the information about the spectral curve [BDKS20].
• In the BKP case, we have observed that both Kontsevich and BGW models have the same simple pattern of -deformation while generating functions for spin Hurwitz numbers have the other one in accordance with [AS21].Reasons for the difference between BGW and spin Hurwitz numbers can be understood by analogy with the KP case.However, the case of the Kontsevich model suggests that the theory of [BDKS20] can be further generalized to include τ -functions of similar to the Kontsevich form.Now let us discuss several questions which are related to the further research: • For both KP and BKP cases the theory of [BDKS20; AS21] can be generalized to include τ -functions of form similar to the Kontsevich model.That is, to the functions r(n) defined up to n mod k.For example, it would be useful to investigate the Hermitian matrix model with cubic potential and a proper star-like choice of contour.
• The recent increasing growth of interest in W -representations of matrix models [MM22; WLZZ22; MMM + 23a; MMM + 23b] revealed a large set of new matrix models.Firstly, it would be useful to make the genus expansion of these models and see how do they fit in the picture of -deformed hierarchies.Secondly, it would be interesting to insert directly into the W -operators, which may uncover their additional structure.
• Finally, the same picture of -deformation should be extended to the CKP and DKP hierarchies, corresponding to the rest infinite-dimensional algebras.Even though these hierarchies are less studied and do not possess well-known matrix model solutions, their investigation is important in the search for the generalization to the (q, t)-KP hierarchy.