Codex @codex  0

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

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

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

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

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

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

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

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

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

Cooling makes the oceanic lithosphere denser and contracts it; isostatic adjustment therefore lowers the seafloor. The cooling layer includes substantial mantle, not just the approximately fixed-thickness igneous oceanic crust.

New plate near a spreading ridge is hot and relatively buoyant. As it moves away, conductive cooling penetrates more deeply, and material formerly in the hot asthenosphere becomes part of the colder thermal boundary layer. To first order,where $\alpha$ is the volumetric thermal expansion coefficient, ${\rho }_m$ a reference mantle density, $T_m$ the hot reference temperature, and $t$ plate age. A positive depth-integrated density anomaly makes the column negatively buoyant relative to the ridge column. Isostasy is restored by downward displacement and mantle redistribution; cold rock also occupies less volume.

At a fixed sufficiently deep datum, deepening the seafloor by $\Delta d$ removes an equivalent thickness of rock from the column and replaces it with less dense seawater. In the usual constant-property approximation this giveswhere ${\rho }_w$ is seawater density and $z_c$ lies below the relevant thermal anomaly. This explains why contraction alone and water-loaded seafloor subsidence are not identical numerical quantities.

For a half-space cooling model,so integration yieldsHere $\kappa$ is thermal diffusivity and $z$ is measured downward from the surface of the thermal model. The leading oceanic thermal subsidence is proportional to $\sqrt{t}$. A finite oceanic plate cooling model limits this growth at old ages. Sediment loading, dynamic topography, anomalous crust thickness and hydrothermal heat transport can produce departures from the simple relation.^[ref]

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

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

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

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

At a convergent margin, the incoming oceanic lithosphere is a chemically and mechanically layered package: sediments, oceanic crust and lithospheric mantle. Subduction does not deliver this whole package unchanged to one destination. Solids can be accreted, underplated, exhumed or recycled deeply, while released fluids and melts carry selected components to the overriding plate. The outcome depends on incoming composition and hydration, slab age, convergence rate, thermal structure, interface coupling and subsequent tectonic history.

The incoming material already contains a history. Sediments include terrigenous debris, biogenic carbonates, siliceous deposits, clay minerals, pore water and organic matter. Oceanic crust consists mainly of basalt, sheeted dikes and gabbro; hydrothermal alteration can add water, carbonates and mobile elements before subduction. Beneath it, residual peridotite forms the lithospheric mantle. Faulting at a spreading environment, a transform boundary or an outer-rise bending region can let seawater enter and cause serpentinization. The lithospheric mantle is not uniformly hydrated to its full thickness. Consequently, sediment thickness alone does not measure the incoming volatile inventory.

Sediments: separation starts near the trench. Compaction and increasing effective pressure expel pore fluids, which can escape through the forearc, transport dissolved constituents and influence pore pressure on the plate interface. Some sediments are scraped off into an accretionary prism, shortened by thrusts and retained in the upper plate. Others pass beneath the leading edge, then detach and undergo tectonic underplating beneath the forearc. Deeper underplated material may return to the surface during tectonic exhumation, preserving a record of burial and deformation. At a sediment-starved erosive margin, little prism may form, and tectonic erosion can entrain material from the overriding plate as well. Thus net continental gain and net continental loss are both possible; they can even alternate at one margin.

Sediments that remain attached to the slab undergo metamorphism, dehydration, and potentially decarbonation or carbonate dissolution. The released fluids can transport boron, barium, strontium, lead and other soluble constituents. Carbon can leave as dissolved species, remain in carbonate, or occur as reduced carbon depending on redox conditions and reaction history. Along sufficiently hot paths, sediment-derived melts or mixtures of sediment and altered rock can transfer material into the mantle wedge; buoyant mixed material may rise as a mélange diapir. Cold paths can instead retain more solid material and some volatile-bearing phases to greater depths. These are alternative reaction and transport pathways, not a single compulsory sequence. A carbon-budget study explicitly includes fluid dissolution, sediment diapirism and storage in the upper plate, showing why volcanic carbon dioxide alone cannot close the whole subduction carbon budget.^[ref]

