Topological recursion and uncoupled BPS structures II: Voros symbols and the $\tau$-function

We continue our study of the correspondence between BPS structures and topological recursion in the uncoupled case, this time from the viewpoint of quantum curves. For spectral curves of hypergeometric type, we show the Borel-resummed Voros symbols of the corresponding quantum curves solve Bridgeland's"BPS Riemann-Hilbert problem". In particular, they satisfy the required jump property in agreement with the generalized definition of BPS indices $\Omega$ in our previous work. Furthermore, we observe the Voros coefficients define a closed one-form on the parameter space, and show that (log of) Bridgeland's $\tau$-function encoding the solution is none other than the corresponding potential, up to a constant. When the quantization parameter is set to a special value, this agrees with the Borel sum of the topological recursion partition function $Z_{\rm TR}$, up to a simple factor.


Introduction
This paper continues our study [IK] of the relationship between the theory of BPS structures and the formalism of topological recursion, in the special case where the former is uncoupled. In the present work, we approach this based on the relationship that each of these two theories has with exact WKB analysis -in particular, the theory of quantum curves. As before, our testing ground will be concrete examples related to the Gauss hypergeometric equation and its confluent degenerations, the spectral curves of hypergeometric type, summarized in Table 1 below. In our previous work, we gave a formula for the topological recursion free energies F g as a sum over BPS cycles, weighted by BPS indices (γ ) whose definition we introduced generalizing [GMN2,BrS]. In the present paper we build on this result (which is formal in ), following the topological recursion construction of the quantum curve in [IKoT1,IKoT2] to solve Bridgeland's BPS Riemann-Hilbert problem [Br1] (which is analytic in nature). Our solution is given by the Borel-resummed Voros symbols of Table 1. Spectral curve • : y 2 − Q • (x) = 0 of hypergeometric type, where • ∈ {HG, dHG, Kum, Leg, Bes, Whi, Web, dBes, Ai} labels the equation Bessel (Bes) x + 4m 2 0 4x 2 m 0 = 0.
the quantum curve, and we furthermore obtain the corresponding "τ -function", together with its interpretation in terms of a potential for the Voros coefficients. We will state our main result in detail below after giving a brief background of these two theories, together with the results of our previous paper [IK]. Bridgeland [Br1], are an axiomatization of the structure of Donaldson-Thomas (DT) invariants of a Calabi-Yau 3 triangulated category with a stability condition. It consists of a triple ( , Z , ) of a "charge lattice" with a skew-symmetric pairing, a "central charge" homomorphism Z : → C, and a collection of numbers ("BPS indices") : → Q, associated to each γ ∈ , satisfying some conditions. This is a special case of the stability structure in [KS1].

BPS structures. BPS structures, introduced by
Given a BPS structure, Bridgeland also associated a Riemann-Hilbert problem on C * = C * with the coordinate , which we call the BPS Riemann-Hilbert problem, seeking a collection of piecewise holomorphic (or meromorphic; see §1.4.2 below and [Ba]) twisted-torus-valued functions X ( ), labeled by a ray ⊂ C * , with a collection of BPS rays on which the functions have discontinuity. More precisely, X is defined on a half-plane H ⊂ C * whose bisecting ray is given by , and when we vary so that it crosses over a BPS ray BPS , then X has a discontinuity described by a cluster-like birational transformation, called the BPS automorphism, encoded by the BPS invariants (γ ) of γ ∈ satisfying Z (γ ) ∈ BPS . The BPS Riemann-Hilbert problem is closely related to the one studied by Gaiotto-Moore-Neitzke in [GMN1,GMN2], in which the BPS spectra for "class S theories" were investigated. These Riemann-Hilbert problems allows us to understand the Kontsevich-Soibelman wall-crossing formula, a non trivial relation among the BPS invariants under variations of stability conditions, as a kind of iso-Stokes property around the "irregular singular point" = 0 (see also [BrTL,FFS]). This wall-crossing structure is axiomatized in the notion of variation of BPS structures.
If a variation of BPS structures over a complex manifold M is nicely parametrized ("miniversal"), the solution to the associated BPS Riemann-Hilbert problem is naturally regarded as a collection of piecewise-holomorphic functions on M × C * , and can be used to define a certain -dependent vector field on M. In this setting, Bridgeland also defined a function on M × C * , which we call the BPS τ -function and denote by τ BPS , by the condition that log τ BPS be a Hamiltonian for this vector field (see §2.5).
If the BPS structure is uncoupled, it was shown in [Br1] that the solutions to the BPS Riemann-Hilbert problem can be solved explicitly in terms of gamma functions, and the associated BPS τ -function is expressed in terms of a product of Barnes G-functions. Although non-trivial wall-crossing does not occur in the uncoupled situation, it was observed in [Br1,Br2] that the BPS τ -function has an alternative interesting interpretation. Namely, τ BPS for the BPS Riemann-Hilbert problem associated with the natural variation of (uncoupled) BPS structures arising from resolved conifold through Donaldson-Thomas theory provides a "non-perturbative partition function" for the resolved conifold; that is, the asymptotic expansion of τ BPS when → 0 coincides with the generating series of all-genus Gromov-Witten invariants. This surprising discovery was also generalized in [Sto].
The study of BPS Riemann-Hilbert problems associated to spectral curves of hypergeometric type through our previous work [IK] -in particular, the construction of their solutions and τ -functions via topological recursion -will be the main theme of our present work. We will recall how the BPS structures arise from such spectral curves in §1.3.

Topological recursion and spectral curves of hypergeometric type.
Topological recursion (TR) is an algorithm generalizing the construction of correlation functions and partition functions in matrix models, whose input data is an algebraic curve enhanced with some additional data, called a spectral curve ( [EO1,CEO]). Mathematically, to a given spectral curve, TR constructs a doubly-indexed sequence {W g,n } g≥0, n≥1 of meromorphic multi-differentials on the spectral curve, and a sequence {F g } g≥0 of complex numbers.
These outputs of TR are expected to encode information of enumerative geometry such as Gromov-Witten invariants, Hurwitz numbers, etc. (essentially due to various dualities in mathematical physics), and there is a large class of spectral curves for which this expectation is known to hold true (see [EO2,E1,DMNPS,DN] etc.). As represented by the celebrated Kontsevich-Witten theorem [Kon,W], generating series of such enumerative invariants are expected to be related to τ -functions of some corresponding integrable systems. In fact, at the level of formal power series of , KdV and Painlevé τ -functions have already been constructed via (the discrete Fourier transform of) the TR partition function Z TR = exp( g≥0 2 g−2 F g ) for certain spectral curves; see [EO1,IS,IM,IMS,MO1,I2,EG,MO2] for example. The actual τ -function should be given as the Borel sum of those formal series, and thus, in this paper, we denote by τ TR the Borel sum of the TR partition function. Although the Borel summability of TR partition functions is not proved in full generality (see [E2] for the growth estimate of the coefficients), those which appear in this paper are Borel summable along all but finitely many rays.
In [IKoT1,IKoT2], the first named author together with collaborators applied the TR formalism to the spectral curves of hypergeometric type, and obtained a result on quantum curves which states the following: after the principal specialization (i.e., setting z 1 = · · · = z n = z), a certain generating function of a primitive F g,n of W g,n (i.e., a function satisfying 1 d z 1 · · · d z n F g,n = W g,n ) gives a WKB solution to a Schrödinger-type ODE, which we denote by E, whose classical limit recovers the spectral curve. Indeed, we may add a tuple of parameters ν which parametrizes the integration contour to define the primitive F g,n , and hence, we obtain a family E = E(ν) of quantum curves. We note that they also depend on the tuple m = (m s ) of parameters contained in the definition of Q(x), which is one of Q • (x) in Table 1. We also observe that the resulting quantum curve is equivalent to an ODE satisfied by classical special functions; for example, the quantum curve of HG is equivalent to the Gauss hypergeometric differential equation. This allows us to study the TR invariants in terms of various properties of special functions. As one of the main results in [IKoT1,IKoT2], an explicit formula in terms of Bernoulli numbers was given for the F g , as well as the Voros coefficients of the quantum curve E(ν), defined as a term-wise integral of the formal series given as the logarithmic derivative of the WKB solution along cycles on the spectral curve. It turns out that the Borel sum of generating series of free energy and Voros coefficients can be interpreted through the language of BPS structures.

Brief review of results of Part I.
To formulate the results of [IK], let us first recall how BPS structures are constructed from spectral curves of hypergeometric type. There is a canonical way to construct (a variation of) BPS structures from such (a family of) spectral curves (see [Br1,§ 7] for the general case, and see also [GMN1] as the origin of this construction). For a given spectral curve = • in Table 1, we take the charge lattice and pairing to be = {γ ∈ H 1 ( , Z) | ι * γ = −γ } and the intersection pairing on it, respectively. Here, is a partial compactification of (see §3.1), and ι is the covering involution ι : → . The central charge is given by the period map Z (γ ) = γ ydx = γ √ Q(x) dx, and the BPS indices (γ ) were defined by a certain weighted counting of degenerations of spectral networks (Stokes graphs) for the meromorphic quadratic differential Q(x) dx 2 . We recall these details more precisely in §3.3.3. The list of BPS structures which we obtained in [IK] following this construction is summarized concisely in Table 3. In all our examples, since the curves in Table 1 are of genus 0, has the trivial intersection pairing so that the resulting BPS structure is a fortiori uncoupled. The rank of is equal to the number of mass parameters m = (m s ) attached to even order poles s of Q(x). If we vary m, we have a variation of BPS structure parametrized by m ∈ M, where M is the space of mass parameters satisfying the assumptions in Table 1. We can describe the central charge Z (γ ) explicitly as a function of m, and it is easy to check that the variation is miniversal.
We recall that our definition of (γ ) generalizes the one given in [GMN1,BrS,Br1] in order to include degenerations involving a saddle connection with endpoint a simple pole of Q(x) dx 2 , and degenerate ring domains (see [IK,§ 3.6]). With this generalization, we have shown in [IK] that the TR free energy F g for the spectral curve is described as a sum over all BPS cycles whose central charge lie in a generically chosen half-plane ( [IK,Theorem 5.3]). The following formula for TR free energy with g ≥ 2 is one of the main results in [IK]: (1.1) Here, H is any half-plane whose boundary rays are not BPS. This universal formula is valid for all our examples, and shows that each term in the free energy is a contribution of a BPS cycle associated with a degenerate spectral network; furthermore the coefficient of each term knows the BPS index of the BPS cycle. We conjectured that the formula (1.1) is still valid for higher degree spectral curves if the associated BPS structure is uncoupled, and checked this expectation numerically for a few examples of degree 3 spectral curves ( [IK,§ 5.3]). However, in our previous work, we did not give any explanation for the "naturality" of our generalized definition of (γ ); namely, we defined them so that the universal formula (1.1) is valid. Therefore, in the present paper, we will try to fill this gap 2 and show the validity of our generalization of (γ ) by showing they control the analytic behaviour of (Borel-resummed) Voros coefficients. We will make this more precise below.
1.4. Results of Part II. The BPS indices (γ ) are defined by weighted counting of degenerate spectral networks, but they have another interpretation in Gaiotto-Moore-Neitzke [GMN1]. That is, the BPS indices should appear in the exponents in the formula which describe the mutation of the Fock-Goncharov cluster coordinates (on the moduli space of flat SL 2 (C)-connections) caused by the degeneration of the spectral network. On the other hand, as is shown in [I1,All1,Ku], the Fock-Goncharov cluster coordinates can be identified (on the oper locus, where this makes sense) with the Borel-resummed Voros symbols (exponential of Voros coefficients) in the exact WKB analysis of a certain Schrödinger-type equation with the spectral curve as its classical limit (see also [IN1,IN2]), where the mutation formula (BPS automorphism) is a consequence of the Stokes phenomenon for the Voros symbols.
Keeping these previous works in mind, to investigate the relation between BPS structures and TR, we study the Stokes structures of the Voros symbols for the quantum curves constructed through TR from the spectral curves of hypergeometric type obtained in [IKoT1,IKoT2], and compare those to the BPS Riemann-Hilbert problems associated with the BPS structure studied in our previous work [IK]. We note here that such analysis of the Stokes phenomenon for the Voros symbols of the (confluent) hypergeometric differential equations are not new; some special cases (e.g., for specific ν) have already been studied in the exact WKB literature (see [AT,ATT1,AIT,T,KoT,KKKoT,KKT] for example). We also note that Stokes structures of the Voros symbols for Schrödinger-type equations were also studied in [DDP1,IN1] in a general setting.
Let us summarize our results more precisely below.

