Stars of $\sim$ 8 - 10 $M_{\odot}$ on the main-sequence form strongly electron-degenerate O+Ne+Mg core and become super-AGB stars. If such an O+Ne+Mg core grows to 1.38 $M_\odot$, electron captures on $^{20}$Ne$(e,\nu_e)^{20}$F$(e,\nu_e)^{20}$O take place and ignite O-Ne deflagration around the center. In this paper, we perform two-dimensional hydrodynamics simulations of the propagation of the O-Ne flame to see whether such a flame induces a collapse of the O+Ne+Mg core due to subsequent electron capture behind the flame or triggers a thermonuclear explosion. We present a series of models to explore how the outcome depends on model parameters for the central density in the range from $10^{9.80}$ to $10^{10.20}$ g cm$^{-3}$, flame structure of both centered and off-centered ignition kernels, special and general relativistic effects, turbulent flame speed formula and the treatments of laminar burning phase. We find that the O+Ne+Mg core obtained from stellar evolutionary models has a high tendency to collapse into a neutron star. We obtain the bifurcation between the electron-capture collapse and thermonuclear explosion. We discuss the implication in nucleosynthesis and the possible observational signals of this class of supernovae.

Recent observations of Type Ia supernovae (SNe Ia) have shown diversified properties of the explosion strength, light curves and chemical composition. To investigate possible origins of such diversities in SNe Ia, we have presented multi-dimensional hydrodynamical study of explosions and associated nucleosynthesis in the near Chandrasekhar mass carbon-oxygen (CO) white dwarfs (WDs) for a wide range of parameters (Leung and Nomoto 2018 ApJ). In the present paper, we extend our wide parameter survey of models to the explosions of sub-Chandrasekhar mass CO WDs. We take the double detonation model for the explosion mechanism. The model parameters of the survey include the metallicity of $Z = 0 - 5~Z_\odot$, the CO WD mass of $M = 0.90 - 1.20~M_\odot$, and the He envelope mass of $M_{\rm He} = 0.05 - 0.20~M_\odot$. We also study how the initial He detonation configuration, such as spherical, bubble, and ring shapes, triggers the C detonation. For these parameters, we derive the minimum He envelope mass necessary to trigger the C detonation. We then examine how the explosion dynamics and associated nucleosynthesis depend on these parameters, and compare our results with the previous representative models. We compare our nucleosynthesis yields with the unusual abundance patterns of Fe-peak elements and isotopes observed in SNe Ia 2011fe, 2012cg and 2014J, as well as SN Ia remnant 3C 397 to provide constraints on their progenitors and environments. We provide the nucleosynthesis yields table of the sub-Chandrasekhar mass explosions, to discuss their roles in the galactic chemical evolution and archaeology.

Recent observations have shown a remarkable diversity of the observational behaviors and explosion mechanisms in thermonuclear supernovae (SNe). An emerging class of peculiar type Ia SNe, called type Iax, show distinct photometric and spectroscopic behaviors. Their origin remains highly controversial, but pure turbulent deflagration of white dwarfs (WDs) has been regarded as the leading formation theory. The large population of type Iax indicates the existence of unidentified Galactic type Iax supernova remnants (SNRs). We report evidence that SNR Sgr A East in the Galactic center resulted from a pure turbulent deflagration of a Chandrasekhar-mass carbon-oxygen WD, an explosion mechanism used for type Iax SNe. Our X-ray spectroscopic study of Sgr A East using 3 Ms Chandra data shows a low ratio of intermediate-mass elements to Fe, and large Mn/Fe and Ni/Fe ratios. This abundance pattern does not accord with the core-collapse or normal type Ia models. Sgr A East is thus the first Galactic SNR for which a likely type Iax origin has been proposed, and the nearest target to study this peculiar class. We compared Sgr A East with the Fe-rich SNRs 3C 397 and W49B, which also have high Mn and Cr abundances and were claimed to result from deflagration-to-detonation explosions of Chandrasekhar-mass WDs (although with disputes). Our study shows that they have distinct abundance patterns. The X-ray spectroscopic studies of thermonuclear SNRs provide observational evidence for the theories that there are diverse explosion channels and various metal outputs for Chandrasekhar-mass WDs.

