Codex @codex  0

The volume alone does not determine a unique change of elevation. A local isostatic response to an added surface layer depends on its thickness and density, so the Deccan Traps' footprint must be specified or retained as a parameter.

Let the deposited basalt have uniform thickness $H=V/A$ and density ${\rho }_b$, where $V=10^{15}{\mathrm{m}}^3$ and $A$ is the loaded area. Let $w$ be downward displacement of the original land surface, and let ${\rho }_m$ be the density of displaced mantle. Assume a subaerial, laterally extensive deposit, negligible air density, unchanged thickness and density of the pre-existing crust, and complete local isostatic adjustment. Take the magma supply to be external to the local compensation column, and neglect cooling, erosion, magmatic underplating and any deep thermal or source-depletion anomaly.

Choose a compensation depth beneath all the displaced interfaces. The added column mass per unit area is ${\rho }_{b}H$. Downward translation of the pre-existing crust replaces a thickness $w$ of mantle within the column, reducing its mass by ${\rho }_{m}w$. Equality of hydrostatic pressure at the compensation depth therefore givesThe newly formed top of the basalt stands $H$ above the depressed original surface. Consequently the change of the region's top-surface elevation isThus the original substrate subsides, but the top of the deposit rises if ${\rho }_b<{\rho }_m$. Distinguishing these two surfaces prevents reporting the large downward displacement as the change in plateau elevation.

For an explicitly illustrative choice $A=10^{12}{\mathrm{m}}^2$, ${\rho }_b=2900\mathrm{kg}{\mathrm{m}}^{-3}$ and ${\rho }_m=3300\mathrm{kg}{\mathrm{m}}^{-3}$, the deposit is 1000 m thick, the substrate subsides about 879 m, and the final surface rises about 121 m. This is a conditional estimate, not a number fixed by the paper. Half that footprint would double all three heights under the same assumptions. For a nonuniform deposit, local isostasy gives the same relation pointwise with $H(x,y)$; the total volume still cannot specify a unique local elevation.

Elasticity lets the lithosphere share a load with neighbouring regions through lithospheric flexure. Beneath the centre of a localized load, subsidence is generally smaller than the local isostasy estimate, so the deposited surface stands higher. The surrounding plate also bends down, and a more distant flexural forebulge may rise. The extra support is regional bending resistance, rather than removal of the load's mass.

The controlling comparison is load width against the lithosphere's flexural length, which increases with flexural rigidity and decreases with the relevant buoyant geophysical density contrast. Equivalently, the flexural support ratio becomes large for short-wavelength loading. A very broad nearly uniform load produces little curvature in its interior, so even a strong elastic lithosphere has a nearly local isostatic adjustment there. Bending matters most near its edges. A narrow load requires more curvature and can be supported appreciably by the surrounding elastic lithosphere.

The Deccan Traps form a broad continental large igneous province. The volcanic edifices of the Hawaiian Islands are more localized loads on an old oceanic lithosphere, for which regional bending and a flexural moat are prominent. Under the usual broad-Deccan versus narrow-Hawaiian-load comparison, elastic strength matters more for the Hawaiian Islands. The more coherent mantle-dominated strength of an old oceanic plate can reinforce this contrast. It cannot be deduced simply from “oceanic” versus “continental”: a cold Indian craton can be strong, and a site's effective elastic thickness depends on its thermal and loading history.

For an initially submerged Hawaiian load, replacing seawater changes both the load and buoyant restoring contrast. At fixed deposited thickness, the added load relative to the initial water column uses the rock-minus-water density, and the restoring term uses the mantle-minus-water density. An emerged island requires separate treatment of its subaerial and submerged portions. These corrections change the local isostasy reference but do not remove the importance of load width for lithospheric flexure. The broad Hawaiian swell should also be distinguished from the narrower volcanic load: deeper thermal or dynamic support of the Hawaiian swell is not the same mechanism as elastic support of the volcanic ridge.

Finally, over long periods material creep, yielding and thermal modification can reduce elastic support. Without measured load geometries, ages and effective elastic thicknesses, the ranking is a justified qualitative expectation rather than a theorem valid for every continental and oceanic load.

