Elementary crystals

Note

Each of these crystals will work with any Cartan matrix input (with weights from the weight lattice corresponding to the Cartan matrix given).

T-crystal

Let $\lambda$ be a weight. As defined in [Kashiwara1993] (see, also, [Kashiwara1995]) the crystal $T_{\lambda }=\{t_{\lambda }\}$ is a single element crystal with the crystal structure defined by

$$\mathrm{wt}(t_{\lambda })=\lambda ,e_i{t}_{\lambda }=f_i{t}_{\lambda }=0,{\epsilon }_i(t_{\lambda })={\phi }_i(t_{\lambda })=-\mathrm{\infty }.$$

The crystal $T_{\lambda }$ shifts the weights of the vertices in a crystal $B$ by $\lambda$ when tensored with $B$, but leaves the graph structure of $B$ unchanged. That is, for all $b\in B$, we have $\mathrm{wt}(t_{\lambda }\otimes b)=\mathrm{wt}(b)+\lambda$:

Sage

sage: B = crystals.Tableaux(['A',2],shape=[2,1])
sage: T = crystals.elementary.T(['A',2], B.Lambda()[1] + B.Lambda()[2])
sage: V = crystals.TensorProduct(T,B)
sage: for x in V:
....:     print(x.weight())
....:
(4, 2, 0)
(3, 3, 0)
(3, 2, 1)
(3, 1, 2)
(2, 2, 2)
(4, 1, 1)
(3, 2, 1)
(2, 3, 1)
sage: for x in B:
....:     print(x.weight() + T[0].weight())
....:
(4, 2, 0)
(3, 3, 0)
(3, 2, 1)
(3, 1, 2)
(2, 2, 2)
(4, 1, 1)
(3, 2, 1)
(2, 3, 1)

Python

>>> from sage.all import *
>>> B = crystals.Tableaux(['A',Integer(2)],shape=[Integer(2),Integer(1)])
>>> T = crystals.elementary.T(['A',Integer(2)], B.Lambda()[Integer(1)] + B.Lambda()[Integer(2)])
>>> V = crystals.TensorProduct(T,B)
>>> for x in V:
...     print(x.weight())
....:
(4, 2, 0)
(3, 3, 0)
(3, 2, 1)
(3, 1, 2)
(2, 2, 2)
(4, 1, 1)
(3, 2, 1)
(2, 3, 1)
>>> for x in B:
...     print(x.weight() + T[Integer(0)].weight())
....:
(4, 2, 0)
(3, 3, 0)
(3, 2, 1)
(3, 1, 2)
(2, 2, 2)
(4, 1, 1)
(3, 2, 1)
(2, 3, 1)

Sage Live

B = crystals.Tableaux(['A',2],shape=[2,1])
T = crystals.elementary.T(['A',2], B.Lambda()[1] + B.Lambda()[2])
V = crystals.TensorProduct(T,B)
for x in V:
    print(x.weight())
for x in B:
    print(x.weight() + T[0].weight())

Warning

Sage uses the opposite convention for the tensor product rule to Kashiwara’s definition, so care must be taken when comparing the examples here with Kashiwara’s papers.

Here is an example using a hyperbolic Cartan matrix:

Sage

sage: A = CartanMatrix([[2,-4],[-4,2]])
sage: La = RootSystem(A).weight_lattice().fundamental_weights()
sage: La
Finite family {0: Lambda[0], 1: Lambda[1]}
sage: T = crystals.elementary.T(A,La[1])
sage: T
The T crystal of type [ 2 -4]
[-4  2] and weight Lambda[1]

Python

>>> from sage.all import *
>>> A = CartanMatrix([[Integer(2),-Integer(4)],[-Integer(4),Integer(2)]])
>>> La = RootSystem(A).weight_lattice().fundamental_weights()
>>> La
Finite family {0: Lambda[0], 1: Lambda[1]}
>>> T = crystals.elementary.T(A,La[Integer(1)])
>>> T
The T crystal of type [ 2 -4]
[-4  2] and weight Lambda[1]

Sage Live

A = CartanMatrix([[2,-4],[-4,2]])
La = RootSystem(A).weight_lattice().fundamental_weights()
La
T = crystals.elementary.T(A,La[1])
T

C-crystal

Defined in [Kashiwara1993], the component crystal $C=\{c\}$ is the single element crystal whose crystal structure is defined by