We conduct a rigorous examination of the nearby red supergiant $\alpha$ Orionis, or Betelgeuse, by drawing on the synthesis of new observational data and three different modeling techniques. Our observational results include the release of new, processed photometric measurements collected with the space-based SMEI instrument prior to Betelgeuse's recent, unprecedented dimming event. Our theoretical predictions include self-consistent results from multi-timescale evolutionary, oscillatory, and hydrodynamic simulations conducted with the Modules for Experiments in Stellar Astrophysics (MESA) software suite. Significant outcomes of our modeling efforts include a precise prediction for the star's radius: $750^{+62}_{-30}\,R_{\odot}$. In concert with additional constraints, this allows us to derive a new, independent distance estimate of $165^ {+16}_{-8}$\,pc and a parallax of $\pi=6.06^{+0.31}_{-0.52}$\,mas, in good agreement with \textit{Hipparcos} but less so with recent radio measurements. Seismic results from both perturbed hydrostatic and evolving hydrodynamic simulations constrain the period and driving mechanisms of Betelgeuse's dominant periodicities in new ways. Our analyses converge to the conclusion that Betelgeuse's $388$ day period is the result of pulsation in the fundamental mode, driven by the $\kappa$-mechanism. Grid-based hydrodynamic modeling reveals that the behavior of the oscillating envelope is mass-dependent, and likewise suggests that the non-linear pulsation excitation time could serve as a mass constraint. Our results corroborate recent conclusions that Betelgeuse is the outcome of a past merger. We place it definitively in the core helium-burning phase near the base of the red supergiant branch, and we report a present-day mass of $16.5$--$19 ~M_{\odot}$---slightly lower than typical literature values.

We present 3D full-sphere supernova simulations of non-rotating low-mass (~9 Msun) progenitors, covering the entire evolution from core collapse through bounce and shock revival, through shock breakout from the stellar surface, until fallback is completed several days later. We obtain low-energy explosions [~(0.5-1.0)x 10^{50} erg] of iron-core progenitors at the low-mass end of the core-collapse supernova (LMCCSN) domain and compare to a super-AGB (sAGB) progenitor with an oxygen-neon-magnesium core that collapses and explodes as electron-capture supernova (ECSN). The onset of the explosion in the LMCCSN models is modelled self-consistently using the Vertex-Prometheus code, whereas the ECSN explosion is modelled using parametric neutrino transport in the Prometheus-HOTB code, choosing different explosion energies in the range of previous self-consistent models. The sAGB and LMCCSN progenitors that share structural similarities have almost spherical explosions with little metal mixing into the hydrogen envelope. A LMCCSN with less 2nd dredge-up results in a highly asymmetric explosion. It shows efficient mixing and dramatic shock deceleration in the extended hydrogen envelope. Both properties allow fast nickel plumes to catch up with the shock, leading to extreme shock deformation and aspherical shock breakout. Fallback masses of <~5x10^{-3} Msun have no significant effects on the neutron star (NS) masses and kicks. The anisotropic fallback carries considerable angular momentum, however, and determines the spin of the newly-born NS. The LMCCSNe model with less 2nd dredge-up results in a hydrodynamic and neutrino-induced NS kick of >40 km/s and a NS spin period of ~30 ms, both not largely different from those of the Crab pulsar at birth.

