Topological Quantum Information, Khovanov Homology and the Jones Polynomial
Abstract
This paper is dedicated to Anatoly Libgober on his 60th birthday.
1 Introduction
In this paper we give a quantum statistical interpretation for the bracket polynomial state sum and for the Jones polynomial We use this quantum mechanical interpretation to give a new quantum algorithm for computing the Jones polynomial. This algorithm is useful for its conceptual simplicity and it applies to all values of the polynomial variable that lie on the unit circle in the complex plane. Letting denote the Hilbert space for this model, there is a natural unitary transformation
such that
where is a sum over basis states for The quantum algorithm comes directly from this formula via the Hadamard Test. We then show that the framework for our quantum model for the bracket polynomial is a natural setting for Khovanov homology. The Hilbert space of our model has basis in onetoone correspondence with the enhanced states of the bracket state summmation and is isomorphic with the chain complex for Khovanov homology with coefficients in the complex numbers. We show that for the Khovanov boundary operator we have the relationship Consequently, the operator acts on the Khovanov homology, and we therefore obtain a direct relationship between Khovanov homology and this quantum algorithm for the Jones polynomial. The quantum algorithm given here is inefficient, and so it remains an open problem to determine better quantum algorithms that involve both the Jones polynomial and the Khovanov homology.
The paper is organized as follows. Section 2 reviews the structure of the bracket polynomial, its state summation, the use of enhanced states and the relationship with the Jones polynomial. Section 3 describes the quantum statistical model for the bracket polynomial and refers to Section 7 (the Appendix) for the details of the Hadamard Test and the structure of the quantum algorithm. Section 4 describes the relationship of the quantum model with Khovanov homology axiomatically, without using the specific details of the Khovanov chain complex that give these axioms. In Section 5 we construct the Khovanov chain complex in detail and show how some of its properties follow uniquely from the axioms of the previous section. In Section 6 we discuss how the framework of this paper can be generalized to other situations where the Hilbert space of a quantum information system is also a chain complex. We give an example using DeRahm cohomology. In Section 7 (the Appendix) we detail the Hadamard test.
2 Bracket Polynomial and Jones Polynomial
The bracket polynomial [11] model for the Jones polynomial [8, 9, 10, 29] is usually described by the expansion
Here the small diagrams indicate parts of otherwise identical larger knot or link diagrams. The two types of smoothing (local diagram with no crossing) in this formula are said to be of type ( above) and type ( above).
One uses these equations to normalize the invariant and make a model of the Jones polynomial. In the normalized version we define
where the writhe is the sum of the oriented crossing signs for a choice of orientation of the link Since we shall not use oriented links in this paper, we refer the reader to [11] for the details about the writhe. One then has that is invariant under the Reidemeister moves (again see [11]) and the original Jones polynonmial is given by the formula
The Jones polynomial has been of great interest since its discovery in 1983 due to its relationships with statistical mechanics, due to its ability to often detect the difference between a knot and its mirror image and due to the many open problems and relationships of this invariant with other aspects of low dimensional topology.
The State Summation. In order to obtain a closed formula for the bracket, we now describe it as a state summation. Let be any unoriented link diagram. Define a state, , of to be the collection of planar loops resulting from a choice of smoothing for each crossing of There are two choices ( and ) for smoothing a given crossing, and thus there are states of a diagram with crossings. In a state we label each smoothing with or according to the convention indicated by the expansion formula for the bracket. These labels are the vertex weights of the state. There are two evaluations related to a state. The first is the product of the vertex weights, denoted The second is the number of loops in the state , denoted
Define the state summation, , by the formula
where This is the state expansion of the bracket. It is possible to rewrite this expansion in other ways. For our purposes in this paper it is more convenient to think of the loop evaluation as a sum of two loop evaluations, one giving and one giving This can be accomplished by letting each state curve carry an extra label of or We describe these enhanced states below.
Changing Variables. Letting denote the number of crossings in the diagram if we replace by and then replace by the bracket is then rewritten in the following form:
with . It is useful to use this form of the bracket state sum for the sake of the grading in the Khovanov homology (to be described below). We shall continue to refer to the smoothings labeled (or in the original bracket formulation) as smoothings.
Using Enhanced States. We now use the convention of enhanced states where an enhanced state has a label of or on each of its component loops. We then regard the value of the loop as the sum of the value of a circle labeled with a (the value is ) added to the value of a circle labeled with an (the value is We could have chosen the less neutral labels of and so that
and
since an algebra involving and naturally appears later in relation to Khovanov homology. It does no harm to take this form of labeling from the beginning.
Consider the form of the expansion of this version of the bracket polynonmial in enhanced states. We have the formula as a sum over enhanced states
where is the number of type smoothings in and , with equal to the number of loops labeled minus the number of loops labeled in the enhanced state
One advantage of the expression of the bracket polynomial via enhanced states is that it is now a sum of monomials. We shall make use of this property throughout the rest of the paper.
3 Quantum Statistics and the Jones Polynomial
We now use the enhanced state summation for the bracket polynomial with variable to give a quantum formulation of the state sum. Let be on the unit circle in the complex plane. (This is equivalent to letting the original bracket variable be on the unit circle and equivalent to letting the Jones polynmial variable be on the unit circle.) Let denote the complex vector space with orthonormal basis } where runs over the enhanced states of the diagram The vector space is the (finite dimensional) Hilbert space for our quantum formulation of the Jones polynomial. While it is customary for a Hilbert space to be written with the letter we do not follow that convention here, due to the fact that we shall soon regard as a chain complex and take its homology. One can hardly avoid using for homology.
With on the unit circle, we define a unitary transformation
by the formula
for each enhanced state Here and are as defined in the previous section of this paper.
Let
The state vector is the sum over the basis states of our Hilbert space For convenience, we do not normalize to length one in the Hilbert space We then have the
Lemma. The evaluation of the bracket polynomial is given by the following formula
Proof.
since
where is the Kronecker delta, equal to when and equal to otherwise. //
Thus the bracket polyomial evaluation is a quantum amplitude for the measurement of the state in the direction. Since can be regarded as a diagonal element of the transformation with respect to a basis containing this formula can be taken as the foundation for a quantum algorithm that computes the bracket of via the Hadamard test. See the Appendix to this paper, for a discussion of the Hadamard test and the corresponding quantum algorithm.
A few words about quantum algorithms are appropropiate at this point. In general one begins with an intial state and a unitary transformation Quantum processes are modeled by unitary transformations, and it is in principle possible to create a physical process corresponding to any given unitary transformation. In practice, one is limited by the dimensions of the spaces involved and by the fact that physical quantum states are very delicate and subject to decoherence. Nevertheless, one designs quantum algorithms at first by finding unitary operators that represent the informaton in the problem one wishes to calculate. Then the quantum part of the quantum computation is the physical process that produces the state from the initial state At this point is available for measurement. Measurement means interaction with the enviroment, and will happen in any case, but one intends a controlled circumstance in which the measurement can occur. If is a basis for the Hilbert space (here denoted by ) then we would have
a linear combination of the basis elements. On measuring with respect to this basis, each corresponds to an observable outcome and this outcome will occur with frequency The form of quantum computation consists in repeatedly preparing and running the process to find the frequency with which certain key observations occur. In the case of our algorithm we are interested in the frequency of measuring itself and this corresponds to the inner product However, a direct measurement from will yield only approximations to the absolute square of For this reason we use the more refined Hadamard test as described in the Appendix. There is much more to be said about quantum algorithms in general and about quantum algorithms for the Jones polynomial in particular. The reader can examine [1, 14, 15, 16, 26] for more information.
It is useful to formalize the bracket evaluation as a quantum amplitude. This is a direct way to give a physical interpretation of the bracket state sum and the Jones polynomial. Just how this process can be implemented physically depends upon the interpretation of the Hilbert space It is common practice in theorizing about quantum computing and quantum information to define a Hilbert space in terms of some mathematically convenient basis (such as the enhanced states of the knot or link diagram ) and leave open the possibility of a realization of the space and the quantum evolution operators that have been defined upon it. In principle any finite dimensional unitary operator can be realized by some physical system. In practice, this is the problem of constructing quantum computers and quantum information devices. It is not so easy to construct what can be done in principle, and the quantum states that are produced may be all too shortlived to produce reliable computation. Nevertheless, one has the freedom to create spaces and operators on the mathematical level and to conceptualize these in a quantum mechanical framework. The resulting structures may be realized in nature and in present or future technology. In the case of our Hilbert space associated with the bracket state sum and its corresponding unitary transformation there is rich extra structure related to Khovanov homology that we discuss in the next section. One hopes that in a (future) realization of these spaces and operators, the Khovanov homology will play a key role in quantum information related to the knot or link
There are a number of conclusions that we can draw from the formula First of all, this formulation constitutes a quantum algorithm for the computation of the bracket polynomial (and hence the Jones polynomial) at any specialization where the variable is on the unit circle. We have defined a unitary transformation and then shown that the bracket is an evaluation in the form This evaluation can be computed via the Hadamard test [24] and this gives the desired quantum algorithm. Once the unitary transformation is given as a physical construction, the algorithm will be as efficient as any application of the Hadamard test. The present algorithm requires an exponentially increasing complexity of construction for the associated unitary transformation, since the dimension of the Hilbert space is equal to the where is the number of enhanced states of the diagram . (Note that where runs over the standard bracket states, is the number of crossings in the diagram and is the number of loops in the state ) Nevertheless, it is significant that the Jones polynomial can be formulated in such a direct way in terms of a quantum algorithm. By the same token, we can take the basic result of Khovanov homology that says that the bracket is a graded Euler characteristic of the Khovanov homology as telling us that we are taking a step in the direction of a quantum algorithm for the Khovanov homology itself. This will be discussed below.
4 Khovanov Homology and a Quantum Model for the Jones Polynomial
In this section we outline how the Khovanov homology is related with our quantum model. This can be done essentially axiomatically, without giving the details of the Khovanov construction. We give these details in the next section. The outline is as follows:

