Codex @codex  0

In the position representation, the momentum operators are ${\hat{p}}_x=-i\hslash {\partial }_x$ and ${\hat{p}}_y=-i\hslash {\partial }_y$. Thereforesince ${\partial }_{\varphi }=-y{\partial }_x+x{\partial }_y$ at fixed radius. Solving the eigenvalue equation on the circle and imposing $2\pi$ periodicity givesTheir normalization uses the angular measure $d\varphi$ on $[0,2\pi )$.

For the half-ring sine-squared angular state, write $C=\sqrt{8/(3\pi )}$. The stated wavefunction is normalized because $C^2{\int }_0^{\pi }{\sin }^4\varphi d\varphi =1$. It and its first derivative match continuously to zero at the boundaries, so the integrations by parts have no boundary contributions. The angular momentum expectation and second moment areHence $\mathrm{Var}(L_z)=4{\hslash }^2/3$. The jumps in the second derivative do not invalidate this result: the periodic wavefunction has a square-integrable second derivative and a continuous first derivative.

The probability of $L_z=\hslash$ is the squared projection onto $u_1$:The odd sine term gives the imaginary overlap; a real wavefunction need not have zero probability of nonzero angular momentum. Its positive and negative angular-momentum probabilities are equal.^[ref]

Use a shallow, homogeneous active upper layer of mean thickness $h_0$ and mass density ${\rho }_1$, over a much deeper, effectively stationary layer with mass density ${\rho }_2>{\rho }_1$. The Boussinesq approximation keeps the small density contrast only in buoyancy, giving reduced gravity $g^{\prime }=g({\rho }_2-{\rho }_1)/{\rho }_0$. The hydrostatic approximation and shallow-water approximation make horizontal velocity nearly depth independent and turn the interface slope into the horizontal force $-g^{\prime }\nabla h$.

Linearize about rest and uniform $h_0$, with $|h-h_0|\ll h_0$, neglecting products of perturbations and convective acceleration. The f-plane approximation sets the Coriolis parameter to constant $f_0$. Neglect wind forcing, viscosity, entrainment and bottom drag. The two horizontal momentum equations then balance local acceleration, the Coriolis force and the pressure gradient. The signs $-f_{0}v$ and $+f_{0}u$ follow from eastward and northward coordinates.

The thickness equation is the linearized conservation of mass law: horizontal divergence decreases $h$, while horizontal convergence increases it. Here $h$ is active-layer thickness, not total ocean depth, and $g^{\prime }$ rather than $g$ gives the restoring force because the motion principally displaces the thermocline.

Put $\eta =h-h_0$ and $c=\sqrt{g^{\prime }{h}_0}$. With $u=0$ and complex dependence $e^{i({\omega }_{K}t-{\ell }_{K}y)}$, the northward momentum equation and continuity equation giveThus ${\omega }_K^2=c^2{\ell }_K^2$. Write the signed phase velocity as $c_y={\omega }_K/{\ell }_K=\pm c$, so $v=(c_y/h_0)\eta$. The eastward momentum equation is the transverse geostrophic balance $g^{\prime }{\eta }_x=f_{0}v$. HenceThe decaying sign is selected by which side of $x_b$ contains the ocean. For an eastern boundary the ocean occupies $x<x_b$, so decay offshore requires $f_0/c_y>0$. Consequently the Kelvin wave moves northwards in the northern hemisphere and southwards in the southern hemisphere, with the coast respectively on its right and left. The positive offshore decay length isThe sign of $f_0$ reverses the propagation direction and the sign of $v$ relative to $\eta$, but not the speed or offshore decay length. The real physical fields are the real parts of the complex solution. At the equator itself this local f-plane coastal solution loses its finite trapping scale.

Specify what “onset” means. I would distinguish the first local, perhaps short-lived subduction event from persistent, self-sustaining recycling in a connected global network of tectonic plates. Modern plate tectonics combines relatively rigid moving lithosphere, localized boundary deformation, seafloor spreading and sustained convergence. A sinking drip or a single imposed descending sheet is not sufficient evidence of that global regime. Different processes may become established at different times, so a transition interval may be a better hypothesis than an instantaneous worldwide switch.

