# 3.3: Effective Potential

- Page ID
- 43226

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\id}{\mathrm{id}}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\kernel}{\mathrm{null}\,}\)

\( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\)

\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\)

\( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)

\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)

\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vectorC}[1]{\textbf{#1}} \)

\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

In the previous section, we applied the geodesic equation to three of the four Schwarzschild coordinates. Applying it to the r-coordinate yields a result that is complicated enough to warrant its own section. The result is

\[\left(\frac{dr}{d\tau}\right)^2=\left(\frac{E}{m}\right)^2-\left (1-\frac{2M}{r}\right)\left (1+\frac{L^2}{m^2r^2}\right).\label{eq:SchwarzschildEffective}\]

Now let's look at this term by term. On the left, there is something that looks like a velocity squared. Next is a term that is constant and represents the total energy. Finally, there is a term that depends only on position and other constants. This is very similar to a familiar equation from Newtonian Mechanics:

\[KE=\frac{1}{2}mv^2=E_\text{total}-PE\]

If you make a graph of potential energy, then you can use the graph to come to conclusions about the motion. For example, let's graph \(PE=\frac{1}{2}kx^2\). The total energy is a constant, which is a horizontal line on the graph. The kinetic energy is then given by the difference of the two curves, as shown in Fig. 3.3.1.

Since kinetic energy must always be greater than or equal to zero, we see that the object is forced to oscillate between the points of intersection (called **turning points**) of the two curves. Furthermore, the smaller the total energy the smaller the range of oscillation. The point, though, is that we can obtain qualitative information about the velocity by plotting the other two parts. We can repeat the process with Equation \ref{eq:SchwarzschildEffective} by plotting \(\frac{E}{m}\) (as a horizontal line) and

\[\frac{PE_\text{eff}(r)}{m}=\sqrt{\left(1-\frac{2M}{r}\right)\left(1+\frac{L^2}{m^2r^2}\right)},\]

where \(PE_\text{eff}(r)\) is called the **effective potential energy**.

An effective potential energy is not technically a potential energy but can still be used to obtain qualitative information about motion.

Fig. 3.3.2 plots this effective potential energy for \(\frac{L}{m}=3.75M\) with dimensionless axes.

Describe the motion depicted in Fig. 3.3.2. How would the motion be different if \(\frac{E}{m}=0.98\)?

**Answer**-
There are three turning points (i.e. points of intersection) on the graph, so the motion depends on which region the particle is in. If \(6.5\le \frac{r}{M}\le 15.5\), then the particle is in a trapped orbit. We know this because the particle has angular momentum, which mean it is orbiting around, but there is a minimum and a maximum r-value to that orbit (given by the turning points). This could, for example, be an elliptical orbit, though we can't say for sure that it is

*exactly*an ellipse. The other possibility is if \(\frac{r}{M}\le 3.5\), in which case there is a maximum orbital r-value but no minimum. That particle may swing around a bit, but it is doomed to fall to the event horizon.If the total energy were \(\frac{E}{m}=0.98\), then there would only be one turning point instead of three (it would intersect way off to the right). In that case, the orbit would have maximim r-value but no minimum, which again means that it is doomed to fall to the event horizon.

Fig. 3.3.2 shows the effective potential energy for \(\frac{L}{m}=0\). Describe the resulting motion. Is that what you would expect?

**Answer**-
Zero angular momentum corresponds to a particle moving radially. The fact that there is only one turning point makes sense, since it shows that a particle initially moving radially away will reach a maximum r-value then fall back toward the event horizon. Note that the particle never passes the event horizon, which is consistent with what we learned from looking at light cones.

Circular orbits are possible wherever the effective potential energy has zero slope. Take the derivative of the effective potential energy with respect to *r*, then set it equal to zero. Use it to show that

\[r^2-\frac{L^2}{Mm^2}r+3\frac{L^2}{m^2}=0.\nonumber\]

Solve the quadratic equation from Box 3.3.1 to show that

\[r=\frac{L^2}{2m^2M}\left [1\pm \sqrt{1-\frac{12M^2m^2}{L^2}}\right].\nonumber\]

Using the result from the previous box, what is the minimum value of \(\frac{L}{m}\) for which you get a real number for \(r\)? Express your answer in terms of \(M\).

Use the minimum value of \(\frac{L}{m}\) from the previous box to determine the r-value of the orbit (in terms of \(M\)). This is called the **Innermost Stable Circular Orbit** (ISCO).