The Highly Accurate Relation Between the Radius and Mass of the White Dwarf Star From Zero to Finite Temperature

In this research, first considering the electron–electron interaction in the high-density Fermi electron gas at T = 0 K, this interaction causes the pressure 2/137 time less than the original value. However, the pressure of the Fermi electron gas should have something to do with temperature. Then, we estimate the temperature effect using statistical mechanics and find that the complicated form of the pressure p depends on temperature at the given particle number N and volume V. According to this, the central density–mass (ρ c -M), central density–radius (ρ c -R), and mass–radius (M-R) relations of the white dwarf star are obtained by considering the equation of state (EOS). Traditional formula gives the problematic mass–radius relation R ∝ M − 1 / 3 for the low-density white dwarf stars because it leads to R→∞ and p→0 when M→0. We correct this relation and obtain two reasonable relations in the relativistic and nonrelativistic regions. In our EOS calculations, the central density is divided into the high-, middle-, and low-density regions. All three relations are almost unchanged until 108 K in the high-density region. The temperature effect mainly affects the middle- and low-density regions, and it becomes explicitly above 107 K. Our calculations can explain Sloan Digital Sky Survey observations where some white dwarf stars with a radius of more than 8 × 103 km have larger mass than the predictions by the relativistic EOS at T = 0 K. This result tells us that the temperature effect is important for the low and middle central-density white dwarf star and also useful to estimate the inner temperature of a white dwarf star.