In our quest to identify the progenitors of Type Ia supernovae (SNe Ia), we first update the nucleosynthesis yields both for near-Chandrasekhar (Ch) and sub-Ch mass white dwarfs (WDs), for a wide range of metallicity, with our two-dimensional hydrodynamical code and the latest nuclear reaction rates. We then include the yields in our galactic chemical evolution code to predict the evolution of elemental abundances in the solar neighborhood and dwarf spheroidal (dSph) galaxies: Fornax, Sculptor, Sextans, and Carina. In the observations of the solar neighborhood stars, Mn shows an opposite trend to $\alpha$ elements, showing an increase toward higher metallicities, which is very well reproduced by deflagration-detonation transition of Ch-mass WDs, but never by double detonations of sub-Ch-mass WDs alone. The problem of Ch-mass SNe Ia was the Ni over-production at high metallicities. However, we found that Ni yields of Ch-mass SNe Ia are much lower with the solar-scaled initial composition than in previous works, which keeps the predicted Ni abundance within the observational scatter. From the evolutionary trends of elemental abundances in the solar neighborhood, we conclude that the contribution of sub-Ch-mass SNe Ia in chemical enrichment is up to 25\%. In dSph galaxies, however, larger enrichment from sub-Ch-mass SNe Ia than in the solar neighborhood may be required, which causes a decrease in [(Mg, Cr, Mn, Ni)/Fe] at lower metallicities. The observed high [Mn/Fe] ratios in Sculptor and Carina may also require additional enrichment from pure deflagrations, possibly as Type Iax supernovae. Future observations of dSph stars will provide more stringent constraints on the progenitor systems and explosion mechanism of SNe Ia.

All ten LIGO/Virgo binary black hole (BH-BH) coalescences reported from the O1/O2 runs have near zero effective spins. There are only three potential explanations of this fact. If the BH spin magnitudes are large then (i) either both BH spin vectors must be nearly in the orbital plane or (ii) the spin angular momenta of the BHs must be oppositely directed and similar in magnitude. Or, (iii) the BH spin magnitudes are small. We test the third hypothesis within the framework of the classical isolated binary evolution scenario of the BH-BH merger formation. We test three models of angular momentum transport in massive stars: a mildly efficient transport by meridional currents (as employed in the Geneva code), an efficient transport by the Tayler-Spruit magnetic dynamo (as implemented in the MESA code), and a very-efficient transport (as proposed by Fuller et al.) to calculate natal BH spins. We allow for binary evolution to increase the BH spins through accretion and account for the potential spin-up of stars through tidal interactions. Additionally, we update the calculations of the stellar-origin BH masses, include revisions to the history of star formation and to the chemical evolution across cosmic time. We find that we can match simultaneously the observed BH-BH merger rate density, BH masses, and effective spins. Models with efficient angular momentum transport are favored. The updated stellar-mass weighted gas-phase metallicity evolution now used in our models appears to be a key in better reproducing the LIGO/Virgo merger rate estimate. Mass losses during the pair-instability pulsation supernova phase are likely overestimated if the merger GW170729 hosts a BH more massive than 50 Msun. We also estimate rate of BH-NS mergers from recent LIGO/Virgo observations. Our updated models of BH-BH, BH-NS and NS-NS mergers are now publicly available at www.syntheticuniverse.org.

Recently observed pulsars with masses $\sim 1.1 ~M_{\odot}$ challenge the conventional neutron star (NS) formation path by core-collapse supernova (CCSN). Using spherically symmetric hydrodynamics simulations, we follow the collapse of a massive white dwarf (WD) core triggered by electron capture, until the formation of a proto-NS (PNS). For initial WD models with the same central density, we study the effects of a static, compact dark matter (DM) admixed core on the collapse and bounce dynamics and mass of the PNS, with DM mass $\sim 0.01 ~M_{\odot}$. We show that increasing the admixed DM mass generally leads to slower collapse and smaller PNS mass, down to about 1.0 $M_{\odot}$. Our results suggest that the accretion-induced collapse of dark matter admixed white dwarfs can produce low-mass neutron stars, such as the observed low-mass pulsar J0453+1559, which cannot be obtained by conventional NS formation path by CCSN.

