Graphical representation of products of quaternion units as 90°-rotations in the planes of 4D-space, spanned by two of {1, i, j, k}. The left factor can be taken as rotated by the right. For example: in blue: 1/i-plane: 1 ⋅ i = i, i/k-plane: i ⋅ j = k in red: 1/j-plane: 1 ⋅ j = j, j/k-plane: j ⋅ i = -k

Here as he walked byon the 16th of October 1843Sir William Rowan Hamiltonin a flash of genius discoveredthe fundamental formula forquaternion multiplicationi2 = j2 = k2 = ijk = −1& cut it on a stone of this bridge

Hamilton knew that the complex numbers could be interpreted as points in a plane, and he was looking for a way to do the same for points in three-dimensional space. Points in space can be represented by their coordinates, which are triples of numbers, and for many years he had known how to add and subtract triples of numbers. However, Hamilton had been stuck on the problem of multiplication and division for a long time. He could not figure out how to calculate the quotient of the coordinates of two points in space.

The great breakthrough in quaternions finally came on Monday 16 October 1843 in Dublin, when Hamilton was on his way to the Royal Irish Academy where he was going to preside at a council meeting. As he walked along the towpath of the Royal Canal with his wife, the concepts behind quaternions were taking shape in his mind. When the answer dawned on him, Hamilton could not resist the urge to carve the formula for the quaternions,

On the following day, Hamilton wrote a letter to his friend and fellow mathematician, John T. Graves, describing the train of thought that led to his discovery. This letter was later published in a letter to a science magazine;[12] Hamilton states:

And here there dawned on me the notion that we must admit, in some sense, a fourth dimension of space for the purpose of calculating with triples ... An electric circuit seemed to close, and a spark flashed forth.[12]

Hamilton called a quadruple with these rules of multiplication a quaternion, and he devoted most of the remainder of his life to studying and teaching them. Hamilton's treatment is more geometric than the modern approach, which emphasizes quaternions' algebraic properties. He founded a school of "quaternionists", and he tried to popularize quaternions in several books. The last and longest of his books, Elements of Quaternions,[13] was 800 pages long; it was edited by his son and published shortly after his death.

After Hamilton's death, his student Peter Tait continued promoting quaternions. At this time, quaternions were a mandatory examination topic in Dublin. Topics in physics and geometry that would now be described using vectors, such as kinematics in space and Maxwell's equations, were described entirely in terms of quaternions. There was even a professional research association, the Quaternion Society, devoted to the study of quaternions and other hypercomplex number systems.

From the mid-1880s, quaternions began to be displaced by vector analysis, which had been developed by Josiah Willard Gibbs, Oliver Heaviside, and Hermann von Helmholtz. Vector analysis described the same phenomena as quaternions, so it borrowed some ideas and terminology liberally from the literature of quaternions. However, vector analysis was conceptually simpler and notationally cleaner, and eventually quaternions were relegated to a minor role in mathematics and physics. A side-effect of this transition is that Hamilton's work is difficult to comprehend for many modern readers. Hamilton's original definitions are unfamiliar and his writing style was wordy and difficult to understand.

The finding of 1924 that in quantum mechanics the spin of an electron can be described using quaternions furthered their revival; they helped to understand how rotations of electrons by 360° can be discerned from those by 720° (the “Plate trick”).[22] As of 2018[update], their use has not overtaken rotation groups.[a]

where a, b, c, d, are real numbers, and i, j, k, are symbols that can be interpreted as unit-vectors pointing along the three spatial axes. In practice, if one of a, b, c, d is 0, the corresponding term is omitted; if a, b, c, d are all zero, the quaternion is the zero quaternion, denoted 0; if one of b, c, d equals 1, the corresponding term is written simply i, j, or k.