Voros solution to the (almost-doubled) BPS Riemann-Hilbert problem.
A natural BPS Riemann-Hilbert problem, which is compared to the Stokes structure on the TR side, is closely related to the one associated with the doubled BPS structure which has a charge lattice D = ⊕ ∨ . In what follows, we identify a relative homology class in H 1 ( , D ∞ , Z), where D ∞ is the divisor supported on the poles of Q(x)dx 2 of order ≥ 2, with an element in ∨ through the Poincaré-Lefschetz duality. For the quantum curves E = E(ν) constructed in [IKoT1,IKoT2], first we show that the Voros coefficients for relative homology classes ("path Voros coefficients") have an explicit expression as a sum over all BPS cycles (which is analogous to the formula (1.1) for the TR free energy): for a given β ∈ ∨ satisfying ι * β = −β, we have with the coefficients being given in terms of the Bernoulli polynomials: (see §4.5 for the precise notation). We also have an explicit expression of Voros coefficients V γ for homology classes γ ∈ ("cycle Voros coefficients") in terms of the parameters m and ν. Together, we find that • The Voros symbols e V γ , e V β are Borel summable as a formal series of along any ray which is not the BPS rays for the corresponding BPS structure given in [IK], and • The Borel-resummed Voros symbols jumps when crosses a BPS ray due to the Stokes phenomenon, and the action of Stokes automorphism on the Voros symbols are given in the exactly same manner as the BPS automorphism (see Proposition 5.2) with the BPS indices (γ ) defined in [IK]. This justifies our definition of (γ ).
The second observation tells us that the Borel-resummed Voros symbols of the quantum curve gives a certain solution to the BPS Riemann-Hilbert problem. To describe the underlying BPS structure, let us introduce the almost-doubled BPS structure ( D , Z D , D ), which has the charge lattice D = ⊕ * where * = {β ∈ ∨ | ι * β = −β}, (1.3) with trivially extended central charge Z D and BPS indices D (we note that ( D , Z D , D ) is also uncoupled). The constant term ξ in the small asymptotics of the solution is specified by the aforementioned parameter ν together with a specific quadratic refinement (which introduces an appropriate sign to give a twisted-torus-valued function). The reason for considering the sublattice * is essentially to avoid a square root type singularity in the Borel-resummed Voros symbols (i.e., so we can obtain a meromorphic solution to the BPS Riemann-Hilbert problem), although the reader can take the full lattice if they are happy to accept such a singularity.
We call the above solution the Voros solution to the almost-doubled BPS Riemann-Hilbert problem, and denote it by X Vor . We give its explicit expression in (5.7). [All2,Br1,Ba] 3 . For the special choice of ξ where ξ(γ ) = 1 if (γ ) = 0, Bridgeland gave the explicit general piecewise-holomorphic solution, which we will denote by X hol for all uncoupled BPS structures in [Br1]. This was generalized by [Ba] who noted that, even in the uncoupled case, holomorphic solutions do not usually exist for other values of ξ , but solved the problem meromorphically in this case for arbitrary ξ . [All2] provides a (meromorphic) solution for the general (coupled) set up of a quadratic differential with essentially arbitrary pole configuration, but at a specific value of ξ , and in terms of Fock-Goncharov coordinates. Finally, [Br1] solves the problem for arbitrary ξ in the simplest nontrivial coupled case of the cubic oscillator, also using Fock-Goncharov coordinates.

Relation between the Voros solution and other solutions. Similar problems have been considered in
Among these, our work is closely related to the problem considered by Barbieri [Ba]. Like her work, we provide a solution for arbitrary ξ , at the cost of considering only uncoupled examples. We note that our solution is not the so-called minimal solution X min in the sense of [Ba], which is reduced to the holomorphic solution in the original case ξ(γ ) = 1 considered by Bridgeland [Br1], but in §5.2 we show that the minimal solution and Voros solution differ only by a simple factor and provide an explicit relation between the two. We will make this more precise in §5.2.

Voros potential, BPS τ -function and TR partition function function
We observe that Voros coefficients of the quantum curve E(ν) naturally define a formal power series valued one-form on the space M of parameters m which is closed (in fact exact), and arises from a potential function φ; in particular, we have a formal power series valued function on M satisfying where {β s } is a natural basis of * . We call φ the Voros potential. The following theorem tell us that the Voros potential (more precisely, exponential of its -derivative) is a topological recursion / quantum curve analogue of the BPS τ -function.
Theorem 1.2 (Theorem 5.10). Let ( D , Z D , D ) be the almost-doubled variation of BPS structures corresponding to one of the families of spectral curves of hypergeometric type, and be any non-BPS ray. The BPS τ -function τ Vor BPS, associated to the Voros solution X Vor is given by the Borel sum of the -derivative of the Voros potential as: where c is an explicit constant factor (see Theorem 5.10).
We can also express τ Vor BPS, as a product of modified Barnes G-functions, which is closely related to the BPS τ -function obtained by Alexandrov-Pioline [AP1].
In all examples, we furthermore find there exists a choice (not unique) of the parameter ν = ν * so that we have ξ D,ν * (γ BPS ) = 1 for all BPS cycles γ BPS . If we specialize the parameter to this value, we have Corollary 1.3 Let ( D , Z D , D ). be the almost-doubled variation of BPS structures corresponding to one of the families of spectral curves of hypergeometric type, and be any non-BPS ray. If we specialize the parameter value ν = ν * , the BPS τ -function

6)
where κ is either 1 or a simple explicit function independent of .
This BPS τ -function τ hol BPS, was originally obtained by Bridgeland in [Br1], where its explicit expression in terms of the Barnes G-function was given. Through a well-known formula for the asymptotic expansion of the Barnes G-function and the formula (1.1) of our previous work [IK], we can relate τ hol BPS to the Borel-resummed topological recursion partiton function τ TR . More precisely, we have: (1.7) Here τ >0 TR, denotes the Borel-resummed TR partition function without the genus 0 part. These results establish a precise link between objects arising from the topological recursion (quantum curves and their Voros coefficients/potential), and objects arising from the BPS structures (solution to the BPS RHP, τ -function, etc.) coming from the same initial data; i.e., the spectral curves of hypergeometric type. We summarize our results as a dictionary between uncoupled BPS structures and topological recursion, depicted in Table 2 below, which holds precisely in all our examples (we believe these examples to be exhaustive of finite uncoupled cases when the spectral curve is of degree 2).
Although the TR is applicable to higher degree spectral curves ( [BHLMR,BoE1,BoE2]), the existence of the associated BPS structure remains conjectural: we expect that Gaiotto-Moore-Neitzke's approach in [GMN3] allows us to define a natural BPS structure for higher degree spectral curves. We will briefly verify our picture continues to hold in this case based on numerical experiments in our previous work, following computations by Y. M. Takei [Ta].

Remarks.
(1) Based on our result, we would like to emphasize the following point. For the special parameter value ν = ν * we have shown that the Borel sum of the (g > 0) TR partition function is equal to τ hol BPS , so that in particular τ TR encodes the solution to the BPS Riemann-Hilbert problem. That is, τ TR is characterized by its discontinuity structure and asymptotics, so that one may approach the partition function of TR from the perspective of Riemann-Hilbert-type methods. While the correspondence is literal for all examples we have considered, we hope that this perspective may be useful to deepen our understanding of the topological recursion in more general cases.
(2) We also note a similar structure to our result appears in the solution to a different (related) problem, but in a certain limit of quadratic differentials with four secondorder poles ("nodal limit of the four-punctured sphere"). In particular, we note the resemblance to the expression for the 1-loop part of the twisted effective superpotential 4 W eff in the conjecture of Nekrasov-Rosly-Shatashvili [NRS,TV,HK,JN], identified with the generating function of ( -)opers. There, a hypergeometric spectral curve also appears, essentially as the four-punctured sphere breaks up into two three-punctured spheres. However, that calculation corresponds to the 1 → 0 limit, whereas the approach of this paper corresponds to the 2 = − 1 = limit. It is somewhat puzzling that these two regimes produce a similar structure. (3) While τ BPS differs slightly from our τ TR , the expression for τ TR in the Weber case (a single BPS state) coincides exactly with the special function γ appearing in the work of Nekrasov-Okounkov [NO] as the perturbative part of the Nekrasov partition function. While the difference may be attributed to possible variants of the definition of τ BPS and choices of renormalization schemes defining the partition function, it would be desirable to understand this relationship better. (4) Finally, we mention that our results can also be interpreted in the context of socalled Joyce structures [Br3]. This is a recently proposed geometric structure arising naturally on the space of stability conditions from the theory of wall-crossing and Donaldson-Thomas invariants (see also [BrSt,AP1,AP2,AST] for the associated hyperkähler structure). Our solution to the BPS Riemann-Hilbert problem provides a canonical example of one such geometry, though we leave its detailed computation and further investigation to future work.
1.6. Organization. The paper is structured as follows. In §2 we recall the notion of BPS structures and the define the associated BPS Riemann-Hilbert problem and τ -function. In §3, we focus on the BPS structures arising from spectral curves of hypergeometric type, recall our previous results, and prepare for the proof of the main results. In §4 we recall the basics of topological recursion, and the relevant results on quantum curves and associated the Voros coefficients, describe them in terms of BPS structures and express their Stokes phenomena. Finally, in §5, we assemble the pieces and solve the BPS Riemann-Hilbert problem, relate τ BPS to the TR side, and remark on evidence in higher rank examples.

Review of BPS Structures and the BPS Riemann-Hilbert Problem
Here we recall the basic facts on BPS structures, including the results of our previous work [IK], define the BPS Riemann-Hilbert problem [Br1], and give some preparatory lemmas about the variations of BPS structures we will work with. • a free abelian group of finite rank equipped with an antisymmetric pairing ·, · : × → Z, (2.1) • a homomorphism of abelian groups Z : → C, and • a map : → Q, satisfying the conditions • Support property: for some (equivalently, any) choice of norm · on ⊗ R, there is some C > 0 such that We call the charge lattice, and the homomorphism Z is called the central charge. The rational numbers (γ ) are called the BPS indices or BPS invariants.
Let us recall some useful terminology [Br1]: • An active class or BPS state is an element γ ∈ which has nonzero BPS index, (γ ) = 0. • A BPS ray is a subset of C * which can be written as γ = Z (γ ) R >0 for some active class γ .
We note that the central charge Z (γ ) for an active class γ is nonzero due to the support property. We can also consider certain classes of BPS structures with especially nice properties. A BPS structure ( , Z , ) is said to be • finite if there are only finitely many active classes, • uncoupled if γ 1 , γ 2 = 0 holds whenever γ 1 , γ 2 are both active classes • integral if the BPS invariant takes values in Z.
Though we will stick to such examples in this paper, in order to formulate the BPS Riemann-Hilbert problem in the next section, we may weaken these conditions. We call a BPS structure ray-finite if there are finitely many active classes γ with Z (γ ) ∈ for a given BPS ray , and generic if γ 1 , γ 2 = 0 whenever γ 1 , γ 2 are active and Z (γ 1 ) and Z (γ 2 ) lie on the same BPS ray.
If the BPS structure is not uncoupled, we call it coupled. This is in fact the more typical situation, but the analysis of the BPS structure and the corresponding Riemann-Hilbert problem given in the next section is much more difficult in this case (see [Br1,All2]). All our considerations in this paper will be for finite, uncoupled and integral BPS structures.

