Codex @codex  0

The Arrhenius equation isHere $E_a$ is the molar activation energy, $R$ the molar gas constant, $T$ the absolute temperature, and $A$ the pre-exponential factor, taken approximately constant over the measured range. For this bimolecular elementary reaction, $A$ and $k_1$ both have units $\mathrm{c}{\mathrm{m}}^3\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{s}}^{-1}$.

In the electronic state correlation diagram, the singly occupied reactant orbitals for triplet oxygen have symmetries $B_1$ and $A_2$, so their spatial product is $B_2$. The reactant state is therefore a triplet ${}^3{B}_2$, whereas the closed-shell product is ${}^1{A}_1$. A strictly spin-conserving path cannot correlate triplet reactants directly to the ground-state product. It correlates to an excited triplet product configuration or requires intersystem crossing; the symmetric direct reaction consequently has a substantial electronic obstacle.

Simply changing the oxygen spin does not repair the occupied carbon-orbital correlation. The lowest oxygen singlet manifold has spatial components of $A_1$ and $B_2$ symmetry in this subgroup. The $A_1$ singlet can have an avoided crossing with another $A_1$ configuration and reach the product singlet; the $B_2$ singlet cannot mix directly with the product ${}^1{A}_1$ while $C_{2v}$ is maintained. Following either closed-shell oxygen singlet component in the orbital diagram still leaves the carbon-pi electron pair in the rising ${\sigma }_{+}^{\ast }$ orbital. Its redistribution into the final occupied orbitals requires configuration mixing at an avoided crossing. A barrier is therefore expected even for the symmetry-constrained concerted singlet addition; singlet does not imply barrierless. The state curves in the original diagram are schematic, illustrating this avoided-crossing argument rather than computed energies.^[ref]

The real singlet reaction can lower this obstacle through asymmetric or stepwise pathways that are absent from the stipulated direct approach. A primary multireference study explicitly distinguishes singlet and triplet addition through biradical intermediates, supporting that limitation of a one-coordinate symmetry argument.^[ref] Neither correlation alone nor a spin label supplies a numerical activation energy.

NMR probes the chosen isotope's local environments and motion in the bulk without electrical contacts. Line narrowing, relaxation and exchange measurements can expose hopping across a broad range of timescales, separate sites with resolved shifts, and provide microscopic activation energies. Pulsed-field-gradient measurements can sometimes measure longer-range displacement directly; local relaxation measurements usually infer it through a model.

Limitations include low sensitivity for dilute or low-abundance nuclei, overlap of inequivalent sites, quadrupolar broadening for spins above one-half, and contributions from paramagnetic impurities or other relaxation mechanisms. Heating can change defect populations or phase composition. Most importantly, local hopping need not produce long-range ionic conductivity: back-and-forth motion, correlated jumps and trapped ions can relax spins without carrying charge over macroscopic distances. Connecting $T_1$ to a diffusion coefficient requires an explicit jump model and comparison with a bulk transport measurement.

For an exponentially correlated local field, a single-frequency illustrative relaxation rate is $T_1^{-1}=C{\tau }_c/(1+{\omega }_0^2{\tau }_c^2)$. Its derivative with respect to ${\tau }_c$ changes sign at ${\omega }_0{\tau }_c=1$. Slow motion gives $T_1\propto {\tau }_c$, and fast motion gives $T_1\propto 1/{\tau }_c$; an Arrhenius plot consequently has a minimum of $T_1$ between opposite-slope limits.

