Calabi–Yau operators of degree two

We show that the solutions to the equations, defining the so-called Calabi–Yau condition for fourth-order operators of degree two, define a variety that consists of ten irreducible components. These can be described completely in parametric form, but only two of the components seem to admit arithmetically interesting operators. We include a description of the 69 essentially distinct fourth-order Calabi–Yau operators of degree two that are presently known to us.


Introduction
The hypergeometric series became famous after in the paper [19] its dual interpretation was discovered: on the one hand the series encodes enumerative information on rational curves on the general quintic Calabi-Yau threefold in P 4 and on the other hand it can be identified as a normalised period of the mirror quintic Calabi-Yau threefold. The power series is the unique (up to scalar) series solution of the hypergeometric operator L := θ 4 − 5 5 x(θ + 1 5 )(θ + 2 5 )(θ + 3 5 )(θ + 4 5 ), θ = x d dx .
The operator L is the first of the family of fourteen hypergeometric fourth order operators related to mirror symmetry for complete intersections in weighted projective spaces. These operators are all of the form where the integer N is chosen as to make the coefficients of the normalised holomorphic solution integral. The properties of these operators have been well-studied from various points of view [1], [40], [25], [18], [36], [37], [30].
The operator L is also the first member of the ever growing list of Calabi-Yau operators [7], [39], a notion that was introduced by the first author and W. Zudilin in [5] by abstracting the properties of L. Calabi-Yau operators in this sense are essentially self-adjoint fourth order fuchsian operators with a point of maximal unipotent monodromy for which strong integrality properties are supposed to hold. In a sense the term Calabi-Yau operator is somewhat of a misnomer, as there exist families of Calabi-Yau varieties, whose Picard-Fuchs operator have no MUM-point, hence are not Calabi-Yau operators in the above sense of the word. For a recent account, see [22] Let us consider a general Nth order differential operator written in θ-form: where the P i are polynomials of degree N. We call the largest r with P r = 0 the degree of the operator. The differential equation a n x n translates into the recursion relation P 0 (n)a n + P 1 (n − 1)a n−1 + . . . + P r (n − r)a n−r = 0 on the coefficients a n of the series φ(x), so the degree of L is equal to the length of this recursion.
The roots of P 0 (θ) are the exponents of the operator at 0. By translation one can define exponents of P at each point of the Riemann sphere P 1 . For the point ∞ we have to apply inversion x → 1/x; the exponents then are the roots of P r (−θ). The Riemann symbol is a table that contains for each singular point the corresponding exponents, that encode information about the local ramification of the solutions. Logarithmic terms may and usually do appear at points with exponents with integral difference. This leads to Jordan blocks in the local monodromy, but these are usually left implicit in the Riemann symbol.
We will usually assume that N = 4 and all exponents at 0 vanish, P 0 (θ) = θ 4 . The most salient feature of this situation is that at 0 we have a canonical Frobenius basis of solution φ 0 , φ 1 , φ 2 , φ 3 , where φ k contains terms with log k (x). The local monodromy has a Jordan block of maximal size, so 0 is a point of maximal unipotent monodromy: MUM. The 14 hypergeometric operators have degree 1 and have a Riemann-symbol of the form exhibiting a conifold singularity at the point 1/N, where the exponents are 0, 1, 1, 2 and the local monodromy around that point is a symplectic reflection; there is a Jordan block of size two. The exponents at ∞ are α 1 ≤ α 2 ≤ α 3 ≤ α 4 and add up in pairs to one: Currently we know over 500 Calabi-Yau operators with degrees running up to 40, but in no way do we expect the current list to be complete. An update of the AESZ-list is in preparation, [8].
At present we know 69 essentially distinct Calabi-Yau operators of order two, to which all further such degree two operators can be related by simple transformations. The operators for which the instanton numbers are 0 are obtained as Sym 3 of a second order operator and do not count as a proper Calabi-Yau operator and do not appear in [7] But contrary to the hypergeometric case, it appears that these degree two operators do not all fall in a single family. First of all, an operator of degree two can have four singularities, with a Riemann symbol of the form exhibiting two distinct conifold points. But in many cases the points a and b coalesce, producing an operator with only three singular points.
As one can observe in the above two examples, the exponents at infinity are symmetrically centered around the value 1 and this seemed to be the case for all further examples we found. It was M. Bogner who first found an example of an Calabi-Yau operator for which this is not the case: it is the operator θ 4 − x(216 θ 4 + 396 θ 3 + 366 θ 2 − 168 θ + 30) + 36 x 2 (3 θ + 2) 2 (6 θ + 7) 2 with Riemann symbol              0 1/108 ∞ 0 0 2/3 0 1/6 2/3 0 1 7/6 0 7/6 7/6              and later two more such operators were found. The existence of such operators was initially a surprise to us, as it appears that within the group of operators with three singular points further distinctions can be made. It is this circumstance that led to this paper.
The structure of the paper is as follows. In the first section we recall some basic facts about the Calabi-Yau condition and formulate the main result of this paper. In the second section we describe how the result can be obtained from a series of simple computations. In the third section we summarise the properties of the operators corresponding to each of the components that we find. In a final section we review the list all degree two Calabi-Yau operators that are presently known to us. In two appendices we include some further monodromy and modular form information on these examples, together with some properties of the lower order operators used in many constructions.