Twisted torus and BPS automorphism. For a lattice with an antisymmetric pairing
·, · , we define the corresponding twisted torus T − as the set of twisted homomorphisms into C * : We also have the usual torus T + := Hom( , C * ). We denote by x γ : T − → C * the twisted character naturally associated to γ ∈ , defined by x γ (ξ ) = ξ(γ ), which generate the coordinate ring of T − . Similarly, the ordinary characters of the usual torus are denoted y γ . The twisted and usual characters satisfy Since T + acts on T − in the obvious way, a point ξ ∈ T − can be understood as a point in T + together with a "base point" σ ∈ T − . There is a particularly natural set of possible base points, intertwining usual inversion with the map x γ → x −γ , those taking values in Z 2 . These are the so-called quadratic refinements σ : → {±1}.
Since T + carries a Poisson structure given by We may transport this to T − using any such base point, giving it a Poisson structure too (independent of the choice). The twisted torus T − will essentially be the space in which the solution to an associated Riemann-Hilbert problem takes values. The Riemann-Hilbert jump contour will be a finite union of rays in C * whose coordinate is denoted by . Here, a ray is a subset of the form = e iϑ R >0 ⊂ C * . Given a ray , we associate a BPS automorphism S( ) as the jump factor in the Riemann-Hilbert problem. In [Br1], the BPS automorphisms were defined as a time 1 Hamiltonian flow of a certain function which contains the information of BPS/DT invariants. Note that these operators were also considered in [KS1,GMN1].
In general, defining the BPS automorphism needs a careful analysis of convergence problems, and one may need to work with an appropriate completion of the twisted torus. However, thanks to [Br1,Proposition 4.2], the BPS automorphisms for finite, uncoupled, integral BPS structures are given by the explicit birational automorphisms (2.6) Here we define S( ) via its pullback S( ) * acting on the coordinate ring of T − , also called the Kontsevich-Soibelman transform in [GMN2]. Note that the operator acts trivially on twisted characters x γ corresponding to "null" elements γ ∈ with zero intersection γ, γ for all γ ∈ . We will take (2.6) as the definition of the BPS automorphism, and consider the associated Riemann-Hilbert problem below.

BPS Riemann-Hilbert problem.
Here we formulate the Riemann-Hilbert problem associated with a ray-finite, generic, integral BPS structure, following [Br1]. We will call it the BPS Riemann-Hilbert problem for short. Roughly speaking, in the BPS Riemann-Hilbert problem, we seek a twisted-torus-valued, piecewise-meromorphic map on C * which jumps when crossing a BPS ray, with the jump given precisely by the BPS automorphism. Note that the BPS Riemann-Hilbert problem was defined for more general BPS structures (i.e., for "convergent" BPS structure) in [Br1]; see [Br1,Appendix A,B] for technical details in defining the BPS automorphisms in general setting. Given a ray = e iϑ R >0 ⊂ C * for some ϑ ∈ R, we denote by the half plane centred on . Then, the BPS Riemann-Hilbert problem is formulated as follows (without loss of generality we may consider H as the domain): for every γ ∈ , with the following properties: (RH1) Jumping. Suppose that two non-BPS rays 1 , 2 ⊂ C * form the boundary rays of an acute sector ⊂ C * taken in the clockwise order. Then, whenever ∈ H 1 ∩ H 2 and 0 < | | 1. Here S( ) is the birational automorphism of T − given by the product of all BPS automorphisms S with ⊂ . More explicitly, whenever is a non-BPS ray. (RH3) Polynomial growth at ∞. For every γ ∈ and non-BPS ray , there exists some We denote the collection of all X ,γ by X , which should then be understood as meromorphic maps where denotes the set of zeroes of all X ,γ , for γ ∈ .
Remark 2.3. Note that the BPS Riemann-Hilbert problem is closely related to the infinitedimensional Riemann-Hilbert problem encountered by Gaiotto-Moore-Neitzke [GMN2] in constructing special Darboux coordinate functions on the moduli space of Higgs bundles (see also [BrTL,FFS]). The solution is connected to exact WKB analysis in [Br1,All2], and is expected to furthermore satisfy a kind of thermodynamic Bethe ansatz equation [G,GMN2,N]. Many interesting properties remain to be understood, even conjecturally. Solutions to the BPS Riemann-Hilbert problem should furthermore provide a map between cluster varieties and spaces of stability structures, proved in certain cases by [All2]. Apart from the uncoupled cases and the works [All2,Br1], very little is currently known about the general solution.
For any finite, uncoupled, and integral BPS structure, Bridgeland [Br1] gave the explicit and unique piecewise-holomorphic solution to this problem for a specific choice of constant term (ξ(γ ) = 1 for every active γ ∈ ) in terms of the gamma function. Barbieri observed in [Ba] that for other values of ξ , there will usually be no holomorphic solution, but generalized the solution of [Br1] as follows. To write it down, we use a variant of the gamma function, (2.14) and will later sometimes write (w) = (w, 0) = (w, 1). It is a multivalued function due to the factor w w+η+1/2 = exp((w + η + 1/2) log w), and we will restrict the domain of the function to Re w > 0 and take the principal branch of log w. The asymptotic expansion of (w, η) is given explicitly (see e.g. [Ba]) as valid when w → ∞ in the complement of any closed sector containing R <0 .
Theorem 2.4 ( [Ba,Br1]). Let ( , Z , ) be a finite, uncoupled, integral BPS structure, and fix ξ ∈ T − . Then the associated BPS Riemann-Hilbert problem has a solution, which associates to a non-BPS ray the map X min whose value on γ ∈ is given by 5 Here, we fix the branch of the logarithm appearing in the formula so that Im log ξ(γ ) ∈ [0, 2π) for γ ∈ . When ξ is chosen so that ξ(γ ) = 1 for all active γ ∈ , this becomes the unique solution which is holomorphic on the half planes H , denoted (2.17) Note that the solution has the structure of a product over BPS states, since the exponent will vanish whenever does. The main observations of [Ba] are • There is an explicit solution X min given by (2.16).
• The solution is not unique, but any two solutions are related by multiplication by a meromorphic function f ( ) of satisfying f (0) = 1. • The solution X min is distinguished by the property that all poles and zeroes must be simple, may lie only at Z (γ ) log ξ +2πik , with γ ∈ and k ∈ Z, and for each γ ∈ the points Z (γ ) log ξ +2πik , k ∈ Z lie on the same side of γ . In our case, we will solve (almost) the same Riemann-Hilbert problem using Voros symbols, but our solution will differ from X min . We will use the expression (2.16) to relate the two precisely in §5.2.

Variation of BPS structures.
BPS structures become even more interesting when considered in families. In this case, we require additional axioms ensuring the behaviour is the same as the Donaldson-Thomas invariants, in particular the wall-crossing phenomenon. is covariantly constant as m varies, whenever the boundary rays of are non-BPS throughout.
The precise meaning of condition (VBPS4) for general variation of BPS structure is given in [Br1,Appendix A]. For us, since we will only consider finite, uncoupled, integral BPS structures, the operator S m ( ) can be realized as a birational automorphism of T −,m , so we do not have to consider any completion.
In all of our examples, the wall-crossing formula will be automatically satisfied as a result of the computation of [IK] -in particular, the BPS indices do not vary with m at all, so the identity is trivially satisfied.
2.5. BPS τ -function. Given a sufficiently nice variation of BPS structures, and a family of solutions to the corresponding BPS Riemann-Hilbert problems, Bridgeland [Br1] defined a notion of τ -function encoding the solution in a single piecewise meromorphic function. To recall this, let us first introduce some terminology.
A variation of BPS structures ( m , Z m , m ) over a complex manifold M is said to be framed if the local system ( m ) m∈M is trivial. For a framed variation of BPS structure, we can identify all the lattices m with a fixed lattice through the parallel transport. Using the identification, we will denote a framed variation of BPS structures by ( , Z m , m ).

Definition 2.6. A framed variation of BPS structures
is a local isomorphism.
Now, suppose we are given a framed, miniversal variation of BPS structure ( , Z m , m ) over M. The miniversality property ensures we have the isomorphism and dually we may identify T * m M with ⊗ C. Then, the pairing . . . , . . . on induces a Poisson structure on M. If we take a (C-)basis γ 1 , . . . , γ r of and use (z 1 , . . . , z r ) = (Z m (γ 1 ), . . . , Z m (γ r )) as a coordinate on M, then the Poisson bracket is explicitly given as For a family of solutions {X (m, )} m∈M to the BPS Riemann-Hilbert problem, we denote by Y (m, ) ∈ T + the "non-twisted part" This is a collection of piecewise meromorphic functions (one for each γ ∈ ) with discontinuities at m ∈ M such that the ray is a BPS ray for the corresponding BPS structure ( m , Z m , m ). Thanks to the isomorphism (2.20), for any non-BPS ray , defines a vector field on M with discontinuities at this locus. Then, we can associate the τ -function corresponding to the family as follows.
Definition 2.7 ([Br1, § 4.6]). Suppose X (m, ) is a family of solutions to the BPS Riemann-Hilbert problem associated to a framed miniversal variation of BPS structures over M. A BPS τ -function for the family of solutions X (m, ), associated to a non-BPS ray , is a possibly-multivalued meromorphic function defined so that the Hamiltonian vector field of the function 2πi log τ BPS, coincides with the vector field ∂ log Y , where the Poisson structure on M is given by (2.21).
In terms of the coordinates above (for a chosen basis γ 1 , . . . , γ r ∈ ), we may write the condition for being a τ -function as In general, it is unknown if there exists a τ -function, and uniqueness is not guaranteed. However, if we additionally impose that τ BPS, is invariant under simultaneous rescaling of the coordinates z i and , and if the antisymmetric form ·, · is nondegenerate, then the τ -function is uniquely specified up to a multiplication by an element in C * .
Recall that Barbieri [Ba] provided a canonical piecewise-meromorphic solution to the BPS Riemann-Hilbert problem for arbitrary uncoupled BPS structures. The τ -function was not computed in that work, but considered in a recent paper [AP1] by Alexandrov-Pioline. The formula is given as follows (see also [BaBrS]).
Proposition 2.8 (cf. [AP1,§ 4]). Let ( , Z , ) be any miniversal variation of finite, uncoupled, integral BPS structures. A BPS τ -function for the family of minimal solutions to the BPS Riemann-Hilbert problem is given by where ϒ is given explicitly in terms of the Barnes G-function ([Bar]): ( 2.28) The formula (2.27) is a consequence of the relation d dw log ϒ(w, η) = w d dw log (w, η). We will also use the following formula which describes the asymptotic expansion of (2.28): for w → ∞, valid when w → ∞ in the complement of any open sector containing R <0 . See [AP1,Appendix B.4] for these properties of ϒ-function. We also note that, when ξ(γ ) = 1 holds for any active class γ , (2.27) is reduced to the BPS τ -function associated with the holomorphic solution (2.17) which was originally studied by Bridgeland in [Br1]: . (2.30) Here we write ϒ(w) = ϒ(w, 0).
Remark 2.9. In [Br2], Bridgeland proved that the τ -function for a BPS Riemann-Hilbert problem arising from the Donaldson-Thomas theory for the resolved conifold gives a "non-perturbative partition function" for the resolved conifold. This means that the asymptotic expansion of the BPS τ -function coincides with the generating series of all genus Gromov-Witten invariants. We expect there should be an analogue of our story in that setting as well.