Oceanic crust: metamorphic transformation and chemical filtration. Altered basalt and gabbro are progressively buried. Depending on their pressure-temperature paths and reaction kinetics, they can pass through blueschist facies or warmer assemblages and eventually form eclogite, dominated by garnet and omphacite. Eclogitization increases density and can strengthen negative slab buoyancy, but a kinetic delay can leave metastable rock. It is incorrect to assign every slab the same sequence or the same fixed depths for these changes.

During burial, breakdown of hydrous phases releases water-rich fluids; different mineral assemblages retain different quantities of water. The fluids migrate through or along the interface and into the mantle wedge, where they cause metasomatism and lower the solidus of peridotite. Flux melting of the mantle wedge generates much of the primary magma feeding a volcanic arc. This is different from saying that the entire descending oceanic crust melts beneath the arc. Most of its refractory mass can remain solid. Partial melting of oceanic crust is more plausible for unusually hot or young slabs, or other favorable thermal conditions, and produces a distinct crustal contribution.

The depth of dehydration is controlled by thermal structure and composition. Global thermal-petrological calculations, for example, predict more efficient deep water transport in old, fast, cold slabs than in hot ones, with substantial differences between volcanic crust, gabbro and hydrated peridotite. Such a calculation supports conditional deep retention; its modeled proportions should not be imposed on every margin.^[ref]

Mechanical detachment adds more branches. Upper-crustal slices, seamounts or oceanic plateaus may be accreted or underplated instead of continuing down the slab. A complete or partial ophiolite can be emplaced onto another plate by obduction in suitable convergence histories. Most intact oceanic crust is not saved this way. High-pressure fragments may be entrained in an interface channel or mélange and exhumed; their small preserved volume should not be mistaken for the dominant global mass flux. Magma generated after slab-derived transfer may also pond and crystallize beneath the upper crust by magmatic underplating, which is a different process from mechanical tectonic underplating of solid slices.

Lithospheric mantle: residual rock and a variable water carrier. Most incoming lithospheric mantle is depleted peridotite left after extraction of ridge melts. Its cold, dense, mechanically coherent portion contributes to slab pull and generally sinks with the plate. Hydrated portions contain serpentine and possibly other hydrous minerals. Breakdown of these phases releases water at conditions that depend on the thermal path, feeding further metasomatism of the mantle wedge and sometimes dehydration embrittlement. The latter is a proposed contributor to intermediate-depth earthquakes, not an exclusive explanation for all such events. Very cold paths may transfer some water onward in high-pressure hydrous phases or as hydrogen incorporated in nominally anhydrous minerals.

Do not confuse the incoming hydrated lithospheric mantle with newly serpentinized forearc mantle in the overriding plate. The latter receives slab-derived fluids and can temporarily store their water. Hydrated, relatively buoyant interface material can facilitate return flow and tectonic exhumation, while the colder deeper bulk of the lithospheric mantle continues down. Some peridotite may also be accreted or obducted with an ophiolite. These exchanges mean that even an exposed serpentinite requires structural and chemical context before it can be assigned to its original plate.

Deep recycling and eventual return. The residual eclogite, remaining sediments and lithospheric mantle can enter the deeper mantle. A slab can temporarily stagnate around the mantle transition zone or penetrate into the lower mantle, depending on buoyancy, phase transformations, geometry and mantle viscosity. Crustal material can remain as compositionally distinct domains, mix, or react with surrounding peridotite; it need not retain the near-surface eclogite assemblage indefinitely as pressure rises. Subsequent mantle convection can return recycled components toward spreading regions or mantle plumes, where partial melting incorporates them into new basalt. A primary study of Hawaiian basalt illustrates a specific model in which melts derived from recycled crust react with peridotite to form a pyroxenite source. It is an interpretation of combined chemical evidence, not a unique label attached to every enriched basalt.^[ref]

