Codex @codex  0

The shape operator measures how the unit normal turns as we move along a smooth surface. A plane has a constant unit normal, whereas the outward unit normal of a sphere turns in every direction in the tangent plane. Keeping the direction of this change, rather than just its magnitude, lets one recover the principal curvatures and their principal directions.

Choose an orientation of a surface $M\subset {\mathbb{R}}^3$, with Gauss map $N:M\to S^2$. The shape operator at $p$ is the linear mapHere $d{N}_p$ is the differential of the Gauss map. Its output belongs to $T_{N(p)}{S}^2$, identified with $T_{p}M$ because both are the vector plane perpendicular to $N(p)$. The minus sign is our convention throughout this branch. ^[ref]

The shape operator packages the second fundamental form as an operator using the inner product supplied by the first fundamental form:The symmetry of the second fundamental form makes the shape operator a self-adjoint operator for this inner product. Thus an orthonormal basis of principal directions diagonalizes it. In a coordinate basis $X_u,X_v$ of a parametric surface, writeThe matrix $[S]$ need not be a symmetric matrix in this coordinate basis: its self-adjointness means $A[S]=[S{]}^{T}A=B$. Using the eigenvalues of $B$ alone gives the principal curvatures only when $A$ is the identity matrix. ^[ref]

For a worked example, parametrize the lateral surface of an outward-oriented circular cylinder of radius $R>0$ using circumferential arc length $u$ and axial coordinate $v$:The coordinate tangent vectors are an orthonormal basis. Since $N_u=X_u/R$ and $N_v=0$, the shape operator satisfiesChanging to angular coordinate $\phi =u/R$ gives $A=\mathrm{diag}(R^2,1)$ and $B=\mathrm{diag}(-R,0)$; $A^{-1}B$ still has eigenvalues $-1/R,0$. This is why the first fundamental form must participate in the calculation.

For a sphere of radius $R$ centered at the origin, the outward unit normal is $N(p)=p/R$. Hence $S_{p}v=-v/R$, so both principal curvatures are $-1/R$. For a plane, $dN=0$ and the shape operator vanishes. Reversing the unit normal negates the shape operator and both principal curvatures; the Gaussian curvature remains unchanged while the mean curvature changes sign, as in orientation reversal of surface curvature.

For two constraints $F(x,y,z)=G(x,y,z)=0$, the implicit function theorem gives differentiable functions $y(x),z(x)$ when the Jacobian determinant $\partial (F,G)/\partial (y,z)$ is nonzero. At its zeros this particular projection can fail even when a geometric curve remains. A nonzero tangent vector with zero $x$ component gives a vertical tangent in that projection. Singular constraint representations require further care: a determinant of zero alone does not prove an infinite derivative or even a regular intersection.

An ideal contact mechanical constraint supplying a normal reaction and no tangential friction. For a stationary constraining surface, the permitted velocity is a tangent vector, so the normal force has zero power and does no work. For a moving surface, normality alone does not imply zero power in the laboratory frame.

Spherical geometry studies lengths and angles measured along the surface of a sphere, rather than through its interior. It is a useful model whenever travel is constrained to a round surface. Its straightest paths are great circle arcs: they are geodesics, even though they appear curved in the surrounding Euclidean space. Two distinct great circles meet at an antipodal pair, so the parallel-line picture of Euclidean geometry does not apply.

For a sphere of radius $R>0$, write two positions as $Ru,Rv$, with $u,v$ unit vectors. Their shorter surface route has great-circle distancewhere the dot product gives the central angle in radians. The straight chord through the surrounding Euclidean space instead has length $2R\sin (d/(2R))$. For a quarter-circle separation these are $\pi R/2\approx 1.571R$ and $\sqrt{2}R\approx 1.414R$: specifying which distance we mean has a measurable effect. On the unit sphere, $R=1$ and great-circle distance equals the central angle.

Consider the spherical triangle with vertices $(R,0,0)$, $(0,R,0)$ and $(0,0,R)$, joined by the shorter great circle arcs. Each side has length $\pi R/2$. At each vertex the two outgoing tangent vectors point along the other coordinate axes and are orthogonal, so all three interior angles are $\pi /2$. This is the positive coordinate octant of the sphere, not a flat triangle drawn inside it.

