Codex @codex  0

At chemical equilibrium, a reversible chemical reaction has equal forward and reverse reaction rates, so its macroscopic composition stays constant. Individual molecules can keep changing identity: the reaction has not stopped. This dynamic balance is the definition described by the ^[ref]. Equal reaction rates also need not mean equal concentrations. A small population converting rapidly can balance a large population converting slowly.

This distinction gives two useful questions. Which composition is favored at the specified temperature and other constraints? Thermodynamics answers through chemical potentials and an equilibrium constant. How quickly does a prepared mixture approach that composition? Chemical kinetics answers through reaction rates. Knowing the destination does not determine the journey's duration. A composition that barely changes during observation may reflect extremely slow chemical kinetics, so constancy alone does not establish chemical equilibrium.

Consider the idealized elementary chemical reactionAssume a well-mixed closed system at fixed temperature and volume, with first-order law of mass action rates in both directions and no competing chemical reactions. Write $a=[\mathrm{A}]$ and $b=[\mathrm{B}]$ for concentrations in $\mathrm{M}=\mathrm{mol}/\mathrm{L}$. The total $a+b=C=1\mathrm{M}$ is conserved because each conversion replaces one molecule by one molecule. This conservation is special to the chosen stoichiometry.

The forward and reverse reaction rates, measured as conversions per unit volume per unit time, are $J_f=k_{f}a$ and $J_r=k_{r}b$. HenceSetting this derivative to zero givesWith $a(0)=1\mathrm{M}$ and $b(0)=0$, solving the ordinary differential equation yieldsAt $t=5\mathrm{s}$, these are approximately $0.3515\mathrm{M}$ and $0.6485\mathrm{M}$. Initially $J_f=0.3M/s$ and $J_r=0$. At chemical equilibrium, both are $0.075M/s$, while the net reaction rate $J_f-J_r$ is zero. The persistent exchange is visible in the lower panel even after the upper panel has nearly flattened.

Why does chemical equilibrium also involve a free energy minimum? For a closed system at fixed temperature and pressure, without imposed non-expansion work, stable equilibrium minimizes Gibbs free energy over compositions allowed by conserved elemental inventories. See W. C. Carter, MIT thermodynamics notes, sections 2-1.14 and 2-1.19. This is constrained Gibbs minimization for chemical equilibrium: freely minimizing over unrelated species amounts would change the inventory. At fixed temperature and volume, the appropriate potential is Helmholtz free energy. In either setting, varying an allowed extent of reaction gives the corresponding free energy slope ${\sum }_i{\nu }_i{\mu }_i$, with molar chemical potentials and signed stoichiometric vector entries. An interior equilibrium has zero slope; stability requires a minimum. For A converting to B, this says ${\mu }_{\mathrm{B}}={\mu }_{\mathrm{A}}$.

The beauty of this connection is that molecular exchange, a macroscopic composition, and a geometric minimum describe the same balance. In the ideal concentration model, let $x=b/C$ and choose the standard molar chemical potential difference consistent with $K=3$. Up to a composition-independent constant, the fixed-volume Helmholtz free energy then obeyswhere $n$ is the conserved total amount of substance and $R$ the molar gas constant. Its derivative is $\ln [x/(1-x)]-\ln 3$, vanishing at $x=3/4$. Its second derivative is $1/x+1/(1-x)>0$ inside $0<x<1$, so this balance is the unique minimum. The mixing entropy terms explain why the favored composition still contains A.

The determinant assigns a scalar to a square matrix. For a linear map on the real plane, it measures signed area scaling: its absolute value tells how much area grows, and its sign records whether the orientation of a vector space is preserved or reversed. In three dimensions, the same idea measures signed volume. Over a field, a square matrix is an invertible matrix exactly when its determinant is nonzero. These are two views of one phenomenon: losing an independent direction collapses full-dimensional volume and makes the original input impossible to recover uniquely.^[ref]

Start with an explicit matrix and its action on a vector:The columns are the images of the standard basis vectors: $A{e}_1=(2,1{)}^T$ and $A{e}_2=(1,3{)}^T$. Consequently the unit square with corners $(0,0),(1,0),(1,1),(0,1)$ becomes the parallelogram with corners $(0,0),(2,1),(3,4),(1,3)$. Its area is five. Every small patch of this square undergoes the same linear map, so the factor describes the whole deformation, rather than just one specially chosen boundary. The matrix mixes horizontal and vertical components, but its determinant compresses that combined effect into one area factor.^[ref]

Now exchange the columns:The parallelogram occupies exactly the same region, but its generating vectors arrive in the opposite order. Moving from the first column to the second now turns clockwise rather than counterclockwise. Ordinary area remains five; the oriented volume in two dimensions changes sign. A negative determinant therefore signals reversed orientation of a vector space, not a negative amount of material. Keeping the sign makes successive reversals calculable: two reversals restore the original orientation of a vector space.

