<p>Whether or not the object is just above or below the threshold (about 13 Jupiter masses) at which the core will become hot enough to burn deuterium is not really relevant to the calculation you wish to understand. The object has cooled beyond the point at which D burning has ceased. It has also pretty much contracted to a minimum sized configuration where electron degeneracy pressure supports it. Therefore very little gravitational potential energy can contribute to subsequent thermal evolution.</p>

<p>All that is left is to calculate the thermal energy content (roughly $3kT/2 \times N_i$ - the number of hydrogen and helium ions, where $T$ is the <em>internal</em> temperature; the degenerate electrons contribute very little to the heat capacity) and then divide this by the luminosity radiated from the surface. The calculation is a little more complex than that because the luminosity fades as the brown dwarf cools.</p>

<p>Just for kicks - here is a back of the envelope calculation. The central temperature of a contracting ball of gas is given by the virial theorem. The brown dwarf you are talking about will contract until D burning (nearly) starts at about $10^{6}$ K - let's assume that is an upper limit to the current core temperature. The average temperature in the interior will be a bit lower than that (let's divide by 2), so $\langle T \rangle \sim 5\times 10^{5}$ K.</p>

<p>If the mass is about $0.012 M_{\odot}$ and it is made almost totally of hydrogen, then the number of protons is about $N_i \sim 1.4 \times 10^{55}$.</p>

<p>Thus the total thermal energy is $E \sim 3N_i k\langle T \rangle/2 \sim 1.5\times 10^{38}$ J.</p>

<p>The <em>current</em> luminosity $L \sim 4 \pi \sigma R^2 T_s^{4}$, where $T_s \simeq 1100$ K is the <em>surface</em> temperature and $R$ is the radius, which will be about the same as Jupiter, i.e. $R\sim 70,000$ km. Thus $L \sim 5 \times 10^{21}$ W. </p>

<p>The current cooling timescale is $E/L \sim 10^{9}$ years - but this is a <em>lower limit</em>.</p>

<p>I am not surprised that this is somewhat lower than your estimate because the cooling time gets much longer as the object cools. This is clearly seen when looking at a proper model calculation; for example see the picture below taken from <a href="https://arxiv.org/abs/astro-ph/9705201" rel="nofollow noreferrer">Burrows et al. (1997)</a>.</p>

<p>The lines show the cooling tracks of surface temperature versus (log) time. The red curves are "planets", the green are "brown dwarfs" - the dividing line is the D-burning threshold - seen as a plateau in the the green brown dwarf curves at early times. Your object is around this threshold. I have marked two points on the plot - the first corresponding to a 13 Jupiter mass object with a surface temperature of 1100K (you can see it must already be at least a hundred million years old) , the second marks the same object when it has cooled to about 273 K (zero Celsius). It does indeed occur after about 10 billion years.</p>

<p><a href="https://i.sstatic.net/eRXBi.jpg" rel="nofollow noreferrer"><img src="https://i.sstatic.net/eRXBi.jpg" alt="Burrows cooling curves"></a></p>

CC BY-SA 3.0
2017-05-28T17:19:02.347Z
Also once the electron degeneracy is at its limit, the electrons are tightly confined to their orbitals so they can no longer transition to higher orbitals and back to release photons.  This would effectively stop heat radiation, thereby drastically slowing down thermal losses.  Thus substantially increasing the cool down period.
0tyranny0poverty
CC BY-SA 3.0
2017-05-28T18:41:08.720Z
@0tyranny0poverty This is quite incorrect. The material near the surface where the radiation arises is not in a degenerate state and is not even ionised. The effect of electron degeneracy in the interior is to reduce the thermal energy because the electrons have an extremely low heat capacity, which decreases the cooling time. That is why I only mentioned heat stored by the ions in my answer.
ProfRob
CC BY-SA 3.0
2017-05-28T17:09:18.400Z
2017-05-28T19:38:16.677Z
2017-05-28T19:38:16.677Z
ProfRob
<p>Brown dwarves usually exhibit deuterium burning in their core, which is where this omnious limit of 13 Jupiter masses comes from. It is the limit at which the core is dense and hot enough (as others have pointed out) enough to make deuterium nuclei fuse together. Deuterium has a relative abundance compared with Hydrogen of $10^{-4}$, but is easier to fuse than regular hydrogen. This is why brown dwarves can exist, after all.</p>

