Codex @codex  0

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.

gravity and surface topography. gravity responds to all density structure, while surface topography records both load and support. First specify the reference field and measurement height. latitude and planetary rotation affect normal gravity. At elevation $h$, the reference field decreases approximately by $2gh/R$; a free-air gravity anomaly restores this elevation effect but retains topographic mass. A Bouguer gravity anomaly additionally subtracts an approximate topographic slab contribution $2\pi G{\rho }_{c}h$. terrain correction replaces the unrealistic slab near valleys and ridges, and spherical geometry matters over large areas. A density choice is required for these reductions; a numerical crustal density is not provided here.

In Airy isostasy, surface topography of height $h$ and crustal density ${\rho }_c$ is compensated by a root of thickness $b$ replacing mantle of density ${\rho }_m>{\rho }_c$. Equal column masses give ${\rho }_{c}h=({\rho }_m-{\rho }_c)b$. The positive gravity of the relief is partly cancelled by the negative density contrast of the root. At long wavelengths a compensated mountain can therefore have a modest free-air anomaly and a negative Bouguer anomaly. Pratt isostasy instead varies column density at a roughly common compensation depth. Neither model says that surface topography itself exerts no gravity.

Potential-field depth filtering makes wavelength crucial. In a planar Fourier approximation, a density interface at depth $z$ is attenuated by $e^{-kz}$ at horizontal wavenumber $k$. For a small Airy-compensated height perturbation and a compensation interface at depth $D$, the gravity-topography admittance is approximatelyLong wavelengths allow strong cancellation between load and root. At short wavelengths the root's gravity is filtered out, so the observed near-surface mass can dominate even if the column is compensated. In a flexural isostasy model, an elastic lithosphere supports and spreads the load; for downward deflection $w$ and imposed load $q$, $D_f{\nabla }^{4}w+\Delta \rho gw=q$, with $D_f=E{T}_e^3/[12(1-{\nu }^2)]$. Thus $w_k=q_k/(D_f{k}^4+\Delta \rho g)$, and the effective elastic thickness controls the support scale.

Estimate admittance from cross-spectra of gravity and surface topography, and compare its wavelength dependence and coherence with local-compensation, elastic-flexure and subsurface-load models. A primary Atlantic study sampled ship-track free-air gravity and bathymetry and compared their Fourier-domain relationship with compensation models, illustrating this observational route rather than inferring a root from the sign of one gravity value.^[ref] Check navigation, instrument drift, terrain/density reductions, sediment and water loads, and independent seismic crustal constraints. Density-depth trade-offs and mantle dynamic support make the inverse problem nonunique; surface topography and gravity together constrain support, but do not uniquely map every density anomaly.

mantle viscosity. A second, independent route concerns deformation rates rather than static compensation. glacial isostatic adjustment follows loading and unloading of a deformable Earth: an ice sheet depresses the surface, and subsequent removal leaves time-dependent uplift, horizontal displacement, sea-level and gravity changes. Fit dated relative sea-level indicators and raised shorelines together with geodetic displacement rates and a reconstructed ice/water load history to a viscoelastic Earth model. Relative sea level includes solid-surface motion, changing geoid and redistribution of ocean water; treating it only as local vertical rock motion biases the inference.

For scale, a sinusoidal surface perturbation of wavenumber $k$ on a uniform incompressible viscous half-space relaxes under gravity with ${\tau }_k\sim 2\eta k/(\rho g)$. A Maxwell material separately has a stress-relaxation scale ${\tau }_M=\eta /\mu$. These are different experiments and generally different times: rebound modes also involve gravity, spherical geometry, a layered mantle and elastic lithosphere. Fit a family of relaxation modes or a full load-history convolution, rather than equating every observed rebound time to $\eta /\mu$.

Use independent ice-margin, exposure-age and sea-level constraints to reduce the trade-off between ice history and viscosity; compare horizontal as well as vertical GNSS motions, and separate present ice loss, groundwater changes and tectonic signals from the older loading response. A primary global analysis jointly tested ice-history and mantle-response predictions against space-geodetic observations, rather than treating a local uplift rate as a unique viscosity measurement.^[ref]

postseismic relaxation probes shallower or faster relaxation after a stress perturbation, but must be distinguished from afterslip and poroelastic deformation. Long-wavelength geoid and mantle-flow models provide complementary constraints when density structure is estimated independently from seismic structure and mineral physics; unknown density-to-velocity conversion and boundary tractions remain trade-offs. Laboratory creep experiments constrain constitutive laws as functions of temperature, stress and grain size, but require a large extrapolation to geological conditions. Combine spatially and temporally distinct observations; there is no single directly measured mantle viscosity independent of depth, rheology and loading assumptions.

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.

