Genus Two Pe-Functions in Hyperelliptic Curves

• Genus two Pe-functions are multivariable Abelian functions defined as logarithmic derivatives of the sigma function on hyperelliptic curves.
• They are constructed using derivatives of theta functions with half-integer characteristics, yielding explicit addition laws and inversion formulas.
• Their intricate differential relations, including trilinear and quartic identities, link period matrices and branch point data, with applications in integrable systems and modular forms.

A genus two Pe-function is a multivariable Abelian function arising as a logarithmic derivative of the genus-two sigma function, itself expressed in terms of the genus-two Riemann theta function with (possibly even) half-integer characteristic. These functions generalize the classical Weierstrass ℘-function to the context of hyperelliptic curves of genus two, and provide the algebraic and analytic toolkit underpinning addition laws, inversion theorems, and explicit formulas for solutions to integrable systems, correlation functions, and modular forms in the genus 2 setting. Their structure is controlled by period matrices, theta constants, and branch point data, and their differential-algebraic relations (including trilinear and quartic identities) encode fundamental properties of the associated complex curves and their Jacobians.

1. Analytic and Algebraic Structure

On a compact Riemann surface $X$ of genus $g=2$, a canonical basis of holomorphic differentials $(\omega_1, \omega_2)$ defines the period matrix $\Omega$: $\Omega_{jk} = \int_{b_j} \omega_k, \quad j,k=1,2,$ with the normalization $\int_{a_j}\omega_k = 2\pi i\,\delta_{jk}$. The genus-two Riemann theta function is

$$\theta(z\mid\Omega) = \sum_{n\in\mathbb{Z}^2} \exp\left\{ \pi i\, n^T \Omega\, n + 2\pi i\, n^T z \right\},$$

where $z\in\mathbb{C}^2$ and $\Omega$ lies in the Siegel upper half-space (symmetric, positive-definite imaginary part) (Tutiya, 2017, Tsuchiya, 13 Jan 2026, Korotkin et al., 2012). One extends this to theta functions with half-integer characteristic $[\alpha;\beta]$,

$$\theta\bigl[{\alpha\atop \beta}\bigr](z|\Omega) = \sum_{n\in\mathbb{Z}^2} \exp\Bigl\{ i\pi (n+\alpha)^T \Omega (n+\alpha) + 2\pi i (n+\alpha)^T (z+\beta) \Bigr\},$$

with $\alpha,\beta\in(1/2)\mathbb{Z}^2$.

On top of this, the genus-two sigma function can be built via

$$\sigma(u) = C \exp\left( \frac{1}{2} u^T \eta \omega^{-1} u \right) \theta[\beta](2 \omega^{-1} u \mid \tau),$$

where $u$ is a vector of Abel-Jacobi coordinates, $\omega$ and $\eta$ are period matrices, and $\beta$ selects the characteristic (Korotkin et al., 2012, Nakayashiki, 2015). The Pe-functions of genus two are defined as derivatives of the logarithm of sigma: $$P_{IJ}(a) = -\partial_I \partial_J \log \sigma(a), \quad P_{IJK}(a) = -\partial_I \partial_J \partial_K \log \sigma(a), \dots$$ These objects are multivariable analogues of the classical Weierstrass $\wp$-function and its derivatives.

2. Pe-Functions: Definition, Properties, and Half-Periods

The Pe-functions $P_{IJ}(a)$ (and higher $P_{I_1\ldots I_k}(a)$) depend on the period matrix and the point $a\in\mathbb{C}^2$ modulo the period lattice. For specific arguments such as non-singular even half-periods—labeled $a_{ij}$, constructed from the images under the Abel map of branch points plus the Riemann constant (Tsuchiya, 13 Jan 2026)—the Pe-functions evaluate to central values governing addition formulas and identities for genus-two objects.

On a genus-2 hyperelliptic curve $y^2 = \prod_{k=1}^5(x-e_k)$, the Abel map images of branch points $Q_j$ provide half-periods, and their combinations $a_{ij}=Q_i+Q_j+R$ (with $R$ the Riemann constant) are even and non-singular: $$\theta\bigl[{\alpha_{ij}\atop\beta_{ij}}\bigr](0|\Omega) \neq 0.$$

3. Differential Equations, Trilinear Relations, and Kleinian Decomposition

The genus-two Pe-functions satisfy a system of 15 fundamental PDEs—ten cubic and five quartic equations—encoding their algebraic structure: $$P_{IJK} P_{LMN} = \sum_{(G,H),(R,S)} C^{IJK,LMN}_{GH,RS} P_{GH} P_{RS} + \text{lower order terms},$$

$$P_{IJKL} = \sum_{(E,F),(G,H)} D^{IJKL}_{EF,GH} P_{EF} P_{GH} + \dots$$

These relations enable reduction of higher powers to at most quadratic in the Pe-functions, and when evaluated at half-periods they yield trilinear relations sufficient to expand arbitrary symmetric products in a finite basis (Tsuchiya, 13 Jan 2026).

Any product of Szegő kernels (cyclic product of genus-2 fermion correlation functions) can thus be expressed as

