"Root systems and Dynkin diagrams(mathematica)"의 두 판 사이의 차이
(피타고라스님이 이 페이지의 이름을 root systems and Dynkin diagrams(mathematica)로 바꾸었습니다.) |
|
(차이 없음)
| |
2010년 3월 14일 (일) 13:22 판
- Root Systems and Dynkin diagrams
- http://en.wikipedia.org/wiki/root_systems
- [1][2][3][4]http://en.wikipedia.org/wiki/Dynkin_diagram
A_n root systems
(* A_n type Cartan matrix *)
r := 2
rt[i_] := UnitVector[r + 1, i] - UnitVector[r + 1, i + 1]
a[i_, j_] := (2 Dot[rt[i], rt[j]])/Dot[rt[j], rt[j]]
A := Table[a[i, j], {i, 1, r}, {j, 1, r}]
A // MatrixForm
B_n root systems
Clear[rt]
(*B_r type Cartan matrix*)
r := 4
rt[i_] :=
If[i < r, UnitVector[r, i] - UnitVector[r, 1 + i], UnitVector[r, r]]
a[i_, j_] := (2 Dot[rt[i], rt[j]])/Dot[rt[j], rt[j]]
A := Table[a[i, j], {i, 1, r}, {j, 1, r}]
Print["root vectors"]
Table[rt[i], {i, 1, r}] // TableForm
Print["Cartan matrix"]
A // MatrixForm
C_n root systems
Clear[rt]
(*C_r type Cartan matrix*)
r := 4
rt[i_] :=
If[i < r, UnitVector[r, i] - UnitVector[r, 1 + i], 2*UnitVector[r, r]]
a[i_, j_] := (2 Dot[rt[i], rt[j]])/Dot[rt[j], rt[j]]
A := Table[a[i, j], {i, 1, r}, {j, 1, r}]
Print["root vectors"]
Table[rt[i], {i, 1, r}] // TableForm
Print["Cartan matrix"]
A // MatrixForm
D_n root systems
Clear[rt]
(*D_r type Cartan matrix*)
r := 6
rt[i_] :=
If[i < r, UnitVector[r, i] - UnitVector[r, 1 + i],
UnitVector[r, r - 1] + UnitVector[r, r]]
a[i_, j_] := (2 Dot[rt[i], rt[j]])/Dot[rt[j], rt[j]]
A := Table[a[i, j], {i, 1, r}, {j, 1, r}]
Print["root vectors"]
Table[rt[i], {i, 1, r}] // TableForm
Print["Cartan matrix"]
A // MatrixForm
G2
E_6 root system
Clear[rt]
(*E_ 6 type Cartan matrix*)
r := 6
rt[1] := {0, 0, 0, 0, 0, 0, 0, 1, \[Minus]1}
rt[2] := {0, 0, 0, 0, 0, 0, 1, -1, 0}
rt[3] := {1, -2, 1, -2, 1, 1, -2, 1, 1}/3
rt[4] := {0, 0, 0, 1, -1, 0, 0, 0, 0}
rt[5] := {0, 0, 0, 0, 1, -1, 0, 0, 0}
rt[6] := {0, 1, -1, 0, 0, 0, 0, 0, 0}
a[i_, j_] := (2 Dot[rt[i], rt[j]])/Dot[rt[j], rt[j]]
A := Table[a[i, j], {i, 1, r}, {j, 1, r}]
Print["root vectors"]
Table[rt[i], {i, 1, r}] // TableForm
Print["Cartan matrix"]
A // MatrixForm
http://en.wikipedia.org/wiki/E6_(mathematics)
E_7 root system
Clear[rt]
(*E_ 7 type Cartan matrix*)
r := 7
alp := Sum[UnitVector[r + 1, i]/2, {i, 1, 4}] -
Sum[UnitVector[r + 1, i]/2, {i, 5, 8}]
rt[i_] :=
Piecewise[{{UnitVector[r + 1, i + 2] - UnitVector[r + 1, 1 + i],
i < 7}, {alp, i == 7}}, {i, 1, r}]
a[i_, j_] := (2 Dot[rt[i], rt[j]])/Dot[rt[j], rt[j]]
A := Table[a[i, j], {i, 1, r}, {j, 1, r}]
Print["root vectors"]
Table[rt[i], {i, 1, r}] // TableForm
Print["Cartan matrix"]
A // MatrixForm
E_8 root system
Clear[rt]
(*E_ 8 type Cartan matrix*)
alp := -Sum[UnitVector[r, i]/2, {i, 1, r}]
r := 8
rt[i_] :=
Piecewise[{{UnitVector[r, i] - UnitVector[r, 1 + i],
i < 7}, {UnitVector[r, i] + UnitVector[r, i - 1], i == 7}, {alp,
i == 8}}, {i, Range[r]}]
a[i_, j_] := (2 Dot[rt[i], rt[j]])/Dot[rt[j], rt[j]]
A := Table[a[i, j], {i, 1, r}, {j, 1, r}]
Print["root vectors"]
Table[rt[i], {i, 1, r}] // TableForm
Print["Cartan matrix"]
A // MatrixForm