For a convex minor spherical triangle with interior angles $\alpha ,\beta ,\gamma$ in radians, the spherical excess formula statesOur octant therefore has area $R^2(3\pi /2-\pi )=\pi R^2/2$. Independently, eight congruent octants partition the sphere, giving $4\pi R^2/8$ again. Thus angles alone determine area when $R$ is fixed, unlike Euclidean geometry, where similar triangles can have arbitrarily different areas. The Gauss-Bonnet theorem explains the connection: the sphere's constant Gaussian curvature $1/R^2$ times the area equals this excess. For spherical triangles whose area is small compared with $R^2$, the excess is small, recovering the familiar flat picture.

For a nondegenerate convex minor spherical triangle, the spherical law of cosines relates two angular sides $b,c$, their included interior angle $A$, and the opposite angular side $a$:Relabelling the vertices gives the corresponding formulas for $b,B$ and $c,C$. As with the spherical law of sines, physical side lengths must first be divided by the sphere's radius. The spherical law of cosines determines $a\in (0,\pi )$ uniquely from $b,c,A$ when the data describe such a spherical triangle. It also recovers each interior angle from three sides, for exampleUsing the cosine and the range $0<A<\pi$ selects a unique angle, avoiding the supplementary-angle ambiguity of an inverse sine calculation.

To derive the spherical law of cosines, let $u,v,w$ be the vertex unit vectors at $A,B,C$, so $u\cdot v=\cos c$, $u\cdot w=\cos b$ and $v\cdot w=\cos a$. The outgoing unit tangent vectors at $u$ are $t_B=(v-\cos cu)/\sin c$ and $t_C=(w-\cos bu)/\sin b$. Their dot product is both $\cos A$ and $(\cos a-\cos b\cos c)/(\sin b\sin c)$, proving the identity. This gives a geometric explanation: spherical trigonometry compares angles in the tangent plane using vectors in the surrounding Euclidean space.

For a worked example on the unit sphere, take the verticesAt $A$, the outgoing unit tangent vectors are $(1,0,0)$ and $(0,1,0)$, so the interior angle is a right angle. The vertex dot products give $b=c=\pi /3=60^{\circ }$. The spherical law of cosines then reduces to the right-angle identity $\cos a=\cos b\cos c$, givingTo find $B$, use the relabelled spherical law of cosines:The physical side lengths on a sphere of radius $R$ are ${\ell }_a\approx 1.318116R$ and ${\ell }_b={\ell }_c=\pi R/3\approx 1.047198R$. For example, on a toy sphere with $R=10$ metres they are approximately $13.18116$, $10.47198$ and $10.47198$ metres; the interior angles stay the same. Applying the planar Pythagorean theorem directly to the angular sides would instead give $\sqrt{60^2+60^2}\approx 84.852814^{\circ }$. That flat right triangle calculation does not give the great-circle distance between $B$ and $C$.

The familiar law of cosines reappears as a local approximation. With $a,b,c$ all small in radians, substitute the small-angle approximations $\sin x\approx x$ and $\cos x\approx 1-x^2/2$ into the spherical law of cosines and keep terms through second order:Multiplication by $R^2$ gives the planar law of cosines for the physical side lengths. For a fixed small surface triangle on a progressively larger sphere, the angular sides shrink and this approximation improves. Thus the curved and flat formulas fit together: the Pythagorean theorem is recovered locally when $A$ is a right angle.

For a nondegenerate convex minor spherical triangle, with angular sides $a,b,c$ opposite interior angles $A,B,C$, the spherical law of sines statesThe sines of sides refer to central angles, so a physical side length ${\ell }_a$ on a sphere of radius $R$ must enter as $\sin ({\ell }_a/R)$. All six angles lie in $(0,\pi )$, and the denominators are positive. The reciprocal form, $\sin a/\sin A=\sin b/\sin B=\sin c/\sin C$, is equivalent. This is Tatum's sine formula, equation 3.5.1, in Celestial Mechanics, section 3.5: spherical sine and cosine formulas.