$$\prod_{i=1}^N S_\delta(z_i, z_{i+1}) = \left[ \prod_{i=1}^N \omega_1(z_i) \right] \sum_{I_1, \dots, I_k} H_{I_1 \dots I_k}(B) P_{I_1 \dots I_k}(a_{ij}),$$

where the $H_{I_1 \dots I_k}$ are meromorphic in the relative Abelian integrals $B_i$ (Tsuchiya, 13 Jan 2026). Ultimately, using Kleinian theory, such Abelian functions admit a decomposition in terms of the $P_{IJ\cdots}$ and their algebraic relations.

4. Connection to Theta-, Sigma-, and Addition Formulas

All genus-two Pe-functions are logarithmic derivatives of the sigma function, itself a composition of the Riemann theta function and an exponential of a certain bilinear form (Korotkin et al., 2012): $$P_{IJ}(a) = \frac{\partial_I \partial_J \theta_R(a)}{ \theta_R(a) } - \frac{ \partial_I \theta_R(a) \partial_J \theta_R(a) }{ \theta_R(a)^2 }.$$ This allows rewriting products or relations originally posed in terms of Pe-functions purely as algebraic combinations of (shifted) derivatives of the genus-two theta function at the relevant half-periods.

Genus-two addition formulas and shift identities generalize the elliptic addition formulas, e.g.,

$$\sigma(u+v)\sigma(u-v)\sigma(0)^2 = \sum_{I\leq J} [P_{IJ}(u)P_{IJ}(v)] \prod_{K<M}(u_K - v_K)(u_M+v_M),$$

giving a purely Abelian function-theoretic expansion in terms of Pe-functions (Tsuchiya, 13 Jan 2026).

5. Thomae-Type Formulas and Branch Point Expressions

The relation of theta constants and their derivatives at half-periods to branch point data is encoded by Thomae-type formulas for genus two (Bernatska, 2019). The even non-singular theta constants at zero,

$$\theta\bigl[\varepsilon(I_0)\bigr](0;\tau) = \epsilon \Bigl( \frac{\det\omega}{\pi^2} \Bigr)^{1/2} \Delta(I_0)^{1/4} \Delta(J_0)^{1/4}$$

(recall $I_0$, $J_0$ partition the six branch points), and derivatives at odd singular half-periods are given by

$$\partial_{v_n} \theta[\varepsilon^{(k)}](0;\tau) = \epsilon \Bigl( \frac{\det\omega}{\pi^2} \Bigr)^{1/2} \Delta(\{k\})^{1/4} \Delta(J_1)^{1/4} (\omega_{1n} - e_k \omega_{2n}),$$

linking algebraic geometry of the hyperelliptic model to the behavior of genus-two Pe- and theta-functions. These relations, together with the classical Bolza formula, allow for recovery of explicit branch points from ratios of theta-derivatives (Bernatska, 2019).

6. Applications: Integrable Systems, Fermion Correlators, and Moduli

Genus-two Pe-functions appear naturally:

• In the construction of explicit genus-two algebro-geometric solutions (e.g., the ILW equation (Tutiya, 2017)) via Baker-Akhiezer functions and the machinery of the Krichever/Dubrovin school;
• In decomposition of cyclic fermion correlators, where the spin structure dependence is expressible entirely in terms of Pe-function values at non-singular even half-periods, enabling algebraic summation over spin structures in superstring amplitudes (Tsuchiya, 13 Jan 2026);
• Through their trilinear and quartic relations, which control reduction of multivariable Abelian products in moduli, and their replacement by shifted theta-derivatives enables rewriting in terms of a single, generically shifted, genus-two theta function.

A crucial aspect is the ability to relate all such Abelian function-theoretic expressions to rational functions in the $(x,y)$-coordinates of points on the underlying curve, utilizing the Jacobi inversion theorem and its modifications for difference arguments (Shigemoto, 2016, Tsuchiya, 13 Jan 2026).

7. Summary Table: Pe-Functions, Sigma Functions, and Theta Expressions

Object	Definition / Construction	Genus-Two Specialization
Sigma function	Combination: exponential quadratic × theta with characteristic	$$\sigma(u) = C \exp(\frac{1}{2} u^T \eta \omega^{-1} u) \theta[\beta](2\omega^{-1}u|\tau)$$
Pe-function	Multiderivatives of $\log\sigma(u)$	$P_{IJ}(a) = -\partial_I \partial_J \log \sigma(a)$
Theta function	Riemann theta function (possibly with characteristic)	$$\theta\bigl[{\alpha\atop\beta}\bigr](z|\Omega) = \sum_{n\in\mathbb{Z}^2} \ldots$$
Abelian relations	Trilinear/quartic identities among Pe-functions	E.g., $$P_{222}^2 = 4P_{11} + 4 e_3 P_{12} + 4 e_2 P_{22} + 4 e_1 P_{22} + \cdots$$
Thomae/Bolza	Evaluation at (derivatives at) half-periods in branch point data	See explicit branch point formulas above, relating theta/Pe values/derivatives to elementary symmetric polynomials

These objects constitute the foundation of explicit Abelian function theory, theta-functional addition laws, and higher-genus extensions of elliptic function identities, providing the algebraic scaffolding for both theoretical and applied considerations in mathematical physics and algebraic geometry (Korotkin et al., 2012, Nakayashiki, 2015, Tsuchiya, 13 Jan 2026).