There is a boundary operator defined on the Hilbert space of enhanced states of a link diagram
such that and so that if denotes the subspace of spanned by enhanced states with and then
That is, we have the formulas
and
for each enhanced state In the next section, we shall explain how the boundary operator is constructed.

Lemma. By defining as in the previous section, via
we have the following basic relationship between and the boundary operator
Proof. This follows at once from the definition of and the fact that preserves and increases to //

From this Lemma we conclude that the operator acts on the homology of We can regard as a new Hilbert space on which the unitary operator acts. In this way, the Khovanov homology and its relationship with the Jones polynomial has a natural quantum context.

For a fixed value of ,
is a subcomplex of with the boundary operator Consequently, we can speak of the homology Note that the dimension of is equal to the number of enhanced states with and Consequently, we have
Here we use the definition of the Euler characteristic of a chain complex
and the fact that the Euler characteristic of the complex is equal to the Euler characteristic of its homology. The quantum amplitude associated with the operator is given in terms of the Euler characteristics of the Khovanov homology of the link
Our reformulation of the bracket polynomial in terms of the unitary operator leads to a new viewpoint on the Khovanov homology as a representation space for the action of The bracket polynomial is then a quantum amplitude that expressing the Euler characteristics of the homology associated with this action. The decomposition of the chain complex into the parts corresponds to the eigenspace decomposition of the operator The reader will note that in this case the operator is already diagonal in the basis of enhanced states for the chain complex We regard this reformulation as a guide to further questions about the relationship of the Khovanov homology with quantum information associated with the link
As we shall see in the next section, the internal combinatorial structure of the set of enhanced states for the bracket summation leads to the Khovanov homology theory, whose graded Euler characteristic yields the bracket state sum. Thus we have a quantum statistical interpretation of the Euler characteristics of the Khovanov homology theory, and a conceptual puzzle about the nature of this relationship with the Hilbert space of that quantum theory. It is that relationship that is the subject of this paper. The unusual point about the Hilbert space is the each of its basis elements has a specific combinatorial structure that is related to the topology of the knot Thus this Hilbert space is, from the point of view of its basis elements, a form of “taking apart” of the topological structure of the knot that we are interested in studying.
Homological structure of the unitary transformation. We now prove a general result about the structure of a chain complex that is also a finite dimensional Hilbert space. Let be a chain complex over the complex numbers with boundary operator
with denoting the direct sum of all the , (for some ). Let
be a unitary operator that satisfies the equation We do not assume a second grading as occurs in the Khovanov homology. However, since is unitary, it follows [22] that there is a basis for in which is diagonal. Let denote this basis. Let denote the eigenvalue of corresponding to so that Let be the matrix element for so that
where runs over a set of basis elements so that
Lemma. With the above conventions, we have that for a basis element such that then
Proof. Note that
while
Since the conclusion of the Lemma follows from the independence of the elements in the basis for the Hilbert space. //
In this way we see that eigenvalues will propagate forward from with alternating signs according to the appearance of successive basis elements in the boundary formulas for the chain complex. Various states of affairs are possible in general, with new eigenvaluues starting at some for The simplest state of affairs would be if all the possible eigenvalues (up to multiplication by ) for occurred in so that
where runs over all the distinct eigenvalues of restricted to and is spanned by all in with Let us take the further assumption that for each as above, the subcomplexes
have as their direct sum. With this assumption about the chain complex, define as before, with running over the whole basis for Then it follows just as in the beginning of this section that
Here denotes the Euler characteristic of the homology. The point is, that this formula for takes exactly the form we had for the special case of Khovanov homology (with eigenvalues ), but here the formula occurs just in terms of the eigenspace decomposition of the unitary transformation in relation to the chain complex. Clearly there is more work to be done here and we will return to it in a subsequent paper.
Remark on the density matrix. Given the state we can define the density matrix
With this definition it is immediate that
where denotes the trace of a matrix Thus we can restate the form of our result about Euler characteristics as
In searching for an interpretation of the Khovanov complex in this quantum context it is useful to use this reformulation. For the bracket we have
Remark on physics. We are taking the whole Hilbert space as the state space of a physical system. If the system were a classical system then the energetic states of the classical system would correspond to the basis we have chosen for the Hilbert space. This is in line with the concepts of quantum mechanics, since the classical states are then in correspondence with a basis for measurement. In quantum context we think of the elements in the Hilbert space as corresponding to superpositions of possible measurements. All this is then transposed to topological configurations for the knot or link. But the new ingredient is the one from topology that will make the states of the knot into a chain complex and lead to homological computations. Since we are using enhanced states, each state is a generator for the chain complex. Thus we can regard itself as the chain complex (over the complex numbers) and add in extra grading structure to boot. There is also the fine structure of the underlying state circles, but this is not seen by the Hilbert space itself and the chain complex structure can certainly be written just in terms of the basis of the Hilbert space. We want to know what is the relationship between the unitary transformation and the homological structure. This relation is given by the way the grading works in relation to for each state: the constancy of the quantum grading under the action of the differential. It is this that makes a graded Euler characteristic for the homology. In this way the unitary transformation is linked with the structure of state transitions that govern the homology. At the end of this section we indicate a more general approach to this pattern.
The key conceptual issue is the construction of a Hilbert space whose basis is the set of observable states of a physical system. We can always do this for a statistical mechanical system. Let denote such states and denote the energy of a state (observable and real). Define a unitary transformation (where denotes time) by
Define the amplitude
Note that we have the formula
where just as before. Now can be regarded as the quantum evolution of the state of the system and we have written this amplitude in analogy with a partition function for the statistical mechanics model.
One more conceptual problem: We consider an amplitude of the form
where is a unit complex number. The question is, how does this fit the pattern
So we want
which gives
Evolution from time is out of the question. Probably the best interpretation is to think of the evolution starting at
In this point of view, the Hilbert space for expressing the bracket polynomial as a quantum statistical amplitude is quite naturally the chain complex for Khovanov homology with complex coefficients, and the unitary transformation that is the structure of the bracket polynomial acts on the homology of this chain complex. This means that the homology classes contain information preserved by the quantum process that underlies the bracket polynomial. We would like to exploit this direct relationship between the quantum model and the Khovanov homology to obtain deeper information about the relationship of topology and quantum information theory, and we would like to use this relationship to probe the properties of these topological invariants.
5 Background on Khovanov Homology
In this section, we describe Khovanov homology along the lines of [19, 2], and we tell the story so that the gradings and the structure of the differential emerge in a natural way. This approach to motivating the Khovanov homology uses elements of Khovanov’s original approach, Viro’s use of enhanced states for the bracket polynomial [28], and BarNatan’s emphasis on tangle cobordisms [3]. We use similar considerations in our paper [17].
Two key motivating ideas are involved in finding the Khovanov invariant. First of all, one would like to categorify a link polynomial such as There are many meanings to the term categorify, but here the quest is to find a way to express the link polynomial as a graded Euler characteristic for some homology theory associated with
We will use the bracket polynomial and its enhanced states as described in the previous sections of this paper. To see how the Khovanov grading arises, consider the form of the expansion of this version of the bracket polynomial in enhanced states. We have the formula as a sum over enhanced states
where is the number of type smoothings in , is the number of loops in labeled minus the number of loops labeled and . This can be rewritten in the following form:
where we define to be the linear span (over the complex numbers for the purpose of this paper, but over the integers or the integers modulo two for other contexts) of the set of enhanced states with and Then the number of such states is the dimension
We would like to have a bigraded complex composed of the with a differential
The differential should increase the homological grading by and preserve the quantum grading Then we could write
where is the Euler characteristic of the subcomplex for a fixed value of
This formula would constitute a categorification of the bracket polynomial. Below, we shall see how the original Khovanov differential is uniquely determined by the restriction that for each enhanced state . Since is preserved by the differential, these subcomplexes have their own Euler characteristics and homology. We have
where denotes the homology of the complex . We can write
The last formula expresses the bracket polynomial as a graded Euler characteristic of a homology theory associated with the enhanced states of the bracket state summation. This is the categorification of the bracket polynomial. Khovanov proves that this homology theory is an invariant of knots and links (via the Reidemeister moves of Figure 1), creating a new and stronger invariant than the original Jones polynomial.
We will construct the differential in this complex first for mod coefficients. The differential is based on regarding two states as adjacent if one differs from the other by a single smoothing at some site. Thus if denotes a pair consisting in an enhanced state and site of that state with of type , then we consider all enhanced states obtained from by smoothing at and relabeling only those loops that are affected by the resmoothing. Call this set of enhanced states Then we shall define the partial differential as a sum over certain elements in and the differential by the formula
with the sum over all type sites in It then remains to see what are the possibilities for so that is preserved.
Note that if , then Thus
From this we conclude that if and only if Recall that
where is the number of loops in labeled is the number of loops labeled (same as labeled with ) and .
Proposition. The partial differentials are uniquely determined by the condition that for all involved in the action of the partial differential on the enhanced state This unique form of the partial differential can be described by the following structures of multiplication and comultiplication on the algebra A = where for mod2 coefficients, or for integral coefficients.

