2-BKP Tau Functions & Integrable Hierarchies

• 2-BKP tau functions are core objects in integrable hierarchies, providing Pfaffian solutions to Hirota bilinear equations and linking algebraic geometry with random matrix theory.
• They expand naturally in terms of Schur Q-functions and serve as square roots of determinantal tau functions in corresponding two-component KP/Toda hierarchies.
• Constructed via neutral fermions and characterized by Pfaffian addition formulae, these tau functions underpin applications in matrix models, combinatorics, and enumerative geometry.

Two-component BKP (2-BKP) tau functions are central objects in the theory of integrable hierarchies associated with infinite-dimensional orthogonal and symplectic groups, providing Pfaffian solutions to systems of Hirota bilinear equations and encoding quantities in algebraic geometry, random matrix theory, and quantum field theory. The 2-BKP hierarchy generalizes the classical BKP hierarchy by introducing two independent sets of variables (“times”) and is intrinsically linked to representation theory, Clifford algebras, fermionic Fock spaces, and symmetric functions. 2-BKP tau functions play a pivotal role as square roots of determinantal tau functions in the corresponding two-component KP/Toda hierarchies, and have remarkable applications in matrix models, Hurwitz enumerative geometry, and singularity theory.

1. Fermionic Construction and the Hirota Bilinear Equations

The 2-BKP hierarchy is defined via neutral (B-type) fermions $\phi^{(1)}(z)$, $\phi^{(2)}(z)$ and their respective modes, with anti-commutation relations

$$[\phi^{(i)}_m, \phi^{(j)}_n]_+ = \delta_{ij} (-1)^m \delta_{m,-n}.$$

The tau function is constructed as a fermionic vacuum expectation value: $$\tau(t^{(1)}, t^{(2)}) = \langle 0|\, g\, |0\rangle$$ where $g$ is an exponential of quadratic forms in the neutral fermions, possibly mixing the two flavors. The key property is the satisfaction of an infinite system of Hirota bilinear equations, expressible as residue identities in Miwa-shifted times: $$\mathrm{Res}_{z=0} \big\{ \tau(t^{(1)} + \varepsilon^{(1)} - [z^2],\, t^{(2)} + \varepsilon^{(2)})\, \tau(t^{(1)} + [z^2],\, t^{(2)}) - \tau(t^{(1)} + \varepsilon^{(1)},\, t^{(2)} + [z^2])\, \tau(t^{(1)},\, t^{(2)} - [z^2]) \big\} = 0,$$ where shifts are in the sets of odd times $t^{(a)}_n$, $n = 1,3,5,\ldots$ (Antonov et al., 2023, Chen et al., 19 Nov 2025).

2. Symmetric Function Expansion: Schur Q-Functions and Pfaffians

A fundamental fact of BKP theory is that tau functions expand naturally in terms of projective Schur Q-functions indexed by strict partitions. For 2-BKP, this expansion is diagonal and multiplicative: $$\tau(t^{(1)}, t^{(2)}) = \sum_{\lambda\, \text{strict}} c_\lambda\, Q_{\lambda}(t^{(1)})\, Q_{\lambda}(t^{(2)}),$$ with explicit combinatorial coefficients often given by the model under consideration (e.g., $$c_{2\alpha} = c^{|\alpha|} \prod_{i=1}^{\ell(\alpha)} (\alpha_i)!$$ for the modified ZJZ model (Antonov et al., 2023)). Wick’s theorem ensures all such tau functions admit a single Pfaffian formula in Miwa variables, and in bosonic language, the tau function is the square root of a determinant in KP variables (Leur et al., 2014, Leur, 2021): $$\tau_{\mathrm{KP}}(\text{odd times}) = \left(\tau_{\mathrm{BKP}}(\text{odd times})\right)^2.$$

3. Relation to Two-Component KP/Toda Hierarchies

2-BKP tau functions are canonically the square root of corresponding 2-component KP (2-Toda) tau functions. Under the fermionic bosonization, the embedding arises via splitting charged fermions into two neutral sectors, and the group element $g$ generating the KP tau becomes a product of $h$ and its involutive image in the Clifford algebra: $$(\tau^{2\text{–BKP}}(t))^2 = \tau^{2\text{–KP}}(t^{(1)} = t,\ t^{(2)} = -t)\quad \text{with all even times set to zero}.$$ This structural identity underpins the Pfaff = $\sqrt{\det}$ phenomenon for ensemble partition functions (e.g. orthogonal, symplectic random matrices) (Lee, 2018, Leur et al., 2014).

