[go: up one dir, main page]

login
Search: a159842 -id:a159842
     Sort: relevance | references | number | modified | created      Format: long | short | data
Number of inequivalent sublattices of index n in hexagonal lattice, where two sublattices are equivalent if they are related by a rotation or reflection preserving the hexagonal lattice.
(Formerly M0420)
+10
24
1, 1, 2, 3, 2, 3, 3, 5, 4, 4, 3, 8, 4, 5, 6, 9, 4, 8, 5, 10, 8, 7, 5, 15, 7, 8, 9, 13, 6, 14, 7, 15, 10, 10, 10, 20, 8, 11, 12, 20, 8, 18, 9, 17, 16, 13, 9, 28, 12, 17, 14, 20, 10, 22, 14, 25, 16, 16, 11, 34, 12, 17, 21, 27, 16, 26, 13, 24, 18, 26, 13, 40, 14
OFFSET
1,3
COMMENTS
The hexagonal lattice is the familiar 2-dimensional lattice in which each point has 6 neighbors. This is sometimes called the triangular lattice.
From Andrey Zabolotskiy, Mar 10 2018: (Start)
If only primitive sublattices are considered, we get A003050.
Here only rotations and reflections preserving the parent hexagonal lattice are allowed. If reflections are not allowed, we get A145394. If any rotations and reflections are allowed, we get A300651.
In other words, the parent lattice of the sublattices under consideration has Patterson symmetry group p6mm, and two sublattices are considered equivalent if they are related via a symmetry from that group [Rutherford]. For other 2D Patterson groups, the analogous sequences are A000203 (p2), A069734 (p2mm), A145391 (c2mm), A145392 (p4), A145393 (p4mm), A145394 (p6).
Rutherford says at p. 161 that his sequence for p6mm differs from this sequence, but it seems that with the current definition and terms of this sequence, this actually is his p6mm sequence, and the sequence he thought to be this one is actually A300651. Also, he says that a(n) != A300651(n) only when A002324(n) > 2 (first time happens at n = 49), but actually these two sequences differ at other terms, too, for example, at n = 42 (see illustration). (End)
REFERENCES
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
LINKS
A. Altshuler, Construction and enumeration of regular maps on the torus, Discrete Math. 4 (1973), 201-217.
A. Altshuler, Construction and enumeration of regular maps on the torus, Discrete Math. 4 (1973), 201-217. [Annotated and corrected scanned copy]
M. Bernstein, N. J. A. Sloane and P. E. Wright, On Sublattices of the Hexagonal Lattice, Discrete Math. 170 (1997) 29-39 (Abstract, pdf, ps).
Amihay Hanany, Domenico Orlando, and Susanne Reffert, Sublattice counting and orbifolds, High Energ. Phys., 2010 (2010), 51, arXiv.org:1002.2981 [hep-th] (see Table 3).
Daejun Kim, Seok Hyeong Lee, and Seungjai Lee, Zeta functions enumerating subforms of quadratic forms, arXiv:2409.05625 [math.NT], 2024.
W. Kurth, Enumeration of Platonic maps on the torus, Discrete Math. 61 (1986), 71-83.
Andrey Zabolotskiy, Sublattices of the hexagonal lattice (illustrations for n = 1..7, 14)
Andrey Zabolotskiy, Coweight lattice A^*_n and lattice simplices, arXiv:2003.10251 [math.CO], 2020.
FORMULA
a(n) = Sum_{ m^2 | n } A003050(n/m^2).
a(n) = (A000203(n) + 2*A002324(n) + 3*A145390(n))/6. [Rutherford] - N. J. A. Sloane, Mar 13 2009
a(n) = Sum_{ d|n } A112689(d+1). - Andrey Zabolotskiy, Aug 29 2019
a(n) = Sum_{ d|n } floor(d/6) + 1 - 1*[d == 2 or 6 (mod 12)] + 1*[d == 4 (mod 12)]. [Kurth] - Brahadeesh Sankarnarayanan, Feb 24 2023
MATHEMATICA
max = 73; A145390 = Drop[ CoefficientList[ Series[ Sum[(1 + Cos[n*Pi/2])*x^n/(1 - x^n), {n, 1, max}], {x, 0, max}], x], 1]; A002324[n_] := (dn = Divisors[n]; Count[dn, _?(Mod[#, 3] == 1 & )] - Count[dn, _?(Mod[#, 3] == 2 & )]); a[n_] := (DivisorSigma[1, n] + 2 A002324[n] + 3*A145390[[n]])/6; Table[a[n], {n, 1, max}] (* Jean-François Alcover, Oct 11 2011, after given formula *)
KEYWORD
nonn,nice,easy
STATUS
approved
Number of inequivalent sublattices of index n in a square lattice, where two sublattices are considered equivalent if one can be rotated to give the other.
+10
8
1, 1, 2, 2, 4, 3, 6, 4, 8, 7, 8, 6, 14, 7, 12, 10, 16, 9, 20, 10, 18, 16, 18, 12, 30, 13, 20, 20, 28, 15, 30, 16, 32, 24, 26, 20, 46, 19, 30, 26, 38, 21, 48, 22, 42, 33, 36, 24, 62, 29, 38, 34, 46, 27, 60, 30, 60, 40, 44, 30, 70, 31, 48, 52, 64, 33, 72, 34, 60, 48, 60
OFFSET
0,3
COMMENTS
If reflections are allowed, we get A054346. If only rotations that preserve the parent square lattice are allowed, we get A145392. The analog for a hexagonal lattice is A054384.
LINKS
Daejun Kim, Seok Hyeong Lee, and Seungjai Lee, Zeta functions enumerating subforms of quadratic forms, arXiv:2409.05625 [math.NT], 2024. See section 6.2 for the Dirichlet g.f. zeta^SL_{x^2+y^2}(s).
Andrey Zabolotskiy, Sublattices of the square lattice (illustrations for n = 1..6, 15, 25).
EXAMPLE
For n = 1, 2, 3, 4 the sublattices are generated by the rows of:
[1 0] [2 0] [2 0] [3 0] [3 0] [4 0] [4 0] [2 0] [2 0]
[0 1] [0 1] [1 1] [0 1] [1 1] [0 1] [1 1] [0 2] [1 2].
PROG
(SageMath)
# see A159842 for the definitions of dc, fin, u, N
def ff(m, k1, minus = True):
def f(n):
if n == 1: return 1
r = 1
for (p, k) in factor(n):
if p % 4 != m or k1 and k > 1: return 0
if minus: r *= (-1)**k
return r
return f
f1, f2, f3 = ff(1, True), ff(1, True, False), ff(3, False)
def a_SL(n):
return (dc(u, N, f1)(n) + dc(u, f3)(n)) / 2
print([a_SL(n) for n in range(1, 100)]) # Andrey Zabolotskiy, Sep 22 2024
CROSSREFS
KEYWORD
nonn,easy,nice
AUTHOR
N. J. A. Sloane, May 06 2000
STATUS
approved
Number of inequivalent sublattices of index n in square lattice, where two sublattices are considered equivalent if one can be rotated or reflected to give the other.
+10
7
1, 1, 2, 2, 4, 3, 5, 3, 7, 5, 7, 4, 11, 5, 8, 8, 12, 6, 13, 6, 15, 10, 11, 7, 21, 9, 13, 12, 18, 9, 21, 9, 21, 14, 16, 13, 29, 11, 17, 16, 28, 12, 28, 12, 25, 21, 20, 13, 39, 16, 24, 20, 29, 15, 34, 18, 36, 22, 25, 16, 47, 17, 26, 29, 38, 21, 40, 18, 36, 26, 36, 19, 58, 20
OFFSET
0,3
COMMENTS
If we count sublattices as equivalent only if they are related by a rotation, we get A054345 instead of this sequence. If we only allow rotations and reflections that preserve the parent (square) lattice, we get A145393; the first discrepancy is at n = 25 (see illustration), the second is at n = 30. If both restrictions are applied, i.e., only rotations preserving the parent lattice are allowed, we get A145392. The analog for the hexagonal lattice is A300651. - Andrey Zabolotskiy, Mar 12 2018
LINKS
Daejun Kim, Seok Hyeong Lee, and Seungjai Lee, Zeta functions enumerating subforms of quadratic forms, arXiv:2409.05625 [math.NT], 2024. See section 6.2 for the Dirichlet g.f. zeta^GL_{x^2+y^2}(s).
Andrey Zabolotskiy, Sublattices of the square lattice (illustrations for n = 1..6, 15, 25)
EXAMPLE
For n = 1, 2, 3, 4 the sublattices are generated by the rows of:
[1 0] [2 0] [2 0] [3 0] [3 0] [4 0] [4 0] [2 0] [2 0]
[0 1] [0 1] [1 1] [0 1] [1 1] [0 1] [1 1] [0 2] [1 2].
PROG
(SageMath)
# See A159842 and A054345 for the definitions of functions used here
def a_GL(n):
return (a_SL(n) + dc(fin(1, 0, 0, 1), u, u, f2)(n)) / 2
print([a_GL(n) for n in range(1, 100)]) # Andrey Zabolotskiy, Sep 22 2024
CROSSREFS
KEYWORD
nonn,easy,nice
AUTHOR
N. J. A. Sloane, May 06 2000
STATUS
approved
Number of inequivalent sublattices of index n in hexagonal lattice, where two sublattices are considered equivalent if one can be rotated to give the other.
+10
5
1, 1, 1, 2, 3, 2, 4, 3, 5, 5, 6, 4, 10, 5, 7, 8, 11, 6, 13, 7, 14, 10, 12, 8, 20, 11, 13, 14, 17, 10, 24, 11, 21, 16, 18, 14, 31, 13, 19, 18, 30, 14, 28, 15, 28, 26, 24, 16, 42, 17, 31, 24, 31, 18, 40, 24, 35, 26, 30, 20, 56, 21, 31, 31, 43, 26, 48, 23, 42, 32, 42, 24, 65
OFFSET
0,4
COMMENTS
The hexagonal lattice is the familiar 2-dimensional lattice in which each point has 6 neighbors. This is sometimes called the triangular lattice.
If reflections are allowed, we get A300651. If only rotations that preserve the parent hexagonal lattice are allowed, we get A145394. The analog for square lattice is A054345. - Andrey Zabolotskiy, Mar 10 2018
LINKS
M. Bernstein, N. J. A. Sloane and P. E. Wright, On Sublattices of the Hexagonal Lattice, Discrete Math. 170 (1997) 29-39 (Abstract, pdf, ps).
Daejun Kim, Seok Hyeong Lee, and Seungjai Lee, Zeta functions enumerating subforms of quadratic forms, arXiv:2409.05625 [math.NT], 2024. See section 6.1 for the Dirichlet g.f. zeta^SL_{x^2+xy+y^2}(s).
Andrey Zabolotskiy, Sublattices of the hexagonal lattice (illustrations for n = 1..7, 14)
PROG
(SageMath)
# see A159842 for the definitions of dc, fin, u, N
def gg(m, k1, minus = True):
def f(n):
if n == 1: return 1
r = 1
for (p, k) in factor(n):
if p % 3 != m or k1 and k > 1: return 0
if minus: r *= (-1)**k
return r
return f
g1, g2, g3 = gg(1, True), gg(1, True, False), gg(2, False)
def a_SL(n):
return (dc(u, N, g1)(n) + 2 * dc(u, g3)(n)) / 3
print([a_SL(n) for n in range(1, 100)]) # Andrey Zabolotskiy, Sep 22 2024
CROSSREFS
KEYWORD
nonn,nice,easy
AUTHOR
N. J. A. Sloane, May 08 2000
STATUS
approved
Number of inequivalent sublattices of index n in hexagonal lattice, where two sublattices are considered equivalent if they are related by any rotation or reflection.
+10
5
1, 1, 2, 3, 2, 3, 3, 5, 4, 4, 3, 8, 4, 5, 6, 9, 4, 8, 5, 10, 8, 7, 5, 15, 7, 8, 9, 13, 6, 14, 7, 15, 10, 10, 10, 20, 8, 11, 12, 20, 8, 17, 9, 17, 16, 13, 9, 28, 11, 17, 14, 20, 10, 22, 14, 25, 16, 16, 11, 34, 12, 17, 20, 27, 16, 26, 13, 24, 18, 24, 13, 40
OFFSET
1,3
COMMENTS
If we count sublattices as equivalent only if they are related by a rotation, we get A054384 instead of this sequence. If we only allow rotations and reflections that preserve the parent (hexagonal) lattice, we get A003051; the first discrepancy is at n = 42 (see illustration), the second is at n = 49. If both restrictions are applied, i.e., only rotations preserving the parent lattice are allowed, we get A145394. The analog for square lattice is A054346.
Although A003051 has its counterpart A003050 which counts primitive sublattices only, this sequence has no such counterpart sequence because a primitive sublattice can turn to a non-primitive one via a non-parent-lattice-preserving rotation, so the straightforward definition of primitiveness does not work in this case.
LINKS
Daejun Kim, Seok Hyeong Lee, and Seungjai Lee, Zeta functions enumerating subforms of quadratic forms, arXiv:2409.05625 [math.NT], 2024. See section 6.1 for the Dirichlet g.f. zeta^GL_{x^2+xy+y^2}(s).
Andrey Zabolotskiy, Sublattices of the hexagonal lattice (illustrations for n = 1..7, 14)
PROG
(SageMath)
# See A159842 and A054384 for the definitions of functions used here
def a_GL(n):
return (a_SL(n) + dc(fin(1, -1, 0, 2), u, u, g2)(n)) / 2
print([a_GL(n) for n in range(1, 100)]) # Andrey Zabolotskiy, Sep 22 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Andrey Zabolotskiy, Mar 10 2018
STATUS
approved
Three-dimensional simplices of determinant n.
+10
4
1, 2, 3, 7, 5, 10, 7, 20, 14, 18, 11, 41, 15, 28, 31, 58, 21, 60, 25, 77, 49, 54, 33, 144, 50, 72, 75, 123, 49, 158, 55, 177, 97, 112, 99, 268, 75, 136, 129, 286, 89, 268, 97, 249, 218, 190, 113, 496, 146, 280
OFFSET
0,2
COMMENTS
Two simplices are considered equal if some integer affine automorphism sends the first to the second.
Probably the same as A159842. - Andrey Zabolotskiy, Mar 12 2018
CROSSREFS
Cf. A159842.
KEYWORD
nice,nonn
AUTHOR
Jacques-Olivier Moussafir (msfr(AT)ceremade.dauphine.fr)
STATUS
approved
Number of symmetrically distinct sublattices (supercells, superlattices, HNFs) of the simple cubic lattice of index n.
+10
4
1, 3, 3, 9, 5, 13, 7, 24, 14, 23, 11, 49, 15, 33, 31, 66, 21, 70, 25, 89, 49, 61, 33, 162, 50, 81, 75, 137, 49, 177, 55, 193, 97, 123, 99, 296, 75, 147, 129, 312, 89, 291, 97, 269, 218, 203, 113, 534, 146, 302, 203, 357, 141, 451, 207, 508, 247, 307, 171, 789
OFFSET
1,2
PROG
(Python)
# see A159842 for the definition of dc, fin, per, u, N, N2
def a(n): # from DeCross's slides
return (dc(u, N, N2)(n) + 6*dc(fin(1, -1, 0, 4), u, u, N)(n)
+ 3*dc(fin(1, 3), u, u, N)(n)
+ 8*dc(fin(1, 0, -1, 0, 0, 0, 0, 0, 3), u, u, per(0, 1, -1))(n)
+ 6*dc(fin(1, 1), u, u, per(0, 1, 0, -1))(n))//24
print([a(n) for n in range(1, 300)])
# Andrey Zabolotskiy, Sep 02 2019
KEYWORD
nonn
AUTHOR
Andrey Zabolotskiy, Mar 12 2018
EXTENSIONS
Terms a(11) and beyond from Andrey Zabolotskiy, Sep 02 2019
STATUS
approved
Number of symmetrically distinct sublattices (supercells, superlattices, HNFs) of the 3D hexagonal lattice of index n.
+10
4
1, 3, 5, 11, 7, 19, 11, 34, 23, 33, 19, 77, 25, 53, 55, 104, 37, 115, 45, 143, 91, 105, 61, 272, 90, 139, 137, 235, 91, 309, 103, 331, 183, 219, 185, 516, 141, 267, 245, 544, 169, 529, 185, 485, 411, 375, 217, 952, 278, 550, 389, 647, 271, 829, 397, 922, 477
OFFSET
1,2
LINKS
Gus L. W. Hart and Rodney W. Forcade, Algorithm for generating derivative structures, Phys. Rev. B 77, 224115 (2008), DOI: 10.1103/PhysRevB.77.224115 [see Table IV].
Materials Simulation Group, Derivative structure enumeration library
Kohei Shinohara, Atsuto Seko, Takashi Horiyama, Masakazu Ishihata, Junya Honda and Isao Tanaka, Enumeration of nonequivalent substitutional structures using advanced data structure of binary decision diagram, J. Chem. Phys. 153, 104109 (2020); preprint: Derivative structure enumeration using binary decision diagram, arXiv:2002.12603 [physics.comp-ph], 2020.
PROG
(Python)
# see A159842 for the definitions of dc, fin, per, u, N, N2
def a(n):
return (dc(u, N, N2)(n) + 6*dc(fin(1, -1, 0, 4), u, u, N)(n)
+ dc(fin(1, 3), u, u, N)(n)
+ 4*dc(fin(1, 0, 1), u, u, per(0, 1, -1))(n)) // 12
print([a(n) for n in range(1, 100)])
# Andrey Zabolotskiy, Feb 03 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Andrey Zabolotskiy, Mar 12 2018
EXTENSIONS
Terms a(11) and beyond from Andrey Zabolotskiy, Feb 03 2020
STATUS
approved
Number of four-dimensional simplical toric diagrams with hypervolume n.
+10
3
1, 2, 4, 10, 8, 19, 13, 45, 33, 47, 30, 129, 43, 96, 108, 226, 78, 264, 102, 357, 226, 277, 163, 813, 260, 425, 436, 780, 297, 1092, 355, 1281, 678, 856, 712, 2215, 569, 1155, 1050, 2537, 752, 2544, 856, 2447, 2048, 1944, 1093, 5388, 1447, 3083, 2150, 3827
OFFSET
1,2
COMMENTS
Also gives the number of distinct abelian orbifolds of C^5/Gamma, Gamma in SU(5).
LINKS
J. Davey, A. Hanany and R. K. Seong, Counting Orbifolds, J. High Energ. Phys. (2010) 2010: 10; arXiv:1002.3609 [hep-th], 2010.
A. Hanany and R. K. Seong, Symmetries of abelian orbifolds, J. High Energ. Phys. (2011) 2011: 27; arXiv:1009.3017 [hep-th], 2010-2011. Table 5 gives a(1)-a(80), but the terms a(36) and a(65) there are apparently erroneous.
Andrey Zabolotskiy, Coweight lattice A^*_n and lattice simplices, arXiv:2003.10251 [math.CO], 2020.
PROG
(Sage)
# see Python in A159842 for the definition of dc, fin, per, u, N, N2
def fin_d(d):
return fin(*(d.get(n+1, 0) for n in range(max(d))))
def a(n): # see Hanany & Seong 2011, Table 1 row D=5 and Table 9
return (dc(u, N, N2, lambda n: n**3)(n) +
10 * dc(u, u, N, N2, fin(1, -1, 0, 8))(n) +
15 * dc(u, u, N, N, fin_d({1: 1, 2: -3, 4: 14, 8: -12, 16: 16}))(n) +
20 * dc(u, u, N, per(0, 1, -1), fin(1, 0, -1, 0, 0, 0, 0, 0, 9))(n) +
20 * dc(u, u, u, per(0, 1, -1), fin(1, -1, 0, 2), fin(1, 0, -1, 0, 0, 0, 0, 0, 3))(n) +
30 * dc(u, u, u, per(0, 1, 0, -1), fin_d({1: 1, 2: -2, 4: 3, 16: 6, 32: -8, 64: 8}))(n) +
24 * dc(u, per(0, 1, -1, -1, 1), per(0, 1, I, -I, -1), per(0, 1, -I, I, -1))(n)) / 120
print([a(n) for n in range(1, 100)])
CROSSREFS
Cf. A003051 (No. of two-dimensional triangular toric diagrams of area n), A045790 (No. of three-dimensional tetrahedral toric diagrams of volume n), A173877, A173878.
KEYWORD
nonn
AUTHOR
Rak-Kyeong Seong (rak-kyeong.seong(AT)imperial.ac.uk), Feb 25 2010
EXTENSIONS
a(16) corrected, terms a(31) and beyond added from Hanany & Seong 2011 by Andrey Zabolotskiy, Jun 30 2019
a(36) corrected from 2202 to 2215 by Andrey Zabolotskiy, Sep 20 2022
STATUS
approved
Number of symmetrically distinct sublattices (supercells, superlattices, HNFs) of the tetragonal lattice of index n.
+10
3
1, 5, 5, 17, 9, 29, 13, 51, 28, 53, 25, 115, 33, 81, 73, 153, 51, 176, 61, 219, 121, 161, 85, 403, 126, 213, 188, 353, 129, 473, 145, 487, 257, 335, 261, 776, 201, 405, 345, 815, 243, 801, 265, 731, 584, 569, 313, 1407, 398, 838, 559, 975, 393, 1256, 573, 1375
OFFSET
1,2
PROG
(Python)
# see A159842 for the definition of dc, fin, per, u, N, N2
def a(n):
return (dc(u, N, N2)(n) + 2*dc(fin(1, -1, 0, 4), u, u, N)(n)
+ 3*dc(fin(1, 3), u, u, N)(n)
+ 2*dc(fin(1, 1), u, u, per(0, 1, 0, -1))(n)) // 8
print([a(n) for n in range(1, 300)])
# Andrey Zabolotskiy, Jan 31 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Andrey Zabolotskiy, Mar 12 2018
EXTENSIONS
Terms a(11) and beyond from Andrey Zabolotskiy, Jan 31 2020
STATUS
approved

Search completed in 0.012 seconds