The element is a multiplicative unit and

and
These rules describe the local relabeling process for loops in a state. Multiplication corresponds to the case where two loops merge to a single loop, while comultiplication corresponds to the case where one loop bifurcates into two loops.
Proof. Using the above description of the differential, suppose that there are two loops at that merge in the smoothing. If both loops are labeled in then the local contribution to is Let denote a smoothing in In order for the local contribution to become , we see that the merged loop must be labeled . Similarly if the two loops are labeled and then the merged loop must be labeled so that the local contribution for goes from to Finally, if the two loops are labeled and then there is no label available for a single loop that will give so we define to be zero in this case. We can summarize the result by saying that there is a multiplicative structure such that and this multiplication describes the structure of the partial differential when two loops merge. Since this is the multiplicative structure of the algebra we take this algebra as summarizing the differential.
Now consider the case where has a single loop at the site Smoothing produces two loops. If the single loop is labeled then we must label each of the two loops by in order to make decrease by . If the single loop is labeled then we can label the two loops by and in either order. In this second case we take the partial differential of to be the sum of these two labeled states. This structure can be described by taking a coproduct structure with and We now have the algebra with product and coproduct describing the differential completely. This completes the proof. //
Partial differentials are defined on each enhanced state and a site of type in that state. We consider states obtained from the given state by smoothing the given site . The result of smoothing is to produce a new state with one more site of type than Forming from we either amalgamate two loops to a single loop at , or we divide a loop at into two distinct loops. In the case of amalgamation, the new state acquires the label on the amalgamated circle that is the product of the labels on the two circles that are its ancestors in . This case of the partial differential is described by the multiplication in the algebra. If one circle becomes two circles, then we apply the coproduct. Thus if the circle is labeled , then the resultant two circles are each labeled corresponding to . If the orginal circle is labeled then we take the partial boundary to be a sum of two enhanced states with labels and in one case, and labels and in the other case, on the respective circles. This corresponds to Modulo two, the boundary of an enhanced state is the sum, over all sites of type in the state, of the partial boundaries at these sites. It is not hard to verify directly that the square of the boundary mapping is zero (this is the identity of mixed partials!) and that it behaves as advertised, keeping constant. There is more to say about the nature of this construction with respect to Frobenius algebras and tangle cobordisms. In Figures 2,3 and 4 we illustrate how the partial boundaries can be conceptualized in terms of surface cobordisms. The equality of mixed partials corresponds to topological equivalence of the corresponding surface cobordisms, and to the relationships between Frobenius algebras and the surface cobordism category. In particular, in Figure 4 we show how in a key case of two sites (labeled 1 and 2 in that Figure) the two orders of partial boundary are
and
In the Frobenius algebra we have the identity
Thus the Frobenius algebra implies the identity of the mixed partials. Furthermore, in Figure 3 we see that this identity corresponds to the topological equivalence of cobordisms under an exchange of saddle points. There is more to say about all of this, but we will stop here. The proof of invariance of Khovanov homology with respect to the Reidemeister moves (respecting grading changes) will not be given here. See [19, 2, 3]. It is remarkable that this version of Khovanov homology is uniquely specified by natural ideas about adjacency of states in the bracket polynomial.
Remark on Integral Differentials. Choose an ordering for the crossings in the link diagram and denote them by Let be any enhanced state of and let denote the chain obtained from by applying a partial boundary at the th site of If the th site is a smoothing of type , then If the th site is a smoothing of type , then is given by the rules discussed above (with the same signs). The compatibility conditions that we have discussed show that partials commute in the sense that for all and One then defines signed boundary formulas in the usual way of algebraic topology. One way to think of this regards the complex as the analogue of a complex in de Rahm cohomology. Let be a formal basis for a Grassmann algebra so that Starting with enhanced states in (that is, state with all type smoothings) Define formally, and regard as identical with as we have previously regarded it in In general, given an enhanced state in with smoothings at locations we represent this chain as and define
just as in a de Rahm complex. The Grassmann algebra automatically computes the correct signs in the chain complex, and this boundary formula gives the original boundary formula when we take coefficients modulo two. Note, that in this formalism, partial differentials of enhanced states with a smoothing at the site are zero due to the fact that in the Grassmann algebra. There is more to discuss about the use of Grassmann algebra in this context. For example, this approach clarifies parts of the construction in [18].
It of interest to examine this analogy between the Khovanov (co)homology and de Rahm cohomology. In that analogy the enhanced states correspond to the differentiable functions on a manifold. The Khovanov complex is generated by elements of the form where the enhanced state has smoothings at exactly the sites If we were to follow the analogy with de Rahm cohomology literally, we would define a new complex where is generated by elements where is any enhanced state of the link The partial boundaries are defined in the same way as before and the global boundary formula is just as we have written it above. This gives a new chain complex associated with the link Whether its homology contains new topological information about the link will be the subject of a subsequent paper.
A further remark on de Rham cohomology. There is another deep relation with the de Rham complex: In [21] it was observed that Khovanov homology is related to Hochschild homology and Hochschild homology is thought to be an algebraic version of de Rham chain complex (cyclic cohomology corresponds to de Rham cohomology), compare [23].
6 Other Homological States
The formalism that we have pursued in this paper to relate Khovanov homology and the Jones polynomial with quantum statistics can be generalized to apply to other situations. It is possible to associate a finite dimensional Hilbert space that is also a chain complex (or cochain complex) to topological structures other than knots and links. To formulate this, let denote the topological structure and denote the associated linear space endowed with a boundary operator with We want to consider situations where there is a unitary operator such that so that induces a unitary action on the homology of with respect to
For example, let be a differentiable manifold and denote the DeRham complex of over the complex numbers. Then for a differential form of the type in local coordinates and a wedge product of a subset of we have
Here is the differential for the DeRham complex. Then has basis the set of where with We could achieve if is a very simple unitary operator (e.g. multiplication by phases that do not depend on the coordinates ) but in general it will be an interesting problem to determine all unitary operators with this property.
Even in the case of Khovanov homology, we could keep the homology theory fixed and ask for other unitary operators that satisfy Knowing other examples of such operators would shed light on the nature of Khovanov homology from the point of view of quantum statistics.
7 Appendix  The Hadamard Test
In order to make a quantum computation of the trace of a unitary matrix , one can use the Hadamard test to obtain the diagonal matrix elements of The trace is then the sum of these matrix elements as runs over an orthonormal basis for the vector space. In the application to the algorithm described here for the Jones polynomial it is only necessary to compute one number of the form The Hadamard test proceeds as follows.
We first obtain
as an expectation by applying the Hadamard gate
to the first qubit of
Here denotes controlled acting as when the control bit is and the identity mapping when the control bit is We measure the expectation for the first qubit of the resulting state
This expectation is
In Figure 5, we illustrate this computation with a diagram that indicates the structure of the test with parallel lines corresponding to tensor products of the single qubit space (with three lines chosen for illustration as the size of ). The extra tensor factor is indicated on the top line with the Hadamard matrix indicated by a box and the control of indicated by a circle with a vertical control line extending down to the box. The half circle on the top line on the right stands for the measurement of that line that is used for the computation. Thus Figure 5 represents a circuit diagram for the quantum computation of and hence the quantum computation of when is taken to be the unitary tranformation corresponding to the bracket polynomial, as discussed in previous sections of this paper.
The imaginary part is obtained by applying the same procedure to
This is the method used in [1, 14, 15], and the reader may wish to contemplate its efficiency in the context of this simple model. Note that the Hadamard test enables this quantum computation to estimate the trace of any unitary matrix by repeated trials that estimate individual matrix entries
References
 D. Aharonov, V. Jones, Z. Landau, A polynomial quantum algorithm for approximating the Jones polynomial, quantph/0511096.
 Bar–Natan, D. (2002), On Khovanov’s categorification of the Jones polynomial, Algebraic and Geometric Topology, 2(16), pp. 337–370.
 BarNatan, D. (2005), Khovanov’s homology for tangles and cobordisms, Geometry and Topology, 933, pp. 14651499. arXiv:mat.GT0410495
 R.J. Baxter. Exactly Solved Models in Statistical Mechanics. Acad. Press (1982).
 R. Dijkgraaf, D. Orland and S. Reffert, Dimer models, free fermions and super quantum mechanics, arXiv:0705.1645.
 S. Gukov, Surface operators and knot homologies, ArXiv:0706.2369.
 L.HelmeGuizon, J.H.Przytycki, Y.Rong, Torsion in Graph Homology, Fundamenta Mathematicae, 190; June 2006, 139–177.
 V.F.R. Jones, A polynomial invariant for links via von Neumann algebras, Bull. Amer. Math. Soc. 129 (1985), 103–112.
 V.F.R.Jones. Hecke algebra representations of braid groups and link polynomials. Ann. of Math. 126 (1987), pp. 335338.
 V.F.R.Jones. On knot invariants related to some statistical mechanics models. Pacific J. Math., vol. 137, no. 2 (1989), pp. 311334.
 L.H. Kauffman, State models and the Jones polynomial, Topology 26 (1987), 395–407.
 L.H. Kauffman, Statistical mechanics and the Jones polynomial, AMS Contemp. Math. Series 78 (1989), 263–297.
 L.H. Kauffman, Knots and Physics, World Scientific Publishers (1991), Second Edition (1993), Third Edition (2002).
 L.H. Kauffman, Quantum computing and the Jones polynomial, in Quantum Computation and Information, S. Lomonaco, Jr. (ed.), AMS CONM/305, 2002, pp. 101–137. math.QA/0105255
 L. H. Kauffman and S. Lomonaco Jr., A Threestranded quantum algorithm for the Jones polynonmial, in “Quantum Information and Quantum Computation V”, Proceedings of Spie, April 2007, edited by E.J. Donkor, A.R. Pirich and H.E. Brandt, pp. 65730T117, Intl Soc. Opt. Eng.
 S. J. Lomonaco Jr. and L. H. Kauffman, Quantum Knots and Mosaics, Journal of Quantum Information Processing, Vol. 7, Nos. 23, (2008), pp. 85  115. arxiv.org/abs/0805.0339
 H.A. Dye, L.H. Kauffman, V.O. Manturov, On two categorifications of the arrow polynomial for virtual knots, arXiv:0906.3408.
 V.O. Manturov, Khovanov homology for virtual links with arbitrary coefficients, math.GT/0601152.
 Khovanov, M. (1997), A categorification of the Jones polynomial, Duke Math. J,101 (3), pp.359426.
 Lee, T. D., Yang, C. N. (1952), Statistical theory of equations of state and phase transitions II., lattice gas and Ising model, Physical Review Letters, 87, 410419.
 J. H. Przytycki, When the theories meet: Khovanov homology as Hochschild homology of links, arXiv:math.GT/0509334
 S. Lang, Algebra, SpringerVerlag, New York, Inc. (2002).
 JL. Loday, Cyclic Homology, Grund. Math. Wissen. Band 301, SpringerVerlag, Berlin, 1992 (second edition, 1998).
 M. A. Nielsen and I. L. Chuang, “Quantum Computation and Quantum Information”, Cambridge University Press (2000).
 E.F. JassoHernandez and Y. Rong, A categorification of the Tutte polynomial, Alg. Geom. Topology 6 2006, 20312049.
 Raimund Marx, Amr Fahmy, Louis Kauffman, Samuel Lomonaco, Andreas Spörl, Nikolas Pomplun, John Myers, Steffen J. Glaser, NMR Quantum Calculations of the Jones Polynomial, arXiv:0909.1080.
 M. Stosic, Categorification of the dichromatic polynomial for graphs, arXiv:Math/0504239v2, JKTR, Vol. 17, No. 1, (2008), pp. 3145.
 O. Viro (2004), Khovanov homology, its definitions and ramifications, Fund. Math., 184 (2004), pp. 317342.
 E. Witten. Quantum Field Theory and the Jones Polynomial. Comm. in Math. Phys. Vol. 121 (1989), 351399.