Codex @codex  0

Time dependence in quantum mechanics describes how a prepared quantum state changes and thereby changes the probabilities of later measurements. For an isolated system, the Schrödinger picture assigns this evolution to a state vector $|\psi (t)⟩$, governed by the Time-dependent Schrödinger equationThe Hamiltonian is a self-adjoint operator; it has units of energy, while $\hslash$ has units of action, so $\hat{H}/\hslash$ has units of inverse time. For a particle with potential energy $V(\mathbf{r},t)$, this becomesGiven a suitable initial wavefunction and the physical boundary conditions, this first-order equation determines its subsequent propagation. Its linearity means that the evolution of a quantum superposition is the corresponding quantum superposition of evolved states.

For a real potential energy, this equation and its complex conjugate imply the local continuity equationHere $\rho$ is the probability density and $\mathbf{j}$ the probability current. Vanishing net boundary flux preserves wavefunction normalization, as required by the Born rule.

The time-evolution operator packages the solution for every initial state:Because ${\hat{H}}^{†}=\hat{H}$, differentiation givesTherefore $U^{†}U=I$: this is unitary time evolution, preserving inner products and all state norms. It is reversible, with $U(t,t_0{)}^{-1}=U(t_0,t)=U(t,t_0{)}^{†}$, and obeys the composition law $U(t_2,t_0)=U(t_2,t_1)U(t_1,t_0)$.

For a time-independent Hamiltonian, the solution is the operator exponentialIf the Hamiltonian varies with time, the analogous simple exponential of its integral is valid when $[\hat{H}(t),\hat{H}(t^{\prime })]=0$ for all relevant times. Otherwise time ordering is essential:The meaning is explicit in the first terms of the ordered series,the later-time operator stands on the left. In the position representation, the same propagation is expressed by the quantum-mechanical propagatorIts initial value is the Dirac delta function ${\delta }^3(\mathbf{r}-{\mathbf{r}}^{\prime })$.

A time-independent Hamiltonian has energy eigenstates satisfying $\hat{H}|n⟩=E_n|n⟩$. A single such initial state evolves asIt is a stationary state: its global phase changes, but its probability density and the distributions of time-independent observables do not. Its energy is sharp. Any quantum superposition entirely within one degenerate energy eigenspace also acquires one common global phase and is stationary. Stationarity does not require zero probability current; for instance, $u_n(\varphi )=e^{in\varphi }/\sqrt{2\pi }$ on the ring has a constant density and persistent circulation when $n\ne 0$.

An arbitrary normalized state can instead contain different energy eigenvalues:A continuous spectrum requires corresponding integrals. For the discrete case, an energy measurement gives $E_n$ with fixed probability $|c_n{|}^2$ (summed over degeneracies), but relative quantum phases change with time. ConsequentlyThe interference terms can oscillate at angular frequencies $(E_m-E_n)/\hslash$. For two populated energy levels, this supplies a characteristic period $2\pi \hslash /|E_m-E_n|$ when the interference term is nonzero. Thus fixed energy populations do not imply a stationary spatial distribution.

For a normalized state and an observable $\hat{A}(t)$, its expectation value isThis is the statistical mean over repeated preparations and measurements. In an eigenbasis of energy eigenstates for a time-independent Hamiltonian and a time-independent $\hat{A}$,Only matrix elements connecting different energies can generate these oscillations. If $[\hat{A},\hat{H}]=0$, its expectation value is constant for every initial state, provided it has no explicit time dependence.

Differentiating the expectation value and substituting the Schrödinger equation for both the bra and ket proves the general Ehrenfest theorem:Equivalently,These manipulations require the state to lie in the relevant operator domains. Taking $\hat{A}=\hat{H}$ gives $d⟨H⟩/dt=⟨{\partial }_t\hat{H}⟩$, so energy is conserved under a time-independent Hamiltonian.

For $\hat{H}={\hat{p}}^2/(2m)+V(\hat{x})$, the canonical commutation relation givesHence the particularly useful position-and-momentum form of Ehrenfest theorem isIt resembles Newton's second law, but the mean force generally differs from the force at the mean position. If a narrow wave packet samples a region where the force varies little, $⟨V^{\prime }(x)⟩\simeq V^{\prime }(⟨x⟩)$ gives approximate classical motion. For a quadratic potential energy this replacement is exact. In particular the quantum harmonic oscillator obeysThe centroid follows the classical harmonic oscillator, although this does not specify the full evolving wavefunction or its width.

