R-Matrix Invariants

2009-05-22T11:25:18Z

RacroCdelt:

http://www.textcaelracrelle.com
{{Manual TOC Sidebar}}

KnotTheory` can compute knot and link invariants associated with matrix solutions of the Yang-Baxter equation, using the program <code>REngine</code>:

{{Startup Note}}


{{HelpAndAbout|
n = 1 |
n1 = 2 |
in = <nowiki>REngine</nowiki> |
out= <nowiki>REngine[K, Rp, Rn, Mcupl, Mcupr, Mcapl Mcapr] returns the invariant associated with the given R-matrices (Rp for positive crossings, Rn for negative crossings) and oriented creation and annihilation M matrices, of the oriented knot or link K. See the Manual for details of convention. Note that REngine does not verify that the given matrices actually define an invariant, use TestRMatrix[..] for this purpose.</nowiki> |
about= <nowiki>REngine was written by Siddarth Sankaran at the University of Toronto, in the summer of 2005</nowiki>}}


In more detail, let <math>V</math> be a free module of finite rank over a commutative ring, and let <math> R, \bar{R} \in End(V \otimes V) </math> be invertible. Further, let there be invertible endomorphisms <math>M_\leftarrow, M_\rightarrow, M^\leftarrow, M^\rightarrow \in End(V) </math>; these correspond to McupL, McupR, McapL, and McapR respectively. Fixing a basis of V, we can regard all of these endomorphisms as matrices.

The following notation is useful: for <math>n = dim(V)</math>, <math>R_{ab}^{cd}</math> refers to the element of '''''R''''' in row <math>n(c-1)+d</math> and column <math>n(a-1) + b</math>.

For an oriented knot or link <math>K</math>, REngine returns the product <math> \tau(K) </math>, which is computed as follows:
#Find a [[MorseLink Presentations| Morse link presentation]] of <math>K</math> such that it is composed only of crossings of the <code>X[k, Over/Under, Down, Down]</code> variety; this is accomplished by adding caps and cups where necessary.
#Label the edges in between events in the Morse link presentation of <math>K</math>.
#Associate each event with one of the matrices as shown, with the indices given by the labelling:
#*<code>X[k, Over, Down, Down]</code>:[[Image:REposcr.png|Positive crossing]]<math>\Rightarrow R_{ab}^{cd}</math> (positive crossing)
#*<code>X[k, Under,Down, Down]</code>:[[Image:REnegcr.png|Negative crossing]] <math>\Rightarrow\bar{R}_{ab}^{cd}</math> (negative crossing)
#*<code>Cup[k, k+1]</code>: [[Image:REmcupr.png|Cup right]] <math>\Rightarrow {M_\rightarrow}_{ab}</math>
#*<code>Cup[k, k-1]</code>: [[Image:REmcupl.png|Cup left]] <math>\Rightarrow {M_\leftarrow}_{ab}</math>
#*<code>Cap[k, k+1]</code>: [[Image:REmcapr.png|Cap right]] <math>\Rightarrow {M^\rightarrow}_{ab}</math>
#*<code>Cap[k, k-1]</code>: [[Image:REmcapl.png|Cap left]] <math>\Rightarrow {M^\leftarrow}_{ab}</math>
#Define <math>\tau(K)</math> as the result of taking the product of the matrices associated with the elements of ''K'', and summing over repeated indices; each sum runs from 1 to <math>n=dim(V)</math>.

[[Image:REtref.png|left|thumb|150px|A trefoil knot]]
As an example, let ''K'' be the trefoil as shown, with the chosen axis upwards, and the strands going down through all the crossings. Then <math>\tau(K) = \sum {M_\leftarrow}_{ab} {M_\rightarrow}_{cd} R_{bc}^{ef} R_{ef}^{gh} R_{gh}^{ij} {M^\rightarrow}_{ai}{M^\leftarrow}_{jd}</math>

where the summation is carried out on all the indices ''a'' through ''i'', each ranging from 1 to n.

Determining that <math>\tau(K)</math> is actually an invariant is simply a matter of checking the equality of the outcomes of Reidemeister moves in various configurations and orientations, as well as a few 'topological' moves. Note that <code>REngine</code> does not perform any checks, while the utility [[TestRMatrix]] performs the tests required to determine regular isotopy.

Let's use the REngine to find the Jones polynomial of the trefoil. The '''''R''''' and '''''M''''' matrices for the Jones polynomial, in the indeterminate t, are as follows:

:<math>R = \begin{pmatrix} t^{1/2} & 0&0&0 \\ 0&0&t&0 \\ 0&t&t^{1/2} - t^{3/2}&0\\0&0&0&t^{1/2} \end{pmatrix} , \bar{R} = R^{-1} </math>

:<math>M^\leftarrow = M^\rightarrow = M_\leftarrow = M_\rightarrow = \begin{pmatrix} 0&t^{-1/4} \\ t^{1/4}&0 \end{pmatrix} </math>

Putting these definitions into Mathematica, we run the REngine on the [[3_1|trefoil]]. Note that the Jones polynomial is usually normalized so that the polynomial of the unknot [[0_1]] is 1, hence the term on the denominator.


{{In|
n = 3 |
in = <nowiki>r={{Sqrt[t], 0, 0, 0},{0, 0, t, 0},{0,t, Sqrt[t] - t^(3/2), 0}, {0,0,0,Sqrt[t]}};
rb = Inverse[r]; mcupl=mcapl=mcupr=mcapr={{0, t^(-1/4)},{t^(1/4), 0}};</nowiki>}}




{{InOut|
n = 4 |
in = <nowiki>jones = REngine[Knot[3,1], r, rb, mcupl, mcapl, mcupr, mcapr] / REngine[Knot[0,1], r, rb, mcupl, mcapr, mcupr, mcapr] //Apart</nowiki> |
out= <nowiki> -4 -3 1
-t + t + -
t</nowiki>}}


===See also===
See also [[QuantumGroups`]], a Mathematica package by Scott Morrison, now included in KnotTheory`, which computes arbitrary quantum knot invariants. (For now, it only uses braid representations, not morse link representations.)

Knot Atlas - User contributions [en]

R-Matrix Invariants