Calculus III: Spherical Coordinates

When integrating over three-dimensional regions, Cartesian coordinates often turn simple volumes into computational nightmares. For spheres, spherical shells, and regions with central symmetry, spherical coordinates provide a far more elegant and powerful framework. Mastering this system is not just a mathematical exercise; it is essential for solving real-world problems in engineering and physics, from calculating the gravitational field of a planet to determining the electric potential around a charged sphere.

Defining the Spherical Coordinate System

In the spherical coordinate system, a point $P$ in space is described by three values: its radial distance from the origin, and two angles. The convention we will use is:

• $\rho$ (rho): The radial distance from the origin to the point ($\rho \ge 0$).
• $\theta$ (theta): The polar angle, measured from the positive z-axis downward ($0\le \theta \le \pi$). Think of it as the "colatitude."
• $\varphi$ (phi): The azimuthal angle, measured from the positive x-axis in the xy-plane, just like the angle in 2D polar or cylindrical coordinates ($0\le \varphi <2\pi$).

Understanding this geometry is the first step. The conversion formulas between spherical and Cartesian coordinates are critical: $$x=\rho \sin \theta \cos \varphi$$ $$y=\rho \sin \theta \sin \varphi$$ $$z=\rho \cos \theta$$ Conversely, to go from Cartesian to spherical coordinates, you use: $$\rho =\sqrt{x^2+y^2+z^2}$$ $$\theta =\arccos \left(\frac{z}{\rho }\right)$$ $$\varphi =\arctan \left(\frac{y}{x}\right)\text{ (with quadrant awareness)}$$

The coordinate surfaces are key to visualization: constant $\rho$ defines a sphere, constant $\theta$ defines a cone, and constant $\varphi$ defines a half-plane. This geometric insight directly informs how you set up integrals.

The Volume Element and Jacobian Derivation

You cannot simply replace $dV$ with $d\rho d\theta d\varphi$. The transformation from Cartesian to spherical coordinates distorts space, stretching it more near the origin than far away. To account for this distortion, we must multiply by the Jacobian determinant of the transformation, $|J|$.

The Jacobian is the determinant of the matrix of all first-order partial derivatives $(x,y,z)$ with respect to $(\rho ,\theta ,\varphi )$. Starting from our conversion formulas: $$J=\frac{\partial (x,y,z)}{\partial (\rho ,\theta ,\varphi )}=\left|\begin{array}{ccc}\sin \theta \cos \varphi & \rho \cos \theta \cos \varphi & -\rho \sin \theta \sin \varphi \\ \sin \theta \sin \varphi & \rho \cos \theta \sin \varphi & \rho \sin \theta \cos \varphi \\ \cos \theta & -\rho \sin \theta & 0\end{array}\right|$$

Computing this determinant (a standard exercise) yields: $$|J|={\rho }^2\sin \theta$$

Therefore, the volume element in spherical coordinates is: $$dV={\rho }^2\sin \theta d\rho d\theta d\varphi$$ This ${\rho }^2\sin \theta$ factor is non-negotiable; forgetting it is the single most common error. Geometrically, ${\rho }^{2}d\rho$ accounts for the volume of a thin spherical shell, while $\sin \theta d\theta d\varphi$ accounts for the area element on that shell's surface.

Setting Up and Evaluating Triple Integrals

To evaluate a triple integral ${\iiint }_{E}f(x,y,z)dV$ in spherical coordinates, you follow a systematic process:

1. Describe the region $E$ in spherical coordinates. Substitute $x,y,z$ using the conversion formulas. The bounds will typically be constants.
2. Replace $dV$ with ${\rho }^2\sin \theta d\rho d\theta d\varphi$.
3. Determine the order and limits of integration. The most common order is $d\rho$ first, then $d\theta$, then $d\varphi$.

• $\rho$ limits: From the origin ($\rho =0$) to the radial boundary of the region (e.g., $\rho =R$ for a solid sphere, or $\rho =a\sec \theta$ for a region bounded by a cone and a plane).
• $\theta$ limits: Sweep from the positive z-axis ($\theta =0$) down to the boundary, often a cone ($\theta =\alpha$) or the xy-plane ($\theta =\pi /2$).
• $\varphi$ limits: Rotate fully around the z-axis from $0$ to $2\pi$ for full symmetry, or a subset for wedge-shaped regions.