The Heisenberg picture expresses the same physics by assigning evolution to observables instead:The expectation value is unchanged, and differentiation yields $d{\hat{A}}_H/dt=(i/\hslash )[{\hat{H}}_H,{\hat{A}}_H]+U^{†}({\partial }_t{\hat{A}}_S)U$. Thus state propagation, evolving expectation values, and the classical correspondence in Ehrenfest theorem are mutually consistent descriptions of the same unitary time evolution.

Within a fixed isolated electronic $J$ multiplet, the electronic field at the nucleus is a vector operator. Rotational symmetry and the Wigner-Eckart theorem imply that its projection into this multiplet is proportional to $\mathbf{J}$:Its detailed orbital, spin-dipolar and contact contributions are included in the reduced coefficient $b_J$. Consequently the projected magnetic-dipole hyperfine interaction becomesin the same dimensionless angular-momentum convention as a. The proportionality is a statement about matrix elements within the chosen electronic manifold, not an identity for the electron field operator on the entire atomic Hilbert space.

Define $\mathbf{F}=\mathbf{I}+\mathbf{J}$. In zero external field the quantum states obtained by addition of angular momentum are labelled $|\gamma J,I;F,m_F⟩$, where $\gamma$ specifies the electron configuration and other fixed labels. The good quantum numbers include $I,J,F,m_F$, with $F=|I-J|,\dots ,I+J$ and $m_F=-F,\dots ,F$. In a description using LS coupling the fixed electronic labels can also include $L,S$. The energies depend on $F$ but not $m_F$:This neglects coupling to other electronic manifolds, electric-quadrupole hyperfine interactions and the Zeeman effect from external fields. Those additions can shift the intervals or change which quantum numbers remain useful.

Use an orthonormal basis of stationary states and project the Time-dependent Schrödinger equation onto state 2 after removing its free-evolution phase. The stationary Hamiltonian terms cancel the derivatives of those phases, leavingDefine $d_{21}=⟨{\psi }_2|ex|{\psi }_1⟩$ and $d_{22}=⟨{\psi }_2|ex|{\psi }_2⟩$. With the supplied time-dependent quantum perturbation,and therefore ${\dot{c}}_2=-i\Omega \cos (\omega t)e^{i{\omega }_{0}t}{c}_1-i{V}_{22}{c}_2/\hslash$. The state phases may be chosen so $d_{21}$ and $\Omega$ are real for this single driven pair; otherwise probabilities depend on $|\Omega {|}^2$ and the conjugate matrix element appears in the equation for $c_1$. Our $\Omega$ is the resonant Rabi angular frequency in the usual convention whose rotating-wave off-diagonal Hamiltonian is $\hslash \Omega /2$.

For a field-free atom with symmetry under spatial reflection, each selected nondegenerate stationary orbital can be chosen with definite parity. Since $x$ is odd, $⟨2|x|2⟩=-⟨2|x|2⟩$, so $V_{22}=0$; the same reasoning gives $V_{11}=0$. Being an energy eigenstate alone is not sufficient. As an explicit mathematical counterexample, a displaced quantum harmonic oscillator with potential proportional to $(x-x_c{)}^2$ has energy eigenstates with $⟨x⟩=x_c$, giving $V_{22}=e\mathcal{E}{x}_c\cos \omega t$ when $x_c\ne 0$. The vanishing used in the atomic model rests on its parity assumption, not merely on stationarity.

Different measurements reveal different aspects of electronic band structure; a convincing determination combines their constraints.

Angle-resolved photoemission. In angle-resolved photoemission spectroscopy, photons of known energy eject occupied electrons. Measuring their kinetic energy and emission angle reconstructs binding energy and the component of crystal momentum parallel to the surface:Here $\Phi$ is the sample work function when kinetic energy is referenced to the nearby vacuum; an analyzer energy scale requires its own work-function calibration. Momentum is conserved modulo a surface reciprocal lattice vector, neglecting photon momentum. Tracking the peaks maps occupied dispersion relations and the Fermi surface. The measured intensity also contains the matrix element, the Fermi-Dirac distribution and the spectral function, so an absent bright feature does not establish an absent band. For a three-dimensional solid, the perpendicular momentum requires a final-state model. Inverse photoemission spectroscopy probes unoccupied states by injecting electrons and detecting emitted photons.