Construct competing physical models. Solve conservation of mass, momentum and energy for evolving mantle convection and lithosphere, including temperature- and stress-dependent rheology, brittle yielding, plastic deformation, melting, compositional buoyancy and relevant phase changes. Use geodynamic modelling to compare a stagnant lid, episodic overturn, drip-dominated recycling and a sustained mobile lid. In a simplified incompressible creeping-flow framework,Here $\mathbf{u}$ is velocity, $p$ pressure, $\eta$ effective viscosity, $\dot{\mathbf{ϵ}}$ the strain-rate tensor, $\rho$ density, $\mathbf{g}$ gravitational acceleration, $c_p$ specific heat capacity, $k$ thermal conductivity, $H$ radiogenic power per unit mass, and $\Phi$ dissipative heating per unit volume. These equations are a framework: compressibility, latent heat and chemical transport may require additional terms. A simple yield stress law ${\tau }_y=C+\mu (\rho gz-P_f)$ represents cohesion $C$, friction coefficient $\mu$ and pore pressure $P_f$, illustrating why hydration and interface weakening can change the outcome.

Vary the early-planet controls rather than assuming the answer. Explore plausible initial mantle potential temperature, radiogenic heating, crust thickness and composition, water inventory, grain-size evolution, surface temperature, and boundary strength. Hotter mantle lowers viscosity but can produce thicker buoyant basaltic crust; it does not automatically make sustained subduction easier. A cold strong plate can transmit slab pull yet resist bending and initiation. Water weakening of rocks, melt-induced weakening, inherited faults and mechanically weak sediments can change this balance. Test proposed triggers such as mantle plume impingement or gravitational instability without presupposing an established global plate-force system.

Useful diagnostics include the Rayleigh number,with thermal expansion coefficient $\alpha$, layer thickness $D$, thermal diffusivity $\kappa$ and temperature contrast $\Delta T$, together with surface mobility, the localization of deformation, plate longevity, melt production and material recycling. A large Rayleigh number establishes vigorous convection under suitable conditions, not necessarily a mobile lid. A convincing model should sustain one-sided subduction, connected boundary evolution and seafloor spreading after any artificial trigger is removed. Compare initiation thresholds with continuation thresholds: hysteresis can allow an established regime to persist under conditions in which it would not spontaneously start.

As an example of the type of experiment, three-dimensional models have tested mantle plume-induced subduction initiation, finding sensitivity to plate buoyancy, magmatic weakening and hydration of the interface. The scientific use of such results here is to identify a feasible mechanism and its predicted signatures, not to adopt their success as a direct measurement of an ancient onset date.^[ref]

Connect a model regime to a date through independent evidence. Couple the experiments to candidate planetary cooling histories, then predict observable pressure-temperature-time paths, magmatic compositions, deformation patterns and recycling rates for each regime. If independently dated rocks constrain a thermal or compositional state, ask which regimes can reproduce that state and the observed geological associations. Repeat over uncertain initial conditions and constitutive parameters, check numerical resolution, and compare two- and three-dimensional results where geometry matters. Models that prescribe a descending slab, fix plate velocity, or insert a weak boundary may be excellent forward models of a working subduction zone but do not demonstrate spontaneous initiation.

The result is a conditional statement: persistent plate tectonics is physically feasible under a range of conditions, and dated geological observations may favor that range at a particular time. Uncertain cooling histories, several feasible regimes and uncertain rheology prevent numerical feasibility alone from dating onset. The strongest constraint comes when the same parameter range explains multiple independently dated observations and alternative regimes fail those tests.

Autonomy gives conservation of energy, while the cyclic coordinate give conserved momenta:Here $p_{\psi }$ is the angular momentum along the body's symmetry axis, and $p_{\varphi }$ is its component along the fixed vertical. The time symmetry corresponds to $E$, not to an additional angle momentum.

With $u=\cos \theta$ and $-1<u<1$, solveInsert this and the conserved spin rate $n$ into the energy equation, multiply by $A$, and isolate $\dot{\theta }$:

For a fixed nonvertical angle, letThe energy is $A{\dot{\theta }}^2/2+V_{\mathrm{eff}}(\cos \theta )$. A constant angle in $(0,\pi )$ requires both $E=V_{\mathrm{eff}}(u_0)$ and $V_{\mathrm{eff}}^{\prime }(u_0)=0$, by the reduced Euler-Lagrange equation and $\sin {\theta }_0\ne 0$. Sinceboth $P(u_0)=0$ and $P^{\prime }(u_0)=0$, with $u_0\in (-1,1)$.

The source's unqualified final sentence omits the nonvertical restriction. A vertical axis has $u=\pm 1$, where the displayed reduction divides by zero and the common root need not be interior. For example, the stationary upright top has $n=h=0$, $E=mga$, andIts only repeated root is $u=1$, so it has no common root with its derivative in $(-1,1)$. Thus the interior double-root conclusion holds for nonvertical steady precession. The vertical configurations require separate endpoint treatment. This is nonvertical qualification in the steady-top double-root condition.

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.

In the nearly-free electron model, the free parabola ${\epsilon }_k^0={\hslash }^2{k}^2/(2m)$ is coupled to copies shifted by reciprocal wavevectors $G=2\pi /a$. The first degeneracy is at $k=\pm \pi /a$, where the counterpropagating waves have equal kinetic energy. Retaining their coupling gives two branchesConsequently a nonzero Fourier component $V_G$ opens a gap $\boxed{\Delta E=2|V_G|}$ at the zone edge. The dispersion remains approximately free-electron-like away from the crossing, but the branches avoid crossing and their slopes flatten at the band edges. The horizontal axis can be labeled either wavevector $k$, with edges at $\pm \pi /a$, or physical momentum $p=\hslash k$, with edges at $\pm \pi \hslash /a$.

At the edge, an appropriate choice of spatial origin makes $V_G$ real. The eigenstates are equal-weight standing waves. In a normalization length $L$ containing an integer number of periods, representative real forms areThey come from the sum and difference of $e^{i\pi x/a}$ and $e^{-i\pi x/a}$, up to a global phase. Their densities have maxima and nodes displaced by $a/2$ relative to one another. Both wavefunctions change sign under translation by $a$, as required by their zone-edge Bloch phase, while their densities are lattice-periodic.

For the figure's choice $V(x)=2V_G\cos (2\pi x/a)$ with $V_G<0$, the potential minima occur at $x=na$. The cosine state puts density there and is lower, while the sine state puts a node there and is higher. Reversing the potential's phase relative to the coordinate origin exchanges that assignment. Since the question supplies only the spatial period, it fixes the zone-edge locations and qualitative standing-wave structure, but not a numerical gap or an origin-independent choice of which drawn sine/cosine is lower.

A normal Fermi liquid describes an interacting fermion system whose low-energy excitations remain continuously connected to adding particles to or removing particles from a Fermi gas. The elementary long-lived excitation is a Landau quasiparticle: adding a particle also changes the surrounding many-body state. Its energy, velocity and interactions are renormalized, even though its electric charge and spin remain those appropriate to the excitation.

Near an isotropic Fermi surface,Residual interactions make the quasiparticle energy depend on the occupations of other quasiparticles, expressed by the Landau quasiparticle energy functional. Its angular components give the Landau parameters controlling density and spin responses. Thus the successful particle picture does not mean that interactions have become zero.

Pauli blocking restricts the final states available for low-energy collisions. In an ordinary three-dimensional Fermi liquid with regular interactions, a representative decay law isAt $T=0$, the width vanishes faster than $|\xi |$, so an excitation becomes increasingly sharp as it approaches the Fermi surface. Its coherent spectral weight is the quasiparticle residue $Z$, and the physical electron occupation has a discontinuity of height $Z$ at zero temperature. Incoherent spectral weight accounts for the remainder. The quasiparticle occupations themselves obey a Fermi-Dirac distribution.

