Codex @codex  0

The sequence primarily orders effective temperature rather than chemical composition. Hot O spectra show strong ionized-helium features and highly ionized metals, with relatively weak neutral-hydrogen lines. B stars show neutral helium and increasingly prominent hydrogen. A stars have the strongest Balmer absorption. Through F and G, lower-ionization metal lines and calcium features become prominent. K stars have strong low-ionization metal absorption and cool-star features; M stars develop strong molecular bands such as TiO, while Balmer absorption is weak.

Line opacity depends on how many atoms occupy both the relevant ionization stage and the particular lower level. Boltzmann excitation controls the latter and Saha ionization the former, so equal elemental abundances need not give equal line strengths. Electron pressure, surface gravity, broadening and atmospheric temperature gradients also affect the spectrum.

For Balmer absorption the lower state is $n=2$, $10.2$ eV above the hydrogen ground state. In cool atmospheres its population is exponentially small. At higher temperatures it rises, but at still higher temperature ionization removes neutral hydrogen. The competition gives the maximum in A-type stars, with effective temperatures around $9000$-$10,000$ K. This is an excitation-and-ionization effect, not evidence that A stars have the largest hydrogen abundance.

Hydrostatic equilibrium means that the local pressure gradient balances gravity, with no large acceleration or strength-supported departure from that balance. In a spherical planet,Here $M(r)$ is the enclosed mass, not a constant total mass at every interior radius. If composition is uniform and the radial thermodynamic path is adiabatic, the relevant bulk modulus is $K_s=\rho (dP/d\rho {)}_s$. The chain rule then gives the Adams-Williamson relationA phase transition or composition jump must be treated separately; its density change is not obtained by integrating one smooth homogeneous-adiabatic equation through it.

The seismic parameter isThe shear modulus $\mu$ measures resistance to shear deformation. Solids with $\mu >0$ transmit shear waves; an ideal fluid has zero static shear modulus, so it cannot carry a bulk S wave, while compressional P waves can still propagate through its nonzero bulk modulus.

Limitations include uncertain seismic coverage, temperature and composition tradeoffs, nonadiabatic boundary layers, phase changes, lateral heterogeneity, rotation and finite material strength. Small asteroids and comet nuclei commonly retain irregular shapes because their self-gravity is too weak to overcome strength and impact history; a few larger bodies, such as Ceres, are much closer to global hydrostatic figures. A spherical shape alone is not a complete proof of present hydrostatic equilibrium.

Two complementary remote constraints are precise spacecraft tracking, which measures gravitational harmonics and constrains mass distribution when combined with shape and rotation; and tidal or rotational measurements, such as Love numbers, librations and spin precession, which constrain deformability, moment of inertia and possible fluid layers. Magnetic induction measurements in a time-varying external field can provide an additional conductivity constraint, but are not a unique density inversion.

For the June 2018 lunar picture, start with a regolith and fractured crust, a largely anorthositic crust tens of kilometres thick, a silicate mantle, a deep partially molten region and a small iron-rich core. GRAIL and Apollo-based models gave a global mean crustal thickness around $34$-$43$ km.^[ref] A representative seismic model has a partial-melt boundary near radius $480$ km, a liquid outer-core boundary near $330$ km and a solid inner-core radius near $240$ km; these are model-dependent radii from the centre, not depths below the surface.^[ref]

The long-running Apollo network had only four principal stations on the near side, with poor global coverage and strong scattering in the fractured near-surface material. Far-side structure and weak deep reflected phases were therefore difficult to constrain. The historical date matters: the answer does not import later Mars seismology or later lunar-core refinements into the 2018 premise.

Set $\lambda \simeq d$ with the convention of part b. ThenThus ordinary solid-state temperatures are far below this rough overlap scale. The electrons need quantum statistics and quantum states throughout that regime. The number is an order-of-magnitude wavelength criterion, not an exact Fermi temperature: electron density, electronic band structure, effective mass and the numerical convention all matter.

