This MODULE implements the Chandrasekhar's degenerate \(\mathrm{EOS}\) for white dwarfs.
See Benz W., Bowers R.L., Cameron A.G.W., Press W.H., 1990, APJ, 348, 647., eqs.(2.4)-(2.5)
FT 19.12.2022
| Type | Visibility | Attributes | Name | Initial | |||
|---|---|---|---|---|---|---|---|
| double precision, | private, | parameter | :: | a_wd | = | 6.02D+21/(c_light2*cm2m)*dens_si2cu |
Constant with dimensions of a pressure in code units |
| double precision, | private, | parameter | :: | a_wd_cgs | = | 6.02D+21*press_si2cgs |
Constant with dimensions of a pressure in CGS units |
| double precision, | private, | parameter | :: | b_wd | = | mu_e*9.82D+8*dens_si2cu |
Constant with dimensions of a density in code units |
| double precision, | private, | parameter | :: | b_wd_cgs | = | mu_e*9.82D+8*kg2g/(m2cm**3) |
Constant with dimensions of a density in CGS units |
| double precision, | private, | parameter | :: | dens_si2cu | = | 1.618654158231174D-21 |
Conversion factor for the baryon mass density, from SI units to code units |
| double precision, | private, | parameter | :: | mu_e | = | 2.D0 |
Mean molecular weight per electron |
Function used in the computation of the degenerate pressure
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| double precision, | intent(in) | :: | x |
Function used in the computation of the degenerate specific internal energy
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| double precision, | intent(in) | :: | x |
Degenerate pressure as a function of density
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| double precision, | intent(in) | :: | rho |
Degenerate pressure as a function of density, all in CGS units
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| double precision, | intent(in) | :: | rho |
Degenerate density as a function of pressure
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| double precision, | intent(in) | :: | pr | |||
| double precision, | intent(in) | :: | rho_left_bracket | |||
| double precision, | intent(in) | :: | rho_right_bracket |
Degenerate specific internal energy as a function of density
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| double precision, | intent(in) | :: | rho |
Test that the implementation is correct. Takes the density in code units as input, and assigns CGS values to the pressure, and the density recomputed from the pressure. It also computed the dimensionless specific internal energy.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| double precision, | intent(in) | :: | rho_input | |||
| double precision, | intent(out) | :: | rho | |||
| double precision, | intent(out) | :: | pr | |||
| double precision, | intent(out) | :: | u |