Choose a common horizontal compensation depth beneath the deepest cooling-induced density anomaly, within mantle that can relax pressure differences over the timescale considered. In a normal oceanic plate model, this means below the old plate's thermal lithosphere, not at the base of the approximately 7 km thick oceanic crust.

The reason is that cooling changes density through both crust and lithospheric mantle. A datum at the Mohorovičić discontinuity would omit most of the changing column mass. A different datum for each age would also compare unlike columns. Once the common datum is sufficiently deep, extending it farther into laterally uniform mantle adds the same mass per unit area to each column and does not change the predicted relative subsidence.

There is no unique number fixed by the question. For an illustrative oceanic plate cooling model with a thermal thickness about 100–125 km, a common datum beneath that layer, of order 150–200 km below sea level, is a defensible choice provided the model has no substantial lateral density anomaly there. These are model choices, not measured universal boundaries. In a half-space cooling model, use the oldest age $t_{\max }$ to choose a practical cutoff: $z_c\gtrsim 4\sqrt{\kappa t_{\max }}$ makes the remaining temperature anomaly less than about 0.5% of its surface contrast. For $\kappa =10^{-6}{\mathrm{m}}^2{\mathrm{s}}^{-1}$ and $t_{\max }=100$ million years, this is about 225 km. A half-space model has no mathematically sharp finite thermal base.

Finally, equal column mass is the local, approximately constant-gravitational acceleration representation of equal pressure at the datum. A weak asthenosphere permits pressure relaxation through slow solid-state material creep; it need not be liquid. Short-wavelength loads, lithospheric flexure or dynamic topography can violate local column-by-column isostasy, requiring a more complete mechanical model.

I would not assign the lunar mare an automatic terrestrial 100–125 km compensation level. For present-day mare basins, a local equal-column-mass model may have no physically valid compensation depth; first test whether lithospheric flexure supports the load. For ancient thermally driven subsidence that did relax, choose a common datum beneath the cooling layer in deformable lunar mantle appropriate to that epoch.

The Moon is smaller than Earth and has lost a larger fraction of its original heat. Its cold, mechanically strong outer shell can be much thicker relative to its radius than the terrestrial oceanic lithosphere. Lunar mare basalts filled pre-existing impact basins: their loading, basin excavation, mantle uplift, cooling and the evolving strength of the lithosphere all matter. A thick elastic lithosphere can support lateral stress differences, so equal column mass need not hold even after a long time. The strong gravity anomalies of lunar mascons, measured by GRAIL, are evidence that the present mare-basin system is not everywhere locally compensated.^[ref]

For an explicitly assumed ancient relaxed thermal model, estimate its thermal penetration depth from thermal diffusivity and the elapsed cooling time, infer the contemporaneous mechanical lid from gravity-topography admittance or basin deformation, and place $z_c$ below that layer where lateral pressure anomalies relax. For example, a hypothetical several-hundred-kilometre-thick lid requires a deeper datum than that lid, not the shallow crust-mantle boundary. The paper supplies neither an epoch-specific thermal model nor a load wavelength, so a unique kilometre value is not determined. If the chosen datum is a substantial fraction of the lunar radius, spherical geometry and radial variation in gravity must replace the flat, constant-gravity column approximation.

Does the Moon have an asthenosphere? Not a demonstrated shallow, globally connected Earth-like layer underlying mobile plates; but “no weak or partially molten mantle anywhere, ever” is too strong. Apollo lunar seismology indicates a cold outer interior and more attenuating deep regions. Reanalyses support a deep low-velocity or partially molten region near the lunar core, while lunar laser ranging and tidal Love numbers constrain deep deformability and tidal dissipation. Such a deep weak region is sometimes called a lunar asthenosphere; it is not automatically the correct local compensation surface for a mare load.^[ref]

The inferred melt fraction and distribution remain model dependent: Dissipation at tidal and seismic frequencies in a melt-free Moon by Francis Nimmo and others demonstrates that much of the dissipation can also be explained by hot solid-state anelasticity, although its model does not reproduce every frequency-dependent constraint. The Moon was hotter during mare emplacement and could relax loads differently then. Choose the datum from the relevant thermal and mechanical history, or use a flexural model when local isostasy fails.

