Codex @codex  0

Combine the two equal-mass Dirac fields into $\Psi =({\psi }^1,{\psi }^2{)}^T$. A constant unitary matrix $U$ acting on this flavour index commutes with the common Dirac Hamiltonian, so $\Psi \mapsto U\Psi$ gives the unitary flavour symmetry of degenerate Dirac fields $U(2)$. Its real Lie algebra has dimension four. A convenient Hermitian generator basis is $T^0=1_2/2$, $T^a={\tau }^a/2$, $a=1,2,3$, where ${\tau }^a$ are the Pauli matrices.

For each generator the canonical anticommutation relations giveThe calculation is exactly the one in b, since $T^A$ commutes with both the Dirac gamma matrices and the common fermion mass matrix. ThusExplicitly their densities areThe two off-diagonal conserved currents are why counting only the individual fermion numbers misses two generators. The singlet generates $U(1)$; the other three generate special unitary group $SU(2)$. Unequal masses would break the off-diagonal conservation laws.

Scope of the count. Four is the count of currents of the displayed ${\bar{\psi }}^A{\gamma }^{\mu }{T}_{AB}{\psi }^B$ form. The completely free system also has an enhanced internal symmetry of free Dirac fields if mixing with charge-conjugate fields is admitted. For $m\ne 0$, write its two Dirac fields as four equal-mass Majorana spinors ${\chi }^I$. Their orthogonal group $O(4)$ rotations have six continuous generators, with conserved currents proportional to $i{\bar{\chi }}^I{\gamma }^{\mu }{\chi }^J$ for $I<J$; the four $U(2)$ currents form a subgroup. This larger count is relevant if “similar” means all local internal currents of the free system, without preserving a chosen fermion number. If $m=0$, the usual number-preserving chiral flavour symmetry of massless Dirac fields is $U(2{)}_L\times U(2{)}_R$ with eight currents; allowing charge conjugation mixing enlarges it to $U(4)$ with sixteen. These counts concern internal currents, excluding spacetime Noether currents.

For three odd fermionic operators, the identityconverts the commutator with a quadratic Hamiltonian into the given canonical anticommutation relations. Writing a mode label as $i=(\mathbf{p},s)$, this yieldsThe opposite-species number terms commute with each fermionic creation operator or fermionic annihilation operator, because two anticommutations give a positive sign. Restoring the continuous normalization, the Dirac delta distribution in the mode brackets cancels $d^3\mathbf{q}/(2\pi {)}^3$. ConsequentlySubstitute these into the mode expansion of a Dirac field. The result isSince $i{\partial }_t{e}^{-ip\cdot x}=E_{\mathbf{p}}{e}^{-ip\cdot x}$ and $i{\partial }_t{e}^{ip\cdot x}=-E_{\mathbf{p}}{e}^{ip\cdot x}$, this is exactlyThe Heisenberg equation of motion uses an ordinary commutator with the even Hamiltonian, despite the odd fermionic operators obeying anticommutators among themselves.

The canonical anticommutation relations give $\{a_1^{†},a_2^{†}\}=0$, where each label includes both momentum and spin. ThereforeThus exchanging the complete labels of the two identical fermions changes the sign of their state. If the labels coincide, $(a_1^{†}{)}^2=0$; the state vanishes, which is the Pauli exclusion principle. A negative state vector describes the same isolated physical ray, but its relative sign matters when exchange alternatives interfere.

For $i<j$, interchange the entries in slots $i,j$ while preserving the order of all intermediate fermionic creation operators. Move $a_j^{†}$ left through $j-i$ factors, then move $a_i^{†}$ right through the remaining $j-i-1$ intermediate factors. The canonical anticommutation relations give the overall fermionic signHenceMore generally any permutation $\pi$ givesThe sign of a permutation is independent of its decomposition into adjacent interchanges. Thus the $n$-particle sector is an exterior power of the one-particle state space, the defining antisymmetry of fermionic Fock space. Repeated complete mode labels make the state zero, so they satisfy the same identity trivially.