BPS Structures from Spectral Curves of Hypergeometric Type
In this section, we will recall the definition and computation of BPS structures associated to spectral curves of hypergeometric type, as described in detail in our previous work [IK].

Spectral curves of hypergeometric type.
Let us summarize the geometric structure (i.e., properties of cycles and paths on) the spectral curves of hypergeometric type which are collected in Table 1. We sometimes borrow terminology from exact WKB analysis; see [KT] for background.
In the rest of this section, we denote by Q(x) one of the rational functions Q • (x) for • ∈ {HG, dHG, Kum, Leg, Bes, Whi, Web, dBes, Ai}, and by ϕ = Q(x)dx 2 the associated meromorphic quadratic differential on X = P 1 . We also introduce the following. 6 • P ⊂ X is the set of poles of ϕ. It has a decomposition P = P ev ∪ P od into the subsets of even order poles and odd order poles of ϕ. Note that all our examples in Table 1 satisfy ∞ ∈ P. • T ⊂ X is the set of zeros and simple poles of ϕ. Borrowing a terminology in WKB literature, we call elements in T turning points since they will be turning points of the quantum curves E = E • constructed in §4 below. • The set Crit = P ∪ T is called the set of critical points of ϕ. Its subset Crit ∞ = P\T consists of infinite critical points of ϕ; that is, the poles of ϕ of order ≥ 2. The points in T = Crit\Crit ∞ are also called finite critical points in this context.
The quadratic differential ϕ is parametrized by a tuple of parameters m = (m s ) s∈P ev which we call mass parameters (temperatures in TR language). We will write ϕ = ϕ(m) when we emphasize the dependence on the mass parameters m. In what follows, we assume the following condition on the mass parameters. (3.1) ( 3.7) where we set For degenerate Bessel curve (dBes) and Airy curve (Ai), we set M dBes = M Ai = ∅ since they have no mass parameters.
It is not hard to check that ϕ has only simple zeros, and no collision of zeros and poles occurs under the Assumption 3.1. Let us consider the Riemann surface associated with the quadratic differential ϕ. This is a branched double cover of X \ P with the projection map π : (x, y) → x. We denote by ι : (x, y) → (x, −y) the covering involution of . We will mainly use x as a local coordinate 7 on \ π −1 (T ). We note that, the Riemann surface for all the examples in Table 1 are of genus 0, and hence topologically are punctured spheres, with punctures corresponding to poles of ϕ. We denote by the compactification of , obtained by adding preimages of poles of ϕ, on which the square root √ ϕ = √ Q(x) dx gives a meromorphic differential. More precisely, to obtain , we add to a single point for each s ∈ P od , and add a pair s ± of two points for each s ∈ P ev : (3.10) Here we identify s ∈ P od with its unique lift on . We keep using the same notations π : → X and ι : → for the projection map and the covering involution. The label of two points s ± , both of which are mapped to s ∈ P ev by π , is chosen so that the following residue formula holds: ( 3.11) For later purposes, we also introduce a partial compactification of obtained by filling punctures corresponding to simple poles of ϕ: (3.12) Here we set If ϕ has no simple poles, then = .
Example 3.2 (Gauss hypergeometric curve [IKoT2,§2.3.1]). The Gauss hypergeometric curve HG is defined from the meromorphic quadratic differential ϕ HG = Q HG (x)dx 2 where the rational function (3.14) Under the assumption m = (m 0 , m 1 , m ∞ ) ∈ M HG , the poles are P = P ev = {0, 1, ∞} and the turning points are T = {b 1 , b 2 }, where b 1 , b 2 are the two simple zeros of Q HG (x) (since there are no simple poles, we have P ∩ T = ∅). Topologically, HG is a sphere with six punctures, and the compactification is given by We note that HG = HG holds, and we have 7 Technically speaking, we must fix branch cuts and identify the coordinate x of the base P 1 with a coordinate of \ π −1 (T ) on the first sheet (then, ι(x) gives a coordinate in the second sheet).

3.2.
Cycles on the spectral curve. Let us fix notation for (relative) homology classes on which will be used to describe the BPS structure in the rest of this section, and to define the Voros coefficients in §4. For spectral curves of hypergeometric type, since is a punctured sphere, the homology group H 1 ( , Z) is generated by the residue cycles γ a (i.e., the class represented by a positively oriented small circle) around the puncture a ∈ D ∞ satisfying the relation The covering involution of also acts on the homology group, and in particular, we have ι * γ s = γ s for s ∈ P od and ι * γ s ± = γ s ∓ for s ∈ P ev . On the other hand, for each s ∈ P ev , we may associate a relative homology class β s ∈ H 1 ( , D ∞ , Z) which is represented by a path from s − to s + on . The action of covering involution on this class is given by ι * β s = −β s . Among various relative homology classes, we are particularly interested in these classes.
Note that the intersection pairing ·, · on H 1 ( , Z) is trivial. There is another intersection pairing defined through Poincaré-Lefschetz duality: Z). The intersection pairing between the above (relative) homology classes is given by, 8 where δ in the right hand side is the Kronecker symbol.

BPS structures from quadratic differentials.
In the previous work [IK], we constructed BPS structures from the quadratic differentials in Table 1. Here we briefly recall the construction based on the spectral network associated with these quadratic differentials. We note that the construction of BPS structure given here works more generally (see [Br1,§ 7]).

Lattice and central charge.
Recall that a BPS structure requires a charge lattice and a central charge Z . Following [GMN2,BrS], we define: equipped with the intersection pairing ·, · : × → Z. The central charge Z : → C is the period integral of √ ϕ: The lattice is called the hat-homology group in [BrS]. It will sometimes be convenient to extend the domain of definition of Z from to the whole lattice H 1 ( , Z) in a natural manner. Since any element γ ∈ can be written as a sum of γ s ± in H 1 (˜ , Z), we may use the formula which follows from (3.11), to express the central charge Z (γ ) for a general element. We will write = m and Z = Z m when we discuss the dependence on the mass parameters.
As we have seen in previous subsection, the intersection pairing ·, · is trivial in our example. Thus the BPS structures we will obtain are automatically uncoupled.

Spectral networks.
To define the BPS indices, we use the notion of a WKB spectral network [GMN3] associated to the quadratic differential ϕ considered in the previous subsection. 9 This can also be understood as the Stokes graph (c.f., [KT]) of the quantum curve E(ν) which will be constructed in §4 through the topological recursion.
Let us define a trajectory of the quadratic differential to be any maximal curve on X along which ( 3.22) is satisfied. There is a well-constrained possible local behaviour of these trajectories, but a rich global behaviour. The general structure of trajectories of quadratic differentials is described in detail in the classic book by Strebel [Str] (see also [BrS]). We are primarily interested in trajectories with at least one endpoint on a turning point -we call these critical trajectories. The set of (oriented) critical trajectories is called the spectral network W ϑ (ϕ) of ϕ at phase ϑ. In the WKB literature, the critical trajectories are called Stokes curves, and form the locus where the Borel-resummed WKB solutions have a discontinuity.
In this paper, we will only use the end result of the previous work [IK], so we will not give details and figures of spectral networks, but we refer the reader to that paper, and proceed to write down the BPS structure that was obtained there.

Degenerate spectral networks and saddle trajectories.
We say the spectral network W ϑ (ϕ) is degenerate if it contains a saddle trajectory (i.e., a critical trajectory both of whose endpoints are turning points) of phase ϑ. These degenerate spectral networks are of paramount importance in this paper and many applications. In the physics of 4d N = 2 QFTs, they correspond to BPS states in the spectrum of the theory [GMN1,GMN3]. From a mathematical point of view, they correspond to stable objects in a 3-Calabi-Yau category associated with a quiver with potential determined by ϕ [BrS].
Let us summarize the possible types of saddle trajectories which can appear in degenerate spectral networks arising from our examples 10 : • A type I saddle connects two distinct simple zero of ϕ.
• A type II saddle connects a simple zero and a simple pole of ϕ.
• A type III saddle connects two distinct simple poles of ϕ.
• A type IV saddle is a closed curve which forms the boundary of a degenerate ring domain (i.e., a maximal region foliated by closed trajectories around a second order pole of ϕ).