The Calabi-Yau condition
The adjoint L * of a differential operator L of the form L = d n dx n + a n−1 is obtained by reading the operator backwards with alternating signs: so that L * (y) = y (n) − (a n−1 y) (n−1) + . . . + (−1) n a 0 y The operator P is called essentially self-adjoint, if there exists a function α = 0 (in an extension of Q(x)), such that Lα = αL * .
It is easy to see that any such α has to satisfy the first order differential equation This essentially adjointness is equivalent to the existence of an invariant pairing on the solution space, symmetric if n is odd and alternating if n is even, [15], [16]. The condition of essential self-adjointness can be expressed by the vanishing of certain differential polynomials in the a i described in [5], [15], [16], [6], called the Calabi-Yau condition.
For an operator of order two the essentially self-adjontness does not impose any conditions, whereas for an differential equation of order three of the form y ′′′ + a 2 y ′′ + a 1 y ′ + a 0 y the Calabi-Yau condition comes down to the vanishing of the quantity whereas for a fourth order operator y (iv) + a 3 y ′′′ + a 2 y ′′ + a 1 y ′ + a 0 y = 0 the quantity has to vanish. In that case the associated differential equation satisified by the 2 × 2 Wronskians reduces its order from 6 to 5.
For operators of higher order, one finds that more than one differential polynomial condition has to be satisfied; here we will be mainly concerned with operators of order four.
It is rather easy to fulfill the condition Q = 0. For example, there exists one particularly nice family of operators of degree n, that we call the main family.

Proposition
If P, Q, R are polynomials of degree n, then the operator has degree n and satisfies the Calabi-Yau condition. If the constant term of P is 1 and of Q and R is zero, then the exponents of L at 0 are all zero. If the roots of P are all distinct, L has n singular points with exponents 0, 1, 1, 2 and a further singular point at infinity. (Here P ′ := θ(P) = x∂P/∂x, etc.) In this paper we take a closer look at differential operators of the form If f = 0, then the numbers α, β, γ, δ are then the characteristic exponents of the operator at ∞.
For this operator the quantity Q is a rational function with as denominator; the singular points of the operator are the roots of this polynomial, together with the point ∞.
The numerator is a polynomial Q of degree 5 in the variable x, so we can write where the coefficients Q i are rather complicated polynomials in the parameters a, b, c, d, e, f and the exponents α, β, γ, δ: So the Calabi-Yau condition is equivalent to six polynomial equations which define a certain affine algebraic set We will give a complete description of the set X, in particular we determine the irreducible components of X. As we have six equations in 10 variables, a first trivial remark is that all irreducible components of X have dimension at least four. We will give explicit parametrisations of all the irreducible components. Note also that the polynomials Q i do not depend on the accessory parameter e, so the solution set X contains the e-line as trivial factor.

Theorem
The algebraic set X defined by the condition Q(a, b, c, d, e, f , α, β, γ, δ) = 0 is the union of ten irreducible components, of which there are seven with f = 0.
Note that if f = 0 we are dealing with an operator of lower degree, so will from now on assume that f = 0 and we will ignore the three components with f = 0. The seven remaining components differ in the behaviour of the average exponents at infinity: There is a single component for which σ is not fixed, that we will call the transverse component. On each of the remaining six components σ has a fixed value. The dimensions and values of σ that appear are summarised in the following diagram. The two five dimensional components are related by a simple transformation, as are the three four dimensional components.
All Calabi-Yau operators known to us, belong to the six-dimensional main component M: or the transverse component T that we describe later.
We remark furthermore that the ideal is, from an algebraic point of view, rather complicated. It is not radical and the components appear with non-trivial multiplicities in the primary decomposition.
For third order operators of degree two the Calabi-Yau condition W = 0 leads in a similar way to an algebraic set Y ⊂ C 8 defined by a system of polynomial equations where now the polynomials are in eight variables. The analysis of this ideal is very similar to the ideal for the fourth order operators, but hardly simpler. The result is also very similar to the case of order four:

Theorem
The algebraic set Y ⊂ C 8 defined by the condition W(a, b, c, d, f , α, β, γ) = 0 is the union of ten irreducible components, of which there are seven with f = 0.
The seven remaining components differ in the behaviour of the (doubled) average exponent Again there is a single transverse component on which the value of σ is not fixed. On each of the remaining six components σ has a fixed value. The dimensions and values of σ that appear are summarised in the following diagram.
The main component (m) is the one with σ = 2, dim = 4 and is given by whereas on the transverse component (t) we have the operator Again, the operators on the two three dimensional components are related by a simple transformation, as are the three two dimensional components. One may speculate that the above patterns extends to higher order operators, but we will not try to pursue it here.

Proof of the theorem
The proof of the above theorem is purely computational and consists of systematically solving of the equations using a computer algebra system like Maple. The miracle here is that this procedure actually works. We will describe the most important steps in the process.

Preliminary reductions
The value of b: The polynomial Q 0 is found to be −4b − 8d + 8c and this allows us directly to eliminate the variable b and thus we will put from now on Symmetry of the exponents: The polynomial Q 5 factors nicely as As we assume that f = 0, this only vanishes if one of the three other factors vanish, meaning that the exponents come in two pairs that add up to the same value that we shall call σ. So without loss of generality we can assume to have the pairs α, σ − α and β, σ − β and we will take the operator to be of the form with b = 2(c − d) as our starting point.
After solving the equations Q 0 = 0, Q 5 = 0, one is left to solve the polynomial equations For sake of concreteness, we give the explicit form of the polynomials: +σ β 2 a 2 + 6 σ cd + 3 σ ad − β 2 ac + β 2 ad − α 2 ac + α 2 ad + σ 2 ac+ −σ 2 ad − 3 σ d 2 − 3 σ c 2 + α σ a 2 + σ β a 2 − σ 2 β a 2 + α 2 σ a 2 , Somewhat to our surprise, it turns out to be rather easy to give a complete solution to these equations. It will be useful to introduce the parameters and use these instead of α and β. One has We will see that most expressions only depend on

The polynomial Q 4
The polynomial Q 4 is an expression that is f 2 times an expression linear in the variable a. As we will only consider cases with f = 0, we can solve for a and find Note that this value for a introduces the denominator So if we work further with this operator, we are implicitly assuming that these two factors do not vanish. If one of them does vanish, we should go back to the previous step, impose these conditions and compute further. For now, we will assume that these factors are non-zero and come back to these cases later.

The case σ = 2
It appears that with this value of a the quantity Q becomes divisible by σ − 2. So in this remarkable σ = 2 case we obtain an operator family with Q = 0. The value of a simplifies to and one obtains what we call the main component: depending c, d, e, f and the exponents α and β as free parameters. Almost all operators of degree two from the AESZ-list are of this type. We obtain our first component of dimension 6.

2.4
The case σ = 2 If however σ = 2 we have to analyse further. Still using the above value for a, the coefficient of x of Q is a complicated expression in c, d, e, f and the exponents, that however is factored by Maple instantly into four factors. First, there are the factors We assume first that these are non-zero, so one is forced put to zero the fourth factor, which leads to a specific value for f : Here a new factor 2σ + ∆ is introduced in the denominator, so we have to come back later to the case 2σ + ∆ = 0.

2.4.1
The subcase σ = 3 Using the above values for a and f , the quantity Q only contains terms with x 2 and x 3 . It factors out the factor σ − 3, which leads to another remarkable operator family with Q = 0: It has c, d, e and α, β as free parameters. This makes up a second component, of dimension 5.

The subcase σ = 3
Still with the given values for a and f , if σ = 3, one has to analyse the quantity Q further. Maple manages to factor the coefficient of x 2 into a product of eight factors, of only three were not supposed to be non-zero at this stage. Each of these factors leads to a linear dependence between c and d.
Upon substitution we obtain an operator that satisfies the condition Q = 0. For this operator one has and has d, e and the exponents α, β, σ as free parameters. So we obtain a further irreducible component of dimension 5. The remarkable thing is that now we have only ∆ apearing in the denominator.
The operator of Bogner mentioned in the introduction is an instance of case A for the parameter choice: Case B: In this case f reduces to 0, so we do not get an operator of degree two.
Case C: In this case the condition Q = 0 reduces to special relations between the exponents. The exceptional cases are: or This concludes the list of all possibilities. We still have to backtrack some of the cases:

Backtracking the remaining cases
Earlier we use the value for a which involved the denominator The introduction of f led to a further denominator 2σ + ∆ The analysis of Case A gave ∆ as denominator, and Case C gave further factor (σ − 1)(σ − 4) + ∆ in the denominator. In each of the cases we have to backtrack and see if we find additional solutions.

2.5.1
The case σ = 1 Looking at the operator with σ = 1, we find from the coefficient of x of Q the value of f . Using this value factoring the coefficient of x 2 of Q we find there is a unique value of c: Here a denominator ∆ appears, so we assume this to be = 0; but the case ∆ = 0 has to be analysed anyway.

2.5.2
The case ∆ = 0 To analyse the further cases, it is useful to write the operator in terms of A and B instead of α and β, so that ∆ = C can be implemented by putting B = C − A. We start with C = 0. The following seven components make up the intersection with ∆ = 0: This is inside the main component.
This is inside the σ = 3 component.

Remaining cases
There are more cases to check and this could be done along the above sketched lines. Instead of this, we used the computer algebra system Singular. There is a library called primdec.lib containing the algorithms for doing primary decomposition. On our computer, Singular could not straight away compute the primary decomposition of the ideal I defined by the vanishing of the x-coefficients of Q. But after taking the ideal quotient of I by the ideals already found, it turned out to be able to find the last component. In total we thus found seven irreducible components of the variety defined by Q = 0.

Summary of the seven components
In total we have found seven components. Membership to a component is defined by certain polynomial relations between the coefficients. Below we give the defining relations for each of the components, the operator with its coefficients together with its generic Riemann symbol.
The big σ = 2 component, also called the main component is defined by the linear equations For a general member of this component the Riemann symbol is: where * are the two roots of 1 + ax + f x 2 = 0. The operator depends on the six parameters a, d, e, f , α, β. The parameters a, f determine the position of the singular fibres.
The parameters α, β determine the four exponents at ∞ which are symmetric around the value 2. The parameters d, e are accessory parameters. However, if the two roots coincide, a 2 = 4 f , then we obtain an operator with three singular points.
The σ = 1 component is defined by the four conditions: The Riemann symbol of the operator is: where the special points x and y are located at: The point y is not really a singular point of the differential operator, as the exponents 0, 1, 2, 3 are those of a regular point of the differential equation. Indeed, as it can be checked that no logarithms occur.
The σ = 3 component is defined by the conditions: The Riemann symbol of the operator is: where the singular points are located at: The operator has five parameters c, d, e, α, β. α, β are parameters that determine the exponents at infinity, symmetric around 3. The parameters c, d determine, together with α and β the position of the two singularities. e is a single accessory parameter.
The σ = 0 component is defined by the conditions: The generic Riemann symbol for this operator is: where the special points x, y are located at the solutions of 1 + ax + f x 2 = 0. The operator depends on four parameters a, e, f , α.
The small σ = 2 component is defined by the four conditions: The Riemann symbol for a general member of this component is: where the special points are The σ = 4 component is defined by the four conditions: The Riemann symbol of this operator is: where the singular points x, y are solutions to 1 + ax + f x 2 = 0.
The transverse component finally is defined by the equations Solving for a, b, c, f we find the family of operators The Riemann symbol of the generic member in this family is: The big and small σ = 2 components do intersect: The transverse component intersects all other components. No other pair of components can intersect, as they have a different σ-value.