We calculate the evolution of massive stars, which undergo pulsational pair-instability (PPI) when the O-rich core is formed. The evolution from the main-sequence through the onset of PPI is calculated for stars with the initial masses of $80 - 140$ $M_{\odot}$ and metallicities of $Z = 10^{-3} - 1.0$ $Z_\odot$. Because of mass loss, $Z \leq 0.5$ $Z_\odot$ is necessary for stars to form He cores massive enough (i.e., mass $>40 ~M_\odot$) to undergo PPI. The hydrodynamical phase of evolution from PPI through the beginning of Fe core collapse is calculated for the He cores with masses of $40 - 62 ~M_\odot$ and $Z = 0$. During PPI, electron-positron pair production causes a rapid contraction of the O-rich core which triggers explosive O-burning and a pulsation of the core. We study the mass dependence of the pulsation dynamics, thermodynamics, and nucleosynthesis. The pulsations are stronger for more massive He cores and result in such a large amount of mass ejection such as $3 - 13$ $M_\odot$ for $40 - 62 ~M_\odot$ He cores. These He cores eventually undergo Fe-core collapse. The $64 ~M_\odot$ He core undergoes complete disruption and becomes a pair-instability supernova. The H-free circumstellar matter ejected around these He cores is massive enough for to explain the observed light curve of Type I (H-free) superluminous supernovae with circumstellar interaction. We also note that the mass ejection sets the maximum mass of black holes (BHs) to be $\sim 50$ $M_{\odot}$, which is consistent with the masses of BHs recently detected by VIRGO and aLIGO.

We perform two-dimensional hydrodynamic simulations for the thermonuclear explosion of Chandrasekhar-mass white dwarfs with dark matter (DM) cores in Newtonian gravity. We include a 19-isotope nuclear reaction network and make use of the pure turbulent deflagration model as the explosion mechanism in our simulations. Our numerical results show that the general properties of the explosion depend quite sensitively on the mass of the DM core M$_{{\rm DM}}$: a larger M$_{{\rm DM}}$ generally leads to a weaker explosion and a lower mass of synthesized iron-peaked elements. In particular, the total mass of $^{56}$Ni produced can drop from about 0.3 to 0.03 $M_{\odot}$ as M$_{{\rm DM}}$ increases from 0.01 to 0.03 $M_{\odot}$. We have also constructed the bolometric light curves obtained from our simulations and found that our results match well with the observational data of sub-luminous Type-Ia supernovae.

Stars with $\sim 8-10~{M}_{\odot}$ evolve to form a strongly degenerate ONeMg core. When the core mass becomes close to the Chandrasekhar mass, the core undergoes electron captures on $^{24}$Mg and $^{20}$Ne, which induce the electron-capture supernova (ECSN). In order to clarify whether the ECSN leads to a collapse or thermonuclear explosion, we calculate the evolution of an $8.4~M_\odot$ star from the main sequence until the oxygen ignition in the ONeMg core. We apply the latest electron-capture rate on $^{20}$Ne including the second forbidden transition, and investigate how the location of the oxygen ignition (center or off-center) and the $Y_e$ distribution depend on the input physics and the treatment of the semiconvection and convection. The central density when the oxygen deflagration is initiated, $\rho_{\rm c,def}$, can be significantly higher than that of the oxygen ignition thanks to the convection, and we estimate $\log_{10}(\rho_{\rm c,def}/\mathrm{g~cm^{-3}})>10.10$. We perform two-dimensional simulations of the flame propagation to examine how the final fate of the ONeMg core depends on the $Y_e$ distribution and $\rho_{\rm c,def}$. We find that the deflagration starting from $\log_{10}(\rho_{\rm c,def}/\mathrm{g~cm^{-3}})>10.01 (< 10.01)$ leads to a collapse (thermonuclear explosion). Since our estimate of $\rho_{\rm c,def}$ exceeds this critical value, the ONeMg core is likely to collapse, although further studies of the convection and semiconvection before the deflagration are important.