How the pathways can be observed. Seismic reflection profiles, drilling, wedge geometry and balanced structural sections constrain accretion, fluid escape and tectonic erosion. High-pressure metamorphism, mineral zoning, radiometric dating and pressure-temperature-time paths in exhumed rocks reveal burial, reaction and return. Their reconstructable histories are especially useful when deformation and fluid inclusions are linked to the dated mineral growth. Earthquake distributions, seismic tomography, seismic attenuation and electrical conductivity constrain the geometry and physical state of the slab and mantle wedge, although fluids, melts and temperature can produce nonunique signals.

Arc magmas carry complementary chemical evidence: relatively mobile barium, boron, lead or strontium can trace fluid-mediated transfer, whereas sediment-sensitive thorium and isotope ratios can help identify other contributions. Element mobility changes with fluid composition and pressure-temperature conditions, so no one element is an infallible tracer. Boron isotope ratios, radiogenic strontium isotope ratios and lead isotope ratios should be interpreted with measured incoming sediments and altered oceanic crust. Arc gases and forearc fluid discharge constrain volatile return; ocean-island basalts, mantle xenoliths and diamond inclusions provide possible evidence of deeper recycling. Each observation samples only part of the system.

For any element $X$, a useful control-volume statement isHere permanent accretion is counted outside the chosen reactive control volume, return includes fluid, melt and exhumed-solid exports, and transient underplating or reservoir growth contributes to storage until exported. For a different boundary choice, the terms must be regrouped consistently. The three incoming layers can have very different transfer efficiencies for bulk rock, water, carbon and a particular trace element. A small mass of escaping fluid can dominate an element's return while most of the solid slab continues into the deep mantle.

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

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

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

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

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

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

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

For a real potential, the Time-dependent Schrödinger equation isMultiply this equation by ${\psi }^{\ast }$, multiply its conjugate by $\psi$, and subtract. The potential terms cancel and the product rule converts the remaining terms into a divergence. ConsequentlyThis is the probability current and its continuity equation.

Put $R(r)=\sin (kr)/r$. The supplied radial Laplacian formula gives ${\nabla }^{2}R=-k^{2}R$ for $r>0$. Also $i\hslash {\partial }_t\psi =E\psi$. Therefore the free-particle equation is satisfied whenThe apparent singularity at zero is removable: $R(r)=k-k^3{r}^2/6+\cdots$ extends smoothly there and satisfies the same equation. Since $R$ is real and the time dependence is only a common phase, ${\psi }^{\ast }\nabla \psi =R\nabla R$ is real, so $\mathbf{j}=0$. The density $R^2$ is time-independent, proving the continuity equation explicitly. This is a standing spherical wave with zero net radial current.

The momentum of a material volume $V(t)$ is ${\int }_{V(t)}\rho \mathbf{u}dV$. Newton's second law equates its time derivative to the pressure force and gravity:Apply the given Reynolds transport theorem componentwise to the left side, and the divergence theorem to the pressure term. ThenSince this holds for arbitrary material volumes and the fields are continuous, the integrand vanishes pointwise. Dividing by positive density gives the Euler equations for an inviscid fluidAt a fixed physical time $t_0$, a streamline is a curve tangent to $\mathbf{u}(\cdot ,t_0)$, which may be parametrized by $d\mathbf{x}/ds=\mathbf{u}(\mathbf{x}(s),t_0)$. The chain rule gives $df(\mathbf{x}(s),t_0)/ds=\mathbf{u}\cdot \nabla f=0$, proving constancy along it. The time is held fixed; a time-dependent streamline statement does not by itself assert material conservation along particle trajectories.

Define $L_{\xi }(f)={\lim }_N{A}_N(f)$. Limits preserve linearity and positivity, and $|L_{\xi }(f)|\le \Vert f{\Vert }_{\infty }$, $L_{\xi }(1)=1$. The Riesz-Markov-Kakutani representation theorem supplies a unique Radon probability measure ${\mu }^{\xi }$. By its definition,which is exactly weak-star convergence in $C(\mathbb{T}{)}^{\prime }$.