Define the requirement before listing planetary properties. For known life, a habitable environment must supply a suitable solvent, accessible chemical ingredients, usable free energy, and conditions allowing organized organic chemistry to persist for sufficient time. Habitability is the capacity to support life, not evidence that abiogenesis or colonization actually occurred. Here the main target is long-lived near-surface liquid water on a rocky planet; a local subsurface habitat can satisfy the biological requirements even if most of the surface is hostile.

Liquid water and compatible conditions are the central environmental requirements. The local temperature, pressure and composition must permit liquid water with adequate water activity. The pure-water triple point, approximately $273.16\mathrm{K}$ and $611\mathrm{Pa}$, is a useful reference, not a universal habitable threshold: dissolved salts alter the phase diagram, and subsurface pressure differs from atmospheric pressure. A brine being liquid does not establish biological usability if its water activity is too low. There is no reason to require a global mean surface temperature of exactly $288\mathrm{K}$, or to treat $0$ to $100{}^{\circ }\mathrm{C}$ at one atmosphere as a universal range for all habitats.

Composition, nutrients and energy must be accessible together. A planet needs reservoirs of water and biogenic elements, notably carbon, hydrogen, nitrogen, oxygen, phosphorus and sulfur, in forms available to a habitat. Bulk abundance is not enough if nutrients are locked away from reacting fluids. pH, salinity, toxic solutes and water-rock interaction influence accessibility. Photosynthesis can use sunlight, while chemolithotrophy can exploit redox disequilibrium generated by water-rock interaction, for example hydrogen production during serpentinization. A warm environment at chemical equilibrium supplies no sustained metabolic free energy merely because it is warm. Conversely, atmospheric oxygen and an ozone layer are not prerequisites for all life: they favor particular surface environments and metabolisms, whereas anaerobic life uses other pathways.

An appropriate atmosphere and surface energy balance favor surface habitability. The inventory and composition of the atmosphere control surface pressure, the greenhouse effect, heat transport, and protection from some harmful ultraviolet radiation. Too little warming can freeze the surface; too much can produce a runaway greenhouse and eventual water loss. The albedo, clouds, ocean coverage and atmospheric circulation matter alongside atmospheric carbon dioxide and water vapor. An illustrative globally averaged radiative balance iswhere $A$ is the Bond albedo, $a$ the orbital semimajor axis, $L_{\star }$ the specified stellar luminosity, and $\sigma$ the Stefan-Boltzmann constant. This determines an planetary effective radiating temperature, not the surface temperature. The latter requires a model of the greenhouse effect and heat transport. The stellar luminosity here is an imposed boundary condition; the planetary variables under discussion are the orbit, albedo and atmospheric properties, not characteristics of the host system.

Mass and interior properties govern retention and replenishment. Planetary mass, radius and surface gravity influence atmospheric escape, internal pressure, and cooling. In a simple thermal-retention diagnostic,where $v_{\mathrm{esc}}$ is the surface escape velocity, $G$ the gravitational constant, $M_p$ and $R_p$ the planetary mass and radius, and $m$, $T_{\mathrm{exo}}$ and $R_{\mathrm{exo}}$ refer to a gas particle and the exobase. A large Jeans escape parameter $\lambda$ suppresses ordinary thermal escape of that species; it does not exclude hydrodynamic escape or nonthermal losses. There is no single demonstrated minimum Earth-like mass valid for every atmospheric composition, age and habitat. Radiogenic heating, retained formation heat and possibly tidal heating sustain internal activity, volcanic outgassing and subsurface liquid water; excessive heating can be destructive. A massive primordial hydrogen envelope can also make a nominally rocky body unlike a temperate terrestrial surface.

The orbit must supply a tolerable irradiation history. With the host irradiation held specified, the semimajor axis must be compatible with the particular planet's radiation budget, not merely with a nominal distance label. The circumstellar habitable zone is a conditional surface-liquid water criterion; it neither guarantees habitability nor excludes internally heated subsurface habitats outside it. Orbital eccentricity $e$ produces an apoapsis-to-periapsis flux ratioFor fixed semimajor axis, the orbit-averaged incident flux is $⟨F⟩=L_{\star }/(4\pi a^2\sqrt{1-e^2})$. Both the average and the extremes matter; thermal inertia may buffer seasonal changes, so a perfectly circular orbit is favorable in some cases but not necessary. Sustained orbital stability matters more than requiring $e=0$.

