Emden–Chandrasekhar equation explained

In astrophysics, the Emden–Chandrasekhar equation is a dimensionless form of the Poisson equation for the density distribution of a spherically symmetric isothermal gas sphere subjected to its own gravitational force, named after Robert Emden and Subrahmanyan Chandrasekhar.^{[1] [2]} The equation was first introduced by Robert Emden in 1907.^[3] The equation^[4] readswhere

\xi is the dimensionless radius and \psi is the related to the density of the gas sphere as \rho = \rho _{c} e ^{-\psi} , where \rho _{c} is the density of the gas at the centre. The equation has no known explicit solution. If a polytropic fluid is used instead of an isothermal fluid, one obtains the Lane–Emden equation. The isothermal assumption is usually modeled to describe the core of a star. The equation is solved with the initial conditions,\psi =0 , (d\psi)/(d\xi) =0 at \xi =0 .

The equation appears in other branches of physics as well, for example the same equation appears in the Frank-Kamenetskii explosion theory for a spherical vessel. The relativistic version of this spherically symmetric isothermal model was studied by Subrahmanyan Chandrasekhar in 1972.^[5]

Derivation

For an isothermal gaseous star, the pressure

p is due to the kinetic pressure and radiation pressurep = \rho (k _{B})/(W H) T + (4 \sigma)/(3 c) T ^{4}

where

\rho is the densityk _{B} is the Boltzmann constantW is the mean molecular weightH is the mass of the protonT is the temperature of the star\sigma is the Stefan–Boltzmann constantc is the speed of light

The equation for equilibrium of the star requires a balance between the pressure force and gravitational force

(1)/(r ^{2}) (d)/(\left( dr) (r ^{2})/(\rho) (dp)/(dr) \right) = - 4 \pi G \rho

where

r is the radius measured from the center and G is the gravitational constant. The equation is re-written as(k _{B} T)/(WH) (1)/(r ^{2}) (d)/(dr) \left(r ^{2 (d ln \rho)/(dr)} \right) = - 4 \pi G \rho Introducing the transformation\psi = ln (\rho _{c})/(\rho) , \xi = r \left( (4 \pi G \rho _{c} W H)/(k _{B} T) \right) ^{1/2}

where

\rho _{c} is the central density of the star, leads to(1)/(\xi ^{2}) (d)/(d\xi) \left(\xi ^{2} (d\psi)/(d\xi) \right) = e ^{-\psi}

The boundary conditions are

\psi =0 , (d\psi)/(d\xi) =0 at \xi =0

For

\xi\ll 1 , the solution goes like\psi = (\xi ^{2})/(6) - (\xi ^{4})/(120) + (\xi ^{6})/(1890) + …

Limitations of the model

Assuming isothermal sphere has some disadvantages. Though the density obtained as solution of this isothermal gas sphere decreases from the centre, it decreases too slowly to give a well-defined surface and finite mass for the sphere. It can be shown that, as

\xi \gg 1 ,(\rho)/(\rho _{c}) = e ^{-\psi} = (2)/(\left[ 1+ \xi ^{2}) (A)/(\cos\left( \xi ^{1/2}) (\sqrt 7)/(2) ln \xi + \delta\right) + O(\xi ^{-1} )\right]

where

A and \delta are constants which will be obtained with numerical solution. This behavior of density gives rise to increase in mass with increase in radius. Thus, the model is usually valid to describe the core of the star, where the temperature is approximately constant.^[6]

Singular solution

Introducing the transformation

x =1 /\xi transforms the equation tox ^{4} (d ^{2\psi})/(dx ^{2}) = e ^{-\psi}

The equation has a singular solution given by

e ^{- \psi _{s}} = 2 x ^{2,} or - \psi _{s} = 2 ln x + ln 2

Therefore, a new variable can be introduced as

- \psi = 2 ln x + z , where the equation for z can be derived,(d ^{2z})/(- dt ^{2}) (dz)/(dt) + e ^{z} -2 =0 , where t =ln x

This equation can be reduced to first order by introducing

y = (dz)/(dt) = \xi (d\psi)/(d\xi) - 2

then we have

y (dy)/(dz) - y + e ^{z-} 2 = 0

Reduction

There is another reduction due to Edward Arthur Milne. Let us define

u = (\xi e ^{-\psi})/(d\psi/d\xi) , v = \xi (d\psi)/(d\xi)

then

(u)/(v) (dv)/(du) =- (u -1)/(u + v -3)

Properties

• If

\psi(\xi) is a solution to Emden–Chandrasekhar equation, then \psi(A\xi) -2ln A is also a solution of the equation, where A is an arbitrary constant.

• The solutions of the Emden–Chandrasekhar equation which are finite at the origin have necessarily

d\psi/d\xi =0 at \xi =0

See also

• Lane–Emden equation
• Frank-Kamenetskii theory
• Chandrasekhar's white dwarf equation

Notes and References

1. Chandrasekhar, Subrahmanyan, and Subrahmanyan Chandrasekhar. An introduction to the study of stellar structure. Vol. 2. Courier Corporation, 1958.
2. Chandrasekhar, S., and Gordon W. Wares. "The Isothermal Function." The Astrophysical Journal 109 (1949): 551-554.http://articles.adsabs.harvard.edu/cgi-bin/nph-iarticle_query?1949ApJ...109..551C&defaultprint=YES&filetype=.pdf
3. Emden, R. (1907). Gaskugeln: Anwendungen der mechanischen Wärmetheorie auf kosmologische und meteorologische Probleme. B. Teubner.
4. Kippenhahn, Rudolf, Alfred Weigert, and Achim Weiss. Stellar structure and evolution. Vol. 282. Berlin: Springer-Verlag, 1990.
5. Chandrasekhar, S. (1972). A limiting case of relativistic equilibrium. In General Relativity (in honor of J. L. Synge), ed. L. O'Raifeartaigh. Oxford. Clarendon Press (pp. 185-199).
6. Henrich, L. R., & Chandrasekhar, S. (1941). Stellar Models with Isothermal Cores. The Astrophysical Journal, 94, 525.