Electron capture on $^{20}$Ne is critically important for the final stage of evolution of stars with the initial masses of 8 - 10 $M_{\odot}$. In the present paper, we evaluate electron capture rates for a forbidden transition $^{20}$Ne (0$_{g.s.}^{+}$) $\rightarrow$ $^{20}$F (2$_{g.s.}^{+}$) in stellar environments by the multipole expansion method with the use of shell-model Hamiltonians. These rates have not been accurately determined in theory as well as in experiments. Our newly evaluated rates are compared with those obtained by a prescription that treats the transition as an allowed Gamow-Teller (GT) transition with the strength determined from a recent $\beta$-decay experiment for $^{20}$F (2$_{g.s.}^{+}$) $\rightarrow$ $^{20}$Ne (0$_{g.s.}^{+}$) \citep{Kirsebom}. We find that different electron energy dependence of the transition strengths between the two methods leads to sizable differences in the weak rates of the two methods. We also find that the Coulomb effects, that is, the effects of screening on ions and electrons are non-negligible. We apply our e-capture rates on $^{20}$Ne to the calculation of the evolution of high-density O-Ne-Mg cores of 8 - 10 $M_{\odot}$ stars. We find that our new rates affect the abundance distribution and the central density at the final stage of evolution.

We present axisymmetric hydrodynamical simulations of accretion-induced collapse (AIC) of dark matter (DM) admixed rotating white dwarfs (WD) and their burst gravitational-wave (GW) signals. For initial WD models with the same central baryon density, the admixed DM is found to delay the plunge and bounce phases of AIC, and decrease the central density and mass of the proto-neutron star (PNS) produced. The bounce time, central density and PNS mass generally depend on two parameters, the admixed DM mass $M_\mathrm{DM}$ and the ratio between the rotational kinetic and gravitational energies of the inner core at bounce $\beta_\mathrm{ic,b}$. The emitted GWs have generic waveform shapes and the variation of their amplitudes $h_+$ show a degeneracy on $\beta_\mathrm{ic,b}$ and $M_\mathrm{DM}$. We found that the ratios between the GW amplitude peaks around bounce allow breaking the degeneracy and extraction of both $\beta_\mathrm{ic,b}$ and $M_\mathrm{DM}$. Even within the uncertainties of nuclear matter equation of state, a DM core can be inferred if its mass is greater than 0.03 $M_{\odot}$. We also discuss possible DM effects on the GW signals emitted by PNS g-mode oscillations. GW may boost the possibility for the detection of AIC, as well as open a new window in the indirect detection of DM.

We present two-dimensional hydrodynamics simulations of near-Chandrasekhar mass white dwarf (WD) models for Type Ia supernovae (SNe Ia) using the turbulent deflagration model with deflagration-detonation transition (DDT). We perform a parameter survey for 41 models to study the effects of the initial central density (i.e., WD mass), metallicity, flame shape, DDT criteria, and turbulent flame formula for a much wider parameter space than earlier studies. The final isotopic abundances of $^{11}$C to $^{91}$Tc in these simulations are obtained by post-process nucleosynthesis calculations. The survey includes SNe Ia models with the central density from $5 \times 10^8$ g cm$^{-3}$ to $5 \times 10^9$ g cm$^{-3}$ (WD masses of 1.30 - 1.38 $M_\odot$), metallicity from 0 to 5 $Z_{\odot}$, C/O mass ratio from 0.3 - 1.0 and ignition kernels including centered and off-centered ignition kernels. We present the yield tables of stable isotopes from $^{12}$C to $^{70}$Zn as well as the major radioactive isotopes for 33 models. Observational abundances of $^{55}$Mn, $^{56}$Fe, $^{57}$Fe and $^{58}$Ni obtained from the solar composition, well-observed SNe Ia and SN Ia remnants are used to constrain the explosion models and the supernova progenitor. The connection between the pure turbulent deflagration model and the subluminous SNe Iax is discussed. We find that dependencies of the nucleosynthesis yields on the metallicity and the central density (WD mass) are large. To fit these observational abundances and also for the application of galactic chemical evolution modeling, these dependencies on the metallicity and WD mass should be taken into account.