Definition of BPS invariants.
Given a degenerate spectral network with a saddle trajectory of type I, II, or III, we can associate a homology class γ BPS ∈ represented by the closed cycle (up to its orientation) on obtained as the pullback by π : → X of the saddle trajectory. On the other hand, for a degenerate spectral network with a type IV saddle, we define γ BPS = γ s + − γ s − ∈ (up to its orientation) where s is the second order pole of ϕ in the degenerate ring domain surrounded by the type IV saddle. We refer to homology classes γ BPS ∈ defined in this way as BPS cycles, and the collection of all BPS cycles (together with their type) as the BPS spectrum of ϕ. Degenerate spectral networks only appears when the phase ϑ coincides with the argument of the central charge Z (γ BPS ) of a certain BPS cycle.
To define the BPS invariants, we should avoid a locus in the parameter space M of ϕ = ϕ(m) on which a multiple degeneration appears in the spectral network simultaneously. We say the parameter m is generic if it lies in the complement M \ W of a set W = W • defined by (3.23) for • = Leg, and W Leg = ∅ for the Legendre case. Otherwise, we say m is non-generic. We denote by M = M \ W the locus of generic parameters. Now we recall the definition of the BPS indices introduced in [GMN2, GMN3,BrS,IK] (see also the recent work [Ha] which arrives at the same BPS indices for BPS cycles associated with Type II and III saddles via Donaldson-Thomas theory).
Definition 3.4. For each m ∈ M , we define the BPS indices { (γ )} γ ∈ as a collection of integers defined as follows: if γ is a BPS cycle associated with a type I saddle, +2 if γ is a BPS cycle associated with a type II saddle, +4 if γ is a BPS cycle associated with a type III saddle, −1 ifγ is a BPS cycle associated with a degenerate ring domain.
In our previous work, [IK], given any of the examples of ϕ from Table 1, we determined the BPS structure following (3.24) at any m ∈ M . The results are summarized in Table 3 below. We refer the reader to [IK] for additional details.
Theorem 3.5 ( [IK,§ 4]). Let • be the label of any spectral curve of hypergeometric type in Table 1. For each m ∈ M • , the BPS structure determined by ϕ • (m) = Q • (x)dx 2 is the one given in Table 3. In particular, the BPS structure obtained is finite, integral and uncoupled.
As before, for s ∈ P ev , γ s ± is the positively oriented small circle around s ± The uncoupledness of the BPS structures in Table 3 is a consequence of the vanishing of the intersection pairing ·, · on . In particular, this vanishing implies that the associated BPS Riemann-Hilbert problem is trivial (i.e., all BPS automorphisms are identity maps). Therefore, for an arbitrarily chosen constant term ξ ∈ T − , the functions attached to each ray ⊂ C * give a solution to the BPS Riemann-Hilbert problem, and it is the unique holomorphic solution which has no zeros or poles on H . To obtain a non-trivial BPS Riemann-Hilbert problem, which will be compared with the Stokes structure arising from topological recursion and quantum curves, we will consider the almost-doubled BPS structure discussed next.
3.4. Doubled and almost-doubled BPS structure.. Given a BPS structure, we may construct a new doubled BPS structure ( D , Z D , D ) via the so-called doubling construction [Br1,§ 2.8]. The doubled lattice is: Clearly the symmetry property continues to hold, and the support property is also equivalent to the support property for the original BPS structure.
We will in fact use a slightly modified doubling procedure, which we call almostdoubling, in which we take a full-rank sublattice in the second factor, and define D := ⊕ * (3.31) equipped with the restricted pairing (denoted by the same symbol ·, · ). Together with the obvious restrictions Z D and D of Z D and D to the sublattice D , we obtain a BPS structure ( D , Z D , D ) which we call the almost-doubled BPS structure associated to the original. The purpose of this restriction will be to ensure our solution to the BPS Riemann-Hilbert problem is single-valued; see Remark 5.5 below. Note that the relative cycles β s defined for s ∈ P ev in §3.2 are examples of paths satisfying ι * β s = −β s since the covering involution exchanges s + and s − . In fact, we can verify that {β s } s∈P ev generates the sublattice * .
The doubling and almost-doubling construction is useful since it produces nontrivial BPS Riemann-Hilbert problems from BPS structures which may have a trivial Riemann-Hilbert problem. In particular, even though the integers (γ ) γ, γ may be zero for the original BPS structure, it is easy to see that (almost-)doubling will allow the analogous expression to take on nonzero values. This is why our introduction of (γ ) = −1 plays a significant role, but did not affect previous results on the subject -any choice of attributed to such BPS cycles would be killed by the intersection pairing in the original BPS structure.
3.5. Miniversality of almost-doubled BPS structure. In [IK] and above, the BPS structures corresponding to hypergeometric type spectral curves were considered as depending on a set of parameters m, but we did not consider them as a family. However, they naturally fit together to form a variation of BPS structures: Proof. Since no collision of zeros and poles of ϕ occurs on the parameter space M, the collection { m } m∈M of lattices forms a local system on the space. The rank of m are computed by using [BrS,Lemma 2.2], and we may verify that the rank of each lattice m coincides with the dimension of M (i.e., the number of even order poles of ϕ) by case-by-case checking. The triviality of this local system follows from the isomorphism described below in (3.37). Moreover, the explicit expression of central charges in Table  3 shows that the period map is holomorphic and local isomorphism.
Although m was not defined for m ∈ W since a multiple degeneration of spectral network can occur there, the definition of m can be extended by continuity across W due to our observation in [IK] which tells the set of BPS cycles and their BPS invariants are common for all m ∈ M as is shown in Table 3 (that is, Table 3 can be taken as the definition of the BPS structure for all m ∈ M). The family of BPS structure on M thus obtained satisfies the support property and wall-crossing formula since our BPS structures are finite and uncoupled (all BPS automorphisms commute). Thus we have a miniversal variation of finite, integral, uncoupled BPS structures over M.
Thanks to the miniversality of the original BPS structure, we have an identification between the tangent space T m M and * m ⊗ C ( Note that the same statement holds for the fully doubled family as well. 3.6. Almost-doubled BPS Riemann-Hilbert problem and the BPS τ -function. We will be interested in the BPS Riemann-Hilbert problem (Problem 2.2) associated with the almost-doubled BPS structure constructed in the previous subsection. We will call it the almost-doubled BPS Riemann-Hilbert problem below.
First of all, it is easy to observe that, the almost-doubled BPS structure ( D , Z D , D ) constructed above is still uncoupled. This is because the original BPS structure ( , Z , ) is uncoupled, and the corresponding BPS indices D are supported on the BPS cycles in the original BPS structure. Therefore, we may happily use the general formula in Theorem 2.4 for the minimal solution even in the (almost)-doubled BPS Riemann-Hilbert problem, once we have fixed an extended constant term ξ D which takes value in the almost-doubled twisted torus (3.32) whose elements may be written as ξ D = (ξ, ξ ∨ ) through the splitting (3.31). Applying Theorem 2.4, we have the following expression of the minimal solution to the almostdoubled BPS Riemann-Hilbert problem associated with ( D , Z D , D ) evaluated at μ ∈ D : or more explicitly in terms of the splitting D = ⊕ * : In §5, we will compare this minimal solution to our Voros solution which arises from the topological recursion and quantum curves. Note from the expression the jumping structure is clear: X min ,γ do not jump at all and are very simple, whereas the X min ,β jump across BPS rays of the original ( , Z , ).
To describe the BPS τ -function, let us use the mass parameters (m s ) s∈P ev which is a natural global coordinate on M (as we have seen, dim M coincides with the number of even order poles of ϕ.) To realize the coordinate as the central charges of a basis, it is convenient to use cycles with rational coefficients; that is, if we define a lattice can = s∈P ev Zλ s , (3.36) together with a linear map Z can : can → C given by Z can (λ s ) = 2πim s , then we have an isomorphism which identifies Z can with Z m (after ⊗Q) for any m ∈ M. We need such rational coefficients since {γ s + − γ s − } s∈P ev is not a Z-basis of m in general. In view of β s , γ s + − γ s − = 2 δ s,s , we may also canonically identify the (almost) dual lattices ∨ can * m by identifying the dual element of λ s with β s . Thanks to the miniversality, we have the following identification: Therefore, the BPS τ -function does not depend on m ∨ s . In particular, the BPS τ -function associated with the minimal solution (3.33) has the same expression given in Proposition 2.8; that is, it is written as a product of ϒ-functions over the original lattice . We will set Z ∨ = 0 (i.e., m ∨ s = 0) in Sect. 5, but this restriction evidently does not change the BPS τ -function.

Topological Recursion, Quantum Curves and Voros Coefficients
In the works [IKoT1,IKoT2], the first named author and collaborators studied the quantization of spectral curves of hypergeometric type via topological recursion (TR). See also [BoE2] for generalities on the quantization of spectral curves, i.e. quantum curves. Moreover, an explicit expression for the Voros symbols of the associated quantum curves was obtained. In this section, we recall the topological recursion and the construction of the quantum curves, as well as the definition and form of their Voros symbols from those works, which we will use to solve to the Riemann-Hilbert problem in §5.1.

Definition of correlators and partition function.
Here we briefly recall several facts from Eynard-Orantin's theory of topological recursion ( [EO1]) and define the objects playing a central role on the in the present paper.
Definition 4.1. A spectral curve is a tuple (C, x, y, B) of the following data: • A compact Riemann surface C, • A pair of non-constant meromorphic functions x, y on C such that dx and dy never vanish simultaneously, and • A symmetric meromorphic bidifferential 11 B on C having a second order pole with biresidue 1 along the diagonal, and holomorphic elsewhere.
We usually denote by z a local coordinate of C, and by z i a copy of the coordinate when we consider multidifferentials. Here, a meromorphic multidifferential is a meromorphic section of the line bundle π * 1 (T * C)⊗· · ·⊗π * n (T * C) on C n , where π j : C n → C is the jth projection map [DN]. We often drop the symbol ⊗ when we express multidifferentials.
We denote by R the set of ramification points of x, that is, R consists of zeros of dx; here we consider x as a branched covering map C → P 1 , and we call the images of ramification points by x the branch points. We also assume that all ramification points of x are simple so that the local conjugate map z →z near each ramification point is well-defined. Then, the topological recursion (TR) is formulated as follows.

Definition 4.2 ([EO1, Definition 4.2]
). The Eynard-Orantin correlator W g,n (z 1 , . . . , z n ) for g ≥ 0 and n ≥ 1 is defined as a meromorphic multidifferential on C by the recursive relation (4.1) for 2g + n ≥ 2 with initial conditions given by W 0,1 (z 0 ) := y(z 0 )dx(z 0 ), W 0,2 (z 0 , z 1 ) := B(z 0 , z 1 ). (4.2) Here we set W g,n ≡ 0 for a negative g, K r denotes the so-called "recursion kernel" given by defined near a ramification point r ∈ R, denotes the disjoint union, and the prime on the summation symbol in (4.1) means that we exclude terms for (g 1 , I 1 ) = (0, ∅) and (g 2 , I 2 ) = (0, ∅) (so that W 0,1 does not appear) in the sum. We have also used the multi-index notation: for One of the central objects in Eynard-Orantin's theory is the genus g free energy F g (g ≥ 0) associated to the spectral curve: . For g ≥ 2, the genus g free energy F g is defined by where (z) is any primitive of y(z)dx(z). The free energies F 0 and F 1 for g = 0 and 1 are also defined, but in a different manner (see [EO1,§ 4

.2.2 and § 4.2.3] for the definition).
We note that the right-hand side of (4.4) indeed does not depend on the choice of (z).

Definition 4.4. The generating series
of F g is called the free energy of the spectral curve. Its exponential is called the (topological recursion) partition function of the spectral curve.
The free energy and the partition functions are regarded as a formal series of . It is known that usually these formal series are divergent (see [E2] for the growth estimate for the coefficients). The objects arising in our examples are, nonetheless, Borel summable in all but finitely many "singular" directions.

Applying TR to spectral curves of hypergeometric type.
In the rest of §4, we denote by one of the curves • for • ∈ {HG, dHG, Kum, Leg, Bes, Whi, Web, dBes, Ai} in Table 1. To apply TR to , we should regard them as spectral curves in the sense of Definition 4.1. This was done in [IKoT2,§ 2.3], where an explicit rational parametrization of was given. That is, for each • in Table 1, there exists a pair of rational functions (x(z), y(z)) = (x • (z), y • (z)) such that we have an isomorphism of punctured Riemann surfaces, with the choice C = P 1 and the set P of poles of x(z) and y(z). Under the notations fixed in §3.2, the set P is of the form where, for each s ∈ P, either x( p s ) = s or x( p s ± ) = s holds depending on the parity of pole order of ϕ at s. The set P is in bijection with the set \ , and hence, (4.7) can be extended to an isomorphism C ∼ = of compact Riemann surfaces. We often use this isomorphism to identify the two in what follows. In the case of HG , the above construction is realized as follows. (See [IKoT2,§ 2

.3] for the others examples).
Example 4.5. A rational parametrization (4.7) of HG is given by the pair of explicit rational functions where HG (m) is given in (3.8) (which is non-zero under Assumption 3.1), and the set , (4.10) The points p s ± ∈ C are the two preimages of s ∈ {0, 1, ∞} by x HG (z), and the labels are chosen so that the residue formula (3.11) hold. The point p s ± is mapped to s ± through the isomorphism C HG . The set of ramification points is given by R = {±1} ⊂ C, and these two points are mapped to the turning points b 1 , b 2 ∈ P 1 by x HG (z). The conjugate map is given by z = 1/z, which corresponds to the covering involution of HG .
Together with the canonical choice of the bidifferential B when C = P 1 , we obtain a spectral curve (C, x, y, B) in the sense of Definition 4.1. We will identify defined in (3.9) (rather, its compactification ), together with the projection on the first coordinate π : → P 1 , with the spectral curve (C, x, y, B) through the isomorphism and the canonical choice of B. Under the isomorphism, the set R ⊂ C of ramification points is mapped to π −1 (T ) ⊂ bijectively.

Free energies of spectral curves of hypergeometric type and BPS structure.
We have now introduced all of the notions needed to recall the result of our previous work [IK]. In [IKoT1,IKoT2], a relationship between the free energy of the spectral curves of hypergeometric type and the Voros coefficients of the associated quantum curves was found (see Remark 5.14 for the precise relation). As a corollary, there are explicit formulas of gth free energy F g = F • g , which we can write uniformly in terms of BPS structures thanks to the results in [IK] for all hypergeometric type examples: Theorem 4.6 ( [IKoT2,IK]). The gth free energy of the spectral curves of hypergeometric type are given explicitly as follows: (4.14) Here the formula for F 0 (resp., F 1 ) is valid up to quadratic polynomials of m i s (resp., additive constants), B k is the kth Bernoulli number given by B k = B k (0), and H is any half-plane whose boundary rays are not BPS.
We note that the prefactor in (4.14) given by the Bernoulli number has an interpretation as the Euler characteristic of the moduli space of Riemann surfaces ( [HZ,P]). The explicit expression in each example can be recovered from the BPS spectrum computed in [IK], summarized in Table 3 (the Airy and degenerate Bessel cases are excluded since the free energy is trivial). These expressions are effectively used to compute the Borel sum of their generating series F TR = g≥0 2g−2 F g . Below, we will introduce the Voros coefficients of the corresponding quantum curves, and write an analogous explicit formula for them.