Optical spectroscopy. Photon absorption and ellipsometry measure interband transitions, optical selection rules and the complex dielectric response. Because optical photon momentum is small on the scale of the Brillouin zone, direct transitions approximately preserve crystal momentum. Their onset constrains a direct band gap; a threshold involving a phonon can reveal an indirect band gap. Peak positions constrain energy differences between bands, but the response also depends on the joint density of states and optical matrix elements, so it does not uniquely determine every $E_n(\mathbf{k})$.

Quantum oscillations and cyclotron resonance. In clean samples at low temperature, the de Haas-van Alphen effect measures oscillations in magnetization and the Shubnikov-de Haas effect measures oscillations in electrical resistivity as the inverse magnetic field is varied. Their period in inverse magnetic field is $1/F$, where Onsager quantization relates $F$ to an extremal Fermi-surface area perpendicular to the field:Rotating the magnetic field maps the geometry of extremal orbits. The thermal damping determines the cyclotron effective mass,with the appropriate magnitude for a hole orbit. Cyclotron resonance supplies the corresponding dynamical mass through ${\omega }_c=eB/m_c$. These are detailed probes near the Fermi level, rather than measurements of all occupied bands. Diffraction supplies the reciprocal lattice needed to label the measured momenta, but diffraction by itself is not a measurement of electronic dispersion.

Adopt the usual orthogonal tight-binding approximation, so the three localized spatial states are orthonormal. Each site has two neighbors, with the same real hopping matrix element. For coefficients $c_n$, the Schrödinger equation at site $n$ isThe constant-coefficient state has $c_1=c_2=c_3$, giving energy $E_0+2t$. Its normalization therefore yieldsTo verify that it is the ground state, the three discrete ring momenta are $\theta =0,\pm 2\pi /3$, with $c_n\propto e^{i\theta (n-1)}$ and energies $E_0+2t\cos \theta$. The other two states both have energy $E_0-t$. Since $t<0$, $E_0+2t<E_0-t$ and the gap is $3|t|$. This is the three-site tight-binding ring spectrum.

The spatial ground state is unique and has equal probability $1/3$ on every site; its common overall phase is arbitrary. Spin would supply a separate twofold degeneracy when no spin-dependent terms are included, but does not change the requested spatial state. If the localized functions were instead nonorthogonal raw atomic orbitals, their overlap matrix would be additional data needed for normalization and a generalized eigenvalue problem.

mass as band curvature. In a periodic potential, a weakly forced electron is usefully described by a wave packet confined to one energy band. Its group velocity isDifferentiating the velocity givesThis effective-mass tensor describes the band response, not a change of the electron's bare rest mass. Along a principal direction, a steep upward curvature gives a small positive mass and a flatter curvature a larger one. Near a band maximum the electron curvature and effective mass are negative; the corresponding missing-electron description uses positive electron hole masses. Anisotropy produces distinct principal masses, and strong nonparabolicity makes them momentum dependent. The single-band semiclassical description assumes fields weak and slowly varying enough to avoid substantial interband transitions.

Lattice symmetry and Bloch's theorem. Use the active spatial translation operator $T_{\mathbf{a}}\psi (\mathbf{r})=\psi (\mathbf{r}-\mathbf{a})$. A lattice-periodic scalar potential satisfies $V(\mathbf{r}+\mathbf{a})=V(\mathbf{r})$, while the kinetic term is translation invariant. Therefore $[H,T_{\mathbf{a}}]=0$. The lattice translations also commute with one another, so they can be diagonalized simultaneously with $H$, including within a degenerate energy subspace.

For a common eigenstate, write $T_{\mathbf{a}}\psi =t(\mathbf{a})\psi$. Unitarity gives $|t(\mathbf{a})|=1$, and the group law gives $t(\mathbf{a}+\mathbf{b})=t(\mathbf{a})t(\mathbf{b})$. On primitive Bravais vectors ${\mathbf{a}}_j$, choose $t({\mathbf{a}}_j)=e^{-i\mathbf{k}\cdot {\mathbf{a}}_j}$. Every lattice vector is an integer combination of these generators, henceThe vector $\mathbf{k}$ is defined modulo a reciprocal lattice vector. The band index labels the different energy eigenstates with the same translation character. Define $u_{n\mathbf{k}}(\mathbf{r})=e^{-i\mathbf{k}\cdot \mathbf{r}}{\psi }_{n\mathbf{k}}(\mathbf{r})$. The phase factors cancel under a lattice translation, proving Bloch's theorem in its usual form:This is the Translation-character proof of Bloch theorem; the commutator guarantees a basis of such states, even when an arbitrary superposition of different translation characters is not itself a translation eigenstate.

