Skip to content

Tensor Products

The tensor product V⊗WV\otimes W converts bilinear dependence on two vector spaces into ordinary linear dependence on one larger vector space. Its distinguished elements are the simple tensors v⊗wv\otimes w, but the construction also contains arbitrary linear combinations of simple tensors.

For finite-dimensional spaces,

dim⁡(V⊗W)=dim⁡V dim⁡W.\dim(V\otimes W) =\dim V\,\dim W.

That multiplicative dimension is the first clue that a tensor product is not an ordered pair or a direct sum. This page develops the mathematical construction. The physical rule assigning tensor-product Hilbert spaces to distinguishable composite systems belongs to Tensor Products of Hilbert Spaces.

Suppose a construction should combine v∈Vv\in V and w∈Ww\in W while remaining linear in either input when the other is fixed. It must obey

(v+v′)⊗w=v⊗w+v′⊗w,v⊗(w+w′)=v⊗w+v⊗w′,(αv)⊗w=α(v⊗w)=v⊗(αw).\begin{aligned} (v+v')\otimes w &=v\otimes w+v'\otimes w,\\ v\otimes(w+w') &=v\otimes w+v\otimes w',\\ (\alpha v)\otimes w &=\alpha(v\otimes w)\\ &=v\otimes(\alpha w). \end{aligned}

These relations have immediate consequences:

0⊗w=0,v⊗0=0.0\otimes w=0, \qquad v\otimes0=0.

They also show that a scalar can be moved from one factor to the other:

(αv)⊗w=α(v⊗w)=v⊗(αw).(\alpha v)\otimes w =\alpha(v\otimes w) =v\otimes(\alpha w).

Consequently, the factors of a nonzero simple tensor are not unique. For α≠0\alpha\ne0,

(αv)⊗(α−1w)=v⊗w.(\alpha v)\otimes(\alpha^{-1}w) =v\otimes w.

The symbol ⊗\otimes therefore does not package two vectors as a literal ordered pair. It records their bilinear combination subject to these relations.

One concrete construction begins with the free vector space generated by formal symbols [v,w][v,w] for all (v,w)∈V×W(v,w)\in V\times W. Quotient by the subspace generated by the bilinearity relations

[v+v′,w]−[v,w]−[v′,w],[v,w+w′]−[v,w]−[v,w′],[αv,w]−α[v,w],[v,αw]−α[v,w].\begin{aligned} [v+v',w]-[v,w]-[v',w],\\ [v,w+w']-[v,w]-[v,w'],\\ [\alpha v,w]-\alpha[v,w],\\ [v,\alpha w]-\alpha[v,w]. \end{aligned}

The equivalence class of [v,w][v,w] is denoted v⊗wv\otimes w. This quotient is the algebraic tensor product V⊗WV\otimes W.

The construction is characterized without reference to formal symbols by its universal property. Let

τ:V×W⟶V⊗W,τ(v,w)=v⊗w.\begin{aligned} \tau&:V\times W\longrightarrow V\otimes W,\\ \tau(v,w)&=v\otimes w. \end{aligned}

For every bilinear map β:V×W→X\beta:V\times W\to X, there is a unique linear map β~:V⊗W→X\widetilde\beta:V\otimes W\to X such that

β=β~∘τ.\beta=\widetilde\beta\circ\tau.

Equivalently,

β~(v⊗w)=β(v,w)\widetilde\beta(v\otimes w)=\beta(v,w)

on simple tensors, and linearity determines its value on every tensor. This property is what makes the tensor product canonical: any two constructions with it are uniquely isomorphic in a way that preserves simple tensors.

Let

E=(e1,…,em),F=(f1,…,fn)\mathcal E=(e_1,\ldots,e_m), \qquad \mathcal F=(f_1,\ldots,f_n)

be bases of VV and WW. Then

E⊗F={ei⊗fj},1≤i≤m,1≤j≤n.\begin{aligned} \mathcal E\otimes\mathcal F &=\{e_i\otimes f_j\},\\ 1\le i&\le m, \qquad 1\le j\le n. \end{aligned}

is a basis of V⊗WV\otimes W. It has mnmn elements, so

dim⁡(V⊗W)=mn.\dim(V\otimes W)=mn.

If

v=∑ixiei,w=∑jyjfj,v=\sum_i x_i e_i, \qquad w=\sum_j y_j f_j,

bilinearity gives

v⊗w=∑i,jxiyj ei⊗fj.v\otimes w =\sum_{i,j}x_i y_j\, e_i\otimes f_j.

Every tensor T∈V⊗WT\in V\otimes W has a unique coordinate expansion

T=∑i,jcij ei⊗fj.T=\sum_{i,j}c_{ij}\, e_i\otimes f_j.

The coefficients cijc_{ij} depend on the chosen product basis, just as vector coordinates depend on a basis. The tensor itself does not.

A nonzero tensor is simple, decomposable, or rank one when it can be written as v⊗wv\otimes w. Most tensors are not simple.

In chosen bases, collect the coefficients of TT into the m×nm\times n matrix

CT=(cij).C_T=(c_{ij}).

If T=v⊗wT=v\otimes w, then

cij=xiyj,CT=xyT,c_{ij}=x_i y_j, \qquad C_T=xy^T,

so CTC_T has matrix rank one. Conversely, every rank-one coefficient matrix factors as xyTxy^T and therefore defines a simple tensor. Thus, for nonzero TT,

T is simple⟺rank⁡CT=1.T\text{ is simple} \quad\Longleftrightarrow\quad \operatorname{rank}C_T=1.

For example,

T=e1⊗f1+2e1⊗f2+3e2⊗f1+6e2⊗f2\begin{aligned} T &=e_1\otimes f_1 +2e_1\otimes f_2\\ &\quad +3e_2\otimes f_1 +6e_2\otimes f_2 \end{aligned}

has coefficient matrix

(1236)=(13)(12).\begin{pmatrix} 1&2\\ 3&6 \end{pmatrix} = \begin{pmatrix} 1\\3 \end{pmatrix} \begin{pmatrix} 1&2 \end{pmatrix}.

Therefore

T=(e1+3e2)⊗(f1+2f2).T=(e_1+3e_2)\otimes(f_1+2f_2).

By contrast, e1⊗f1+e2⊗f2e_1\otimes f_1+e_2\otimes f_2 has coefficient matrix I2I_2 and is not simple. The minimum number of simple tensors needed in a bipartite decomposition equals the rank of CTC_T. This is the entry point to Singular Value Decomposition and Schmidt Decomposition as Linear Algebra.

If VV and WW are finite-dimensional inner-product spaces, define the inner product first on simple tensors by

⟨v1⊗w1,v2⊗w2⟩=⟨v1,v2⟩V⟨w1,w2⟩W,\langle v_1\otimes w_1, v_2\otimes w_2 \rangle = \langle v_1,v_2\rangle_V \langle w_1,w_2\rangle_W,

and extend sesquilinearly. In particular,

∥v⊗w∥=∥v∥ ∥w∥.\lVert v\otimes w\rVert =\lVert v\rVert\,\lVert w\rVert.

If {ei}\{e_i\} and {fj}\{f_j\} are orthonormal, then

⟨ei⊗fj,ek⊗fℓ⟩=δikδjℓ.\langle e_i\otimes f_j, e_k\otimes f_\ell \rangle =\delta_{ik}\delta_{j\ell}.

The product basis is therefore orthonormal. For

T=∑i,jcij ei⊗fj,T=\sum_{i,j}c_{ij}\, e_i\otimes f_j,

one has

∥T∥2=∑i,j∣cij∣2.\lVert T\rVert^2 =\sum_{i,j}\lvert c_{ij}\rvert^2.

In finite dimensions no further completion is needed. In infinite dimensions, the algebraic tensor product is generally incomplete in this norm; its completion is the Hilbert-space tensor product.

Let

A:V→V′,B:W→W′A:V\to V', \qquad B:W\to W'

be linear maps. Their tensor product is the unique linear map

A⊗B:V⊗W⟶V′⊗W′A\otimes B: V\otimes W\longrightarrow V'\otimes W'

defined by

(A⊗B)(v⊗w)=Av⊗Bw.(A\otimes B)(v\otimes w) =Av\otimes Bw.

Well-definedness follows from the universal property because (v,w)↦Av⊗Bw(v,w)\mapsto Av\otimes Bw is bilinear.

Compatible tensor-product maps obey

(A⊗B)(C⊗D)=(AC)⊗(BD),(A⊗B)†=A†⊗B†,IV⊗IW=IV⊗W.\begin{aligned} (A\otimes B)(C\otimes D) &=(AC)\otimes(BD),\\ (A\otimes B)^\dagger &=A^\dagger\otimes B^\dagger,\\ I_V\otimes I_W &=I_{V\otimes W}. \end{aligned}

If AA and BB are invertible, unitary, Hermitian, or positive, then A⊗BA\otimes B inherits the corresponding property, with the qualification that the tensor product of two Hermitian operators is Hermitian but need not be positive unless both spectra have compatible signs.

Local-factor maps commute:

(A⊗IW)(IV⊗B)=A⊗B,(IV⊗B)(A⊗IW)=A⊗B.\begin{aligned} (A\otimes I_W)(I_V\otimes B) &=A\otimes B,\\ (I_V\otimes B)(A\otimes I_W) &=A\otimes B. \end{aligned}

This identity is algebraic; it does not assert that all physical operations on separated subsystems can be implemented independently.

If Avi=aiviAv_i=a_i v_i and Bwj=bjwjBw_j=b_jw_j, then

(A⊗B)(vi⊗wj)=aibj(vi⊗wj).(A\otimes B)(v_i\otimes w_j) =a_i b_j(v_i\otimes w_j).

When eigenbases exist, the eigenvalues of A⊗BA\otimes B are all products aibja_i b_j, counted with multiplicity. In finite dimensions,

rank⁡(A⊗B)=rank⁡A rank⁡B,tr⁡(A⊗B)=(tr⁡A)(tr⁡B).\begin{aligned} \operatorname{rank}(A\otimes B) &=\operatorname{rank}A\, \operatorname{rank}B,\\ \operatorname{tr}(A\otimes B) &=(\operatorname{tr}A) (\operatorname{tr}B). \end{aligned}

If AA is m×mm\times m and BB is n×nn\times n, then

det⁡(A⊗B)=(det⁡A)n(det⁡B)m.\det(A\otimes B) =(\det A)^n(\det B)^m.

These formulas include algebraic multiplicities. They are particularly useful for checking symbolic and numerical constructions.

Choose ordered bases of VV and WW, then choose an ordering of the product basis. With the lexicographic convention

(e1⊗f1,…,e1⊗fn,e2⊗f1,…,em⊗fn),\left( \begin{gathered} e_1\otimes f_1,\ldots,e_1\otimes f_n,\\ e_2\otimes f_1,\ldots,e_m\otimes f_n \end{gathered} \right),

the matrix of A⊗BA\otimes B is the Kronecker product

[A⊗B]=(a11B⋯a1mB⋮⋱⋮am1B⋯ammB).[A\otimes B] = \begin{pmatrix} a_{11}B&\cdots&a_{1m}B\\ \vdots&\ddots&\vdots\\ a_{m1}B&\cdots&a_{mm}B \end{pmatrix}.

For example, if

A=(120−1),B=(0110),A= \begin{pmatrix} 1&2\\ 0&-1 \end{pmatrix}, \qquad B= \begin{pmatrix} 0&1\\ 1&0 \end{pmatrix},

then

A⊗B=(01021020000−100−10).A\otimes B = \begin{pmatrix} 0&1&0&2\\ 1&0&2&0\\ 0&0&0&-1\\ 0&0&-1&0 \end{pmatrix}.

Changing the product-basis ordering conjugates this matrix by a permutation matrix. The abstract operator is unchanged, but coordinate arrays and code must use the same convention.

Tensor products are associative and symmetric up to canonical isomorphism, not literal equality. The associator sends

(u⊗v)⊗w⟼u⊗(v⊗w),(u\otimes v)\otimes w \longmapsto u\otimes(v\otimes w),

and the swap map sends

SV,W(v⊗w)=w⊗v.S_{V,W}(v\otimes w) =w\otimes v.

These maps justify suppressing parentheses when factor order is otherwise clear. They do not justify silently reversing factors. In coordinates, the swap is represented by a nontrivial permutation matrix, often called a commutation matrix.

Tensor products distribute over direct sums:

(V⊕W)⊗X≅(V⊗X)⊕(W⊗X).(V\oplus W)\otimes X \cong (V\otimes X)\oplus(W\otimes X).

But V⊗WV\otimes W is not generally isomorphic to V⊕WV\oplus W: their dimensions are multiplicative and additive, respectively. See Direct Sums.

In finite dimensions there are natural isomorphisms

(V⊗W)∗≅V∗⊗W∗,Hom⁡(V,W)≅W⊗V∗.\begin{aligned} (V\otimes W)^* &\cong V^*\otimes W^*,\\ \operatorname{Hom}(V,W) &\cong W\otimes V^*. \end{aligned}

Under the second isomorphism, a simple tensor w⊗φw\otimes\varphi corresponds to the rank-one map

v⟼φ(v)w.v\longmapsto\varphi(v)w.

This identification explains why outer products, rank-one operators, and two-index coefficient arrays are manifestations of the same tensor structure. In infinite dimensions the analogous identifications require topological qualifications and are not automatic for arbitrary bounded operators.

When factors are labeled AA and BB, write

HA⊗HB\mathcal H_A\otimes\mathcal H_B

and keep labels on ambiguous vectors and identities:

∣i⟩A⊗∣j⟩B,XA⊗IB.\lvert i\rangle_A\otimes\lvert j\rangle_B, \qquad X_A\otimes I_B.

Compact notation such as

∣ij⟩≡∣i⟩A⊗∣j⟩B\lvert ij\rangle \equiv \lvert i\rangle_A\otimes\lvert j\rangle_B

is safe only after the factor order has been declared. Likewise, writing XAX_A as shorthand for XA⊗IBX_A\otimes I_B is convenient only when the ambient space is unambiguous.

The notation is mathematical. The statement that a physical composite has state space HA⊗HB\mathcal H_A\otimes\mathcal H_B is a composition postulate; see Tensor Products in Core Formalism. Site-wide ordering conventions are collected in Tensor-Product Ordering.

Tensor-product dimensions grow multiplicatively. If kk factors each have dimension dd, the full dimension is dkd^k, and a dense operator has d2kd^{2k} entries. This growth is structural, not an implementation accident.

For reliable computation:

  • record the factor order and product-basis order explicitly;
  • distinguish row-major and column-major reshaping conventions;
  • test a Kronecker implementation on basis vectors;
  • exploit sparse and structured tensor products;
  • apply local operators by reshaping and contracting when possible instead of forming the full matrix;
  • compare abstract identities before comparing flattened arrays;
  • track whether complex conjugation is required by a dual or adjoint.

A reshape is not itself a basis-independent operation. It becomes meaningful only after dimensions and index ordering have been specified.

  • Treating v⊗wv\otimes w as an ordered pair (v,w)(v,w).
  • Assuming every tensor is simple.
  • Expanding (v+v′)⊗(w+w′)(v+v')\otimes(w+w') but omitting the cross terms.
  • Forgetting that scalar factors can move between tensor factors.
  • Confusing the abstract tensor product with one Kronecker coordinate matrix.
  • Reversing factor order without applying the swap isomorphism.
  • Writing AA on a composite space without the required identity factors.
  • Comparing flattened arrays produced by different basis-order conventions.
  • Assuming a direct sum and a tensor product describe the same composition.
  • Using the finite-dimensional dual-space identities without topological qualifications in infinite dimensions.
  1. Let VV have basis (e1,e2)(e_1,e_2) and WW have basis (f1,f2,f3)(f_1,f_2,f_3). Expand

    (2e1−e2)⊗(f1+3f3)(2e_1-e_2)\otimes(f_1+3f_3)

    in the product basis and give its coefficient matrix.

Solution

Set X=(2e1−e2)⊗(f1+3f3)X=(2e_1-e_2)\otimes(f_1+3f_3). Bilinearity gives

X=2e1⊗f1+6e1⊗f3−e2⊗f1−3e2⊗f3.\begin{aligned} X &=2e_1\otimes f_1 +6e_1\otimes f_3\\ &\quad -e_2\otimes f_1 -3e_2\otimes f_3. \end{aligned}

With rows indexed by (e1,e2)(e_1,e_2) and columns by (f1,f2,f3)(f_1,f_2,f_3), the coefficient matrix is

(206−10−3).\begin{pmatrix} 2&0&6\\ -1&0&-3 \end{pmatrix}.

It has rank one, as expected for a simple tensor.

  1. Determine which of the following tensors are simple:

    T1=e1⊗f1+e1⊗f2+2e2⊗f1+2e2⊗f2,T2=e1⊗f1+e2⊗f2.\begin{aligned} T_1 &=e_1\otimes f_1 +e_1\otimes f_2\\ &\quad +2e_2\otimes f_1 +2e_2\otimes f_2,\\ T_2 &=e_1\otimes f_1 +e_2\otimes f_2. \end{aligned}
Solution

Their coefficient matrices are

C1=(1122),C2=(1001).C_1= \begin{pmatrix} 1&1\\ 2&2 \end{pmatrix}, \qquad C_2= \begin{pmatrix} 1&0\\ 0&1 \end{pmatrix}.

The first matrix has rank one and factors as

C1=(12)(11).C_1= \begin{pmatrix} 1\\2 \end{pmatrix} \begin{pmatrix} 1&1 \end{pmatrix}.

Hence

T1=(e1+2e2)⊗(f1+f2).T_1=(e_1+2e_2)\otimes(f_1+f_2).

The second matrix has rank two, so T2T_2 is not simple.

  1. Prove that both local-factor products equal A⊗BA\otimes B:

    (A⊗IW)(IV⊗B)=A⊗B,(IV⊗B)(A⊗IW)=A⊗B.\begin{aligned} (A\otimes I_W)(I_V\otimes B) &=A\otimes B,\\ (I_V\otimes B)(A\otimes I_W) &=A\otimes B. \end{aligned}

    by checking simple tensors. Why is that sufficient?

Solution

For every v∈Vv\in V and w∈Ww\in W, set x=v⊗wx=v\otimes w and abbreviate

L=(A⊗IW)(IV⊗B),R=(IV⊗B)(A⊗IW).\begin{aligned} L&=(A\otimes I_W)(I_V\otimes B),\\ R&=(I_V\otimes B)(A\otimes I_W). \end{aligned}

Then

Lx=(A⊗IW)(v⊗Bw)=Av⊗Bw.\begin{aligned} Lx &=(A\otimes I_W)(v\otimes Bw)\\ &=Av\otimes Bw. \end{aligned}

In the reverse order,

Rx=(IV⊗B)(Av⊗w)=Av⊗Bw.\begin{aligned} Rx &=(I_V\otimes B)(Av\otimes w)\\ &=Av\otimes Bw. \end{aligned}

The two linear maps agree on every simple tensor. Simple tensors span V⊗WV\otimes W, so the maps agree everywhere.

  1. Let AA have eigenvalues 22 and −1-1, and let BB have eigenvalues 33, 44, and 55, with eigenbases in both spaces. List the eigenvalues of A⊗BA\otimes B, compute its trace, and verify the trace-product identity.
Solution

The product eigenvectors have eigenvalues

6, 8, 10, −3, −4, −5.6,\ 8,\ 10,\ -3,\ -4,\ -5.

Their sum is

tr⁡(A⊗B)=24−12=12.\begin{aligned} \operatorname{tr}(A\otimes B) &=24-12\\ &=12. \end{aligned}

Separately,

tr⁡A=1,tr⁡B=12,\operatorname{tr}A=1, \qquad \operatorname{tr}B=12,

so

tr⁡(A⊗B)=(tr⁡A)(tr⁡B)=12.\operatorname{tr}(A\otimes B) =(\operatorname{tr}A)(\operatorname{tr}B) =12.
  • S. Axler, Linear Algebra Done Right, 3rd ed., Springer, 2015.
  • S. Roman, Advanced Linear Algebra, 3rd ed., Springer, 2008.
  • W. H. Greub, Multilinear Algebra, 2nd ed., Springer, 1978.
  • R. A. Horn and C. R. Johnson, Topics in Matrix Analysis, Cambridge University Press, 1991.
  • M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, Cambridge University Press, 2010.