j-invariant

In mathematics, Klein's j-invariant, regarded as a function of a complex variable τ, is a modular function of weight zero for SL(2, Z) defined on the upper half-plane of complex numbers. It is the unique such function which is holomorphic away from a simple pole at the cusp such that

\( j\left(e^{\frac{2}{3}\pi i}\right) = 0, \quad j(i) = 1728. \)

Klein's j-invariant in the complex plane

Rational functions of j are modular, and in fact give all modular functions. Classically, the j-invariant was studied as a parameterization of elliptic curves over C, but it also has surprising connections to the symmetries of the Monster group (this connection is referred to as monstrous moonshine).

Definition

Real part of the j-invariant as a function of the nome q on the unit disk

Phase of the j-invariant as a function of the nome q on the unit disk

Further information: Elliptic curve § Elliptic curves over the complex numbers and Modular forms

While the j-invariant can be defined purely in terms of certain infinite sums (see g2, g3 below), these can be motivated by considering isomorphism classes of elliptic curves. Every elliptic curve E over C is a complex torus, and thus can be identified with a rank 2 lattice; i.e., two-dimensional lattice of C. This is done by identifying opposite edges of each parallelogram in the lattice. It turns out that multiplying the lattice by complex numbers, which corresponds to rotating and scaling the lattice, preserves the isomorphism class of the elliptic curve, and thus we can consider the lattice generated by 1 and some τ in H (where H is the Upper half-plane). Conversely, if we define

\( \begin{align} g_2 &= 60\sum_{(m,n) \neq (0,0)} (m + n\tau)^{-4}, \\ g_3 &= 140\sum_{(m,n) \neq (0,0)} (m + n\tau)^{-6}, \end{align} \)

then this lattice corresponds to the elliptic curve over C defined by y2 = 4x3 − g2x - g3 via the Weierstrass elliptic functions. Then the j-invariant is defined as

\( j(\tau) = 1728 \frac{g_2^3}{\Delta} \)

where the modular discriminant Δ is

\( \Delta = g_2^3 - 27g_3^2 \)

It can be shown that Δ is a modular form of weight twelve, and g2 one of weight four, so that its third power is also of weight twelve. Thus their quotient, and therefore j, is a modular function of weight zero, in particular a meromorphic function H → C invariant under the action of SL(2, Z). As explained below, j is surjective, which means that it gives a bijection between isomorphism classes of elliptic curves over C and the complex numbers.
The fundamental region
The fundamental domain of the modular group acting on the upper half plane.

The two transformations τ → τ + 1 and τ → τ−1 together generate a group called the modular group, which we may identify with the projective special linear group PSL(2, Z). By a suitable choice of transformation belonging to this group,

\( \tau \mapsto \frac{a\tau + b}{c\tau +d}, \qquad ad-bc =1, \)

we may reduce τ to a value giving the same value for j, and lying in the fundamental region for j, which consists of values for τ satisfying the conditions

\( \begin{align} |\tau| &\ge 1 \\ -\tfrac{1}{2} &< \mathfrak{R}(\tau) \le \tfrac{1}{2} \\ -\tfrac{1}{2} &< \mathfrak{R}(\tau) < 0 \Rightarrow |\tau| > 1 \end{align} \)

The function j(τ) when restricted to this region still takes on every value in the complex numbers C exactly once. In other words, for every c in C, there is a unique τ in the fundamental region such that c = j(τ). Thus, j has the property of mapping the fundamental region to the entire complex plane.

As a Riemann surface, the fundamental region has genus 0, and every (level one) modular function is a rational function in j; and, conversely, every rational function in j is a modular function. In other words the field of modular functions is C(j).

Class field theory and j

The j-invariant has many remarkable properties:

If τ is any of the singular moduli, that is, any element of an imaginary quadratic field with positive imaginary part (so that j is defined) then j(τ) is an algebraic integer.[1]

The field extension Q[j(τ), τ]/Q(τ) is abelian, that is, it has an abelian Galois group.