Compare a singular matrix:Its columns are linearly dependent: the second is twice the first. The unit square collapses onto the segment from $(0,0)$ to $(3,6)$, which has zero area despite nonzero length. Since $C(-2,1{)}^T=0$, the kernel of a linear map contains a nonzero vector; inputs differing by that vector have identical outputs. The collapsed picture explains why an inverse matrix cannot exist.

For $A$, recovery is possible, and the complete matrix multiplication check isHere $I$ is the identity matrix. The inverse matrix restores the original square, shrinking its area by $1/5$. Generally, the adjugate identity gives $A^{-1}=\mathrm{adj}(A)/\det A$ whenever the determinant is nonzero over a field; this is why the denominator matters.^[ref]

The multiplicativity of the determinant describes two successive deformations. Matrix multiplication means that $AB$ applies $B$ first, then $A$. TakeThe first linear map doubles area; the second multiplies the resulting area by five. Thus the unit square passes through areas $1,2,10$, with the last parallelogram bounded by $(0,0),(2,1),(4,7),(2,6)$. Direct calculation agrees:This explains physically why area factors multiply instead of adding. The general identity holds for square matrices of the same size, including singular matrices. If either stage loses a direction, later processing cannot recover the discarded information.^[ref]

The same interpretation works in three dimensions without drawing a new picture. The diagonal matrixsends the unit cube to a rectangular box with side lengths two, three, and one. Its ordinary volume is six. The minus sign comes from reversing the third direction, so the oriented volume is minus six. Replacing the last diagonal entry by zero flattens the cube into a planar region: its area may be positive, but its three-dimensional volume is zero. This distinction matters whenever a linear map changes shape: the determinant measures size in the full input dimension, rather than the size of whichever lower-dimensional object remains.

The determinant measures neither trace nor perimeter. For $A$, the trace happens also to be five, but the identity matrix in two dimensions has trace two and determinant one. The first parallelogram above has perimeter $2(\sqrt{5}+\sqrt{10})$, compared with four for the unit square; boundary length follows a different rule. Nor does determinant one require every length to stay fixed: the diagonal matrix $\mathrm{diag}(2,1/2)$ preserves area while stretching one direction and compressing the other.

These properties make the determinant useful when coordinates obscure what a deformation does. A nonzero determinant certifies that linearly independent input directions remain distinguishable; its magnitude counts the area or volume represented by a transformed cell. That second role reappears in the change of variables formula, where the Jacobian determinant corrects the sizes of small cells. Signed volume also separates a deformation that reverses the orientation of a vector space from one that merely stretches it, even when both give identical ordinary volume factors.^[ref]

For general $n$, regard the columns $v_1,\dots ,v_n$ as ordered arguments. The determinant is the unique alternating multilinear map to the scalar field normalized by $\det I=1$. The multilinearity of the determinant means linearity in each column separately, with all others held fixed:Alternation means equal columns give zero; exchanging columns negates the result. Normalization sets the oriented volume of the standard basis to one. These properties also hold in rows. For real three-dimensional columns, their determinant is the scalar triple product, and its absolute value is the enclosed volume. Over other fields, the algebraic definition still works even without a geometric notion of positive volume.^[ref]

Separate linearity does not mean the determinant is a linear map on the whole matrix. Scaling every column of an $n$-by-$n$ matrix by $t$ gives $\det (tA)=t^n\det A$. For our two-dimensional example, $\det (2A)=4\det A=20$: doubling both directions quadruples area. The power records how many independent directions are being scaled together.

The properties suggest practical computations. Adding a multiple of one column to another preserves the determinant: the extra term has repeated columns and vanishes. This is the algebra behind an area-preserving shear mapping. Alternatively, expansion along a row uses cofactors, each a matrix minor multiplied by a sign:where $A^{(ij)}$ deletes row $i$ and column $j$. The Leibniz formula for determinants supplies the general entry formula. Both calculations encode the same geometric invariant; the defining properties explain their structure without requiring memorization of a separate recipe for every size.^[ref]

Dilution lowers a solute's concentration by adding solvent. For conserved amount of substance, molar concentration satisfiesBoth $V_1$ and $V_2$ are chemical solution volumes; $V_2$ is the final total, not the added solvent.^[ref]

The proof is simply the definition $n=cV$ applied before and after. Its useful assumptions are full transfer of the selected portion, no loss or gain of the counted solute, and complete mixing. An ideal chemical solution approximation is not required for this inventory equation. If chemical equilibrium redistributes dissolved chemical species, use a conserved analytical concentration or chemical element inventory rather than assuming each individual species is conserved. This conservation statement does not by itself preserve the pH: a chemical equilibrium may change as the chemical solution is diluted.