A short proof explains the symmetry of the spherical law of sines. Let $u,v,w$ be the vertex unit vectors at $A,B,C$. At $u$, the outgoing unit tangent vectors towards $v,w$ are $t_B=(v-\cos cu)/\sin c$ and $t_C=(w-\cos bu)/\sin b$. Their angle is $A$, so the magnitude of their cross product is $\sin A$. Since both tangent vectors are orthogonal to $u$, the scalar triple product obeysDividing by $\sin a\sin b\sin c$ proves the spherical law of sines. Its beauty is that the same three-dimensional scalar triple product can be calculated at any vertex, producing all three ratios at once.

For the spherical triangle worked out under the spherical law of cosines, $A=90^{\circ }$, $b=c=60^{\circ }$ and $\cos a=1/4$. Hence $\sin a=\sqrt{15}/4$, and the spherical law of sines givesThe inverse sine gives $B\approx 63.434949^{\circ }$, but the supplementary angle $116.565051^{\circ }$ has the same sine. The spherical law of sines alone therefore leaves two candidates for this angle. The spherical law of cosines, using the known three sides, gives $\cos B=1/\sqrt{5}>0$ and selects $63.434949^{\circ }$; likewise $C=B$. When only two sides and an opposite angle are given, two distinct valid spherical triangles can remain. A cosine relation resolves an inverse sine branch only when enough additional data are known; it does not make every ambiguous data set unique.

A spherical triangle is bounded by three great circle arcs. For the usual convex minor spherical triangle, the vertices are neither antipodal nor on one great circle, and its angular sides and interior angles lie in $(0,\pi )$. An interior angle is measured between outgoing tangent vectors to the two sides at a vertex, rather than between straight chords inside the sphere. The spherical law of sines and spherical law of cosines determine its side-angle relations; the spherical excess formula determines its area.

A subdivision surface is the limit of repeated mesh refinement and vertex averaging. Exact interrogation evaluates this limit and, where defined, its tangent vectors and normal vectors; a finite refined control mesh is generally only an approximation to the limit.

// Target: analysis.bigb

For the radial profile $(\rho (t),z(t))$, with $\rho >0$, use the parametrized surface $\mathbf{R}(t,\theta )=(\rho \cos \theta ,\rho \sin \theta ,z)$, $0\le \theta <2\pi$. The two tangent vectors are orthogonal and have Euclidean norms $\sqrt{{\rho }^{\prime 2}+z^{\prime 2}}$ and $\rho$. Their cross product gives the displayed surface area. Geometric surface area assumes a single covering of the radial profile. A general three-dimensional generating curve must first be replaced by its actual radius $\rho (t)=\sqrt{x(t{)}^2+y(t{)}^2}$; its spatial arc length cannot simply replace the radial profile's arc length.

A tangential boundary derivative differentiates the boundary restriction of a function along a tangent vector. Prescribing it fixes the boundary value only up to a constant on each connected boundary component. This differs from a Neumann boundary condition, which specifies a normal derivative.

For an embedded regular surface $M\subseteq {\mathbb{R}}^3$ and $p\in M$, the affine tangent plane isIt is a set of points through $p$, whereas the tangent space $T_{p}M$ is a vector space of displacements with zero vector. Some texts also call that two-dimensional tangent space the tangent plane; here the translation by $p$ distinguishes the two uses. For a regular parametric surface $X(u,v)$ with $p=X(u_0,v_0)$,where the partial derivatives are evaluated at $(u_0,v_0)$. Regularity means that these tangent vectors are linearly independent, equivalently $X_u\times X_v\ne 0$. Their cross product is a normal vector. ^[ref]

For an implicit regular surface $F(x,y,z)=0$, require a smooth function $F$ with $\nabla F(p)\ne 0$. ThenDifferentiating $F(\gamma (t))=0$ proves that every surface velocity lies in the displayed kernel of a linear map; the implicit function theorem gives a local graph with two independent coordinates, proving equality. The gradient is the normal vector with respect to the ambient dot product. ^[ref]