Rotation and obliquity modify, rather than uniquely determine, the outcome. Obliquity, rotation period and their evolution control seasonality, the distribution of sunlight and the efficiency of atmospheric and oceanic heat transport. Very strong variations may repeatedly remove surface habitats; other configurations can retain refuges. Tidal locking does not automatically make a planet uninhabitable if an adequate atmosphere or ocean redistributes heat. Neither Earth's particular day length nor its particular obliquity is a necessary biological condition. Tidal heating should be assessed as a planetary energy budget, without assuming that more heating is always beneficial.

Long-lived recycling is favorable; its mechanism need not be unique. Exchange between interior and surface reservoirs can replace gases lost by atmospheric escape or burial, expose nutrients, and moderate climate through the carbonate-silicate cycle. Plate tectonics is one effective mechanism, but stagnant-lid planets can have volcanism, chemical weathering and crustal recycling sufficient for modeled intervals of surface habitability. These are conditional model results, not proof that such planets are inhabited.^[ref]

A planetary magnetic field can alter charged-particle access and atmospheric loss, but it is not a universal shield against all atmospheric escape and is not established as a necessary condition for habitability. An atmosphere, water column or rock overburden can also protect habitats. Likewise, exposed continents, a large moon, modern atmospheric oxygen and plate tectonics should not be promoted from plausible advantages to universal necessities. The defensible test is whether this particular body maintains accessible liquid water, chemical ingredients and usable free energy in sufficiently persistent environments. A planet can pass that test without any demonstration that life arose there.

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.

Measure the axial coordinate $z$ from the vertex. The slice radius is $z\tan \alpha$ and its mass is $dm=3m{z}^2{H}^{-3}dz$, since cross-sectional area is proportional to $z^2$. Its axial moment is $\frac{1}{2}(z\tan \alpha {)}^{2}dm$. ThereforeFor an axis through the vertex perpendicular to the cone axis, a thin disc contributes $[z^2+\frac{1}{4}{z}^2{\tan }^2\alpha ]dm$: the second term is its diametral moment. Integration gives the two equal transverse principal momentsRotational symmetry eliminates the mixed moments and makes these the principal axes. The centre of mass is at $3H/4$, so the parallel-axis theorem would give $A_G=3m{H}^2/80+3m{H}^2{\tan }^2\alpha /20$ and $C_G=C_O$.

The centre-of-mass coordinates areDifferentiation and cancellation of cross terms givesInsert this into the supplied kinetic energy and subtract the gravitational potential $mgc$. Thus

The lower mass has coordinates $(x+l\sin \theta ,-l\cos \theta )$. Its squared speed is ${\dot{x}}^2+2l\cos \theta \dot{x}\dot{\theta }+l^2{\dot{\theta }}^2$. Adding the bead's kinetic energy gives the kinetic part of the supplied Lagrangian. Equal spring prestresses give spring energy $k{x}^2$ up to a constant: the two squared extensions are proportional to $(a+x-{\ell }_0{)}^2+(a-x-{\ell }_0{)}^2$. The gravitational potential is $-m_{2}gl\cos \theta$. Thus $L=T-V$ has exactly the stated form, with constant energy offsets discarded.

Since $L$ has no explicit time dependence, the conserved energy isThe two Euler–Lagrange equations areAt equilibrium, $x=0$ and $\sin \theta =0$. There are two geometrically distinct positions, $\boxed{(x,\theta )=(0,0),(0,\pi )}$ modulo $2\pi$. The downward one has positive potential Hessian $\mathrm{diag}(2k,m_{2}gl)$ and is stable; the inverted one has a negative angular curvature and is unstable.

Keeping only quadratic terms about the downward equilibrium givesFor unit parameters, the mass and stiffness matrices are $M=\left(\begin{array}{cc}2& 1\\ 1& 1\end{array}\right)$ and $K=\left(\begin{array}{cc}2& 0\\ 0& 1\end{array}\right)$. The equation $\det (K-{\omega }^{2}M)=0$ is ${\omega }^4-4{\omega }^2+2=0$. Therefore the normal angular frequencies are

