OFFSET
4,1
COMMENTS
It is strongly believed that this sequence is never zero. If this were to hold true, these conformal blocks divisors would form a basis of Pic(M_{0,n})^{S_n}.
REFERENCES
A. Beauville, Conformal blocks, fusion rules and the Verlinde formula, Bar-Ilan University, Ramat Gan, Israël, 1996, pages 75-99.
K. Ueno, Conformal field theory with gauge symmetry, American Mathematical Society, 2008
LINKS
A. Beauville, Conformal blocks, fusion rules and the Verlinde formula, arXiv:alg-geom/9405001, 1994.
N. Fakhruddin, Chern classes of conformal blocks on M_{0,n}, arxiv:0904.2918 [math.AG], 2009-2011.
EXAMPLE
For n=6 the intersection matrix is [[0,3],[2,0]] giving a determinant of -6. In general these matrices have a high degree of symmetry and periodicity.
MATHEMATICA
slncbIntersection[n_, i_, j_] := Module[{k = n-Mod[i j, n]}, If[Floor[i j/ n] + Floor[j (n-i-2)/n] == j-2, If[1 <= k <= j, Return[k], If[j <= k <= 2j-1, Return[2j-k]], Return[0]]]; 0];
slnIntersectionMatrix[n_] := Module[{matdict}, Do[r = slncbIntersection[n, i, j]; matdict[i-1, j-2] = r, {i, 1, n/2-1}, {j, 2, n/2}]; Return[ Table[ matdict[i, j], {i, 0, n/2-2}, {j, 0, n/2-2}]]];
Table[Det[slnIntersectionMatrix[n]], {n, 4, 31}] (* Jean-François Alcover, Aug 10 2018, translated from Sage *)
PROG
(Macaulay 2 with package ConfBlocks) slnl1intmat = (n) -> ( if even(n) then g=lift(n/2-1, ZZ) else g=lift((n-1)/2-1, ZZ); cu:=sn1curves(n); M = apply(#cu, i -> apply(g, j -> (1/1)*CdotCBslml1(cu_i, n, 1, apply(n, k -> j+2)))); matrix M )
(Sage)
# Based on a formula in "Conformal Blocks Divisors on M_{0, n}"
# with Maxim Arap, Angela Gibney and David Swinarski.
def sln_cb_intersection(n, i, j):
"""Gives the intersection on M_{0, n} of the sl_n level 1, omega_j^n conformal blocks divisor with the F-curve F_{1, 1, i, n-i - 2}"""
k = n - (i*j % n)
if (i*j/n).floor() + (j*(n-i-2)/n).floor() == j-2:
if k in [1..j]:
return k
elif k in [j..2*j-1]:
return 2*j - k
else:
return 0
def sln_intersection_matrix(n):
"""gives the g X g intersection matrix of sln level 1 omega_j^n conformal blocks divisors with 1, 1, i F-curves"""
matdict = dict()
for j in [2..(n/2).floor()]:
for i in [1..(n/2).floor()-1]:
r = sln_cb_intersection(n, i, j)
if r != 0:
matdict[(i-1, j-2)] = r
return matrix(matdict)
[sln_intersection_matrix(n).determinant() for n in [4..104]]
CROSSREFS
KEYWORD
sign
AUTHOR
Jim Stankewicz (stankewicz(AT)gmail.com), Mar 19 2010
STATUS
approved