View source for Aharonov–Jones–Landau algorithm

In [[computer science]], the '''Aharonov–Jones–Landau algorithm''' is an efficient [[quantum algorithm]] for obtaining an additive [[approximation algorithm|approximation]] of the [[Jones polynomial]] of a given link at an arbitrary [[root of unity]]. Finding a multiplicative approximation is a [[Sharp-P|#P]]-hard problem,<ref>{{cite arXiv |last=Kuperberg |first=Greg |date=2009 |title=How hard is it to approximate the Jones polynomial? |eprint=0908.0512 }}</ref> so a better approximation is considered unlikely. However, it is known that computing an additive approximation of the Jones polynomial is a [[BQP]]-complete problem.<ref>{{cite arXiv |first1=Michael | last1=Freedman |first2=Michael | last2=Larsen |first3=Zhenghan | last3=Wang |date=2000 |title=A modular functor which is universal for quantum computation |eprint=quant-ph/0001108 }}</ref>

The algorithm was published in 2009 in a paper written by [[Dorit Aharonov]], [[Vaughan Jones]] and [[Zeph Landau]]. 

== History ==
In the early 2000s, a series of papers by [[Michael Freedman]], [[Alexei Kitaev]], [[Michael J. Larsen]], and [[Zhenghan Wang]] demonstrated that [[Topological quantum computer|topological quantum computers]] based on [[topological quantum field theory]] had the same computational power as [[Quantum circuit|quantum circuits]]. In particular, they showed that the braiding of anyons in the [[Fibonacci category]] could be used to approximate the Jones polynomial evaluated at a [[Root of unity|primitive 5th root of unity]]. They then showed that this problem was [[BQP]]-complete.

Putting these results together, this implies that there is a polynomial length quantum circuit which approximates the Jones polynomial at 5th roots of unity. This algorithm was completely inaccessible to ordinary quantum computer scientists, however, since the papers by Freedman-Kitaev-Larsen-Wang used heavy machinery from [[Manifold|manifold topology]]. The contribution of Aharanov-Jones-Landau was to simplify this complicated implicit algorithm in such a way that it would be palatable to a larger audience.

== The Markov trace ==

The first idea behind the algorithm is to find a more tractable description for the operation of evaluating the Jones polynomial. This is done by means of the Markov trace.

The "Markov trace" is a trace operator defined on the [[Temperley–Lieb algebra]] <math>TL_n(d)</math> as follows: given a <math>T\in TL_n(d)</math> which is a single [[Kauffman diagram]], let <math>tr(T)=d^{a-n}</math> where <math>a</math> is the number of loops attained by identifying each point in the bottom of <math>T</math>'s Kauffman diagram with the corresponding point on top. This extends linearly to all of <math>TL_n(d)</math>.

The Markov trace is a trace operator in the sense that <math>\operatorname{tr}(1)=1</math> and <math>\operatorname{tr}(XY)=\operatorname{tr}(YX)</math> for any <math>X,Y\in TL_n(d)</math>. It also has the additional property that if <math>X</math> is a Kauffman diagram whose rightmost strand goes straight up then <math>\operatorname{tr}(XE_{n-1}) = \frac{1}{d}\operatorname{tr}(X)</math>.

A useful fact exploited by the AJL algorithm is that the Markov trace is the unique trace operator on <math>TL_n(d)</math> with that property.<ref>{{cite journal |first=V.F.R |last=Jones |title=Index for subfactors |journal=Invent Math |volume=1 |issue=72  |year=1983 |doi=10.1007/BF01389127|bibcode=1983InMat..72....1J }}</ref>

== Representing <math>B_n</math> in <math>TL_n(d)</math> ==

For a complex number <math>A</math> we define the map <math>\rho_A:B_n\to TL_n(d)</math> via <math>\sigma_i\mapsto AE_i+A^{-1}I</math>. It follows by direct calculation that if <math>A</math> satisfies that <math>d= -A^2-A^{-2}</math> then <math>\rho_A</math> is a representation.

Given a braind <math>B\in B_n</math> let <math>B^{tr}</math> be the link attained by identifying the bottom of the diagram with its top like in the definition of a Markov trace, and let <math>V_{B^{tr}}</math> be the result link's Jones polynomial. The following relation holds:

:<math>V_{B^{tr}}(A^{-4})=(-A)^{3w(B^{tr})}d^{n-1}\operatorname{tr}(\rho_A(B))</math>

where <math>w</math> is the [[writhe]]. As the writhe can be easily calculated classically, this reduces the problem of approximating the Jones polynomial to that of approximating the Markov trace.

== The path model representation of <math>TL_n(d)</math> ==

We wish to construct a complex representation <math>\tau</math> of <math>TL_n(d)</math> such that the representation <math>\tau\circ\rho_A</math> of <math>B_n</math> will be unitary. We also wish that our representation will have a straightforward encoding into qubits.

Let 
:<math>Q_{n,k}=\left\{ q\in\left\{ 1,\ldots,k-1\right\} ^{n+1}\mid q(1)=1,\left|q(i)-q(i+1)\right|=1\right\} </math>
and let <math>V_{n,k}=\mathbb{C}[Q_{n,k}]</math> be the vector space which has <math>Q_{n,k}</math> as an orthonormal basis.

We choose define a linear map <math>\tau:TL_n(d)\to V_{n,k}</math> by defining it on the base of generators <math>\{1,E_1,\ldots,E_{n-1}\}</math>. To do so we need to define the matrix element <math>\tau(E_i)_{q,q'}</math> for any <math>q,q'\in Q_{n,k}</math>.

We say that <math>q</math> and <math>q'</math> are 'compatible' if <math>q(j)=q'(j)</math> for any <math>j\ne i+1</math> and <math>q(i)=q(i+2)</math>. Geometrically this means that if we put <math>q</math> and <math>q'</math> below and above the Kauffman diagram in the gaps between the strands then no connectivity component will touch two gaps which are labeled by different numbers.

If <math>q</math> and <math>q'</math> are incompatible set <math>\tau(E_i)_{q,q'}=0</math>. Else, let <math>l</math> be either <math>q(i)</math> or <math>q(i+2)</math> (at least one of these number must be defined, and if both are defined they must be equal) and set

:<math>\tau\left(E_{i}\right)_{q,q'}=\begin{cases}
\frac{\lambda_{l+1}}{\lambda_{l}} & q\left(i+1\right)=q'(i+1)>l\\
\frac{\sqrt{\lambda_{l-1}\lambda_{l+1}}}{\lambda_{l}} & q\left(i+1\right)\ne q'(i+1)\\
\frac{\lambda_{l-1}}{\lambda_{l}} & q\left(i+1\right)=q'(i+1)<l
\end{cases}</math>

where <math>\lambda_l=\sin(\pi l/k)</math>. Finally set <math>d=2\cos(\pi/k)</math>.

This representation, known as the ''path model representation'', induces a unitary representation of the braid group.<ref>{{cite journal |first=V.F.R |last=Jones |title=A polynomial invariant for knots via von Neumann algebras |journal=Bull. Amer. Math. Soc. |volume=12  |year=1985 |pages=103–111 |doi=10.1090/s0273-0979-1985-15304-2|doi-access=free }}</ref><ref>{{cite journal |first=V.F.R |last=Jones |title=Braid groups, Hecke Algebras and type II factors|journal=Geometric methods in Operator Algebras |volume=123  |year=1986 |pages=242-273 }}</ref> Moreover, it holds that <math>d=-A^2-A^{-2}</math> for <math>A=ie^{-i\pi/2k}</math>.

This implies that if we could approximate the Markov trace in this representation this will allow us to approximate the Jones polynomial in <math>A^{-4}=e^{2\pi i/k}</math>.

== A quantum version of the path model representation ==

In order to be able to act on elements of the path model representation by means of quantum circuits, we need to encode the elements of <math>Q_{n,k}</math> into qubits in a way which allows us to easily describe the images of the generators <math>E_i</math>.

We represent each path as a sequence of moves, where <math>0</math> indicates a move to the right and <math>1</math> indicates a move to the left. For example, the path <math>1,2,3,2</math> will be represented by the state <math>\left|001\right\rangle </math>.

This encodes <math>\mathbb{C}[Q_{n,k}]</math> as a subspace of the state space on <math>k-1</math> qubits.

We define the operators <math>\varphi_i</math> within this subspace we define

:<math>\Phi_{i}\left|w\right\rangle =\begin{cases}
0 & w_{i}=w_{i+1}\\
\frac{\lambda_{z_{i}-\left(-1\right)^{w_{i}}}}{\lambda_{z_{i}}}\left|w\right\rangle +\frac{\sqrt{\lambda_{z_{i}-1}\lambda_{z_{i}+1}}}{\lambda_{z_{i}}}X_{i}X_{i+1}\left|w\right\rangle  & w_{i}\ne w_{i+1}
\end{cases}</math>

where <math>X_i</math> is the [[Pauli matrices|Pauli matrix]] flipping the <math>i</math>th bit and <math>z_i</math> is the position of the path represented by <math>\left|w\right\rangle</math> after <math>i</math> steps.

We arbitrarily extend <math>\varphi_i</math> to be the identity on the rest of the space.

We note that mapping <math>E_i\mapsto\varphi_i</math> retains all the properties of the path model representation. Specifically, it induces a unitary representation <math>\varphi</math> of <math>B_n</math>. It is fairly straightforward to show that <math>\varphi_i</math> can be implemented by <math>\text{poly}\left(n,k\right)</math> gates, so it follows that <math>\varphi(B)</math> can be implemented for any <math>B\in B_n</math> using <math>\text{poly}\left(m,n,k\right)</math> where <math>m</math> is the number of crossings in <math>B</math>.

== A quantum version of the Markov trace ==

The benefit of this construction is that it gives us a way to represent the Markov trace in a way which can be easily approximated.

Let <math>\mathcal{H}_{n,k}</math> be the subspace of paths we described in the previous clause, and let <math>\mathcal{H}_{n,k,l}</math> be the subspace spanned by basis elements which represent walks which end on the <math>l</math>-th position.

Note that each of the operators <math>\varphi_i</math> fix <math>\mathcal{H}_{n,k,l}</math> setwise, and so this holds for any <math>W\in\text{Im}\Phi\left(TL_{n}\left(d\right)\right)</math>, hence the operator <math>W|_l:=W|_{\mathcal{H}_{n,k,l}}</math> is well defined.

We define the following operator:

:<math>\operatorname{Tr}_{n}\left(W\right)=\frac{1}{N}\sum_{l=1}^{k-1}\lambda_{l}\operatorname{Tr}\left(W|_{l}\right)</math>

where <math>\operatorname{Tr}</math> is the usual matrix trace.

It turns out that this operator is a trace operator with the Markov property, so by the theorem stated above it has to be the Markov trace. This finishes the required reductions as it establishes that to approximate the Jones polynomial it suffices to approximate <math>\operatorname{Tr}_{n}</math>.

== The algorithm ==

 '''algorithm''' Approximate-Jones-Trace-Closure '''is'''
     '''input''' <math>B\in B_n</math> with <math>m</math> crossings
                 An integer <math>k</math>
     '''output''' a number <math>s</math> such that <math>|s-V_{B^{tr}}(e^{2\pi/k})|<1/poly(n,k,m)</math> 
                  with all but exponentially small probability
     '''repeat''' for <math>j=1</math> to <math>poly(n,m,k)</math>
         1. Pick a random <math>l\in\{1,\ldots,k-1\}</math> such that the probability
            to choose a particular <math>l</math> is proportional to <math>\lambda_l</math>
         2. Randomly pick <math>q\in Q_{n,k}</math> which ends in position <math>l</math>
         3. Using the [[Hadamard test (quantum computation)|Hadamard test]] create a random variable <math>x_j</math> with
            <math>\mathbb{E}\left[x_{j}\right]=\mathcal{Re}\left\langle \alpha\mid Q\left(B\right)\mid\beta\right\rangle </math>
     Do the same to create <math>y_j</math> with <math>\mathbb{E}\left[y_{j}\right]=\mathcal{Im}\left\langle \alpha\mid Q\left(B\right)\mid\beta\right\rangle </math>
     let <math>r</math> be the average of <math>x_j+iy_j</math>
     '''return''' <math>(-A)^{3w(B^{tr})}d^{n-1}r</math>

Note that the parameter <math>d</math> used in the algorithm depends on <math>k</math>.

The correctness of this algorithm is established by applying the [[Hoeffding's inequality|Hoeffding bound]] to <math>x_j</math> and <math>y_j</math> separately.

==Notes==
{{reflist|30em}}

==References==

# D. Aharonov, V. Jones, Z. Landau - [https://arxiv.org/abs/quant-ph/0511096 A Polynomial Quantum Algorithm for Approximating the Jones Polynomial]


{{DEFAULTSORT:Aharonov-Jones-Landau algorithm}}
[[Category:Quantum algorithms]]