$$\mathrm{wt}(c)=0,e_{i}c=f_{i}c=0,{\epsilon }_i(c)={\phi }_i(c)=0.$$

Note $C\cong B(0)$, where $B(0)$ is the highest weight crystal of highest weight $0$.

The crystal $C\otimes T_{\mu }$ is useful when finding subcrystals inside irreducible highest weight crystals $B(\lambda )$ where $\lambda$ is larger than $\mu$ in the lexicographic order. For example:

Sage

sage: P = RootSystem("C2").weight_lattice()
sage: La = P.fundamental_weights()
sage: h = P.simple_coroots()
sage: T = crystals.elementary.T("C2", 2*La[1])
sage: C = crystals.elementary.Component(P)
sage: B = crystals.TensorProduct(C,T)
sage: b = B(C[0],T[0])
sage: for i in B.index_set(): print(b.epsilon(i))
-2
0
sage: for i in B.index_set(): print(b.phi(i))
0
0
sage: for i in B.index_set(): print(b.f(i))
None
None
sage: for i in B.index_set(): print(b.e(i))
None
None

Python

>>> from sage.all import *
>>> P = RootSystem("C2").weight_lattice()
>>> La = P.fundamental_weights()
>>> h = P.simple_coroots()
>>> T = crystals.elementary.T("C2", Integer(2)*La[Integer(1)])
>>> C = crystals.elementary.Component(P)
>>> B = crystals.TensorProduct(C,T)
>>> b = B(C[Integer(0)],T[Integer(0)])
>>> for i in B.index_set(): print(b.epsilon(i))
-2
0
>>> for i in B.index_set(): print(b.phi(i))
0
0
>>> for i in B.index_set(): print(b.f(i))
None
None
>>> for i in B.index_set(): print(b.e(i))
None
None

Sage Live

P = RootSystem("C2").weight_lattice()
La = P.fundamental_weights()
h = P.simple_coroots()
T = crystals.elementary.T("C2", 2*La[1])
C = crystals.elementary.Component(P)
B = crystals.TensorProduct(C,T)
b = B(C[0],T[0])
for i in B.index_set(): print(b.epsilon(i))
for i in B.index_set(): print(b.phi(i))
for i in B.index_set(): print(b.f(i))
for i in B.index_set(): print(b.e(i))

This new crystal can be summarized into the R-crystal below.

R-crystal

For a fixed weight $\lambda$, the crystal $R_{\lambda }=\{r_{\lambda }\}$ is a single element crystal with the crystal structure defined by

$$\mathrm{wt}(r_{\lambda })=\lambda ,e_i{r}_{\lambda }=f_i{r}_{\lambda }=0,{\epsilon }_i(r_{\lambda })=-⟨h_i,\lambda ⟩,{\phi }_i(r_{\lambda })=0,$$

where $\{h_i\}$ are the simple coroots. See page 146 [Joseph1995], for example, for more details. (Note that in [Joseph1995], this crystal is denoted by $S_{\lambda }$.)

Tensoring $R_{\lambda }$ with a crystal $B$ results in shifting the weights of the vertices in $B$ by $\lambda$ and may also cut a subset out of the original graph of $B$.

Warning

Sage uses the opposite convention for the tensor product rule to Kashiwara’s definition, so care must be taken when comparing the examples here with some of the literature.

For example, suppose $\mu \le \lambda$ in lexicographic ordering on weights, and one wants to see $B(\mu )$ as a subcrystal of $B(\lambda )$. Then $B(\mu )$ may be realized as the connected component of $R_{\mu -\lambda }\otimes B(\lambda )$ containing the highest weight $r_{\mu -\lambda }\otimes u_{\lambda }$, where $u_{\lambda }$ is the highest weight vector in $B(\lambda )$:

Sage

sage: La = RootSystem(['B',4]).weight_lattice().fundamental_weights()
sage: Bla = crystals.NakajimaMonomials(['B',4], La[1]+La[2])
sage: Bmu = crystals.NakajimaMonomials(['B',4], La[1])
sage: R = crystals.elementary.R(['B',4], -La[2])
sage: T = crystals.TensorProduct(R,Bla)
sage: mg = mg = T(R[0], Bla.module_generators[0])
sage: S = T.subcrystal(generators=[mg])
sage: G = T.digraph(subset=S)
sage: Bmu.digraph().is_isomorphic(G, edge_labels=True)
True
sage: view(G, tightpage=True) # optional - dot2tex graphviz, not tested (opens external window)

