Let $M_n$ be the "affine monoid scheme" of $n \times n$ matrices under multiplication (like an affine group scheme but no inverses).

**Claim:** Every polynomial monoid homomorphism $M_n \to M_1$ is a non-negative integer power of the determinant; in particular the determinant is the unique generator of the monoid of such homomorphisms.

The same is true for the more familiar affine group scheme $GL_n$ (with "integer power") but then we need to distinguish between the determinant and its inverse, which we can do using the fact that the determinant extends to $M_n$.

*Proof.* Let $\varphi : M_n \to M_1$ be such a homomorphism. Passing to any algebraically closed field of characteristic zero, the diagonalizable matrices are Zariski dense and $\varphi$ must be invariant under conjugation, so $\varphi$ is determined by its restriction to diagonal matrices and therefore must be a polynomial in the eigenvalues. Since diagonal matrices are conjugate to all their diagonal permutations, $\varphi$ must be a symmetric polynomial in the eigenvalues. By restricting to each diagonal entry separately $\varphi$ restricts to a polynomial homomorphism $M_1 \to M_1$ which must be the same for each diagonal entry (because each of the corresponding copies of $M_1$ is conjugate in $M_n$), and it's not hard to show with a direct computation that such a homomorphism (concretely, a polynomial $f(x) \in \mathbb{Z}[x]$ satisfying $f(xy) = f(x) f(y)$ and $f(1) = 1$) must be of the form $f(x) = x^k$ for some non-negative integer $k$. By writing a diagonal matrix $\text{diag}(\lambda_1, \dots \lambda_n)$ as a product of diagonal matrices of the form $\text{diag}(1, \dots \lambda_i, \dots 1)$ we conclude that $\varphi = (\lambda_1 \dots \lambda_n)^k$ for some non-negative integer $k$ as desired. (Unfortunately this argument doesn't construct the determinant; we need to show that it exists separately.) $\Box$

An alternative way to state the conclusion is that the bialgebra $\mathcal{O}_{M_n}$ of functions on $M_n$, which concretely is $\mathbb{Z}[x_{ij}], 1 \le i, j \le n$ equipped with the comultiplication $\Delta x_{ij} = \sum_k x_{ik} \otimes x_{kj}$ encoding matrix multiplication, has the property that its monoid of grouplike elements (elements satisfying $\Delta x = x \otimes x$) is generated by the determinant.

This is indirect, plus it doesn't immediately tell you that the determinant is invariant under change of basis (although that's easy to prove once you know it's a homomorphism), but it has the pleasant property of specifying a single universal polynomial over $\mathbb{Z}$ which specializes to the determinant over every field (even every commutative ring), without needing to work over any particular field.

Instead of talking about generators we can also single out the determinant using the fact that, among all the powers $\det^k$, the determinant itself is the only one which is multilinear.