×
近期发现有不法分子冒充我刊与作者联系,借此进行欺诈等不法行为,请广大作者加以鉴别,如遇诈骗行为,请第一时间与我刊编辑部联系确认(《中国物理C》(英文)编辑部电话:010-88235947,010-88236950),并作报警处理。
本刊再次郑重声明:
(1)本刊官方网址为cpc.ihep.ac.cn和https://iopscience.iop.org/journal/1674-1137
(2)本刊采编系统作者中心是投稿的唯一路径,该系统为ScholarOne远程稿件采编系统,仅在本刊投稿网网址(https://mc03.manuscriptcentral.com/cpc)设有登录入口。本刊不接受其他方式的投稿,如打印稿投稿、E-mail信箱投稿等,若以此种方式接收投稿均为假冒。
(3)所有投稿均需经过严格的同行评议、编辑加工后方可发表,本刊不存在所谓的“编辑部内部征稿”。如果有人以“编辑部内部人员”名义帮助作者发稿,并收取发表费用,均为假冒。
                  
《中国物理C》(英文)编辑部
2024年10月30日

New insights into the limit of the magnetic monopole flux and the heating source in white dwarfs

Figures(2) / Tables(11)

Get Citation
Jing-Jing Liu, Dong-Mei Liu and Liang-Huan Hao. New insights on the limit of the magnetic monopole flux and the heating source in white dwarfs[J]. Chinese Physics C. doi: 10.1088/1674-1137/acdc8b
Jing-Jing Liu, Dong-Mei Liu and Liang-Huan Hao. New insights on the limit of the magnetic monopole flux and the heating source in white dwarfs[J]. Chinese Physics C.  doi: 10.1088/1674-1137/acdc8b shu
Milestone
Received: 2023-04-09
Article Metric

Article Views(975)
PDF Downloads(31)
Cited by(0)
Policy on re-use
To reuse of subscription content published by CPC, the users need to request permission from CPC, unless the content was published under an Open Access license which automatically permits that type of reuse.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

New insights into the limit of the magnetic monopole flux and the heating source in white dwarfs

  • College of Science, Hainan Tropical Ocean University, Sanya 572022, China

Abstract: Based on the magnetic monopole (MM) catalytic nuclear decay (Rubakov-Callan (RC) effect), we propose five new models to discuss the limit of the MM flux and the heating energy resources of white dwarfs (WDs) based on observations of 13 red giant branch (RGB) stars. We find that the number of MMs captured can reach a maximum value of $ 9.1223\times10^{24} $ when $ m=10^{17} $ GeV, $n_{\rm B}=5.99\times10^{31}~\rm{cm^{-3}}$, $\phi=7.59\times 10^{-26} \rm{cm^{-2}s^{-1}sr^{-1}}$. The good agreement of our calculated luminosities for WDs with observation provides support for our model based on the RC effect by MMs. We obtain a new limit of the MM flux of $\xi=\phi\langle\sigma_m v_{\rm{T}}\rangle_{-28}\leq $$ 9.0935\times10^{-13}\rm{cm^{-2}s^{-1}sr^{-1}}$, and $ \xi\leq 4.9950\times10^{-13}\rm{cm^{-2}s^{-1}sr^{-1}} $ at $n_{\rm B}=5.99\times10^{31}\rm{cm^{-3}}$ when $m=10^{15}\rm{GeV}$, $ \beta=9.4868\times10^{-3} $, and $ m=10^{17}\rm{GeV}, $$ ~\beta=10^{-3} $, respectively. Our results show that the RC effect could cause heating that prevents white dwarfs from cooling down into a stellar graveyard. Our results will also provide a new idea for further research on the upper limit of MM flow (note: $n_{\rm B}, \sigma_m, m, \phi, \xi$ are the baryon number density, reaction cross section, mass, MM flux, and the new limit of the MM flux, respectively, and $ \beta=v_T/c $ is the ratio of the speed of MMs to that of light).

    HTML

    I.   INTRODUCTION
    • White dwarfs (hereafter WDs) are usually made up of C+O. However, it is also possible for their cores to be hot enough to burn carbon but not hot enough to burn Ne, forming a WD with a core of O+Ne+Mg. At the later stage of WDs, the star ejects large quantities of matter. After great mass loss, if the remaining core mass is less than 1.44 solar masses, the star may evolve into a WD.

      WDs form at very high temperatures. Because they have no source of energy, they will therefore gradually give off heat and cool down, whereupon its radiation will decrease over time from its initial high color temperature to red. This surface temperature is defined in astronomy as the effective temperature $ T_{\rm{eff}} $ as per Stefan's law, so that

      $ \begin{array}{*{20}{l}} L_{\rm{rad}}=4\pi R^2\sigma T^4_{\rm{eff}}, \end{array} $

      (1)

      where R the radius of the star, $ L_{\rm{rad}} $ is the radiation luminosity, and $ \sigma=5.6704\times10^{-5}\rm{erg s^{-1}cm^{-2}K^{-4}} $ is the radiation constant from Stefan's law. $ T_{\rm{eff}} $ is a measure of the energy flux at the surface and not a real temperature, but it nevertheless constitutes a useful measure of the atmospheric temperature of the star.

      As is well known, the effective temperature of WDs is mostly in the range 5500–40000 K, while a few are outside this range, and the internal temperature of WDs is on the order of $ \sim10^6-10^7 $K, with a total thermal energy less than $ 10^{47} $ergs. Mestel. [1] discussed the energy sources of WDs. Avakian. [2] also studied the configurations of hot WDs with nuclear sources of energy. Bildsten & Hall. [3] discussed the sources of white dwarfs, suggesting that there are small amounts of $ ^{22} $Ne in some WDs that may constitute an extra source of heat in carbon-oxygen WDs. Single-particle $ ^{22} $Ne sedimentation may be considered a possible heat source [4, 5]. However, some work suggests that $ ^{22} $Ne must separate into clusters, enhancing diffusion, in order for sedimentation to provide heating on the observed timescale. Recently, the sources of ultra-high-energy photons for WD pulsars have been discussed by Lobato et al. [6]. Cheng et al. [7] discuss the cooling anomaly of high-mass WDs, pointing out that $ ^{22} $Ne settling in C/O-core WDs could account for this extra cooling delay. Caplan et al. [8] studied this topic using molecular dynamics methods and phase diagrams, from which they ruled out the isotope $ ^{22} $Ne as a possible cause of the extra heating. Therefore, the problem of additional heat sources for WDs remains a challenging topic.

      In this paper, we selected 13 red giant branch (RGB) stars to present a new model of the number of magnetic monopoles (hereafter MMs) captured, and seek to solve the energy source problem for WDs based on MM catalytic nuclear decay (the Rubakov-Callan (RC) effect) [9, 10]. MMs are hypothetical magnetic particles with a single north or south magnetic pole, which have been proposed in string theory. Research on MMs has long been a hot topic among physicists and astronomers. Some papers have discussed the issues of MMs, (e.g., Callan [9], Detrixhe et al. [11], Frank et al. [12], Fujii & Pierre [13], Kain [14], Rajantie [15]). Recently, Mavromatos & Mitsou. [16] discussed the developments in both theory and experimental searches for MMs in past, current, and future colliders and in the cosmos. We also are interested in the problem of MMs and other related issues ((e.g., Liu. [17, 18], Liu & Gu [19], Liu et al. [20], Liu [2126], Peng et al. [27, 28]).

      The arrangement of this paper is as follows. In the next section, we discuss the number of possible magnetic monopoles in space and the luminosity due to the RC effect by MMs. In Section III, we describe our models and the luminosity function due to magnetic monopole catalytic nuclear decay. In Section IV, some results and discussions are presented. Finally, our conclusions are summarized in Section V.

    II.   THE NUMBER OF MMs CAPTURED IN SPACE AND THE LUMINOSITY DUE TO THE RC EFFECT BY MMs
    • According to some research, the interaction of MMs with neutral hydrogen atoms is very weak. Therefore, during the process of formation of celestial bodies, very few MMs are captured in the collapse of a neutral hydrogen cloud and collect in the core of a star or planet. MMs typically may be contained within stars and planets, and they are mostly captured from space during their lifetime after formation. One type of interaction that MMs undergo in stars and planets might be the RC effect, through which MMs may catalyze nucleon decay, as expressed by $p+M\rightarrow e^+\pi^0+ M+ {\rm debris}$ (x%) and $p+M\rightarrow e^+ \mu^{\pm} + M+ {\rm debris}$ (y%). The ratio of the cross sections of the above reactions $ x/y $ may be a few percent, which is on the order of $ \sim10^{-4} $ [10]. Bernreuther and Craigie [29] discussed the cross section of monopole-induced proton decay in SU(5) for the above reactions, showing that the ratio of cross sections of the above reactions $ x/y $ may be $ (2.5\sim2.8)\times10^{-4}/(0.1\sim0.3) $.

      The number (flow) of MMs captured in space at the surface of stars and planets (including earth) is estimated as follows [30]:

      $ N_{m}=4\pi^2 R^2\eta\phi t\left[1+\left(\frac{v_{\rm{esc}}}{v_m}\right)^2\right], $

      (2)

      where ϕ is the flux of MMs intercepted in space and $ v_{\rm{esc}}=(2GM/R)^{1/2} $, t are the escape velocity from the star and cooling age of the star, respectively. $ v_m $ represents the velocity of MMs in space, For a planet, $ \dfrac{v_{\rm{esc}}}{v_m}\ll1 $, and $ \dfrac{v_{\rm{esc}}}{v_m}\approx1 $ for ordinary stars, but $ \dfrac{v_{\rm{esc}}}{v_m}\gg1 $ for compact objects such as WDs, neutron stars, galactic nuclei, and quasars. Therefore, the ratio between the number of MMs captured and the number of nuclei in the MMs accumulation area is [30]

      $ \frac{N_m}{N_B}=5.10\times10^{-34}\eta R^2_\ast\left(\frac{t_9}{M_\ast}\right)\left(\frac{\phi}{\phi_0}\right)\left[1+4.256v^{-2}_{-3}\left(\frac{M_\ast}{R_\ast}\right)\right], $

      (3)

      where $ R_\ast=R/R_{\odot} $, $ M_\ast=M/M_{\odot} $, $ v_{-3}=v_m/10^{-3}c $, $ t_{9}=t/10^9 $Yr, and c is the speed of light. $ \phi_0\approx 10^{-12} \rm{cm^{-2}s^{-1}Sr^{-1}} $ [31]. η is the probability of the capture of MMs by stars, which depends on the ratio of the penetration distance $ l_{\rm{pd}} $ of an MM in a star to the star's radius. In general, we have $ l_{\rm{pd}}\approx1.2\times10^{30}v_{-3}n_e^{-1}T_e^{1/2} $ for plasma [30]. For example, with $n_e\sim10^{22} {\rm cm^{-3}}, T_e\sim10^6$K for the sun, we have $ l_{\rm{pd}}\sim10^{11}, \eta\sim0.7 $; however, for WDs and neutron stars, $ n_e\geq10^{30}, n_e\geq10^{35}\rm{cm^{-3}} $, respectively, and $ \eta\sim1 $. For quasars and active galactic nuclei, $ M_\ast\sim10^8, T_e\sim10^5 $K, $ \eta\sim1 $.

      The velocity of MMs is also determined as a function of monopole mass by $ \beta_{\rm{T}} $ as follows [32].

      $ \begin{array}{*{20}{l}} v_m=c\beta_{\rm{T}}=\left\lbrace \begin{array}{ll}\; 3\times10^{-3}c\left(\dfrac{10^{16}\rm{GeV}}{m}\right)^{0.5}\; \; ( m<10^{17}\rm{GeV})\\ \; 10^{-3}c \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; (\rm{otherwise}). \end{array} \right. \end{array} $

      (4)

      According to Eq. (3), for WDs we have $ \eta\sim1 $; thus, the total number of magnetic monopoles trapped in space after the formation of stars (or planets) is estimated to be

      $ N_{m}=7.18\times10^{11}n_{\rm{B}}R^5_\ast\phi\left(\frac{t_9}{M_\ast}\right)\left[1+4.256v^{-2}_{-3}\left(\frac{M_\ast}{R_\ast}\right)\right]. $

      (5)

      For astronomy, the most important property of a magnetic monopole is that it can trigger the RC effect, as independently proposed by Rubakov and Callen [9, 10]. The reaction cross section is about $ \sigma_m\approx10^{-25}\sim10^{-26}\rm{cm}^2 $, almost reaching the Thomson cross section ($6.665\times 10^{-25}\rm{cm}^2$). The luminosity of various types of celestial bodiesdue to the RC effect (i.e., RC luminosity) can be estimated as follows. In the core area, where the magnetic monopole is concentrated, the nuclear decay reaction is catalyzed by the magnetic monopoles and the total luminosity produced is [30]

      $ L_m\approx\frac{4\pi}{3}r_c^3n_mn_{\rm{B}}\langle\sigma_m v_{\rm{T}}\rangle m_{\rm{B}}c^2=N_m n_{\rm B}\langle\sigma_m v_{\rm{T}}\rangle m_{\rm{B}}c^2, $

      (6)

      where $ r_c $, and $n_m, n_{\rm B}$ are the radius of the stellar central region and the number densities of MMs and nucleons, respectively.

      In Eq. (6), $ v_{\rm{T}} $ is the thermal movement speed of the nucleus relative to the magnetic monopole. We will ignore the thermal velocity of the magnetic monopole due to its great mass. We thus only consider the contributions from the thermal velocity of the nucleus. According to $1/2mv_{\rm T}^2=3/2kT$, we have $v_{\rm{T}}=\sqrt{3kT/m_B} \approx 1.5745\times 10^{7}T_6^{1/2}$cm/s, where T is the temperature, $ T_6=T/10^6 $K, $ k=1.38\times10^{-16}\rm{erg/s} $ is the Boltzmann constant, and $m_{\rm B}\approx1.67\times10^{-24}$g is the nucleon mass. As the central temperature of WDs is about $ \sim10^6 $K, we have $ v_T\sim10^{-3}c $.

      As a general rule, the reaction cross section of the RC effect $ \sigma_m $ is $10^{-26}\sim10^{-24}\;\rm{cm}^2$. Ma & Tang [33] gave a value $ \sigma_m\approx4.28676\times10^{-24}\rm{cm}^2 $ for the cross section of the proton with different channels using SU(5) grand unification theory. In the RC process, MMs induced nucleon decay, followed by nucleon decay into $ \pi^0 $ mesons, $ \mu^\pm $ leptons, and positrons $ e^+ $, and $ \mu^\pm $ and $ \pi^0 $ again decay into photons and electron-positron pairs $ e^\pm $. The positrons then undergo annihilation with electrons to form photons. The net effect is that the rest mass energy of nucleons ($ m_{\rm{B}}c^2 $) is entirely converted to radiant energy with 100% efficiency ($ 1m_{\rm{B}}c^2\approx1\rm{GeV}\approx1.6\times10^{-3}\rm{ergs} $).

    III.   THE MAGNETIC MONOPOLE MODEL AND RC LUMINOSITY INSIDE WHITE DWARFS

      A.   The mass and radius for the red giant branch (RGB) phase

    • It is well known that WDs evolve from red giants with mass less than $ 8M_{\odot} $. The MM content of the trapping accumulation set is mainly from the red giant branch (RGB) phase. The $ M_{\ast} $ and $ R_{\ast} $ for RGB stars is given by [34, 35]

      $ M^{\rm{r}}_{\ast}=\frac{M_{\rm{RG}}}{M_\odot}=\left(\frac{\nu_{\rm{max}}}{\nu_{\rm{max},\odot}}\right)^3\left(\frac{\Delta\nu}{f_{\Delta\nu}\Delta\nu_{\odot}}\right)^{-4}\left(\frac{T_{\rm{eff}}}{T_{\rm{eff},\odot}}\right)^{1.5}, $

      (7)

      $ R^{\rm{r}}_{\ast}=\frac{R_{\rm{RG}}}{R_\odot}=\left(\frac{\nu_{\rm{max}}}{\nu_{\rm{max},\odot}}\right)\left(\frac{\Delta\nu}{f_{\Delta\nu}\Delta\nu_{\odot}}\right)^{-2}\left(\frac{T_{\rm{eff}}}{T_{\rm{eff},\odot}}\right)^{0.5}, $

      (8)

      where $ \nu_{\rm{max},\odot}=3090\mu $Hz, $ \Delta\nu_{\odot}=135.1\mu $Hz, $ T_{\rm{eff},\odot}=5777 $K, and $ \Delta\nu\simeq \alpha(\nu_{\rm{max}/\mu{\rm{Hz}}})^{\beta} $ (here $\alpha=0.268\mu$Hz, $ \beta=0.758 $)[36]. $ f_{\Delta\nu} $ is a typical asteroseismic correction factor, which is between 0.98 and 1.02 [37]. The predicted $ \nu_{\rm{max}} $ values are calculated using the scaling relation [38]

      $ \frac{\nu_{\rm{max}}}{\nu_{\rm{max}\odot}}=\left(\frac{g}{g_{\odot}}\right)\left(\frac{T_{\rm{eff}}}{T_{\rm{eff},\odot}}\right)^{-0.5}. $

      (9)

      where g is the surface gravity of the star and $ g_\odot=27487\;\rm{cm/s^2} $.

    • B.   The magnetic monopole RC effect model in white dwarfs

    • According to Eqs. (5)–(9), the number of magnetic monopoles captured from space and the total luminosity due to the RC effect by the MMs are given by

      $ N_{m}=7.18\times10^{11}n_{\rm{B}}(R^{r}_\ast)^5\phi\left(\frac{t_9}{M^{r}_\ast}\right)\left[1+4.256v^{-2}_{-3}\left(\frac{M^{r}_\ast}{R^{r}_\ast}\right)\right]. $

      (10)

      $ \begin{aligned}[b] L_m\approx&\frac{4\pi}{3}r_c^3n_mn_{\rm{B}}\langle\sigma_m v_{\rm{T}}\rangle m_{\rm{B}}c^2=N_mn_{\rm B}\langle\sigma_m v_{\rm{T}}\rangle m_{\rm{B}}c^2,\\ =&1.15\times10^{9}n^2_{\rm{B}}(R^{r}_\ast)^5\phi t_9(M^{r}_\ast)^{-1} \\ &\times\langle\sigma_m v_{\rm{T}}\rangle\left[1+4.256v^{-2}_{-3}\left(\frac{M^{r}_\ast}{R^{r}_\ast}\right)\right]. \end{aligned} $

      (11)

      Defining $ \xi=\phi\langle\sigma_m v_{\rm{T}}\rangle_{-28}=\phi\langle\sigma_m v_{\rm{T}}\rangle/10^{-28} $, Eq. (11) may be rewritten

      $ \begin{aligned}[b] L_m\approx&\frac{4\pi}{3}r_c^3n_mn_{\rm{B}}\langle\sigma_m v_{\rm{T}}\rangle m_{\rm{B}}c^2=N_mn_B\langle\sigma_m v_{\rm{T}}\rangle m_{\rm{B}}c^2,\\ =&1.15\times10^{37}n^2_{\rm{B}}(R^{\rm{r}}_\ast)^5t_9(M^{\rm{r}}_\ast)^{-1}\xi\left[1+4.256v^{-2}_{-3}\left(\frac{M^{\rm{r}}_\ast}{R^{\rm{r}}_\ast}\right)\right].\\ \end{aligned} $

      (12)

      From Eqs. (1) and (12), and $ L_m=L_{\rm{rad}} $, we obtain

      $ \xi=5.29\times10^{-15}\sigma M^{r}_\ast t_9^{-1} n^{-2}_{\rm{B}}(R^{r}_\ast)^{-3}T^4_{\rm{eff}}\left[1+4.256v^{-2}_{-3}\left(\frac{M^{r}_\ast}{R^{r}_\ast}\right)\right]^{-1}. $

      (13)
    IV.   RESULTS AND DISCUSSIONS
    • The study of MMs has held considerable interest since MMs were found to be a generic feature of grand unified gauge theories in the physical fields. The theoretical predictions of monopole abundance are problematic in the standard cosmology, as far too many monopoles would have survived annihilation for the universe to have reached its present state. For example, the galactic field that yields the Parker bound is $ \phi(\sigma v)_{-28}\leq 10^{-16} \rm{cm^{-2}s^{-1}sr^{-1}} $ [39]. Due to MM RC decay, another limit of the flux may be $ \phi(\sigma v)_{-28}\leq10^{-21}\rm{cm^{-2}s^{-1}sr^{-1}} $[31]. A better-understood limit in WDs may be $ \phi(\sigma v)_{-28}\leq10^{-18} \rm{cm^{-2}s^{-1}sr^{-1}} $ [32]. Then the bound has been stated as $ \phi(\sigma v)_{-28}\leq10^{-28} \;\rm{cm^{-2}s^{-1}sr^{-1}} $ by Freese and Krasteva [40]. In this paper, we study MMs and their numbers in space and discuss our MM model and the luminosity due to RC effect. We select the following typical parameters: $ m=10^{15}, 10^{17} $GeV, $ \sigma_m=10^{-24}\rm{cm^2} $, $\phi=5.59\times10^{-28}, 10^{-26}, 7.59\times 10^{-26}\rm{cm^{-2}s^{-1}sr^{-1}}$, and $n_{\rm B}=5.99\times10^{31}, 1.89\times 10^{32} \rm{cm^{-3}}$, $ t_9=1, 10 $.

      Everyone knows that WDs evolve from red giants with mass less than 8 M$ _\odot $, we focus on 13 typical RGB stars listed in Tables 1 and 2 [41] and discuss their mass, radius, and cooling age. Depending on an asteroseismic correction factor, we propose five models (I$-$V) to examine the problem of the energy source of WDs.

      starTIC IDGaia IDGania magDistance (pc)Galactic substructureSpectral source
      1TIC20897763Gaia 23656494710338280969.41484457.879 ± 9.435Gaia Enceladus SausageAPOGEE
      2TIC341816936Gaia 142177604633572300811.6231547.85 ± 45.63Gaia Enceladus SausageAPOGEE
      3TIC393961551Gaia 15063876279179368969.57166500.533 ± 6.312Gaia Enceladus SausageAPOGEE
      4TIC453888381Gaia 523025673034745715210.9714788.614 ± 15.999HaloGALAH
      5TIC279510617Gaia 548055045064301721610.7551933.263 ± 22.520HaloGALAH
      6TIC300938910Gaia 527067501829784422410.5629607.156 ± 7.7075HaloGALAH
      7TIC198204598Gaia 162989868534727385610.9455952.885 ± 37.512Halo···
      8TIC1008989Gaia 37896392809526103049.72882370.56 ± 5.85Gaia Enceladus Sausage···
      9TIC91556382Gaia 506500965033314739210.0855870.289 ± 34.565Halo···
      10TIC159509702Gaia 170919509028171827212.15421595.67 ± 55.605Halo···
      11TIC25079002Gaia 46693160657002229769.91465716.804 ± 15.5995DiskAPOGEE
      12TIC177242602Gaia 526229539536721228810.1451532.831 ± 14.294DiskAPOGEE
      13TIC9113677Gaia 324548565060765158410.1791491.438 ± 10.040Thick Disk···

      Table 1.  Information of the 13 red giant stars selected [41].

      star$ \nu_{\rm{max}} $$ \Delta\nu $ /μHzTeff /K[Fe/H][α/Fe]
      161.31383 ± 1.217686.81739 ± 0.249904988 ± 127−1.274 ± 0.0190.219 ± 0.021
      236.34 ± 0.764.2969 ± 0.07155068 ± 100−1.873 ± 0.1070.248 ± 0.023
      361.34 ± 1.756.68 ± 0.415121 ± 105−1.0751 ± 0.01230.156 ± 0.014
      450.37 ± 1.595.953 ± 0.0494741 ± 100−0.728 ± 0.070.32 ± 0.021
      528.57 ± 0.163.566 ± 0.0154450 ± 100−0.49 ± 0.050.281 ± 0.017
      6106.30 ± 0.929.464 ± 0.1314908 ± 100−0.792 ± 0.050.2566 ± 0.0165
      745.86 ± 0.315.132 ± 0.0324979 ± 100······
      8104.33059 ± 1.466189.80317 ± 0.153364893 ± 100······
      931.38665 ± 1.030974.189 ± 0.1865192 ± 100······
      1045.37 ± 0.535.090 ± 0.0274724 ± 100······
      1145.238 ± 0.624.967 ± 0.1214797 ± 830.1636 ± 0.006−0.010 ± 0.006
      1254.66 ± 0.335.663 ± 0.0314603 ± 100−0.1176 ± 0.0060.125 ± 0.007
      1387.22437 ± 0.870428.365 ± 0.1984764 ± 100······

      Table 2.  Information of the 13 RGB stars [41].

      In this paper, we select the asteroseismic correction factors $ f_{\bigtriangleup\nu}= $ 0.98, 0.99, 1.00, 1.01, 1.02 for our typical models under study, which correspond to models (I $-$ V). Based on Eqs. (8) and (9), we can calculate the mass and the radii of RGB stars $ M^{r}_\ast $ (I $-$ V) and $ R^{r}_\ast $ (I $-$ V) given in Table 3. We determine the stellar (RGB) ages using the packages BASTA [42], isochrones [43], isoclassify [44], PARAM [45, 46], and scaling-giants [47]. These ages correspond to $ t_9 $ (I $ \sim $ V) in Table 4.

      star$ M_{\ast} $(I)$ M_{\ast} $(II)$ M_{\ast} $(III)$ M_{\ast} $(IV)$ M_{\ast} $(V)$ R_{\ast} $(I)$ R_{\ast} $(II)$ R_{\ast} $(III)$ R_{\ast} $(IV)$ R_{\ast} $(V)
      10.89160.92860.96671.00591.04646.95417.09677.24087.38637.5333
      21.20471.25461.30611.35911.413810.45810.67210.88911.10811.329
      31.00751.04931.09231.13671.18247.34217.49277.64497.79857.9537
      40.78790.82060.85420.88890.92467.30457.45437.60567.75857.9129
      51.01541.05751.10091.14561.191611.18611.41611.64711.88212.118
      61.22111.27171.32391.37761.43306.20576.33306.46166.59146.7226
      71.15871.20671.25621.30721.35989.17039.35849.54849.74049.9342
      80.98961.04041.08301.12701.17235.66995.78625.90376.02246.1422
      90.89100.92800.96601.00521.04579.61939.816610.01610.21710.421
      100.92210.96030.99971.04031.08218.54488.72018.89729.07609.2566
      111.19861.24831.29951.35231.40669.47889.67329.869610.06810.268
      121.17621.22491.27521.32701.38038.63088.80788.98679.16739.3497
      131.05871.10261.14781.19451.24256.42506.55686.68996.82446.9602

      Table 3.  The mass and radii of the 13 RGB stars selected [41]. The asteroseismic correction factor $ f_{\bigtriangleup\nu}= $ 0.98, 0.99, 1.00, 1.01, 1.02 defines models I $-$ V, respectively.

      starscaling-giants Age (Gyr)isochrones Age (Gyr)isoclassify Age (Gyr)PARAM Age (Gyr)BASTA Age (Gyr)
      $ t_9 $(I)$ t_9 $(II)$ t_9 $(III)$ t_9 $(IV)$ t_9 $(V)
      15.8 $ \pm $ 3.08.77 $ \pm $ 2.95.68$ ^{+0.79}_{-0.83} $9.29$ ^{+2.82}_{-3.23} $9.0$ ^{+2.7}_{-2.5} $
      22.9 $ \pm $ 1.89.16 $ \pm $ 2.757.44$ ^{+1.36}_{-1.34} $6.52$ ^{+4.46}_{-3.80} $7.9$ ^{+2.8}_{-1.9} $
      34.7 $ \pm $ 6.05.93 $ \pm $ 3.055.63$ ^{+0.84}_{-0.69} $5.72$ ^{+4.61}_{-3.14} $7.5$ ^{+3.3}_{-2.2} $
      410.5 $ \pm $ 3.59.72 $ \pm $ 2.5012.78$ ^{+0.78}_{-1.59} $11.66$ ^{+1.50}_{-2.42} $14.8$ \pm $2.8
      56.4 $ \pm $1.17.99 $ \pm $ 3.3710.23$ ^{+0.95}_{-2.34} $7.68$ ^{+4.55}_{-4.45} $7.6$ ^{+0.9}_{-0.8} $
      62.7 $ \pm $ 0.83.693 $ \pm $ 0.7023.34$ ^{+0.40}_{-0.33} $2.11$ ^{+0.61}_{-0.49} $2.9$ ^{+0.6}_{-0.4} $
      73.3 $ \pm $ 0.88.749 $ \pm $ 2.1611.54$ ^{+0.66}_{-2.83} $1.98$ ^{+7.13}_{-0.77} $5.8$ ^{+1.6}_{-1.2} $
      85.8 $ \pm $ 2.06.843 $ \pm $ 3.336.84$ ^{+1.50}_{-0.90} $7.27$ ^{+2.95}_{-2.09} $10.6$ ^{+3.5}_{-2.9} $
      97.1 $ \pm $ 5.06.841 $ \pm $ 3.596.43$ ^{+1.89}_{-0.96} $5.27$ ^{+2.84}_{-1.77} $8.9$ ^{+3.8}_{-2.9} $
      104.8 $ \pm $ 1.29.344 $ \pm $ 2.3319.98$ ^{+1.47}_{-1.18} $7.92$ ^{+3.47}_{-2.49} $10.0$ ^{+3.2}_{-2.5} $
      112.1 $ \pm $ 3.27.07 $ \pm $ 3.651.92$ ^{+0.17}_{-0.11} $1.84$ ^{+0.18}_{-0.14} $3.6$ ^{+0.3}_{-0.4} $
      124.4 $ \pm $ 2.77.49 $ \pm $ 3.456.54$ ^{+0.44}_{-0.64} $7.60$ ^{+4.02}_{-3.74} $6.4$ ^{+0.6}_{-1.0} $
      135.1$ \pm $2.07.07 $ \pm $ 3.656.77$ ^{+1.9}_{-1.35} $8.24$ ^{+3.24}_{-2.87} $9.7$ ^{+3.3}_{-2.3} $

      Table 4.  The cooling age of the 13 RGB stars selected [41]. The five cooling ages correspond to models I $-$ V. The discussions of the cooling ages are based on the packages BASTA [42], isochrones [43], isoclassify [44], PARAM [45, 46], and scaling-giants [47].

      Based on Eqs. (7)–(9), there are two global seismic parameters, the frequency of maximum power, $ \nu_{\rm{max}} $, and the mean large frequency separation, $ \bigtriangleup \nu $ , to describe the oscillations of solar-like stars and radius. The surface gravity and temperature strongly determine the value of the frequency of maximum power, which is given by $ \nu_{\rm{max}}\propto gT^{-0.5}_{\rm{eff}}\propto MR^{-2}T^{-0.5}_{\rm{eff}} $ [48]. On the other hand, the travel time of sound from the center to the surface of a star will directly influence $ \bigtriangleup \nu $, which is sensitive to the mean stellar density and is given by $ \bigtriangleup \nu\propto \rho^{0.5}\propto M^{0.5}R^{-1.5} $ [49]. In order to reduce systematic errors, by considering the effect on mass, [Fe/H] and $ T_{\rm{eff}} $ interpolation over a grid of models, a modification strategy was adopted by Sharma et al. [37]. Therefore, the asteroseismic correction factor in this paper is selected as $ f_{\bigtriangleup \nu}= $ 0.98, 0.99, 1.00, 1.02, we propose five model (I$-$V). Detailed information on the mass and radius of RGB stars is shown in Table 3.

      Based on Eqs. (10) and (12), the age of RGB stars is a key parameter for estimating the number of MMs captured and the luminosity of WDs due to the RC effect by MMs. By combining masses inferred from the asteroseismic parameters with stellar atmospheric parameters and using stellar isochrones, the cooling ages of RGB stars were determined as shown in Table 4 by $ t_9 $. Using the packages BASTA [42], isochrones [43], isoclassify [44], PARAM [45, 46], and scaling-giants [47], estimates are obtained of the cooling ages of RGB stars, which are the stellar ages since zero-age main sequence stars are considered.

      The scaling-giants package accepts asteroseismic parameters, metallicity, and temperature as inputs for model I. In the isochrones, isoclassify, PARAM, and BASTA packages, we take asteroseismic, photometric, and spectroscopic parameters as inputs for models II−V. From the SYD pipeline [50], along with effective temperatures determined through the direct method of isoclassification and metallicity determined by either the APOGEE or GALAH surveys, and using measured seismic $ \nu_{\rm{max}} $ and $ \bigtriangleup\nu $ values, we can determine the ages with scaling-giants for model I. The SYD pipeline is an automated pipeline to estimate global oscillation parameters, such as the frequency of maximum power ($ \nu_{\rm{max}} $) and the large frequency spacing ($ \triangle\nu $), for a large number of time series. By using 2MASS K magnitudes, asteroseismic $ \bigtriangleup\nu $, Gaia parallaxes, and temperatures and metallicities from spectroscopy, we can determine the ages with isochrones, isoclassify, and PARAM, and BASTA packages for models II−V, respectively. Some useful parameters are presented in Tables 13. Detailed discussions are given in Grunblatt et al. [41].

      Based on Eq. (10), we give a computational assessment for the number of MMs captured for WDs in Tables 5 and 6 for $ m=10^{15}, 10^{17} $GeV, respectively. Our results in Table 5 show that the maximum number of MMs captured in models I $ \sim $V are $ 4.8126\times10^{22} $ and $ 6.5345\times10^{24} $ when $ \phi=5.59\times10^{-28} \;\rm{cm^{-2}s^{-1}sr^{-1}} $, and $ \phi=7.59\times 10^{-26} \rm{cm^{-2}s^{-1}sr^{-1}} $, respectively. Table 6 also shows that the number of MMs captured has a maximum of $ 9.1223\times10^{24} $ (e.g., when $m=10^{17}{\rm GeV}, n_{\rm B}=5.99\times 10^{-31} \rm{cm^{-3}}, \phi=7.59\times10^{-26} \rm{cm^{-2}s^{-1} sr^{-1}}$, for model III). One can see that the number of MMs captured increases as the flux of MMs increases due to $ N_m\propto \phi $ according to Eq. (10). On the other hand, when the flux and mass of MMs are certain, there is no significant difference found in the number of MMs captured for same stars among the five models. However, for different models, there is a small difference. The reasons for these differences can be from the differences of the cooling ages for RGB stars due to differences in the stellar parameters selected, such as the asteroseismic correction factor and $ [\rm{Fe/H}] $. From the above analysis, it may be seen that the number of MMs captured can reach the maximum value of $ 9.1223\times10^{24} $ when $m=10^{17}{\rm GeV}, n_{\rm B}= 5.99\times10^{-31}\rm{cm^{-3}}, \phi= 7.59\times10^{-26} \rm{cm^{-2}s^{-1}sr^{-1}}$.

      $ \phi=5.59\times10^{-28}\rm{cm^{-2}s^{-1}sr^{-1}} $$ \phi=7.59\times10^{-26}\rm{cm^{-2}s^{-1}sr^{-1}} $
      star$ N_m $(I)$ N_m $(II)$ N_m $(III)$ N_m $(IV)$ N_m $(V)$ N_m $(I)$ N_m $(II)$ N_m $(III)$ N_m $(IV)$ N_m $(V)
      12.5599e214.1144e212.8307e214.9152e215.0524e213.4758e235.5864e233.8435e236.6738e236.8601e23
      27.2823e212.4449e222.7645e221.9626e222.5231e229.8877e233.3197e243.7536e242.6648e243.4258e24
      32.4095e213.2314e213.2591e213.5154e214.8906e213.2716e234.3876e234.4251e234.7731e236.6404e23
      46.6993e216.5919e219.2068e218.9176e211.2010e229.0962e238.9503e231.2501e241.2108e241.6307e24
      52.6667e223.5386e224.8126e223.8356e224.0272e223.6207e244.8046e246.5345e245.2079e245.4680e24
      64.9402e207.1828e206.9014e204.6290e206.7508e206.7077e229.7527e229.3705e226.2851e229.1661e22
      74.4692e211.2594e221.7647e221.4790e229.9908e216.0681e231.7100e242.3961e242.0081e241.3565e24
      88.2517e201.0349e211.0989e211.2401e211.9185e211.1204e231.4051e231.4921e231.6837e232.6049e23
      91.5855e221.6237e221.6212e221.4106e222.5275e222.1527e242.2046e242.2012e241.9152e243.4318e24
      105.7326e211.1862e221.3458e221.1338e221.5189e227.7836e231.6105e241.8273e241.5395e242.0624e24
      113.2438e211.1608e223.3488e213.4071e217.0729e214.4044e231.5761e244.5469e234.6261e239.6034e23
      124.3369e217.8473e217.2788e218.9801e218.0238e215.8885e231.0655e249.8830e231.2193e241.0895e24
      131.2784e211.8839e211.9164e212.4764e213.0932e211.7359e232.5579e232.6020e233.3624e234.1999e23

      Table 5.  The number of MMs captured in the five models I$-$V (corresponding to typical asteroseismic correction factor values $ f_{\bigtriangleup\nu}=0.98, 0.99, 1.00, 1.01, 1.02 $ ) when $m=10^{15}{\rm GeV}, n_{\rm B}=5.99\times10^{31}\rm{cm^{-3}}, \phi=5.59\times10^{-28}, 7.59\times10^{-26}\rm{cm^{-2}s^{-1}sr^{-1}}$.

      $ \phi=5.59\times10^{-28}\rm{cm^{-2}s^{-1}sr^{-1}} $$ \phi=7.59\times10^{-26}\rm{cm^{-2}s^{-1}sr^{-1}} $
      star$ N_m $(I)$ N_m $(II)$ N_m $(III)$ N_m $(IV)$ N_m $(V)$ N_m $(I)$ N_m $(II)$ N_m $(III)$ N_m $(IV)$ N_m $(V)
      13.9330e216.3662e214.4113e217.7145e217.9867e215.3401e238.6439e235.9895e231.0475e241.0844e24
      21.0794e223.6479e224.1522e222.9675e223.8405e221.4656e244.9531e245.6377e244.0292e245.2146e24
      33.7921e215.1235e215.2058e215.6570e217.9291e215.1488e236.9565e237.0683e237.6810e231.0766e24
      49.7253e219.6300e211.3536e221.3195e221.7884e221.3205e241.3075e241.8379e241.7915e242.4282e24
      53.6811e224.9121e226.7186e225.3850e225.6862e224.9981e246.6696e249.1223e247.3117e247.7206e24
      68.9937e201.3196e211.2795e218.6609e201.2747e211.2212e231.7918e231.7373e231.1760e231.7307e23
      76.8316e211.9388e222.7358e222.3091e221.5709e229.2759e232.6324e243.7146e243.1352e242.1330e24
      81.4320e211.8114e211.9401e212.2082e213.4458e211.9443e232.4595e232.6342e232.9982e234.6787e23
      92.2008e222.2668e222.2762e221.9920e223.5898e222.9883e243.0778e243.0906e242.7046e244.8742e24
      108.3230e211.7331e221.9789e221.6778e222.2622e221.1301e242.3532e242.6868e242.2781e243.0715e24
      114.9599e211.7875e225.1930e215.3209e211.1124e226.7345e232.4270e247.0510e237.2247e231.5105e24
      126.8084e211.2410e221.1597e221.4414e221.2975e229.2444e231.6851e241.5746e241.9571e241.7617e24
      132.1582e213.2067e213.2891e214.2856e215.3977e212.9304e234.3540e234.4659e235.8190e237.3289e23

      Table 6.  The number of MMs captured in the five models I$-$V (corresponding to typical asteroseismic correction factor values $ f_{\Delta\nu}=0.98, 0.99, 1.00, 1.01, 1.02 $ ) when $m=10^{17}{\rm GeV}, n_{\rm B}=5.99\times10^{31}\rm{cm^{-3}}, \phi=5.59\times10^{-28}, 7.59\times10^{-26}\rm{cm^{-2}s^{-1}sr^{-1}}$.

      Figures 1 and 2 display the luminosities as a function of $ t_9 $ of WDs of the five models under different astronomical conditions. From the two figures, it can be seen that the calculated luminosities agree well with observations. One can also conclude the same from Table 7. For example, the ranges of our calculated luminosities are $ 1.3336\times10^{34}\sim 7.1985\times10^{36}\rm{erg \; s^{-1}} $ and $ 1.8224\times10^{34}\sim 1.0871\times10^{36}\rm{erg\; s^{-1}} $ for models I and V, respectively, in Table 7.

      Figure 1.  (color online) The luminosities as a function of $ t_9 $ for the 13 WDs [41] for the five models I $-$V when $n_{\rm B}=5.99\times10^{31}\rm{cm^{-3}}$, $ m=10^{15} $GeV, $ \sigma_m=10^{-24}\rm{cm^2}, \phi=10^{-26}\rm{cm^{-2}s^{-1}sr^{-1}} $ at the temperature$ T_6=1 $.

      Figure 2.  (color online) The luminosities as a function of $ t_9 $ of the 13 WDs [41] for models I $-$ V when $n_{\rm B}=5.99\times10^{31}\rm{cm^{-3}}$, $ m=10^{15} $GeV, $ \sigma_m=10^{-24}\rm{cm^2}, \phi=10^{-26}\rm{cm^{-2}s^{-1}sr^{-1}} $ at the temperature $ T_6=10 $.

      starour resultsobservations
      $L_{m }$(I)$L_{m }$(II)$L_{m }$(III)$L_{m }$(IV)$L_{m }$(V)$ L_{\rm{rad}} $(I)$ L_{\rm{rad}} $(II)$ L_{\rm{rad}} $(III)$ L_{\rm{rad}} $(IV)$L_{m }$(V)
      16.9103e341.1107e357.6414e341.3269e351.3639e351.0318e351.0746e351.1187e351.1641e351.2109e35
      21.9658e356.6000e357.4626e355.2980e356.8110e352.4869e352.5900e352.6962e352.8057e352.9185e35
      36.5043e348.7231e348.7978e349.4896e341.3202e351.2778e351.3308e351.3854e351.4417e351.4996e35
      41.8085e351.7795e352.4853e352.4073e353.2420e359.2913e349.6763e341.0073e351.0482e351.0904e35
      57.1985e359.5522e351.2992e361.0354e361.0871e361.6913e351.7614e351.8336e351.9081e351.9848e35
      61.3336e341.9390e341.8630e341.2496e341.8224e347.7022e348.0214e348.3505e348.6895e349.0388e34
      71.2064e353.3998e354.7637e353.9925e352.6970e351.7814e351.8552e351.9313e352.0097e352.0905e35
      82.2275e342.7936e342.9665e343.3475e345.1790e346.3515e346.6147e346.8861e347.1657e347.4537e34
      94.2799e354.3831e354.3763e353.8078e356.8228e352.3176e352.4137e352.5127e352.6147e352.7198e35
      101.5475e353.2020e353.6329e353.0607e354.1003e358.3978e348.7458e349.1046e349.4742e349.8551e34
      118.7565e343.1336e359.0399e349.1974e341.9093e351.6398e351.7078e351.7779e351.8501e351.9244e35
      121.1707e352.1183e351.9649e352.4242e352.1660e351.1526e351.2004e351.2496e351.3004e351.3526e35
      133.4511e345.0855e345.1732e346.6849e348.3500e347.3599e347.6649e347.9794e348.3034e348.6371e34

      Table 7.  Comparisons of the calculated luminosities ($L_{m}$) due to the RC effect by MMs with observations ($ L_{\rm{rad}} $) for the five models I $-$ V at $n_{\rm B}=5.99\times10^{31}{\rm cm^{-3}}, m=10^{15}$GeV, $\sigma_m=10^{-24}{\rm cm^2}, \phi=10^{-26}{\rm cm^{-2}s^{-1}sr^{-1}}, T_6=1$ ($ \xi=1.5745\times10^{-15}\rm{cm^{-2}s^{-1}sr^{-1}} $).

      In order to facilitate comparison of our results with the observed data, we define the scale factor $ k_i $as the ratio of our results due to the RC effect to the observed luminosity. We note that the largest differences between our results and the observed values are one order of magnitude. For example, based on Table 7 and Table 9 for star 5, the maximal scale factors are $ k_1=4.2562 $, $ k_2=5.4231 $, $ k_3=7.0851 $, $ k_4=5,4264 $, and $ k_5=5.4772 $ for models I−V, respectively. However, in Table 8 and Table 9, the maximal scale factors for star 5 are $ k_1=13.459 $, $ k_2=17.149 $, $ k_3=22.405 $, $ k_4=17.160 $ and $ k_5=17.321 $ for models I−V, respectively.

      starTable 7Table 8
      $ k_1 $$ k_2 $$ k_3 $$ k_4 $$ k_5 $$ k_1 $$ k_2 $$ k_3 $$ k_4 $$ k_5 $
      10.669731.03360.683091.13981.12642.11793.26852.16013.60453.5619
      20.790482.54832.76781.88832.33382.49978.05858.75275.97147.3801
      30.50900.65550.63500.65820.88041.60962.07282.00822.08152.7840
      41.94641.83902.46732.29652.97336.15515.81537.80227.26239.4024
      54.25625.42317.08515.42645.477213.45917.14922.40517.16017.321
      60.17310.24170.22310.14380.20160.547530.76440.70550.45470.6376
      70.67731.83262.46661.98661.29012.14165.79517.80006.28214.0797
      80.35070.422330.43080.46720.69481.10901.33551.36231.47732.1972
      91.84671.81591.74171.45632.50855.83965.74255.50764.60527.9327
      101.84273.66123.99023.23064.16065.827311.57812.61810.21613.157
      110.53401.83480.50850.49710.99221.68865.80231.60791.57213.1374
      121.01571.76471.57241.86421.60133.21205.58064.97235.89525.0638
      130.46890.66350.64830.80510.96681.48282.09812.05022.54593.0572

      Table 9.  The comparisons of the scale factor $ k_i $ ($ i=(1-5) $ indexes models I−V), which is the ratio of our results ($L_{m}$) due to RC effect by MMs to the observed luminosities ($ L_{\rm{rad}} $) for the five models I−V from Tables 7 and 8.

      starour resultsobservations
      $L_{m }$(I)$L_{m }$(II)$L_{m }$(III)$L_{m }$(IV)$L_{m }$(V)$ L_{\rm{rad}} $(I)$ L_{\rm{rad}} $(II)$ L_{\rm{rad}} $(III)$ L_{\rm{rad}} $(IV)$L_{m }$(V)
      12.1852e353.5122e352.4164e354.1959e354.3130e351.0318e351.0746e351.1187e351.1641e351.2109e35
      26.2165e352.0871e362.3599e361.6754e362.1538e362.4869e352.5900e352.6962e352.8057e352.9185e35
      32.0568e352.7585e352.7821e353.0009e354.1749e351.2778e351.3308e351.3854e351.4417e351.4996e35
      45.7189e355.6271e357.8593e357.6125e351.0252e369.2913e349.6763e341.0073e351.0482e351.0904e35
      52.2764e363.0207e364.1083e363.2743e363.4378e361.6913e351.7614e351.8336e351.9081e351.9848e35
      64.2172e346.1316e345.8913e343.9515e345.7628e347.7022e348.0214e348.3505e348.6895e349.0388e34
      73.8151e351.0751e361.5064e361.2625e368.5286e351.7814e351.8552e351.9313e352.0097e352.0905e35
      87.0440e348.8342e349.3808e341.0586e351.6377e356.3515e346.6147e346.8861e347.1657e347.4537e34
      91.3534e361.3861e361.3839e361.2041e362.1576e362.3176e352.4137e352.5127e352.6147e352.7198e35
      104.8936e351.0126e361.1488e369.6788e351.2966e368.3978e348.7458e349.1046e349.4742e349.8551e34
      112.7691e359.9092e352.8587e352.9085e356.0377e351.6398e351.7078e351.7779e351.8501e351.9244e35
      123.7022e356.6988e356.2135e357.6659e356.8495e351.1526e351.2004e351.2496e351.3004e351.3526e35
      131.0913e351.6082e351.6359e352.1140e352.6405e357.3599e347.6649e347.9794e348.3034e348.6371e34

      Table 8.  The comparisons of the luminosity ($L_{m}$) due to RC effect by MMs with the observed values ($ L_{\rm{rad}} $) for the five models I $-$ V at $n_{\rm B}=5.99\times10^{31}{\rm cm^{-3}}, m=10^{15}$GeV, $\sigma_m=10^{-24}{\rm cm^2}, \phi=10^{-26}{\rm cm^{-2}s^{-1}sr^{-1}}, T_6=10$ ($ \xi=4.9790\times10^{-15}\rm{cm^{-2}s^{-1}sr^{-1}} $).

      According to the analysis above, our results agree well with observations at lower relativistic densities and temperatures, but the greatest difference is about two orders of magnitude at higher relativistic densities and temperatures in WDs. The good agreement of our calculatedresults with observations shows that our model with the RC effect by MMs yields realistic results, as well as suggesting that the energy source of WDs is the RC effect by MMs.

      The monopole flux problem is well known to be a highly challenging and interesting issue. Some scholars have made pioneering work on this subject, such as Freese [32], Parker [39], Kolb & Turner [31], and Freese & Krasteva [40]. In this paper, we discuss this problem by considering the RC effect of MMs. Tables 10 ($ m=10^{15}\rm{GeV}, \beta=9.4868\times10^{-3} $) and 11 ($ m=10^{17}\rm{GeV}, \beta=1.00\times10^{-3} $) show the flux of MMs for the four models I$ \sim $III and V, which correspond to typical asteroseismic correction factor values $ f_{\bigtriangleup\nu}=0.98, 0.99, 1.00, 1.02 $ when $n_{\rm B}=5.99\times10^{31}, 1.89\times10^{32}\rm{cm^{-3}}$. One can see that the maximum of the monopole flux values are $ 9.0935\times10^{-13}\rm{cm^{-2}s^{-1}sr^{-1}} $ and $ 5.8519\times10^{-13} \rm{cm^{-2}s^{-1}sr^{-1}} $ in Tables 10 and 11, respectively.

      star$ n_B=5.99\times10^{31}\rm{cm^{-3}} $$ n_B=1.89\times10^{32}\rm{cm^{-3}} $
      ξ(I)ξ(II)ξ(III)ξ(IV)ξ(V)ξ(I)ξ(II)ξ(III)ξ(IV)ξ(V)
      12.3509e-131.5233e-132.3050e-131.3813e-131.3978e-132.3614e-141.5301e-142.3152e-141.3875e-141.4041e-14
      21.9918e-136.1786e-145.6885e-148.3380e-146.7465e-142.0007e-146.2061e-155.7139e-158.3752e-156.7765e-15
      33.0933e-132.4021e-132.4794e-132.3920e-131.7884e-133.1070e-142.4128e-142.4904e-142.4026e-141.7964e-14
      48.0892e-148.5618e-146.3815e-146.8559e-145.2954e-148.1252e-158.5999e-156.4099e-156.8865e-155.3190e-15
      53.6993e-142.9033e-142.2223e-142.9015e-142.8746e-143.7157e-152.9162e-152.2322e-152.9145e-152.8874e-15
      69.0935e-136.5135e-137.0573e-136.0949e-137.8094e-139.1340e-146.5425e-147.0887e-147.0998e-147.8442e-14
      72.3248e-138.5916e-146.3833e-147.9257e-141.2204e-132.3352e-148.6299e-156.4117e-157.9610e-151.2259e-14
      84.4895e-133.7281e-133.6549e-133.3703e-132.2661e-134.5095e-143.7447e-143.6711e-143.3854e-142.2761e-14
      98.5262e-148.6703e-149.0402e-141.0812e-136.2765e-148.5642e-158.7090e-159.0804e-151.0860e-146.3044e-15
      108.5442e-144.3005e-143.9459e-144.8737e-143.7843e-148.5823e-154.3196e-153.9635e-154.8954e-153.8011e-15
      112.9486e-138.5810e-143.0965e-133.1671e-131.5870e-132.9617e-148.6192e-153.1103e-143.1812e-141.5940e-14
      121.5501e-138.9220e-141.0013e-138.4458e-149.8324e-141.5570e-148.9617e-151.0058e-148.4834e-159.8762e-15
      133.3578e-132.3731e-132.4285e-131.9557e-131.6286e-133.3727e-142.3837e-142.4394e-141.9644e-141.6359e-14

      Table 10.  The flux of MMs for the five models I−V (corresponding to typical asteroseismic correction factor values $ f_{\Delta\nu}= $ 0.98, 0.99, 1.00, 1.01, 1.02) when $m=10^{15}{\rm GeV}, \beta=9.4868\times10^{-3}, n_{\rm B}=5.99\times10^{31}, 1.89\times10^{32}\rm{cm^{-3}}$.

      star$n_{\rm B}=5.99\times10^{31}\rm{cm^{-3} }$$n_{\rm B}=1.89\times10^{32}\rm{cm^{-3} }$
      ξ(I)ξ(II)ξ(III)ξ(IV)ξ(V)ξ(I)ξ(II)ξ(III)ξ(IV)ξ(V)
      11.5302e-139.8450e-141.4791e-138.8011e-148.8428e-141.5370e-149.8889e-151.4857e-148.8403e-158.8822e-15
      21.3438e-134.1410e-143.7874e-145.5146e-144.4323e-141.3498e-144.1595e-153.8042e-155.5391e-154.4520e-15
      31.9655e-131.5150e-131.5522e-131.4864e-131.1031e-131.9742e-141.5218e-141.5591e-141.4930e-141.1080e-14
      45.5723e-145.8607e-144.3406e-144.6336e-143.5561e-145.5971e-155.8868e-154.3599e-154.6543e-153.5720e-15
      52.6798e-142.0914e-141.5918e-142.0667e-142.0359e-142.6918e-152.1008e-151.5989e-152.0759e-152.0450e-15
      64.9950e-133.5454e-133.8064e-135.8519e-134.1359e-135.0173e-143.5612e-143.8234e-145.8779e-144.1543e-14
      71.5209e-135.5812e-144.1175e-145.0764e-147.7618e-141.5276e-145.6060e-154.1359e-155.0990e-157.7963e-15
      82.5870e-132.1299e-132.0702e-131.8927e-131.2617e-132.5985e-142.1394e-142.0794e-141.9011e-141.2673e-14
      96.1421e-146.2105e-146.4385e-147.6561e-144.4190e-146.1695e-156.2382e-156.4672e-157.6902e-154.4387e-15
      105.8850e-142.9433e-142.6835e-143.2935e-142.5410e-145.9112e-152.9564e-152.6955e-153.3081e-152.5523e-15
      111.9284e-135.5727e-141.9968e-132.0279e-131.0090e-131.9370e-145.5975e-152.0057e-142.0370e-141.0135e-14
      129.8741e-145.6415e-146.2849e-145.2620e-146.0805e-149.9181e-155.6666e-156.3129e-155.2854e-156.1076e-15
      131.9890e-131.3942e-131.4150e-131.1301e-139.3331e-141.9979e-141.4004e-141.4213e-141.1351e-149.3746e-15

      Table 11.  The flux of MMs for the five models I−V (corresponding to typical asteroseismic correction factor values $ f_{\Delta\nu}= $ 0.98, 0.99, 1.00, 1.01, 1.02) when $m=10^{17}{\rm GeV}, \beta=1.0\times10^{-3}, n_{\rm B}=5.99\times10^{31}, 1.89\times10^{32}\rm{cm^{-3}}$.

      It is very interesting to note that the monopole flux decreases as $n_{\rm B}$ increases from $ 5.99\times10^{31} $ to $ 1.89\times 10^{32}\rm{cm^{-3}} $ in Tables 10 and 11. This is not hard to understand according to Eq. (10). Based on our calculations above from Tables 10, due to the RC effect by MMs, we obtain new limits on the MM flux of $\xi\leq 9.0935\times 10^{-13}\rm{cm^{-2}s^{-1}sr^{-1}}$ and $\xi\leq 9.1340\times 10^{-14} \rm{cm^{-2}s^{-1}sr^{-1}}$ when $n_{\rm B}=5.99\times10^{31} \rm{cm^{-3}}$ and $ 1.899\times 10^{32}\rm{cm^{-3}} $, respectively. In Table 11, the new limits on the MM flux are $\xi\leq 4.9950 \times10^{-13}\rm{cm^{-2}s^{-1}sr^{-1}}$ and $\xi\leq 5.0173\times10^{-14} \rm{cm^{-2}s^{-1}sr^{-1}}$ when $n_B=5.99\times10^{31}$ and $ 1.899\times10^{32} $, respectively.

      Based on the above analysis, we obtain new limits on the MM flux of $\xi\leq 9.0935\times10^{-13} \rm{cm^{-2}s^{-1}sr^{-1}}$, and $ \xi\leq 4.9950\times10^{-13}\rm{cm^{-2}s^{-1}sr^{-1}} $ at $n_{\rm B}=5.99\times10^{31}\rm{cm^{-3}}$ when $ m=10^{15}\rm{GeV}, \beta=9.4868\times10^{-3} $, and $ m=10^{17}\rm{GeV}, \beta=10^{-3} $, respectively. When we estimate the number of MMs captured, the MM luminosity and the limit of the MM flux, as samples, we consider the 13 RGB stars in our MM model for the following reasons. First, WDs originate from RGB stars. Second, compared with WDs, RGB stars have a very large surface area. According to Eq. (10), we expect RGB phases to capture more MMs during their evolution period. Third, since the MM is a superheavy particle, when MMs are captured by an RGB, they will be deposited in the star core. If all MMs captured by RGBs remain, their number will be much larger than that of those captured by WDs. Thus, the number of MMs calculated inside an WD will be more accurate than for MMs captured only during the WD phase. For example, Freese et al. [51] showed that if the MMs captured by stars in the main sequence stage all survive, the MM flow due to neutron star catalysis can be strengthened by up to 7 orders of magnitude. Finally, based on Schwarzschid [52], the nuclear energy generation rates of the proton-proton and CNO cycle are $ \epsilon_{\rm{pp}}\approx10\rho_{100}T_7^4 \rm{erg\; g^{-1}s^{-1}} $and $ \epsilon_{\rm{CNO}}\approx8\rho_{100}T_7^{16} \rm{erg\; g^{-1}s^{-1}} $, respectively (where $ T_7=T/10^7\rm{K} $ is the temperature, and $ \rho_{100}=\rho/100 $ is the density). Based on the discussions fof Bjork et al. [53], we may select the mass of the outer layer of the RGB as being from $0.005\sim0.02 {M}_{\odot}$ (the main component is hydrogen). Thus, when $T_7=0.1-1$, $ \rho_{100}\sim10^{-4} $, and we obtain a proton-proton nuclear energy generation rate of $10^{24} - 10^{28} \rm{ergs\; s^{-1}}$, which is $\ll L_m=10^{34} - 10^{36}\rm{ergs\; s^{-1}}$ in our calculations. Based on the above analysis, and the fact that RGB stars are the origins of WDs, we therefore have $ L_{\rm{m}}\approx L_{\rm{rad}} $.

      One can also conclude that with the increasing number of MMs captured, the luminosity due to the RC effect by MMs increases linearly with time until it becomes the main contribution to the total luminosity. One can even observe that for some of the oldest white dwarfs, the luminosity may have passed its minimum and some reheating may have occurred.

      It may be suggested that the annihilation of magnetic and antimagnetic monopoles could result in a significant reduction in the number of monopoles and the catalytic luminosity of the monopoles in the WDs. Dicus et al. [54] calculated the annihilation cross sections of magnetic monopoles and antimonopoles caused by two-body and three-body recombination. Their results show that this kind of annihilation has little effect on the flux and luminosity of the monopole. On the other hand, some WDs may have magnetic fields of up to $ 10^5 $G according to observations. The forces generated by the magnetic field inside the white dwarf must balance the gravitational and Coulomb interactions.The magnetic field may keep the monopole and antimonopole distributions far enough apart for annihilation to be negligible.

      On the other hand, neutrino emission in WDs is a very interesting issue. Based on the discussions of Althaus et al. [55], when WD is very hot, neutrino emissions could be a major source of cooling. However, based on the relatively low temperature environment of WDs (e.g., $ T_6=1, 10 $ in our paper) we discussed the heating resource problem with our MMs model. The neutrino processes inside WDs at such low temperatures (e.g., $ T_6=1 $) may not be the main cooling process (see discussions by Itoh et al. [56]). On the other hand, Izawa [57] also discussed the neutrinos emitted according to Eqs. (23), (24) in his paper and calculated the neutrino emitted per one nucleon decay at low and high energy components, finding that these neutrino losses did not affect the structure or evolution of Rubakov stars because the energy lost through the neutrino emission is smaller than 100 MeV per one nucleon decayn although about two neutrinos are emitted furing the decay of one nucleon.

      According to the above analysis, one can see that MMs pass through space to be captured by WDs. MMs trapped inside a WD can catalyze the decay of nuclei, which can function as an energy source to keep the WDs hot.

    V.   CONCLUSIONS AND OUTLOOKS
    • We have presented five MMs models of WD energy resources to discuss their cooling based on certain observations of 13 RGB stars. We find that the number of MMs captured can reach a maximum value of $ 9.1223\times10^{24} $ when $m=10^{17}\;\;{\rm GeV}, \; n_{\rm B}=5.99\times10^{31}\; \rm{cm^{-3}}, \phi=7.59\times10^{-26}\;\;\rm{cm^{-2}s^{-1}sr^{-1}}$. The good agreement with observations of our luminosities due to the RC effect by MMs calculated for WDs shows that our model is reasonable. We conclude that the energy source of WDs may be the RC effect. Due to the RC effect by MMs, we obtain a new limit of MM flux of $\xi\leq 9.0935\times 10^{-13} \;\; \rm{cm^{-2}s^{-1}sr^{-1}}$ and $\xi\leq 4.9950\times10^{-13}\;\;\rm{cm^{-2}s^{-1}sr^{-1}}$ at $n_{\rm B}=5.99\times 10^{31} \rm{cm^{-3}}$ when $m=10^{15}\rm{GeV}, \;\beta=9.4868\times10^{-3}$, and $m=10^{17}\rm{GeV}, \;\beta=10^{-3}$, respectively.

      In this paper, the main highlights may be given as follows. First, we created detailed estimates of the cooling ages of 13 RGB stars using the packages BASTA [42], isochrones [43], isoclassify [44], PARAM [45,46], and scaling-giants [47]. Second, we proposes five new models to discuss the energy resources and the cooling of WDs and compare the luminosities with observations for 13 RGB stars due to the RC effect. Finally, the new limit of the MM flux is obtained based on our models.

      As is widely known, research on MMs haa always been a hot frontier topic in the fields of nuclear physics and astrophysics. The search for MMs remians a difficult and challenging problem, and the flux of magnetic monopoles in the universe remains uncertain. The neutrino emissivity rates due to the RC effect also may play a key role in the process of WD and neutron star evolution. These challenging problem will be our future issues.

Reference (57)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return