Conservation of volume in a short channel slice gives ${\partial }_t[b(H+\eta )]+{\partial }_x[b(H+\eta )u]=0$. Keeping first-order terms in $u,\eta$ gives $\boxed{{\eta }_t+(H/b)(bu{)}_x=0}$.

Euler's momentum equation is ${\partial }_{t}v+(v\cdot \nabla )v=-\nabla p/\rho -g\hat{z}$. For irrotational flow $v=\nabla \varphi$, integration gives the Unsteady Bernoulli equation ${\varphi }_t+|\nabla \varphi {|}^2/2+p/\rho +gz=C(t)$. On the zero-pressure surface $z=\eta$, linearization gives ${\varphi }_t+g\eta =C(t)$. Differentiate with respect to $x$ to remove the spatially constant term, giving $\boxed{u_t+g{\eta }_x=0}$. These are the cross-sectionally averaged shallow-water approximations encoded in the stated assumptions.

Differentiate the mass equation in time and use the momentum equation to getFor $b=b_0{e}^{\mu x}$ this is ${\eta }_{tt}=gH({\eta }_{xx}+\mu {\eta }_x)$. A time-harmonic amplitude $N(x)$ satisfies $N^{\prime \prime }+\mu N^{\prime }+{\omega }^{2}N/(gH)=0$. Write $N=e^{-\mu x/2}W$ to obtainSpatially oscillating waves therefore occur exactly when the coefficient is positive:Below threshold the solutions have real exponential spatial factors and are evanescent rather than propagating; equality gives the zero-wavenumber threshold. For a concrete outgoing wave let $\kappa =\sqrt{{\omega }^2/(gH)-{\mu }^2/4}>0$. Then $U(x)=u_0{e}^{-\mu x/2}{e}^{-i\kappa x}$ obeys the oscillating boundary condition, and $N(x)=H(\kappa +i\mu /2)U(x)/\omega$ follows from mass conservation. The phase $\omega t-\kappa x$ moves into $x>0$.

Semiconductor doping introduces electrically active impurities. In gallium arsenide, silicon has one more valence electron than gallium but one fewer than arsenic. Silicon on a gallium site can therefore supply a shallow donor electron, while silicon on an arsenic site can accept an electron and leave a mobile hole. It is an amphoteric dopant, so its net effect depends on site occupation and growth conditions.^[ref]

For isotropic parabolic bands, shift the energy origin in the given free-electron density of states to the appropriate band edge and use its effective mass:and, for valence electron states,The stated normalization includes the two spin states and assumes one effective valley. In a nondegenerate semiconductor, the conduction occupation is approximately $e^{-(E-\mu )/(k_{B}T)}$ and the probability of a valence hole is approximately $e^{(E-\mu )/(k_{B}T)}$. Substituting $x=(E-{ϵ}_c)/(k_{B}T)$ in the first integral, and the corresponding hole energy in the second, giveswhere the effective density of states values are $N_c=2[m_e{k}_{B}T/(2\pi {\hslash }^2){]}^{3/2}$ and $N_v=2[m_h{k}_{B}T/(2\pi {\hslash }^2){]}^{3/2}$. The factor $\sqrt{\pi }/2$ comes from the supplied Gaussian-weighted integral. Multiplying cancels the shared thermodynamic equilibrium chemical potential and proves the semiconductor mass-action law:For the given sample $N_D=2\times 10^{21}{\mathrm{m}}^{-3}$ and $N_A=10^{21}{\mathrm{m}}^{-3}$. With complete donor ionization and occupied acceptors, charge neutrality gives $n-p=N_D-N_A=10^{21}{\mathrm{m}}^{-3}$. At $293\mathrm{K}$ the given masses and gap giveSolving the compensated semiconductor equations, or using $n_i\ll N_D-N_A$, yieldsThe very small minority density is consistent with the wide gap and the strongly n-type balance. Also $n/N_c\simeq 2.6\times 10^{-3}\ll 1$ and $p/N_v\ll 1$, checking the nondegenerate approximation. The acceptor levels are occupied in this n-type regime; complete acceptor ionization in the compensation equation is appropriate to the model.

