Hayashi limit

A Hertzsprung–Russell diagram displays a plot of a star's surface temperature against the luminosity. On this diagram, the Hayashi limit forms a nearly vertical line at about 2,500 K. The outer layers of low temperature stars are always convective, and models of stellar structure for fully convective stars do not provide a solution to the right of this line. Thus in theory, stars are constrained to remain to the left of this limit during all periods when they are in hydrostatic equilibrium, and the region to the right of the line forms a type of "forbidden zone". Note, however, that there are exceptions to the Hayashi limit. These include collapsing protostars, as well as stars with magnetic fields that interfere with the internal transport of energy through convection.

Red giants are stars that have expanded their outer envelope in order to support the nuclear fusion of helium. This moves them up and to the right on the H-R diagram. However, they are constrained by the Hayashi limit not to expand beyond a certain radius. Stars that find themselves across the Hayashi limit have large convection currents in their interior driven by massive temperature gradients. Additionally, those stars states are unstable so the stars rapidly adjust their states, moving in the Hertzprung-Russel diagram until they reach the Hayashi limit.

When lower mass stars in the main sequence start expanding and becoming a red giant the stars revisit the Hayashi track. The Hayashi limit constrains the asymptotic giant branch evolution of stars which is important in the late evolution of stars and can be observed, for example, in the ascending branches of the Hertzsprung–Russell diagrams of globular clusters, which have stars of approximately the same age and composition.

The Hayashi limit is named after Chūshirō Hayashi, a Japanese astrophysicist.

Despite its importance to protostars and late stage main sequence stars, the Hayashi limit was only recognized in Hayashi's paper in 1961. This late recognition may be because the properties of the Hayashi track required numerical calculations that were not fully developed before.

Derivation of the limit

We can derive the relation between the luminosity, temperature and pressure for a simple model for a fully convective star and from the form of this relation we can infer the Hayashi limit. This is an extremely crude model of what occurs in convective stars, but it has good qualitative agreement with the full model with less complications. We follow the derivation in Kippenhahn, Weigert, and Weiss in Stellar Structure and Evolution.

Nearly all of the interior part of convective stars has an adiabatic stratification (corrections to this are small for fully convective regions), such that

${\displaystyle \frac{\delta \ln T}{\delta \ln P}={\mathrm{\nabla }}_{\text{adiabatic}}=0.4}$, which holds for an adiabatic expansion of an ideal gas.

We assume that this relation holds from the interior to the surface of the star—the surface is called photosphere. We assume ${\displaystyle {\mathrm{\nabla }}_{\text{adiabatic}}}$ to be constant throughout the interior of the star with value 0.4. However, we obtain the correct distinctive behavior.

For the interior we consider a simple polytropic relation between P and T:

${\displaystyle P=C{T}^{(1+n)}}$

With the index ${\displaystyle n=3/2}$.

We assume the relation above to hold until the photosphere where we assume to have a simple absorption law

${\displaystyle \kappa ={\kappa }_0{P}^a{T}^b}$

Then, we use the hydrostatic equilibrium equation and integrate it with respect to the radius to give us

${\displaystyle P_0=\text{Constant}\ast {\left(\frac{M}{R^2}{T}_{\text{eff}}^{-b}\right)}^{\frac{1}{1+a}}}$

For the solution in the interior we set ${\displaystyle P=P_0}$ ; ${\displaystyle T=T_{\text{eff}}}$ in the P-T relation and then eliminate pressure of this equation. Luminosity is given by the Stephan-Boltzmann law applied to a perfect black body:

${\displaystyle L=4\pi R^2\sigma T_{\text{eff}}^4}$.

Thus, any value of R corresponds to a certain point in the Hertzsprung–Russell diagram.

Finally, after some algebra this is the equation for the Hayashi limit in the Hertzsprung–Russell diagram:

${\displaystyle \log (T_{\text{eff}})=A\log (L)+B\log (M)+\text{constant}}$

With coefficients

${\displaystyle A=\frac{0.75a-0.25}{b-5.5a+1.5}}$, ${\displaystyle B=\frac{0.5a-1.5}{b-5.5a+1.5}}$

Takeaways from plugin in ${\displaystyle a\approx 1}$ and ${\displaystyle b\approx 3}$ for a cool hydrogen ion dominated atmosphere opacity model (${\displaystyle T<5000\text{ K}}$):

• The Hayashi limit must be far to the right in the Hertzsprung–Russell diagram which means temperatures have to be low.
• The Hayashi limit must be very steep. The gradient of Luminosity with respect to temperature has to be large.
• The Hayashi limit shifts slightly to the left in the Hertzsprung–Russell diagram for increasing M.

These predictions are supported by numerical simulations of stars.