<p><strong>TL; DR Somewhere between now and a few hundred billion years time. (For a co-moving volume)</strong> Now read on.</p>

<p>If stellar remnants are included, then the answer is very far in the future indeed, if and when the constituents of baryons begin to decay. So let's assume  that "stars" means those things that are undergoing nuclear fusion reactions to power their luminosity. Let's further assume that the stellar mass function, (<span class="math-container">$N(m)$</span> is the number of stars per unit mass) we see in the neighbourhood of the Sun is representative of populations in all galaxies at all times (difficult to make a start, without assuming this).</p>

<p>The number of stars that have been born is equal to the sum over time (the integral) and over mass of <span class="math-container">$N(m)$</span> multiplied by the rate at which mass is turned into stars in a comoving volume of the universe <span class="math-container">$\Phi(t)$</span>.</p>

<p>We then need to subtract a sum over time and mass of the rate of stellar death in the same comoving volume. The rate of stellar death is the rate of stellar birth at a time <span class="math-container">$t-\tau(m)$</span>, where <span class="math-container">$\tau(m)$</span> is the mass-dependent stellar lifetime. We ignore mass transfer in binary systems and assume that multiples can be treated as independent stellar components. </p>

<p>Thus the number of stars at time <span class="math-container">$t$</span> is approximately
<span class="math-container">$$N_*(t) = \int_0^{t} \int_m N(m) \Phi(t') - N(m)\Phi(t'-\tau(m))\ dm\ dt'\ .$$</span>
To find where this is a maximum, we differentiate with respect to time and then equate to zero. i.e. We look for the time when the stellar birth and death rates are the same.</p>

<p>I was going to (and possibly still will) attempt some sort of analytic approximation, but <a href="https://arxiv.org/abs/1403.0007" rel="noreferrer">Madau &amp; Dickinson (2014)</a> have done it better and taken into account the metallicity dependence of stellar lifetimes and the chemical evolution of galaxies. The star formation rate peaked about 10 billion years ago, is more than an order of magnitude lower now and is exponentially decreasing with a time constant of 3.9 billion years.</p>

<p>The integrated stellar mass is shown in their Fig 11 (shown below). It is still increasing today, but at a very low rate and has <em>not</em> passed through a maximum. The reason for this is that most stars have masses of 0.2-0.3 solar masses and lifetimes much longer than the age of the universe. Even if these stars are added at a very slow rate, their death rate is zero at present.</p>

<p><a href="https://i.sstatic.net/vHOvT.png" rel="noreferrer"><img src="https://i.sstatic.net/vHOvT.png" alt="Integrated stellar mass density (from Madau &amp; Dickinson 2014)"></a></p>

<p>If star formation did continue at a low-level then the number of stars would only begin to significantly diminish once the stars near the peak of the stellar mass function, that were born at the earliest times, start to die. The lifetime of a 0.25 solar mass star is around a trillion years (<a href="https://iopscience.iop.org/article/10.1086/304125" rel="noreferrer">Laughlin et al. 1997</a>). </p>

<p>On the other hand if star formation ceased <em>now</em> then the number of stars would immediately begin to diminish.</p>

<p>Perhaps we could argue that the current exponential decline will continue and the peak 
will come in another few billion years when stars of 0.8-0.9 solar masses begin dying off. However, that is futurology given that we have no first principles theory that explains the time-dependence of star formation, so I believe the best answer that can be given is somewhere between now and a few hundred billion years time.</p>

<p>Note that this answer assumes a co-moving volume. If the question asked is phrased in terms of the observable universe then because the number of stars has nearly reached a plateau, then the answer becomes close to whatever age the volume of the observable universe is maximised. I say "close to" because you have to factor in that the observable universe includes stars in distance slices at all cosmic epochs. I am unwilling to undertake this horrendous calculation, but note that the current concordance cosmological model has our observable universe slowly increasing from around a radius of 45 billion light years now, to about 60 billion light years in the far future <a href="https://arxiv.org/abs/astro-ph/0310808" rel="noreferrer">Davis &amp; Lineweaver 2005</a>, and this may compensate for a slow decline in the number of stars in a co-moving volume.</p>

CC BY-SA 4.0
2019-05-06T08:14:05.387Z
okay this will take some substantial time to read and think about, a few billion years at least. Thank you!
uhoh
CC BY-SA 4.0
2019-05-06T16:53:23.890Z
"that have (are, or will) undergone nuclear fusion reactions to power their luminosity" This confuses me - does this not include "stellar remnants"?
Keith McClary
CC BY-SA 4.0
2019-05-06T17:27:14.737Z
@KeithMcClary rewrite underway...
ProfRob
CC BY-SA 4.0
2019-05-06T17:51:27.270Z
So we're doomed, *Doomed* ..... ! :-)
StephenG - Help Ukraine
CC BY-SA 4.0
2019-05-07T08:03:16.540Z
This answer fits with what I got when I did a Monte Carlo simulation using the SFR and a lifespan model, ending up with peak star in about 10 billion years. It varies a lot between galaxies and neighbourhoods of course. The main uncertainty is infalling gas mass.
Anders Sandberg
CC BY-SA 4.0
2019-05-06T08:11:19.717Z
2019-05-06T17:56:41.453Z
2019-05-06T17:56:41.453Z
ProfRob
<p>There are very roughly a "mole" of stars in the universe. <a href="https://en.wikipedia.org/wiki/Star" rel="noreferrer">Wikipedia</a> quotes an estimate of <span class="math-container">$3 \times 10^{23}$</span> though the number is associated with some debate and uncertainly.</p>