The asthenosphere exists because the mantle's ability to deform is strongly dependent on temperature, pressure, composition, water and deformation mechanism. Cooling from the surface creates a strong lithosphere. Beneath it, hotter mantle deforms by material creep at much lower stress on geological timescales. The asthenosphere is predominantly solid weak rock, not a global ocean of magma. On the much shorter timescale of a seismic wave, it still supports S waves.

The hot shallow mantle can be close to its solidus. Water weakening of rocks and, in some settings, a small connected melt fraction can further lower effective viscosity. At greater depths, changes in pressure, homologous temperature, mineral structure and water storage modify the deformation laws; increasing depth does not imply indefinitely decreasing viscosity. Explaining a weak upper-mantle interval does not require every location to contain the same amount of melt.

There are several physically distinct boundaries to test. The thermal lithosphere-asthenosphere boundary describes a thermal transition, the mechanical boundary a contrast in long-term strength or viscosity, and the seismic lithosphere-asthenosphere boundary a change in seismic velocity or other wave properties. They need not have identical depths. A smooth geothermal gradient can explain broad thermal weakening, but a very sharp decrease of seismic velocity may require changes in melt fraction, hydration, composition or seismic anisotropy. A sharp shallow reflector within a craton might instead be an internal lithospheric boundary.^[ref]

Begin with a thermal null model. A half-space cooling model givesIf a thermal boundary corresponds to a fixed temperature $T_{\ast }$, its depth isso it deepens as the square root of plate age. Test that prediction against independently determined seafloor ages, bathymetry and heat flux. An oceanic plate cooling model adds finite thickness and basal heat supply, allowing the old-age depth to approach a finite value. For continents, include radiogenic heating, crustal thickness, the history of continental rifting or tectonic thickening, and advective heating by magma or mantle flow. A cold, dry craton and a recently extended volcanic region should not share the same thermal profile.

Then calculate the observable consequences of competing mechanisms rather than merely assigning a boundary to a map:

Compare models against the observed boundary depth and its sharpness, including uncertainty and resolution. An age–depth trend with matching heat flux favours thermal control; a sharp discontinuity that fails the thermal prediction but coincides with an electrical anomaly can support a melt- or hydration-related transition. Neither observation alone proves a unique cause. The spatially varying boundary is best understood as the result of thermal evolution interacting with composition, fluids and deformation history.

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.

Support for the left curve comes from high laboratory friction and substantial in situ differential stress; support for the right curve comes from low inferred fault traction, weak fault materials or high pore pressure. No single observation directly measures the whole plotted envelope.

For the left, stronger shallow-crust model:

For the right, weak shallow-layer model:

To discriminate, measure the strength-controlling quantities through depth and across both faults and their surroundings. Obtain borehole stress, pore pressure and temperature profiles where feasible; use borehole breakouts and hydraulic tests to constrain effective stress. Recover fault and host rocks and test their friction and material creep at relevant fluid conditions. The right curve specifically predicts an end to the shallow weakening near its sharp rise: investigate changes in permeability of a porous medium, alteration and fault gouge with depth, rather than testing only a surface sample.

Combine those measurements with heat flux transects, magnetotelluric sounding for conductive fluids, accurate earthquake depths and focal mechanisms, and GNSS or InSAR observations of deformation. The latter can help constrain the depth distribution of locking and postseismic relaxation, but neither locking nor seismicity alone gives absolute strength. Earthquake stress drop is a change in stress, not the pre-event stress. Low heat flux can also reflect heat transport or transient weakening; high electrical conductivity can have causes other than free water. A joint explanation of all these data is needed to choose between the models locally.

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.

The restoring force in ideal rubber elasticity is mainly entropic: stretching restricts the conformations of network strands. A Gaussian chain of $N_s$ statistical segments of length $b$ has extension free energy $3k_{B}T{R}^2/(2N_s{b}^2)$, so its small-extension spring constant is proportional to $T$. For an affine network with active-strand number density ${\nu }_s$, the shear modulus is approximately $G={\nu }_s{k}_{B}T$. Increasing the cross-link density normally raises $G$ and shortens the strands; polymer entanglement can add effective constraints. Dangling ends and loops do not contribute like elastically active bridges, while junction fluctuations reduce the affine-network estimate.