Let Λ be the lattice in C generated by {1, τ}. It is easy to see that all of the elements of Q(τ) which fix Λ under multiplication form a ring with units, called an order. The other lattices with generators {1, τ′}, associated in like manner to the same order define the algebraic conjugates j(τ′) of j(τ) over Q(τ). Ordered by inclusion, the unique maximal order in Q(τ) is the ring of algebraic integers of Q(τ), and values of τ having it as its associated order lead to unramified extensions of Q(τ).

These classical results are the starting point for the theory of complex multiplication.
Transcendence properties

In 1937 Theodor Schneider proved the aforementioned result that if τ is a quadratic irrational number in the upper half plane then j(τ) is an algebraic integer. In addition he proved that if τ is an algebraic number but not imaginary quadratic then j(τ) is transcendental.

The j function has numerous other transcendental properties. Kurt Mahler conjectured a particular transcendence result that is often referred to as Mahler's conjecture, though it was proved as a corollary of results by Yu. V. Nesterenko and Patrice Phillipon in the 1990s. Mahler's conjecture was that if τ was in the upper half plane then exp(2πiτ) and j(τ) were never both simultaneously algebraic. Stronger results are now known, for example if exp(2πiτ) is algebraic then the following three numbers are algebraically independent, and thus at least two of them transcendental:

\( j(\tau), \frac{j^\prime(\tau)}{\pi}, \frac{j^{\prime\prime}(\tau)}{\pi^2} \)

The q-expansion and moonshine

Several remarkable properties of j have to do with its q-expansion (Fourier series expansion), written as a Laurent series in terms of q = exp(2πiτ), which begins:

\( j(\tau) = {1 \over q} + 744 + 196884 q + 21493760 q^2 + 864299970 q^3 + 20245856256 q^4 + \cdots \)

Note that j has a simple pole at the cusp, so its q-expansion has no terms below q−1.

All the Fourier coefficients are integers, which results in several almost integers, notably Ramanujan's constant:

\( e^{\pi \sqrt{163}} \approx 640320^3 + 744. \)

The asymptotic formula for the coefficient of qn is given by

\( \frac{e^{4\pi\sqrt{n}}}{\sqrt{2}n^{3/4}}, \)

as can be proved by Hardy–Littlewood circle method.[2][3]

Moonshine

More remarkably, the Fourier coefficients for the positive exponents of q are the dimensions of the graded part of an infinite-dimensional graded algebra representation of the monster group called the moonshine module – specifically, the coefficient of qn is the dimension of grade-n part of the moonshine module, the first example being the Griess algebra, which has dimension 196,884, corresponding to the term 196884q. This startling observation was the starting point for moonshine theory.

The study of the Moonshine conjecture led J.H. Conway and Simon P. Norton to look at the genus-zero modular functions. If they are normalized to have the form

\( q^{-1} + {O}(q) \)

then Thompson showed that there are only a finite number of such functions (of some finite level), and Cummins later showed that there are exactly 6486 of them, 616 of which have integral coefficients.[4]
Alternate Expressions

We have

\( j(\tau) = \frac{256(1-x)^3}{x^2} \)

where x = λ(1−λ) and λ is the modular lambda function

\( \lambda(\tau) = \frac{\theta_2^4(0,\tau)}{\theta_3^4(0,\tau)} = k^2(\tau) \)

a ratio of Jacobi theta functions \theta_{m}, and is the square of the elliptic modulus k(\tau).[5] The value of j is unchanged when λ is replaced by any of the six values of the cross-ratio:[6]

\( \left\lbrace { \lambda, \frac{1}{1-\lambda}, \frac{\lambda-1}{\lambda}, \frac{1}{\lambda}, \frac{\lambda}{\lambda-1}, 1-\lambda } \right\rbrace \)

The branch points of j are at {0, 1, ∞}, so that j is a Belyi function.[7]
Expressions in terms of theta functions

