# Solving the Tolman–Oppenheimer–Volkoff (TOV) equation

The Pressure of a static spherical object (say star), which has the Schwarzchild metric outside it, satisfies the following differential equation called the TOV equation.

$$\frac{\mbox{d}P}{\mbox{d}r}=-\left(P+\rho\right)\frac{m(r)+4\pi r^3P }{r\left(r-2m(r)\right)}$$

Here $\rho$ the mass density is a function of $r$, and $m(r)=4 \pi\int_0^r\rho(r^{\prime}){r^{\prime}}^2 d{r^{\prime}}$.

How do I solve this equation? Wald states the solution as follows; but Please could someone suggest how I solve this analytically?

($\rho_0$ is the density at the radius of the star $R$, and $M$ is the total mass $M=4 \pi\int_0^R\rho(r^{\prime}){r^{\prime}}^2 d{r^{\prime}}$)

$$P(r)=\rho_0\left(\frac{\left(1-\frac {2M}{R}\right)^\frac12 - \left(1-\frac {2Mr^2 }{R^3 }\right)^\frac12 }{\left(1-\frac {2Mr^2 }{R^3 }\right)^\frac12 -3\left(1-\frac{2M}{R}\right)^\frac12 }\right)$$

• Note that there is an additional important assumption between your first equation and your second. In order to solve TOV you need an equation of state $\rho = \rho(P)$. In the case at hand you assume an incompressible equation of state $\rho = \rho_0 = \text{constant}$, which is a highly idealised case. With $\rho$ a constant you can explicitly find $m(r)$ and simplifications follow... Jul 4, 2013 at 1:41
• @ramanujan_dirac : You may note that $\frac{dm}{dr} = 4 \pi r^2\rho$, so you can replace $\rho$ in the first equation, so getting a equation with $P, \frac{dP}{dr},m, \frac{dm}{dr}$. Unfortunately, this is a complex equation.. Jul 4, 2013 at 6:31
• @Trimok It's pretty simple really. Put in $m = 4\pi r^3 \rho_0/ 3$ (Yes, the density is constant here. There is no analytic solution for the general equation of state.) Cancel out some factors of $r$ and notice that you can seperate variables. You get $\mathrm{d}P/(\text{quadratic in }P) = f(r) \mathrm{d}r$ and both integrals can be done by elementary techniques. Jul 4, 2013 at 7:04
• @MichaelBrown Yes, but I was not considering the constant case... Jul 4, 2013 at 7:08

Summarizing and expanding on the comments...

First, the given solution only holds in the case $\rho \equiv \rho_0$ is constant (though Wald's wording could stand to be clearer on this point). You won't find a nice general expression for an arbitrary equation of state, which is why numerical results are key for modeling real compact objects like neutron stars.

Proceeding with our assumption, the mass inside a sphere is given by $m = (4\pi/3) \rho_0 r^3$, with the total mass related to the total radius in the same way: $M = (4\pi/3) \rho_0 R^3$. Then the differential equation can be rewritten $$\frac{\mathrm{d}P}{\mathrm{d}r} = -(P + \rho_0) \frac{(4\pi/3)\rho_0r^3 + 4\pi r^3P}{r(r - 2(4\pi/3)\rho_0r^3)} = -\frac{4\pi}{3} (P + \rho_0) \frac{(\rho_0 + 3P)r}{1-8\pi\rho_0r^2/3}.$$

We now proceed with separation of variables. The only physics left is in the boundary condition: $P = 0$ at $r = R$.1 With this in mind, we can write $$\int_P^0 \frac{\mathrm{d}P}{(P+\rho_0)(3P+\rho_0)} = \int_r^R \frac{r\ \mathrm{d}r}{2\rho_0r^2-3/4\pi}.$$ Both of these integrals are quite solvable using cute tricks from high school calculus, nothing more advanced required.2

The result of integrating is $$\frac{1}{2\rho_0} \log\left(\frac{P+\rho_0}{3P+\rho_0}\right) = \frac{1}{4\rho_0} \log\left(\frac{8\pi\rho_0R^2/3-1}{8\pi\rho_0r^2/3-1}\right).$$ Using the earlier definition of $M$ this is clearly equivalent to $$\frac{P+\rho_0}{3P+\rho_0} = \left(\frac{1-2M/R}{1-2Mr^2/R^3}\right)^{1/2},$$ which can be solved for $P$ to yield the desired result.

Finally, it should be noted that ODE solving is not the most interesting thing going on here. What is far more intriguing is how one gets the differential equation to begin with - it is not trivial. Remarkably, Karl Schwarzschild did much of this computation in his 1916 paper "Über das Gravitationsfeld einer Kugel aus inkompressibler Flüssigkeit nach der Einsteinschen Theorie" ("On the Gravitational Field of a Sphere of Incompressible Fluid according to Einstein’s Theory") just a few months after Einstein published his GR field equations.

1 Alternatively one could say $P = P_\mathrm{c}$ at $r = 0$ and integrate the other way, finding a formula for pressure at any radius in terms of the central pressure and mass. The solution can then be inverted to find the radius in terms of the pressure, and the value of this at $P = 0$ would give the size of the object. This procedure make sense if one has a better handle on the microphysics responsible for pressure than on the actual sizes of objects, as is often the case.

2 For the pressure integral, you can complete the square to show that the integrand looks like the derivative of the inverse hyperbolic tangent (or equivalently the inverse tangent with $\sqrt{-1}$ in appropriate places). Then it's a simple matter to invert the exponential formula for $\tanh$ to find the logarithmic version of $\tanh^{-1}$.

• +1 Not only did Schwarzschild do it so soon after Einstein's publication, he did it while suffering from a horribly painful autoimmune disease he developed while serving in the trenches of the first world war. They don't make em like they used to! Jul 5, 2013 at 7:27
• @Chris: Thanks. I was under the impression $\rho$ is not constant, and hence didn't attempt with $\rho \equiv \rho_0$, the integral in this case as you show is quite trivial.
– user7757
Jul 10, 2013 at 4:20
• @Elnur What you suggested isn't wrong, but it's actually equivalent to what is already there ;)
– user10851
May 2, 2016 at 19:40