At the same temperature, momenta increase as the square root of mass, so wavelengths decrease as its inverse square root. For a nucleus of mass $M_N\simeq A{m}_p$,At the electron overlap threshold, ${\lambda }_e=d$, the nuclear wavelength is about $7.0\mathrm{pm}/\sqrt{A}$, far below the $0.3$ nm spacing. Nuclear exchange effects are consequently much less important than electron exchange effects. This statement about nuclear overlap does not mean lattice vibrations are always classical.

At room temperature, electrons require quantum mechanics, including their occupied orbitals or bands and Fermi statistics. Heavy nuclei can often be treated as distinguishable particles near classical equilibrium lattice sites, because their exchange wavelengths are short. Their coupled displacements form phonons, however, and the occupation of a vibrational mode is classical only if $k_{B}T\gg \hslash \omega$ for that mode. Even at room temperature some high-frequency modes remain quantum mechanical.

At $1$ K, quantum lattice vibrations, zero-point motion and the scarcity of thermally excited phonons are essential. Many heavy nuclei still do not have overlapping exchange wavefunctions; quantum vibration and quantum exchange are distinct criteria. For the same spacing, the nuclear overlap threshold from the rough estimate is $T_N\simeq 211/A$ K, illustrating why the nuclear mass cannot be omitted from a blanket low-temperature claim.

At a separation of $1$ fm, the proton Coulomb energy is approximately $V=e^2/(4\pi {\epsilon }_{0}r)\simeq 1.44$ MeV, whereas $k_{B}T\simeq 0.862$ keV at the stated solar-core temperature. As a deliberately crude rectangular-barrier estimate, take a width of $1$ fm and use the two-proton reduced mass $\mu =m_p/2$. ThensoTreating one proton as fixed and using $m_p$ instead gives about $0.59$. Both estimates only say that this arbitrarily shortened barrier would be fairly transparent; the opaque-barrier exponential is itself uncontrolled for such a small exponent.

There is an important conditioning issue. If the protons have already reached a separation of one femtometre, most of the long-range Coulomb suppression has already been overcome. The rectangular estimate cannot represent the probability for thermal protons starting far apart. Their classical turning separation at this energy is about $1.7\times 10^3$ fm, not one femtometre.

Two major suppressions are missing from the one-femtometre rectangular estimate. First, the Coulomb barrier extends far beyond nuclear distances and is not rectangular. A semiclassical integral through it produces the Gamow factor $\exp (-2\pi \eta )$, with $\eta ={\alpha }_{\mathrm{fs}}\sqrt{\mu c^2/(2E)}$; at the illustrative thermal energy it is already of order $10^{-11}$. A thermal average additionally weights the high-energy tail and is not the same as evaluating a single thermal energy.

Second, two protons cannot simply form a stable bound diproton. The proton–proton chain starts with a weak-interaction conversion producing a deuteron, a Positron and a neutrino. The weak transition probability supplies a further very large suppression. Thus the quoted tiny fusion probability is not a pure one-dimensional barrier transmission probability with universal energy-independent meaning.

The Sun contains of order $10^{57}$ nucleons in a hot dense core and offers enormously many encounters over astronomical times. A very small reaction probability per encounter can therefore supply its luminosity while consuming only a small fraction of the fuel per unit time. The small weak reaction rate is also crucial for the long stellar lifetime; a larger core temperature increases burning and triggers expansion and cooling.

The Shockley-Queisser calculation treats a single-junction solar cell with a single band gap, in steady state, with one collected electron-hole pair per absorbed above-gap photon, an ideal absorption edge, thermalized carriers and radiative recombination as the limiting recombination process. The source spectrum, cell temperature, angular acceptance and optical concentration must be specified. Below-gap photons are transmitted; energy above the gap is lost during carrier thermalization; the cell must emit photons by detailed balance, producing a recombination current and an entropy-related voltage loss. Delivering power at finite current also puts the operating voltage below the open-circuit voltage. Nonradiative recombination, resistance and imperfect collection worsen a real cell but are not required for the ideal limit.