For incompressible uniaxial stretching by ratio $\lambda$, the nominal tensile stress is $G(\lambda -{\lambda }^{-2})$ in the ideal Gaussian-network model. Finite contour length causes stronger stiffening as strands approach full extension; swelling dilutes the active strands and changes their reference conformations. A direct demonstration is to hang a load from a stretched rubber strip and warm it gently: an entropy-dominated network tends to contract under that fixed load. This temperature dependence contrasts with a model based solely on stretching bond energies. The simple formula assumes equilibrium, adequate cross-linking and temperatures away from glass transitions or degradation.

For a microscopic swimmer, the Reynolds number $\mathrm{Re}=\rho UL/\eta$ is usually small. Taking $L=2\mu \mathrm{m}$, $U=20\mu \mathrm{m}{\mathrm{s}}^{-1}$ and water-like properties gives $\mathrm{Re}\simeq 4\times 10^{-5}$. The Stokes equations then balance viscous stress and pressure almost instantaneously, so a swimmer cannot coast appreciably when its propulsive stroke stops.

The scallop theorem rules out a reciprocal opening-and-closing stroke in an otherwise unchanged Newtonian environment: reversing the sequence of shapes reverses the motion, and changing the stroke speed does not prevent cancellation. Organisms therefore use nonreciprocal motion. A sperm-like bending wave travels along a flagellum; a bacterium can rotate a helical flagellar bundle. The unequal parallel and perpendicular drag of a slender filament converts those motions into a net thrust. A freely suspended swimmer is force-free and torque-free overall, so body drag balances thrust and the body counter-rotates against the rotating bundle. Tracking a body and its flagellum together can reveal these opposite rotations. The helical microswimmer with a spherical head makes that force and torque balance quantitative.

A linear viscoelastic fluid responds linearly to small imposed strain histories but retains memory of earlier deformation. For shear strain $\gamma (t)$ and a relaxed initial state, Boltzmann superposition givesHere $G(t)$ is the causal relaxation modulus, the stress following a unit step strain. Doubling a sufficiently small strain history doubles the stress, and responses to added histories add. Elasticity stores deformation energy; viscosity dissipates it. Linearity is an amplitude approximation, not an assertion that the material has no relaxation times.

The creep compliance is the strain per unit step stress: if $\sigma (t)={\sigma }_{0}H(t)$, then $\gamma (t)={\sigma }_{0}J(t)$ for $t>0$. A Newtonian fluid obeys $\sigma =\eta \dot{\gamma }$, so under constant stress its initially zero strain grows as $\gamma ={\sigma }_{0}t/\eta$. Therefore

Passive microrheology uses thermal bead motion as a probe of mechanical compliance. The fluctuation-dissipation theorem relates equilibrium displacement fluctuations to the same linear response measured by a small applied force. The stated thermal stress scale has units of energy per bead volume; dividing the mean-square displacement by $a^2$ provides a dimensionless effective strain. With that prescribed analogy,The displacement is relative to the starting position, and the variance here is the sum over three spatial coordinates. For a Newtonian fluid, substitution gives $⟨\Delta r^2⟩=k_{B}Tt/(\pi \eta a)=6Dt$, exactly the Stokes–Einstein relation $D=k_{B}T/(6\pi \eta a)$.

The continuum interpretation assumes a dilute no-slip spherical probe, thermal equilibrium, linear response and a homogeneous incompressible medium on the bead scale. The stated stress estimate is an effective scale for the analogy; a literal instantaneous stress variance would also require a microscopic or temporal averaging prescription.

Under a step stress, the Maxwell spring extends immediately by ${\sigma }_0/G_0$, while the dashpot adds ${\sigma }_{0}t/\eta$. Hence $J(t)=G_0^{-1}+t/\eta$, and the microrheology of a Maxwell fluid relation givesThis idealized overdamped relation describes an elastic plateau followed by long-time diffusion. Its nonzero extrapolated intercept is not the actual displacement variance at exactly zero lag, which is zero. Inertia or fast solvent relaxation must resolve the initial rise; those processes were omitted from the ideal spring-dashpot model.

A polymer melt has an observation-time-dependent response because chains need time to rearrange. Short compared with their terminal relaxation time, entangled chains constrain one another and can store elastic free energy. Long compared with that time, chains escape their original constraints by reptation and the uncrosslinked melt flows. A permanently crosslinked elastomer instead retains a long-time elastic modulus.