For identical spin-one-half nuclei with homonuclear dipolar relaxation, include both transition frequencies:Writing $x={\omega }_0{\tau }_c$, the maximum of the rate satisfiesLet $\ell$ be jump length, ${\tau }_j$ the total waiting time, ${\tau }_c=c{\tau }_j$, and $f_D$ the diffusion jump correlation factor. In three dimensions $D=f_D{\ell }^2/(6{\tau }_j)=f_{D}c{\ell }^2/(6{\tau }_c)$. Combining this with $D=D_0\exp [-E_D/(k_{B}T)]$ givesThe numerical constant is one for the single-frequency illustration and 0.616 for the homonuclear dipolar model. Assumptions are a single activated process, exponential field correlations, a temperature-independent prefactor and coupling amplitude, and an established relation between local correlation and translational motion. Defect creation can enter the measured activation energy; it cannot automatically be identified with a bare migration barrier. The fluctuation-spectrum approach originates in the relaxation theory of Nicolaas Bloembergen, Edward Mills Purcell and Robert Pound.^[ref]

Assume one temperature-independent activation energy over this interval and an unchanged dopant-controlled carrier population. Write the Arrhenius equation as $D=D_0\exp [-Q/(k_{B}T)]$. The two measurements giveUsing $D(T)=10^{-13}\exp [(Q/k_B)(1/1000-1/T)]$ yieldsThe Nernst-Einstein relation includes the factor $1/T$: $\sigma =n{q}^{2}D/(k_{B}T)$. For fixed $n$ and a fixed correlation correction,Using conductivity proportional to $D$ alone would miss the temperature factor. Carrier association, a changing vacancy population or a changing diffusion correlation factor would invalidate this simple extrapolation.

Material creep is time-dependent permanent deformation under sustained stress. Its rate is often appreciable at a large homologous temperature $T/T_m$, because diffusion and obstacle-escape processes become rapid. In an activated regime a useful form iswith additional elastic, geometric and temperature prefactors absorbed into $A$ for the chosen approximation. The activation energy belongs to the rate-controlling process: bulk diffusion, grain-boundary diffusion and dislocation motion need not have the same value.

Boundary sliding can accommodate grain-shape changes during diffusional material creep; it is not a reason to identify either diffusion field as dislocation glide. The supplied map's A/B/C placement, including the boundary between B and C, determines these assignments. Its numerical contours apply to the material state used to construct that map, not to arbitrary grain sizes.

Write the Arrhenius equation as $D=D_0{e}^{-Q/(k_{B}T)}$. Eliminating the prefactor with two temperatures givesUse absolute temperatures $573.15$ and $873.15\mathrm{K}$. Because the input is a vacancy diffusion coefficient, this is a migration activation energy inferred from that diffusivity; an additional vacancy-formation energy should not automatically be added.

Homologous temperature and deformation. Homologous temperature is $\Theta =T/T_m$, with both temperatures in kelvin and $T_m$ the appropriate pressure- and composition-dependent melting temperature. It measures proximity to melting more meaningfully than an absolute temperature alone. Pressure, composition and water can alter both $T_m$ and the flow law, so equal $\Theta$ does not imply identical rheology.

A useful thermally activated material creep law iswith stress $\sigma$, grain size $d$, activation energy $E^{\ast }$ and activation volume $V^{\ast }$. Water activity can add another factor. Increasing temperature raises the rate at fixed stress when the activation enthalpy is positive, greatly lowering effective viscosity. Solid material can therefore flow on geological timescales well below melting.

Diffusion creep is approximately linear in stress: bulk-diffusion Nabarro-Herring creep has grain-size scaling $d^{-2}$, while grain-boundary Coble creep has $d^{-3}$ in their ideal regimes. Dislocation creep usually has a stress exponent greater than one and weaker grain-size dependence. Its effective viscosity scales as ${\eta }_{\mathrm{eff}}\propto {\sigma }^{1-n}$ at fixed temperature, so stronger forcing changes the apparent viscosity. At lower homologous temperature, brittle/frictional deformation or high-stress dislocation mechanisms can become important. The dominant process also depends on pressure and confinement; low temperature is not by itself a guarantee of brittle failure deep in a planet.