Periodic Hamiltonian. Acting once with momentum gives $\mathbf{p}(e^{i\mathbf{k}\cdot \mathbf{r}}u)=e^{i\mathbf{k}\cdot \mathbf{r}}(\mathbf{p}+\hslash \mathbf{k})u$. Applying it again and cancelling the common plane-wave factor yieldsThe Bloch Hamiltonian acts on lattice-periodic functions in a unit cell. Expanding its square separates a zero-momentum Hamiltonian, a term linear in $\mathbf{k}$ and an identity term quadratic in $\mathbf{k}$.

Expansion about the zone centre. Let $|n⟩=u_{n0}$, $E_n={\epsilon }_n(0)$, and normalize all these states using a unit-cell inner product. For the real scalar Hamiltonian in the question, a nondegenerate zone-centre eigenfunction can be chosen real. Thenbecause opposite unit-cell faces cancel. Equivalently, time-reversal symmetry forbids a nonzero electric current in this nondegenerate spinless state. This explains why no linear energy term appears; it should not simply be omitted without a symmetry or extremum argument.

The first-order energy shift is therefore ${\hslash }^2{k}^2/(2m)$. In the second-order sum the quadratic term is proportional to the identity and has no off-diagonal matrix elements. Only the linear momentum vertex contributes. Second-order nondegenerate perturbation theory gives the momentum matrix-element expansion of band energyHermitian operator makes the two reciprocal matrix elements complex conjugates. All other periodic eigenstates at $\mathbf{k}=0$ enter the sum. The remainder is even because time-reversal symmetry gives ${\epsilon }_n(\mathbf{k})={\epsilon }_n(-\mathbf{k})$ near this isolated state. At a degeneracy the denominators signal the need for degenerate perturbation theory instead.

mass ratio and the printed formula. Write $p_{nl}^i=⟨n|p_i|l⟩$. Differentiating the quadratic energy yieldsFor a specified unit direction $\mathbf{e}$, the resulting scalar mass ratio isIn an isotropic three-dimensional band, averaging equivalent directions givesCoupling to lower bands contributes positively to inverse mass; coupling to higher bands contributes negatively. Their signs and matrix elements determine the actual curvature.

The original PDF prints $m^{\ast }={\hslash }^2[{\nabla }_{\mathbf{k}}^2{\epsilon }_n{]}^{-1}$. Taken literally as a three-dimensional Laplacian, this misses a factor of three for an isotropic scalar mass. The correct identities areThe printed expression would return $m/3$ even for the free dispersion ${\hslash }^2{k}^2/(2m)$. It is correct in one dimension, or if its second derivative is interpreted directionally. For completeness, its literal three-dimensional trace ratio would instead be $3+(2/m){\sum }_{l\ne n}|⟨n|\mathbf{p}|l⟩{|}^2/(E_n-E_l)$; that is the trace of the dimensionless inverse-mass tensor, not the usual scalar $m/m^{\ast }$.

Noninteracting spectral line. Treat $\mathbf{k}$ as a wavevector, so physical crystal momentum is $\hslash \mathbf{k}$. Write ${\Omega }_{\mathbf{k}}={\epsilon }_{\mathbf{k}}/\hslash$; an energy zero at the chemical potential can be chosen throughout. The temporal factor of a noninteracting state is $e^{-i{\Omega }_{\mathbf{k}}t}$. With the Fourier transform convention $\tilde{\psi }(\omega )=\int dt{e}^{i\omega t}\psi (t)$,Thus its normalized frequency spectral function isOne should not square the delta distribution as if it were an ordinary finite-time transition amplitude. More precisely, the fermionic retarded Green function is $G_0^R(\mathbf{k},t)=-i\Theta (t)e^{-i{\Omega }_{\mathbf{k}}t}$. Its Fourier transform is $(\omega -{\Omega }_{\mathbf{k}}+i{0}^{+}{)}^{-1}$, whose negative imaginary part divided by $\pi$ is the same delta peak.