Quantum curves of hypergeometric type.
It is known since [GS,Z,MS,DM,BoE2] etc. that the Eynard-Orantin correlators are closely related to the expansion coefficients of the WKB solution of a certain Schrödinger-type equation, called a quantum curve, whose classical limit gives the spectral curve. Here we review the statement for the Gauss hypergeometric spectral curve HG (see [IKoT2,§ 2.3] for other examples).
First we introduce the following notations for various subsets which exclude turning and ramification points: (4.16) on C which depends on z ∈ C and a tuple ν = (ν p ) p∈P of complex parameters, which we call quantization parameters, satisfying p∈P ν p = 1. (4.17) Let W g,n (z 1 , . . . , z n ) be the correlator defined by TR from the spectral curve (C, x, y, B) described in §4.2 and [IKoT2]. Let us define a "primitive" of W g,n by 12 F g,n (z 1 , . . . , z n ) Restricting this to a diagonal set z 1 = · · · = z n = z, we have a function on C . Taking their generating series, we introduce 12 For a second kind meromorphic differential ω(z) on C which is holomorphic at p ∈ P , we shall define its integral with the divisor D(z; ν) for z ∈ C\{poles of ω} by ζ ∈D(z;ν) ω(ζ ) = p∈P ν p ζ =z ζ = p ω(ζ ). The condition (4.17) implies that d ζ ∈D(z;ν) ω(ζ ) = ω(z). We also generalize the definition for multidifferentials in a straightforward way. as a formal (Laurent) series in whose coefficients are functions 13 on C .
Roughly speaking, the framework of quantum curves claims that the exponential of S TR (z, ) satisfies a Schrödinger-type ODE (after combining this with the inverse map of the isomorphism C \P ∼ → (which is a restriction of (4.7))), which has the original spectral curve as its classical limit. We will not give a general statement of quantum curves in [BoE2,IKoT1] here since it requires considerable notational preparation. Instead, we focus on the main example, the quantization of the Gauss hypergeometric curve HG .

Then, its exponential (for any choice of local inverse z(x))
gives a formal solution of the following differential equation: The equation E HG given in (4.21) is a 2nd order linear ODE with three regular singular points at 0, 1 and ∞, and hence, it is equivalent to the Gauss hypergeometric differential equation (see [IKoT2,§ 2.3.1]). We note that the Gauss hypergeometric curve (3.14) is nothing but the classical limit → 0 of the equation E HG ; that is, E HG is a quantization of the Gauss curve (3.14). We call E HG the quantum Gauss curve (for the divisor D(z; ν)).
In [IKoT1,IKoT2], it was proved that this quantization procedure works for all spectral curves = • shown in Table 1, and the resulting quantum curves E are members of the confluent family of hypergeometric type differential equations (there is some overlap with the results of [BoE2]). We refer to [IKoT2,§ 2.3] for explicit descriptions of rational parametrizations of hypergeometric type curves, choice of divisors, and the associated quantum curves for the other eight cases. We call these equations E the quantum curves of hypergeometric type (associated to each example).
Remark 4.8. It is well-known that a Schrödinger-type ODE with a small parameter like (4.21) has a formal solution of the form (4.22) called the WKB solution. For a construction of the WKB solution, we refer [KT,§ 2] for example. The derivative d S WKB,k (z(x))/dx of the coefficients in the WKB solutions are meromorphic functions on the compactified classical limit which are determined by a certain recursion relation. Theorem 4.7 claims that S WKB coincides with the formal series S TR in (4.19) obtained from the topological recursion, up to additive constants at each order of .

Voros coefficients of quantum curves of hypergeometric type.
Here we recall the definition of the Voros symbols for the quantum curve E constructed above. See [V,DDP1,IN1,IN2] for the definition for general Schrödinger-type ODEs, and their roles in the exact WKB analysis.

Definition of Voros coefficients.
We define two kinds of Voros coefficients; the first type are associated with homology classes γ ∈ H 1 ( , Z), while the second type are associated with relative homology classes β ∈ H 1 ( , D ∞ , Z). Following [IN1], we refer to these as cycles and paths, respectively. We will keep using notations like γ for cycles in H 1 ( , Z) and β for paths in H 1 ( , D ∞ , Z) throughout.
Consider the formal series S TR (z, ) defined in (4.19), and consider its term-wise derivative d S TR (z, ) = k≥−1 k d S TR,k (z) (= d S WKB (z, )) with respect to z. We take its "odd part" d S odd TR (z, ) = k≥−1 k S odd TR,k (z) defined by 14 where ι * means the pullback by the covering involution ι of . This is a formal series whose coefficients are meromorphic on satisfying Although the coefficients d S odd TR,k have poles on π −1 (T ), the residue of d S odd TR,k at those points vanish due to the property (4.24). Therefore, the integral of d S odd TR,k along a cycle γ ∈ H 1 ( , Z) is well-defined for any k ≥ −1. On the other hand, among the coefficients, d S odd TR,−1 and d S odd TR,0 have poles on π −1 (T ) ∪ D ∞ , while d S odd TR,k for k ≥ 1 have poles only at π −1 (T ). This is because the correlators W g,n with 2g −2+n ≥ 1 have poles only at the ramification points. Moreover, W g,n with 2g −2+n ≥ 1 have no residue at ramification points, and hence, the integral of d S odd TR,k along a path β ∈ H 1 ( , D ∞ , Z) is well-defined for k ≥ 1. Thus, we define (4.25) by subtracting off the k = 0, 1 terms. Then, we define the Voros coefficients for the quantum curves as follows.
• For any cycle γ ∈ H 1 ( , Z), the (cycle) Voros coefficient of the quantum curve E for γ is a formal series of defined by the term-wise integral: (4.26) • For any β ∈ H 1 ( , D ∞ , Z), the (path) Voros coefficient of the quantum curve E for β is a formal series of defined by the term-wise integral 15 : • The exponentiated formal series e V γ and e V β are called the Voros symbol of E for γ ∈ and β ∈ * , respectively. We refer to these as cycle Voros symbols and path Voros symbols, respectively.
For our purpose, we will be mainly interested in the Voros coefficients for cycles and paths in our previously chosen sublattices; that is, for γ ∈ and β ∈ * . We will recall the explicit form of the Voros coefficients in these cases.
Remark 4.10. In [IKoT1,IKoT2], the path Voros coefficients were defined by using d S TR instead of d S odd TR . Since the even part d S even TR := d S TR − d S odd TR is invariant under the covering involution ι * , its integral along β ∈ * vanishes. Hence, our path Voros coefficients (4.27) for β ∈ * agree with those computed in [IKoT1,IKoT2].

Cycle Voros coefficients for γ ∈ .
In what follows, we will use the same notation for cycles and relative cycles on the spectral curve defined in §3.2. Recall that the homology group H 1 ( , Z) is generated by the residue cycles around the punctures D ∞ . In particular, any element γ ∈ can be decomposed as a sum of residue cycles around s + or s − for s ∈ P ev . Through the residue computation, we have the following.
Proposition 4.11. For any quantum curve of hypergeometric type and any γ ∈ , we have with ν : H 1 ( , Z) → C the linear map given on generators by ν(γ s ± ) = ±ν s , where we set ν s := ν p s+ − ν p s − .
Proof. It is enough to show the formula for γ s ± with any s ∈ P ev . Since W g,n with 2g − 2 + n ≥ 1 have poles only at ramification points, we can verify that the residue of d S TR at s ± only comes from those of d S (4.29) The residue of the first term is given by due to the residue formula (3.11). The integrand in the second term is expressed as . (4.31) The summands for p = p s ∓ in (4.31) are holomorphic at z = p s ± , but the summand for p = p s ∓ has a simple pole at z = p s ± . Consequently, we have Res z= p s ± ζ ∈D(z;ν) (4.32) Since Res ) holds, we have the desired formula for γ s ± : (4.33) The general formula follows immediately from (4.33).
In particular, we note V γ is just a finite series in with two terms.

Path
Voros coefficients for β ∈ * . As we mentioned in §3.4, * is generated by β s , s ∈ P ev , which is a class represented by a path from s − to s + . The path Voros coefficients of the hypergeometric differential equations for those paths were studied in [SS,T,KoT,AT,KKT,ATT1,AIT,IKoT1,IKoT2] etc. Let us recall and rewrite the formula for their explicit expression, following [IKoT2,IK]. The corresponding quadratic differentials Q • (x)dx 2 have a unique even order pole s; thus m = (m s ), and we have the Voros coefficient for the canonical relative homology class β s ∈ * • associated with s. According to [IKoT2,Theorem 3.1 (iv)], the coefficients V β s ,k = V β s ,k (m s , ν s ) are explicitly given as follows: (4.37) Here B k (t) is the k-th Bernoulli polynomial defined through its generating series we tw (4.38) These explicit formulas are derived by finding and solving the difference equations (with respect to the -shift of parameters m) satisfied by the Voros coefficients (see [IKoT2] and the above references for details).
These two cases are of fundamental importance, since the Voros coefficients for other examples are obtained as superpositions. In fact, following our interpretation of the topological recursion free energies via BPS structures [IK] we may write the Voros coefficients computed in [IKoT2] uniformly as follows: Proposition 4.13. For any quantum curve of hypergeometric type, the coefficient of k in the Voros coefficient for any relative homology class β ∈ * is expressed as: where H is any half-plane whose boundary rays are not BPS.
for Bernoulli polynomials shows the expression (4.39) is independent of the choice of half-plane. Taking the sign convention (3.29) into account, the claim is proved immediately by comparing with the expression of Voros coefficients obtained in [IKoT2, Theorem 3.1 (iv)] (written explicitly in Appendix A), together with Table 3.
The expression (4.39) is crucially important for our purpose.