Define the nome \( q=e^{\pi i \tau} and the Jacobi theta function,

\( \vartheta(0; \tau) = \vartheta_{00}(0; \tau) = 1 + 2 \sum_{n=1}^\infty \left(e^{\pi i\tau}\right)^{n^2} = \sum_{n=-\infty}^\infty q^{n^2} \)

from which one can derive the auxiliary theta functions. Let,

\( \begin{align} a &= \theta_{2}(0; q) = \vartheta_{10}(0; \tau) \\ b &= \theta_{3}(0; q) = \vartheta_{00}(0; \tau) \\ c &= \theta_{4}(0; q) = \vartheta_{01}(0; \tau) \end{align} \)

where \( \theta_{m} \) and \(\vartheta_{n} \) are alternative notations, and \( a^4 - b^4 + c^4 = 0 \). Then,

\( \begin{align} g_2(\tau) &= \tfrac{2}{3}\pi^4 \left(a^8 + b^8 + c^8\right) \\ g_3(\tau) &= \tfrac{4}{27}\pi^6 \sqrt{\frac{(a^8+b^8+c^8)^3-54(abc)^8}{2}} \\ \Delta &= g_2^3-27g_3^2 = (2\pi)^{12} \left(\tfrac{1}{2}a b c\right)^8 = (2\pi)^{12}\eta(\tau)^{24} \end{align} \)

for Weierstrass invariants g2, g3, and Dedekind eta function η(τ). We can then express j(τ) in a form which can rapidly be computed.

\( j(\tau) = 1728\frac{g_2^3}{g_2^3-27g_3^2} = 32 {(a^8 + b^8 + c^8)^3 \over (a b c)^8}

Algebraic definition

So far we have been considering j as a function of a complex variable. However, as an invariant for isomorphism classes of elliptic curves, it can be defined purely algebraically. Let

\( y^2 + a_1 xy + a_3 y = x^3 + a_2 x^2 + a_4 x + a_6 \)

be a plane elliptic curve over any field. Then we may define

\( b_2 = a_1^2 + 4a_2,\quad b_4 = a_1a_3 + 2a_4 \)
\( b_6 = a_3^2 + 4a_6,\quad b_8 = a_1^2a_6 - a_1a_3a_4 + a_2a_3^2 + 4a_2a_6 - a_4^2 \)
\( c_4 = b_2^2 - 24b_4,\quad c_6 = -b_2^3 + 36b_2b_4 - 216b_6 \)

and

\( \Delta = -b_2^2b_8 + 9b_2b_4b_6 - 8b_4^3 - 27b_6^2 \)

the latter expression is the discriminant of the curve.

The j-invariant for the elliptic curve may now be defined as

\( j = {c_4^3 \over \Delta} \)

In the case that the field over which the curve is defined has characteristic different from 2 or 3, this definition can also be written as

\( j = 1728{c_4^3 \over c_4^3-c_6^2}

Inverse function

The inverse function of the j-invariant can be expressed in terms of the hypergeometric function 2F1 (see also the article Picard–Fuchs equation). Explicitly, given a number N, to solve the equation j(τ) = N for τ can be done in at least four ways.

Method 1: Solving the sextic in λ,

\( j(\tau) = \frac{256(1-\lambda(1-\lambda))^3}{(\lambda(1-\lambda))^2} \)

where λ is the modular lambda function. Let x = λ(1−λ) and the sextic can be solved as a cubic in x. Then,

\( \tau = i \ \frac{{}_2F_1 \left (\tfrac{1}{2},\tfrac{1}{2},1,1 - \lambda \right )}{{}_2F_1 \left (\tfrac{1}{2},\tfrac{1}{2},1,\lambda \right)} \)

for any of the six values of λ.

Method 2: Solving the quartic in γ,

\( j(\tau) = \frac{27(1 + 8\gamma)^3}{\gamma(1 - \gamma)^3} \)

then for any of the four roots,

\( \tau = \frac{i}{\sqrt{3}} \frac{{}_2F_1 \left (\tfrac{1}{3},\tfrac{2}{3},1,1-\gamma \right)}{{}_2F_1 \left(\tfrac{1}{3},\tfrac{2}{3},1,\gamma \right )} \)

Method 3: Solving the cubic in β,

\(j(\tau) = \frac{64(1+3\beta)^3}{\beta(1-\beta)^2} \)

then for any of the three roots,

\( \tau = \frac{i}{\sqrt{2}} \frac{{}_2F_1 \left (\tfrac{1}{4},\tfrac{3}{4},1,1-\beta \right)}{{}_2F_1 \left(\tfrac{1}{4},\tfrac{3}{4},1,\beta \right )} \)

Method 4: Solving the quadratic in α,

\( j(\tau)=\frac{1728}{4\alpha(1-\alpha)} \)

then,

\( \tau = i \ \frac{{}_2F_1 \left (\tfrac{1}{6},\tfrac{5}{6},1,1-\alpha \right)}{{}_2F_1 \left(\tfrac{1}{6},\tfrac{5}{6},1,\alpha \right )} \)

One root gives τ, and the other gives 1/τ, but since j(τ) = j(1/τ), then it doesn't make a difference which α is chosen. The latter three methods can be found in Ramanujan's theory of elliptic functions to alternative bases.

The inversion is highly relevant to applications via enabling high-precision calculations of elliptic functions periods even as their ratios become unbounded. A related result is the expressibility via quadratic radicals of the values of j at the points of the imaginary axis whose magnitudes are powers of 2 (thus permitting compass and straightedge constructions). The latter result is hardly evident since the modular equation of level 2 is cubic.
Pi formulas

The Chudnovsky brothers found in 1987,[8]

\( \frac{1}{\pi} = \frac{12}{640320^{3/2}} \sum_{k=0}^\infty \frac{(6k)! (163 \cdot 3344418k + 13591409)}{(3k)!(k!)^3 (-640320)^{3k}} \)

which uses the fact that \( j\big(\tfrac{1+\sqrt{-163}}{2}\big) = -640320^3 \). For similar formulas, see the Ramanujan–Sato series.
Special values

The j-invariant vanishes at the "corner" of the fundamental domain at

\( \tfrac{1}{2}\left(1 + i \sqrt{3}\right). \)

Here are a few more special values (only the first four of which are well known; in what follows, j means J/1728 throughout):

\( \begin{align} j(i) &= j \left( \tfrac{1 + i}{2} \right) = 1 \\ j\left(\sqrt{2}i\right) &= \big(\tfrac{5}{3}\big)^3 \\ j(2i) &= \big(\tfrac{11}{2}\big)^3 \\ j\left(2\sqrt{2}i\right) &= \tfrac{125}{216} \left(19 + 13\sqrt{2} \right)^3\\ j(4i) &= \tfrac{1}{64} \left(724 + 513\sqrt{2} \right)^3\\ j\left( \tfrac{1 + 2i}{2} \right) &= \tfrac{1}{64} \left(724 - 513\sqrt{2} \right)^3\\ j\left( \tfrac{1 + 2\sqrt{2}i}{3} \right) &= \tfrac{125}{216} \left(19 - 13\sqrt{2} \right)^3\\ j(3i) &= \tfrac{1}{27} \left(2 + \sqrt{3}\right)^2 \left(21 + 20\sqrt{3}\right )^3 \\ j\left(2\sqrt{3}i\right) &= \tfrac{125}{16} \left(30 + 17\sqrt{3}\right)^3\\ j\left( \tfrac{1 + 7\sqrt{3}i}{2} \right) &= -\tfrac{64000}{7} \left(651 + 142\sqrt{21} \right)^3\\ j\left(\tfrac{1 + 3\sqrt{11}i}{10} \right) &= \tfrac{64}{27} \left(23 - 4\sqrt{33}\right)^2 \left(-77 + 15\sqrt{33} \right)^3\\ j\left(\sqrt{21}i\right) &= \tfrac{1}{32} \left(5 + 3\sqrt{3}\right)^2 \left(3 + \sqrt{7} \right)^2 \left( 65 + 34\sqrt{3} + 26\sqrt{7} + 15\sqrt{21}\right)^3\\ j\left( \tfrac{\sqrt{30}i}{1} \right) &= \tfrac{1}{16} \left(10 + 7\sqrt{2} + 4\sqrt{5} + 3\sqrt{10} \right)^4 \left( 55 + 30\sqrt{2} + 12\sqrt{5} + 10\sqrt{10} \right)^3\\ j\left( \tfrac{\sqrt{30}i}{2} \right) &= \tfrac{1}{16} \left(10 + 7\sqrt{2} - 4\sqrt{5} - 3\sqrt{10} \right)^4 \left( 55 + 30\sqrt{2} - 12\sqrt{5} - 10\sqrt{10} \right)^3\\ j\left( \tfrac{\sqrt{30}i}{5} \right) &= \tfrac{1}{16} \left(10 - 7\sqrt{2} + 4\sqrt{5} - 3\sqrt{10} \right)^4 \left( 55 - 30\sqrt{2} + 12\sqrt{5} - 10\sqrt{10} \right)^3\\ j\left( \tfrac{\sqrt{30}i}{10} \right) &= \tfrac{1}{16} \left(10 - 7\sqrt{2} - 4\sqrt{5} + 3\sqrt{10} \right)^4 \left( 55 - 30\sqrt{2} - 12\sqrt{5} + 10\sqrt{10} \right)^3\\ j\left(\tfrac{1+\sqrt{31}i}{2}\right)&=\left(1-\left(1+\tfrac{\sqrt{19}}{2}\left(\left(\tfrac{13-\sqrt{93}}{13+\sqrt{93}}\right)^{1/2}\left(\tfrac{\sqrt{31}+\sqrt{27}}{\sqrt{31}-\sqrt{27}}\right)^{1/3}+\left(\tfrac{13+\sqrt{93}}{13-\sqrt{93}}\right)^{1/2}\left(\tfrac{\sqrt{31}-\sqrt{27}}{\sqrt{31}+\sqrt{27}}\right)^{1/3}\right)\right)^2\right)^3\\ j(\sqrt{70}i) &= \left(1 + \tfrac{9}{4}\left(303 + 220\sqrt{2} + 139\sqrt{5} + 96\sqrt{10}\right)^2 \right)^3\\ j(\sqrt{94}i) &= \left(1 + \tfrac{9}{64}\left(2454 + 1736\sqrt{2} + \left(546 + 384\sqrt{2}\right)\sqrt{9 + 8\sqrt{2}} + \left(527 + 373\sqrt{2} + \left(117 + 83\sqrt{2}\right)\sqrt{9 + 8\sqrt{2}}\right)\sqrt{3 + 4\sqrt{2} + 3\sqrt{9 + 8\sqrt{2}}}\right)^2\right)^3\\ j(7i) &= \left( 1 + \tfrac{9}{4}\sqrt{21+8\sqrt{7}} \left(30 + 11\sqrt{7} + \left (6+\sqrt{7} \right )\sqrt{21+8\sqrt{7}}\right)^2 \right)^3\\ j(8i) &= \left( 1 + \tfrac{9}{4} \sqrt[4]{2} \left (1 + \sqrt{2} \right) \left(123 + 104\sqrt[4]{2} + 88\sqrt{2} + 73\sqrt[4]{8}\right)^2 \right)^3\\ j(10i) &= \left(1 + \tfrac{9}{8}\left(2402 + 1607\sqrt[4]{5} + 1074\sqrt[4]{25} + 719\sqrt[4]{125}\right)^2 \right)^3\\ j \left( \frac{5 \, i}{2} \right) &= \left(1 + \tfrac{9}{8}\left(2402 - 1607\sqrt[4]{5} + 1074\sqrt[4]{25} - 719\sqrt[4]{125}\right)^2 \right)^3\\ j(2\sqrt{58}i) &= \left(1+\tfrac{9}{256}\left(1+\sqrt{2}\right)^5\left(5+\sqrt{29}\right)^5\left(793+907\sqrt{2}+237\sqrt{29}+103\sqrt{58}\right)^2\right)^3\\ j\left( \tfrac{1 + \sqrt{1435}i}{2} \right) &= \left( 1 - 9 \left ( 9892538 + 4424079\sqrt{5} + 1544955\sqrt{41} + 690925\sqrt{205} \right )^2 \right)^3\\ j\left( \tfrac{1 + \sqrt{1555}i}{2} \right) &= \left( 1 - 9 \left ( 22297077 + 9971556\sqrt{5} + \left ( 3571365 + 1597163\sqrt{5} \right ) \sqrt{\tfrac{31 + 21\sqrt{5}}{2}} \right)^2 \right)^3\\ \end{align} \)

Several special values were calculated in 2014:[9]

\( \begin{align} j \left( \frac{5 \, i + 1}{2} \right) &= \left( \frac{2927 - 1323 \, \sqrt{5}}{2} \right)^3,\\ j \left( 5 \, i \right) &= \left( \frac{2927 + 1323 \, \sqrt{5}}{2} \right)^3,\\ j \left( \frac{5 \, i + 2}{4} \right) &= \Bigg( \frac{\left( 1 + \sqrt{5} \, \right)^{37}}{2^{39}} \Bigg( 1190448488 - 858585699 \, \sqrt{2} + 540309076 \, \sqrt{5} - 374537880 \, \sqrt{10} \, - \, \sqrt[4]{5} \left( 693172512 - 595746414 \, \sqrt{2} + 407357424 \, \sqrt{5} - 240819696 \, \sqrt{10} \, \right) \Bigg) \Bigg)^3,\\ j \left( \frac{10 \, i + 1}{2} \right) &= \Bigg( \frac{\left( 1 + \sqrt{5} \, \right)^{37}}{2^{39}} \Bigg( 1190448488 - 858585699 \, \sqrt{2} + 540309076 \, \sqrt{5} - 374537880 \, \sqrt{10} \, + \, \sqrt[4]{5} \left( 693172512 - 595746414 \, \sqrt{2} + 407357424 \, \sqrt{5} - 240819696 \, \sqrt{10} \, \right) \Bigg) \Bigg)^3,\\ j \left( \frac{5 \, i}{4} \right) &= \Bigg( \frac{\left( 1 + \sqrt{5} \, \right)^{37}}{2^{39}} \Bigg( 1190448488 + 858585699 \, \sqrt{2} + 540309076 \, \sqrt{5} + 374537880 \, \sqrt{10} \, - \, \sqrt[4]{5} \left( 693172512 + 595746414 \, \sqrt{2} + 407357424 \, \sqrt{5} + 240819696 \, \sqrt{10} \, \right) \Bigg) \Bigg)^3,\\ j(20 \, i) &= \Bigg( \frac{\left( 1 + \sqrt{5} \, \right)^{37}}{2^{39}} \Bigg( 1190448488 + 858585699 \, \sqrt{2} + 540309076 \, \sqrt{5} + 374537880 \, \sqrt{10} \, + \, \sqrt[4]{5} \left( 693172512 + 595746414 \, \sqrt{2} + 407357424 \, \sqrt{5} + 240819696 \, \sqrt{10} \, \right) \Bigg) \Bigg)^3. \end{align} \)

All preceding values are real. A complex conjugate pair might be inferred exploiting the symmetry described in the reference, along with the values for j(10 i) and j(5 i/2), given above:

\( \begin{align} j \left( \frac{5 \, i \pm 1}{4} \right) &= \left(1 - \tfrac{9}{8}\left((2402 - 1074\sqrt{5}) \, i \pm (1607 - 719\sqrt{5}) \sqrt[4]{5} \right)^2 \right)^3. \end{align} \)

Four more special values are given as two complex conjugate pairs:[10]

\( \begin{align} j \left( \frac{4 \left( 5 \, i \pm 1 \right)}{13} \right) = \Bigg( \frac{\left( 1 - \sqrt{5} \, \right)^{37}}{2^{39}} \Bigg( 1190448488 - 858585699 \, \sqrt{2} - 540309076 \, \sqrt{5} + 374537880 \, \sqrt{10} \, \pm \, \textit{i} \, \sqrt[4]{5} \left( 693172512 - 595746414 \, \sqrt{2} - 407357424 \, \sqrt{5} + 240819696 \, \sqrt{10} \, \right) \Bigg) \Bigg)^3,\\ j \left( \frac{5 \left( 4 \, i \pm 1 \right)}{17} \right) = \Bigg( \frac{\left( 1 - \sqrt{5} \, \right)^{37}}{2^{39}} \Bigg( 1190448488 + 858585699 \, \sqrt{2} - 540309076 \, \sqrt{5} - 374537880 \, \sqrt{10} \, \pm \, \textit{i} \, \sqrt[4]{5} \left( 693172512 + 595746414 \, \sqrt{2} - 407357424 \, \sqrt{5} - 240819696 \, \sqrt{10} \, \right) \Bigg) \Bigg)^3 \end{align} \)

References

Silverman, Joseph H. (1986). The Arithmetic of Elliptic Curves. Graduate Texts in Mathematics 106. Springer-Verlag. p. 339. ISBN 0-387-96203-4. Zbl 0585.14026.
Petersson, Hans (1932). Über die Entwicklungskoeffizienten der automorphen Formen. Acta Mathematica 58 (1). pp. 169–215. doi:10.1007/BF02547776. MR 1555346.
Rademacher, Hans (1938). The Fourier coefficients of the modular invariant j(τ). American Journal of Mathematics 60 (2) (The Johns Hopkins University Press). pp. 501–512. doi:10.2307/2371313. JSTOR 2371313. MR 1507331.
Cummins, C.J. (2004). "Congruence subgroups of groups commensurable with PSL(2,Z)$ of genus 0 and 1". Exp. Math. 13 (3): 361–382. ISSN 1058-6458. Zbl 1099.11022.
Chandrasekharan (1985) p.108
Chandrasekharan, K. (1985), Elliptic Functions, Grundlehren der mathematischen Wissenschaften 281, Springer-Verlag, p. 110, ISBN 3-540-15295-4, Zbl 0575.33001
Girondo, Ernesto; González-Diez, Gabino (2012), Introduction to compact Riemann surfaces and dessins d'enfants, London Mathematical Society Student Texts 79, Cambridge: Cambridge University Press, p. 267, ISBN 978-0-521-74022-7, Zbl 1253.30001
Chudnovsky, David V.; Chudnovsky, Gregory V. (1989), "The Computation of Classical Constants", Proceedings of the National Academy of Sciences of the United States of America 86 (21): 8178–8182, doi:10.1073/pnas.86.21.8178, ISSN 0027-8424, JSTOR 34831, PMC 298242, PMID 16594075.
Adlaj, Semjon. "Multiplication and division on elliptic curves, torsion points and roots of modular equations" (PDF). Retrieved 17 October 2014.

Adlaj, Semjon (2014). "Torsion points on elliptic curves and modular polynomial symmetries" (PDF). The joined MSU-CCRAS Computer Algebra Seminar. Moscow, Russia.

Apostol, Tom M. (1976), Modular functions and Dirichlet Series in Number Theory, Graduate Texts in Mathematics 41, New York: Springer-Verlag, MR 0422157. Provides a very readable introduction and various interesting identities.
Apostol, Tom M. (1990), Modular functions and Dirichlet Series in Number Theory (2nd ed.), ISBN 0-387-97127-0, MR 1027834
Berndt, Bruce C.; Chan, Heng Huat (1999), "Ramanujan and the modular j-invariant" (PDF), Canadian Mathematical Bulletin 42 (4): 427–440, doi:10.4153/CMB-1999-050-1, MR 1727340. Provides a variety of interesting algebraic identities, including the inverse as a hypergeometric series.
Cox, David A. (1989), Primes of the Form x^2 + ny^2: Fermat, Class Field Theory, and Complex Multiplication, New York: Wiley-Interscience Publication, John Wiley & Sons Inc., MR 1028322 Introduces the j-invariant and discusses the related class field theory.
Conway, John Horton; Norton, Simon (1979), "Monstrous moonshine", Bulletin of the London Mathematical Society 11 (3): 308–339, doi:10.1112/blms/11.3.308, MR 0554399. Includes a list of the 175 genus-zero modular functions.
Rankin, Robert A. (1977), Modular forms and functions, Cambridge: Cambridge University Press, ISBN 0-521-21212-X, MR 0498390. Provides a short review in the context of modular forms.
Schneider, Theodor (1937), "Arithmetische Untersuchungen elliptischer Integrale", Math. Annalen 113: 1–13, doi:10.1007/BF01571618, MR 1513075.

Mathematics Encyclopedia