<p>The energy flux they generate in the core is slowly diffused outwards until it reaches the shell from which it can escape. The virial theorem then tells us, that exactly half of their luminosity we see on the surface is retained in the dwarf. Thus, is we know an energy content $E$ and we know the energy loss rate $\dot E$ (which is the luminosity) astronomers usually go ahead to construct an estimate of the e-folding cooling &amp; contraction time as $t_{\text{cool}} = E/\dot E$.</p>

<p>This can be used as a first estimate of the future fate of the brown dwarf. The estimate will now change based on whether the object is actually a brown dwarf that fuses deuterium, or not. This is easily understandable because above the deuterium burning limit the energy loss rate will be something like $\dot E = \dot E_{\text{contraction}} + \dot E_{\text{deuterium burning}}$, while below that limit the deuterium burning term is zero, and thus the estimate for $t_{\text{cool}}$ increases.</p>

CC BY-SA 3.0
2017-05-26T19:00:42.593Z
Is deuterium formed in the Big Bang along with Hydrogen, Helium, and Lithium?  I thought deuterium formed as a product of regular Protium (Hydrogen) fusion.  Doesn't fusion require the deuterium to be 99.9% pure?  How would it work if only 1 of 10,000 hydrogen atoms are deuterium?
0tyranny0poverty
CC BY-SA 3.0
2017-05-27T03:13:55.280Z
@0tyranny0poverty: So many questions, I'll answer them in order: Yes, and that at roughly the  fraction that I've cited. Also. No, only in low-density, high temperature plasmas that humans use on Earth. The wiki-article on deuterium fusion gives a bit more detail on that.
AtmosphericPrisonEscape
CC BY-SA 3.0
2017-05-28T17:14:27.040Z
How do you equate half the luminosity to an energy content? They are different quantities measured in different units. Both terms on the RHS of the energy loss equation you write are actually *heating* terms. The only cooling term is radiation  from the surface.
ProfRob
CC BY-SA 3.0
2017-05-29T12:53:18.053Z
@RobJeffries: Errr you're right. Seems I wrote this in a half-asleep state. I'll do corrections.
AtmosphericPrisonEscape
CC BY-SA 3.0
2017-05-26T14:13:53.297Z
2017-05-28T05:33:44.650Z
2017-05-28T05:33:44.650Z
AtmosphericPrisonEscape
<p>The surface temperature of a brown dwarf (or a star, for that matter) is considerably less than the core temperature.  As <a href="http://www.scholarpedia.org/article/Brown_dwarfs" rel="nofollow noreferrer">scholarpedia</a> writes, </p>

<blockquote>
  <p>Depending on the mass of the brown dwarf, core temperatures can reach
  10^4 to 6×10^6 K, while core densities range from 10 to 10^3 g/cm^3</p>
</blockquote>

<p>There's a lot more energy in the body than if the entire dwarf were at 1100 K. Even after downgrading to planet status, the core temperature is maybe a decade greater than the surface temperature.</p>

CC BY-SA 3.0
2017-05-26T12:08:12.913Z
2017-05-26T12:08:12.913Z
Carl Witthoft
<p>I was reading <a href="http://newatlas.com/brown-dwarf-is-a-planet/49445/" rel="nofollow noreferrer">this</a> article and it said something that didn't sit right with me.</p>

<p>"SIMP0136 still gives off light mostly in the infrared wavelengths, as its temperature is now approximately 830° C (1,526° F)," Gagné told us. "It will take it approximately <strong>10 billion years more before it reaches 0° C</strong>"</p>

<p>SIMP0136 was classified as a brown dwarf star but when they just recently figured out its mass to be only 12.7 times that of Jupiter they downgraded it to a gas giant planet.  The threshold to qualify as a brown dwarf is 13 Jupiter masses.</p>