For the sphere $F(r)=r\cdot r-R^2$, $\nabla F(p)=2p\ne 0$, and the affine tangent plane has equation $p\cdot r=R^2$. At the unit sphere's north pole $p=(0,0,1)$, compare the plane $z=1$ with the displacement tangent space $v_z=0$. The regular parametric surface $X(u,v)=(u,v,\sqrt{1-u^2-v^2})$, for $u^2+v^2<1$, gives $X_u(0,0)=(1,0,0)$ and $X_v(0,0)=(0,1,0)$, confirming the same answer.

This analytical illustration shows a unit sphere and a vertical cross-section. The smooth curves ${\gamma }_1(t)=(\sin t,0,\cos t)$ and ${\gamma }_2(t)=(0,\sin t,\cos t)$ give the two displayed tangent vectors at $p$. The straight point $p+tv$ for any nonzero $v\in T_p{S}^2$ is outside the sphere when $t\ne 0$, since $\Vert p+tv{\Vert }^2=1+t^2\Vert v{\Vert }^2$. A tangent plane therefore captures first-order motion, not an exact patch of the regular surface.

The same idea makes local approximation useful. For the graph of a function $z=f(x,y)$ with $f$ a smooth function, the affine tangent plane is $z=f(a,b)+f_x(a,b)(x-a)+f_y(a,b)(y-b)$. ^[ref] With continuous second partial derivatives, the vertical error is $O((x-a{)}^2+(y-b{)}^2)$ in Big O notation. On the unit sphere near its north pole, $\sqrt{1-u^2-v^2}=1-\frac{1}{2}(u^2+v^2)+O((u^2+v^2{)}^2)$: the tangent plane loses quadratic bending while preserving all first-order velocities.

A failed regularity test needs interpretation. Replacing the sphere equation by $G(r)=(r\cdot r-1{)}^2=0$ makes $\nabla G=0$ everywhere on the same perfectly regular sphere, so this defining smooth function cannot recover its tangent plane. Similarly, the usual spherical parametric surface loses matrix rank at the poles, while the graph parametrization above remains regular there. A vanishing gradient or degenerate parametrization invalidates that calculation; it alone does not prove that the underlying geometry has no tangent plane.

At a fixed point $p$ of an $n$-dimensional smooth manifold $M$ without boundary, the tangent space $T_{p}M$ is the real vector space of all tangent vectors at $p$. In manifold charts $x=(x^1,\dots ,x^n)$ and $y=(y^1,\dots ,y^n)$, a tangent vector has componentsThe coordinate basis transforms inversely, so the tangent vector is unchanged. This invertible Jacobian matrix also preserves vector addition and scalar multiplication, giving a vector space of dimension $n$ independently of coordinates. ^[ref]

For example, at $(x,y)=(1,0)$ in ${\mathbb{R}}^2$, use the alternative manifold chart $u=x$, $w=y+x^2$. The Jacobian matrix sends $(v_x,v_y)=(2,3)$ to $(v_u,v_w)=(2,7)$. Indeed, the smooth curve $(x(t),y(t))=(1+2t,3t)$ gives $w(t)=1+7t+4t^2$. Conversely, ${\partial }_u{|}_p={\partial }_x{|}_p-2{\partial }_y{|}_p$ and ${\partial }_w{|}_p={\partial }_y{|}_p$, so $2{\partial }_u{|}_p+7{\partial }_w{|}_p=2{\partial }_x{|}_p+3{\partial }_y{|}_p$. The numbers change; their action on every smooth function does not.

For an embedded submanifold $M\subseteq {\mathbb{R}}^N$, the tangent space can be identified with the vector subspaceA proof outline explains why every coordinate vector is possible: in a manifold chart $x$, lift the short straight path $x(p)+ta$ using $x^{-1}$. Its tangent vector has components $a$. For a regular local parametrization $X:U\subseteq {\mathbb{R}}^n\to M$, the same chain rule gives $T_{p}M=\mathrm{im}D{X}_q$, where $X(q)=p$. This is the tangent space from a local parametrization.

For the sphere $S_R^2=\{p\in {\mathbb{R}}^3:p\cdot p=R^2\}$, differentiating that equation gives $p\cdot v=0$. Conversely, for any nonzero $v$ with $p\cdot v=0$, the smooth curvelies on the sphere and has ${\gamma }_v^{\prime }(0)=v$; use a constant smooth curve for $v=0$. Thus $T_p{S}_R^2=\{v:p\cdot v=0\}$, not merely a subset of that plane. At $p=(0,0,R)$ this vector subspace is $v_z=0$, while the affine tangent plane drawn through $p$ is $z=R$.