Three ways to change those assumptions are tandem solar cells, which divide the spectrum among gaps and reduce both transmission and thermalization losses; hot-carrier solar cells, which aims to collect excess carrier energy before thermalization; and carrier multiplication or spectral conversion, which aims to use a high-energy photon for more than one excitation or combine otherwise unusable sub-gap photons. These concepts have additional practical losses and do not automatically outperform an optimized conventional device. The usual unconcentrated single-junction optimum is approximately one third under a terrestrial solar spectrum, rather than a universal efficiency independent of illumination.

For the geometrical estimate, take the stated horizontal equinox irradiance to set a constant direct-normal reference. The equinox solar-noon elevation is ${\alpha }_e=90^{\circ }-51.75^{\circ }=38.25^{\circ }$, so this reference is $I_n=1/\sin {\alpha }_e\simeq 1.62\mathrm{kW}{\mathrm{m}}^{-2}$. A south-facing panel with tilt $\beta =41.1^{\circ }$ has noon incidence factor $\cos (51.75^{\circ }-\delta -41.1^{\circ })$, where $\delta =\pm 23.44^{\circ }$. Relative to the $1\mathrm{kW}{\mathrm{m}}^{-2}$ plane-of-array rating,These are the results of that idealized geometrical normalization, not credible seasonal production forecasts. Its inferred direct-normal irradiance already exceeds the approximately $1.36\mathrm{kW}{\mathrm{m}}^{-2}$ solar irradiance outside the atmosphere, so the supplied horizontal benchmark cannot describe a purely direct noon beam literally. Real direct and diffuse fractions, atmospheric attenuation at low winter elevation, cell temperature and inverter clipping are unspecified. An alternative assumption of $1\mathrm{kW}{\mathrm{m}}^{-2}$ direct-normal irradiance gives approximately $4.88$ and $4.14$ kW. A physically unique maximum cannot be inferred without an irradiance model.

The stated annual output corresponds to $4100/5=820$ full-rated hours, an average power of $0.468$ kW and a capacity factor of $9.4\%$. Night, cloud, seasonal day length and other losses explain the large difference from an instantaneous sunlit estimate.

For a national scale estimate, take UK area $2.44\times 10^{11}{\mathrm{m}}^2$, population $6.6\times 10^7$, and annual mean solar input $100\mathrm{W}{\mathrm{m}}^{-2}$. Covering five percent with twenty-percent-efficient modules givesThis is comparable with an order-$100$ kWh daily per-person primary-energy budget, and much larger than electricity use alone, of order $10$-$15$ kWh per person per day. It is a gross area-and-flux estimate: access land, shading, seasonal storage, curtailment and conversion between primary energy and electrical end use all affect its interpretation.

Fission divides a heavy nucleus into more tightly bound fragments, releasing energy through the increase in binding energy per nucleon. A neutron-induced chain reaction can sustain heat production in a controlled reactor. Fusion combines light nuclei into more tightly bound products and needs sufficiently high temperature and confinement to obtain a useful reaction rate despite Coulomb repulsion. Both ultimately supply heat that must be converted to electricity, unless a different energy-conversion scheme is explicitly introduced.

A typical fission releases about $200$ MeV per event, corresponding to roughly $8\times 10^{13}$ J per kilogram of fully fissioned heavy fuel. The deuterium-tritium fusion energy of $17.6$ MeV per reaction corresponds to about $3.4\times 10^{14}$ J per kilogram of reacting mixture. These are nuclear reaction energies, not electricity per kilogram of mined or engineered fuel. Actual fuel use, thermal efficiency, parasitic power and incomplete burn change plant output.

Fission has substantial operating experience and can provide high-capacity-factor generation, but faces capital, waste, decommissioning, fuel-cycle and safety constraints. Fusion avoids a self-sustaining neutron chain reaction in the fuel, but its systems still face neutron damage, activated materials, tritium management and demanding heat extraction. Demonstrating a hot plasma or a favorable reaction-energy balance is not the same as demonstrating economical net electricity. The useful comparison distinguishes scientific feasibility, engineering availability and whole-system cost.

A battery stores electrochemical free energy. Representative pack specific energy is of order $0.1$-$0.25$ kWh/kg, with electrical round-trip efficiency roughly $85$-$95\%$. It responds rapidly and can be modular, but energy capacity requires more active material, and ageing depends on time, temperature and cycling.