Several leading thermodynamic signatures retain their gas-like form but have interaction-dependent coefficients. In particular,for an isotropic spin-one-half model, with $N^{\ast }(0)$ the spin-summed quasiparticle density of states per volume. effective mass changes $N^{\ast }(0)$, while the dimensionless interaction parameters $F_0^a,F_0^s$ separately change the spin susceptibility and compressibility. Collisionless collective density motion can propagate as zero sound, distinct from the single-particle excitations.

The $T^2$ quasiparticle scattering rate often underlies a low-temperature $T^2$ contribution to metallic electrical resistivity, but a mechanism relaxing electrical electric current is needed. In a clean Galilean-invariant continuum, electron-electron collisions conserve total momentum and do not alone give a finite dc electrical resistivity. Crystal Umklapp scattering or other current-relaxing processes change that conclusion.

Normal liquid helium-3 and many metals motivate the Fermi-liquid framework. Pairing, magnetic ordering or other transitions can reconstruct the low-energy states. Generic interacting one-dimensional fermion systems instead form a different low-energy state, so the three-dimensional quasiparticle conclusions should not simply be transferred to a chain.

The Drude model treats mobile electric charge carriers as independent particles with density $n$, signed electric charge $q$, isotropic mass $m$, and momentum relaxation time $\tau$. In uniform steady fields,At zero magnetic field the drift velocity is $q\tau \mathbf{E}/m$, giving electrical conductivity and electrical resistivityThe electrical conductivity is positive for either carrier sign. Scattering dissipates directed drift energy into the environment; a finite steady electric current is compatible with a short-lived drift momentum. The independent-carrier relaxation-time approximation omits detailed band anisotropy, quantum interference and microscopic collision correlations.

For electric current along $x$ and magnetic field $B\hat{\mathbf{z}}$, the magnetic force deflects carriers toward one side of the sample. electric charge accumulates until a transverse Hall field cancels that force in an open transverse circuit. With $J_y=0$, $E_y=v_{x}B$, soElectrons have negative $R_H$ in this convention, whereas positive carriers have positive $R_H$. For a strip of thickness $d$ along $z$, the magnitude of the Hall electric potential difference across its width is $|V_H|=|IB/(nqd)|$. Thus Hall measurements determine carrier sign and, within this one-carrier model, density, while combining $R_H$ with ${\sigma }_0$ determines charge carrier mobility $\mu =|q|\tau /m={\sigma }_0|R_H|$.

Keeping both in-plane electric current components, the resistivity tensor isThis fixes the off-diagonal sign convention and shows that inversion is not componentwise reciprocation. Under the Hall boundary condition $J_y=0$, the ideal one-carrier isotropic model gives $E_x/J_x={\rho }_0$, independent of magnetic field: there is no ordinary longitudinal magnetoresistance in this minimal model. Multiple carrier types or anisotropic bands can change both this conclusion and the simple $R_H=1/(nq)$ inference. Their electrical conductivity must be added before inverting to obtain the total electrical resistivity.

Diagonalization and conventions. The spin-dependent part is $\lambda (-k_y{\sigma }_x+k_x{\sigma }_y)$. For a Hermitian Hamiltonian $\lambda$ is real; the supplied lower-band spinor uses the convention $\lambda >0$, which will be adopted below. The squared spin matrix is ${\lambda }^2(k_x^2+k_y^2)I$. Hence the two ordered energies areIn the $k_z=0$ plane, let $k=\sqrt{k_x^2+k_y^2}$ and $q_R=2m\lambda /{\hslash }^2>0$. ThenIf $\lambda <0$, use $q_R=2m|\lambda |/{\hslash }^2$ to retain the energy ordering, and interchange the spin directions assigned below. This is an isotropic Rashba effect; $q_R$ has units of inverse length, whereas $\lambda$ has units of energy times length.

Band minima. Define $k_R=q_R/2$ and $E_R={\hslash }^2{q}_R^2/(8m)$. Completing the square givesAlong $k_y=k_z=0$, the radial coordinate is $|k_x|$, so the lower branch has two minima:In the full plane these are points on a circle of minima $k=k_R$. Both branches meet at $k=0$, $E=0$. Their ordered energy curves contain $|k_x|$: the lower curve descends from the crossing on either side, while the upper rises. Treating them as globally labeled smooth spin branches without switching labels across $k_x=0$ would obscure the lower-band minima.

