Academy
PhysicsElectromagnetismElectrostaticsAdvanced Integration: Variable Charge Densities

Advanced Integration: Variable Charge Densities

Mastering Coulomb's Law when charge distribution is non-uniform across 1D structures.

Advanced Integration: Variable Charge Densities

In our previous analysis of continuous bodies, we assumed the charge was spread uniformly, allowing us to treat the linear charge density λ=Q/L\lambda = Q/L as a constant. JEE Advanced routinely shatters this assumption.

When charge density is a function of position—such as λ=f(x)\lambda = f(x) or λ=f(θ)\lambda = f(\theta)—the setup of the differential element dqdq must rigorously incorporate this function before any integration can occur.

Case 1: The Non-Uniform Rod

The Problem: A rod of length LL lies on the x-axis from x=ax = a to x=a+Lx = a + L. Its linear charge density varies proportionally with the distance from the origin, given by λ=λ0x\lambda = \lambda_0 x, where λ0\lambda_0 is a positive constant. Find the exact electrostatic force exerted on a positive point charge qq placed at the origin.

Step 1: Define the Differential Element

Consider an infinitesimal element of the rod of length dxdx at a position xx. The charge dqdq on this element is no longer a simple fraction of the total charge. It is dictated by the local density:

dq=λdx=(λ0x)dxdq = \lambda \cdot dx = (\lambda_0 x) dx

Step 2: Apply Coulomb's Law

Treating dqdq as a point charge, the infinitesimal repulsive force dFdF on the test charge qq (directed along the negative x-axis) is:

dF=14πϵ0qdqx2dF = \frac{1}{4\pi\epsilon_0} \frac{q \cdot dq}{x^2}

Substitute our expression for dqdq:

dF=14πϵ0q(λ0xdx)x2=qλ04πϵ0dxxdF = \frac{1}{4\pi\epsilon_0} \frac{q (\lambda_0 x \cdot dx)}{x^2} = \frac{q \lambda_0}{4\pi\epsilon_0} \frac{dx}{x}

Notice the Mathematical Shift: Because the charge density grows linearly with xx, it cancels one power of xx in the 1/x21/x^2 denominator of Coulomb's Law. This fundamentally changes the nature of the integral.

Step 3: Integration

We integrate dFdF over the physical boundaries of the rod, from x=ax = a to x=a+Lx = a + L:

F=aa+Lqλ04πϵ0dxxF = \int_{a}^{a+L} \frac{q \lambda_0}{4\pi\epsilon_0} \frac{dx}{x}

Pulling the constants out of the integral:

F=qλ04πϵ0aa+L1xdxF = \frac{q \lambda_0}{4\pi\epsilon_0} \int_{a}^{a+L} \frac{1}{x} dx

The integral of 1/x1/x is the natural logarithm, lnx\ln|x|:

F=qλ04πϵ0[ln(x)]aa+LF = \frac{q \lambda_0}{4\pi\epsilon_0} \left[ \ln(x) \right]_{a}^{a+L}

F=qλ04πϵ0(ln(a+L)ln(a))F = \frac{q \lambda_0}{4\pi\epsilon_0} \left( \ln(a+L) - \ln(a) \right)

Applying logarithm properties (ln(A)ln(B)=ln(A/B)\ln(A) - \ln(B) = \ln(A/B)), we arrive at the final force:

F=qλ04πϵ0ln(1+La)F = \frac{q \lambda_0}{4\pi\epsilon_0} \ln\left(1 + \frac{L}{a}\right)


Case 2: The Semi-Circular Ring with Angular Variation

The Problem: A thin semi-circular ring of radius RR lies in the upper half of the xy-plane, centered at the origin. Its linear charge density varies with the angle θ\theta (measured from the positive x-axis) according to the function λ=λ0cosθ\lambda = \lambda_0 \cos\theta. Find the force on a positive point charge qq located at the origin.

Step 1: The Geometry of the Differential Element

Let's take a small angular slice dθd\theta at an angle θ\theta. The arc length of this slice is ds=Rdθds = R d\theta. The charge on this slice is:

dq=λds=(λ0cosθ)(Rdθ)dq = \lambda \cdot ds = (\lambda_0 \cos\theta) (R d\theta)

Step 2: The Force Vector and Symmetry

The force dFdF from this element acts radially inward (since λ\lambda is positive for 0<θ<π/20 < \theta < \pi/2 and qq is positive). The magnitude is:

dF=14πϵ0qdqR2=qλ0cosθRdθ4πϵ0R2=qλ04πϵ0RcosθdθdF = \frac{1}{4\pi\epsilon_0} \frac{q \cdot dq}{R^2} = \frac{q \lambda_0 \cos\theta R d\theta}{4\pi\epsilon_0 R^2} = \frac{q \lambda_0}{4\pi\epsilon_0 R} \cos\theta d\theta

We must resolve this into components:

  • dFx=dFcosθdF_x = -dF \cos\theta
  • dFy=dFsinθdF_y = -dF \sin\theta

The Symmetry Argument: The charge density is λ0cosθ\lambda_0 \cos\theta. For θ\theta between 00 and π/2\pi/2 (first quadrant), λ\lambda is positive. For θ\theta between π/2\pi/2 and π\pi (second quadrant), cosθ\cos\theta is negative, so λ\lambda is negative. Because the first quadrant pushes the charge qq away, and the second quadrant pulls it in, the y-components will completely cancel out. We only need to integrate the x-component.

Step 3: Integrating the X-Component

Fx=0π(dFcosθ)F_x = \int_{0}^{\pi} (-dF \cos\theta)

Substitute our expression for dFdF:

Fx=0π(qλ04πϵ0Rcosθ)cosθdθF_x = -\int_{0}^{\pi} \left( \frac{q \lambda_0}{4\pi\epsilon_0 R} \cos\theta \right) \cos\theta d\theta

Fx=qλ04πϵ0R0πcos2θdθF_x = -\frac{q \lambda_0}{4\pi\epsilon_0 R} \int_{0}^{\pi} \cos^2\theta d\theta

Using the half-angle identity cos2θ=1+cos(2θ)2\cos^2\theta = \frac{1 + \cos(2\theta)}{2}:

Fx=qλ04πϵ0R0π(12+cos(2θ)2)dθF_x = -\frac{q \lambda_0}{4\pi\epsilon_0 R} \int_{0}^{\pi} \left( \frac{1}{2} + \frac{\cos(2\theta)}{2} \right) d\theta

Fx=qλ08πϵ0R[θ+sin(2θ)2]0πF_x = -\frac{q \lambda_0}{8\pi\epsilon_0 R} \left[ \theta + \frac{\sin(2\theta)}{2} \right]_{0}^{\pi}

Evaluating from 00 to π\pi, the sine terms vanish, leaving just π\pi:

Fx=qλ0π8πϵ0RF_x = -\frac{q \lambda_0 \pi}{8\pi\epsilon_0 R}

The force is entirely along the negative x-axis. This problem demonstrates that understanding the physical symmetry dictated by the density function is just as critical as executing the integral.