If $\xi /(2\pi )=p/q$ in lowest terms, the orbit has exactly $q$ points. Averaging a periodic sequence givesThis includes $q=1$. If $\xi /\pi$ is irrational, then $e^{ik\xi }\ne 1$ for every nonzero integer $k$. The preceding calculation gives $L_{\xi }(z^k)=0$ for $k\ne 0$ and $L_{\xi }(1)=1$, exactly the integrals against normalized Haar measure $m$. Density of trigonometric polynomials implies equality on all continuous functions, and uniqueness givesThese are the periodic-orbit and irrational-rotation cases of equidistribution.

Use normalized Haar measure $dm=dt/(2\pi )$. A homogeneous Banach space on the circle is a Banach subspace $B\subseteq L^1(\mathbb{T})$ with $\Vert f{\Vert }_1\le \Vert f{\Vert }_B$, such that translations ${\tau }_{s}f(t)=f(t-s)$ preserve $B$, satisfy $\Vert {\tau }_{s}f{\Vert }_B=\Vert f{\Vert }_B$, and obey $\Vert {\tau }_{s}f-f{\Vert }_B\to 0$ as $s\to 0$.

For continuous $g$, the map $s\mapsto g(s){\tau }_{s}f$ is continuous with values in the complete space $B$. Define its Banach-valued Riemann integral, equivalently its Bochner integral, byThe triangle inequality for Riemann sums and passage to their norm limit give the canonical estimateThe continuous inclusion into $L^1$ carries the same Riemann sums to the ordinary convolution, so this definition agrees almost everywhere with $\int f(t-s)g(s)dm(s)$. Density of continuous functions in $L^1$, completeness of $B$, and the estimate extend the operation uniquely to all $g\in L^1$. The ordinary $L^1$ convolution estimate identifies that extension with the same integral formula.

Definitions and pressure convention. $R(t)$ is the scale factor, $t$ cosmic time and dots time derivatives. $G$ is the gravitational constant, $c$ the speed of light, $\rho$ the mass-equivalent density of the material fluid and $k=+1,0,-1$ the spatial-curvature sign in this normalization. $\Lambda$ here has units of inverse time squared, rather than inverse length squared. Comparison of the acceleration equation with $\ddot{R}/R=\Lambda /3-4\pi G(\rho +3P/c^2)/3$ fixesThis $ϵ$ is a pressure parameter, unrelated to Question 1's radiative efficiency. pressureless matter has $P=0$ and $ϵ=0$. When solving the generalized model below, $ϵ$ is assumed constant.

Evidence for a positive vacuum component. Standardized Type Ia supernovae measure the luminosity-distance–redshift relation. Their corrected apparent magnitudes at cosmological distances imply greater distances than a decelerating matter-dominated model of comparable parameters, favoring late accelerated expansion. Light-curve standardization, interstellar extinction, selection bias and possible population evolution must be controlled rather than treating every supernova as an identical candle.^[ref][ref]

The cosmic microwave background and baryon acoustic oscillation measurements constrain geometry and matter density in combination. Near spatial flatness together with a subcritical matter density leaves room for an additional smooth component; cosmic ages and structure growth supply consistency tests. A positive cosmological constant fits these observations within the standard framework. They do not establish that every allowed dark-energy theory must be a strictly constant vacuum energy; the inference depends on the cosmological and gravitational model.

pressureless matter evolution. Set $k=0$ and $ϵ=0$. Since $\dot{H}=\ddot{R}/R-H^2$, eliminating $\rho$ givesFor $H=b\coth (at)$, $\dot{H}=-ab{\mathrm{csch}}^2(at)$ and ${\coth }^{2}x=1+{\mathrm{csch}}^{2}x$. ThereforeOn the expanding Big-Bang branch with positive $\Lambda$, chooseThe origin $t=0$ is the Big Bang; an arbitrary time translation would shift this origin. Integrating gives $R\propto {\sinh }^{2/3}(at)$. At early times $H\sim 2/(3t)$, while at late times $H\to \sqrt{\Lambda /3}$.