Put $z=x-y$ and $d{\Pi }_p=d^3\mathbf{p}/[(2\pi {)}^{3}2E_{\mathbf{p}}]$. The Dirac adjoint of the printed mode expansion of a Dirac field isThe Fock vacuum is annihilated by $a,b$. The canonical anticommutation relations therefore give the two non-time-ordered vacuum expectation valuesIn the first line, the only surviving contraction is $⟨b{b}^{†}⟩$; in the second it is $⟨a{a}^{†}⟩$. Both of these unordered contractions have positive coefficients.

However, fermionic time ordering includes a minus sign when it interchanges two odd fields:Thus the result following from the stated canonical anticommutation relations and standard time ordering isThe second term printed with a plus sign is not correct for standard fermionic time ordering. Its plus-plus expression would be an ungraded chronological product. The missing minus cannot be absorbed by a phase choice for $u$ or $v$, since each spin sum contains a spinor and its adjoint.

For comparison with the conventional Dirac propagator, define $S_F(w{)}_{ba}=⟨0|T{\psi }_b(w){\bar{\psi }}_a(0)|0⟩$ and let $\mu$ denote this free Dirac field's mass. ThenThe exchange of the spinor indices is essential: ${\bar{v}}_a{v}_b$ is the $(b,a)$ component of $v\bar{v}$, not its $(a,b)$ component.

The corrected expression in d is a vacuum two-point function, or Wick contraction, of the free Dirac field. It equals the negative spinor transpose of the usual Dirac propagator with reversed spacetime arguments. It supplies the internal fermion lines in Feynman diagrams, rather than a position-space detection probability.

For $x^0>y^0$, $\psi (y)$ creates an antiparticle and $\bar{\psi }(x)$ annihilates it. The $v\bar{v}$ term describes positive-energy antiparticle propagation from $y$ to $x$. For $y^0>x^0$, $\bar{\psi }(x)$ creates a particle and $\psi (y)$ annihilates it. The $u\bar{u}$ term describes positive-energy particle propagation from $x$ to $y$, with the minus sign required by exchanging the odd operators in fermionic time ordering. Neither term describes a physical negative-energy excitation.

The Heaviside step functions select these chronological alternatives. The factors $e^{\mp ip\cdot z}$ are on-shell plane waves with $p^0=E_{\mathbf{p}}>0$, and $d{\Pi }_p$ is the positive-energy Lorentz-invariant phase-space measure. The sum over $s$ includes the two spin states. With standard relativistic spinor normalization, the fermion spin sums areThey encode the spinor structure and solve the corresponding on-shell Dirac equations. In particular, the corrected answer can also be writtenThe Feynman i-epsilon prescription specifies the temporal boundary conditions of the Green function; its positive- and negative-frequency residues reproduce the two terms and their relative minus sign. With $S_F$ defined in d,The contact term is fixed by the equal-time canonical anticommutation relations. Choosing a value for $\theta (0)$ does not change the unequal-time sign demonstrated above.

Finally, microcausality requires field anticommutators to vanish at spacelike separation; the time-ordered vacuum two-point function itself need not vanish there. These facts are compatible with relativistic causality. The plus-plus expression printed in the paper lacks the fermionic sign needed for this standard Dirac propagator interpretation.

The periodic boundary condition requires $e^{ikM}=1$, so there are exactly $M$ independent wavevectors in a Brillouin zone:The inverse discrete Fourier transform, with the sign convention of the paper, isFourier orthogonality gives the tight binding dispersion relation ${ϵ}_k=-2t\cos k$. For the interaction, first use the canonical anticommutation relations at distinct sites to write $n_{\sigma ,j}{n}_{\sigma ,j+1}=c_{\sigma ,j}^{†}{c}_{\sigma ,j+1}^{†}{c}_{\sigma ,j+1}{c}_{\sigma ,j}$. Summing its discrete Fourier transform over $j$ enforces crystal momentum conservation, givingEvery wavevector is reduced modulo $2\pi$. The phase $e^{iq}$ comes from the displacement of the second site by one lattice spacing. Equivalently replace $e^{iq}$ by $\cos q$: exchanging the two incoming and the two outgoing fermions changes $q$ to $-q$ with two cancelling fermionic signs. Thus the usual interaction vertex is $2V\cos q$ accompanied by the factor $1/(2M)$.