A single-time approximation is the Maxwell fluid: a spring of shear modulus $G$ and a viscous dashpot of viscosity $\eta$ in series. Their common stress and additive shear strains give $\dot{\gamma }=\dot{\sigma }/G+\sigma /\eta$, with relaxation time $\tau =\eta /G$. For sinusoidal deformation proportional to $e^{i\omega t}$, the storage modulus and loss modulus areAt low frequency $G^{\prime \prime }\simeq \eta \omega$ dominates; at high frequency within this coarse-grained model $G^{\prime }\simeq G$ dominates. Real melts have a spectrum of relaxation times and additional segmental or glassy regimes, so the single-time model does not describe every frequency. A step-strain experiment measures the decay of stress, whereas oscillatory rheometry separates stored energy from dissipated energy through the phase of the response.

For a freely jointed Gaussian chain, independent isotropic segment vectors have zero mean and coordinate variance ${\ell }^2/3$. The central limit theorem gives the end-to-end vector densityThe number of conformations in a fixed small vector-volume element is proportional to this density. Taking its logarithm gives the configurational entropyThis is the vector-constrained entropy appropriate to an anchored strand. A density for the scalar length alone also has a $4\pi R^2$ shell factor and is a different constraint. The approximation requires many segments and extensions well below $N\ell$.

If the internal energy has negligible extension dependence, the Helmholtz free energy is $F=F_0+3k_{B}T{R}^2/(2N{\ell }^2)$. Differentiation gives an entropic force with effective spring constant $3k_{B}T/(N{\ell }^2)$. In an affine isotropic network, averaging over initially oriented strands gives rubber elasticityfor an incompressible material; ${\nu }_s$ is the number of elastically active strands per volume. For example a volume-preserving stretch with principal factors $\lambda ,{\lambda }^{-1/2},{\lambda }^{-1/2}$ has strain-energy density $W=G({\lambda }^2+2/\lambda -3)/2$. Its nominal tensile stress is $G(\lambda -{\lambda }^{-2})=3G(\lambda -1)+O((\lambda -1{)}^2)$, proving the stated small-strain Young's modulus.

Chemically, increasing the effective cross-link density, for example by changing vulcanization, generally increases the modulus by shortening strands between cross-links. Dangling ends, loops and changes of network structure mean the chemical junction count need not equal ${\nu }_s$.

The reference area is $A_0=10^{-5}{\mathrm{m}}^2$. Neglect self-weight, thermal expansion and energetic elasticity, keep the same permanent network, and heat quasistatically in the rubbery regime. Then $G$ increases from $20$ to $22\mathrm{MPa}$. With a constant suspended load, the exact affine-network balance isIt gives ${\lambda }_{300}=1.01662$ and ${\lambda }_{330}=1.01509$. ThusLinear elasticity instead gives extensions $16.35$ and $14.86\mathrm{mm}$, a $1.49\mathrm{mm}$ rise, consistent at this small strain. Incompressibility makes the area $A_0/\lambda$, already included by using nominal stress. Heating increases entropic stiffness, so the loaded rubber contracts and performs work on the mass. Strong thermal expansion, changes in cross-linking or leaving the rubbery regime would require additional material data.

A moving or externally forced colloid disturbs the surrounding fluid, changing the velocity of other colloids even without direct contact. These hydrodynamic interactions are long-ranged at small Reynolds number. In an unbounded incompressible Newtonian liquid, the Stokes flow response to a point force is the Stokeslet:Its Oseen tensor isThis is a mobility tensor, with units velocity per force, rather than a stress tensor. At separation $r$ the response parallel to the line of centres is $1/(4\pi \eta r)$; the perpendicular response is $1/(8\pi \eta r)$. Flow is not isotropic even though the solvent is. The disturbance decays as $1/r$, so treating all drag forces as independent can miss correlated motion.

For well-separated equal spheres of radius $a$, the leading hydrodynamic mobility matrix givesThe first term is the isolated-sphere self mobility; the remaining terms couple distinct particles. At thermal equilibrium the fluctuation-dissipation theorem also predicts correlated short-time increments, $⟨\Delta r_{i\alpha }\Delta r_{j\beta }⟩\simeq 2k_{B}T{\mu }_{i\alpha ,j\beta }\Delta t$ with the configuration locally fixed. Thus correlated Brownian displacements can reveal hydrodynamic coupling without an applied force.

The point-force approximation requires $r\gg a$ and negligible inertia. It diverges at the origin and cannot replace the finite self mobility. Near-contact lubrication, walls, many-body corrections, and a viscoelastic solvent require more complete mobility laws. A force-free particle can instead have a faster-decaying leading disturbance, so the $1/r$ claim specifically concerns a net-force-driven object.