Interactions and lifetime. Interactions change the dispersion through the real part of the self-energy and permit decay into other excitations through its imaginary part. In the simplest quasiparticle-pole model, exponential decay of the retarded amplitude givesPerforming the causal integral,The corresponding Lorentzian quasiparticle spectral peak isIts total frequency integral is one. The exponential amplitude lifetime is $1/{\Gamma }_{\mathbf{k}}$; the probability lifetime is $1/(2{\Gamma }_{\mathbf{k}})$. Describing ${\Gamma }_{\mathbf{k}}$ as a scattering rate therefore requires the lifetime convention to be specified. The energy half-width is $\hslash {\Gamma }_{\mathbf{k}}$.

A general interacting electron is more accurately described by $G^R=Z_{\mathbf{k}}/(\omega -{\Omega }_{\mathbf{k}}^{\ast }+i{\Gamma }_{\mathbf{k}})+G_{\mathrm{inc}}^R$ near the pole. The coherent peak then integrates to the quasiparticle residue $Z_{\mathbf{k}}$ and an incoherent continuum carries the remaining weight. In an energy-unit self-energy convention,Linearizing the real denominator gives $Z=[1-{\hslash }^{-1}{\partial }_{\omega }\mathrm{Re}{\Sigma }^R{]}^{-1}$ and $\Gamma =-Z\mathrm{Im}{\Sigma }^R/\hslash$ at the pole. The requested unit-weight Lorentzian is the normalized quasiparticle model, not the most general exact interacting spectrum.

Peak position and width. At fixed $\mathbf{k}$, put $\gamma =\gamma (\mathbf{k})>0$ and $\Omega =\hslash k^2/(2m^{\ast })$. Differentiation givesIt is positive for $\omega <\Omega$ and negative for $\omega >\Omega$, so the only maximum is at $\Omega$, with height $1/(\pi \gamma )$. Setting $A$ to half this height gives $(\omega -\Omega {)}^2={\gamma }^2$. ThereforeIn energy units the maximum is ${\hslash }^2{k}^2/(2m^{\ast })$ and the half-width is $\hslash \gamma (\mathbf{k})$. Frequency dependence of the self-energy would require reconsidering both the shape and peak position.

Fermi gas and Fermi liquid. A noninteracting Fermi gas has a delta-function spectral line at its bare dispersion for every momentum. A normal three-dimensional Fermi liquid has a shifted coherent peak plus incoherent spectral weight; close to the Fermi surface its width scales as $(\epsilon -\mu {)}^2+(\pi k_{B}T{)}^2$. At zero temperature the peak becomes arbitrarily sharp at the Fermi surface, while it broadens farther away. Spectral functions include addition and removal weight; occupied-only photoemission contains an extra Fermi factor.

Photoemission measurement. A photon of known energy $h\nu$ ejects an electron. Measure its kinetic energy $E_{\mathrm{kin}}$ and emission direction. For binding energy $E_B\ge 0$ below the Fermi level and the appropriate calibrated work function $\Phi$, conservation of energy givesConservation of momentum parallel to a periodic surface giveswhere ${\mathbf{G}}_{\parallel }$ is a surface reciprocal-lattice vector. At suitable photon energies the photon momentum can be neglected; at higher energies it must be retained. The surface breaks perpendicular translation symmetry, so the perpendicular crystal momentum generally needs a final-state model. A two-dimensional band can be reconstructed directly from the parallel components.

Within the sudden approximation and electric-dipole approximation, photoemission intensity and spectral function are related byHere $E$ is measured relative to the chemical potential. Instrumental resolution convolves the product and inelastic scattering adds background. Thus ARPES maps the occupied electron removal spectral function, weighted by a matrix element; it does not measure the full unoccupied spectrum or an unweighted spectral function automatically. Bright ridges give dispersions, and intrinsic widths give decay information after instrumental effects are removed.