Pumped storage stores gravitational potential energy. At a $400$ m head with $89\%$ generation efficiency, useful energy is $\eta gh\simeq 0.97$ Wh/kg of water. Its low mass-specific energy is acceptable when inexpensive water and large reservoirs are available. Storage volume and generator rating can be chosen substantially independently, making multi-hour service practical, but suitable sites and civil works are major constraints.

Hydrogen stores chemical energy, with a lower heating value around $33$ kWh/kg of hydrogen before allowing for tank mass. Electrolysis, compression or other storage processing, and conversion back to electricity can reduce electrical round-trip efficiency to roughly $25$-$45\%$, depending on the route. Separate tanks or caverns can make long-duration energy capacity less expensive than adding electrochemical cells, but containment, infrastructure, leakage and conversion equipment matter.

The appropriate comparison includes power rating, energy duration, cycle count, system mass or area, and whether heat or electricity is the final useful product. A large fuel heating value does not imply a high electricity-storage efficiency.

Use the squared three-dimensional endpoint displacements, not their arithmetic mean. The mean squared displacement estimate and corresponding diffusion coefficient areThe factor six is for three-dimensional Brownian motion. For an isolated no-slip sphere, the Stokes–Einstein relation gives the Brownian radius estimate from endpoint displacements:Here room temperature is taken as $298\mathrm{K}$; choosing $293\mathrm{K}$ gives $55\mathrm{nm}$, a much smaller difference than the sampling uncertainty. Under independent increments with a normal distribution, $\mathrm{sd}(\hat{D})/D=\sqrt{2/(3n)}\simeq 0.37$ for $n=5$. The sampling uncertainty of an endpoint diffusivity estimate therefore makes this a rough radius estimate. It assumes negligible drift and hydrodynamic coupling between probes, the supplied bulk viscosity, and actual three-dimensional tracking. Localization noise adds to the squared displacements and would bias this inferred radius downward if uncorrected.

For the three-dimensional worm-like chain, define the bending coefficient ${\kappa }_b$ by $F_b=({\kappa }_b/2)\int |{\partial }_s\mathbf{t}{|}^{2}ds$. Its persistence length is $l_p={\kappa }_b/(k_{B}T)$. Holding that coefficient and the polymer contour length $L$ fixed givesThe temperatures must be absolute. The end-to-end distance of a worm-like chain obeysFor a long flexible chain, $L\gg l_p$, the typical size is $R_{\mathrm{rms}}\simeq \sqrt{2L{l}_p}$, so warming contracts the chain, with $R_{\mathrm{rms}}(T_2)/R_{\mathrm{rms}}(T_1)\simeq \sqrt{T_1/T_2}$. The finite-length formula also decreases as $l_p$ decreases, but the rod limit $L\ll l_p$ changes very little. A large segment count alone does not guarantee $L\gg l_p$. If the material's bending coefficient itself depends on temperature, replace the first ratio by $[{\kappa }_b(T_2)/{\kappa }_b(T_1)]T_1/T_2$; the sign then needs that additional material information. This calculation concerns a generic physical filament, with no biomolecular application.

For monovalent mobile ions, the linearized Debye–Hückel screening length isHere $n_i$ is a number density in ${\mathrm{m}}^{-3}$, not a molarity. Use the question's idealized neutral-water convention, $c_{+}=c_{-}=10^{-7}\mathrm{mol}{\mathrm{L}}^{-1}$, and $T=293.15\mathrm{K}$. Both ion species contribute:For added salt concentration $c_s$ in $\mathrm{mol}{\mathrm{L}}^{-1}$, the ionic strength in this model is $I=c_s+10^{-7}\mathrm{mol}{\mathrm{L}}^{-1}$. With the numerical molarity $I_M=I/(1\mathrm{mol}{\mathrm{L}}^{-1})$,The water contribution is negligible at that salt level. The factor $10^3$ converts molarity into $\mathrm{mol}{\mathrm{m}}^{-3}$; omitting it or counting only one ion species gives a wrong screening length. Exact water autoionization and activities vary with temperature, so the first estimate uses the simplified neutral $10^{-7}\mathrm{M}$ concentrations implied by the question. At the higher salt level, ion-size effects, nonideal activities and changes in permittivity limit the numerical accuracy of this continuum estimate.