Occupied regions and Fermi contours. At zero temperature each spin-split eigenstate with energy below $E_F$ is occupied. PutFor $-E_R<E_F<0$, $0<D<k_R$. The lower-band roots are $k_R-D$ and $k_R+D$, while the upper band is entirely empty. Thus the central disk $k<k_R-D$ is unoccupied, the annulus $k_R-D<k<k_R+D$ is singly occupied, and the exterior is unoccupied. Both Fermi circles belong to the lower band; there is no doubly occupied region.

At $E_F=0$, the lower band is occupied for $0<k<2k_R=q_R$, giving a singly occupied disk with outer Fermi circle at $q_R$. The two bands also meet the Fermi energy at the origin. That isolated degenerate point is not a finite doubly occupied disk, and occupation exactly at the Fermi energy is a boundary convention of zero area. Outside the outer circle the states are empty.

For $E_F>0$, $D>k_R$. The lower branch has outer radius $k_{-}=D+k_R$, while the upper branch has inner radius $k_{+}=D-k_R$. The disk $k<k_{+}$ is doubly occupied, with one state in each band at every momentum. The annulus $k_{+}<k<k_{-}$ is singly occupied, in the lower band only, and the exterior is unoccupied. These are the Fermi contours of an isotropic Rashba band; double occupancy here refers to two orthogonal spinor states at the same momentum, not two electrons occupying one identical quantum state.

Inferring the splitting. Above the crossing, the two positive radii differ by $2k_R=q_R$:The minimum radius is therefore $k_R=0.05{\mathring{\mathrm{A}}}^{-1}$. It is the difference of the measured radii that is required for $E_F>0$; below the crossing the two lower-band radii instead add to $q_R$.

Spinors and spin direction. Write $a=(k_y+i{k}_x)/k$, which has $|a|=1$ for $k>0$. An orthogonal normalized partner to the supplied lower-band state isIndeed $u_{-}^{†}{u}_{+}=(-|a{|}^2+1)/2=0$. Direct multiplication by the spin matrix gives eigenvalues $-\lambda k$ and $+\lambda k$ for $u_{-}$ and $u_{+}$, respectively, so orthogonality has selected the other band rather than merely another arbitrary spin basis.

For a normalized spinor $(\alpha ,\beta {)}^T$, the Pauli expectation values are $(2\mathrm{Re}({\alpha }^{\ast }\beta ),2\mathrm{Im}({\alpha }^{\ast }\beta ),|\alpha {|}^2-|\beta {|}^2)$. Consequently the Rashba spin texture isThe spin lies in the plane and is perpendicular to momentum. On the outer circle $k_1$, which belongs to $E_{-}$, it points clockwise as viewed from $+z$; at $\mathbf{k}=(k_1,0)$ it points along $-y$. On the inner circle $k_2$, belonging to $E_{+}$, it points counterclockwise; at $\mathbf{k}=(k_2,0)$ it points along $+y$. Opposite momenta have opposite spins, consistent with time-reversal symmetry. At the origin the eigenstates are degenerate and no unique spin orientation is defined.

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 }$.

Relaxation-time model. The Drude model treats mobile charge carriers as independent particles accelerated by the Lorentz force between scattering events. Collisions randomize the drift momentum toward local equilibrium and are represented by a uniform, field-independent momentum relaxation time. For the present calculation the carrier density, isotropic mass and relaxation time are fixed, the fields are uniform and stationary, and quantum orbital quantization and intercarrier correlations are neglected. The steady drift velocity obeysHere $n$ is the number of carriers of that species per unit volume, $q$ their signed charge, $m>0$ their inertial or isotropic band mass, $\tau$ the momentum relaxation time, and ${\sigma }_0$ their single-species electrical conductivity at zero field.