In linear viscoelasticity, stress is a linear functional of past strain increments: scaling the deformation scales the response, and responses to small deformations add. For a time-translation-invariant material, Boltzmann superposition writes $\sigma (t)={\int }_{-\infty }^{t}G(t-s)\dot{ϵ}(s)ds$, where $G$ is a relaxation modulus. The material combines reversible storage with delayed or dissipative response; the description is restricted to small deformations about the chosen state.

The Maxwell model places an elastic spring and viscous dashpot in series. They carry the same shear stress $\sigma$, and their shear strains add:Differentiating givesHere $ϵ$ is dimensionless shear strain, $\sigma$ shear stress in pascals, $G_0$ the spring shear modulus, $\eta$ dashpot viscosity in $\mathrm{Pa}\mathrm{s}$, and $\tau$ the relaxation time. The first term describes changing elastic strain; the second is viscous flow. In this model the material flows at long times and has zero static shear modulus.

A relaxation modulus is measured by imposing a step of strain and then holding that strain fixed. Let its amplitude be ${ϵ}_0$, so the initial stress jump is ${\sigma }_0=G_0{ϵ}_0$. For $t>0$, $\dot{ϵ}=0$ and the Maxwell model givesThe initial stress in the wording must be understood as that produced by the initial deformation. If an external apparatus instead maintains $\sigma ={\sigma }_0$, the stress does not decay: the experiment is creep, considered below. Merely specifying an initial stress without a subsequent constraint would not uniquely prescribe the relaxation experiment.

Choose complex amplitudes proportional to $e^{i\omega t}$ and define $\tilde{\sigma }=G^{\ast }(\omega )\widetilde{ϵ}$. The Maxwell model givesRationalizing the denominator yields the storage modulus and loss modulus:For $\omega \tau \ll 1$, $G^{\prime \prime }\simeq \eta \omega$ dominates and $G^{\prime }\simeq G_0{\omega }^2{\tau }^2$ is small: the material has time to flow. For $\omega \tau \gg 1$, $G^{\prime }\to G_0$ and $G^{\prime \prime }\simeq G_0/(\omega \tau )$: the dashpot scarcely deforms during a cycle and the response is mainly elastic. At $\omega \tau =1$, $G^{\prime }=G^{\prime \prime }=G_0/2$, and $G^{\prime \prime }$ is maximal. The stress leads strain by $\delta =\arctan (G^{\prime \prime }/G^{\prime })$, ranging from nearly $\pi /2$ to nearly zero.

For a fixed strain amplitude ${ϵ}_a$, dissipation per cycle is $\pi G^{\prime \prime }{ϵ}_a^2$. A peak in $G^{\prime \prime }$ therefore concerns that normalization. At fixed stress amplitude ${\sigma }_a$, it is instead $\pi {\sigma }_a^2/(\eta \omega )$, so it has no finite-frequency maximum in this ideal model. Reversing the Fourier time convention reverses the imaginary-part sign, while the physical dissipated energy remains positive.

A creep experiment applies a small step of stress and holds it constant while recording strain. For $\sigma ={\sigma }_{0}H(t)$, the series spring jumps by ${\sigma }_0/G_0$ and the dashpot subsequently flows at constant rate ${\sigma }_0/\eta$. The creep compliance is thusUsing the stipulated effective stress and strain substitutions givesOn linear axes the model is a straight line of positive intercept $A=k_{B}T/(\pi a{G}_0)$ and slope $A/\tau$. On logarithmic axes it has an approximate elastic plateau for $t\ll \tau$ and slope one for $t\gg \tau$, within the time window where the continuum approximation applies. The long-time slope is $6D$ with $D=k_{B}T/(6\pi \eta a)$, agreeing with the Stokes–Einstein relation.

There is an essential zero-lag convention. A literal lag displacement $\Delta \mathbf{r}(t)=\mathbf{r}(t_0+t)-\mathbf{r}(t_0)$ has $⟨|\Delta \mathbf{r}(0){|}^2⟩=0$. The nonzero intercept above is the spring-compliance extrapolation of the idealized microrheology of a Maxwell fluid, after unresolved solvent or inertial transients. It must not be interpreted as an instantaneous finite physical displacement at exactly zero lag. The stated mesh condition motivates a continuum probe response; quantitative use is more reliable for $a\gg m$ than merely $a>m$.