How can the Hertzsprung-Russell diagram be used in star building?

I've grabbed an observational Hertzsprung"“Russell diagram from the wikipedia article about the same.

In this case it's a plot of 22,000 stars and you can see that there's a lot more flexibility in where the stars lie, in practice a not insignificant number lie in some of the "black bits" of the simplified diagram you've posted. However there are still clear black areas.

posted over 6 years ago

Separatrix‭

11 reputation 5 1 1 0

To accurately answer your question, you might need to use a stellar evolution code, either doing your own modeling or looking up existing data tables. I'd recommend the MESA code for the former approach, and the Geneva grids for the latter (see Eggenberger et al. (2008) for details). Numerical simulations are excellent, and as the results are often available to the public, you can get some nice results with a bit of effort.

That said, simply using these doesn't tell you anything about the population as a whole - in other words, how the masses of stars are distributed. This can be computed easily using an initial mass function, or IMF. An IMF is something of the form $\xi(m)\Delta m$, which tells you how many stars were born with masses between $m$ and $m+\Delta m$.¹ This is usually in the form of a power law, i.e. $$\xi(m)\Delta m=\xi_0m^{-a}\Delta m$$ where $m$ is in solar masses and $\xi_0$ is a constant. Essentially, if you want to find out the number of stars which have masses between $m_1$ and $m_2$, you integrate: $$N(m_1,m_2)=\xi_0\int_{m_1}^{m_2}m^{-a}dm$$ where $\xi_0$ is a normalization constant such that $N(m_1,m_2)=N$ for a sample of $N$ stars. A common IMF is the Kroupa IMF, represented as a broken power law - in other words, $a$ is different for different mass ranges. This is really the answer to your question, if you're talking about the likelihood of finding a star in a given mass range.

I've written some code for this answer (see on Github) that generates a few plots. The first one is the Kroupa IMF, for a sample where $N=100$:

Here is the cumulative distribution function, the number of stars of mass less than $m$:

However, stars are born and die at different rates. We can divide up their lifetimes into four stages:

1. Pre-main sequence evolution²
2. Main sequence evolution
3. Post-main sequence evolution
4. Stellar remnant

If you want to find the time it takes a star to pass through each phase, you can use some simple analytical approximations. We can estimate timescales for pre-main sequence evolution ($\tau_{\text{PMS}}$) and main sequence evolution ($\tau_{\text{MS}}$). These give the amount of time a star spends in a stage as a function of its initial mass: $$\tau_{\text{PMS}}=10^7\left(\frac{M}{M_{\odot}}\right) ^{-2.5}\text{ years}, \quad\tau_{\text{MS}}=10^{10}\left(\frac{M}{M_{\odot}}\right)^{-2.5}\text{ years}$$ The second exponent is often picked between $-2.5$ and $-3$; I've chosen $-2.5$ as a conservative estimate. Then, seeing that $\tau_{\text{PMS}}\ll\tau_{\text{MS}}$, we can see that at time $t$ stars of a given minimum mass won't have reached the main sequence yet, and stars of a given maximum mass will have already evolved off of it. Therefore, the total number of main sequence stars at a time $t$ is $$N_{\text{MS}}(t)=\xi_0\int_{m_\text{min}}^{m_\text{max}}m^{-a}dm$$ where $$m_{\text{min}}(t)=\left(\frac{t}{10^7\text{ yrs}}\right)^{-2/5},\quad m_{\text{max}}(t)=\left(\frac{t}{10^{10}\text{ yrs}}\right)^{-2/5}$$ I've used these to create plots of the number of stars on the pre-main sequence, main sequence, and post-main sequence tracks over 100 billion years. Notice that even after 100 billion years, many of the stars still have not left the main sequence; they're low-mass red dwarfs.

The answer to your second question - can stars be found in the black regions - is, essentially, a yes. All of the H-R diagrams here show the two large areas above the main sequence populated by giants and supergiants, but stars evolve back and forth along that area. For instance, a star may change from being a red supergiant to a blue supergiant (or vice versa), often including periods of violent activity. The evolved stars above the main sequence fall in a number of places.

You're not likely to see quite as many stars below the main sequence - the notable exception being white dwarfs - but some stars, called subdwarfs, do lie in that area. However, they're usually closer to the main sequence than to the white dwarfs. I don't think any are clearly visible on the H-R diagrams above, but they do exist. Kapteyn's star is a good example of a fairly cool subdwarf.

If you want to get a better idea of which gaps are filled and which ones aren't, you can try looking at a star catalogue. Gaia recently released a slew of new data this spring, and you should be able to determine temperature and luminosity for many stars. Assuming you can query that (or another) database, perhaps you can find some stars in unusual places.

¹ Bear in mind that a star may lose mass during the course of its lifetime, often through stellar winds (or, in a select few cases, violent eruptions).
² You've heard of a protostar; a protostar is a pre-main sequence star that is still enshrouded by material from the molecular cloud that formed it, and is likely not visible.

posted over 6 years ago

HDE 226868‭

601 reputation 61 463 80 49

1. Yes, stars can exist in the "black" area of your diagram. Here is a more detailed diagram, and this only plots 22 000 stars. There are more than a million times more stars just in the Milky Way galaxy.

2. That depends entirely on what region of space you are looking at.

   And why is this even important? Are you going to be quoting statistics at your readers? You will make them fall asleep, just sayin'...

posted over 6 years ago

MichaelK‭

21 reputation 0 2 2 0