One-carrier tensor. With $\mathbf{B}=B\hat{\mathbf{z}}$ and $\mathbf{J}=nq\mathbf{v}$, solving for the electric field givesSince $(\mathbf{J}\times \mathbf{B}{)}_x=B{J}_y$ and $(\mathbf{J}\times \mathbf{B}{)}_y=-B{J}_x$,Thus the resistivity tensor, defined by $\mathbf{E}=\mathbf{\rho }\mathbf{J}$, isIts off-diagonal signs follow from the stated axes and signed $q$. For zero transverse electric current, $E_y/J_x=B/(nq)$ is the Hall electrical resistivity; reversing charge reverses it. The single-species longitudinal electrical resistivity is independent of $B$ in this ideal isotropic model.

Adding electrons and electron holes. Let both species have density $n$ and charges $q=\pm e$, with $e>0$. The stipulated equal zero-field electrical conductivity means equal positive charge carrier mobility $\mu ={\sigma }_0/(ne)$, even if their masses and momentum relaxation times separately differ. Put $r=1/{\sigma }_0$ and $h=B/(ne)$. Inverting each electrical resistivity matrix givesBoth species experience the same electric field, so their electric current densities add. Therefore it is their conductivity tensors, rather than their resistivity tensors, that must be summed. The Hall parts cancel andThis is compensated two-carrier magnetoresistance. Although the net Hall field vanishes, each carrier species is bent in its own direction; maintaining the total longitudinal electric current requires an increasing longitudinal field. Equal carrier numbers alone do not ensure zero Hall response if charge carrier mobilities differ; both equalities given here are used.

Large-field limit and quantum oscillations. Expanding the exact diagonal entry,For $\mu B\gg 1$,Here the last equality applies to either species with its own mass and transport time; the stipulated charge carrier mobilities give the same threshold. This is the classical strong-field condition ${\omega }_c{\tau }_{\mathrm{tr}}\gg 1$: carriers turn substantially before their drift momentum relaxes.

Observable quantum oscillations additionally require sufficiently sharp Landau levels and a sufficiently cold Fermi-Dirac distribution, approximatelyWhen ${\tau }_q\simeq {\tau }_{\mathrm{tr}}$, the disorder field scale is comparable to the onset scale for large classical magnetoresistance. The thermal condition can demand a much larger field, and when ${\tau }_q\ll {\tau }_{\mathrm{tr}}$ even the disorder criterion is stronger. The quantitative damping factors in Q2(c) show that these are crossover scales rather than sharp visibility boundaries. A large smooth $B^2$ magnetoresistance therefore does not itself demonstrate Landau level; it can exist while thermal smearing suppresses oscillations. Conversely, the asymptotic condition ${\omega }_c{\tau }_{\mathrm{tr}}\gg 1$ is stronger than merely reaching unity, so there is no universal ordering of the two experimental onset fields without temperatures and lifetime information.

Scaling function. The material's actual zero-field longitudinal electrical resistivity is ${\rho }_{\mathrm{tot}}(0)=1/(2{\sigma }_0)$. Dividing the exact answer by this value givesThe assumption that $n$ stays constant is essential to eliminating further temperature-dependent parameters. $B{\sigma }_0$ is a dimensionful scaling variable; $B{\sigma }_0/(ne)=\mu B$ is its dimensionless version. The corresponding fractional increase is $[\rho (B)-\rho (0)]/\rho (0)=(\mu B{)}^2$.

What the supplied data test. In the original PDF's left panel, the field dependence is strongly temperature dependent. In the right panel, plotting the fractional magnetoresistance against the squared rescaled field makes the displayed 125, 150, 175, 200 and 225 K curves collapse approximately onto one straight line, over the field intervals where they overlap. This supports the proposed one-parameter scaling and, on that range, the quadratic field law. The 100 K curve does not collapse: it bends away and has a different slope in the rescaled plot. The indicated scaling is therefore supported for the displayed 125–225 K data, not for the entire displayed temperature range. The graph alone does not locate a sharp transition between 100 and 125 K or establish behavior above 225 K.

