The problem of finding the total charge from a variable volume charge density is a classic application of Gauss's Law and integration in electrostatics. It tests your ability to set up a volume integral in spherical coordinates and execute it flawlessly. Let's dive into the elegant mathematics behind this charge distribution.
The Anatomy of the Charge Distribution
Imagine a solid sphere of radius R. Unlike a simple conductor where charge resides only on the surface, here the charge is distributed throughout the entire volume. However, it's not spread evenly. The volume charge density is given by:
This function tells us two things. First, the 1/r2 term means the charge is highly concentrated near the center and thins out as you move away. Second, the exponential term e−2r/a causes a rapid, exponential decay of the charge density.
Slicing the Sphere
The Volume Element
To find the total charge Q, we cannot simply multiply density by volume. We must integrate. We slice the sphere into infinitely many thin, concentric spherical shells.
Consider one such shell at a distance r from the center, with an infinitesimally small thickness dr. If we were to unroll this shell, it would resemble a flat sheet with a surface area of 4πr2 and a thickness of dr. Therefore, the volume of this tiny shell is:
The Elegance of the Integral
The total charge Q is the sum of the charges in all these tiny shells, from the center (r=0) to the surface (r=R). We set up the integral:
Substituting our expressions for ρ(r) and dV:
Q=∫0R(r2Ae−2r/a)(4πr2dr)
Here is where the mathematical elegance of the problem reveals itself. The r2 in the denominator of the density perfectly cancels the r2 in the volume element!
We are left with a straightforward exponential integral.
Executing the Integration
Integrating e−2r/a with respect to r yields:
Now, we apply the limits from 0 to R:
Since e0=1, we can simplify the expression:
To make the equation neater, we absorb the negative sign into the bracket:
Rescuing the Radius
Our ultimate goal is to find the radius R. We need to isolate the exponential term. First, divide both sides by 2πaA:
Rearranging to isolate the exponential:
To bring R down from the exponent, we take the natural logarithm (ln or loge) of both sides:
Finally, multiply by −2a:
Using the logarithmic property −log(x)=log(x−1)=log(1/x), we can rewrite the final answer to match the given options:
This beautiful result shows exactly how the radius of the sphere depends on the total charge it contains and the parameters of its internal charge distribution.