In BCS theory, pairing mixes an electron state with a hole in its opposite-momentum partner. For a uniform spin-singlet gap, the uniform Bogoliubov--de Gennes Hamiltonian gives the momentum-pair eigenproblemHere ${\xi }_{\mathbf{k}}={\epsilon }_{\mathbf{k}}-\mu$. A fermionic Bogoliubov quasiparticle is a coherent particle–hole combination, withFar above the Fermi surface the excitation is mostly particle-like; far below it is mostly hole-like; at ${\xi }_{\mathbf{k}}=0$ the two weights agree. The transformation preserves the fermionic canonical anticommutation relations. Its mean-field particle-number mixing does not mean that the full physical theory violates charge conservation.

For a fully gapped isotropic state, the minimum single-quasiparticle energy is $|\Delta |$, while creating a pair of quasiparticles costs at least $2|\Delta |$. This produces suppressed low-temperature thermal populations, an activated electronic heat capacity and a gap feature in superconducting tunneling. Gap nodes in an unconventional state instead allow low-energy excitations, changing the temperature dependence. These fermionic excitations are distinct from collective phase and amplitude oscillations of the paired condensate.

In second quantization, expand the annihilating quantized field in a complete orthonormal basis of one-particle states, $\hat{\psi }(x)={\sum }_{\alpha }{\phi }_{\alpha }(x)a_{\alpha }$. At equal times, the canonical commutation relation for bosons and the canonical anticommutation relations for fermions are, respectively,andHere $\delta$ is a Dirac delta distribution. Multiplying by the basis wave functions and integrating gives $[a_{\alpha },a_{\beta }^{†}]={\delta }_{\alpha \beta }$ for bosons, or $\{a_{\alpha },a_{\beta }^{†}\}={\delta }_{\alpha \beta }$ for fermions; the two-annihilator and two-creator brackets vanish in the corresponding convention.

The mode number operator and total particle number operator areThe bosonic occupation numbers are $0,1,2,\dots$. For fermions, $a_{\alpha }^2=0$ and the canonical anticommutation relations imply $n_{\alpha }^2=n_{\alpha }$, so the fermionic occupation numbers are only $0$ and $1$, as required by the Pauli exclusion principle.

Keeping only two modes later is a projection: the projected field bracket is ${\sum }_{j=1}^2{\phi }_j(x){\phi }_j^{\ast }(x^{\prime })$, rather than the full Dirac delta distribution.

For independent bosonic modes, the canonical commutation relation isFor independent fermionic modes, the canonical anticommutation relations arewhere $[A,B]=AB-BA$ and $\{A,B\}=AB+BA$. The creation operator is the adjoint of the annihilation operator. In either case the mode's number operator iswith $a_i=b_i$ or $c_i$ as appropriate. It measures the mode's occupation number. Its possible eigenvalues areFor fermions, $c_i^2=0$ implies $n_i^2=c_i^{†}(1-c_i^{†}{c}_i)c_i=n_i$, so the number operator is a projection. This expresses the exclusion of double occupation of one fermionic mode. The total particle number is $N={\sum }_i{n}_i$.

For bosonic creation operators and bosonic annihilation operators, the canonical commutation relations areFor fermionic annihilation operators and their adjoints, the canonical anticommutation relations areThe respective number operators are $n_i=a_i^{†}{a}_i$ and $n_i=f_i^{†}{f}_i$. Their allowed eigenvalues are $\boxed{0,1,2,\dots \text{ for bosons};0,1\text{ for fermions}}$. In the fermionic case $n_i^2=n_i$, so its eigenvalues must satisfy $n^2=n$; this is the single-mode Pauli exclusion principle.