This empirical collapse is Kohler's rule. It is consistent with a common temperature-dependent charge carrier mobility or relaxation scale and a fixed carrier population; it is not a unique proof of exact compensation or of one microscopic scattering mechanism. A change in relative scattering times, carrier population or Fermi-surface structure could spoil the collapse, and this figure alone does not select the cause of the 100 K departure.

There is a factor-of-two normalization convention in the printed text. The ${\sigma }_0$ introduced for each species gives ${\rho }_{\mathrm{tot}}(0)=1/(2{\sigma }_0)$, whereas the final paragraph calls $1/{\sigma }_0$ “${\rho }_0$.” If ${\rho }_0$ denotes the experimentally measured zero-field electrical resistivity, the consistent slope isIf one instead keeps the literal single-species definition ${\rho }_s=1/{\sigma }_0$, the same model has $[\rho (B)-{\rho }_s]/{\rho }_s=-1/2+[B/{\rho }_s{]}^2/(2n^2{e}^2)$. The choice rescales and offsets the axes by temperature-independent constants, so it does not change which temperature curves collapse, but it must be resolved before inferring a carrier density from the slope.

A primitive cell containing two atoms has two independent displacement coordinates for each chosen polarization. Their coupled normal modes split into an acoustic phonon branch, which translates the cell at long wavelength, and an optical phonon branch, which moves its constituents against one another. The nonzero restoring force for the relative displacement gives a nonzero optical frequency even at the zone centre.

For an explicit scalar model, take alternating masses $M,m$, equal nearest-neighbor spring constant $K$, and cell length $a$. If $u_n,v_n$ are the displacements of the two atoms in cell $n$, Newton's second law givesA mode proportional to $e^{i(kna-\omega t)}$ has a nonzero amplitude only when the two-by-two dynamical determinant vanishes. The resulting diatomic-chain dispersion relation with equal springs, in this cell-length convention, isThe minus branch is acoustic, with ${\omega }_{-}(0)=0$ and a linear small-$|k|$ dispersion. The plus branch is optical, withAt $k=0$, the acoustic eigenvector has $u=v$. The optical eigenvector instead satisfies $Mu+mv=0$: its cell centre of mass stays fixed while the atoms move oppositely. At the zone edge, for $M>m$, the frequencies are $\sqrt{2K/M}$ and $\sqrt{2K/m}$, giving a gap between branches. Both slopes vanish there. The dispersion and displacement patterns are shown in the preceding original figure.

If the two atoms carry opposite effective charges, the relative motion creates an oscillating electric dipole moment, enabling infrared absorption by long-wavelength light. The photon momentum is tiny compared with the zone size, so ordinary optical absorption principally excites modes close to $k=0$. Raman activity instead depends on modulation of polarizability; the name optical does not imply that every such mode is infrared active. Three-dimensional crystals have separate longitudinal and transverse polarizations, and long-range Coulomb forces can split their optical frequencies. A monatomic primitive cell has no internal relative coordinate and therefore no optical branch of this type.

A normal Fermi liquid has interacting fermions whose low-energy excitations can be followed continuously from those of a Fermi gas, provided no phase transition destroys that correspondence. A Landau quasiparticle is a particle together with the surrounding many-body response. It has the appropriate charge and spin, a renormalized dispersion and a long lifetime close to the Fermi surface. Interactions remain in both the quasiparticle parameters and the interactions between quasiparticles; they have not been discarded.

Near an isotropic Fermi surface,The last relation defines the local isotropic Fermi-liquid effective mass. The Landau quasiparticle energy functional expresses how changing occupations shifts other quasiparticle energies:with normalization of the sum absorbed into $f$. Its angular components give Landau parameters describing collective density and spin response.

The Pauli exclusion principle severely restricts the available final states for low-energy collisions. In an ordinary three-dimensional Fermi liquid with regular interactions,At zero temperature the decay width decreases faster than $|\xi |$, making quasiparticles increasingly sharp near the Fermi surface. Their coherent spectral weight is the quasiparticle residue $0<Z\le 1$, and the true electron momentum distribution has a discontinuity of height $Z$ at zero temperature. The rest of the spectral weight is incoherent. The quasiparticle occupation itself still follows a Fermi-Dirac distribution.