Age and the density fractions. Dividing the flat constraint by $H^2$, with ${\rho }_{\mathrm{crit}}=3H^2/(8\pi G)$, givesThe solution implies $\sqrt{{\Omega }_{\Lambda }}=\tanh (at)$, soHere $H_0=2.3331\times 10^{-18}{\mathrm{s}}^{-1}$, using the printed megaparsec conversion. Taking one year as $365.25$ days givesFor the same $H_0$, an Einstein-de Sitter universe hasA cosmic model must be older than its oldest objects, with additional time allowed for their formation. Old stellar populations are dated through stellar-evolution fitting, particularly the globular cluster main-sequence turnoff, with uncertainties from distances, composition and input stellar physics. A pre-examination globular-cluster age constraint on cosmology obtained a $95\%$ lower cosmic-age limit of $11\mathrm{Gyr}$ after allowing for formation time. This conflicts with the $9.05\mathrm{Gyr}$ model but is compatible with the $13.09\mathrm{Gyr}$ calculation.^[ref] Individual stellar age estimates are not exact clocks, and the numerical $13.09\mathrm{Gyr}$ is the result for the specified $H_0$ and matter fraction, not a claim to be a contemporary best-fit age.

Constant nonzero pressure parameter. Eliminating $\rho$ before setting the pressure to zero instead givesHenceFor constant $ϵ>-3$, the positive-density expanding Big-Bang branch isEquivalently use $t^{\prime }=(3+ϵ)t/3$ in the pressureless matter solution. The scale factor obeys $R\propto {\sinh }^{2/(3+ϵ)}[(3+ϵ)H_{\Lambda }t/2]$, and the early limit is $H\simeq 2/[(3+ϵ)t]$. The age of a flat barotropic matter-vacuum universe with the same present density fractions and $H_0$ isThusThe small positive pressure makes the early universe expand with a smaller age at the same late expansion parameters.

The acceleration transition. Normalize the present scale factor to one. Constant $ϵ$ and the continuity equation give $\rho ={\rho }_0(1+z{)}^{3+ϵ}$. Dividing the acceleration equation by $H_0^2$ givesAt the barotropic matter-vacuum acceleration transition its right-hand side vanishes. For the parameters here, $1+ϵ>0$, so the past transition isAt greater redshift the material term dominates and expansion decelerates; at smaller redshift vacuum energy dominates and expansion accelerates. The pressure factor $1+ϵ$ belongs in the numerator of the gravitating material term; omitting it would confuse this model with pressureless matter.

Three tracers give complementary estimates of a cluster's gravitating mass. First, galaxy redshifts measure a line-of-sight velocity dispersion ${\sigma }_{\mathrm{los}}$. For a roughly stationary bound cluster, the virial theorem gives a cluster velocity-dispersion mass estimatewith a coefficient depending on the density profile and orbital anisotropy. The galaxy luminosity and plausible stellar mass-to-light ratio supply too little mass to bind typical rich clusters. Foreground interlopers, substructure, velocity anisotropy and a recent merger can bias the dynamical estimate; a large dispersion is not by itself proof of equilibrium.

Second, diffuse X-ray emission measures a hot intracluster plasma, not merely unresolved starlight. Spectral continuum and lines constrain its temperature, while surface-brightness profiles constrain its density after line-of-sight projection and emissivity modeling. If the gas has constant composition and is approximately spherical and in hydrostatic equilibrium,This cluster hydrostatic mass estimator uses $\mu m_p$ as the mean mass per pressure-supporting particle. The inferred potential requires substantially more mass than the stars plus measured gas, although the gas contains much of the visible baryonic mass. Turbulence, bulk motion, temperature structure and departures from hydrostatic equilibrium are important systematic errors.

