Determination of normalized extremal quasimodular forms of depth 1 with integral Fourier coefficients

The main purpose of this paper is to determine all normalized extremal quasimodular forms of depth 1 whose Fourier coefficients are integers. By changing the local parameter at infinity from $q=e^{2\pi i \tau}$ to the reciprocal of the elliptic modular $j$-function, we prove that all normalized extremal quasimodular forms of depth 1 have a hypergeometric series expression and that integrality is not affected by this change of parameters. Furthermore, by transforming these hypergeometric series expressions into a certain manageable form related to the Atkin(-like) polynomials and using the lemmas that appeared in the study of $p$-adic hypergeometric series by Dwork and Zudilin, the integrality problem can be reduced to the fact that a polynomial vanishes modulo a prime power, which we prove. We also prove that all extremal quasimodular forms of depth 1 with appropriate weight-dependent leading coefficients have integral Fourier coefficients by focusing on the hypergeometric expression of them.