Borel sum of Voros coefficients.
Since the values of Bernoulli polynomials grow factorially as k tends to infinity, the Voros coefficients for relative homology classes are divergent series of . However, we can see that the Voros coefficients in our examples are Borel summable as a formal series of along all but finitely many rays. Here we briefly recall the notion of the Borel sum for our specific examples -for a general theory of the Borel summation, we refer [C,Sa]. We note that our computations presented here were essentially done in previous works [T,KoT,AT,ATT1,AIT] etc. which discussed the Voros coefficients for hypergeometric-type differential equations in slightly different presentations (with or without quantization parameters ν). First, for the Voros coefficient V β s ( ) = k≥1 k V β s ,k , let us define its Borel transform V β s (ζ ) as follows: Namely, V β s (ζ ) is the term-wise inverse Laplace transform, where ζ is the variable which is Laplace dual to 1/ . Although the Voros coefficient is divergent series of , V β s has a finite radius of convergence. For a given ray of the form = e iϑ R >0 , the Voros coefficient V β s is said to be Borel summable along if the Laplace integral . 16 This requires that V β s (ζ ) has an analytic continuation to a sectorial neighborhood of the integration contour e iϑ R ≥0 , and growth at most exponentially when |ζ | → ∞. The resulting integral (4.42) is called the Borel sum of the Voros coefficient, if it is well-defined. Watson's lemma (e.g., [BH,Chapter4]) implies that the asymptotic expansion of the Borel sum recovers the original divergent series when approaches 0: For the Borel sum of the Voros coefficients of our examples, the following is known: Proposition 4.14. For any example of hypergeometric type and any s ∈ P ev , the Voros coefficient V β s ( ) with the expression (4.39) and the corresponding Voros symbol are Borel summable along any ray which is not BPS in the corresponding BPS structure in Table 3. Moreover, the Borel sum S e V βs ( ) of the Voros symbol can be extended as a meromorphic function of on the whole half-plane H .
Proof. The statement was essentially shown (without the language of BPS structures) in the previous works [T,KoT,AT,KKT,ATT1,AIT], but we include a proof here for the convenience of reader. Note that, since the Voros coefficients in our examples are obtained by superposing V Web β ∞ and V Bes β 0 , it suffices to prove the statement for the two examples. For the Weber case, using the expression (4.37) of the coefficients, we have (4.44) Hence, V Web β ∞ is not Borel summable along the rays ±2πim ∞ R ≥0 , which are exactly the BPS rays for the BPS structure constructed from Web , since the Borel transform V Web β ∞ has a sequence of simple poles on the set 2πim ∞ Z =0 . The Borel sum is well-defined for any other ray = e iϑ R >0 , and explicitly computed as follows: (4.45) where (z, η) is the explicit function given in (2.14). The formula (4.45) is derived with the aid of the following Binet's integral formula for the logarithm of the -function (see [WW,p. 249] for example): Since the Borel summation commutes with taking exponential (c.f., [Sa,§ 13]), we have (4.46) Thus we have proved the claim in the Weber case. We can perform a similar computation in the Bessel case. We can verify that V Bes β 0 is Borel summable except for the BPS rays ±2πim 0 R >0 , and obtain: arg(2πim 0 ) < ϑ < arg(2πim 0 ) + π. (4.47) The other cases follow from the formula (4.39) On the other hand, since the Voros coefficient V γ = V • γ for γ ∈ is a finite series of (c.f., Proposition 4.11), its Borel sum along any ray is identical to itself in all examples: (4.48) 4.7. Stokes phenomenon for Voros symbols. The equality (4.48) shows that the cycle Voros coefficients do not exhibit any jumping behaviour when we vary the ray for the Borel summation. On the other hand, as we have seen in the proof of Proposition 4.14, the Borel sum of the path Voros coefficients have a discontinuity across the non-Borel summable directions when we vary . The jump of the Borel sum can be understood as the Stokes phenomenon for the divergent series V β s around = 0. The discontinuity follows from the presence of singular points of the Borel transform V β s (ζ ) in the ζ -plane. We have the following behaviour for the two fundamental cases of Weber and Bessel: Theorem 4.15 (c.f., [T,AIT]). Let • ∈ {Web, Bes}, and take a BPS ray BPS = ±2πim s R >0 of the corresponding BPS structure. We denote by 1 and 2 the rays obtained by perturbing BPS in counter-clockwise and clockwise direction, respectively (i.e., k = e iϑ k R >0 with ϑ 1 = arg(±2πim s ) + δ and ϑ 2 = arg(±2πim s ) − δ with a sufficiently small δ > 0). Then, we have the following relation (as meromorphic functions) between the Borel-resummed Voros symbols on H 1 ∩ H 2 : for • = Web and s = ∞, for • = Bes and s = 0. (4.49) Proof. We only give a proof for the Weber case based on the idea of [T] where the case ν ∞ = 0 was discussed. It follows from the expression (4.44) that the Borel transform has simple poles only on BPS rays. Through the residue calculus, we have Since the overall factor ∓1 is identical to minus the intersection pairing (3.18), we obtain (4.49) for the Weber case by taking the exponential. The same method can be applied to Bessel.
The jump formula (4.49) for the Borel-resummed Voros symbols is closely related to the one imposed in the BPS Riemann-Hilbert problem in §2.3. Indeed, when ν ∞ = 1, the Borel resumed Voros symbol S ϑ e V Web for the Weber curve, which is expressed by the function (2.14), coincides with the solution of the doubled A 1 BPS Riemann-Hilbert problem discussed in [Br1,§ 5] after identifying the mass parameter with the central charge (recall that (w, 0) = (w, 1) = (w)). For generic ν ∞ , the Voros symbol also gives the solution to the "meromorphic" BPS Riemann-Hilbert problem discussed by Barbieri [Ba,§ 4]. We will make this more precise in §5.1.
Remark 4.16. In the language of the resurgent analysis (c.f., [Sa]), the Voros coefficients for our examples are simple resurgent series, and the computation (4.50) in the above proof is closely related to the so-called alien calculus. See [T,KoT,KKKoT,KKT,AIT] for more discussion. Note that (4.49) is a consequence of the Euler reflection formula of the gamma function, as shown in [Br1,Ba]. The above computation gives a resurgencetheoretic proof of the formula.
4.8. The Voros potential. Since Voros coefficients are integrals over homology classes, we may consider V ·,k as a map V ·,k : ( ⊕ * ) ⊗ C → C, where we complexify in the obvious way. Recall that we are interested in the variation of BPS structures ( D , Z D , D ), over T * M • which by Corollary 3.7 is in fact miniversal. The miniversality of the original ( , Z , ) provides an isomorphism: (4.51) and dually we may identify T * m M and ⊗ C. In particular, this means the "interesting part" (namely, the path Voros symbols) V ·,k : ∨ ⊗ C → C (abusing notation) of the maps V ·,k may be interpreted as a one-form on M • .
Fixing generators {γ s − ι * γ s } for ⊗ C and taking as local coordinates m s = Z (γ s )/4πi on M • we may write these one-forms as It turns out these forms are closed. Thus, there is a formal series in of functions φ = k≥1 φ k k such that where d M • denotes the exterior derivative on M • only, and ω := k≥1 ω k k is a formal series of one-forms. Each term in φ is a function on the space M • depending on the quantization parameters ν. More precisely: Proposition 4.17. The Voros coefficients of quantum curves of hypergeometric type are generated by a potential 2πiφ, where φ = k≥1 φ k k is called the Voros potential. That is, the forms ω k , k ≥ 1 are closed, and the Voros potential is given by where H is any half-plane whose boundary rays are not BPS. Whenever is not a BPS ray, φ is Borel summable along , and 2πi S φ is a potential for S ω.
Proof. To determine closedness and compute φ k , we calculate directly from (4.53), caseby-case using the explicit formula (4.39) or Appendix A. In the Gauss hypergeometric and degenerate hypergeometric cases, we use the identity for Bernoulli polynomials . The Borel summability follows from the known asymptotic expansion (see e.g. [BaBrS]) of the Barnes double gamma function in a sufficiently large sector.
We will provide an interpretation of the Voros potential in the next section.

Solution to BPS Riemann-Hilbert Problem and the τ -Function
We may now formulate and prove our main results.