Third, strong-lensing arcs and multiple images, together with coherent weak-lensing shear of background galaxies, measure projected gravitational deflection. In general relativity,This cluster gravitational-lensing mass estimate does not require the cluster galaxies or gas to be in dynamical equilibrium. It does require source-distance information, a lens model and control of point-spread-function-induced shapes, line-of-sight projection, unrelated line-of-sight structure and the mass-sheet degeneracy.

The convergence of dynamical, hydrostatic and lensing masses, exceeding the baryonic inventory, is the principal evidence. The Bullet Cluster adds a spatial test: collisional X-ray gas lags behind the galaxies during the merger, while the inferred lensing concentrations lie nearer the collisionless galaxy component. An additional weakly collisional gravitating mass accounts naturally for that separation.^[ref] The conclusion assumes a gravitational framework and does not identify the microscopic dark-matter particle; the different tracer systematics should be evaluated rather than treating three mass estimates as completely assumption free.

cosmic inflation is an early interval of accelerated expansion, $\ddot{R}>0$. For expansion this means a shrinking comoving Hubble radius, $d[(RH{)}^{-1}]/dt<0$. A nearly constant positive scalar potential can supply it: in $c=1$ units,Potential domination gives ${\rho }_{\varphi }+3P_{\varphi }<0$ and approximately exponential expansion. Sustaining the interval requires sufficiently slow evolution, not just setting the initial field velocity to zero on a steep potential.

The horizon problem is the near equality of conditions in widely separated microwave-background regions that could not equilibrate in the standard decelerating hot Big Bang history. cosmic inflation permits the currently observed patch to descend from an earlier region small enough to have been causally connected, before that patch is stretched to enormous size. The claim concerns a causal history extended before the conventional hot phase; merely accelerating a collection of already disconnected regions would not explain their common temperature.

For the flatness problem, the Friedmann constraint givesA growing $RH$ drives this curvature fraction toward zero, whereas ordinary matter- or radiation-dominated deceleration makes it grow. cosmic inflation therefore suppresses observable curvature without requiring an equally extreme initial tuning within the later hot evolution.

Pre-existing monopoles and other unwanted relics are also diluted by the volume increase: a conserved particle number density falls as $e^{-3N}$ over $N$ e-folds. The transition to the hot universe must reheat efficiently, and its temperature must not regenerate an excluded relic population. The required number of e-folds is model and reheating dependent; “about sixty” is a useful scale for many scenarios rather than a universal exact requirement.

cosmic inflation also supplies a mechanism for perturbations. quantum fluctuations are stretched beyond the Hubble radius, and nearly scale-invariant comoving curvature perturbations can later seed cosmic structure. Simple single-field slow-roll models commonly predict approximately adiabatic, nearly Gaussian fluctuations, with a small spectral tilt; tensor perturbations would provide an additional test. Acoustic coherence and the observed scalar spectrum are compatible with this account, but they do not single out one inflaton potential. cosmic inflation addresses causal homogeneity, curvature and relic dilution while motivating the initial perturbations; it does not automatically solve the cosmological constant problem or determine reheating.

Use signature $(+---)$, coordinates $x^0=ct$ and four-velocity $u^{\mu }=d{x}^{\mu }/d\tau$ normalized by $u_{\mu }{u}^{\mu }=c^2$. The derivative is the metric-compatible Levi-Civita connection. Fix the curvature sign by $[{\nabla }_{\alpha },{\nabla }_{\beta }]V^{\mu }=R^{\mu }{}_{\nu \alpha \beta }{V}^{\nu }$. Then the required contractions areWith $G^{\mu \nu }=R^{\mu \nu }-g^{\mu \nu }R/2$, the contracted Bianchi identity gives ${\nabla }_{\mu }{G}^{\mu \nu }=0$. The Einstein field equations $G^{\mu \nu }+\Lambda g^{\mu \nu }=(8\pi G/c^4)T^{\mu \nu }$, with constant $\Lambda$, therefore require covariant conservation of stress-energy. This is a consistency relation between geometry and the total matter/field source. It does not mean ${\partial }_{\mu }{T}^{\mu \nu }=0$ in arbitrary coordinates, nor guarantee a global conserved energy in an arbitrary curved spacetime. The geometric identity alone, without the field equations or matter dynamics, is not an equation of motion for a freely specified material source.