4. Matrix Models and Partition Function Formulas

2-BKP tau functions arise naturally as grand partition functions of two-matrix integrals, such as modifications of the Zinn-Justin–Zuber model, providing exactly solvable cases: $$J_N(p^{(1)}, p^{(2)}) = \int_{U_N \times U_N} e^{\mathrm{Tr}(U_1 T U_1^{-1} U_2)}\, \exp\left[\sum_{m\,\mathrm{odd}} \left(p^{(1)}_m\, \mathrm{Tr}\,U_1^m + p^{(2)}_m\, \mathrm{Tr}\,U_2^m \right)\right]\, d\mu_1(U_1)\, d\mu_2(U_2),$$ which expands as a sum over doubled strict partitions in Schur Q-functions. In various limits, these formulas produce explicit multi-integral and Pfaffian solutions linked to classical ensembles and their deformations (Antonov et al., 2023, Orlov et al., 2016, Orlov et al., 2012).

5. Algebraic Structure: Grassmannians, Lax Operators, and Reductions

The geometric underpinning is the Sato or Shiota isotropic Grassmannian, where a maximal isotropic subspace $$U \subset V = \mathbb{C}((z_1^{-1}))e_1 \oplus \mathbb{C}((z_2^{-1}))e_2$$ is associated to each tau-function, and the tau itself is reconstructed via a Plücker coordinate expansion. The Lax formalism gives two scalar Lax operators $L_1$, $L_2$ coupled by a pseudo-differential $H$, whose evolution is encoded in the time variables. Reductions by degree $(M_1, M_2)$ correspond to constrained hierarchies (e.g., Novikov–Veselov, Kac–Wakimoto type D) and fix the form of quadratic and Virasoro constraints (Chen et al., 19 Nov 2025, Cheng et al., 2018).

6. Pfaffian Addition Formulae and Darboux Transformations

The addition (Fay) formulae for 2-BKP are expressed by recursive Pfaffian relations generating the entire hierarchy: $$\frac{\tau\big(\mathbf t - 2\sum_{i=1}^{N_1}[\lambda_i^{-1}]_1 - 2\sum_{j=1}^{N_2}[\mu_j^{-1}]_2\big)}{\tau(\mathbf t)}\, \prod_{i < k}\frac{\lambda_i - \lambda_k}{\lambda_i + \lambda_k} \prod_{j < \ell} \frac{\mu_\ell - \mu_j}{\mu_\ell + \mu_j} = \mathrm{Pf}\, K$$ where $K$ is an explicit block matrix kernel comprised of shifted tau functions. Iterated Darboux transformations provide systematic methods for constructing multi-soliton tau functions and generating Virasoro-type symmetries (Chen et al., 19 Nov 2025).

7. Applications to Geometry, Combinatorics, and Connected Correlators

The 2-BKP tau functions encode enumerative invariants in algebraic geometry, notably spin Hurwitz numbers and the descendant potential of type $D$ singularities (miniversal unfoldings). Connected $(n,m)$-point bosonic correlators are recovered via cycle-sum formulas involving the tau function's affine coordinate kernel, yielding closed expressions for higher-genus Gromov–Witten invariants, generating functions for random partitions, and combinatorial models for symmetric functions. In matrix models, the partition function structure is governed by 2-BKP Pfaffians; for spin double Hurwitz numbers, the associated generating functions are explicit Pfaffian tau functions (Wang et al., 2022, Lee, 2018, Wang et al., 2022).

---

In summary, 2-BKP tau functions provide an algebraic, combinatorial, and geometric framework for integrable hierarchies underlying a spectrum of phenomena from mathematical physics to algebraic geometry, with rich Pfaffian structures, deep connections to 2-KP/Toda determinants, and explicit realizations in matrix models, symmetric function theory, and enumerative geometry (Antonov et al., 2023, Leur et al., 2014, Leur, 2021, Chen et al., 19 Nov 2025, Cheng et al., 2018, Wang et al., 2022, Orlov et al., 2016, Orlov et al., 2012, Wang et al., 2022).