<p>I'd like to know if there are estimates of when the number of stars in the universe will maximize. Is it expected to increase asymptotically to some maximum, or will it peak and then decrease.</p>

<p>I suppose this could depend on what the definition of "star" is taken to be, if brown or black dwarf objects are counted or not. I don't want to pre-specify, it's more likely that a good, well-informed answer will include this information.</p>

CC BY-SA 4.0
2019-05-06T04:00:22.027Z
Good question! The definition will also need to consider whether to include stellar remnants – and it's not a trivial consideration, firstly by proportion (e.g. of the 100 closest stars, 8 are white dwarfs = 8%), and secondly by cumulative effect (most stars end up as stellar remnants). If you include them, then finding a maximum is likely to be surprisingly complex, as it will need to take into account the number of binaries that end up merging, the proportion of SNEs that leave no compact body remnant, the number of stars and compact bodies that fall into SMBHs...
Chappo Hasn't Forgotten
CC BY-SA 4.0
2019-05-06T04:07:36.463Z
... and the trajectory of accelerating expansion, i.e. at what far-distant point in time does space expand fast enough to prevent a molecular cloud fragment from collapsing, and thus no new stars can be born?
Chappo Hasn't Forgotten
CC BY-SA 4.0
2019-05-06T07:48:42.230Z
It actually mainly depends on what you define as the "universe" and "when". I have written an answer for the number of stars in a co-moving volume as a function of cosmic epoch. The answer is considerably trickier for the number of stars in the observable universe, which is what wikipedia is talking about (with order of magnitude uncertainties actually).
ProfRob
CC BY-SA 4.0
2019-05-06T08:07:07.997Z
@RobJeffries if there is an estimate out there using any supportable definition of "universe" and "when", that would be just fine. As an aside, I'd thought there was at least some handle on the size/mass of the *whole universe* based on what is observable and models, but apparently not.
uhoh
CC BY-SA 4.0
2019-05-06T21:34:34.547Z
"I suppose this could depend on what the definition of "star" is taken to be" Actually, I think that the question of what "universe" means may be more of an issue. Does it means "everything in existence"? Does it mean "everything that is, as of time t, in the observable universe"? Does it mean "everything that is in space that is *now* in the observable universe"?
Acccumulation
CC BY-SA 4.0
2019-05-07T00:42:56.717Z
@Acccumulation hmm... I've just asked [Are there ways to estimate size of the "whole universe"?](https://astronomy.stackexchange.com/q/31793/7982)
uhoh
<cosmology><stellar-evolution><galactic-dynamics><fate-of-universe><modeling>
When will the number of stars be a maximum?
uhoh
# When will the number of stars be a maximum?

(Asked by: username_0 on 2019-05-06)