For relativistic dust, $\rho$ is proper rest-mass density, measured in the fluid rest frame, and $u^{\mu }$ is its timelike four-velocity with the $c^2$ normalization above. In a local Cartesian inertial frame $u^{\mu }=\gamma (c,\mathbf{v})$, soAccordingly $T^{0i}/c$ is momentum density and $c{T}^{0i}$ is energy flux in direction $i$. The equality of momentum-density and energy-flux components after the appropriate factors of $c$ reflects the symmetry of the stress-energy tensor. In the comoving frame the spatial stress is zero, as required for pressureless matter.

Expand the dust conservation equation by the product rule:Contract with $u_{\nu }$. metric compatibility and constant normalization give $u_{\nu }{\nabla }_{\mu }{u}^{\nu }={\nabla }_{\mu }(u_{\nu }{u}^{\nu })/2=0$, leaving $c^2{\nabla }_{\mu }(\rho u^{\mu })=0$. Substitute this rest-mass continuity law back and divide by nonzero $\rho$:This proves geodesic motion of uncharged dust along each fluid flow line. pressure forces or another exchange of momentum would change the argument; at zero density there is no dust worldline to constrain.

The electromagnetic field tensor is antisymmetric and has six independent components. In the stated local inertial coordinates choose $F^{0i}=-E_i/c$ and $F^{ij}=-{ϵ}_{ijk}{B}_k$, givingThe fields are those measured by that local observer; changing observer mixes them through the Lorentz transformation of $F$. With ${ϵ}^{0123}=+1$ and ${\tilde{F}}^{\mu \nu }={ϵ}^{\mu \nu \alpha \beta }{F}_{\alpha \beta }/2$, two independent quadratic contractions areThe second is a pseudoscalar: it is invariant under proper Lorentz transformations and changes sign under an orientation reversal. Opposite choices for the orientation tensor change its written sign. A nonzero null electromagnetic field can have both contractions zero, so the two numbers do not determine whether the tensor itself vanishes.

For charged relativistic dust, insert $j^{\alpha }=(q/m)\rho u^{\alpha }$ into the force-density equation:Contraction with $u_{\nu }$ again gives the continuity law, because $u_{\nu }{F}^{\nu \alpha }{u}_{\alpha }=0$ by antisymmetry. After removing the continuity term,In a local inertial frame the spatial components reduce to $d(\gamma m\mathbf{v})/dt=q(\mathbf{E}+\mathbf{v}\times \mathbf{B})$ and the temporal component gives $d(\gamma m{c}^2)/dt=q\mathbf{E}\cdot \mathbf{v}$. This is the relativistic Lorentz force, with gravity represented by the connection. The force is orthogonal to $u$ and preserves the rest-mass normalization. Matter alone now exchanges energy and momentum with the electromagnetic field; adding the field tensor's stress-energy restores zero divergence of the combined source.

A thin spherical shell has outward pressure force $-4\pi r^2(dP/dr)dr$ and inward gravitational force $GM(r)4\pi r^2\rho (r)dr/r^2$. Hydrostatic equilibrium requiresA uniform-density solar model gives $M(r)=M_{\odot }(r/R_{\odot }{)}^3$, $\bar{\rho }=3M_{\odot }/(4\pi R_{\odot }^3)$, and after integration from $R_{\odot }$ to zero,Thus the appropriate orders of magnitude are $10^{14}$–$10^{16}\mathrm{Pa}$ and $10^7\mathrm{K}$. The uniform-density approximation underestimates the actual central pressure because the Sun is centrally concentrated. Radiation pressure is neglected in applying the ideal gas law.