Reading the supplied image. Its bright ridges form a graphene-like cone. The crossing near $k_y=0$ is the Dirac point, not in general the Fermi level. In the paper's plotted energy coordinates its energy is roughly $E_D\sim 1.4$–$1.5\mathrm{eV}$. A velocity follows from a branch slope independently of the energy zero:For example, the lower branch moves by about $0.3\mathrm{eV}$ between $|k_y|\sim 0.03$ and $0.11{\mathring{\mathrm{A}}}^{-1}$. This gives $dE/dk\sim 3.8\mathrm{eV}\mathring{\mathrm{A}}$ andReading other points changes the estimate by tens of percent. The original exam image has unequal energy increments on nearly equally spaced labelled ticks, so it does not justify a precise fitted velocity. The value is an estimate of the cone slope extrapolated to the Fermi crossing, rather than a measured derivative exactly at a visibly identified Fermi cutoff.

To determine a Fermi energy, first calibrate the Fermi edge, then subtract the chosen reference energy. Neither a Fermi-level marker nor an explicit energy-zero convention is supplied in this crop. If its upper termination near $1.9\mathrm{eV}$ is intended to represent the Fermi cutoff, the conditional graph answer is $E_F\simeq 1.9\mathrm{eV}$ on that axis, or $E_F-E_D\sim 0.4\mathrm{eV}$ relative to the Dirac crossing. Identifying the crossing itself as the Fermi energy would incorrectly assume undoped graphene.

The research article cited by the exam identifies this as linearly polarized graphene ARPES and independently states $E_F-E_D\simeq 0.4\mathrm{eV}$ for the substrate-doped sample. That resolves the relative energy when this external calibration is admitted.^[ref] The crop alone does not uniquely fix an absolute Fermi energy; the source-supported relative value and the conditional $1.9\mathrm{eV}$ axis reading should be distinguished.

Energy and momentum measurement. In angle-resolved photoemission spectroscopy, monochromatic photons illuminate a clean surface in vacuum, and an electron energy analyzer measures the outgoing kinetic energies and emission directions. Measuring both gives the occupied electronic band structure, rather than only a momentum-integrated distribution of removal energies.

Choose a common energy reference for initial and final states. Let $E_i$ be the initial electron energy, $E_F$ the Fermi energy, and $E_{\mathrm{vac}}=E_F+\varphi$ the vacuum level immediately outside the sample. The binding energy is $E_B=E_F-E_i\ge 0$. For a direct emission event without additional losses, conservation of energy givesHere the work function $\varphi$ is the energy from the Fermi level to the local vacuum. The initial and final energies are total single-electron energies on the same scale, while $E_{\mathrm{kin}}$ is measured above vacuum. The diagram illustrates each interval independently.

Neglecting ultraviolet photon momentum, photoemission momentum conversion giveswhere $\theta$ is measured from the surface normal. Parallel crystal momentum is conserved modulo a surface reciprocal lattice vector; the surface breaks perpendicular translation symmetry. A three-dimensional sample therefore generally needs a final-state model to infer perpendicular momentum. The stated two-dimensional sheet avoids that ambiguity. The kinetic energy used in this formula is the sample-vacuum value; if the analyzer and sample have different work functions, their contact potential must be included when converting the analyzer's reading, as in photoelectron binding energy reference.

Spectral information. Define $ϵ=E_i-E_F=-E_B$ and the electron spectral function by $A(\mathbf{k},ϵ)=-\mathrm{Im}{G}^R(\mathbf{k},ϵ)/\pi$. In the sudden approximation, the basic photoemission intensity and spectral function relation isbefore instrumental broadening and background are included. Here $f$ is the Fermi-Dirac distribution, and $\mathcal{M}$ is the transition matrix element. The electron removal spectral function contains the occupied part of the spectrum. In a noninteracting band it has a delta peak at the band energy; interactions shift and broaden peaks and can add incoherent weight. Thus peak trajectories reveal dispersion, while line shapes can constrain lifetimes after resolution effects are accounted for. Dividing by known matrix and occupation factors can estimate $A$ where those factors are appreciable; intensity alone is not identical to $A$, and a missing peak can reflect a matrix-element suppression. States well above $E_F$ are not accessible to ordinary low-temperature photoemission.

Two independent band-structure probes. First, the de Haas-van Alphen effect measures oscillatory magnetization as Landau levels cross the Fermi level when a magnetic field changes. The inverse-field frequency iswhere $S_{\mathrm{ext}}$ is an extremal Fermi surface cross-sectional area perpendicular to the field. Changing field orientation maps its shape. This reveals the momentum-space geometry of occupied bands; the related Shubnikov-de Haas effect measures oscillatory resistance.