For a clean liquid interface at fixed temperature and composition, surface tension is the reversible free-energy cost per added area, $\gamma =(\partial F/\partial A{)}_{T,V,N}$, so $dF=\gamma dA$. Its units are both $\mathrm{J}{\mathrm{m}}^{-2}$ and $\mathrm{N}{\mathrm{m}}^{-1}$. Molecules at an interface have a different environment from those in the bulk; creating more interface changes their interactions and entropy. For constant $\gamma$, the interfacial free energy is $\gamma A$ up to a reference constant. This is not generally identical to the internal surface energy: the interfacial entropy density is $-d\gamma /dT$, giving internal-energy density $\gamma -Td\gamma /dT$.

A line of length $\ell$ in one interface transmits a tangential force $\gamma \ell$, perpendicular to the line. A soap film has two interfaces, giving force $2\gamma \ell$ on a sliding rim. Quasistatic extension by $dx$ creates total area $2\ell dx$, recovering that force from $dF=2\gamma \ell dx$. Uniform positive surface tension favors minimum area at fixed volume, hence an isolated spherical droplet when gravity is negligible. curvature gives the Young–Laplace equation, with outward normal chosen consistently:A spherical liquid drop has $\Delta p=2\gamma /R$; a thin soap bubble with two similar interfaces has approximately $4\gamma /R$. The capillary length $\sqrt{\gamma /(\Delta \rho g)}$ compares gravity with capillarity. surface tension also restores capillary waves, rather than providing a bulk static shear modulus.

Interfacial contamination and added surfactants can change $\gamma$ substantially; a nonuniform value produces tangential Marangoni stress. A rim-force measurement should therefore control temperature, cleanliness and the number of interfaces. Capillary-rise methods additionally need the contact angle, because the vertical component is proportional to $\gamma \cos \theta$. An independent frame-force measurement can distinguish a contact-angle change from a true change of surface tension.

Use polymer volume fraction $\varphi$ horizontally and the Flory interaction parameter $\chi$ vertically. Lower $\chi$ gives better solvent quality. The simplest Flory–Huggins solution theory assigns the dimensionless mixing free energy per lattice siteThus $\chi <1/2$ gives positive excluded volume and a good-solvent coil, $\chi =1/2$ is the long-chain theta condition, and $\chi >1/2$ is the poor-solvent side. It is the net two-segment interaction that vanishes at theta, not every molecular force.

At fixed good-solvent quality, the polymer overlap volume fraction ${\varphi }^{\ast }\sim N{a}^3/R_g^3$ separates isolated coils from a semidilute polymer solution. Above overlap, many chains interpenetrate, but there is still solvent on the segment scale. As $\varphi$ approaches order unity the solution is concentrated; these concentration-regime boundaries are crossovers, not first-order phase transitions. The overlap fraction increases from order $N^{-4/5}$ in a strong good solvent toward order $N^{-1/2}$ at theta. The crossover can be expressed as ${\varphi }^{\ast }\sim N^{-4/5}(1-2\chi {)}^{-3/5}$ while $1-2\chi \gg N^{-1/2}$; it saturates at the theta scaling when that condition fails.

On the poor-solvent side there can instead be coexistence of polymer-poor and polymer-rich phases. The binodal follows a common tangent of $f$; the spinodal followsThe critical point of a Flory–Huggins polymer solution isFor finite $N$, ${\chi }_c$ lies above $1/2$. Therefore poor solvent does not automatically mean demixing at every concentration, and an overlap line must not be mistaken for a coexistence boundary. The original figure uses the stated ideal model for an illustrative chain length; it does not assert a measured phase boundary. A vertical temperature axis would additionally require a material-specific relation $\chi (T)$.

For a thermally elastic network, estimate one elastically active element per mesh volume. Each contributes an elastic free-energy scale of order $k_{B}T$, giving the entropic network modulusat $T=298\mathrm{K}$. Mesh size alone would not determine the modulus of an arbitrary bending-dominated or athermal network; this is the entropic scaling assumption needed for the estimate.

