Second Quantization

Second quantization is a mathematical framework for describing systems with variable or large numbers of identical quantum particles. Rather than tracking which particle is in which state (as in first quantization), second quantization describes the state of the system by specifying how many particles occupy each single-particle state. This occupation number representation is natural for both bosons and fermions and is the standard language of quantum field theory, condensed matter physics, and quantum chemistry.

First vs. Second Quantization

In first quantization, you label particles (electron 1, electron 2, …) and write wavefunctions that depend on all particle coordinates. For N electrons, the wavefunction psi(r_1, r_2, …, r_N) must be antisymmetric under exchange of any two electrons (Pauli exclusion principle), which requires Slater determinants and becomes unwieldy for large N.

In second quantization, you describe the same system by specifying the occupation of each spatial orbital: how many electrons (0 or 1 for fermions) are in orbital 1, orbital 2, and so on. The state |n_1, n_2, …, n_M> for M orbitals is called a Fock state. Antisymmetry is enforced automatically by the algebra of the creation and annihilation operators, not by the structure of the wavefunction.

Creation and Annihilation Operators

The fermionic creation operator a_i^dagger adds one electron to orbital i, and the annihilation operator a_i removes one electron from orbital i. They satisfy the canonical anticommutation relations:

{a_i, a_j^dagger} = delta_ij    (anticommutator)
{a_i, a_j} = 0
{a_i^dagger, a_j^dagger} = 0

The anticommutation relations enforce Pauli exclusion: applying a_i^dagger twice to a state gives zero, since you cannot put two fermions in the same orbital.

Molecular Hamiltonians in Second Quantization

The electronic Hamiltonian of a molecule in second quantization has the form:

H = sum_pq h_pq a_p^dagger a_q  +  (1/2) sum_pqrs g_pqrs a_p^dagger a_q^dagger a_r a_s

The first term (one-body) captures kinetic energy and electron-nucleus attraction. The second term (two-body) captures electron-electron repulsion. The coefficients h_pq and g_pqrs are molecular integrals that can be computed classically from the nuclear geometry.

This form is exact within a chosen set of orbitals. Computing the ground state of this Hamiltonian on a quantum computer requires mapping the fermionic operators to qubit operators, a step that is not trivial because qubits satisfy different algebra from fermions.

Mapping Fermions to Qubits

Quantum computers operate on qubits (spin-1/2 systems), not fermions. To simulate the molecular Hamiltonian on a quantum computer, the fermionic creation and annihilation operators must be mapped to qubit Pauli operators. The Jordan-Wigner transformation is the most widely used mapping:

a_j = (Z_0 x Z_1 x ... x Z_{j-1}) x (X_j + iY_j) / 2

The string of Z operators encodes the fermionic anticommutation relations using the ordering of qubits. After this mapping, the molecular Hamiltonian becomes a weighted sum of Pauli strings, which can be evaluated on a quantum computer using VQE or phase estimation.

Bosonic Second Quantization

For bosons, the creation and annihilation operators satisfy commutation (rather than anticommutation) relations:

[b_i, b_j^dagger] = delta_ij    (commutator)

Bosons can occupy the same mode without restriction, so occupation numbers can be any non-negative integer. Bosonic second quantization applies to photonic quantum computing, quantum optics simulations, and vibrational modes of molecules.

Relevance to Quantum Computing

Second quantization is the standard language for expressing molecular and materials Hamiltonians that quantum computers are designed to simulate. Libraries like OpenFermion, PySCF, and the quantum chemistry modules of PennyLane and Qiskit Nature all work in the second-quantized representation, automating the conversion from molecular geometry to qubit Hamiltonian.