The small components
It turns out that the components with σ = 0, 2, 4 are closely related to each other and the solutions of the corresponding differential equations can be related in a simple manner. The families of operators in question are: Their relation is most easily understood from the perspective of the operator P 2 (a, d, e, A) and its exponents. The discriminant ∆(x) of the operator P 2 (a, d, e, A) factors as (1 + dx/4)(1 + (4a − d)x/4). The exponents at the first factor are −1, 0, 1, 2 and at the second factor 0, 1, 2, 3. The exponents at infinity are of the form Multiplication of a solution to P 2 (a, d, e, A) by (1 + dx/4) will shift the exponents at the first singular point by +1 and the exponents at infinity by −1, which then leads to an operator with two points with exponents 0, 1, 2, 3 and α, −α, iα, −iα as exponents at infinity, i.e with σ = 0 and ∆ = 0. Similarly, division of a solution to P 2 (a, d, e, A) by (1 + (4a − d)x/4) leads to a shift at the second singularity by −1 and hence produces two singularities with exponents −1, 0, 1, 2. At infinity the exponents shift by +1 and are of the form We remark that the exponents at infinity can only be all real in the very special case that these all coincide.
A precise statement relating these operators is the following: It is a remarkable fact that the local monodromies around the special points = 0 and = ∞ are trivial. The monodromy around 0 is MUM, so has infinite order and thus the monodoromy group is the non-reductive group Z, which clearly prevents these operators from being Picard-Fuchs operators of a family of algebraic varieties.
Something similar happens for the σ = 1 and γ = 3 components, given by the operator families One verifies directly that These operators are also remarkable. One of the roots of the polynomial 1 + ax + f x 2 is in fact a regular point, the other is a conifold point, so the three non-trivial monodromies are precisely as those for the hypergeometric operators. In fact, if we multiply the hypergeometric operator with the linear factor 1 + Mx, the result is the operator for the specific value of the accessory parameter e. If e takes on another value, the operator is no longer equal to linear factor × hypergeometric.

The known Calabi-Yau operators of degree 2
Below we will give an overview of the 69 Calabi-Yau operators of degree two that are currently known to us. We also include some further remarkable operators that are not strictly Calabi-Yau in the sense of [5]. Of course, this is just a small portion of the list of the still growing list Calabi-Yau operators that was started in [7] and now contains more then 500 members. An update of this list in preparation [8].
For operators of degree two there are cases with three and with four singular points.
Most of the known operators can be related to hypergeometric operators or operators that are convolutions, that is, can be obtained via Hadamard product of operators of lower order. Recall that the Hadamard product of two power series φ(x) := ∑ n a n x n and ψ(x) := ∑ n b n x n is the series φ * ψ(x) = ∑ n a n b n x n and by a classical theorem, if φ and ψ satisfy a linear differential equations, then so does φ * ψ. The Hadamard product corresponds to the multiplicative convolution of local systems, so are under very good control, see [23], [24].
Frobenius basis: For a fourth order operator P with at MUM-point at 0, there is a canonical basis of solutions y 0 , y 1 , y 2 , y 3 to Py = 0 on a sufficiently small slit disc around the origin, called the Frobenius basis: The Calabi-Yau operators from [5] are characterised by three integrality conditions: 1) integrality of the solution y 0 : 2) integrality of the q-coordinate q(x): 3) integrality of the instanton numbers n d .
These can be computed in several different ways. The second logarithmic derivative of y 2 /y 0 expressed in the q-coordinate is called the Yakawa-coupling K(q) of P: Now expand K(q) as a Lambert-series The numbers n d are called the (normalised, n 0 = 1) instanton numbers of P and are required to to be integral. It is natural to allow small denominators in f 0 , q(x) and n d , so require only N-integrality, for some denominator N.

Description of the operators
Many of the Calabi-Yau operators of degree 2 that we will describe involve special lower order operators of Calabi-Yau type. These are introduced and discussed in Appendix B.

Operators with three singular points
for which the hypergeometric operator, scaled by N, is a Calabi-Yau operator, [1]. Corresponding to these, there are also 14 hypergeometric fifth order Calabi-Yau operator These operators have a Yifan Yang pull-back to 14 special fourth order operators, called the tilde-operators 1, 2, . . . , 14. These operators replace the more complicated hat-operatorŝ i, i = 1, 2, . . . , 14 that appeared in the list [7] into which they can be transformed. However, operator 3, is not new, as it can be reduced to operator 2.33

II. The 16 operators of type H * µ(H)
There are four special hypergeometric second order operators called A, B, C, D. The Möbius transformation that interchanges the singularity with infinity leads to four Möbius transformed hypergeometric operators where e, h, i, j refer to the names given in [5], see also Appendix B.
Thus we can form the 16 Hadamard product A * µ(B), etcetera and obtain 16 fourth order operators of degree two with three singular points. But there are two surprises: the operator C * µ(A) can be reduced to the hypergeometric operator AESZ#3 = 1.3 and the operator C * µ(B) has vanishing instanton numbers and in fact is the third symmetric power of the hypergeometric second order operator So from this we obtain only 14 Calabi-Yau operators of degree two that are 'new'.