Hamilton describes a quaternion q=a+bi+cj+dk{\displaystyle q=a+b\ \mathbf {i} +c\ \mathbf {j} +d\ \mathbf {k} }, as consisting of a scalar part and a vector part. The quaternion bi+cj+dk{\displaystyle b\ \mathbf {i} +c\ \mathbf {j} +d\ \mathbf {k} } is called the vector part (sometimes imaginary part) of q, and a is the scalar part (sometimes real part) of q. A quaternion that equals its real part (that is, its vector part is zero) is called a scalar, and is identified with the corresponding real number. That is, the real numbers are a subset of the quaternions. A quaternion that equals its vector part is called a vector quaternion.

The set of quaternions is made a 4 dimensional vector space over the real numbers, with {1,i,j,k}{\displaystyle \left\{1,\mathbf {i} ,\mathbf {j} ,\mathbf {k} \right\}} as a basis, by the componentwise addition

The real quaternions commute with all other quaternions, that is aq = qa for every quaternion q and every real quaternion a. In algebraic terminology this is to say that the field of real quaternions are the center of this quaternion algebra.

The product is first given for the basis elements (see next subsection), and then extended to all quaternions by using the distributive property and the center property of the real quaternions. The Hamilton product is not commutative, but associative, thus the quaternions form an associative algebra over the reals.

Additionally, every nonzero quaternion has an inverse with respect to the Hamilton product:

The center of a noncommutative ring is the subring of elements c such that cx = xc for every x. The center of the quaternion algebra is the subfield of real quaternions. In fact, it is a part of the definition that the real quaternions belong to the center. Conversely, if q = a + bi + cj + dk belongs to the center, then

0=iq−qi=2cij+2dik=2ck−2dj,{\displaystyle 0=iq-qi=2cij+2dik=2ck-2dj,}

and c = d = 0. A similar computation with j instead of i shows that one has also b = 0. Thus q = a is a real quaternion.

The noncommutativity of multiplication has some unexpected consequences, among them that polynomial equations over the quaternions can have more distinct solutions than the degree of the polynomial. For example, the equation z2 + 1 = 0, has infinitely many quaternion solutions, which are the quaternions z = bi + cj + dk such that b2 + c2 + d2 = 1. Thus these "roots of –1" form a unit sphere in the three-dimensional space of vector quaternions.

For two elements a1 + b1i + c1j + d1k and a2 + b2i + c2j + d2k, their product, called the Hamilton product (a1 + b1i + c1j + d1k) (a2 + b2i + c2j + d2k), is determined by the products of the basis elements and the distributive law. The distributive law makes it possible to expand the product so that it is a sum of products of basis elements. This gives the following expression:

A quaternion of the form a + 0i + 0j + 0k, where a is a real number, is called scalar, and a quaternion of the form 0 + bi + cj + dk, where b, c, and d are real numbers, and at least one of b, c or d is nonzero, is called a vector quaternion. If a + bi + cj + dk is any quaternion, then a is called its scalar part and bi + cj + dk is called its vector part. Even though every quaternion can be viewed as a vector in a four-dimensional vector space, it is common to define refer to the vector part as vectors in three-dimensional space. With this convention, a vector is the same as an element of the vector space R3.[b]

Hamilton also called vector quaternions right quaternions[24][25] and real numbers (considered as quaternions with zero vector part) scalar quaternions.

If a quaternion is divided up into a scalar part and a vector part, i.e.

Conjugation of quaternions is analogous to conjugation of complex numbers and to transposition (also known as reversal) of elements of Clifford algebras. To define it, let q=a+bi+cj+dk{\displaystyle q=a+bi+cj+dk} be a quaternion. The conjugate of q is the quaternion q∗=a−bi−cj−dk{\displaystyle q^{*}=a-bi-cj-dk}. It is denoted by q∗, q,[7]qt, or q~{\displaystyle {\tilde {q}}}. Conjugation is an involution, meaning that it is its own inverse, so conjugating an element twice returns the original element. The conjugate of a product of two quaternions is the product of the conjugates in the reverse order. That is, if p and q are quaternions, then (pq)∗ = q∗p∗, not p∗q∗.

Unlike the situation in the complex plane,
the conjugation of a quaternion can be expressed entirely with multiplication and addition:

Conjugation can be used to extract the scalar and vector parts of a quaternion. The scalar part of p is (p + p∗) / 2, and the vector part of p is (p − p∗) / 2.

The square root of the product of a quaternion with its conjugate is called its norm and is denoted ||q|| (Hamilton called this quantity the tensor of q, but this conflicts with modern meaning of "tensor"). In formula, this is expressed as follows:

This is always a non-negative real number, and it is the same as the Euclidean norm on H considered as the vector space R4. Multiplying a quaternion by a real number scales its norm by the absolute value of the number. That is, if α is real, then

Using conjugation and the norm makes it possible to define the reciprocal of a non-zero quaternion. The product of a quaternion with its reciprocal should equal 1, and the considerations above imply that the product of q{\displaystyle q} and q∗/‖q‖2{\displaystyle q^{\ast }/\left\Vert q\right\|^{2}} (in either order) is 1. So the reciprocal of q is defined to be

This makes it possible to divide two quaternions p and q in two different ways (when q is non-zero). That is, their quotient can be either p q−1 or q−1p. The notation p/q is ambiguous because it does not specify whether q divides on the left or the right.

Cayley graph of Q8. The red arrows represent multiplication on the right by i, and the green arrows represent multiplication on the right by j.

The set H of all quaternions is a vector space over the real numbers with dimension 4. (In comparison, the real numbers have dimension 1, the complex numbers have dimension 2, and the octonions have dimension 8.) Multiplication of quaternions is associative and distributes over vector addition, but it is not commutative. Therefore, the quaternions H are a non-commutative associative algebra over the real numbers. Even though H contains copies of the complex numbers, it is not an associative algebra over the complex numbers.

Because it is possible to divide quaternions, they form a division algebra. This is a structure similar to a field except for the non-commutativity of multiplication. Finite-dimensional associative division algebras over the real numbers are very rare. The Frobenius theorem states that there are exactly three: R, C, and H. The norm makes the quaternions into a normed algebra, and normed division algebras over the reals are also very rare: Hurwitz's theorem says that there are only four: R, C, H, and O (the octonions). The quaternions are also an example of a composition algebra and of a unital Banach algebra.

Because the product of any two basis vectors is plus or minus another basis vector, the set {±1, ±i, ±j, ±k} forms a group under multiplication. This non-Abelian Group is called the quaternion group and is denoted Q8.[26] The real group ring of Q8 is a ring R[Q8] which is also an eight-dimensional vector space over R. It has one basis vector for each element of Q8. The quaternions are the quotient ring of R[Q8] by the ideal generated by the elements 1 + (−1), i + (−i), j + (−j), and k + (−k). Here the first term in each of the differences is one of the basis elements 1, i, j, and k, and the second term is one of basis elements −1, −i, −j, and −k, not the additive inverses of 1, i, j, and k.

The vector part of a quaternion can be interpreted as a coordinate vector in R3, therefore the algebraic operations of the quaternions reflect the geometry of R3. Operations such as the vector dot and cross products can be defined in terms of quaternions, and this makes it possible to apply quaternion techniques wherever spatial vectors arise. A useful application of quaternions has been to interpolate the orientations of key-frames in computer graphics.[14]

For the remainder of this section, i, j, and k will denote both the three imaginary[27] basis vectors of H and a basis for R3. Notice that replacing i by −i, j by −j, and k by −k sends a vector to its additive inverse, so the additive inverse of a vector is the same as its conjugate as a quaternion. For this reason, conjugation is sometimes called the spatial inverse.

This shows that the noncommutativity of quaternion multiplication comes from the multiplication of vector quaternions. It also shows that two quaternions commute if and only if their vector parts are collinear. Hamilton[28] showed that this product computes the third vertex of a spherical triangle from two given vertices and their associated arc-lengths, which is also an algebra of points in Elliptic geometry.

Unit quaternions can be identified with rotations in R3 and were called versors by Hamilton[28]. Also see Quaternions and spatial rotation for more information about modeling three-dimensional rotations using quaternions.

