login
Theta series of square lattice with respect to deep hole.
(Formerly M3319)
7

%I M3319 #33 Mar 12 2021 22:24:41

%S 4,8,4,8,8,0,12,8,0,8,8,8,4,8,0,8,16,0,8,0,4,16,8,0,8,8,0,8,8,8,4,16,

%T 0,0,8,0,16,8,8,8,0,0,12,8,0,8,16,0,8,8,0,16,0,0,0,16,12,8,8,0,8,8,0,

%U 0,8,8,16,8,0,8,8,0,12,8,0,0,16,0,8,8,0,24,0,8,8,0,0,8,8,0,4,16,8,8,16,0,0

%N Theta series of square lattice with respect to deep hole.

%C Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).

%C In [Jacobi 1829] on page 105 is equation 18: "2 k K / Pi = 4 sqrt(q) + 8 sqrt(q^5) + 4 sqrt(q^9) [...]". - _Michael Somos_, Sep 09 2012

%D J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag, p. 106.

%D N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

%D C. G. J. Jacobi, Fundamenta Nova Theoriae Functionum Ellipticarum, 1829.

%H G. C. Greubel, <a href="/A005883/b005883.txt">Table of n, a(n) for n = 0..1000</a>

%H Michael Somos, <a href="/A010815/a010815.txt">Introduction to Ramanujan theta functions</a>

%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/RamanujanThetaFunctions.html">Ramanujan Theta Functions</a>

%F Expansion of Jacobi theta constant q^(-1/2) * theta_2(q)^2 in powers of q^2. - _Michael Somos_, Oct 31 2006

%F G.f.: 4 * (Product_{k>0} (1 - x^k) * (1 + x^(2*k))^2)^2. - _Michael Somos_, Oct 31 2006

%F From _Michael Somos_, Sep 09 2012: (Start)

%F Expansion of 4 * psi(x)^2 in powers of x where psi() is a Ramanujan theta function.

%F Expansion of q^(-1) * (1/2) * (1 - k') * K / (Pi/2) in powers of q^4 where k', K are Jacobi elliptic functions.

%F Expansion of q^(-1/2) * k * K / (Pi/2) in powers of q^2 where k, K are Jacobi elliptic functions.

%F Expansion of q^(-1/4) * 2 * k^(1/2) * K / (Pi/2) in powers of q where k, K are Jacobi elliptic functions.

%F Expansion of 4 * q^(-1/4) * eta(q^2)^4 / eta(q)^2 in powers of q.

%F a(n) = 4 * A008441(n). (End)

%e 4 + 8*x + 4*x^2 + 8*x^3 + 8*x^4 + 12*x^6 + 8*x^7 + 8*x^9 + 8*x^10 + 8*x^11 + ...

%e 4*q + 8*q^3 + 4*q^5 + 8*q^7 + 8*q^9 + 12*q^13 + 8*q^15 + 8*q^19 + 8*q^21 + ...

%e Theta = 4*q^(1/2) + 8*q^(5/2) + 4*q^(9/2) + 8*q^(13/2) + 8*q^(17/2) + ...

%t a[0] = 4; a[n_] := 4*DivisorSum[4n+1, (-1)^Quotient[#, 2]&]; Table[a[n], {n, 0, 100}] (* _Jean-François Alcover_, Nov 04 2015, translated from PARI *)

%t s = 4*(QPochhammer[q^2]^4/QPochhammer[q]^2)+O[q]^100; CoefficientList[s, q] (* _Jean-François Alcover_, Nov 09 2015 *)

%o (PARI) {a(n) = if( n<0, 0, n = 4*n + 1; 4 * sumdiv(n, d, (-1)^(d\2)))} /* _Michael Somos_, Oct 31 2006 */

%Y Cf. A008441.

%K nonn

%O 0,1

%A _N. J. A. Sloane_