There are very roughly a "mole" of stars in the universe. [Wikipedia](https://en.wikipedia.org/wiki/Star) quotes an estimate of $3 \times 10^{23}$ though the number is associated with some debate and uncertainly.

I'd like to know if there are estimates of when the number of stars in the universe will maximize. Is it expected to increase asymptotically to some maximum, or will it peak and then decrease.

I suppose this could depend on what the definition of "star" is taken to be, if brown or black dwarf objects are counted or not. I don't want to pre-specify, it's more likely that a good, well-informed answer will include this information.

**Tags:** cosmology, stellar-evolution, galactic-dynamics, fate-of-universe, modeling

> **Comments:**
> **username_1**: Good question! The definition will also need to consider whether to include stellar remnants – and it's not a trivial consideration, firstly by proportion (e.g. of the 100 closest stars, 8 are white dwarfs = 8%), and secondly by cumulative effect (most stars end up as stellar remnants). If you include them, then finding a maximum is likely to be surprisingly complex, as it will need to take into account the number of binaries that end up merging, the proportion of SNEs that leave no compact body remnant, the number of stars and compact bodies that fall into SMBHs...  
> **username_1**: ... and the trajectory of accelerating expansion, i.e. at what far-distant point in time does space expand fast enough to prevent a molecular cloud fragment from collapsing, and thus no new stars can be born?  
> **username_2**: It actually mainly depends on what you define as the "universe" and "when". I have written an answer for the number of stars in a co-moving volume as a function of cosmic epoch. The answer is considerably trickier for the number of stars in the observable universe, which is what wikipedia is talking about (with order of magnitude uncertainties actually).  
> **username_0**: @RobJeffries if there is an estimate out there using any supportable definition of "universe" and "when", that would be just fine. As an aside, I'd thought there was at least some handle on the size/mass of the *whole universe* based on what is observable and models, but apparently not.  
> **username_3**: "I suppose this could depend on what the definition of "star" is taken to be" Actually, I think that the question of what "universe" means may be more of an issue. Does it means "everything in existence"? Does it mean "everything that is, as of time t, in the observable universe"? Does it mean "everything that is in space that is *now* in the observable universe"?  
> **username_0**: @username_3 hmm... I've just asked [Are there ways to estimate size of the "whole universe"?](https://astronomy.stackexchange.com/q/31793/7982)  

---

## 1 answers

### Answer by username_2 on 2019-05-06. Score: 21 | Accepted answer

**TL; DR Somewhere between now and a few hundred billion years time. (For a co-moving volume)** Now read on.

If stellar remnants are included, then the answer is very far in the future indeed, if and when the constituents of baryons begin to decay. So let's assume that "stars" means those things that are undergoing nuclear fusion reactions to power their luminosity. Let's further assume that the stellar mass function, ($N(m)$ is the number of stars per unit mass) we see in the neighbourhood of the Sun is representative of populations in all galaxies at all times (difficult to make a start, without assuming this).

The number of stars that have been born is equal to the sum over time (the integral) and over mass of $N(m)$ multiplied by the rate at which mass is turned into stars in a comoving volume of the universe $\Phi(t)$.

We then need to subtract a sum over time and mass of the rate of stellar death in the same comoving volume. The rate of stellar death is the rate of stellar birth at a time $t-\tau(m)$, where $\tau(m)$ is the mass-dependent stellar lifetime. We ignore mass transfer in binary systems and assume that multiples can be treated as independent stellar components.

Thus the number of stars at time $t$ is approximately
$$N\_\*(t) = \int\_0^{t} \int\_m N(m) \Phi(t') - N(m)\Phi(t'-\tau(m))\ dm\ dt'\ .$$
To find where this is a maximum, we differentiate with respect to time and then equate to zero. i.e. We look for the time when the stellar birth and death rates are the same.

I was going to (and possibly still will) attempt some sort of analytic approximation, but [Madau & Dickinson (2014)](https://arxiv.org/abs/1403.0007) have done it better and taken into account the metallicity dependence of stellar lifetimes and the chemical evolution of galaxies. The star formation rate peaked about 10 billion years ago, is more than an order of magnitude lower now and is exponentially decreasing with a time constant of 3.9 billion years.

The integrated stellar mass is shown in their Fig 11 (shown below). It is still increasing today, but at a very low rate and has *not* passed through a maximum. The reason for this is that most stars have masses of 0.2-0.3 solar masses and lifetimes much longer than the age of the universe. Even if these stars are added at a very slow rate, their death rate is zero at present.

[![Integrated stellar mass density (from Madau & Dickinson 2014)](https://i.sstatic.net/vHOvT.png)](https://i.sstatic.net/vHOvT.png)

If star formation did continue at a low-level then the number of stars would only begin to significantly diminish once the stars near the peak of the stellar mass function, that were born at the earliest times, start to die. The lifetime of a 0.25 solar mass star is around a trillion years ([Laughlin et al. 1997](https://iopscience.iop.org/article/10.1086/304125)).

On the other hand if star formation ceased *now* then the number of stars would immediately begin to diminish.

Perhaps we could argue that the current exponential decline will continue and the peak
will come in another few billion years when stars of 0.8-0.9 solar masses begin dying off. However, that is futurology given that we have no first principles theory that explains the time-dependence of star formation, so I believe the best answer that can be given is somewhere between now and a few hundred billion years time.

Note that this answer assumes a co-moving volume. If the question asked is phrased in terms of the observable universe then because the number of stars has nearly reached a plateau, then the answer becomes close to whatever age the volume of the observable universe is maximised. I say "close to" because you have to factor in that the observable universe includes stars in distance slices at all cosmic epochs. I am unwilling to undertake this horrendous calculation, but note that the current concordance cosmological model has our observable universe slowly increasing from around a radius of 45 billion light years now, to about 60 billion light years in the far future [Davis & Lineweaver 2005](https://arxiv.org/abs/astro-ph/0310808), and this may compensate for a slow decline in the number of stars in a co-moving volume.

> **Comments:**
> **username_0**: okay this will take some substantial time to read and think about, a few billion years at least. Thank you!  
> **username_4**: "that have (are, or will) undergone nuclear fusion reactions to power their luminosity" This confuses me - does this not include "stellar remnants"?  
> **username_2**: @username_4 rewrite underway...  
> **username_5**: So we're doomed, *Doomed* ..... ! :-)  
> **username_6**: This answer fits with what I got when I did a Monte Carlo simulation using the SFR and a lifespan model, ending up with peak star in about 10 billion years. It varies a lot between galaxies and neighbourhoods of course. The main uncertainty is infalling gas mass.  

---