III. The four Hadamard products I * µ(H ′ )
There are also four special hypergeometric third order operators A ′ , B ′ , C ′ , D ′ , analogous to second order operators A, B, C, D. Interchanging the singularity with the point at infinity, we obtain third order operators of degree two with three singular points. Taking Hadamard product with the central binomial coefficient produces four fourth order operators of degree two. One of these, I * β, can be transformd by a quadratic transformation to an operator (I * β) * which in turn can be transformed to the hypergeometric operator AESZ#3. The operator I * θ can be transformed in 2.17, so only two of these give rise to new Calabi-Yau operators.

IV. Four sporadic operators
The operator 2.63 = #AESZ84 can not be constructed in this way. Furthermore, there are three operators of Bogner-type that were obtained from transformations of the higher degree operators AESZ#245, AESZ#406 and AESZ#255. As these operators have unusual properties, they probably should better be understood in terms of the higher degree versions. We remark that the local monodromies around the two singularities = 0 are of finite order. So we are dealing here with special Heun operators. They appear in the work of Beukers [12], and Zagier [42] and are related to the rational elliptic surfaces described by Beauville [11], so we call them the BZB-operators. They were denoted a, b, c, d, f , g in [5].

Operators with four singular points
Taking the Hadamard product with the hypergeometric operators A, B, C, D produces 24 fourth order operators of degree two with four singular points.

II. The six Hadamard products I * BZB ′
The BZB-operators have also analogs as third order operators that we will call the BZB'operators, that were called α, γ, δ, ǫ, ζ, η in [5]. Taking I * with these operators produce six new fourth order operators of degree two with four singularities.
III. Six Sporadic operators with four singularities III.a. Three operators I * sporadic third order: There are three sporadic third order operators of degree two, described in Appendix B. These lead to fourth order operators of degree two by taking I * with them.
III.b: Three original operators: There are three operators that apparently can not be constructed starting from simpler operators of lower order; we call them original operators.

The count
So in total we have described Calabi-Yau operators of degree two with three singular points and with four singular points, making up a total of 69. Because these all have different instanton numbers, these are essentially different and there are no algebraic transformations mapping one to any of the other operators.

Miscellaneous operators
There are two further types of operators of degree two that deserve to be mentioned at this place, as they challenge the precise definition of Calabi-Yau operator, as formulated in [5].
Reducible operators: A large scale Zagier-type search for degree two operators with integral solutions has been perfomed by Pavel Metelitsyn. Apart from the operators equivalent to the ones mentioned above, there are some interesting reducible fourth order operators. With regards to the q-coordinate and instanton numbers these examples behave very much in the same way as the irreducible ones. An example is the operator One finds: A Strange Operator: M. Bogner [15] also found an example of an operator of degree two with an integer solution, integral q-coordinate, integral K(q)-series, but for which the instanton numbers are badly non-integral. We start from the third order operator The operator I * L is the fourth order operator We find the following nice integral series: However, if we expand K(q) as a Lambert-series we do not get integral instanton numbers: This is rather puzzling. Conjecturally, already the fact that φ(x) is an integral series implies that the operator is a factor of a Picard-Fuchs equation of a family of varieties defined over Q and indeed our operator is of geometrical origin. It was conjectured in [5] that the integrality of the q-series implies the integrality of the instanton numbers, which in this example is not the case. The integrality of solution and mirror map clearly indicate that we have a rank four Calabi-Yau motive and one would expect the general arguments for the integrality of [41] to be applicable, but apparently they are not. There might exist a different scaling of the coordinate that repairs this defect, but up to now we have been unable to find it.
Acknowledgement This work for this paper was done in 2016 and we had hoped to find many new degree two operators. That did not happen and we never found the time to write the paper. We intended it to be a possible reference for this class of operators. The untimely death of Gert left the paper in a unfinished state and it is only now I found the time to complete it. We thank M. Bogner, J. Hofmann and K. Samol for help with computations, done years ago, and W. Zudilin and O. Gorodetsky for useful comments, corrections and showing interest in this paper.

Appendix A
The following table lists the operators as they appear in the electronic database accessable at https://cydb.mathematik.uni-mainz.de/. (These are ordered by degree, so that the nth operator of degree d gets the code d.n in that list.) The second column lists the AESZ-numbers as appearing in [7]. The last three list the first three instanton numbers (normalised, n 0 = 1). As these do not change under transformations of the coordinate, the pair of numbers |n 1 |, |n 3 | provide a simple 'fingerprint/superseeker' to identify an operator. In Appendix C one finds some monodromy and modular form information for these operators.