Python

>>> from sage.all import *
>>> La = RootSystem(['B',Integer(4)]).weight_lattice().fundamental_weights()
>>> Bla = crystals.NakajimaMonomials(['B',Integer(4)], La[Integer(1)]+La[Integer(2)])
>>> Bmu = crystals.NakajimaMonomials(['B',Integer(4)], La[Integer(1)])
>>> R = crystals.elementary.R(['B',Integer(4)], -La[Integer(2)])
>>> T = crystals.TensorProduct(R,Bla)
>>> mg = mg = T(R[Integer(0)], Bla.module_generators[Integer(0)])
>>> S = T.subcrystal(generators=[mg])
>>> G = T.digraph(subset=S)
>>> Bmu.digraph().is_isomorphic(G, edge_labels=True)
True
>>> view(G, tightpage=True) # optional - dot2tex graphviz, not tested (opens external window)

Sage Live

La = RootSystem(['B',4]).weight_lattice().fundamental_weights()
Bla = crystals.NakajimaMonomials(['B',4], La[1]+La[2])
Bmu = crystals.NakajimaMonomials(['B',4], La[1])
R = crystals.elementary.R(['B',4], -La[2])
T = crystals.TensorProduct(R,Bla)
mg = mg = T(R[0], Bla.module_generators[0])
S = T.subcrystal(generators=[mg])
G = T.digraph(subset=S)
Bmu.digraph().is_isomorphic(G, edge_labels=True)
view(G, tightpage=True) # optional - dot2tex graphviz, not tested (opens external window)

$i$-th elementary crystal

For $i$ an element of the index set of type $X$, the crystal $B_i$ of type $X$ is the set

$$B_i=\{b_i(m):m\in \mathbf{Z}\},$$

where the crystal structure is given by $\mathrm{wt}(b_i(m))=m{\alpha }_i$ and

$$\begin{array}{r}\begin{array}{rlrl}{\phi }_j(b_i(m))& =\{\begin{array}{ll}m& \text{ if }j=i,\\ -\mathrm{\infty }& \text{ if }j\ne i,\end{array}& {\epsilon }_j(b_i(m))& =\{\begin{array}{ll}-m& \text{ if }j=i,\\ -\mathrm{\infty }& \text{ if }j\ne i,\end{array}\\ f_j{b}_i(m)& =\{\begin{array}{ll}{b}_i(m-1)& \text{ if }j=i,\\ 0& \text{ if }j\ne i,\end{array}& e_j{b}_i(m)& =\{\begin{array}{ll}{b}_i(m+1)& \text{ if }j=i,\\ 0& \text{ if }j\ne i.\end{array}\end{array}\end{array}$$

See [Kashiwara1993] or [Kashiwara1995] for more information. Here is an example:

Sage

sage: B = crystals.elementary.Elementary("A2",1)
sage: S = B.subcrystal(max_depth=4, generators=[B(0)])
sage: sorted(s for s in S)
[-4, -3, -2, -1, 0, 1, 2, 3, 4]
sage: G = B.digraph(subset=S)
sage: view(G, tightpage=True) # optional - dot2tex graphviz, not tested (opens external window)

Python

>>> from sage.all import *
>>> B = crystals.elementary.Elementary("A2",Integer(1))
>>> S = B.subcrystal(max_depth=Integer(4), generators=[B(Integer(0))])
>>> sorted(s for s in S)
[-4, -3, -2, -1, 0, 1, 2, 3, 4]
>>> G = B.digraph(subset=S)
>>> view(G, tightpage=True) # optional - dot2tex graphviz, not tested (opens external window)

Sage Live

B = crystals.elementary.Elementary("A2",1)
S = B.subcrystal(max_depth=4, generators=[B(0)])
sorted(s for s in S)
G = B.digraph(subset=S)
view(G, tightpage=True) # optional - dot2tex graphviz, not tested (opens external window)

Warning

To reiterate, Sage uses the opposite convention for the tensor product rule to Kashiwara’s definition. In particular, using Sage’s convention, one has $T_{\lambda }\otimes B_i\cong B_i\otimes T_{s_i\lambda }$, where $s_i$ is the $i$-th simple reflection.

Make live