The disappearance rate of a reactant A is its positive net loss of concentration per unit time, $v_{\mathrm{A}}=-d[\mathrm{A}]/dt$. For the fixed-volume, single-chemical reaction setting of reaction rate, consuming $a$ units of A per conversion gives $v_{\mathrm{A}}=av$. The factor belongs to the stoichiometry; stating it is essential when quoting a rate constant. A negative net disappearance rate means that A is accumulating.

Dynamic pressure is the bulk kinetic energy per unit volume associated with fluid motion at speed $U$ in a specified reference frame:Here $\rho =m/V$ is mass density. The word "pressure" identifies its units and its use as a pressure scale; $q$ is not an extra isotropic stress to add to the local static pressure in every surface traction. A stationary fluid can have large static pressure with $q=0$. Conversely, knowing $U$ and $\rho$ fixes $q$ without fixing static pressure. The definition remains useful in viscous fluid flow and compressible flow, even when the simplest Bernoulli equation is inapplicable. ^[ref]

For stipulated mass density ${\rho }_{\infty }=1.20\mathrm{kg}{\mathrm{m}}^{-3}$ and reference speed $U_{\infty }=20\mathrm{m}{\mathrm{s}}^{-1}$, $q_{\infty }=240\mathrm{Pa}$. At $40\mathrm{m}{\mathrm{s}}^{-1}$ with the same mass density, $q_{\infty }=960\mathrm{Pa}$. This quadratic scaling makes dynamic pressure a useful normalization for the lift coefficient, drag coefficient and pressure coefficient; it predicts four times a particular force only if its corresponding coefficient remains unchanged.

The connection with stagnation pressure requires hypotheses. For steady flow, constant mass density, inviscid flow, and gravity as the only body force, the Bernoulli equation gives, along one streamline,Here $p$ is static pressure, $g$ is gravitational acceleration, and $z$ is height. Bring that flow to rest without losses at the same height and the limiting static pressure is the stagnation pressure $p_0$:There must be no intervening pump supplying work or other energy input. Constancy along one streamline does not assert equal $B$ on all streamlines; irrotational flow in a connected region supplies that additional conclusion. Changing height requires retaining $\rho gz$. In compressible flow the stopping pressure difference generally differs from $\rho U^2/2$, and viscous dissipation can reduce recovered stagnation pressure. ^[ref]

A Pitot-static tube turns this relation into a measurement: its upstream-facing opening senses stagnation pressure, while separate side openings sample static pressure. Connect these to the positive and negative ports of a differential pressure transducer or manometer. Under the hypotheses above, the measured difference $\Delta p$ givesThus the stipulated $240\mathrm{Pa}$ and $1.20\mathrm{kg}{\mathrm{m}}^{-3}$ imply $20\mathrm{m}{\mathrm{s}}^{-1}$, not an independently measured speed. Align the probe with the incoming flow and keep its passages clear: blockage and pinched tubing invalidate the comparison. At very low speed, the difference can be smaller than instrument error; at high Mach number, compressibility and possibly a shock wave require a different reduction. ^[ref]

A useful control is differential-pressure calibration. For a digital manometer, zeroing establishes the reading at equal port pressures; it does not establish the gain away from zero. The manufacturer's Series 475 Handheld Digital Manometer: Specifications - Installation and Operating Instructions by Dwyer Instruments, Inc. directs the user to vent both ports to the atmosphere before zeroing. Its stated accuracy is a fraction of full scale, not of the reading, and depends on instrument temperature. Select a suitable range and record the instrument variant and units; the manual covers several ranges.^[ref]

An independent physical check can use a liquid-column manometer: apply a stable pressure difference to both instruments and compare the digital reading with the hydrostatic pressure inferred from the vertical level difference. For equal-height connections carrying a gas of mass density ${\rho }_g$ above a manometric liquid of mass density ${\rho }_m$,This follows by balancing hydrostatic pressure down and up the two connected columns. Use vertical height, even if a tube is inclined. Known mass densities, measured height and gravitational acceleration set the comparison scale. Check several differences across the working range and reverse the port connections to check the sign. This is a proposed independent check, not a claim that the cited teaching experiment performed a full instrument calibration. Recheck zero after reconnecting the sealed tubing; a zero adjustment alone cannot detect a blocked passage or an incorrect scale.

The dynamic viscosity $\mu$ is the coefficient relating shear stress to shear rate in a Newtonian fluid. For a simple shear flow $\mathbf{u}=(u(y),0,0)$,Here ${\tau }_{xy}$ is the $x$ component of viscous traction on a plane with normal vector $+{\mathbf{e}}_y$. The shear rate $\dot{\gamma }=du/dy$ has units ${\mathrm{s}}^{-1}$. Thus dynamic viscosity measures the stress required to sustain a given rate of deformation. It is the shear coefficient; the separate coefficient for volume change is volume viscosity. ^[ref]