Just as complex numbers can be represented as matrices, so can quaternions. There are at least two ways of representing quaternions as matrices in such a way that quaternion addition and multiplication correspond to matrix addition and matrix multiplication. One is to use 2 × 2 complex matrices, and the other is to use 4 × 4 real matrices. In each case, the representation given is one of a family of linearly related representations. In the terminology of abstract algebra, these are injectivehomomorphisms from H to the matrix ringsM(2, C) and M(4, R), respectively.

Using 2 × 2 complex matrices, the quaternion a + bi + cj + dk can be represented as

There is a strong relation between quaternion units and Pauli matrices. Obtain the eight quaternion unit matrices by taking a, b, c and d, set three of them at zero and the fourth at 1 or −1. Multiplying any two Pauli matrices always yields a quaternion unit matrix, all of them except for −1. One obtains −1 via i2 = j2 = k2 = ijk = −1. E.g. the last equality is

In fact, there exist 48 distinct representations of this form. More precisely, there are 48 sets of quadruples of matrices such that a function sending 1, i, j, and k to the matrices in the quadruple is a homomorphism, that is, it sends sums and products of quaternions to sums and products of matrices.[31]
In this representation, the conjugate of a quaternion corresponds to the transpose of the matrix. The fourth power of the norm of a quaternion is the determinant of the corresponding matrix. As with the 2 × 2 complex representation above, complex numbers can again be produced by constraining the coefficients suitably; for example, as block diagonal matrices with two 2 × 2 blocks by setting c = d = 0.

Each 4x4 matrix representation of quaternions corresponds to a multiplication table of unit quaternions. For example, the last matrix representation given above corresponds to the multiplication table

×

a

d

−b

−c

a

a

d

−b

−c

−d

−d

a

c

−b

b

b

−c

a

−d

c

c

b

d

a

which is isomorphic — through {a↦1,b↦i,c↦j,d↦k}{\displaystyle \{a\mapsto 1,b\mapsto i,c\mapsto j,d\mapsto k\}} — to

×

1

k

−i

−j

1

1

k

−i

−j

−k

−k

1

j

−i

i

i

−j

1

−k

j

j

i

k

1

Constraining any such multiplication table to have the identity in the first row and column and for the signs of the row headers to be opposite to those of the column headers, then there are 3 possible choices for the second column (ignoring sign), 2 possible choices for the third column (ignoring sign), and 1 possible choice for the fourth column (ignoring sign); that makes 6 possibilities. Then, the second column can be chosen to be either positive or negative, the third column can be chosen to be positive or negative, and the fourth column can be chosen to be positive or negative, giving 8 possibilities for the sign. Multiplying the possibilities for the letter positions and for their signs yields 48. Then replacing 1 with a, i with b, j with c, and k with d and removing the row and column headers yields a matrix representation of a + bi + cj + dk.

Quaternions are also used in one of the proofs of Lagrange's four-square theorem in number theory, which states that every nonnegative integer is the sum of four integer squares. As well as being an elegant theorem in its own right, Lagrange's four square theorem has useful applications in areas of mathematics outside number theory, such as combinatorial design theory. The quaternion-based proof uses Hurwitz quaternions, a subring of the ring of all quaternions for which there is an analog of the Euclidean algorithm.

Quaternions can be represented as pairs of complex numbers. From this perspective, quaternions are the result of applying the Cayley–Dickson construction to the complex numbers. This is a generalization of the construction of the complex numbers as pairs of real numbers.

Let C2 be a two-dimensional vector space over the complex numbers. Choose a basis consisting of two elements 1 and j. A vector in C2 can be written in terms of the basis elements 1 and j as

(a+bi)1+(c+di)j.{\displaystyle (a+bi)1+(c+di)j.\ }

If we define j2 = −1 and ij = −ji, then we can multiply two vectors using the distributive law. Writing k in place of the product ij leads to the same rules for multiplication as the usual quaternions. Therefore, the above vector of complex numbers corresponds to the quaternion a + bi + cj + dk. If we write the elements of C2 as ordered pairs and quaternions as quadruples, then the correspondence is

