Radius of curvature and center of curvature

[1] In differential geometry, the radius of curvature, R, is the reciprocal of the curvature. For a curve, it equals the radius of the circular arc which best approximates the curve at that point. For surfaces, the radius of curvature is the radius of a circle that best fits a normal section or combinations thereof.

## Definition

In the case of a space curve, the radius of curvature is the length of the curvature vector.

In the case of a plane curve, then R is the absolute value of

${\displaystyle {\frac {ds}{d\varphi }}={\frac {1}{\kappa }},}$

where s is the arc length from a fixed point on the curve, φ is the tangential angle and ${\displaystyle \scriptstyle \kappa }$ is the curvature.

If the curve is given in Cartesian coordinates as y(x), then the radius of curvature is (assuming the curve is differentiable up to order 2):

${\displaystyle R=\left|{\frac {\left(1+y'^{\,2}\right)^{3/2}}{y''}}\right|,\qquad {\mbox{where}}\quad y'={\frac {dy}{dx}},\quad y''={\frac {d^{2}y}{dx^{2}}},}$

and | z | denotes the absolute value of z.

If the curve is given parametrically by functions x(t) and y(t), then the radius of curvature is

${\displaystyle R=\;\left|{\frac {ds}{d\varphi }}\right|\;=\;\left|{\frac {{\big (}{{\dot {x}}^{2}+{\dot {y}}^{2}}{\big )}^{3/2}}{{\dot {x}}{\ddot {y}}-{\dot {y}}{\ddot {x}}}}\right|,\qquad {\mbox{where}}\quad {\dot {x}}={\frac {dx}{dt}},\quad {\ddot {x}}={\frac {d^{2}x}{dt^{2}}},\quad {\dot {y}}={\frac {dy}{dt}},\quad {\ddot {y}}={\frac {d^{2}y}{dt^{2}}}.}$

Heuristically, this result can be interpreted as

${\displaystyle R={\frac {\left|\mathbf {v} \right|^{3}}{\left|\mathbf {v} \times \mathbf {\dot {v}} \right|}},\qquad {\mbox{where}}\quad \left|\mathbf {v} \right|=\left|({\dot {x}},{\dot {y}})\right|=R{\frac {d\varphi }{dt}}.}$

## Formula

If ${\displaystyle \gamma :\mathbb {R} \rightarrow \mathbb {R} ^{n}}$ is a parameterized curve in ${\displaystyle \mathbb {R} ^{n}}$ then the radius of curvature at each point of the curve, ${\displaystyle \rho :\mathbb {R} \rightarrow \mathbb {R} }$, is given by

${\displaystyle \rho ={\frac {|\gamma '|^{3}}{\sqrt {|\gamma '|^{2}\;|\gamma ''|^{2}-(\gamma '\cdot \gamma '')^{2}}}}}$.

As a special case, if f(t) is a function from ${\displaystyle \mathbb {R} }$ to ${\displaystyle \mathbb {R} }$, then the curvature of its graph, ${\displaystyle \gamma (t)=(t,f(t))}$, is

${\displaystyle \rho (t)={\frac {|1+f'^{2}(t)|^{3/2}}{|f''(t)|}}.}$

### Derivation

Let ${\displaystyle \gamma }$ be as above, and fix ${\displaystyle t}$. We want to find the radius ${\displaystyle \rho }$ of a parameterized circle which matches ${\displaystyle \gamma }$ in its zeroth, first, and second derivatives at ${\displaystyle t}$. Clearly the radius will not depend on the position (${\displaystyle \gamma (t)}$), only on the velocity (${\displaystyle \gamma '(t)}$) and acceleration (${\displaystyle \gamma ''(t)}$). There are only three independent scalars that can be obtained from two vectors v and w, namely v·v, v·w, and w·w. Thus the radius of curvature must be a function of the three scalars ${\displaystyle |\gamma '^{2}(t)|\,\!}$, ${\displaystyle |\gamma ''^{2}(t)|\,\!}$ and ${\displaystyle \gamma '(t)\cdot \gamma ''(t)}$.

The general equation for a parameterized circle in ${\displaystyle \mathbb {R} ^{n}}$ is

${\displaystyle g(u)=A\cos(h(u))+B\sin(h(u))+C\,\!}$

where ${\displaystyle C\in \mathbb {R} ^{n}}$ is the center of the circle (irrelevant since it disappears in the derivatives), ${\displaystyle A,B\in \mathbb {R} ^{n}}$ are perpendicular vectors of length ${\displaystyle \rho }$ (that is, ${\displaystyle A\cdot A=B\cdot B=\rho ^{2}}$ and ${\displaystyle A\cdot B=0}$), and ${\displaystyle h:\mathbb {R} \rightarrow \mathbb {R} }$ is an arbitrary function which is twice differentiable at t.

The relevant derivatives of g work out to be

${\displaystyle {\begin{array}{lll}|g'|^{2}&=&\rho ^{2}(h')^{2}\\g'\cdot g''&=&\rho ^{2}h'h''\\|g''|^{2}&=&\rho ^{2}\left((h')^{4}+(h'')^{2}\right)\end{array}}}$

If we now equate these derivatives of g to the corresponding derivatives of ${\displaystyle \gamma }$ at t we obtain

${\displaystyle {\begin{array}{lll}|\gamma '^{2}(t)|&=&\rho ^{2}h'^{2}(t)\\\gamma '(t)\cdot \gamma ''(t)&=&\rho ^{2}h'(t)h''(t)\\|\gamma ''^{2}(t)|&=&\rho ^{2}(h'^{4}(t)+h''^{2}(t))\end{array}}}$

These three equations in three unknowns (${\displaystyle \rho }$, ${\displaystyle h'(t)}$ and ${\displaystyle h''(t)}$) can be solved for ${\displaystyle \rho }$, giving the formula for the radius of curvature:

${\displaystyle \rho (t)={\frac {|\gamma '^{3}(t)|}{\sqrt {|\gamma '^{2}(t)|\;|\gamma ''^{2}(t)|-(\gamma '(t)\cdot \gamma ''(t))^{2}}}}}$

or, omitting the parameter (t) for readability,

${\displaystyle \rho ={\frac {|\gamma '|^{3}}{\sqrt {|\gamma '|^{2}\;|\gamma ''|^{2}-(\gamma '\cdot \gamma '')^{2}}}}}$.

## Examples

### Semicircles and circles

For a semi-circle of radius a in the upper half-plane

${\displaystyle y={\sqrt {a^{2}-x^{2}}},\quad y'={\frac {-x}{\sqrt {a^{2}-x^{2}}}},\quad y''={\frac {-a^{2}}{(a^{2}-x^{2})^{3/2}}},\quad R=|-a|=a.}$
An ellipse (red) and its evolute (blue). The dots are the vertices of the ellipse, at the points of greatest and least curvature.

For a semi-circle of radius a in the lower half-plane

${\displaystyle y=-{\sqrt {a^{2}-x^{2}}},\quad R=|a|=a.}$

The circle of radius a has a radius of curvature equal to a.

### Ellipses

In an ellipse with major axis 2a and minor axis 2b, the vertices on the major axis have the smallest radius of curvature of any points ${\displaystyle \left(R={\frac {b^{2}}{a}}\right)}$, and the vertices on the minor axis have the largest radius of curvature of any points ${\displaystyle \left(R={\frac {a^{2}}{b}}\right)}$.

## Applications

### Stress in semiconductor structures

Stress in the semiconductor structure involving evaporated thin films usually results from the thermal expansion (thermal stress) during the manufacturing process. Thermal stress occurs because film depositions are usually made above room temperature. Upon cooling from the deposition temperature to room temperature, the difference in the thermal expansion coefficients of the substrate and the film cause thermal stress.[2]

Intrinsic stress results from the microstructure created in the film as atoms are deposited on the substrate. Tensile stress results from microvoids in the thin film, because of the attractive interaction of atoms across the voids.

The stress in thin film semiconductor structures results in the buckling of the wafers. The radius of the curvature of the stressed structure is related to stress tensor in the structure, and can be described by modified Stoney formula.[3] The topography of the stressed structure including radii of curvature can be measured using optical scanner methods. The modern scanner tools have capability to measure full topography of the substrate and to measure both principal radii of curvature, while providing the accuracy of the order of 0.1% for radii of curvature of 90 m and more.[4]