Appendix B
In this appendix we collect some information on operators of lower order that are relevant to the construction of the order four Calabi-Yau operators of degree two.

First order operator
There is a single first order operator that we call I: 2n n x n is the unique holomorphic solution of I, and in fact coincides for x ≤ 1/4 with the algebraic function 1 √ 1 − 4x .

Second order operators The four hypergeometric cases
There are four very remarkable hypergeometric second order operators of the form with Riemann symbol of the form These four operators appear in many different areas of mathematics, like the alternative theories of Ramanujan and are described at many places, see e.g. [21].
Here MW stands for the Mordell-Weil group of sections of the fibration. The matrices A, B and (BA) −1 are monodromy matrices for loops around 0, c = 0, ∞ and ∞ respectively. In the monodromy column, the group generated by A and B is identified. The column H/MP denotes the notation used in [35] and [32]. The last column denotes the type of the Picard-Fuchs operator. We see that there are three surfaces that belong to operator A, but the surface I 1 I 4 I * 1 shows it is a member of a very regular series of four surfaces, which we will say to belong to the operators A, B, C, D. The surfaces belonging to the operators A, B, C, D have a clear relation to the geometry of the simple elliptic singularities. The fibres over ∞ correspond to configurations of rational curves intersecting in the pattern of the affine Dynkin diagrams. The corresponding simple elliptic singularities have a single modulus. These families of elliptic curves lead to elliptic surfaces that cover the corresponding extremal elliptic surfaces.

Case Dynkin Equation Cover
A D 5 The A case is slightly different, as the corresponding singularity is not a hypersurface, but a complete intersection of two quadrics. The last three cases correspond to the euclidian triples (p, q, r) of integers ≥ 2 with the property that

The transformed hypergeometric cases
The Möbius-transformation x → −

The six Beukers-Zagier-Beauville operators
There are six special second order operators of degree two with four singular points of the form We refer to [28] for a representation of these sequences as constant term of a Laurent polynomial.
Unfortunately there is no natural naming for these operators. In the first column of the table we used the "AZ-names" used in [5]. In [13] and [42] the same operators are named A, B, . . . , F (Not to be confused with the hypergeometric operators A, B, C, D mentioned above!); these names appear in the second column as "BZ-names". These operators also appear as Picard-Fuchs operators for rational elliptic surfaces, namely the six semistable families of elliptic curves with four exceptional fibres, that have been studied by Beauville, [11]. The modular level gives a more or less natural way order the operators.
fibres j Pencil Group Operator The coefficients of the holomorphic solution of the operator b are the Apéry-numbers for the irrationality of π 2 . It belongs to the elliptic surface with fibres I 5 , I 5 , I 1 , I 1 over the modular curve X 1 (5), as was first noticed by Beukers in [12].
The operators a, c, g differ by a Möbius-transformations. The base of this family can be identified with the modular curve for the congruence subgroup Γ 0 (6). Over that curve we have the semi-stable elliptic surface with fibres of type I 6 , I 3 , I 2 , I 1 . The operator d belongs to the two isogenous elliptic surfaces with fibres I 8 , I 2 , I 1 , I 1 and I 4 , I 4 , I 2 , I 2 .
Operator f belongs to the two isogenous elliptic surfaces with fibres I 9 , I 1 , I 1 , I 1 and I 3 , I 3 , I 3 , I 3 which is nothing but the Hesse pencil.

Third order operators Hypergeometric operators
Closely related to the operators A, B, C, D there are also four very remarkable hypergeometric third order operators Also important are the symmetric squares of the operators A, B, C, D. These are of the form θ 3 − x(2θ + 1)(aθ 2 + aθ + b) + a 2 x 2 (θ + α)(θ + 1)(θ + β), α + β = 2 and have Riemann-symbol of the form

The four transformed symmetric squares
More important for us are the transformed symmetric squares of the operators A, B, C, D.
These are of the form and have Riemann-symbol of the form The names β, ι, θ, κ stem from [5].

The six BZB' -operators
Corresponding to the six BZB-operators, there are also six special third order operators of degree two with four singular points of the form and hence have a Riemann-symbol of the form a 1 , a 2 , a 3 , . . . a n δ 7 3 81 1, 3, 9, 3, −279 For a Laurent-representation of these sequences we refer to [28].
The correspondence between the BZB-operators and BZB ′ -operators is determined by the following table: There is a general formula expressing the relation between the BZB and BZB' operators: which shows that the third order operators are twisted versions of the symmetric squares of the corresponding second order operators. The factor Q(x) is the discriminant of the second order operator with y 0 as solution.

