Solid angle

Not to be confused with spherical angle.

In geometry, a solid angle (symbol: Ω) is the two-dimensional angle in three-dimensional space that an object subtends at a point. It is a measure of how large the object appears to an observer looking from that point. In the International System of Units (SI), a solid angle is expressed in a dimensionless unit called a steradian (symbol: sr).

A small object nearby may subtend the same solid angle as a larger object farther away. For example, although the Moon is much smaller than the Sun, it is also much closer to Earth. Indeed, as viewed from any point on Earth, both objects have approximately the same solid angle as well as apparent size. This is evident during a solar eclipse.

Definition and properties

An object's solid angle in steradians is equal to the area of the segment of a unit sphere, centered at the angle's vertex, that the object covers. A solid angle in steradians equals the area of a segment of a unit sphere in the same way a planar angle in radians equals the length of an arc of a unit circle. Solid angles are often used in physics, in particular astrophysics. The solid angle of an object that is very far away is roughly proportional to the ratio of area to squared distance. Here "area" means the area of the object when projected along the viewing direction.

Any area on a sphere which is equal in area to the square of its radius, when observed from its center, subtends precisely one steradian.

The solid angle of a sphere measured from a point in its interior is 4π sr, and the solid angle subtended at the center of a cube by one of its faces is one-sixth of that, or 2π/3 sr. Solid angles can also be measured in square degrees (1 sr = (180/π)² square degree), in square minutes and square seconds, or in fractions of the sphere (1 sr = 1/4π fractional area), also known as spat (1 sp = 4π sr).

In spherical coordinates there is a simple formula for the differential,

d\Omega = \sin\theta\,d\theta\,d\varphi

where $\theta$ is the colatitude (angle from the North pole) and $\varphi$ is the longitude.

The solid angle for an arbitrary oriented surface S subtended at a point P is equal to the solid angle of the projection of the surface S to the unit sphere with center P, which can be calculated as the surface integral:

\Omega = \iint\limits_S \frac{ \hat{r} \cdot \hat{n} \,d\Sigma }{r^2}\ = \iint\limits_S \sin\theta\,d\theta\,d\varphi

where $\hat{r} = \vec{r} / r$ is the unit vector corresponding to $\vec{r}$ , the position vector of an infinitesimal area of surface $\, d\Sigma$ with respect to point P, and where $\hat{n}$ represents the unit normal vector to $\, d\Sigma$ . Even if the projection on the unit sphere to the surface S is not isomorphic, the multiple folds are correctly considered according to the surface orientation described by the sign of the scalar product $\hat{r} \cdot \hat{n}$ .

Thus one can approximate the solid angle subtended by a small facet having flat surface area dΣ, orientation $\hat{n}$ , and distance r from the viewer as:

d\Omega = 4 \pi (d\Sigma / A) \, (\hat{r} \cdot \hat{n})

where the surface area of a sphere is $A=4\pi r^2$ .

Practical applications

Defining luminous intensity and luminance, and the correspondent radiometric quantities radiant intensity and radiance.
Calculating spherical excess E of a spherical triangle
The calculation of potentials by using the boundary element method (BEM)
Evaluating the size of ligands in metal complexes, see ligand cone angle.
Calculating the electric field and magnetic field strength around charge distributions.
Deriving Gauss's Law.
Calculating emissive power and irradiation in heat transfer.
Calculating cross sections in Rutherford scattering.
Calculating cross sections in Raman scattering.
The solid angle of the acceptance cone of the optical fiber

Solid angles for common objects

Cone, spherical cap, hemisphere

Section of cone (1) and spherical cap (2) inside a sphere. In this figure θ = A/2 and r = 1.

The solid angle of a cone with apex angle 2θ, is the area of a spherical cap on a unit sphere

\Omega = 2\pi \left (1 - \cos {\theta} \right)\

For small θ such that sin θ ≈ θ, this reduces to the area of a circle πθ².

The above is found by computing the following double integral using the unit surface element in spherical coordinates:

\int_0^{2\pi} \int_0^{\theta} \sin\theta' \, d \theta' \, d \phi = 2\pi\int_0^{\theta} \sin\theta' \, d \theta' \ = 2\pi\left[ -\cos\theta' \right]_0^{\theta}\ = 2\pi\left(1 - \cos\theta \right)\

This formula can also be derived without the use of calculus. Over 2200 years ago Archimedes proved that the surface area of a spherical cap is always equal to the area of a circle whose radius equals the distance from the rim of the spherical cap to the point where the cap's axis of symmetry intersects the cap.^[1] In the diagram this radius is given as:

2r \sin \left( \frac{ \theta}{2} \right)

Hence for a unit sphere the solid angle of the spherical cap is given as:

\Omega = 4\pi \sin^2 \left( \frac{\theta}{2} \right) = 2\pi \left (1 - \cos {\theta} \right)\

When θ = π/2 , the spherical cap becomes a hemisphere having a solid angle 2π.

The solid angle of the complement of the cone (picture a melon with the cone cut out) is clearly:

4\pi - \Omega = 2\pi \left(1 + \cos {\theta} \right)

This is also the solid angle of the part of the celestial sphere that a Terran astronomical observer positioned at latitude θ can see as the earth rotates. At the equator you see all of the celestial sphere, at either pole only one half.

The solid angle subtended by a segment of a spherical cap cut by a plane at angle $\gamma$ from the cone's axis and passing through the cone's apex can be calculated by the formula:^[2]

\Omega = 2 \left[ \arccos \left(\frac{\sin\gamma}{\sin\theta}\right) - \cos\theta \arccos\left(\frac{\tan\gamma}{\tan\theta}\right) \right]\

Tetrahedron

Let OABC be the vertices of a tetrahedron with an origin at O subtended by the triangular face ABC where $\vec a\ ,\, \vec b\ ,\, \vec c$ are the vector positions of the vertices A, B and C. Define the vertex angle $\theta_a \,$ to be the angle BOC and define $\theta_b ,\, \theta_c$ correspondingly. Let $\phi_{ab} \,$ be the dihedral angle between the planes that contain the tetrahedral faces OAC and OBC and define $\phi_{bc} ,\, \phi_{ac}$ correspondingly. The solid angle $\Omega$ subtended by the triangular surface ABC is given by

\Omega = \left(\phi_{ab} + \phi_{bc} + \phi_{ac}\right)\ - \pi\

This follows from the theory of spherical excess and it leads to the fact that there is an analogous theorem to the theorem that "The sum of internal angles of a planar triangle is equal to $\pi$ ", for the sum of the four internal solid angles of a tetrahedron as follows:

\sum_{i=1}^4 \Omega_i = 2 \sum_{i=1}^6 \phi_i\ - 4 \pi\

where $\phi_i \,$ ranges over all six of the dihedral angles between any two planes that contain the tetrahedral faces OAB, OAC, OBC and ABC.

An efficient algorithm for calculating the solid angle $\Omega$ subtended by the triangular surface ABC where $\vec a\ ,\, \vec b\ ,\, \vec c$ are the vector positions of the vertices A, B and C has been given by Oosterom and Strackee:^[3]

\tan \left( \frac{1}{2} \Omega \right) = \frac{\left|\vec a\ \vec b\ \vec c\right|}{abc + \left(\vec a \cdot \vec b\right)c + \left(\vec a \cdot \vec c\right)b + \left(\vec b \cdot \vec c\right)a}

where

\left|\vec a\ \vec b\ \vec c\right|=\vec a \cdot (\vec b \times \vec c)

denotes the scalar triple product of the three vectors;

\vec a

is the vector representation of point A, while

\, a

is the magnitude of that vector (the origin-point distance)

\vec a \cdot \vec b

denotes the scalar product.

When implementing the above equation care must be taken with the atan function to avoid negative or incorrect solid angles. One source of potential errors is that the scalar triple product can be negative if a,b,c have the wrong winding. Computing abs(det) is a sufficient solution since no other portion of the equation depends on the winding. The other pitfall arises when the scalar triple product is positive but the divisor is negative. In this case atan returns a negative value that must be biased by $\pi$ .

Another useful formula for calculating the solid angle of the tetrahedron at the origin O that is purely a function of the vertex angles $\theta_a ,\, \theta_b ,\, \theta_c$ is given by L'Huilier's theorem^[4]^[5] as

\tan \left( \frac{1}{4} \Omega \right) = \sqrt{ \tan \left( \frac{\theta_s}{2}\right) \tan \left( \frac{\theta_s - \theta_a}{2}\right) \tan \left( \frac{\theta_s - \theta_b}{2}\right) \tan \left(\frac{\theta_s - \theta_c}{2}\right)}

where

\theta_s = \frac {\theta_a + \theta_b + \theta_c}{2}

Pyramid

The solid angle of a four-sided right rectangular pyramid with apex angles $\, a$ and $\, b$ (dihedral angles measured to the opposite side faces of the pyramid) is

\Omega = 4 \arcsin \left( \sin \left({a \over 2}\right) \sin \left({b \over 2}\right) \right)

If both the side lengths (α and β) of the base of the pyramid and the distance (d) from the center of the base rectangle to the apex of the pyramid (the center of the sphere) are known, then the above equation can be manipulated to give

\Omega = 4 \arctan \frac {\alpha\beta} {2d\sqrt{4d^2 + \alpha^2 + \beta^2}}

The solid angle of a right n-gonal pyramid, where the pyramid base is a regular n-sided polygon of circumradius (r), with a pyramid height (h) is

\Omega = 2\pi - 2n \arctan\left(\frac {\tan \left({\pi\over n}\right)}{\sqrt{1 + {r^2 \over h^2}}} \right)

The solid angle of an arbitrary pyramid with an n-sided base defined by the sequence of unit vectors representing edges $\{ s_1, s_2, ..., s_n \}$ can be efficiently computed by:^[2]

\Omega = 2\pi - \arg \prod_{j=1}^{n} \left( \left( s_{j-1} s_j \right)\left( s_{j} s_{j+1} \right) - \left( s_{j-1} s_{j+1} \right) + i\left[ s_{j-1} s_j s_{j+1} \right] \right)

where parentheses (* *) is a scalar product and square brackets [* * *] is a scalar triple product, and $i$ is an imaginary unit. Indices are cycled: $s_0 = s_n$ and $s_{n+1} = s_1$ .

Latitude-longitude rectangle

The solid angle of a latitude-longitude rectangle on a globe is $\left ( \sin \phi_N - \sin \phi_S \right ) \left ( \theta_E - \theta_W \,\! \right)\mathrm{sr}$ , where $\phi_N \,\!$ and $\phi_S \,\!$ are north and south lines of latitude (measured from the equator in radians with angle increasing northward), and $\theta_E \,\!$ and $\theta_W \,\!$ are east and west lines of longitude (where the angle in radians increases eastward).^[6] Mathematically, this represents an arc of angle $\phi_N - \phi_S \,\!$ swept around a sphere by $\theta_E - \theta_W \,\!$ radians. When longitude spans 2π radians and latitude spans π radians, the solid angle is that of a sphere.

A latitude-longitude rectangle should not be confused with the solid angle of a rectangular pyramid. All four sides of a rectangular pyramid intersect the sphere's surface in great circle arcs. With a latitude-longitude rectangle, only lines of longitude are great circle arcs; lines of latitude are not.

Sun and Moon

The Sun is seen from Earth at an average angular diameter of about 9.35×10⁻³ radians. The Moon is seen from Earth at an average diameter of 9.22×10⁻³ radians. We can substitute these into the equation given above for the solid angle subtended by a cone with apex angle $2 \theta \,\!$ :

\Omega = 2 \pi \left (1 - \cos {\theta} \right)

The resulting value for the Sun is 6.87×10⁻⁵ steradians. The resulting value for the Moon is 6.67×10⁻⁵ steradians. In terms of the total celestial sphere, the Sun and the Moon subtend fractional areas of 0.000546% (Sun) and 0.000531% (Moon). On average, the Sun is larger in the sky than the Moon even though it is much, much farther away.

Solid angles in arbitrary dimensions

The solid angle subtended by the complete (d-1)-dimensional spherical surface of the unit sphere in d-dimensional Euclidean space can be defined in any number of dimensions d. One often needs this solid angle factor in calculations with spherical symmetry. It is given by the formula

\Omega_{d} = \frac{2\pi^\frac{d}{2}}{\Gamma\left(\frac{d}{2}\right)}\

where $\Gamma$ is the Gamma function. When d is an integer, the Gamma function can be computed explicitly.^[7] It follows that

\Omega_{d} = \begin{cases} \frac{1}{ \left(\frac{d}{2} - 1 \right)!} 2\pi^\frac{d}{2}\ & d\text{ even} \\ \frac{\left(\frac{1}{2}\left(d - 1\right)\right)!}{(d - 1)!} 2^d \pi^{\frac{1}{2}(d - 1)}\ & d\text{ odd} \end{cases}

This gives the expected results of 4π steradians for the 3D sphere bounded by a surface of area 4πr² and 2π radians for the 2D circle bounded by a circumference of length 2πr. It also gives the slightly less obvious 2 for the 1D case, in which the origin-centered 1D "sphere" is the interval [ -r, r ] and this is bounded by 2 limiting points.

References

↑ "Archimedes on Spheres and Cylinders". Math Pages. 2015.
1 2 Mazonka, Oleg (2012). "Solid Angle of Conical Surfaces, Polyhedral Cones, and Intersecting Spherical Caps". arXiv:1205.1396.
↑ Van Oosterom, A; Strackee, J (1983). "The Solid Angle of a Plane Triangle". IEEE Trans. Biom. Eng. BME-30 (2): 125–126. doi:10.1109/TBME.1983.325207.
↑ "L'Huilier's Theorem – from Wolfram MathWorld". Mathworld.wolfram.com. 2015-10-19. Retrieved 2015-10-19.
↑ "Spherical Excess – from Wolfram MathWorld". Mathworld.wolfram.com. 2015-10-19. Retrieved 2015-10-19.
↑ "Area of a Latitude-Longitude Rectangle". The Math Forum @ Drexel. 2003.
↑ Jackson, FM (1993). "Polytopes in Euclidean n-space". Bulletin (Institute of Mathematics and its Applications) 29 (11/12): 172–174.

External links

Arthur P. Norton, A Star Atlas, Gall and Inglis, Edinburgh, 1969.
M. G. Kendall, A Course in the Geometry of N Dimensions, No. 8 of Griffin's Statistical Monographs & Courses, ed. M. G. Kendall, Charles Griffin & Co. Ltd, London, 1961
Weisstein, Eric W., "Solid Angle", MathWorld.

Classical mechanics SI units

Dimensions	—	L	L²	Dimensions	—	—	—
Linear/translational quantities				Angular/rotational quantities
T	time: $t$ s	absement: $A$ m s		T	time: $t$ s
—		distance: $d$ , position: $r$ , $s$ , $x$ , displacement m	area: $A$ m²	—		angle: $θ$ , angular displacement: $θ$ rad	solid angle: $Ω$ rad², sr
T⁻¹	frequency: $f$ s⁻¹, Hz	speed: $v$ , velocity: $v$ m s⁻¹	kinematic viscosity: $ν$ , specific angular momentum: $h$ m² s⁻¹	T⁻¹	frequency: $f$ s⁻¹, Hz	angular speed: $ω$ , angular velocity: $ω$ rad s⁻¹
T⁻²		acceleration: $a$ m s⁻²		T⁻²		angular acceleration: $α$ rad s⁻²
T⁻³		jerk: $j$ m s⁻³		T⁻³		angular jerk: $ζ$ rad s⁻³

M	mass: $m$ kg			ML²	moment of inertia: $I$ kg m²
MT⁻¹		momentum: $p$ , impulse: $J$ kg m s⁻¹, N s	action: $𝒮$ , actergy: $ℵ$ kg m² s⁻¹, J s	ML²T⁻¹		angular momentum: $L$ , angular impulse: $Δ L$ kg m² s⁻¹	action: $𝒮$ , actergy: $ℵ$ kg m² s⁻¹, J s
MT⁻²		force: $F$ , weight: $F g$ kg m s⁻², N	energy: $E$ , work: $W$ kg m² s⁻², J	ML²T⁻²		torque: $τ$ , moment: $M$ kg m² s⁻², N m	energy: $E$ , work: $W$ kg m² s⁻², J
MT⁻³		yank: $Y$ kg m s⁻³, N s⁻¹	power: $P$ kg m² s⁻³, W	ML²T⁻³		rotatum: $P$ kg m² s⁻³, N m s⁻¹	power: $P$ kg m²s⁻³, W

This article is issued from Wikipedia - version of the Sunday, February 14, 2016. The text is available under the Creative Commons Attribution/Share Alike but additional terms may apply for the media files.