Example: Find the volume of a sphere of radius $R$. Here, $f(x,y,z)=1$. $$V={\int }_{\varphi =0}^{2\pi }{\int }_{\theta =0}^{\pi }{\int }_{\rho =0}^R{\rho }^2\sin \theta d\rho d\theta d\varphi$$ Evaluating from the inside out: $$={\int }_0^{2\pi }{\int }_0^{\pi }{\left[\frac{{\rho }^3}{3}\right]}_0^R\sin \theta d\theta d\varphi ={\int }_0^{2\pi }{\int }_0^{\pi }\frac{R^3}{3}\sin \theta d\theta d\varphi$$ $$={\int }_0^{2\pi }{\left[-\frac{R^3}{3}\cos \theta \right]}_0^{\pi }d\varphi ={\int }_0^{2\pi }\left(\frac{2R^3}{3}\right)d\varphi =\frac{2R^3}{3}\cdot 2\pi =\frac{4}{3}\pi R^3$$

Applications in Physics and Engineering

The true power of spherical coordinates is revealed in physical applications. Consider calculating the gravitational field inside a planet of uniform density $\delta$. By symmetry, you only need the field at a distance $r$ from the center. You can set up an integral over the spherical mass interior to $r$, using the inverse-square law. The spherical symmetry makes the angular integrals trivial, leaving a simple radial integral that yields the classic result: the gravitational force increases linearly with $r$ inside the uniform sphere.

Similarly, in electrostatics, finding the electric potential due to a spherical shell of charge is vastly simpler in spherical coordinates. The integral form of Coulomb's law involves integrating $1/|\mathbf{r}-{\mathbf{r}}^{\prime }|$ over the source charge distribution. For a spherical shell, the symmetry allows you to align your coordinate system so the field point lies on the z-axis, dramatically simplifying the distance calculation and the ensuing integrals. This is how one derives the crucial result that the potential outside a spherical shell is identical to that of a point charge, while inside it is constant.

Choosing the Optimal Coordinate System

You now have three primary tools: Cartesian, cylindrical, and spherical coordinates. The choice is a strategic decision that can reduce a triple integral from a page of algebra to three lines. Follow this framework:

• Use Cartesian Coordinates: When the region $E$ is bounded primarily by planes (e.g., a rectangular box, a tetrahedron). The integrand $f(x,y,z)$ should also be simple in these variables.
• Use Cylindrical Coordinates: When there is symmetry or structure around an axis (usually the z-axis). Look for regions involving cylinders, circular cross-sections, or cones. The transformation is $x=r\cos \varphi ,y=r\sin \varphi ,z=z$, with $dV=rdrd\varphi dz$.
• Use Spherical Coordinates: When there is symmetry about a point (the origin). Spheres, spherical shells, and regions between cones are prime candidates. Also use it when the integrand contains terms like $x^2+y^2+z^2$, which simplifies to ${\rho }^2$.

For a solid hemisphere, for instance, spherical coordinates are ideal. For a cylinder capped by a paraboloid, cylindrical coordinates are best. For a wedge cut from a rectangular block, stick with Cartesian.

Common Pitfalls

1. Forgetting the Jacobian (${\rho }^2\sin \theta$): This is the cardinal sin. Without it, your volume and mass calculations will be drastically wrong. Always write $dV={\rho }^2\sin \theta d\rho d\theta d\varphi$ as your first step.
2. Mixing Up $\theta$ and $\varphi$: The convention used here (and in most calculus textbooks and engineering fields) is that $\theta$ is the polar angle from the z-axis. In some physics contexts, these are swapped. Always check the definition: $\theta$ ranges from $0$ to $\pi$, which is the domain of $\sin \theta$ needed for the Jacobian.
3. Incorrect Limits of Integration: Assuming the $\rho$ limit is always a constant is a frequent oversight. If your region is bounded by a cone ($\theta =\alpha$) and a sphere ($\rho =R$), the $\rho$ limit is constant. But if it's bounded by a cone and a plane $z=a$, the upper $\rho$ limit becomes $a\sec \theta$, a function of $\theta$. Always sketch the $\rho$-line from the origin outward through the region.
4. Applying Spherical Coordinates Inappropriately: Don't force spherical coordinates onto a cylindrical region. If the problem has clear axial symmetry (like a long wire or a pipe), cylindrical coordinates are almost always simpler.

Summary

• Spherical coordinates $(\rho ,\theta ,\varphi )$ describe points in 3D space using radial distance and two angles, ideal for regions with spherical symmetry.
• The transformation requires the use of a Jacobian determinant, yielding the critical volume element $dV={\rho }^2\sin \theta d\rho d\theta d\varphi$.
• Setting up triple integrals involves describing the region in spherical coordinates and establishing bounds, typically integrating in the order $d\rho$, $d\theta$, then $d\varphi$.
• This coordinate system is indispensable for solving foundational physics and engineering problems involving gravitational fields and electric potentials of spherical masses and charge distributions.
• Choosing between Cartesian, cylindrical, and spherical coordinates is a strategic skill; select spherical coordinates for point-symmetric regions bounded by spheres and cones.