The formal elastic intercept in the given model is consequentlyThe corresponding three-dimensional root-mean-square displacement is about $0.20\mu \mathrm{m}$. The temperature cancels from this estimate because the assumed modulus is itself proportional to $T$. This is the extrapolated plateau from the preceding part, not a nonzero exact zero-lag displacement. Here $m/a=1/2$, so a continuum estimate is only approximate.

In passive microrheology, track dilute inert probe spheres by video microscopy and average squared lag displacements over starting times and probes. An isotropic two-dimensional image gives $⟨\Delta r^2{⟩}_{3d}=(3/2)⟨\Delta r^2{⟩}_{2d}$. Alternatively, dynamic light scattering probes ensemble displacement statistics; the original optical microrheology study compared inferred storage and loss moduli with mechanical measurements.^[ref]

Remove drift, calibrate spatial scale and frame timing, estimate localization noise using an immobilized-bead control, and account for motion blur. Keep probes away from walls and check that changing probe radius or tracer concentration does not change the inferred bulk response unexpectedly. Noise alone can generate a false short-time intercept, while adsorption or confinement can mimic elasticity. No numerical long-time slope is fixed because a numerical viscosity is not supplied.

The tension is supplied by Stokes drag on the free bead. With bead radius $b=0.5\mu \mathrm{m}$, the stipulated isolated-sphere approximation givesThe trap also balances drag on the trapped bead, but that bead's drag is not an additional force transmitted through the chain. The polymer contour length is $L=Na=6\mu \mathrm{m}$. Use a room-temperature estimate $T\simeq 300\mathrm{K}$, consistent with the water-viscosity approximation imported from part (a); the expressions below retain the dependence on $T$.

For a theta solvent, Gaussian chain elasticity gives $F_{\mathrm{el}}=3k_{B}T{R}^2/(2N{a}^2)$ and henceThis is the mean extension along the force, rather than the unforced end-to-end distance of a polymer chain $a\sqrt{N}=60\mathrm{nm}$. At $300\mathrm{K}$, $fa/(k_{B}T)=0.410$; the finite-link freely jointed chain result $R=L[\coth (fa/k_{B}T)-k_{B}T/(fa)]$ is $0.810\mu \mathrm{m}$, close to the Gaussian estimate.

The good-solvent force-extension scaling changes the elasticity. In the tensile blob regime, a blob has size ${\xi }_t=k_{B}T/f$, contains $g\sim ({\xi }_t/a{)}^{5/3}$ segments, and the aligned string of $N/g$ blobs has lengthOnly the scaling and order of magnitude of the good-solvent value are determined: solvent-dependent swelling amplitudes are not supplied. The nonlinear force law is the essential correction to extrapolating a weak-force coil spring.^[ref]

To expose the unrealistic extrapolation, the unforced good-solvent size estimate is $R_0\sim a{N}^{3/5}=151\mathrm{nm}$. Treating this coil as a linear spring with spring constant $3k_{B}T/R_0^2$ predicts $R\simeq 5.17\mu \mathrm{m}$ at $300\mathrm{K}$, although its required condition $f{R}_0\ll k_{B}T$ fails by a factor of about $103$. Dropping the order-one factor three, as in the rough estimate $k\sim k_{B}T/R_0^2$, even gives $R\sim 15.5\mu \mathrm{m}>L$. An extension beyond the polymer contour length is an artefact of extending the coil's linear response outside its range. The tensile blob estimate avoids that violation but here ${\xi }_t\simeq 1.46\mathrm{nm}$ and $g\sim 4.4$, so it is only marginally within a many-segment-blob regime. finite extensibility of a polymer chain, chain drag and flow disturbance by both large beads further limit quantitative realism. The specified neglect of chain drag is a model assumption, not evidence that it is physically negligible.