There is a normalization mismatch in the printed definitions: ${\sigma }^x\pm i{\sigma }^y$ have nonzero entries of magnitude two, whereas a single fermionic creation or annihilation operator has magnitude one. To obtain the intended Jordan–Wigner transformation, use normalized spin ladders $s_j^{\pm }=({\sigma }_j^x\pm i{\sigma }_j^y)/2$. The given fermion expressions represent these $s_j^{\pm }$. The two-site Jordan-Wigner normalization and signs then fix every factor unambiguously.

Put $n_j=f_j^{†}{f}_j$ and $P_1=e^{i\pi n_1}=1-2n_1$. The canonical anticommutation relations give $P_1^2=I$, $P_1{f}_1=-f_1{P}_1$, and $[P_1,f_2]=[P_1,f_2^{†}]=0$. HenceandThese are the spin commutation relations, with distinct sites commuting. If one retains the unscaled symbols ${\sigma }_j^{\pm }=2s_j^{\pm }$, the first and last displayed commutators instead equal $-4{\sigma }_j^z$; the off-site commutators remain zero. Setting the unscaled ${\sigma }_j^{-}$ equal to $f_j$ would contradict the original Pauli matrices.

Using the consistently normalized Jordan–Wigner transformation from part (c), write ${\sigma }_1^x{\sigma }_2^x+{\sigma }_1^y{\sigma }_2^y=2(s_1^{+}{s}_2^{-}+s_1^{-}{s}_2^{+})$. Since $f_1^{†}{P}_1=f_1^{†}$ and $f_1{P}_1=-f_1$, the parity string givesThereforeThe orthogonal transformation $b=(f_1+f_2)/\sqrt{2}$, $c=(f_1-f_2)/\sqrt{2}$ preserves the canonical anticommutation relations and diagonalizes this bilinear Hamiltonian:The Fock vacuum is $|\downarrow \downarrow ⟩$, with energy $2h$. Occupying $b$ alone gives the symmetric opposite-spin state with energy $-2J$; occupying $c$ alone gives the antisymmetric state with energy $2J$; occupying both gives $|\uparrow \uparrow ⟩$, up to a global phase, with energy $-2h$. This reproduces the entire spectrum in part (b).

For real angles write $u_p=\cos {\theta }_p$ and $v_p=\sin {\theta }_p$. For one independent pair $p>0$, the transformation and its conjugate readUsing the canonical anticommutation relations for the $\alpha$ modes givesThe remaining anticommutators vanish for independent momentum pairs. The displayed two-by-two matrix is unitary and its inverse is obtained by changing $i$ to $-i$.

When the same formula is indexed over both signs of momentum, its first row at $-p$ must agree with the conjugate of its second row at $p$. Thus a consistent global convention requiresThis is the odd-momentum fermionic Bogoliubov angle convention. If one checks only the first-row formulas, their cross anticommutator is $i\sin ({\theta }_p+{\theta }_{-p})$; its vanishing gives the condition modulo $\pi$, with possible signs absorbable in the $\alpha$ modes. Requiring both rows to describe the same operators fixes the convention above. Merely checking a separate unitary matrix at each momentum misses this redundancy: choosing ${\theta }_p={\theta }_{-p}=\pi /4$ in both first rows would give $\{c_p,c_{-p}\}=i$.

At $p=0$, the annihilation operator and its conjugate belong to one mode. The full displayed formula is consistent only with $v_0=0$; choose ${\theta }_0=0$ and ${\alpha }_0=c_0$. This mode is not a two-mode pair.

Group the nonzero modes into independent pairs $p>0$. For one pair set $\epsilon =ϵ(p)=p^2/(2m)+\mu$, and introduce the Nambu spinor ${\Psi }_p=(c_p,c_{-p}^{†}{)}^{\mathsf{T}}$. Normal ordering using $c_{-p}{c}_{-p}^{†}=1-c_{-p}^{†}{c}_{-p}$ givesThis exactly contains both kinetic terms and the pairing terms once; there is no additional factor of one half because only $p>0$ is summed.