In the complex numbers, C, there are just two numbers, i and −i, whose square is −1 . In H there are infinitely many square roots of minus one: the quaternion solution for the square root of −1 is the unit sphere in R3. To see this, let q = a + bi + cj + dk be a quaternion, and assume that its square is −1. In terms of a, b, c, and d, this means

a2−b2−c2−d2=−1,{\displaystyle a^{2}-b^{2}-c^{2}-d^{2}=-1,}

2ab=0,{\displaystyle 2ab=0,}

2ac=0,{\displaystyle 2ac=0,}

2ad=0.{\displaystyle 2ad=0.}

To satisfy the last three equations, either a = 0 or b, c, and d are all 0. The latter is impossible because a is a real number and the first equation would imply that a2 = −1. Therefore, a = 0 and b2 + c2 + d2 = 1. In other words: a quaternion squares to −1 if and only if it is a vector quaternion with norm 1. By definition, the set of all such vectors forms the unit sphere.

Only negative real quaternions have infinitely many square roots. All others have just two (or one in the case of 0).[citation needed]

The identification of the square roots of minus one in H was given by Hamilton[32] but was frequently omitted in other texts. By 1971 the sphere was included by Sam Perlis in his three-page exposition included in Historical Topics in Algebra (page 39) published by the National Council of Teachers of Mathematics. More recently, the sphere of square roots of minus one is described in Ian R. Porteous's book Clifford Algebras and the Classical Groups (Cambridge, 1995) in proposition 8.13 on page 60.

(Note that this is not the same as qs+‖q‖⋅Uq{\displaystyle q_{s}+\lVert q\rVert \cdot \mathbf {U} q}.) The versor of the vector part of q, Uq→v{\displaystyle \mathbf {U} {\vec {q}}_{v}}, is a unit vector quaternion, so its square is −1. Therefore, it determines a copy of the complex numbers by the function

Under this function, q is the image of the complex number qs+‖q→v‖i{\displaystyle q_{s}+\lVert {\vec {q}}_{v}\rVert i}. Thus H is the union of complex planes intersecting in a common real line, where the union is taken over the sphere of square roots of minus one, bearing in mind that the same plane is associated with the antipodal points of the sphere.

The relationship of quaternions to each other within the complex subplanes of H can also be identified and expressed in terms of commutativesubrings. Specifically, since two quaternions p and q commute (i.e., pq = qp) only if they lie in the same complex subplane of H, the profile of H as a union of complex planes arises when one seeks to find all commutative subrings of the quaternion ring. This method of commutative subrings is also used to profile the split-quaternions, which as an algebra over the reals are isomorphic to 2 × 2 real matrices.

Like functions of a complex variable, functions of a quaternion variable suggest useful physical models. For example, the original electric and magnetic fields described by Maxwell were functions of a quaternion variable.

The term "conjugation", besides the meaning given above, can also mean taking an element a to rar−1 where r is some non-zero element (quaternion). All elements that are conjugate to a given element (in this sense of the word conjugate) have the same real part and the same norm of the vector part. (Thus the conjugate in the other sense is one of the conjugates in this sense.)

Thus the multiplicative group of non-zero quaternions acts by conjugation on the copy of R3 consisting of quaternions with real part equal to zero. Conjugation by a unit quaternion (a quaternion of absolute value 1) with real part cos(θ) is a rotation by an angle 2θ, the axis of the rotation being the direction of the vector part. The advantages of quaternions are:

The image of a subgroup of versors is a point group, and conversely, the preimage of a point group is a subgroup of versors. The preimage of a finite point group is called by the same name, with the prefix binary. For instance, the preimage of the icosahedral group is the binary icosahedral group.

The versors' group is isomorphic to SU(2), the group of complex unitary 2×2 matrices of determinant 1.

