Summary: We show that almost all the linear differential operators factors obtained in the analysis of the \(n\)-particle contributions \(\widetilde {\chi}^{(n)}\)’s of the susceptibility of the Ising model for \(n \leq 6\) are linear differential operators associated with elliptic curves. Beyond the simplest differential operators factors which are homomorphic to symmetric powers of the second order operator associated with the complete elliptic integral \(E\), the second and third order differential operators \(Z_{2}, F_{2}, F_{3}, \widetilde L_3\) can actually be interpreted as modular forms of the elliptic curve of the Ising model. A last order-4 globally nilpotent linear differential operator is not reducible to this elliptic curve, modular form scheme. This operator is shown to actually correspond to a natural generalization of this elliptic curve, modular form scheme, with the emergence of a Calabi-Yau equation, corresponding to a selected \(_{4}F_{3}\) hypergeometric function. This hypergeometric function can also be seen as a Hadamard product of the complete elliptic integral \(K\), with a remarkably simple algebraic pull-back (square root extension), the corresponding Calabi-Yau fourth order differential operator having a symplectic differential Galois group \(SP(4,\mathbb C)\). The mirror maps and higher order Schwarzian ODEs, associated with this Calabi-Yau ODE, present all the nice physical and mathematical ingredients we had with elliptic curves and modular forms, in particular an exact (isogenies) representation of the generators of the renormalization group, extending the modular group \(SL(2,\mathbb Z)\) to a \(GL(2,\mathbb Z)\) symmetry group.