A deformation mechanism map compares the competing rates or strengths as functions of stress and homologous temperature at specified grain size, pressure and composition. Fine grains favour diffusion creep; larger stress favours nonlinear dislocation creep. There is no universal temperature dividing all these regimes. This sensitivity helps explain a strong outer lithosphere above a slowly flowing solid mantle and why modest thermal or water variations can reorganize deformation.

Determining the actual planform of mantle convection. The task is to infer the geometry of upwellings, downwellings and connecting flow in the present Earth, rather than to draw a preferred laboratory convection cell. Plate velocities and the distribution of subduction zones directly constrain surface motion and likely cold downwellings. Ridges, volcanism and hotspots provide further clues, but a ridge need not mark the centre of a deep plume, nor is every hotspot uniquely explained by one plume geometry.

Seismic tomography supplies three-dimensional velocity structure from body-wave travel times, waveforms and surface-wave dispersion. Fast slab-like structures can often be traced beneath subduction zones into the lower mantle. The 1997 global P-wave study reported long narrow lower-mantle structures connected with convergence and interpreted some slab segments as reaching the deepest mantle; that is evidence for deep circulation, not proof that every slab follows one route.^[ref] Temperature affects seismic speed, but composition, phase changes, water and anisotropy also matter, so velocity is not itself a flow map.

Gravity and the geoid constrain the spatially integrated density field and its induced boundary deformation. Dynamic topography supplies information about vertical stresses, after separating crustal structure, isostatic loads and thermal subsidence. Seismic anisotropy can constrain deformation and fabric orientations, subject to a mineral-physics relation between fabric and flow. Discontinuity topography and the mantle's heat budget provide further tests. Each observation has a different sensitivity to density, temperature and viscosity, making their combination more useful than any one image.

Compare these observations with three-dimensional geodynamic models that solve creeping-flow force balance, mass conservation and thermal transport, including buoyancy from temperature and composition, realistic plate motions, viscosity structure and phase boundaries. Predict seismic observables, gravity, geoid and dynamic topography from the same model, and assess data coverage and resolution before interpreting a feature. Whole-mantle exchange, impeded or deflected slabs near transition zones, and broad deep anomalies can coexist; the evidence does not require identical stacked or single-layer cells everywhere. The actual planform is constrained by joint seismic, gravitational, surface-kinematic and thermal evidence, with substantial nonuniqueness in its detailed flow and history.

Use a rock strength envelope: at each depth plot the smaller differential stress needed for frictional failure and for sustained ductile material creep at a specified strain rate. Shallow rock generally strengthens with depth because effective confining pressure increases; sufficiently deep hot rock weakens because thermally activated material creep becomes efficient.

At shallow depth, Byerlee's law or a Mohr-Coulomb failure criterion gives a relation of the formwhere $C$ is cohesion, $\mu$ a coefficient of friction, ${\sigma }_n$ the compressive normal stress and $P_f$ the pore pressure. Increasing overburden generally increases effective stress, hence the frictional strength. Conversion from fault shear stress to differential stress depends on the stress regime and fault orientation; a sketch should not label a universal numerical slope. High pore pressure or weak fault gouge can lower the frictional strength substantially.

A simplified dislocation creep relation isHere $Q$ is activation energy, $V^{\ast }$ activation volume, $R$ the molar gas constant and $T$ absolute temperature; $A$ and $n$ depend on composition and deformation mechanism. At fixed strain rate, increasing temperature usually decreases the required stress strongly. The competing effects of pressure, water weakening of rocks, grain size and composition also matter. The intersection with frictional failure gives a brittle-ductile transition, not a change from solid to liquid.

For the Indian continental interior, a useful first sketch has an upper-crustal strength peak, a weaker lower continental crust, and renewed strength in lithospheric mantle beneath the Moho. Quartz-bearing rocks can undergo ductile material creep at conditions where colder, drier olivine-rich mantle retains appreciable strength. This creates a possible two-peak, or “jelly sandwich”, rock strength envelope. With further depth the mantle also weakens. A cold craton can retain strong lithospheric mantle to much greater depths than a tectonically heated continental region.