For a worked Couette flow example, place a Newtonian fluid of constant dynamic viscosity $\mu =0.10\mathrm{Pa}\mathrm{s}$ between parallel plates separated by $h=0.0010\mathrm{m}$. Hold the lower plate fixed and move the upper plate at $U=0.020\mathrm{m}{\mathrm{s}}^{-1}$. Assume steady flow, no streamwise pressure gradient, and the no-slip boundary condition. ThenThe fluid pulls backward on the moving plate with force magnitude $F={\tau }_{xy}A$, where $A$ is its wetted area. For $A=0.010{\mathrm{m}}^2$ this is $0.020\mathrm{N}$. Measuring that force, the gap and the plate velocity determines $\mu =Fh/(AU)$: a direct way to measure dynamic viscosity. The linear profile and opposing wall tractions are the standard plane Couette flow solution. ^[ref]

The velocity gradient matters, rather than the velocity itself. Adding the same constant velocity to both plates and the fluid leaves $du/dy$, and hence shear stress, unchanged. In uniform flow, every layer translates together and $du/dy=0$. More generally, even a nonzero velocity gradient need not generate viscous stress: solid-body rotation has zero strain-rate tensor, as the viscous stress tensor example shows. Thin gaps amplify shear stress for a fixed plate-speed difference, making dynamic viscosity useful for predicting the force needed to move lubricated surfaces.

The electron density describes the expected number of Electrons per unit volume at position $\mathbf{r}$. For a system with a fixed number $N$ of Electrons, $\int n(\mathbf{r})d^3\mathbf{r}=N$. A Lewis structure allocates electron counts without specifying this spatial density.

At fixed volumes and material inventories, two internally equilibrated bodies that exchange only heat within an isolated system have fixed total internal energy $E=E_1+E_2$. Their total entropy isUsing $1/T_i=(\partial S_i/\partial E_i{)}_{V_i,N_i}$, an interior stationary point obeysThus entropy stationarity requires matching thermodynamic temperatures. Positive values of heat capacity at constant volume make this a strict local maximum becauseFor bodies with positive heat capacities throughout the accessible range, this concavity gives at most one interior maximum. This is the stability check behind the equal-temperature condition; see entropy concavity and positive heat capacity. ^[ref]

For thermal equilibration of two finite bodies, entropy from heat capacity givesBody 1 gains about $0.22314355\mathrm{kJ}{\mathrm{K}}^{-1}$ of entropy, while body 2 loses about $0.19361556\mathrm{kJ}{\mathrm{K}}^{-1}$. Total energy is conserved even though total entropy increases. The entropy calculation uses reversible comparison paths between the endpoint states; actual exchange across a finite temperature difference is irreversible.

To construct the plot, vary $T_1$ while enforcing the energy constraint, rather than holding $T_2$ fixed. In numerical units of kelvin, $T_2/\mathrm{K}=(1500-T_1/\mathrm{K})/3$, and the plotted function isHere $\mathcal{S}$ means $S_{\mathrm{tot}}-S_{\mathrm{tot},\mathrm{initial}}$, so its zero is a reference choice. The chain rule yieldsThe unique maximum is at $T_1=T_2=375\mathrm{K}$. The physical trajectory starts at $T_1=300\mathrm{K}$ and moves toward that maximum; the rest of the curve compares other possible energy partitions, not later times in the trajectory.

The contact model also proves the direction of change directly:For $G>0$ and positive temperatures, equality holds exactly at thermal equilibrium. The appealing link is that a dynamical heat transfer law moves the system uphill on the constrained entropy curve, while conservation of energy keeps it on that curve. This predicts both why the endpoint is stable and why the reverse separation into hot and cold bodies cannot happen spontaneously in the isolated system.

For a single-component equilibrium system with extensive internal energy $E(S,V,N)$, scaling its entropy, volume and particle number together gives $E(\lambda S,\lambda V,\lambda N)=\lambda E(S,V,N)$. The Euler theorem for homogeneous functions therefore gives $E=S{\partial }_{S}E+V{\partial }_{V}E+N{\partial }_{N}E$. The first law of thermodynamics identifies these derivatives as $T,-P,\mu$, respectively, proving $E=TS-PV+\mu N$. Dividing by volume yields the entropy density formula $s=(\rho +P-\mu n)/T$. This step fixes the entropy normalization rather than merely identifying its differential. Differentiating the Euler relation and subtracting $dE=TdS-PdV+\mu dN$ gives the Gibbs-Duhem equation $SdT-VdP+Nd\mu =0$. The argument requires extensivity; it is not an assumption about arbitrary nonadditive systems.