INTRODUCTION
The white dwarf star has been investigated for many years, and it was named first in 1922 (Holberg, 2005). It usually has a very high density with mass similar to our Sun, but the volume is small like the Earth. The Chandrasekhar mass limit is the famous restriction for the upper mass of a non-rotating and uncharged white dwarf star (Chandrasekhar, 1935;Chandrasekhar, 2012). The reported largest mass seems to be the one found in 2007, which is 1.33 times as large as the solar mass M ʘ (Kepler et al., 2007). The white dwarf star is thought to be the type of the low to medium mass stars in the final evolution stage (Schutz, 1985;Hans and Ruffini, 1994;Mould, 2002;Misner et al., 2017). The early theory to explain its mass upper limit is based on the ideally degenerate Fermi electron gas (Chandrasekhar, 1935;Kittel and Kroemer, 1980;Schutz, 1985;Huang, 1987;Greiner et al., 1995;Honerkamp, 2002;Schwabl, 2002;Chandrasekhar, 2012). The calculation adopts all electrons like free particles occupying all energy levels until Fermi energy as they are at zero temperature (Koester and Chanmugam, 1990). It is amazing that, even in the high-temperature and highpressure situation, the ideal Fermi gas still works to explain the existence of the white dwarf star. It makes us curious to further discuss the more detailed temperature effect through statistical mechanics.
Actually, the central temperature of the Fermi electron gas in the white dwarf star might be about 10 7 -10 8 K (Soares, 2017), the temperature effect should be further considered to get more accurate results. The problems about surface luminosity and temperature evolution of the white dwarf stars have been discussed (D'antonia and Mazzitelli, 1990).
Here, we care about the central temperature more. According to statistical mechanics, we first build the pressure produced by the Fermi electron gas at given temperature T, the number of electrons N, and the total volume V. Then, we discuss the temperature effect on the electron pressure. Next, the relations between the central density ρ c , mass M, and radius R of the white dwarf star are derived from the equation of state (EOS). Finally, on the basis of these relations, the temperature effect is discussed and the results are compared to the Sloan Digital Sky Survey (SDSS) observations. Although the temperature is maximum at the center and minimum at the surface, we still adopt the uniform temperature approximation as most research did. After all, it is more complicated to consider the temperature varying with position, and our results imply that the uniform temperature approximation may remain valuable. Furthermore, the temperature can change with time because the white dwarf stars cool down by thermal emission. However, the cooling time is much longer than the astronomical observations, and the temperature of the white dwarf star is kept constant at each calculation in our research. If we want to compare the calculation results with the astronomical observations in several decades, then to hold the initial temperature throughout each calculation is a good enough way.

THE DEGENERATE FERMI ELECTRON GAS FOR THE WHITE DWARF STAR
First, we review the calculation of the upper mass limit for the white dwarf star. It adopts the ideally degenerate Fermi electron gas and considers the relativistic kinetic energy in the calculation (Huang, 1987;Greiner et al., 1995). Because the electron has spin s ± 1 2 , each energy state permits two electrons occupied. Each electron has the rest mass m e , and its relativistic kinetic energy E k at momentum p is The Fermi electron gas with the total number N and total volume V has total kinetic energy where h is the Planck's constant and p F is the Fermi momentum defined as Considering a white dwarf star mainly consisting of helium nuclei, then the total mass M in terms of the mass m p of a proton and the mass m n of a neutron is If we define the parameter then Eq. 2 becomes where The pressure produced by the ideal Fermi electron gas is (Huang, 1987) It is almost 1,000 times larger than the pressure of the helium nuclei in the same white dwarf star (Greiner et al., 1995). Further discussions give the relation between the radius R and mass M of the star for the relativistically high-density Fermi electron gas where R 2πm e c h R, and In Eq. 12, G is the gravitational constant and δ is a parameter of pure number. Some considerations (Greiner et al., 1995) give the upper mass limit M 0 in unit of the mass M ʘ of our Sun which is also the upper limit for appearance of the white dwarf star at T 0 K.

THE CORRECTION OF THE ELECTRON-ELECTRON INTERACTION FOR THE WHITE DWARF STAR AT T = 0
The ideal Fermi electron gas has been widely discussed in solid state physics. The ground state energy of non-relativistically highdensity Fermi electron gas has been calculated by the Hartree-Fock approximation (Mattuck, 1976;Mahan, 2000), and the energy per electron at T 0 is where E HF is the total energy of the Fermi electron gas, and r s is defined by using the Bohr radius a B V N The first two terms are dominant terms, and the ratio of the first term to the second one is proportional rs or N 1/3 . As N increases, the first term increases faster than the second one. Actually, the calculation of the first term at the right-hand side in Eq. 14 should use Eq. 6 because of the relativistic electrons. Considering x F >>1 in the relativistic region, then Eq. 6 becomes The second term consider the Feynman diagram of the oyster type, so this correlation energy E 1 is (Mattuck, 1976;Mahan, 2000) where K e is the Coulomb's constant. Using Eqs. 16 and 17, the pressure P HF of the Fermi electron gas at T 0 K is where 2πK e e 2 /hc is the fine structure constant (Gasiorowicz, 1974;Levin, 2002;Roger, 2002;Gottfried and Yan, 2003). It means that the electron-electron interaction causes the pressure to be about 2/137 time less than the original value. Some related discussions can be checked in the early reference (Salpeter, 1961).

THE TEMPERATURE EFFECT ON THE PRESSURE OF THE IDEAL FERMI ELECTRON GAS IN THE WHITE DWARF STAR
The EOS for ideal and non-interacting electron gas has been given in many resources like Ref. 12. In reality, the central temperature of a star is usually about 10 7 -10 8 K, and the upper mass limit in Eq. 13 calculated at T 0 should be improved. Otherwise, it cannot reflect how the relation between the radius and mass of the white dwarf star varies with temperature. Then, we consider the case for T>>0, and the grand partition function in statistical mechanics (Greiner et al., 1995) is where E k is the kinetic energy, β 1/k B T, and z exp (μβ) with μ being the chemical potential of the Fermi electron gas. Because the energy eigenstates are treated as arbitrarily close to each other in a very large volume, the grand partition function becomes Integrating it by parts, then we have (Greiner et al., 1995) ln Z g 4πV h 3 β 3 ∞ 0 p 3 dp dE k dp where g 2s + 1 is the degeneracy factor and Substituting Eq. 22 into Eq. 21 and considering Fermi energy E F >> m e c 2 , it gives Using the Taylor series expansion to the second-order term, then we have Then, the integral gives where we define the function The corresponding Fermi energy E F is roughly 20 MeV (Huang, 1987) and 1/(2m e c 2 β) ∼ 1/1,000 at 10 7 K. The chemical potential μ∼E F , so z exp (βμ)∼exp (20,000). When z>>1, the approximation of Eq. 26 (Greiner et al., 1995) is and the ratio of the first term to the second term is about According to the relation ln Z pV/k B T, the pressure causing by the Fermi electron gas is Then, we calculate the particle number N (T, V, z) using a similar way in statistical mechanics. It gives Using Taylor series expansion to the second-order term, then we have Further calculation gives Combining Eq. 29 with Eq. 32, it gives P electron gas in terms of T, V, and N, that is, Substituting Eq. 27 into Eq. 33 and further rearrangement give In Eq. 34, (ln z) can be obtained by substituting Eq. 27 into Eq. 32. Then, we have In the following, we have to solve the important chemical potential for the ultra-relativistic Fermi electron gas. Most textbooks and references still use the chemical potential for the non-relativistic Fermi electron gas to deal with the ultrarelativistic Fermi electron problem, so we propose the appropriate chemical potential here. Considering (lnz)>>(m e c 2 /k B T) and using (lnz) μ/k B T with the definition of the Fermi momentum in Eq. 3, it further gives Finding the cubic roots at both sides in Eq. 36 and expanding the bracket to the second-order (1/μ) term give where the relativistic Fermi energy by using Eq. 1 is defined as It reasonably reveals that μ→E F when T→0 in Eq. 38. Substituting Eq. 38 into Eq. 34, the pressure of the electron gas in Eq. 34 is It explicitly tells us that the total pressure depends on temperature complicatedly at the given particle number N and volume V. The pressure of the Fermi electron gas in the white dwarf star should have something to do with the temperature as it is in Eq. 40.
After obtaining the pressure of the degenerate Fermi electron gas varying with temperature, then we can estimate the relation between mass and radius of the white dwarf star without considering the equations of the hydrostatic equilibrium. The relation between V and R is Using Eqs. 4 and 10, (11) (35), it gives (Huang, 1987) The equilibrium condition by considering the gravitational self-energy (Huang, 1987) is where In Eq. 38, δ is a parameter of pure number and G is the gravitational constant (Greiner et al., 1995). Substituting Eq. 38 into Eq. 37, then we improved Eq. 9 to be M R Further arrangement gives Then, we adopt the approximation of E F ∼ cp F in Eq. 46, and it gives Substituting Eq. 42 into Eq. 47 and arranging it, then we have Solving R in Eq. 48 and expanding the square root to the second order, the relation between radius, temperature, and mass is M 0 in Eq. 50 is the upper mass limit of the white dwarf star in the extremely relativistic case (Huang, 1987;Greiner et al., 1995). Equation 49 at T 0 K shows the same result as the case of p F >> m e c in some references (Huang, 1987;Greiner et al., 1995). According to Eq. 49 at T 0 K, the radius of the white dwarf star with the solar mass is about 2,700 km (Greiner et al., 1995). The more important thing is that the temperature term appears in the denominator of Eq. 49, which is much more reasonable than it is at T 0 K. It explicitly tells us that the relation between M and R depends on T, m n , and m e . Then, the mass-radius curves described in Eq. 49 at T 10 8 K and 8 × 10 8 K are respectively drawn in Figure 1. The upper mass limit M 0 is 1.44 M ⊙ at δ 0.65 and M ⊙ 1.99 × 1,030 kg (Greiner et al., 1995;Graham, 2000). The T 0 K curve is very close to the T 10 8 K curve. The temperature effect is explicit and the mass-radius curve at the high temperature is different from it at the low temperature. Most early research has derived the relation of the radius and mass of the white dwarf star (Huang, 1987;Greiner et al., 1995). However, the result shows the divergence of the radius when mass goes to zero for the lowmass white dwarf star (Huang, 1987;Greiner et al., 1995). As we know, the density of the white dwarf star is about (Huang, 1987;Greiner et al., 1995) ρ ≡ N V ≈ 10 7 g cm 3 ≈ 10 13 mole m 3 ≈ 6 × 10 36 electrons m 3 When we consider the ultra-relativistic condition p F >> m e c, it means 8π m e c 2 3 3h 3 5.2 × 10 35 electrons m 3 As long as the white dwarf star has the minimal density as shown in Eq. 51, this condition in Eq. 52 is automatically satisfied. In our common environment on the Earth, most things have density much lower than that of the white dwarf star in Eq. 51. When something disappears like the water evaporating, its volume will become larger and larger, and then, its density is closer and closer to zero. A white dwarf star with very large radius as its mass approaches zero is much unreasonable because it directly violates the criterion in Eq. 52 and our knowledge about the white dwarf star. The similar mass-radius relation trend has been shown in Figure 9 of the reference (Bvdard et al., 2017) where pure iron cores surrounded by helium and hydrogen layers are considered. Three demonstrations reveal arc curves whose maximum radii appearing at mass M between 0.40 and 0.55 M ⊙ , and then, the radii in these three cases are all decreased when the mass of the white dwarf star is below maximum (Bédard et al., 2017). Even for the C/O-core model, a core consisting of a uniform mixture of carbon and oxygen in equal proportions, the mass-radius relations of three different temperatures all show gradual decrease below M 1.10 M ⊙ (Bédard et al., 2017). Another similar mass-radius trend at the low-mass region is shown in Figure 2 of the reference (Nunes et al., 2021) where it clearly exhibits all five temperature cases, leading to the radial coordinate R→0 at the mass of the white dwarf star M→0. This result supports our above statements that the radius of the white dwarf star is close to zero when its mass approaches zero. The zeromass thing shall occupy no volume in space, not to be a sphere with a divergent radius. On the other hand, although our calculations do not have relation with the magnetic field, it is still worthy to mention that some recent research also show the radius of the highly magnetized white dwarf star approaching zero at their mass M close to zero (Das and Mukhopadhyay, 2012;Kundu and Mukhopadhyay, 2012;Das and Mukhopadhyay, 2013). This kind of the strongly magnetic white dwarf star has the strong magnetic field about 10 15 -10 17 G at the center (Das and Mukhopadhyay, 2012) and may explain the origin of overluminous peculiar type Ia supernova. The research of the strongly magnetic carbon-oxygen white dwarf star exhibits the mass-radius relation in Figure 2 ( Kundu and Mukhopadhyay, 2012) where the Chandrasekhar FIGURE 1 | The relation between the radius and mass of the white dwarf star calculated by Eq. 49 at the condition p F >> m e c without considering the equations of hydrostatic equilibrium. mass limit is exceeded and the generic mass limit is proposed to be 2.58 M ⊙ . The mass-radius relations in Figures 2 and 3 of the reference (Das and Mukhopadhyay, 2012) are shown for onelevel, two-level, and three-level systems. The curve of the onelevel system is single-valued, but one mass corresponds to three radii at most for the two-level system and four radii at most for the three-level system (Das and Mukhopadhyay, 2012). All curves also exhibit radii close to zero at mass M approaching zero. Usually speaking, the lower level system corresponds to the higher magnetic field at the center, and the larger Fermi energy has the higher magnetic field at the center (Das and Mukhopadhyay, 2012). For example, the central magnetic field is 6.6 × 10 13 G at E F 2m e c 2 for the one-level system, but it is as high as 8.82 × 10 17 G at E F 200 m e c 2 for the same system (Das and Mukhopadhyay, 2012). The higher-level mass-radius relation trend is much different from the non-magnetic white dwarf stars that we discuss in this manuscript. In addition, Figure 3 of the reference (Das and Mukhopadhyay, 2013) also show the mass-radius relations of the magnetic white dwarf stars for one-level systems where all the radii approach zero as their mass are close to zero. By the way, it also points out that the peculiar type Ia supernovae SN 2003fg, SN 2006gz, SN 2007if, and SN 2009dc may have super-Chandrasekhar mass limit of up to 2.4-2.8 M ⊙ (Das and Mukhopadhyay, 2012). Furthermore, the stability of such strongly magnetic white dwarf stars has also been discussed (Xhamel et al., 2013).
Another thing that we have to notice is the radius of the white dwarf star at its mass close to the upper mass limit. The mass-radius relation trends of most research and textbooks (Huang, 1987;Greiner et al., 1995;Carroll and Ostlie, 2006;Parsons et al., 2010;De Carvalho et al., 2014;Hermes et al., 2014;Franzon and Schramm, 2015;Bera and Bhattacharya, 2016;Boshkayev et al., 2016;Parsons et al., 2017;Sahu et al., 2017;Tremblay et al., 2017;Nunes et al., 2021) including ours all predict a zero radius when the mass is equal to the Chandrasekhar mass limit. However, it belongs to the black hole category, so we have to define a lower radius limit when the mass is close to M 0 . Because we consider a non-rotating and non-charged white dwarf star, this lower radius limit is the Schwarzschild radius, which is GM 0 /c 2 4.25 km for M 0 1.44 M ⊙ .
The other consideration is about m e c 2 >>E F for the nonrelativistic case, then the partition function, pressure, and number of the ideal Fermi electron gas become P electron gas ≈ 8π(2m e ) 3/2 (k B T) 5/2 3h 3 Γ 5 2 f 5/2 (z) and Combining Eq. 55 with Eq. 54, it gives the relation between P electron gas , T, V, and N P electron gas ≈ 2Nk B T 3V Here, we use the non-relativistic chemical potential for E F >> k B T (Honerkamp, 2002) μ Further rearrangement of Eq. 56 and ignoring T 2 terms in Eq. 56 and T 4 term in Eq. 57 give After obtaining the pressure of the ideal Fermi electron gas varying with temperature, then we can estimate the relation between mass and radius of the white dwarf star. The equilibrium condition in Eq. 43 becomes 8πm 3 e c 3 3h 3 and then it gives This reasonably shows that the radius disappears as M→0 verifying the previous viewpoint that zero mass represents no occupation or zero radius. However, the previous result shows the divergence of radius at M→0 (Carroll and Ostlie, 2006), and the mass-radius relation is where Z and A are the number of protons and nucleons, respectively. This expression is very problematic because the pressure is close to zero as R→∞ which can be seen directly from Eq. 59 that p is proportional to the inverse volume of the star or ∝ 1/R 3 . As we know, the pressure inside the white dwarf star is much higher than it is inside the Sun. The zero pressure inside the white dwarf star cannot satisfy the criteria in Eq. 52. When the pressure is close to zero everywhere, it cannot support a stable star anymore and the white dwarf star disappears at this zero-density and zero-pressure limit situation. Therefore, it is incorrect that the radius of the white dwarf star is divergent when its mass goes to zero. Besides, in Figure 6 of the reference (De Carvalho et al., 2014), the pressure in 12 C white dwarf star decreases from 10 22 to 10 16 erg/cm 3 when the density decreases from 5 × 10 6 to 2 × 10 2 g/cm 3 . Theoretically speaking, when its density goes to zero, the pressure must approach zero. If the pressure is not zero but the density is zero at some temperatures (De Carvalho et al., 2014), then the calculations shall have some problem. A white dwarf star with a very large radius when its mass goes to zero is much unreasonable because it directly violates the criterion in Eq. 52 and our knowledge about the white dwarf star. In addition, as mentioned before, another similar mass-radius relation trend at the low-mass region clearly exhibits in Figure 2 of the reference (Nunes et al., 2021) that all five temperature cases lead to the radial coordinate R→0 as the mass of the white dwarf star M→0. This result supports our prediction given in Eq. 61 where R→0 at M→0. In summary, both the relativistic and non-relativistic electron gases lead to zero radius at M approaching zero, which correct the divergent result in the old R ∝ M −1/3 .

THE MASS-RADIUS RELATION OBTAINED BY CONSIDERING THE EQUATIONS OF HYDROSTATIC EQUILIBRIUM
Next, we further obtain some relations from the equation of hydrostatic equilibrium for the stellar structure in the Tolman-Oppenheimer-Volkoff (TOV) form (Koester and Chanmugam, 1990;De Carvalho et al., 2014;Boshkayev et al., 2016;Carvalho et al., 2018). There are two equations in TOV considered here without the correction 1/c 2 terms (Boshkayev et al., 2016): and dm(r) dr 4πr 2 ρ(r), where p(r), m(r), and ρ(r) are the distributions of the pressure, mass, and mass density varying with the radial position in the star, respectively. This ignorance of the 1/c 2 terms makes us pay more attention to the dominant terms. Theoretically speaking, once the distribution of ρ(r) is known, the distributions of the pressure and mass can be obtained by substituting ρ(r) into Eqs. 63 and 64. The boundary conditions are m (0) 0, ρ(0) ρ c , m (R + ) ρ(R + ) 0, and dρ(r)/dr 0 at r 0 where ρ c is the central density.
Two boundary conditions of m (0) 0 at the center and p (R + ) 0 have been introduced in the paper (Koester and Chanmugam, 1990;Carvalho et al., 2018;Nunes et al., 2021). Then, defining a parameter where Using this parameter in Eq. 65, the high Fermi-energy pressure in Eq. 40 becomes x 2 (2x − 1) Finally, we can obtain the differential equation of ρ(r) by substituting Eq. 69 into Eq. 68. Similarly, the case for the low Fermi-energy pressure can be obtained by using Eq. 56. For μ≥k B T, where and the expression of μ in Eq. 57 is adopted here. Using Eqs. 64 and 68 with the boundary conditions ρ(0) ρ c , ρ(R + ) 0, and dρ/dr 0 at r 0, we can obtain the central density-mass, central density-radius, and mass-radius relations numerically. The numerical method is the fourth-order Runge-Kutta method (Nakamura, 1995). According to the definition of the parameter x in Eq. 65, the calculations are divided into three regions by considering the central mass density ρ c at different temperature. The temperature of the white dwarf star is considered homogeneously here. The highdensity region is for ρ c > ρ 0 , and we choose ρ c ≥ 5 × 10 10 kg/m 3 in our high-density calculations by using Eqs. 67 and 69. The low-density region is for ρ c < ρ, and we choose ρ c ≤ 10 9 kg/m 3 in the low-density calculations by using Eqs. 56 and 70. Between 10 9 kg/m 3 and 5 × 10 10 kg/m 3 is the middledensity region, where it can be approximated by connecting the high-and low-density regions directly. In Figure 2A, the relation between mass and the central density of the white dwarf star is given in the high-, middle-, and low-density regions at different temperature. In the high-density region, the mass is close to 1.4 M ⊙ after 10 13 kg/m 3 . Those curves are almost the same one from low temperature to 10 8 K in this region, so only two cases at 10 7 K and 10 8 K are shown. In the low-density region, those curves have tiny deviations until T 10 7 K, and especially, they are almost coincident at ρ c ≥ 2 × 10 7 kg/m 3 . At T 5 × 10 7 K, the curve is explicitly changed, and the starting point is at 5 × 10 6 kg/m 3 . It means that the white dwarf star at this temperature has a central density higher than 5 × 10 6 kg/m 3 . As temperature increases, the starting point of the central density also increases. It is about 10 7 kg/m 3 at 10 8 K, 4 × 10 7 kg/m 3 at 2 × 10 8 K, and 7 × 10 7 kg/m 3 at 3 × 10 8 K. In the middle-density region, although the central density is only from 10 9 kg/m 3 to 5 × 10 10 kg/m 3 , the range of mass covers a large interval from 0.4 M ⊙ to 1.1 M ⊙ . It means that a large part of the white dwarf stars are in this region.
In Figure 2B, the relation between the central density and mass is given in the high-, middle-, and low-density regions at different temperature. Both axes are shown in log scale. The central density is from 10 5 kg/m 3 to 10 14 kg/m 3 , and the radius is from 5 × 10 5 to 5 × 10 7 . At 10 6 K, it shows a logarithm relation between the central density and the radius of the white dwarf star. The higher central density is, the smaller the radius is. However, the trend is broken at 10 7 K and above. At 10 7 K, there is a turning point around the central density of 10 6 kg/m 3 , which means that the maximal radius of the white dwarf star at this temperature is less than 3 × 10 7 m or 3 × 10 4 km. This turning point also means that the radius of the white dwarf star cannot increase infinitely as mentioned in Eq. 62. When the mass goes to zero, the radius also approaches zero. If we extend the curve to zeros radius, then it exhibits that the central density is always above 10 5 kg/m 3 at 10 7 K. The turning point increases in the central density as temperature increases, but the maximal radius decreases at the same time. The turning point is roughly 10 7 kg/m 3 at 5 × 10 7 K and 10 8 kg/m 3 at 2 × 10 8 K. In conclusion, the same radius of the white dwarf star would correspond to different central density, and one is in the high-density region and the other is in the lowdensity region.
In Figure 2C, the mass-radius relation for the white dwarf star is shown. As mentioned in Figure 2A, the mass is close to 1.4 M ⊙ after 10 13 kg/m 3 in the high-density region. Those curves are almost the same one from low temperature to 10 8 K in the highdensity region, and only two cases at 10 7 K and 10 8 K are shown in this region. In the low-density region, those curves have tiny deviations until T 10 7 K, and it is almost coincident at ρ c ≥ 2 × 10 7 kg/m 3 . In the middle-density region, the range of mass covers a large interval from 0.4 to 1.1 M ⊙ where the central density is only from 10 9 kg/m 3 to 5 × 10 10 kg/m 3 . It also implies that a large part of the white dwarf stars that we found astronomically belong to the middle-density region. Those results are the same as Figure 2A. Especially, at T 5 × 10 7 K and above, there exist some parts where the mass at the same radius is larger than the curve at 10 7 K and below. Recently, the Sloan Digital Sky Survey Release 4 shows a lot of observations having larger mass compared to the relativistic EOS at T 0 K when the radius is larger than 8 × 10 3 km (De Carvalho et al., 2014). By using our calculations in Figure 2C, this phenomenon can be explained because the higher temperature results in these white dwarf stars with larger mass appearing at the same radius in the low-density region. Usually, the temperature of these white dwarf stars is higher than 10 7 K. Those parts are denoted by the dotted elliptic curve in Figure 2C. This explanation can also extend to the middle region in Figure 2C.
Compare Figure 2C with Figure 1, both mass-radius relation trend is similar for T > 0. We can find that all cases show the maximum radii appearing between 0 and 1.4 M ⊙ . In addition, as mentioned previously, the similar mass-radius relation trend has been shown in Figure 9 of the reference (Bédard et al., 2017) where pure iron cores surrounded by helium and hydrogen layers are considered. Three demonstrations reveal arc curves whose maximum radii appearing at mass M between 0.40 and 0.55 M ⊙ , and then, the radii in these three cases are all decreased when the mass of the white dwarf star is below maximum (Bédard et al., 2017). Even for the C/O-core model, a core consisting of a uniform mixture of carbon and oxygen in equal proportions, the mass-radius relations of three different temperatures all show gradual decrease below M 1.10 M ⊙ (Bédard et al., 2017).

CONCLUSION
In summary, the mass-radius relation of the white dwarf star derived according to statistical mechanics shows that the temperature effect has to be considered at high temperature above 10 7 K. After all, the ideally degenerate Fermi electron gas is described at T 0 K, and the temperature effect would show something difference. The other correction is due to the electron-electron interaction considered at T 0. The calculation considers the relativistic electrons, and the result shows that this effect causes the pressure 2/137 time less than the original value. It means that the many-particle effect appears and causes about 1.5% deviation in pressure. When the temperature effect is considered, the pressure is calculated by statistical mechanics. According to the Fermi energy, two cases are calculated. One is E F >> m e c 2 for the relativistic case, and the other is E F << m e c 2 for the non-relativistic case. Because of the temperature effect, the chemical potential is also temperature-dependent and different expression in these two cases. From the deductions, the pressure produced by the Fermi electron gas depends on temperature complicatedly at the given particle number N and volume V. Traditional formula gives the problematic relation R ∝ M −1/3 because it leads to R→∞ and p→0 as M→0. Therefore, we correct this relation and obtain two much reasonable relations in the relativistic and nonrelativistic regions. The mass-radius relation of the white dwarf star without considering the equations of hydrostatic equilibrium is obtained for the relativistic and non-relativistic cases in which the radius is shown to be zero as the mass of the white dwarf star disappears.
Then, further considering the equations of hydrostatic equilibrium, the central density-mass, the central density-radius, and the mass-radius relations are obtained. The central density is divided into the high-, middle-, and low-density regions where the results can be coincident with the SDSS observations. In the high-density region, three relations are almost unchanged until 10 8 K. The temperature effect mainly affects the low-density region at temperature above 10 7 K. Especially, the mass-radius curves show some parts having larger mass at the same radius when temperature is higher. It gives a way to explain the SDSS observations at the radius more than 8 × 10 3 km that the mass of the white dwarf star is often larger than the prediction by the relativistic EOS at zero temperature. Those white dwarf stars of larger mass just correspond to the low-and middle-density regions. It means that we should consider the temperature effect to get better calculations that can reasonably explain the astronomical observations. Although the temperature is maximum at the center and minimum at the surface in reality, we still adopt the uniform temperature approximation as most research did. Our results imply that the uniform temperature approximation may remain valuable. Furthermore, the temperature can change with time because the white dwarf stars cool down by thermal emission. However, the cooling time is much longer than the astronomical observations, the temperature of the white dwarf star is kept constant at each calculation in our research. It is a good enough way to hold the initial temperature throughout each calculation if we want to compare the calculation results with the astronomical observations in several decades.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.

AUTHOR CONTRIBUTIONS
The single author first reads a lot of references to propose this research idea. Then, he carefully derives all the mathematical equations in this manuscript and checks them several times. He also writes down the programs to calculate some examples to demonstrate the concepts revealing in this manuscript. After obtaining some results, he starts to write this paper and finishes it. Finally, he prepares this paper and submits it.