(${Abridged~version}$) We review how the single degenerate models for Type Ia supernovae (SNe Ia) works. In the binary star system of a white dwarf (WD) and its non-degenerate companion star, the WD accretes either H-rich matter or He and undergoes H and He shell-burning. We summarize how such shell-burning depend on the accretion rate and the WD mass and how the WD blows strong wind. We identify the following evolutionary routes for the accreting WD to trigger a thermonuclear explosion. With a decreasing rate with time, 5 phases are observed: 1. The WD increases its mass by stable H burning and blows a strong wind, which strips a part of the companion star's envelope to control the accretion rate and forms circumstellar matter (CSM). The WD may explode within CSM as an "SN Ia-CSM". 2. The wind stops and an SN Ia is triggered under steady-stable H shell-burning, which is observed as a super-soft X-ray source: "SN Ia-SSXS". 3. H shell-burning becomes unstable and many flashes recur. The WD undergoes recurrent nova (RN) whose mass ejection is smaller than the accreted matter, evolving to an "SN Ia-RN". 4. If the companion is a He star or WD, the accretion of He can trigger He and C double detonations at the sub-Chandrasekhar (sub-Chand) mass or the WD grows to the Chand mass while producing and He-wind: "SN Ia-He CSM". 5. If the accreting WD rotates rapidly, delayed trigger of an SN Ia beyond standard Chand mass can happen. The (super-Chand) WD contracts to become a delayed SN Ia after losing angular momentum. The companion star has become a He WD and CSM has disappeared: "SN Ia-He WD". We update nucleosynthesis yields of the carbon deflagration model W7, delayed detonation model WDD2, and the sub-Chand mass model to constrain the yields (such as Mn) by comparing the observations. We note the important metallicity effects on $^{58}$Ni and $^{55}$Mn.

High-resolution spectroscopy of the core of the Perseus Cluster of galaxies, using the $Hitomi$ satellite above 2 keV and the $XMM$-$Newton$ Reflection Grating Spectrometer at lower energies, provides reliable constraints on the abundances of O, Ne, Mg, Si, S, Ar, Ca, Cr, Mn, Fe, and Ni. Accounting for all known systematic uncertainties, the Ar/Fe, Ca/Fe, and Ni/Fe ratios are determined with a remarkable precision of less than 10%, while the constraints on Si/Fe, S/Fe, and Cr/Fe are at the 15% level, and Mn/Fe is measured with a 20% uncertainty. The average biases in determining the chemical composition using archival CCD spectra from $XMM$-$Newton$ and $Suzaku$ range typically from 15-40%. A simple model in which the enrichment pattern in the Perseus Cluster core and the proto-solar nebula are identical gives a surprisingly good description of the high-resolution X-ray spectroscopy results, with $\chi^2=10.7$ for 10 d.o.f. However, this pattern is challenging to reproduce with linear combinations of existing supernova nucleosynthesis calculations, particularly given the precise measurements of intermediate $\alpha$-elements enabled by $Hitomi$. We discuss in detail the degeneracies between various supernova progenitor models and explosion mechanisms, and the remaining uncertainties in these theoretical models. We suggest that including neutrino physics in the core-collapse supernova yield calculations may improve the agreement with the observed pattern of $\alpha$-elements in the Perseus Cluster core. Our results provide a complementary benchmark for testing future nucleosynthesis calculations required to understand the origin of chemical elements.

Observations of type Ia supernovae include information about the characteristic nucleosynthesis associated with these thermonuclear explosions. We consider observational constraints from iron-group elemental and isotopic ratios, to compare with various models obtained with the most-realistic recent treatment of electron captures. The nucleosynthesis is sensitive to the highest white-dwarf central densities. Hence, nucleosynthesis yields can distinguish high-density Chandrasekhar-mass models from lower-density burning models such as white-dwarf mergers. We discuss new results of post-processing nucleosynthesis for two spherical models (deflagration and/or delayed detonation models) based upon new electron capture rates. We also consider cylindrical and 3D explosion models (including deflagration, delayed-detonation, or a violent merger model). Although there are uncertainties in the observational constraints, we identify some trends in observations and the models. We make a new comparison of the models with elemental and isotopic ratios from five observed supernovae and three supernova remnants. We find that the models and data tend to fall into two groups. In one group low-density cores such as in a 3D merger model are slightly more consistent with the nucleosynthesis data, while the other group is slightly better identified with higher-density cores such as in single-degenerate 1D-3D deflagration models. Hence, we postulate that both types of environments appear to contribute nearly equally to observed SNIa. We also note that observational constraints on the yields of $^{54}$Cr and $^{54}$Fe, if available, might be used as a means to clarify the degree of geometrical symmetry of SNIa explosions.