At low temperature dopant freeze-out reduces the conduction-electron population; a quantitative curve would require dopant binding energies and occupation factors. Once donors are ionized while intrinsic excitation is negligible, $n\simeq N_D-N_A$ gives an extrinsic plateau. At higher temperature intrinsic pairs eventually dominate, so $n\simeq p\simeq n_i(T)\propto T^{3/2}{e}^{-E_g/(2k_{B}T)}$ if the gap and masses are treated as constant. Real gap and mass changes modify that last expression.

Use a representative conduction-electron density $n\sim 10^{29}{\mathrm{m}}^{-3}$, roughly one mobile electron per atom in a dense metal, and the bare electron mass for this estimate. The plasma frequency is thenThis is an order-of-magnitude answer because no particular metal or carrier density is specified. It corresponds to an ordinary frequency $f_p\simeq 2.8\times 10^{15}\mathrm{Hz}$, excitation energy $\hslash {\omega }_p\simeq 12\mathrm{eV}$, and free-space wavelength near $1.1\times 10^{-7}\mathrm{m}$.

The physical oscillation is collective motion of the electron gas against the positive ionic background. Displacing the electrons relative to the ions produces a restoring electric field. In the simplest uniform-displacement picture, $\ddot{x}+{\omega }_p^{2}x=0$: the frequency increases with the density of mobile electric charge and decreases with carrier inertia. Quantizing this longitudinal collective oscillation gives a plasmon.

The same scale controls the ideal collisionless Drude model optical response. Its cold collisionless plasma dielectric response is ${ϵ}_r(\omega )=1-{\omega }_p^2/{\omega }^2$. Below ${\omega }_p$, the negative relative permittivity makes a transverse electromagnetic wave evanescent, producing strong reflection; above it, the ideal model admits propagation. Real-metal collisions, interband photon absorption and background polarization modify this simple optical edge. Thus ${\omega }_p$ is a collective charge-response scale, not a single-electron scattering rate or a sharp transparency threshold for every real metal.

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.

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.

Use renormalized propagators of mass $m$ and tree vertex $-ig{\mu }^{ϵ}$. The interaction counterterms have Feynman factors $-i{\delta }_{m^2}$ for a two-leg mass insertion and $-i{\mu }^{ϵ}g{\delta }_g$ for a four-leg insertion. At one loop the two-point one-particle-irreducible vertex isThe directed tadpole has $S=1$: interchanging its two ends would reverse the field labels, so the real-scalar tadpole's factor $1/2$ does not apply. It is independent of external momentum, so no wave-function counterterm is required at this order.

For the four-point vertex, label the two unstarred incoming fields by $1,2$ and the starred ones by $3,4$, with all momenta incoming. The tree contact, the three oriented bubbles and the contact counterterm giveThe first bubble has two equivalent arrows along the same direction between its vertices, giving $S=2$. In the other two channels the arrows oppose one another, and interchanging the lines does not preserve their oriented attachments, giving $S=1$. This is the oriented complex-scalar quartic graph factors rule.

The mass tadpole is order $g{m}^2$, while the four-point loop is order $g^2$. In the multiplicative coupling convention of the source,The distinction between ${\delta }_g$ and the whole additive four-leg coefficient $g{\delta }_g$ matters for the requested order counting.

To first order, invert $m_0^2=m^2[1+g/(8{\pi }^2ϵ)]$ and use $g={\lambda }_0{\mu }^{-ϵ}+O({\lambda }_0^2)$:The remainder is a perturbative statement at fixed regulator. When differentiating, hold the bare mass $m_0^2$ fixed as well as ${\lambda }_0$; otherwise the mass running would not be defined by a fixed bare theory. Since $\mu {\partial }_{\mu }{\mu }^{-ϵ}=-ϵ{\mu }^{-ϵ}$,Thus, with the source's extra factor of one half,The factor $-ϵ$ from the running dimensionless bare coupling cancels the counterterm pole. Differentiating as though $g$ were independent of $\mu$ would miss this finite result. Some conventions call $\mu d\log m^2/d\mu$ itself the mass anomalous dimension; that convention is twice the one defined here.