This also explains a mechanical constraint without requiring a force calculation. In the ideal mathematical model of a particle confined to a stationary sphere, its allowed velocities at $p$ are exactly $T_p{S}_R^2$: the radial component must vanish, but either independent tangential direction is available. For a general stationary holonomic constraint $F(q)=0$ with independent equations, the regular level set theorem gives $T_{q}M=\ker D{F}_q$. These are instantaneous kinematic possibilities; they do not determine the actual subsequent motion. On a general smooth manifold, tangent spaces at different base points cannot be identified canonically; an affine connection supplies parallel transport when such comparisons are needed.

A tangent vector at a point $p$ of a smooth manifold without boundary records an instantaneous velocity through $p$. Its base point is part of the data; measuring its length additionally requires a Riemannian metric. Intrinsically, a smooth curve $\gamma$ with $\gamma (0)=p$ defines its tangent vector by how every local smooth function changes along it:This is a derivation at a point. Two smooth curves through $p$ represent the same tangent vector if these derivatives agree for every such smooth function; equivalently, their coordinate velocity vectors agree in one manifold chart. The chain rule makes that agreement independent of the chosen manifold chart. ^[ref]

For an embedded submanifold of Euclidean space, the intrinsic tangent vector is represented by the ordinary derivative ${\gamma }^{\prime }(0)$, with its base point $p$ retained. A drawn arrow starting at $p$ represents that tangent vector; its endpoint $p+{\gamma }^{\prime }(0)$ need not lie on the embedded submanifold. Nor is a tangent vector required to be a unit tangent vector: a stationary smooth curve gives the zero vector, which belongs to the tangent space but has no direction to normalize.

Here are explicit calculations. In ${\mathbb{R}}^2$, $\gamma (t)=(t,t^2)$ and $\eta (t)=(t,2t^2)$ have different routes but the same tangent vector $(1,0)$ at $(0,0)$. For $f(x,y)=3x+y$, both give $v(f)=3$. For a circle of radius $R$, $c(t)=(R\cos (\omega t),R\sin (\omega t))$ has tangent vector $(0,R\omega )$ at $(R,0)$: changing $\omega$ changes the speed or reverses travel without changing the circle.

A reparametrization $\beta (s)=\gamma (h(s))$ with $h(0)=0$ gives, by the chain rule, ${\beta }^{\prime }(0)=h^{\prime }(0){\gamma }^{\prime }(0)$. A regular change of curve parameter has $h^{\prime }(0)\ne 0$. The singular choice $h(s)=s^3$ instead givesBoth describe the same parabola, whose tangent line at the origin is still the horizontal axis. The zero velocity diagnoses this singular change of curve parameter, whose inverse is not differentiable at zero; it does not show that the geometric tangent line is absent. The regular curve $\gamma$ supplies its nonzero direction.

The diagram is an analytical illustration of these exact parametric curves, with equal parameter increments marked on each route. It separates geometry from traversal: the tangent line depends on the route, while a velocity vector also depends on how the curve parameter advances.

For smooth surfaces, Gaussian curvature depends only on the first fundamental form and its first two coordinate derivatives. Consequently any smooth local isometry $f:M\to N$ satisfies $K_N(f(p))=K_M(p)$. The surprise is that $K$ is defined using the shape operator, which describes bending in space, yet its value can be recovered from measurements made entirely within the surface.

In his 1827 paper, Carl Friedrich Gauss expressed Gaussian curvature through the metric coefficients and their partial derivatives in Article 11, then stated its invariance under development in Article 12. Article 13 explicitly compares a plane with surfaces developable onto it, including cylindrical surfaces. The paper was presented on 8 October 1827 and printed in 1828; the Morehead–Hiltebeitel English translation appeared in 1902. ^[ref]