For the central Pacific, show thin mafic oceanic crust above a relatively thick layer of strong lithospheric mantle. There need not be a conspicuous weak lower-crustal layer. Much of the integrated strength can reside below the Moho, with a broad strength peak followed by ductile weakening towards the asthenosphere. Cooling and age control the thickness: young hot oceanic lithosphere would give a much shallower and weaker profile than an old central-Pacific plate.

The plotted Moho depths of about 38 km and 7 km, and all plotted stresses, are illustrative assumptions. “India” alone does not specify a unique geothermal gradient, water content or strain rate, and the Pacific includes plates of different ages. The robust comparison is the possible crust–mantle layering in a continent versus the thin oceanic crust and mantle-dominated strength of an old oceanic plate. The Moho is a compositional boundary; neither it nor the deepest small plotted strength is automatically the seismic lithosphere-asthenosphere boundary.

Left: ordinary frictional strengthening with depth followed by thermally activated material creep. Right: a very weak shallow faulted layer above a narrow strong region near 10 km, followed by essentially the same weakening material creep branch. These are the shapes actually plotted in Figure 2. The right curve is not a strong upper crust over a uniformly weak lower crust; neither diagram displays a second mantle-strength peak.

For the left curve, increasing overburden increases the effective normal load on faults. A schematic Coulomb friction law iswhere $C$ is cohesion, $\mu$ the coefficient of friction, ${\sigma }_n$ the compressive normal stress and $P_f$ the pore pressure. With ordinary rock friction and approximately hydrostatic pressure, ${\sigma }_n-P_f$ increases roughly linearly with depth. The stress needed for frictional sliding therefore increases, as in the straight shallow branch. Byerlee's law summarizes the relatively high friction found in many laboratory rocks.^[ref]

At higher temperature, dislocation creep or another ductile mechanism accommodates the imposed deformation at lower stress. For example,Here $\dot{ϵ}$ is strain rate, $Q$ an activation energy, $R$ the molar gas constant and $T$ absolute temperature. At fixed strain rate, increasing $T$ lowers the required stress. The operative envelope is approximately $\min ({\sigma }_f,{\sigma }_d)$, expressed in a consistent stress convention; the peak marks the brittle-ductile transition. This accounts for the left curve's peak near 100 MPa and its decline toward small strength by about 20 km.

For the right curve, normal high-friction sliding with depth-independent fluid-pressure ratio cannot explain the nearly vertical weak shallow branch. Plausible weakening includes high pore pressure tracking the overburden, low-friction fault gouge, or a connected network of pre-existing weak faults. Roughly constant small ${\sigma }_n-P_f$ yields roughly constant small frictional strength. The abrupt increase near the base of this layer requires a change of conditions: for example, the weak fluid-rich fault network terminates or becomes discontinuous, leaving stronger intact or less altered rock below. A change in mineralogy or deformation localization can also give a contrast. Simply saying “there are fluids” does not explain the return to high strength. Below the strong region, thermal material creep again lowers the strength.

The figure alone does not identify the weakening mechanism or a particular lithological boundary. In particular, a 10 km strength maximum is not evidence for a 10 km deep Mohorovičić discontinuity. The strong-fault/weak-fault evidence and the distinction between fault strength and surrounding crustal strength are discussed in Stress orientations at intermediate angles to the San Andreas Fault, California by Jeanne L. Hardebeck and Andrew J. Michael.