BPS indices and the Stokes structure of the Voros coefficients.
In the work of Gaiotto-Moore-Neitzke [GMN2], the BPS indices { (γ )} γ ∈ appeared in a formula which describes the jump of the Fock-Goncharov cluster coordinate on the moduli space of the framed local systems. In fact, it is easy to see this jump formula is in fact nothing but the required property (2.10) in the BPS Riemann-Hilbert problem. The purpose of this subsection is to solve the BPS Riemann-Hilbert problem, and thus show that our definition (3.24) of BPS indices proposed in [IK] is consistent with this point of view if we identify the Fock-Goncharov coordinates with the Borel-resummed Voros symbols (c.f., [I1,All2,Ku]). We consider the almost-doubled BPS structure ( D , Z D , D ) with the almost-doubled lattice D = ⊕ * . For convenience, in what follows we take Z ∨ = 0 so when we consider the variation of our almost-doubled BPS structures, they are regarded as a family over M (i.e., the zero section of T * M).
. denote any of the almost-doubled BPS structures of hypergeometric type. We define a meromorphic map, the Voros solution, as follows. For any non-BPS ray , set 17 for any cycles γ ∈ , and for any paths β ∈ * .
Here, σ : D → {±1} is the quadratic refinement characterized by the following property: and σ (β) = +1 for β ∈ * . Observe that some choice of quadratic refinement was necessary in order to convert the T D,+ -valued (e −V γ , e V β ) into something T D,− -valued. It is shown in [BrS,Lemma 3.2] that an appropriate subset of BPS cycles generates the lattice , and it is not difficult to see that such a quadratic refinement exists. The collection {X Vor ,μ } μ∈ D gives a map X Vor from H to the almost-doubled twisted torus T −,D .
We justify the terminology "solution" in the remainder of the section. First we show the jumping behaviour: . be an almost-doubled BPS structure of hypergeometric type, and ⊂ C * an acute sector whose boundary rays 1 , 2 , taken in the clockwise order, are non-BPS rays. Then we have the following relation on H 1 ∩ H 2 : Proof. The BPS cycles and BPS indices for hypergeometric type spectral curves are summarized in Table 3. Since our BPS structure is finite, it suffices to prove the claim when contains only one BPS ray. Also, it is enough to verify (5.5) for the generators of the almost doubled lattice; that is, for μ = (γ , 0) with γ ∈ and for μ = (0, β) with β ∈ * . We first note that the cycle Voros symbols, for homology classes μ = (γ , 0) with γ ∈ , do not jump as we have seen in §4.7. hence, they automatically satisfy the formula (5.5) since γ, γ = 0 for all γ ∈ , but (μ) = 0 for μ = (0, β).
Since * is generated by β s for s ∈ P ev , the remaining task is to verify (5.5) for μ = (0, β s ). It is easily see that this follows from the expression (4.39) of the Voros coefficients for those paths and the jump property (4.49) for the Weber and Bessel case. Table 3) and the sign convention (3.29) into account, the formula (4.49) can be written as where γ BPS is the unique BPS cycles whose central charge lies on BPS = ±2πim s R >0 (i.e., γ BPS = γ ∞ ± for • = Web, and γ BPS = γ 0 ± − γ 0 ∓ for • = Bes). Thus we have verified (5.5) for • = Web and Bes. Since the path Voros coefficients in the general case (4.39) is a superposition of those of Weber and Bessel, we can conclude that (5.5) is valid for all other cases as well.
Remark 5.3. From the explicit expressions (4.39) for the Voros coefficients, we have the explicit formula of their Borel sum for any non-BPS ray and any μ ∈ D : (5.7) The factor ξ D, ν ∈ T D,− is defined by where ν(γ ) was defined in Proposition 4.11 for γ ∈ , and we extend it to an element in Hom( D , C) by setting ν(β) = 0 for β ∈ * .
Remark 5.5. Although a fractional power appears in the second line of (5.7), we can verify that (5.7) is indeed a meromorphic function of since the pairing of γ s ± −γ s ∓ ∈ (only such cycles can give −1 as BPS invariant) with μ ∈ D is always even. This shows clearly why we had to consider the almost-doubled lattice rather than the (fully) doubled one. Otherwise, the powers of 1 2 appearing in the presence of loop-type BPS cycles give rise to non-meromorphic X Vor , since the intersection pairing on the doubled lattice is not always even.
Having shown that X Vor satisfies the jump property, we may conclude Theorem 5.6 . Let ( , Z , ) denote a BPS structure obtained from any spectral curve of hypergeometric type. Then X Vor is a meromorphic solution to the BPS Riemann-Hilbert problem associated to the almost-doubled BPS structure ( D , Z D , D ), with the constant term ξ D,ν given by (5.8).
Proof. We have already shown (RH1), and (RH2) follows from the explicit formula (5.7). The asymptotics at infinity follow from the explicit expression and the properties of the function (see, e.g. [Ba,§ 3].

Relation to the minimal and holomorphic solutions.
We note that while we have given a solution to the BPS Riemann-Hilbert problem, it does not always agree with the ones given by [Ba,Br1] which we wrote down in §2. However, we may relate them explicitly (for simplicity, we restrict Re ν(γ BPS )), noting in each case that the difference is at worst multiplication by a simple meromorphic function: Corollary 5.7. Fix a non-BPS ray . For all BPS cycles γ BPS ∈ with Z (γ BPS ) ∈ iH , suppose Re ν(γ BPS ) ∈ (−1, 1] whenever (γ BPS ) = −1, and Re ν(γ BPS ) ∈ (−2, 0] whenever (γ BPS ) = −1. Then, in the half-plane H , the Voros solution (5.7) and the minimal solution (2.16) to the BPS Riemann-Hilbert problem associated to ( D , Z D , D ) with constant term ξ D, ν are related via X Vor ,γ = 1 · X min ,γ , γ ∈ (5.9) X Vor ,β = ρ ,β · X min ,β , β∈ * (5.10) where ρ ,β : H → C is meromorphic, and given explicitly by (5.11) In particular, when ν s = 0 for all s ∈ P ev , the Voros solution agrees with the minimal solution -that is, X Vor = X min holds for any non-BPS .
Remark 5.8. Note, the relation for additional values of ν can be obtained by shifting the second argument in the function.
Proof. The first line of (5.9) is immediate. In order to prove the second line, we may restrict our attention to the Weber and Bessel cases thanks to the superposition structure (4.39). In view of (4.46), we have the following in the Weber case: where γ = γ ∞ ± is chosen so that Z (γ ) ∈ iH . Since (γ ) = 1, our assumption ensures that Im (πiν(−γ )+πi) ∈ [0, 2π), so that log ξ ν∞ (−γ ) 2πi . Thus the Voros solution agrees with the minimal solution On the other hand, in the Bessel case, we may rewrite the expression for the Voros coefficient as follows: (5.14) where γ = γ 0 ± − γ 0 ∓ is chosen to satisfy Z (γ ) ∈ i H , and we used the property (w, η + 1) = (w, η) (1 + η w ) to obtain the second line (recall the pairing is even on the almost-doubled lattice, so there is no issue of branching). Under the assumption on ν(−γ ), the first factor in (5.14) agrees with the minimal solution X min,Bes and hence, the second factor in (5.14) is ρ ,β 0 . By the evennness of the pairing on the almost-doubled lattice, this is a meromorphic function of . The general expression follows from the superposition structure, formula (4.39).
Corollary 5.9. For any spectral curve of hypergeometric type, there exists a choice of parameter ν * so that ξ D,ν * (γ ) = 1 (5.16) holds for all active γ , and the Voros and holomorphic solutions for the associated BPS Riemann-Hilbert problem are related, for any non-BPS ray , as: where β = β ( ) is either 1 or a simple meromorphic function independent of .
Proof. We need only check at the level of asymptotic expansions. The asymptotic expansion of X hol ,β is given by [Br1,§ 5.3]: Comparing this with equation (4.39), and recalling that B k+1 (0) = B k+1 (1) = B k+1 , we can proceed case-by-case. For the Weber and Whittaker cases there is only one term in the sum, so choosing ν ∞ = 1 (or ν ∞ = −1) gives exact agreement, and for Bessel ν 0 = 0 similarly gives which also agrees after accounting for the intersection number. In the case of Kummer we may demand B k+1 (γ BPS ) = B k+1 whenever γ BPS has (γ BPS ) = −1 which forces ν 0 , ν ∞ = 0 or 1 (but not both). Setting ν 0 = 0, ν ∞ = 1, we obtain the agreement. It is not hard to see that this choice satisfies (5.16). The simplest case in which a discrepancy between X Vor and X hol appears is Legendre. For concreteness, let us write out (4.39) in this case: where we used the identity B k (1) = B k and B k (ν + 1) = B k (ν) + kν k . Through the comparison to (5.19), we can see that the Borel sum of e V Leg β∞ along differs from X hol ,β ∞ by an additional term − log(1 − 2m ∞ ), which is independent of since the formal series itself is convergent. Its exponential (1 − 2m ∞ ) −1 is meromorphic function of and gives the factor β ∞ = Leg β ∞ . The degenerate Gauss and Gauss hypergeometric cases follow similarly, first choosing ν s = 0 or 1 in order to obtain agreement from terms with (γ ) = −1, and computing the resulting discrepancy appearing for terms with (γ ) = −1. In both cases, we may make an arbitrary choice of one of the second order poles and set the corresponding ν s = 1, with the remaining νs = 0. Thus we obtain a discrepancy of For the other poless, we have βs = 1. Again, we can verify that (5.16) is satisfied even in the cases where there is a non-trivial discrepancy β .

BPS τ -function and the Voros potential.
Having established the Voros solution to the almost-doubled BPS Riemann-Hilbert problem, we can ask for a corresponding BPS τ -function. Such a function exists, and it turns out there is a very natural interpretation on the TR side in terms of the Voros potential: Here, the prefactor c = c ( ) is independent of m and chosen explicitly as follows: (5.25) Proof. Using the definition of the τ -function for the almost-doubled problem, it is not hard to see the condition for being a τ -function is exactly the same as the condition (rightmost in (3.39)) for being a potential for −∂ log Y ·, , where Y ·, is the untwisted part of the solution and · is valued in * . For the Voros solution, this is the equation (4.53) for the Voros potential (up to the sign). Then from the explicit formulas, it is easy to check the equation holds at the level of asymptotic expansions, and since both sides are Borel summable, the relation holds analytically. The constant prefactor c is chosen so that τ Vor BPS is invariant under the rescaling (m s , ) → (λm s , λ ) for any λ ∈ C * . As in the previous section, we can also obtain the relation between τ Vor BPS and τ min BPS (which together with the explicit formula for τ min BPS in Proposition 2.8 may be used to compute the Voros potential explicitly). We have: where κ is a simple explicit function given by Here (5.29) follows from the asymptotic property (2.29) of ϒ-function, and (5.30) is given by the leading term of -derivative of the Voros potential (together with log c ; see Theorem 5.10). As we have seen in the proof of Corollary 5.7, holds under the same assumptions on ν. Therefore, using the properties B 2 (x) = B 2 (1− x) and B 2 (1 + x) − B 2 (x) = 2x of the Bernoulli polynomial, we have This fixes the overall constant factor in κ , and thus we have (5.27).
In particular, we see that the τ -functions do not agree for the Voros and minimal solutions, but differ only by a very simple factor coming from BPS cycles with (γ BPS ) = −1. where κ is either 1 or a simple explicit function independent of .
Proof. Replacing ρ with in the previous proof, we obtain κ = 1 for the Weber, Whittaker, Bessel, and Kummer cases. For Legendre, degenerate Gauss, and Gauss, we obtain where s ∈ P ev is a pole chosen as in the proof of Corollary 5.9.

Borel-resummed TR partition function as a BPS τ -function.
Let us finally show that the Borel-resummed topological recursion partition function essentially agrees with the BPS τ -function which we have constructed so far -this provides the "Borel-resummed version" of our previous result in [IK]. Proof. Using the asymptotic property (2.29) of ϒ-function, we can verify that τ hol BPS, has the following asymptotic expansion when → 0 in H : (5.38) We can verify that the right hand side is nothing but g≥1 2 g−2 F g by using Theorem 4.6. Since the asymptotic expansion of the ϒ-function is valid in a sufficiently large sector, this asymptotic property ensures that F TR and its exponential Z TR are Borel summable along , and the relation (5.37) holds tautologically.
Together with Corollary 5.11, the above formula tells us that the TR partition function is realized as the specialization of the BPS τ -function associated with the Voros solution.
Remark 5.14. We note that the specialization to the special value of ν is necessary, since the BPS τ -function in general clearly depends on ν, but τ TR depends only on free energies (and thus not on the quantization parameter ν ). We may offer a middle ground as follows. In [IKoT1,IKoT2], it was shown that the Voros coefficient can be recovered as the so-called "difference value" of the free energies. To simplify the notation, let us focus only on the Weber case, and refer to [IKoT2] for the others. It was found that (see [IKoT2]). The last two terms in (5.39) may be removed by a suitable redefinition of the Voros coefficient. Thus, at the level of asymptotic expansions, the difference value of τ hol BPS gives the Voros coefficient, i.e. the solutions to the Riemann-Hilbert problem for arbitrary choices of ν are encoded in the topological recursion τ -function τ TR (the BPS τ -function at ν = ν * ). This should be related to the observation of [BaBrS], who obtained a similar difference equation for the minimal solution (see also [A1,A2,AS] where closely related difference equations are discussed). Since it nonetheless encodes the solution for the full family of ν, one may consider this as an alternative way to define τ BPS , depending on which properties (geometric characterization, independence of ν) one desires the most.
Remark 5.15. By looking at the asymptotic expansion of the both sides of (5.24) evaluated at ν = ν * , we have where we use the same normalization (5.35) of F 1 on the right-hand side. Thus, we see there are two different relations between the solution to the BPS Riemann-Hilbert problem (or the Voros coefficients) and the BPS τ -function (or the TR partition function); that is, the difference equation discussed (5.39) and the relation (5.41) with the Voros potential. From the view point of special functions, this reflects the fact that the function and Barnes G-functions are related in two different ways (i.e., difference relation and integro-differential relation) as: 5.5. Higher degree spectral curves. We have so far only discussed variations of BPS structure arising from quadratic differentials. However, due to the work [GMN3] by Gaiotto-Moore-Neitzke, it is natural to expect that a tuple of higher differentials (or a higher degree spectral curve) also defines a BPS structure. Although many examples have been studied in e.g. [GMN3], there is still no mathematically rigorous description of the corresponding BPS structures in these cases, so such a picture remains conjectural. On the other hand, TR is applicable to higher degree spectral curves [BHLMR,BoE1], and its relationship to WKB analysis is also discussed in [BoE2] (under a certain admissibility assumption on spectral curves). The Borel summability and resurgence structure for the WKB solutions / Voros symbols of the associated quantum curve is not fully established, but it is natural to expect that they are controlled by an (appropriately generalized) spectral network (Stokes graph) as in the rank 2 cases (c.f., [BNR,AKT1,AKT2]).
Thus it is reasonable to ask if the results in this paper generalize to higher degree curves. Rather than formulate a precise conjecture, we will simply recall two higher degree examples considered numerically in Part I, and observe that (thanks to results of Y.M. Takei [Ta]) they behave as expected based on the results we have obtained so far.
Example: (1, 4) curve Let us consider the spectral curve arising as the classical limit of the 3rd order hypergeometric differential equation of type (1, 4), studied in [OK,Hi], which is explicitly given as follows: (1,4) : 3y 3 + 2t y 2 + x y − m ∞ = 0. (5.42) Here m ∞ is a parameter assumed to be non-zero, and we regard t as a constant. The curve (1,4) is a genus 0 curve with two punctures at ∞ ± , and we can verify that Res x=∞ ± ydx = ±m ∞ .
A crude numerical experiment performed in our previous work [IK] suggested a single degeneration occurs, the so-called "three-string web" also observed in [GMN3]. According to their result, the associated BPS cycle is simply the residue cycles γ ∞ ± , and the BPS index is assigned so that (γ ∞ ± ) = 1.
According to [IKo] and [Ta,Theorem 4.6], the explicit expression for the Voros coefficient of the quantum curve for the path from ∞ − to ∞ + is given as follows: (5.43) Since the formula is identical to that of the Weber curve, we have an explicit description via (4.46) of the jump property after taking the Borel sum. Taking an appropriate lattice, one may check that this Voros symbol solves the corresponding almost-doubled BPS Riemann-Hilbert problem.
Numerical experiment [IK] together with the general rules for spectral networks and BPS states predicted from physics (e.g. [GMN3]) suggests that (γ ) = 0 for all classes γ . [Ta,Theorem 4.6] computed the Voros coefficient explicitly, which is simply After taking the Borel sum, the Voros symbol does not jump, and provides a solution in accordance with the vanishing of the BPS indices (γ ).
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

Declarations
Data availability Data sharing is not applicable to this article as no datasets were generated or analyzed during the current study.
• For the quantum Kummer curve: • For the quantum Legendre curve: • For the quantum Bessel curve: • For quantum Whittaker curve: • For quantum Weber curve: • The Voros coefficients are not defined for the quantum Airy curve and quantum degenerate Bessel curve, since they have trivial BPS structure (i.e., = 0).