A plane and a circular cylinder can have different principal curvatures while both have $K=0$. Their matching first fundamental forms exhibit the actual local isometry. In contrast, no nonempty open patch of a sphere of radius $R$ can be flattened smoothly into a plane while preserving its first fundamental form: $1/R^2$ cannot become zero. Making the patch small reduces practical distortion, but never produces an exact local isometry. This is a geometric obstruction to a perfectly distance-preserving flat map of a spherical world, independent of where one places a seam.

The following proof outline connects the two descriptions of Gaussian curvature. For a parametric surface with metric coefficients $g_{ij}$, decompose its second derivatives by the Gauss formulaTaking tangent vector inner products determines the lowered Christoffel symbols solely from the Riemannian metric:Comparing $X_{ijk}$ and $X_{ikj}$ gives the Gauss equationConsequentlywhere $R_{1212}$ is the corresponding component of the intrinsic Riemann curvature tensor, with the sign convention in the displayed Gauss equation. The final expression uses only $g_{ij}$ and its first two derivatives.

The timelike cone is the open set of nonzero tangent vectors with timelike squared norm. It has future and past components once a time orientation is chosen. Homogeneity extends an identity for unit timelike vectors to this whole open set. In particular, if a symmetric quadratic form vanishes on every timelike vector, its polynomial vanishes on an open set and therefore vanishes identically. This observation proves the converse from geodesic first integrals to the Killing equation.

A timelike curve has a nonzero tangent vector with negative Minkowski norm in metric signature $(-,+,+,+)$ at every regular point. A physical massive observer's world line is also future-directed. Parametrizing it by coordinate time gives $-c^2+|\mathbf{v}{|}^2<0$, or $|\mathbf{v}|<c$: the local tangent condition is exactly the subluminal speed condition. It is stronger than checking only that two endpoints have timelike separation.

The tractoid is the surface of revolution of a unit tractrix. For $X(\theta ,v)=(h(\theta ),\sin \theta \cos v,\sin \theta \sin v)$ with $0<\theta <\pi /2$, its coordinate tangent vectors are orthogonal and have squared lengths ${\cot }^2\theta ,{\sin }^2\theta$. The coordinates $x=v$, $y=\csc \theta$ transform this first fundamental form into $(d{x}^2+d{y}^2)/y^2$, giving a local isometry to the hyperbolic plane. Thus its Gaussian curvature is minus one. The angular coordinate is periodic globally, while its unwrapped coordinate is a real variable; this distinguishes the local metric identification from a global embedding of the whole hyperbolic plane.

For a parametric surface $F(u,v)$ and a scalar implicit equation $f(P)=0$, set $g=f\circ F$. At a regular transversal intersection, $\nabla g\ne 0$, and $(-g_v,g_u)$ is a parameter tangent vector. Normalize its image under $DF$ to choose physical arc length steps. A predictor followed by Newton method correction of $g=0$ and a transverse plane equation gives a local continuation method. Certified patch bounds are useful for finding seeds; corner signs alone do not detect a closed intersection entirely inside a patch.

For a regular curve $\mathbf{r}(t)$ with nonzero velocity vector, its unit tangent vector isThe Euclidean norm in the denominator removes the traversal speed, leaving a unit vector that records direction. A positive change of traversal speed leaves this direction unchanged; reversing the direction of travel reverses $T$. In arc length coordinates, $T=d\mathbf{r}/ds$.

For the parabola $\mathbf{r}(t)=(t,t^2)$, the velocity vector is $(1,2t)$ and the speed is $\sqrt{1+4t^2}$. At $t=1$, $T=(1,2{)}^T/\sqrt{5}$. The tangent vector $(1,2{)}^T$ and the unit tangent vector point in the same direction, but only the latter has Euclidean norm one. The principal normal vector is a further choice based on how this direction changes, where the curvature is nonzero. ^[ref]

In elliptic cylindrical coordinates, write a vector field as $A=A_u{e}_u+A_v{e}_v+A_z{e}_z$ using the local orthonormal unit vectors. The scale factors of orthogonal coordinates are $(h,h,1)$, so the divergence and curl areApplying the divergence to the gradient of a scalar function gives $\Delta \varphi =({\varphi }_{uu}+{\varphi }_{vv})/h^2+{\varphi }_{zz}$. These formulas hold where $h>0$. They use physical components along unit vectors, rather than components along the unnormalized coordinate tangent vectors.