A two-dimensional hydrodynamics code for Type Ia supernovae (SNIa) simulations is presented. The code includes a fifth-order shock-capturing scheme WENO, detailed nuclear reaction network, flame-capturing scheme and sub-grid turbulence. For post-processing we have developed a tracer particle scheme to record the thermodynamical history of the fluid elements. We also present a one-dimensional radiative transfer code for computing observational signals. The code solves the Lagrangian hydrodynamics and moment-integrated radiative transfer equations. A local ionization scheme and composition dependent opacity are included. Various verification tests are presented, including standard benchmark tests in one and two dimensions. SNIa models using the pure turbulent deflagration model and the delayed-detonation transition model are studied. The results are consistent with those in the literature. We compute the detailed chemical evolution using the tracer particles' histories, and we construct corresponding bolometric light curves from the hydrodynamics results. We also use a Graphics Processing Unit (GPU) to speed up the computation of some highly repetitive subroutines. We achieve an acceleration of 50 times for some subroutines and a factor of 6 in the global run time.

We study the equilibrium structures of white dwarfs with dark matter cores formed by non-self-annihilating dark matter DM particles with mass ranging from 1 GeV to 100 GeV, which are assumed to form an ideal degenerate Fermi gas inside the stars. For DM particles of mass 10 GeV and 100 GeV, we find that stable stellar models exist only if the mass of the DM core inside the star is less than O(10^-3) Msun and O(10^-6) Msun, respectively. The global properties of these stars, and in particular the corresponding Chandrasekhar mass limits, are essentially the same as those of traditional white dwarf models without DM. Nevertheless, in the 10 GeV case, the gravitational attraction of the DM core is strong enough to squeeze the normal matter in the core region to densities above neutron drip, far above those in traditional white dwarfs. For DM with particle mass 1 GeV, the DM core inside the star can be as massive as around 0.1 Msun and affects the global structure of the star significantly. In this case, the radius of a stellar model with DM can be about two times smaller than that of a traditional white dwarf. Furthermore, the Chandrasekhar mass limit can also be decreased by as much as 40%. Our results may have implications on to what extent type Ia supernovae can be regarded as standard candles - a key assumption in the discovery of dark energy.

In [Leung et al., Phys. Rev. D 84, 107301 (2011)], we presented our results on using a general relativistic two-fluid formalism to study the hydrostatic equilibrium configuration of an admixture of degenerate dark matter and normal nuclear matter. In this work, we present more analysis to complement our previous findings. We study the radial oscillation modes of these compact stars in detail. We find that these stars in general have two classes of oscillation modes. For a given total mass of the star, the first class of modes is insensitive to the dark-matter particle mass. They also reduce properly to the oscillation modes of the corresponding ordinary neutron star, with the same total mass, when the mass fraction of dark matter tends to zero. On the other hand, the second class of modes is due mainly to the dark-matter fluid. In the small dark-matter mass fraction limit, these modes are characterized purely by the oscillations of dark matter, while the normal matter is essentially at rest. In the intermediate regime where the mass fractions of the two fluids are comparable, the normal matter oscillates with the dark matter due to their coupling through gravity. In contrast to the first class of modes, the frequencies of these dark-matter dominated modes depend sensitively on the mass of dark-matter particles.

We study the hydrostatic equilibrium configuration of an admixture of degenerate dark matter and normal nuclear matter by using a general relativistic two-fluid formalism. We consider non- self-annihilating dark matter particles of mass ~ 1 GeV. The mass-radius relations and moments of inertia of these dark-matter admixed neutron stars are investigated and the stability of these stars is demonstrated by performing a radial perturbation analysis. We find a new class of compact stars which consists of a small normal matter core with radius of a few km embedded in a ten-kilometer- sized dark matter halo. These stellar objects may be observed as extraordinarily small neutron stars that are incompatible with realistic nuclear matter models.