<p>Like gas giant planets, brown dwarf stars do not gave enough mass for fusion to occur so the term star is a bit misleading.  Since there is no fusion, then it should eventually cool down.  The 10 billion year cool down time to go from its present 830 deg C to 0 deg C seems way too long.  How do they calculate this? </p>

CC BY-SA 3.0
2017-05-26T03:19:55.180Z
Consider how large it is and how slow it takes for something that size to cool off and how much heat was created at it's formation.  It's probably close to a million degrees in it's core and its over 5,000 times the mass of the Earth.
userLTK
<astrophysics><star-formation><gas-giants><brown-dwarf><thermodynamics>
How do we calculate heat flow or cooling rate for no fusion brown dwarf star?
0tyranny0poverty
# How do we calculate heat flow or cooling rate for no fusion brown dwarf star?

(Asked by: username_0 on 2017-05-25)

I was reading [this](http://newatlas.com/brown-dwarf-is-a-planet/49445/) article and it said something that didn't sit right with me.

"SIMP0136 still gives off light mostly in the infrared wavelengths, as its temperature is now approximately 830° C (1,526° F)," Gagné told us. "It will take it approximately **10 billion years more before it reaches 0° C**"

SIMP0136 was classified as a brown dwarf star but when they just recently figured out its mass to be only 12.7 times that of Jupiter they downgraded it to a gas giant planet. The threshold to qualify as a brown dwarf is 13 Jupiter masses.

Like gas giant planets, brown dwarf stars do not gave enough mass for fusion to occur so the term star is a bit misleading. Since there is no fusion, then it should eventually cool down. The 10 billion year cool down time to go from its present 830 deg C to 0 deg C seems way too long. How do they calculate this?

**Tags:** astrophysics, star-formation, gas-giants, brown-dwarf, thermodynamics
> **Comments:**
> **username_1**: Consider how large it is and how slow it takes for something that size to cool off and how much heat was created at it's formation.  It's probably close to a million degrees in it's core and its over 5,000 times the mass of the Earth.  

---

## 3 answers

### Answer by username_2 on 2017-05-28. Score: 5

Whether or not the object is just above or below the threshold (about 13 Jupiter masses) at which the core will become hot enough to burn deuterium is not really relevant to the calculation you wish to understand. The object has cooled beyond the point at which D burning has ceased. It has also pretty much contracted to a minimum sized configuration where electron degeneracy pressure supports it. Therefore very little gravitational potential energy can contribute to subsequent thermal evolution.

All that is left is to calculate the thermal energy content (roughly $3kT/2 \times N\_i$ - the number of hydrogen and helium ions, where $T$ is the *internal* temperature; the degenerate electrons contribute very little to the heat capacity) and then divide this by the luminosity radiated from the surface. The calculation is a little more complex than that because the luminosity fades as the brown dwarf cools.

Just for kicks - here is a back of the envelope calculation. The central temperature of a contracting ball of gas is given by the virial theorem. The brown dwarf you are talking about will contract until D burning (nearly) starts at about $10^{6}$ K - let's assume that is an upper limit to the current core temperature. The average temperature in the interior will be a bit lower than that (let's divide by 2), so $\langle T \rangle \sim 5\times 10^{5}$ K.

If the mass is about $0.012 M\_{\odot}$ and it is made almost totally of hydrogen, then the number of protons is about $N\_i \sim 1.4 \times 10^{55}$.

Thus the total thermal energy is $E \sim 3N\_i k\langle T \rangle/2 \sim 1.5\times 10^{38}$ J.

The *current* luminosity $L \sim 4 \pi \sigma R^2 T\_s^{4}$, where $T\_s \simeq 1100$ K is the *surface* temperature and $R$ is the radius, which will be about the same as Jupiter, i.e. $R\sim 70,000$ km. Thus $L \sim 5 \times 10^{21}$ W.

The current cooling timescale is $E/L \sim 10^{9}$ years - but this is a *lower limit*.