For a symmetric Flory-Huggins polymer blend, the mixing free energy per lattice site isThe critical conditions $f^{\prime \prime }=0$ and $d^{3}f/d{\varphi }^3=0$ give ${\varphi }_c=1/2$ and ${\chi }_c=2/N=0.002$. ThereforeSince $d\chi /dT<0$, heating improves miscibility: this is an upper critical solution temperature in the stipulated model. The very high value is an extrapolation of the supplied interaction-parameter fit, not evidence that an intact experimental blend can be brought reversibly to this temperature; chemical stability and the fitted temperature range must be checked before interpreting it as an accessible transition.

At $293.15\mathrm{K}$, $\chi =0.132421$ and $N\chi \simeq 132$, well into strong segregation. For the polymer-blend interface width, state the convention explicitly. In the equal-segment, long-chain square-gradient approximation,Its first integral gives ${\varphi }^{\prime }=\sqrt{24\chi }\varphi (1-\varphi )/a$, soThe inverse-central-slope width is consequentlyA convention calling the hyperbolic-tangent scale itself the width instead quotes $0.752\mathrm{nm}$; these describe the same profile. The predicted interface spans only a few statistical segments, so the continuum approximation for gradients is quantitatively marginal. Intrinsic composition width and experimentally broadened width also differ because roughness and finite spatial resolution can broaden a measured profile.

A viscoelastic material stores part of the deformation energy and dissipates part of it. In linear viscoelasticity, a small imposed shear strain produces a stress determined by its history, $\tau (t)={\int }_{-\infty }^{t}G(t-t^{\prime })\dot{\gamma }(t^{\prime })d{t}^{\prime }$. A shear step probes stress relaxation; a stress step probes creep compliance. A purely elastic solid has persistent stress after a held strain, while a Newtonian fluid relaxes that stress and creeps continuously under constant stress.

For a Maxwell model, a spring of shear modulus $G_0$ and a dashpot of viscosity $\eta$ are in series:For the convention $\gamma \propto e^{i\omega t}$,The storage modulus dominates at short times or high frequencies; the loss modulus dominates at long times or low frequencies. The energy dissipated per cycle per volume is $\pi G^{\prime \prime }{\gamma }_0^2$. A Kelvin-Voigt model places spring and dashpot in parallel and instead gives bounded creep $J(t)=(1-e^{-G_{0}t/\eta })/G_0$: a simple illustration of a viscoelastic solid. Real materials can require a spectrum of relaxation times rather than one element.

A small-amplitude oscillatory rheometer measures stress amplitude and phase relative to imposed strain, giving $G^{\prime }=|{\tau }_0/{\gamma }_0|\cos \delta$ and $G^{\prime \prime }=|{\tau }_0/{\gamma }_0|\sin \delta$. First sweep strain amplitude to identify the linear range; control temperature, sample ageing, wall slip and instrument inertia. Frequency-dependent phase lag alone is insufficient if these apparatus effects are unresolved.

DLVO theory adds attractive Van der Waals forces and electric double layer repulsion to obtain an effective colloid pair potential. For two identical spheres of radius $b$ at a surface gap $h\ll b$, a useful schematic form isHere $A_H$ is the Hamaker constant. For a symmetric monovalent electrolyte with number density $n$ of each ionic species, ${\kappa }_D^2=2n{e}^2/(ϵk_{B}T)$. The exponential is a simplified screened repulsion; its amplitude depends on electric potential at the surface, charge regulation and geometry. Constant-charge and constant-potential boundary conditions need not produce the same close-approach law.

A primary attractive minimum near contact may be separated from dispersed configurations by an energy barrier; a shallower secondary minimum can produce loose reversible aggregation. A barrier much larger than $k_{B}T$ can give long-lived kinetic stabilization without making the dispersed state the global free-energy minimum. Adding salt shortens the Debye–Hückel screening length and can lower the barrier, allowing colloidal coagulation. Changing pH can alter surface ionization, so screening alone does not determine stability.

One can vary salt concentration while holding particle concentration, temperature and pH fixed, then follow particle aggregation using light scattering or microscopy. A solvent-only control and an independently measured surface-charge indicator help distinguish a change of screening from adsorption or charge reversal. Classical DLVO theory omits steric stabilization, hydration forces and strong ionic correlations; the formal attraction diverging at contact is cut off by molecular-scale physics.