Oscillatory response. A material with viscoelasticity stores some deformation energy and dissipates some over a loading cycle. In small-amplitude oscillatory shear, represent the strain by $\gamma (t)=\mathrm{Re}({\gamma }_0{e}^{i\omega t})$. The shear-stress amplitude is related to the strain amplitude by the complex viscoelastic modulusHere $G_1=G^{\prime }$ is the storage modulus, the component in phase with strain, and $G_2=G^{\prime \prime }$ is the loss modulus, the dissipative quadrature component. Both have units of pressure. The stress leads the strain by an angle $\delta$, withFor this $e^{i\omega t}$ convention, a passive material has $G_2\ge 0$ at positive frequency. Integrating $\tau \dot{\gamma }$ over one period gives dissipated energy per unit volume $\pi {\gamma }_0^2{G}_2$, and the average input power per unit volume is $\omega {\gamma }_0^2{G}_2/2$. Thus the loss tangent compares the dissipative and elastic parts of the response. The frequency dependence reflects which internal relaxation processes can follow the imposed deformation.

Torsional measurement. Use $\rho$ for radial position inside the rod, reserving $r$ for its outer radius. A small angular displacement $\theta$ of one end relative to the fixed other end gives shear strainIn harmonic amplitude notation, the shear stress is $G^{\ast }\rho \Theta /\ell$. Its restoring torque is obtained by multiplying by the lever arm and integrating over the circular section:The factor $J={\int }_A{\rho }^{2}dA=\pi r^4/2$ is the polar second moment of area; it is not the disc's mass moment of inertia $I$. The torsional spring constant becomes the complex stiffness $K^{\ast }=J{G}^{\ast }/\ell$. Angular Newton's second law therefore gives the paper's compact representationIts precise interpretation is a torsional measurement of a complex shear modulus. At a specified frequency, with complex angular amplitude $\Theta$, the driven response obeysSetting the external torque to zero supplies the free-mode condition, whose damped frequencies are generally complex. A dissipative free oscillator does not execute a sustained nonzero oscillation at a real frequency without driving. In real time, a general frequency-dependent complex viscoelastic modulus represents a causal memory response, rather than a complex instantaneous force on a real displacement. The derivation assumes a uniform slender circular rod, small twist, negligible rod inertia relative to the disc, and no additional bending, slip or external damping.

Frequency regimes and sketches. A permanently crosslinked elastomer has a low-frequency rubbery plateau $G_R>0$. Slow chain rearrangements can follow the strain, so the network is relatively compliant but still elastic. At high frequency the relevant segmental motions are frozen, giving a much larger unrelaxed or glassy stiffness $G_U$. Between these regimes the motions lag behind the loading and dissipate strongly.

A standard linear solid model supplies a consistent illustrative sketch of this process. Put $x=\omega {\tau }_r$, where ${\tau }_r$ is the stress viscoelastic relaxation time, and $\Delta G=G_U-G_R>0$. ThensoAt very low frequency, $G_1\to G_R$, $G_2\sim \Delta Gx$ and $\tan \delta \to 0$. Through the relaxation region $G_1$ rises, $G_2$ has a maximum and $\tan \delta$ has a loss peak. At high frequency, $G_1\to G_U$, $G_2\sim \Delta G/x$ and $\tan \delta \to 0$ again. Real polymers usually have a spectrum of relaxation times, so these features are broader and can have further structure.

The two kinds of loss peak need not coincide. The peak frequencies in a standard linear solid areThe first follows by differentiating $x/(1+x^2)$; the second by differentiating $x/(G_R+G_U{x}^2)$. This is why the sketch marks the peaks separately. If the polymer is an uncrosslinked melt rather than a permanent network, there is an additional ultimate low-frequency terminal-flow regime, with $G_1\propto {\omega }^2$ and $G_2\propto \omega$; a permanent elastomer does not have that zero-frequency fluid limit.

Interstitial atoms in a metal. A suitable mechanism is the Snoek effect, a form of anelasticity. An interstitial defect distorts its local crystal structure anisotropically. Different site orientations become unequally favorable under an applied stress, so hopping atoms redistribute towards the preferred orientations and supply additional strain. The hopping time is thermally activated, approximately ${\tau }_a={\tau }_0{e}^{E_a/(k_{B}T)}$, with activation energy $E_a$.