I am not surprised that this is somewhat lower than your estimate because the cooling time gets much longer as the object cools. This is clearly seen when looking at a proper model calculation; for example see the picture below taken from [Burrows et al. (1997)](https://arxiv.org/abs/astro-ph/9705201).

The lines show the cooling tracks of surface temperature versus (log) time. The red curves are "planets", the green are "brown dwarfs" - the dividing line is the D-burning threshold - seen as a plateau in the the green brown dwarf curves at early times. Your object is around this threshold. I have marked two points on the plot - the first corresponding to a 13 Jupiter mass object with a surface temperature of 1100K (you can see it must already be at least a hundred million years old) , the second marks the same object when it has cooled to about 273 K (zero Celsius). It does indeed occur after about 10 billion years.

[![Burrows cooling curves](https://i.sstatic.net/eRXBi.jpg)](https://i.sstatic.net/eRXBi.jpg)

> **Comments:**
> **username_0**: Also once the electron degeneracy is at its limit, the electrons are tightly confined to their orbitals so they can no longer transition to higher orbitals and back to release photons.  This would effectively stop heat radiation, thereby drastically slowing down thermal losses.  Thus substantially increasing the cool down period.  
> **username_2**: @0tyranny0poverty This is quite incorrect. The material near the surface where the radiation arises is not in a degenerate state and is not even ionised. The effect of electron degeneracy in the interior is to reduce the thermal energy because the electrons have an extremely low heat capacity, which decreases the cooling time. That is why I only mentioned heat stored by the ions in my answer.  

---

### Answer by username_3 on 2017-05-26. Score: 3

Brown dwarves usually exhibit deuterium burning in their core, which is where this omnious limit of 13 Jupiter masses comes from. It is the limit at which the core is dense and hot enough (as others have pointed out) enough to make deuterium nuclei fuse together. Deuterium has a relative abundance compared with Hydrogen of $10^{-4}$, but is easier to fuse than regular hydrogen. This is why brown dwarves can exist, after all.

The energy flux they generate in the core is slowly diffused outwards until it reaches the shell from which it can escape. The virial theorem then tells us, that exactly half of their luminosity we see on the surface is retained in the dwarf. Thus, is we know an energy content $E$ and we know the energy loss rate $\dot E$ (which is the luminosity) astronomers usually go ahead to construct an estimate of the e-folding cooling & contraction time as $t\_{\text{cool}} = E/\dot E$.

This can be used as a first estimate of the future fate of the brown dwarf. The estimate will now change based on whether the object is actually a brown dwarf that fuses deuterium, or not. This is easily understandable because above the deuterium burning limit the energy loss rate will be something like $\dot E = \dot E\_{\text{contraction}} + \dot E\_{\text{deuterium burning}}$, while below that limit the deuterium burning term is zero, and thus the estimate for $t\_{\text{cool}}$ increases.

> **Comments:**
> **username_0**: Is deuterium formed in the Big Bang along with Hydrogen, Helium, and Lithium?  I thought deuterium formed as a product of regular Protium (Hydrogen) fusion.  Doesn't fusion require the deuterium to be 99.9% pure?  How would it work if only 1 of 10,000 hydrogen atoms are deuterium?  
> **username_3**: @0tyranny0poverty: So many questions, I'll answer them in order: Yes, and that at roughly the  fraction that I've cited. Also. No, only in low-density, high temperature plasmas that humans use on Earth. The wiki-article on deuterium fusion gives a bit more detail on that.  
> **username_2**: How do you equate half the luminosity to an energy content? They are different quantities measured in different units. Both terms on the RHS of the energy loss equation you write are actually *heating* terms. The only cooling term is radiation  from the surface.  
> **username_3**: @RobJeffries: Errr you're right. Seems I wrote this in a half-asleep state. I'll do corrections.  

---

### Answer by username_4 on 2017-05-26. Score: 0

The surface temperature of a brown dwarf (or a star, for that matter) is considerably less than the core temperature. As [scholarpedia](http://www.scholarpedia.org/article/Brown_dwarfs) writes,

> Depending on the mass of the brown dwarf, core temperatures can reach
> 10^4 to 6×10^6 K, while core densities range from 10 to 10^3 g/cm^3

There's a lot more energy in the body than if the entire dwarf were at 1100 K. Even after downgrading to planet status, the core temperature is maybe a decade greater than the surface temperature.


---