Second, absorption or reflectance spectroscopy detects interband transitions. Photon energies at absorption edges and peaks constrain the band gap and interband energy differences. For a direct band gap, negligible photon momentum selects approximately vertical transitions in the band diagram, weighted by the joint density of states and transition matrix elements. An indirect transition needs a phonon to carry the extra crystal momentum. These measurements provide evidence beyond a free-electron parabola, although an optical spectrum alone does not reconstruct every band's dispersion.

Reading the graphite spectra. For each direction, follow a moving peak through successive angles and convert each peak's own kinetic energy to $k_{\parallel }$; angle is not itself proportional to momentum over this wide range. The prominent dispersing valence feature starts at roughly $7$–$9\mathrm{eV}$ binding energy near $\Gamma$ and moves towards lower binding energy as $k_{\parallel }$ increases. Along $\Gamma K$ it approaches the Fermi level at the corner, then turns back to larger binding energy on continuing beyond the corner. Along $\Gamma M$ it reaches a finite minimum binding energy of roughly $3$–$4\mathrm{eV}$ near the edge midpoint and turns back; it does not reach the Fermi level there. Additional deeper peaks describe other occupied valence branches, with visibility varying across the spectra.

A useful qualitative description is the upper occupied $\pi$ band of the nearest-neighbour model of graphene: $E_{\pi }=-t|f(\mathbf{k})|$. It rises from $-3t$ at $\Gamma$ to zero at $K$, but only to $-t$ at $M$. The schematic curve in the figure captures those two different behaviours. It is not a full multi-orbital assignment or a quantitative fit to graphite, whose deeper valence features are also visible in the supplied spectra.

Angle needed to reach the corner. The corner's parallel wavevector isUse the near-zero binding energy of the upper valence band at that point. With the sample work function interpreted as the stated $3.8\mathrm{eV}$ separation below vacuum,HenceThe energy and angle readings in the PDF are approximate. If a peak at the corner were instead assigned a small binding energy $E_B$, the estimate becomesFor $E_B=0$–$1\mathrm{eV}$ this gives about $57$–$60^{\circ }$, consistent with the near-Fermi feature in that angular range. Substituting the full photon energy for the kinetic energy would omit the work-function cost and underestimate the required angle.

Let $[a,a^{†}]=1$, $n=a^{†}a$, and take $S^{+}=(S^{-}{)}^{†}=\sqrt{2S-n}a$. The physical Holstein–Primakoff occupation constraint limits the Hilbert space, which contains only $n=0,1,\dots ,2S$. Acting in this number basis givesThe upper and lower endpoints are annihilated as required. ConsequentlyThe change in $n$ also gives $[S^z,S^{\pm }]=\pm S^{\pm }$. These matrix elements prove the spin algebra on the full physical spin-$S$ space. Extending the oscillator space to arbitrary occupation would not give a physical finite spin: the occupation constraint and operator ordering in the square root matter. bosons on distinct sites have $[a_m,a_n^{†}]={\delta }_{mn}$ and all other site commutators zero.

For one bosonic mode, define the unnormalized bosonic coherent state $|z⟩=e^{z{a}^{†}}|0⟩$. The commutator $[a,a^{†}]=1$ gives $a|z⟩=z|z⟩$ and $⟨z^{\prime }|z⟩=e^{{\bar{z}}^{\prime }z}$. The normalized state is $e^{-|z{|}^2/2}|z⟩$; it is not orthogonal to a different coherent state. The coherent-state resolution of identity follows from Gaussian moments,For many modes take products of these states and measures. The complex amplitudes have the natural scale of square root occupation; coherent states turn normal-ordered matrix elements into classical symbols without replacing the quantum trace by a single classical configuration.

For fermions, $\{c_i,c_j^{†}\}={\delta }_{ij}$ and a nilpotent annihilation operator cannot have an arbitrary nonzero ordinary complex eigenvalue. Introduce independent Grassmann variables ${\eta }_i,{\bar{\eta }}_i$, anticommuting among themselves and with the fermionic operators, and define $|\eta ⟩=\exp (-{\sum }_i{\eta }_i{c}_i^{†})|0⟩$. Then $c_i|\eta ⟩={\eta }_i|\eta ⟩$ and $⟨{\bar{\eta }}^{\prime }|\eta ⟩=e^{{\sum }_i{\bar{\eta }}_i^{\prime }{\eta }_i}$. With an ordered Berezin integral measure satisfying $\int d\bar{\eta }d\eta \bar{\eta }\eta =-1$ for one pair,for an even operator $O$; this is the fermionic coherent-state trace. In one mode the trace integrand is $e^{-\bar{\eta }\eta }(O_{00}-\bar{\eta }\eta O_{11})$, whose integral is $O_{00}+O_{11}$. Without the sign reversal it would be the supertrace, not the thermal trace. These formal Grassmann-valued coherent states do not represent ordinary complex-valued fermion amplitudes.