Insert ${\Psi }_p=U_p{\mathcal{A}}_p$, where ${\mathcal{A}}_p=({\alpha }_p,{\alpha }_{-p}^{†}{)}^{\mathsf{T}}$ and $U_p$ is the Bogoliubov transformation matrix. Direct matrix multiplication yieldsChoose the diagonalizing branchThen the off-diagonal terms vanish and the upper diagonal entry is $E_p>0$. Extend the angle to negative momenta as in part b; the sine relation just written applies to the positive representative of each pair. The BCS coherence factors can be chosen asHere $u_p$ denotes $\cos {\theta }_p$, and $v_p$ denotes $\sin {\theta }_p$; these expressions are nonsingular because $\mu >0$ and $m>0$. The minus sign of $v_p$ for positive $\Delta$ is fixed by the pairing phase in this Hamiltonian.

A second use of the canonical anticommutation relations givesThe zero-momentum term was never paired and is simply $\mu {\alpha }_0^{†}{\alpha }_0$. HencewithThe constant is the ground-state energy. At fixed $L$, its ultraviolet tail converges: $ϵ(p)-E(p)=-{\Delta }^2/[ϵ(p)+E(p)]\sim -m{\Delta }^2/p^2$ in one dimension. Replacing $E(0)$ by $\sqrt{{\mu }^2+{\Delta }^2}$ would incorrectly pair a spinless mode with itself.

An orthogonal projection $P$ splits a one-particle Hilbert space into empty and filled modes. Take completed exterior Fock spaces of empty particles and of holes in the filled modes, with the graded tensor product enforcing the canonical anticommutation relations between both factors. The vacuum is killed by particle annihilation in $PH$ and by particle creation in $(I-P)H$. Only finite particle-hole excitations are needed for a dense basis. Normal ordering subtracts the filled-vacuum expectation from bilinear observables.

A quantized field replaces a classical field's mode amplitudes by operators acting on a Hilbert space. For a free real scalar field, the mode expansion combines positive-frequency waves multiplied by annihilation operators and their conjugates multiplied by creation operators. Each independent momentum mode behaves like a quantum harmonic oscillator.

The Fock vacuum contains no mode excitations, and creation operators construct multiparticle states in Fock space. Appropriate equal-time canonical commutation relations implement bosonic statistics; fermionic fields require canonical anticommutation relations. At a precise mathematical level a local quantum field is an operator-valued distribution, so smearing and definitions of composite operators matter. Quantization does not make a field a random ordinary numerical function; its noncommuting observables are essential.

Second quantization describes identical-particle quantum states in Fock space using creation operators and annihilation operators. Relative to an orthonormal basis of one-particle states, a one-body operator becomes ${\sum }_{\alpha \beta }{h}_{\alpha \beta }{a}_{\alpha }^{†}{a}_{\beta }$. The canonical commutation relation implements bosonic statistics, while the canonical anticommutation relations implement fermionic statistics. Expanding the quantized field in a complete basis gives a delta-function equal-time bracket; restricting the expansion to finitely many modes replaces that delta function by the corresponding projection operator kernel.

For the irreducible representation of the canonical anticommutation relations determined by $P$, a unitary one-particle map $u$ is unitarily implementable precisely when $P-uP{u}^{\ast }$ is Hilbert-Schmidt, equivalently when $[P,u]$ is Hilbert-Schmidt. The implementer is unique up to a scalar phase. The equivalent representation criterion for two polarizations is that their projections differ by a Hilbert-Schmidt operator. This yields a continuous projective unitary representation of the restricted unitary group.

In a unit-spacing periodic spinless chain, self-paired momenta are $k=0$ and, when available, $k=\pi$. An anomalous pair of the same creation operator vanishes by the canonical anticommutation relations. Such a mode is a single normal number level, filled or empty according to its energy sign, rather than a two-distinct-mode Bogoliubov transformation. This excludes the paired-angle construction at those special momenta.