Celestia icon indicating copy to clipboard operation
Celestia copied to clipboard

Proper WR temperatures

Open Trolligi112477-ww opened this issue 3 years ago • 3 comments

As you may know, I am obsessed with Wolf-Rayet stars. I'm a bit frustrated whenever I play around with WRs and see that they're so much cooler than they should be. For example, HIP 5100 (WR 2) in Celestia has a temperature of 50,000 K, while in reality that temperature is over 140,000 K.

I would like the Wolf-Rayet temperatures to be raised, to reflect their real-life temperatures. Note: This is for the hydrogen-free Wolf-Rayet stars. For the hydrogen-rich ones (WNh stars) I would like a separate class to be created that has temperatures that reflect the WNh stars.

Sometimes, I use the hottest white dwarf (D0) which has a temperature of 100,000 K, to depict hot Wolf-Rayet stars. This works to an extent, but not very well for WRs between 60kK and 100kK, and ones between 110kK and 130kK.

Possible temperatures for WRs: (Note: The temperatures for early and late WCs are mostly for CSPNe with [WC] spectra) static float tempWN[10] = { 185000, 158000, 140000, 115000, 89000, 75000, 63000, 54000, 47000, 39000 };

static float tempWC[12] = { 225000, 190000, 165000, 140000, 125000, 105000, 83000, 70000, 58000, 46000, 40000, 34000 };

static float tempWNh[12] = { 133000, 106000, 92000, 72000, 61000, 51000, 42000, 36000, 31000, 27000, 22000, 15000 };

and because there are the 10 and 11 subclasses,

const char* SubclassNames[11] = { "0", "1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "" };

const char* SpectralClassNames[StellarClass::NormalClassCount] = { "O", "B", "A", "F", "G", "K", "M", "R", "S", "N", "WC", "WN", "WNh", "WO", "?", "L", "T", "Y", "C", };

Trolligi112477-ww avatar Dec 28 '20 12:12 Trolligi112477-ww

An update as now I know more about WRs.

static float tempWN[11] = { 220000, 180000, 140000, 100000, 80000, 70000, 60000, 50000, 42000, 35000, 30000 };

static float tempWC[13] = { 225000, 190000, 165000, 140000, 110000, 85000, 78000, 70000, 60000, 45000, 32000, 25000, 20000 };

static float tempWNh[12] = { 200000, 160000, 140000, 85000, 70000, 60000, 50000, 45000, 40000, 35000, 30000, 25000 };

and because there are the 10 and 11 subclasses,

const char* SubclassNames[11] = { "0", "1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "" };

const char* SpectralClassNames[StellarClass::NormalClassCount] = { "O", "B", "A", "F", "G", "K", "M", "R", "S", "N", "WC", "WN", "WNh", "WO", "?", "L", "T", "Y", "C", };

Note that WC0, WC1, WC2, WC3, WC10, WC11 and WC12 are only for CSPNe and WN0 and WN1 (and consequently WN0h and WN1h) aren't real. Also, WO1 should be 220,000 K.

Trolligi112477-ww avatar Nov 20 '21 14:11 Trolligi112477-ww

So can you update the current pull request or make a new one, with sources?

SevenSpheres avatar Nov 20 '21 18:11 SevenSpheres

My sources are: https://ui.adsabs.harvard.edu/abs/2019A%26A...625A..57H/abstract https://ui.adsabs.harvard.edu/abs/2019A%26A...621A..92S/abstract https://ui.adsabs.harvard.edu/abs/2020A%26A...640A..10W/abstract (for WC10-12)

For the WC0-3, WN0-1 and WO1 temperatures, I just assumed that they would continue to climb, as it would be hard to use the CSPNe as a source for them, because it seems that early-type CSPNe spectral types are at best weakly correlated with temperatures.

Trolligi112477-ww avatar Dec 30 '21 19:12 Trolligi112477-ww