# Definition

To derive the determinant rigorously from the concept of nn -dimensional volume, we define a function

f:(Rn)nRf: (\mathbb{R}^n)^n \to \mathbb{R}

that maps nn vectors (v1,v2,,vn)(\mathbf{v}_1, \mathbf{v}_2, \dots, \mathbf{v}_n) to a scalar representing the signed volume of the parallelotope they span. For this function to behave like a measure of volume in Euclidean space, it must satisfy three fundamental axioms:

  1. Multilinearity: For any ii, ff is linear with respect to the ii-th vector:

    f(,cvi+u,)=cf(,vi,)+f(,u,)f(\dots, c\mathbf{v}_i + \mathbf{u}, \dots) = cf(\dots, \mathbf{v}_i, \dots) + f(\dots, \mathbf{u}, \dots)

  2. Alternating Property: If any two vectors are equal, the volume is zero (the shape collapses):

    f(v1,,vi,,vi,,vn)=0f(\mathbf{v}_1, \dots, \mathbf{v}_i, \dots, \mathbf{v}_i, \dots, \mathbf{v}_n) = 0

  3. Normalization: The volume of the unit hypercube spanned by the standard basis {e1,,en}\{\mathbf{e}_1, \dots, \mathbf{e}_n\} is unity:

    f(e1,e2,,en)=1f(\mathbf{e}_1, \mathbf{e}_2, \dots, \mathbf{e}_n) = 1

From the alternating property, we can immediately derive that swapping any two vectors negates the result. Consider f(,vi+vj,,vi+vj,)=0f(\dots, \mathbf{v}_i + \mathbf{v}_j, \dots, \mathbf{v}_i + \mathbf{v}_j, \dots) = 0. By expanding this using multilinearity and canceling the terms f(,vi,,vi,)f(\dots, \mathbf{v}_i, \dots, \mathbf{v}_i, \dots) and f(,vj,,vj,)f(\dots, \mathbf{v}_j, \dots, \mathbf{v}_j, \dots) which are zero, we arrive at:

f(,vi,,vj,)=f(,vj,,vi,)f(\dots, \mathbf{v}_i, \dots, \mathbf{v}_j, \dots) = -f(\dots, \mathbf{v}_j, \dots, \mathbf{v}_i, \dots)

Now, let AA be an n×nn \times n matrix where the jj-th column is vj\mathbf{v}_j. We express each column vector in the standard basis: vj=i=1naijei\mathbf{v}_j = \sum_{i=1}^n a_{ij} \mathbf{e}_i. Substituting these into the volume function ff, we use multilinearity to expand the expression completely:

f(v1,,vn)=f(i1=1nai11ei1,,in=1nainnein)=i1,,in=1nai11ainnf(ei1,,ein)f(\mathbf{v}_1, \dots, \mathbf{v}_n) = f\left(\sum_{i_1=1}^n a_{i_1 1} \mathbf{e}_{i_1}, \dots, \sum_{i_n=1}^n a_{i_n n} \mathbf{e}_{i_n}\right) = \sum_{i_1, \dots, i_n = 1}^n a_{i_1 1} \dots a_{i_n n} f(\mathbf{e}_{i_1}, \dots, \mathbf{e}_{i_n})

Due to the alternating property, any term in this summation where the indices (i1,,in)(i_1, \dots, i_n) are not distinct vanishes. Therefore, we only sum over sequences that are permutations σ\sigma of the set {1,,n}\{1, \dots, n\}. For any such permutation, the alternating property implies:

f(eσ(1),,eσ(n))=sgn(σ)f(e1,,en)f(\mathbf{e}_{\sigma(1)}, \dots, \mathbf{e}_{\sigma(n)}) = \text{sgn}(\sigma) f(\mathbf{e}_1, \dots, \mathbf{e}_n)

Applying the normalization axiom f(e1,,en)=1f(\mathbf{e}_1, \dots, \mathbf{e}_n) = 1, the expression simplifies to the Leibniz formula, which serves as the formal definition of the determinant:

det(A)=σSnsgn(σ)j=1naσ(j),j\det(A) = \sum_{\sigma \in S_n} \text{sgn}(\sigma) \prod_{j=1}^n a_{\sigma(j), j}

This proof demonstrates that the determinant is the unique function that satisfies the geometric requirements of a signed volume scaling factor.