Problem 10 - Classical Mechanics - An Undergraduate Text

Question

A particle of mass \(m\) moves under the central field $$ \boldsymbol{F}=-m\left(\frac{\gamma e^{-\epsilon r / a}}{r^{2}}\right) \widehat{\boldsymbol{r}}, $$ where \(\gamma, a\) and \(\epsilon\) are positive constants. Find the apsidal angle for a nearly circular orbit of radius \(a\). When \(\epsilon\) is small, show that the perihelion of the orbit advances by approximately \(\pi \epsilon\) on each revolution.

Step-by-Step Solution

Answer

The apsidal angle of the orbit for this central field is \(\pi(1 + \epsilon)\), and under small \(\epsilon\), the perihelion of the orbit advances by approximately \(\pi \epsilon\) on each revolution.

Step 1: Derive equation of motion

In central force problems, the equation of motion can be derived from the equation for radial force by making use of the conservation of angular momentum \(L\), such that \(L^2 = m^2r^4\dot{\theta}^2\). So the equation of motion is given by \(-\frac{d}{dr}(U_{eff}) = -\frac{\gamma e^{-\epsilon r / a}}{r^{2}} + \frac{L^2}{m^2r^3}\), where the radial potential \(U_{eff}\) is associated with the radial force.

Step 2: Solve for radial coordinate

To find the small changes in the radial coordinate \(r\), differentiate \(U_{eff}\) with respect to \(r\) and solve for \(\dot{r}\). The equation becomes \(\dot{r} = \sqrt{ \frac{2}{m} [E - U_{eff}(r)] }\).

Step 3: Determine small changes in angular coordinate

Next, determine the small changes in the angular coordinate \(\theta\) using the conservation of angular momentum again, given by the equation \(\dot{\theta} = \frac{L}{m r^2}\). Divide \(\frac{dr}{d\theta}\) by \(\frac{d^2r}{d\theta^2}\) to obtain the equation of trajectory in polar coordinates, obtaining as result Binet's equation.

Step 4: Solve Binet's equation

At this point, obtain the solution of Binet's equation using small perturbations. By setting \(u = \frac{1}{r}\) and \(\frac{du}{d \theta} = - \epsilon u\), we can integrate with respect to \(\theta\) to find \(u = \frac{1}{a}(1 + \epsilon cos(\theta))\). Solve for \(r\) by taking the reciprocal of \(u\) and use it to find the apsidal angle.

Step 5: Calculate the change in perihelion

Finally, calculate the change in the perihelion, which can be given by \(\Delta \phi = 2 \pi - 2 \arccos(- \frac{1}{\epsilon}) \approx 2 \pi (1 - \frac{1}{\epsilon}) = \pi \epsilon\), for small \(\epsilon\).

Key Concepts

Apsidal Angle Calculation

The apsidal angle, sometimes called the angle of precession, is a crucial concept in central field mechanics. It describes the angular separation between the successive perihelion and aphelion points of an orbital path. For a central force field given by \( \boldsymbol{F}=-m\left(\frac{\gamma e^{-\epsilon r / a}}{r^{2}}\right)\widehat{\boldsymbol{r}}, \), where \( m, \gamma, a, \) and \( \epsilon \) are constants with \( \epsilon \) being small, we are particularly interested in the apsidal angle for a nearly circular orbit of radius \( a \).

To calculate this angle accurately, we analyze the radial motion through the effective potential energy and use the conservation of angular momentum. The trajectory equation is simplified by assuming a nearly circular orbit, which allows using small angular approximations. The calculated apsidal angle comes to light as an outcome of the perturbations caused by the exponential term in the force field. This approach not only simplifies the complex dynamics of the orbit but also illuminates the subtle effects of the force field parameters on orbital paths. In essence, the apsidal angle serves as a measure of the orbit's rotation after each complete revolution.

Binet's Equation

Understanding Binet's equation is integral to solving orbital dynamics problems in a central force field. This elegant mathematical tool encapsulates the motion of a particle under central field forces in polar coordinates. Binet's equation for classical orbits is derived by manipulating the differential equations governing the orbital motion.

For the given central force field, Binet's equation translates the complex motion into a straightforward relation involving the radius \( r \) and the angular coordinate \( \theta \). By introducing the variable \( u = \frac{1}{r} \), we simplify the problem and can identify the radial changes over the angular coordinate. This approach breaks down the problem into minuscule effects which, once accumulated, reveal the overall precession of the orbit. Applying Binet’s equation yields an interpretable pattern of how the particle trails along its path and offers a practical method for analyzing orbits that are affected by slight perturbations.

Perihelion Advance

Perihelion advance, or the precession of the perihelion, is a phenomenon observed in orbital mechanics where the point of closest approach of an orbiting body to the central mass shifts over time. For orbits influenced by non-inverse-square forces, like the one introduced in the exercise, the perihelion advances with each orbit.

Small perturbations in the force can lead to a gradual but consistent precession of the orbit. In our context, the exponential decay term characterizes this perturbation which causes the orbit to rotate progressively. When \( \epsilon \) is slight, the perihelion advances approximately by \( \pi \epsilon \) after each complete orbit. This analytical approximation brings out the influence of the exponential term on the orbital precession and reinforces the understanding of celestial mechanics, further echoing the nuanced interplay between mathematical abstractions and physical reality.

Conservation of Angular Momentum

The principle of conservation of angular momentum is fundamental in the study of orbital motion and central field forces. It states that if no external torque acts on a system, the total angular momentum of that system remains constant. In the scenario provided, the central force field ensures a torque-free environment for the orbiting particle, thus guaranteeing the conservation of angular momentum.

Employing this principle simplifies the complex equations of motion, allowing us to equate the moment of inertia and the angular velocity at different points along the orbital path. This conservation law ties together the radial and angular coordinates of the particle's trajectory, providing a backbone to derive Binet's equation and analyze the system's stability. Essentially, the conservation of angular momentum is not merely a convenient simplification but an intrinsic property that shapes the very fabric of orbital mechanics.