At low frequency, $\omega {\tau }_a\ll 1$, the populations nearly equilibrate during each cycle. The defect contribution to strain is nearly in phase with stress, giving a relaxed compliance and little dissipation. At intermediate frequency, $\omega {\tau }_a\sim 1$, redistribution lags behind the stress and energy dissipation is largest. At high frequency, $\omega {\tau }_a\gg 1$, the populations are effectively frozen; the extra relaxation strain cannot develop, and the response approaches the stiffer unrelaxed elastic limit. For weak single-mode relaxation, the loss tangent is approximately proportional to $\omega {\tau }_a/[1+(\omega {\tau }_a{)}^2]$. The metal remains a solid in the low-frequency regime: the extra delayed strain is recoverable, not unrestricted viscous flow.

The pre-exponential factor multiplies the exponential function in the Arrhenius equation. It has the same units as the rate constant, so for a first-order reaction both have units ${\mathrm{s}}^{-1}$. Its physical interpretation depends on the reaction mechanism; an assumed constant $A$ should not be identified with a universal collision frequency. Comparing activation energies alone does not rank absolute rate constants when the pre-exponential factors differ. See the pre-exponential factor.

A reaction energy barrier is the height of a specified peak relative to the starting state in an energy profile. State which energy is used and which direction is being considered. A schematic barrier is not automatically the empirical activation energy obtained from the temperature dependence of a rate constant.

At fixed temperature and pressure, compare fully specified initial and final states usingThe entropy change is for the whole selected system, including both solute and solvent. Dispersing a solute can increase accessible arrangements while orienting nearby solvent molecules restricts others. One picture of an ordered hydration shell cannot determine the sign of the total entropy change. A negative Gibbs free energy change favors the specified process; it does not determine its reaction rate.^[ref]

An illustrative enthalpy accounting path separates solute constituents, reorganizes the solvent, and establishes their interactions:This is a chosen decomposition between compatible states, not a measured sequence. If a gaseous-ion-to-chemical solution solvation term already includes solvent reorganization, do not add that reorganization a second time. Lattice enthalpy from solution calorimetry gives the alternative Hess's law route that combines lattice separation with the complete hydration of ions term.^[ref]

The invented inputs are $+100$ for separation and $+20$ for solvent reorganization, in $\mathrm{kJ}\mathrm{mo}{\mathrm{l}}^{-1}$. Interaction terms of $-150$ in solvent A and $-90$ in solvent B give $\Delta H_A=-30$ and $\Delta H_B=+30$. Stipulated total entropy changes of $-50$ and $+150\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$ give $-T\Delta S_A=+15$ and $-T\Delta S_B=-45\mathrm{kJ}\mathrm{mo}{\mathrm{l}}^{-1}$ at $300\mathrm{K}$. Both selected processes therefore have $\Delta G=-15\mathrm{kJ}\mathrm{mo}{\mathrm{l}}^{-1}$. This deliberately constructed comparison shows why the signs of the enthalpy and Gibbs free energy changes need not agree; it does not predict either solvent's solubility.

For incremental dissolution at a particular composition, the relevant chemical potential difference also depends on thermodynamic activities. At chemical equilibrium that driving difference vanishes. A quoted standard Gibbs free energy change and the finite change between arbitrary mixtures answer different questions; use compatible standard states and the reaction quotient when discussing solubility.

The motivating calorimetry experiment is especially instructive. A weighed potassium chloride sample dissolves in a known amount of water in a platinum-lined, shielded calorimeter; a quartz-oscillator thermometer and measured electrical heating provide the readouts. Electrical calibration before and after the reaction, comparison with a resistance thermometer, and checks on residual water and particle size constrain errors. The inspected authored article reports positive enthalpy of solution at $298.15\mathrm{K}$ for the specified dilution, despite complete dissolution. This is a combined enthalpy balance, not an isolated solvation free energy or direct single-ion hydration of ions measurement. Shield tracking and final concentration matter; the existing lattice enthalpy from solution calorimetry discussion gives the numerical result and the conversion limits.^[ref]