The Quaternions can be generalized into further algebras called quaternion algebras. Take F to be any field with characteristic different from 2, and a and b to be elements of F; a four-dimensional unitary associative algebra can be defined over F with basis 1, i, j, and ij, where i2 = a, j2 = b and ij = −ji (so (ij)2 = −ab).

Quaternion algebras are isomorphic to the algebra of 2×2 matrices over F or form division algebras over F, depending on
the choice of a and b.

The usefulness of quaternions for geometrical computations can be generalised to other dimensions by identifying the quaternions as the even part Cℓ+3,0(R) of the Clifford algebra Cℓ3,0(R). This is an associative multivector algebra built up from fundamental basis elements σ1, σ2, σ3 using the product rules

In this picture, quaternions correspond not to vectors but to bivectors – quantities with magnitude and orientations associated with particular 2D planes rather than 1D directions. The relation to complex numbers becomes clearer, too: in 2D, with two vector directions σ1 and σ2, there is only one bivector basis element σ1σ2, so only one imaginary. But in 3D, with three vector directions, there are three bivector basis elements σ1σ2, σ2σ3, σ3σ1, so three imaginaries.

This reasoning extends further. In the Clifford algebra Cℓ4,0(R), there are six bivector basis elements, since with four different basic vector directions, six different pairs and therefore six different linearly independent planes can be defined. Rotations in such spaces using these generalisations of quaternions, called rotors, can be very useful for applications involving homogeneous coordinates. But it is only in 3D that the number of basis bivectors equals the number of basis vectors, and each bivector can be identified as a pseudovector.

There are several advantages for placing quaternions in this wider setting:[34]

Rotors are a natural part of geometric algebra and easily understood as the encoding of a double reflection.

In geometric algebra, a rotor and the objects it acts on live in the same space. This eliminates the need to change representations and to encode new data structures and methods, which is traditionally required when augmenting linear algebra with quaternions.

Rotors are universally applicable to any element of the algebra, not just vectors and other quaternions, but also lines, planes, circles, spheres, rays, and so on.

In the conformal model of Euclidean geometry, rotors allow the encoding of rotation, translation and scaling in a single element of the algebra, universally acting on any element. In particular, this means that rotors can represent rotations around an arbitrary axis, whereas quaternions are limited to an axis through the origin.