To construct the many-body thermal path integral, write $\mathcal{Z}=\mathrm{Tr}{e}^{-\beta K}$ with $K=H-\mu \hat{N}$ and divide the imaginary-time interval $\mathcal{T}=\beta \hslash$ into slices of width $\delta \tau$. Insert the appropriate coherent-state resolution of identity between factors. For a normal-ordered Hamiltonian, the adjacent-slice kernel iscoherent-state time slicing makes the overlap and Gaussian measure combine into the first-order time derivative. The grand-canonical coherent-state action in the regulated continuum limit isFor bosons the integration variables are complex, with the conjugate labels on the original integration contour; for fermions the barred and unbarred Grassmann variables are independent. The continuum notation retains the time-sliced ordering prescription. For example, $n(n-1)=a^{†2}{a}^2$ has normal symbol ${\bar{z}}^2{z}^2$, whereas $n^2$ has that symbol plus $\bar{z}z$. Dropping that distinction changes a theory. The coherent-state thermal boundary conditions follow because the bosonic trace closes paths periodically and the fermionic sign in the trace closes them antiperiodically:The barred fields obey the same respective boundary conditions. Hence Matsubara frequencies are $2\pi n/(\beta \hslash )$ and $(2n+1)\pi /(\beta \hslash )$. The derivative term encodes the canonical structure; omitting it would remove the quantum dynamics.^[ref]

As a check, a free diagonal mode with $K=\xi a^{†}a$ gives $Z_B=(1-e^{-\beta \xi }{)}^{-1}$ for $\xi >0$, whereas $K=\xi c^{†}c$ gives $Z_F=1+e^{-\beta \xi }$. Gaussian integration produces an inverse determinant for bosons and a determinant for independent complex fermions, with their respective temporal boundary conditions. Differentiating $\log Z$ recovers Bose-Einstein and Fermi-Dirac occupations. Interactions become higher powers in the normal-ordered symbols, and sources generate time-ordered correlations. The coherent-state construction preserves statistics, ordering and boundary conditions together; a formal continuum expression without its measure and temporal prescription is incomplete.

Define $\mathcal{M}$ by removing the factor $i(2\pi {)}^4{\delta }^{(4)}(p_1+p_2-p_3-p_4)$ from the connected transition matrix element. The single tree-level Feynman diagram in i givesThere is no identical-particle exchange term because the species differ. Each incoming fermion has two equally populated spin angular momentum states in an unpolarized beam. Their probabilities multiply to $1/4$, giving a spin average over the four initial choices. The unobserved final spin angular momentum states are distinct outcomes, so their probabilities are added by a spin sum, with no final averaging:Taking the complex conjugate of ${\bar{u}}_3{u}_1$ gives ${\bar{u}}_1{u}_3$. The fermion spin sum therefore converts each squared bilinear into a gamma-matrix trace:Using $\mathrm{tr}1=4$, $\mathrm{tr}{\gamma }^{\mu }=0$, $\mathrm{tr}({\gamma }^{\mu }{\gamma }^{\nu })=4{\eta }^{\mu \nu }$, and the mass shell conditions $p_j^2=m^2$, each gamma-matrix trace isThe second equality follows from $q^2=2m^2-2p_1\cdot p_3$; momentum conservation gives the identical result for $p_2,p_4$. HenceThe printed prefactor $g^2$ is incorrect for the stated coupling. Each of the two Yukawa interaction vertices contributes one $g$ to the scattering amplitude, so its square is proportional to $g^4$. As a concrete counterexample, take $g=2$, $m>0$, and forward elastic kinematics $p_3=p_1$, $p_4=p_2$. Then $q^2=0$ and the spin average is $16g^4=256$, whereas the printed expression gives $16g^2=64$.