%I M4509 #52 Mar 12 2021 22:24:41

%S 8,32,48,64,104,96,112,192,144,160,256,192,248,320,240,256,384,384,

%T 304,448,336,352,624,384,456,576,432,576,640,480,496,832,672,544,768,

%U 576,592,992,768,640,968,672,864,960,720,896,1024,960,784,1248,816,832,1536

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

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

%C The D_4 lattice is the set of all integer quadruples [a, b, c, d] where a + b + c + d is even. The deep holes are quadruples [a, b, c, d] where each coordinate is half an odd integer and where a + b + c + d is even. - _Michael Somos_, May 23 2012

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

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

%H T. D. Noe, <a href="/A005879/b005879.txt">Table of n, a(n) for n = 0..10000</a>

%H G. Nebe and N. J. A. Sloane, <a href="http://www.math.rwth-aachen.de/~Gabriele.Nebe/LATTICES/D4.html">Home page for this lattice</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>

%H <a href="/index/Da#D4">Index entries for sequences related to D_4 lattice</a>

%F Expansion of Jacobi theta_2(q)^4/(2q) in powers of q^2. - _Michael Somos_, Apr 11 2004

%F Expansion of q^(-1/2) * 8 * (eta(q^2)^2 / eta(q))^4 in powers of q. - _Michael Somos_, Apr 11 2004

%F Expansion of 8 * psi(x)^4 in powers of x where psi() is a Ramanujan theta function. - _Michael Somos_, May 23 2012

%F Expansion of (phi(q)^4 - phi(-q)^4) / (2 * q) in powers of q^2. - _Michael Somos_, May 23 2012

%F G.f.: 8 * (Product_{k>0} (1 - x^k) * (1 + x^k)^2)^4. - _Michael Somos_, Apr 11 2004

%F a(n) = 8 * A008438(n) = 4 * A005880(n) = A000118(2*n + 1) = - A096727(2*n + 1). - _Michael Somos_, Nov 01 2006

%e 8 + 32*x + 48*x^2 + 64*x^3 + 104*x^4 + 96*x^5 + 112*x^6 + 192*x^7 + ...

%e 8*q + 32*q^3 + 48*q^5 + 64*q^7 + 104*q^9 + 96*q^11 + 112*q^13 + ...

%e .

%e For n = 2 the objects counted are the ways to represent the integer 5 = (2*n+1) as a sum of 4 squares, 0 and negative numbers allowed.

%e [-2,-1,0,0], [-2,0,-1,0], [-2,0,0,-1], [-2,0,0,1], [-2,0,1,0], [-2,1,0,0],

%e [-1,-2,0,0], [-1,0,-2,0], [-1,0,0,-2], [-1,0,0,2], [-1,0,2,0], [-1,2,0,0],

%e [0,-2,-1,0], [0,-2,0,-1], [0,-2,0,1], [0,-2,1,0], [0,-1,-2,0], [0,-1,0,-2],

%e [0,-1,0,2], [0,-1,2,0], [0,0,-2,-1], [0,0,-2,1], [0,0,-1,-2], [0,0,-1,2],

%e [0,0,1,-2], [0,0,1,2], [0,0,2,-1], [0,0,2,1], [0,1,-2,0], [0,1,0,-2],

%e [0,1,0,2], [0,1,2,0], [0,2,-1,0], [0,2,0,-1], [0,2,0,1], [0,2,1,0],

%e [1,-2,0,0], [1,0,-2,0], [1,0,0,-2], [1,0,0,2], [1,0,2,0], [1,2,0,0],

%e [2,-1,0,0], [2,0,-1,0], [2,0,0,-1], [2,0,0,1], [2,0,1,0], [2,1,0,0].

%e - _Peter Luschny_, Nov 03 2015

%p S:= series(JacobiTheta2(0,q)^4/(2*q), q, 202):

%p seq(coeff(S,q,2*j),j=0..100); # _Robert Israel_, Nov 03 2015

%t (* a(n) gives the number of ways to represent the integer 2n+1 as a sum of 4 squares *) a[n_] := SquaresR[4, 2n+1]; Table[a[n], {n, 0, 52}] (* _Jean-François Alcover_, Nov 03 2015 *)

%t terms = 53; QP = QPochhammer; s = 8 QP[q^2]^8/QP[q]^4 + O[q]^terms; CoefficientList[s, q] (* _Jean-François Alcover_, Jul 07 2017, after _Michael Somos_ *)

%o (PARI) {a(n) = if( n<0, 0, 8 * sigma(2*n + 1))} /* _Michael Somos_, Apr 11 2004 */

%o (PARI) q='q+O('q^66); Vec(8*(eta(q^2)^2/eta(q))^4) \\ _Joerg Arndt_, Nov 03 2015

%Y Cf. A000118, A005880, A008438, A096727.

%K nonn,easy,nice

%O 0,1

%A _N. J. A. Sloane_