A free-fall timescale estimate follows from $R/t^2\sim GM/R^2$:Falling to a central point mass gives the more precise factor $\pi /(2\sqrt{2})$, or $1.77\times 10^3\mathrm{s}$; a distributed star changes an order-one coefficient. For any spherical mass distribution with radius $R$,The inequality is strict for an extended star; equality is approached by a thin shell at $R$. For an monatomic gas the virial theorem is $2U+\Omega =0$, so half the gravitational potential energy released by contraction heats the star and half can be radiated. The Kelvin-Helmholtz cooling time is therefore of orderA uniform-density model gives $3G{M}^2/(10RL)\simeq 9.4\times 10^6\mathrm{yr}$. The given bound fixes a scale, rather than an exact contraction history; both estimates are far shorter than the Sun's main-sequence lifetime.

The dominant $z^2/{\beta }^2$ factor raises the energy loss as a charged projectile slows, while the logarithm moderates that rise. This produces a Bragg peak close to the end of its particle range. At still lower velocity, changing effective charge, shell effects and nuclear stopping invalidate extrapolation of the elementary Bethe stopping formula to zero energy. A finite initial energy gives a finite continuous-slowing-down range,A higher initial energy increases the particle range; greater electron density increases stopping power per unit length. Absorber $Z/A$ and $I$ also affect the depth. Multiple Coulomb scattering, energy-loss straggling, nuclear reactions and the incident energy spread broaden the measured peak and distal fall-off; nuclear secondaries can produce a small tail. Superposing many incident energies makes a spread-out Bragg peak for an extended target. The schematic curve below illustrates shapes only and is not patient radiation dosimetry.

Put $x=E_I/(k_{B}T)$. With $\zeta =1/2$, the heat capacity formula becomesTaking $E_I=13.6\mathrm{eV}$ gives $\boxed{T=1.02\times 10^4\mathrm{K}}$. At half ionization, $n_e{n}_p/n_H=n_e$, so the Saha ionization equation directly givesThis uses $2g_p/g_H=1$ as specified in part (a); a different internal-state convention must be carried consistently into the numerical electron density. The original nuclei density is $n_0=2n_e\simeq 10^{21}{\mathrm{m}}^{-3}$.

A thin spherical shell has outward pressure force $-4\pi r^2(dP/dr)dr$ and inward gravitational force $GM(r)4\pi r^2\rho (r)dr/r^2$. Hydrostatic equilibrium requiresA uniform-density solar model gives $M(r)=M_{\odot }(r/R_{\odot }{)}^3$, $\bar{\rho }=3M_{\odot }/(4\pi R_{\odot }^3)$, and after integration from $R_{\odot }$ to zero,Thus the appropriate orders of magnitude are $10^{14}$–$10^{16}\mathrm{Pa}$ and $10^7\mathrm{K}$. The uniform-density approximation underestimates the actual central pressure because the Sun is centrally concentrated. Radiation pressure is neglected in applying the ideal gas law.

A free-fall timescale estimate follows from $R/t^2\sim GM/R^2$:Falling to a central point mass gives the more precise factor $\pi /(2\sqrt{2})$, or $1.77\times 10^3\mathrm{s}$; a distributed star changes an order-one coefficient. For any spherical mass distribution with radius $R$,The inequality is strict for an extended star; equality is approached by a thin shell at $R$. For an monatomic gas the virial theorem is $2U+\Omega =0$, so half the gravitational potential energy released by contraction heats the star and half can be radiated. The Kelvin-Helmholtz cooling time is therefore of orderA uniform-density model gives $3G{M}^2/(10RL)\simeq 9.4\times 10^6\mathrm{yr}$. The given bound fixes a scale, rather than an exact contraction history; both estimates are far shorter than the Sun's main-sequence lifetime.