The quaternions are "essentially" the only (non-trivial) central simple algebra (CSA) over the real numbers, in the sense that every CSA over the reals is Brauer equivalent to either the reals or the quaternions. Explicitly, the Brauer group of the reals consists of two classes, represented by the reals and the quaternions, where the Brauer group is the set of all CSAs, up to equivalence relation of one CSA being a matrix ring over another. By the Artin–Wedderburn theorem (specifically, Wedderburn's part), CSAs are all matrix algebras over a division algebra, and thus the quaternions are the only non-trivial division algebra over the reals.

CSAs – rings over a field, which are simple algebras (have no non-trivial 2-sided ideals, just as with fields) whose center is exactly the field – are a noncommutative analog of extension fields, and are more restrictive than general ring extensions. The fact that the quaternions are the only non-trivial CSA over the reals (up to equivalence) may be compared with the fact that the complex numbers are the only non-trivial field extension of the reals.

"I regard it as an inelegance, or imperfection, in quaternions, or rather in the state to which it has been hitherto unfolded, whenever it becomes or seems to become necessary to have recourse to x, y, z, etc." — William Rowan Hamilton [Quoted in a letter from Tait to Cayley].

"Time is said to have only one dimension, and space to have three dimensions. […] The mathematical quaternion partakes of both these elements; in technical language it may be said to be "time plus space", or "space plus time": and in this sense it has, or at least involves a reference to, four dimensions. And how the One of Time, of Space the Three, Might in the Chain of Symbols girdled be." — William Rowan Hamilton [Quoted in R.P. Graves, Life of Sir William Rowan Hamilton].

"Quaternions came from Hamilton after his really good work had been done; and, though beautifully ingenious, have been an unmixed evil to those who have touched them in any way, including Clerk Maxwell." — W. Thompson, Lord Kelvin. (1892).

"I came later to see that, as far as the vector analysis I required was concerned, the quaternion was not only not required, but was a positive evil of no inconsiderable magnitude; and that by its avoidance the establishment of vector analysis was made quite simple and its working also simplified, and that it could be conveniently harmonised with ordinary Cartesian work." — Oliver Heaviside. (1893). Electromagnetic Theory volume I, pp. 134–135. London: The Electrician Printing and Publishing Company.

"Neither matrices nor quaternions and ordinary vectors were banished from these ten [additional] chapters. For, in spite of the uncontested power of the modern Tensor Calculus, those older mathematical languages continue, in my opinion, to offer conspicuous advantages in the restricted field of special relativity. Moreover, in science as well as in every-day life, the mastery of more than one language is also precious, as it broadens our views, is conducive to criticism with regard to, and guards against hypostasy [weak-foundation] of, the matter expressed by words or mathematical symbols." — Ludwik Silberstein. (1924). Notes on preparing the second edition of his Theory of Relativity.

"... quaternions appear to exude an air of nineteenth century decay, as a rather unsuccessful species in the struggle-for-life of mathematical ideas. Mathematicians, admittedly, still keep a warm place in their hearts for the remarkable algebraic properties of quaternions but, alas, such enthusiasm means little to the harder-headed physical scientist." — Simon L. Altmann. (1986).

^A more personal view of quaternions was written by Joachim Lambek in 1995. He wrote in his essay If Hamilton had prevailed: quaternions in physics: "My own interest as a graduate student was raised by the inspiring book by Silberstein". He concluded by stating "I firmly believe that quaternions can supply a shortcut for pure mathematicians who wish to familiarize themselves with certain aspects of theoretical physics." Lambek, J. (1995). "If Hamilton had prevailed: Quaternions in physics". Math. Intelligencer. Vol. 17 no. 4. pp. 7–15. doi:10.1007/BF03024783.

^It is important to note that the vector part of a quaternion is, in truth, an "axial" vector or "pseudovector", not an ordinary or "polar" vector, as was formally proven by S.L. Altmann in Ch. 12 of his 1986 book, Rotations, Quaternions and Double Groups. A polar vector can be represented in calculations (for example, for rotation by a quaternion "similarity transform") by a pure imaginary quaternion, with no loss of information, but the two should not be confused. The axis of a "binary" (180°) rotation quaternion corresponds to the direction of the represented polar vector in such a case.

^Robert E. Bradley, Charles Edward Sandifer (2007). Leonhard Euler: life, work and legacy. p. 193. ISBN0-444-52728-1.CS1 maint: Uses authors parameter (link) . They mention Wilhelm Blaschke's claim in 1959 that "the quaternions were first identified by L. Euler in a letter to Goldbach written on 4 May 1748," and they comment that "it makes no sense whatsoever to say that Euler "identified" the quaternions in this letter ... this claim is absurd."

^Tomb Raider (1996) is often cited as the first mass-market computer game to have used quaternions to achieve smooth three-dimensional rotations. See, for example Nick Bobick (July 1998). "Rotating objects using quaternions". Game Developer.

^"Maths - Transformations using Quaternions". EuclideanSpace. A rotation of q1 followed by a rotation of q2 is equivalent to a single rotation of q2*q1. Note the reversal of order, that is, we put the first rotation on the right hand side of the multiplication.

Crowe, Michael J. (1967), A History of Vector Analysis: The Evolution of the Idea of a Vectorial System, University of Notre Dame Press. Surveys the major and minor vector systems of the 19th century (Hamilton, Möbius, Bellavitis, Clifford, Grassmann, Tait, Peirce, Maxwell, Macfarlane, MacAuley, Gibbs, Heaviside).

For molecules that can be regarded as classical rigid bodies molecular dynamics computer simulation employs quaternions. They were first introduced for this purpose by D.J. Evans, (1977), "On the Representation of Orientation Space", Mol. Phys., vol 34, p 317.