The three sporadic third order operators
We know of three very special third order operators that seem to be not directly to the analogous of the second order operators. These are the following ones: This example has an A-incarnation as K3-surface that is the intersection of four hyperplane sections of type (1, 1) in P 3 × P 3 . If we take I * , we obtain a fourth order operator 2.60 with A-incarnation a CY threefold that appears as intersection of two (1, 1) and one (2, 2) hypersurface in P 3 × P 3 , [9].
Sporadic 2: This example has an A-incarnation as K3-surface that is the intersection of six hyperplane sections of the Grassmanian G(2, 6) in its Plücker embedding. If we take I * , we obtain a fourth order operator 2.61 with A-incarnation a CY threefold that appears as intersection of four linear and a quadratic hypersurface in Grassmanian G(2, 6) in its Plücker embedding, [10].
Sporadic 3: These sporadic operators and sequences were also found by S. Cooper [20], where they are called s 10 , s 7 and s 18 and where the explicit description for a n in the last case also can be found. We refer to [28] for representations as constant term sequence of a Laurent polynomial.

Miscellaneous third order operators
There is a remarkable infinite series of third order operators with an integral solution.

Theorem
The differential operator as solution.
Proof: The Riemann symbol of the operator P is The operator is the symmetric square of a second order operator with Rieman symbol Multiplication by a factor (1 − x) k/2 gives a hypergeometric operator with Riemann symbol If we compare with the Riemann symbol of the standard hypergeometric operator for

Corollary
If k and d are rational numbers, then the only the primes appearing in denominator of k and d appear in the denominator of the coefficients of the holomorphic solution of P. After an appropriate rescaling, the coefficients become integral.

Appendix C
In this Appendix we list all essentially distinct Calabi-Yau operators of order four and degree two that are known to us. We also present some further information. To present the monodromy, it is sometimes convenient to use u i (x) = y i (x)/(2πi) 3 which we call the scaled Frobenius basis. The monodromy transformation around 0 in this basis is given by the matrix     1 0 0 0 In the cases with four singular points, the operator belongs to the main component and the Riemann symbol of the operator has the form where c 1 , c 2 are solutions to the equation We will always assume that the exponents α, β are between 0 and 1. For each operator we give the monodromy transformations around the other singular points the with respect to this basis. This information was computed by J. Hofmann, [29]. A very common transformation is the symplectic reflection v → v− < v, w > w in a vector w = (a, b, c, d) represented by the matrix Such a tranformation is found at conifold points, where the exponents are 0 1 1 2. For the conifold nearest to the origin the vector is of the form (a, 0, c, d) and in that case the characteristic numbers h 3 , c 2 h, c 3 = χ of the operator are determined by the reflection vector via These numbers are formally attached to the differential operator, but have the interpretation as characteristic numbers of the mirror manifold, see [39].

Modular form information.
For conifold points appearing at rational values of the parameters one can determine the coeffients of a weight four modular form as described in the thesis of K. Samol, [33], [34]. We include here the name of modular form from S 4 (Γ 0 (N) as it appears in the list of C. Meyer, [31]. For example, (6/1) denotes the first (and only) weight four newform of level 6, which is the η-product For critical points at real quadratic irrationalities one expects the appearance of Hilbert modular forms; at imaginary quadratic irrationalities one would expect Bianchi modular forms. Another phenomenon that may occur is the of a so-called K-point, see [39]. Here there are two single logarithms appearing (two size two Jordan blocks), and one expects the appearance of a modular form from S 3 (Γ 0 (N)).
Below we list the 90-degree rotated extended Riemann-symbols; each row starts with the singular value, followed by the four exponents at that point, then monodromy information and finally modular form information. At the MUM point the monodromy is always standard, and we give the first few instanton numbers as a substitute for modular form information. If Hilbert-modular, Bianchi or weight three modular forms are expected, we indicate this with h, b or k. In case the monodromy is of finite order we put −. The operators with three singular points which are on the main components can all be written in the form There is a single accessory parameter e.

2.63: #84
for which the hypergeometric operator, scaled by N, is a Calabi-Yau operator, [1]. Corresponding to these, there are also 14 hypergeometric fifth order Calabi-Yau operator These operators have a Yifan Yang pull-back to 14 special fourth order operators, called the tilde-operators 1, 2, . . . , 14. These operators replace the more complicated hat-operatorŝ i, i = 1, 2, . . . , 14 that appeared in the list [7] into which they can be transformed.