Neutron Star Formation Research Articles

TWIN BINARIES: STUDIES OF STABILITY, MASS TRANSFER, AND COALESCENCE

Motivated by suggestions that binaries with almost equal-mass components ("twins") play an important role in the formation of double neutron stars and may be rather abundant among binaries, we study the stability of synchronized close and contact binaries with identical components in circular orbits. In particular, we investigate the dependency of the innermost stable circular orbit on the core mass, and we study the coalescence of the binary that occurs at smaller separations. For twin binaries composed of convective main-sequence stars, subgiants, or giants with low mass cores (M_c <~0.15M, where M is the mass of a component), a secular instability is reached during the contact phase, accompanied by a dynamical mass transfer instability at the same or at a slightly smaller orbital separation. Binaries that come inside this instability limit transfer mass gradually from one component to the other and then coalesce quickly as mass is lost through the outer Lagrangian points. For twin giant binaries with moderate to massive cores (M_c >~0.15M), we find that stable contact configurations exist at all separations down to the Roche limit, when mass shedding through the outer Lagrangian points triggers a coalescence of the envelopes and leaves the cores orbiting in a central tight binary. In addition to the formation of binary neutron stars, we also discuss the implications of our results for the production of planetary nebulae with double degenerate central binaries.

Formation of millisecond pulsars with CO white dwarf companions - I. PSR J1614−2230: evidence for a neutron star born massive

The recent discovery of a 2 M_sun binary millisecond pulsar (Demorest et al. 2010) has not only important consequences for the equation-of-state of nuclear matter at high densities but also raises the interesting question if the neutron star PSR J1614-2230 was born massive. The answer is vital for understanding neutron star formation in core collapse supernovae. Furthermore, this system raises interesting issues about the nature of the progenitor binary and how it evolved during its mass exchanging X-ray phase. In this paper we discuss the progenitor evolution of PSR J1614-2230. We have performed detailed stellar evolution modelling of intermediate-mass X-ray binaries undergoing Case A Roche-lobe overflow (RLO) and applied an analytic parameterization for calculating the outcome of either a common envelope evolution or the highly super-Eddington isotropic re-emission mode. We find two viable possibilities for the formation of the PSR J1614-2230 system: either it contained a 2.2-2.6 M_sun giant donor star and evolved through a common envelope and spiral-in phase or, more likely, it descended from a close binary system with a 4.0-5.0 M_sun main sequence donor star via Case A RLO. We conclude that the neutron star must have been born with a mass of ~1.95 M_sun or 1.7+-0.15 M_sun, respectively - which significantly exceeds neutron star birth masses in previously discovered radio pulsar systems. Based on the expected neutron star birth masses from considerations of stellar evolution and explosion models, we find that the progenitor star of PSR J1614-2230 is likely to have been more massive than 20 M_sun.

Gravitationless black holes

A gravitationless black hole model is proposed in accord with a five-dimensional fully covariant Kaluza-Klein (K-K) theory with a scalar field, which unifies the four-dimensional Einsteinian general theory of relativity and Maxwellian electromagnetic theory. It is shown that a dense compact core of a star, when it collapses to a critical density, suddenly turns off or shields its gravitational field. The core, if its mass exceeds an upper limit, directly collapses into a black hole. Otherwise, the extremely large pressure, as the gravity is turned off, immediately stops the collapse and drives the mantle material of supernova moving outward, which leads to an impulsive explosion and forms a neutron star as a remnant. A neutron star can further evolve into a black hole when it accretes enough matter from a companion star such that the total mass exceeds a lower limit. The black hole in the K-K theory is gravitationless at the surface because the scalar field is infinitely strong, which varies the equivalent gravitational constant to zero. In general, a star, at the end of its evolution, is relatively harder to collapse into a gravitationless K-K black hole than a strong gravitational Schwarzschild black hole. This is consistent with the observation of some very massive stars to form neutron stars rather than expected black holes. In addition, the gravitationless K-K black hole should be easier to generate jets than a Schwarzschild black hole.

Neutron Star Formation Could Awaken the Vacuum

Neutron Star Formation Could Awaken the Vacuum

CONSTRAINTS ON NATAL KICKS IN GALACTIC DOUBLE NEUTRON STAR SYSTEMS

Since the discovery of the first double neutron star (DNS) system in 1975 by Hulse and Taylor, there are currently 8 confirmed DNS in our galaxy. For every system, the masses of both neutron stars, the orbital semi- major axis and eccentricity are measured, and proper motion is known for half of the systems. Using the orbital parameters and kinematic information, if available, as constraints for all system, we investigate the immediate progenitor mass of the second-born neutron star and the magnitude of the supernova kick it received at birth, with the primary goal to understand the core collapse mechanism leading to neutron star formation. Compared to earlier studies, we use a novel method to address the uncertainty related to the unknown radial velocity of the observed systems. For PSR B1534+12 and PSR B1913+16, the kick magnitudes are 150 - 270 km/s and 190 - 450 km/s (with 95% confidence) respectively, and the progenitor masses of the 2nd born neutron stars are 1.3 - 3.4 Msun and 1.4 - 5.0 Msun (95%), respectively. These suggest that the 2nd born neutron star was formed by an iron core collapse supernova in both systems. For PSR J0737-3039, on the other hand, the kick magnitude is only 5 - 120 km/s (95%), and the progenitor mass of the 2nd born neutron star is 1.3 - 1.9 Msun (95%). Because of the relatively low progenitor mass and kick magnitude, the formation of the 2nd born neutron star in PSR J0737-3039 is potentially connected to an electron capture supernova of a massive O - Ne - Mg white dwarf. For the remaining 5 Galactic DNS, the kick magnitude ranges from several tens to several hundreds of km/s, and the progenitor mass of the 2nd formed neutron star can be as low as ~1.5 Msun, or as high as ~8 Msun. Therefore in these systems, it is not clear which type of supernova is more likely to form the 2nd neutron star.

Neutrino Signal of Electron-Capture Supernovae from Core Collapse to Cooling

An 8.8M{⊙} electron-capture supernova was simulated in spherical symmetry consistently from collapse through explosion to essentially complete deleptonization of the forming neutron star. The evolution time (∼9 s) is short because high-density effects suppress our neutrino opacities. After a short phase of accretion-enhanced luminosities (∼200 ms), luminosity equipartition among all species becomes almost perfect and the spectra of ν{e} and ν{μ,τ} very similar, ruling out the neutrino-driven wind as r-process site. We also discuss consequences for neutrino flavor oscillations.

The long and the short of it: modelling double neutron star and collapsar Galactic dynamics

The work presented here examines populations of double compact binary systems and tidally enhanced collapsars. We make use of BINPOP and BINKIN, two components of a recently developed population synthesis package. Results focus on correlations of both binary and spatial evolutionary population characteristics. Pulsar and long duration gamma-ray burst observations are used in concert with our models to draw the conclusions that: double neutron star binaries can merge rapidly on timescales of a few million years (much less than that found for the observed double neutron star population), common envelope evolution within these models is a very important phase in double neutron star formation, and observations of long gamma-ray burst projected distances are more centrally concentrated than our simulated coalescing double neutron star and collapsar Galactic populations. Better agreement is found with dwarf galaxy models although the outcome is strongly linked to the assumed birth radial distribution. The birth rate of the double neutron star population in our models range from 4-160 Myr^-1 and the merger rate ranges from 3-150 Myr^-1. The upper and lower limits of the rates results from including electron capture supernova kicks to neutron stars and decreasing the common envelope efficiency respectively. Our double black hole merger rates suggest that black holes should receive an asymmetric kick at birth.

Neutron star dynamics and gravitational waves

As several large scale interferometers are beginning to take data at sensitivities where astrophysical sources are predicted, the direct detection of gravitational waves may well be imminent. This would open the gravitational-wave window to our Universe and should lead to a much improved understanding of the most violent processes imaginable; the formation of black holes and neutron stars following core collapse supernovae and the merger of compact objects at the end of binary inspiral. Among the best candidate sources for gravitational waves are the oscillations, but mainly the rotational instabilities of neutron stars which can emit quite strong gravitational wave signals via which one may reveal the details of their structure. Magnetars also are neutron stars with ultra strong magnetic field whose periodic flaring activity is associated with starquakes. They are also a potential source of gravitational waves while even the quasi-periodic oscillations (QPOs) observed in the electromagnetic spectrum can reveal the details of their structure.

High space velocities of single radio pulsars versus low orbital eccentricities and masses of double neutron stars: Evidence for two different neutron star formation mechanisms

Radio pulsars tend to be high-velocity objects, which implies that the majority of them received a velocity kick of several hundreds of km s −1 at birth. However, six of the eight known double neutron stars in the galactic disk have quite low orbital eccentricities (0.085–0.27), indicating – taking into account the eccentricity induced by the mass-loss effects in the second supernova – that their second-born neutron stars received hardly any velocity kick at birth. The second-born neutron stars in these systems tend to have low masses (1.25 ± 0.07 M ⊙) and the same is true for the neutron star in the close white dwarf-neutron star system PSRJ 1145-6545, which also is the second-born star in a system of low eccentricity ( e = 0.17). Such a low mass would fit well with the formation of the second-born neutron stars by the electron-capture collapse of a Chandrasekhar-mass degenerate O–Ne–Mg core of stars of initial main-sequence mass between 8 and about 12 M ⊙. Podsiadlowski et al. (2004) have pointed out that neutron-star formation in this mass range is possible only if the star is in a close binary, following loss of the hydrogen envelope, while single stars can leave neutron stars only if they started out with masses above about 12 M ⊙; in the latter case a neutron star forms by the final collapse of an iron core. We point out here that the low masses of the second-born neutron stars in the low-eccentricity double neutron stars, together with the high space velocities of the single radio pulsars find a consistent explanation if one postulates, following ( Van den Heuvel, 2004; Podsiadlowski et al., 2004; Dewi et al., 2005) that neutron stars formed by electron-capture collapse of a degenerate O–Ne–Mg core receive hardly any kick velocity at birth, whereas neutron stars formed by iron-core collapse receive a large space velocity at birth. The existence of a population of low-kick neutron stars in binaries has important implications for understanding the formation of the neutron star population of globular clusters and for the origin of low-mass X-ray binaries and their descendants, the wide circular-orbit binary radio pulsars.

A neutron star with a carbon atmosphere in the Cassiopeia A supernova remnant

The surface of hot neutron stars is covered by a thin atmosphere. If there is accretion after neutron-star formation, the atmosphere could be composed of light elements (H or He); if no accretion takes place or if thermonuclear reactions occur after accretion, heavy elements (for example, Fe) are expected. Despite detailed searches, observations have been unable to confirm the atmospheric composition of isolated neutron stars. Here we report an analysis of archival observations of the compact X-ray source in the centre of the Cassiopeia A supernova remnant. We show that a carbon atmosphere neutron star (with low magnetic field) produces a good fit to the spectrum. Our emission model, in contrast with others, implies an emission size consistent with theoretical predictions for the radius of neutron stars. This result suggests that there is nuclear burning in the surface layers and also identifies the compact source as a very young ( approximately 330-year-old) neutron star.

Effects of mirror dark matter on neutron stars

If dark matter is made of mirror baryons, they presumably are present in all gravitationally bound structures. Here, we investigate some effects of mirror dark matter on neutron stars and discuss possible observational consequences. The general-relativistic hydrostatic equations are generalized to spherical objects with multiple fluids that interact by gravity. We use the minimal parity-symmetric extension of the standard model, which implies that the microphysics is the same in the two sectors. We find that the mass–radius relation is significantly modified in the presence of a few percent mirror baryons. This effect mimics that of other exotica, e.g., quark matter. In contrast to the common view that the neutron-star equilibrium sequence is unique, we show that it depends on the relative number of mirror baryons to ordinary baryons. It is therefore, history dependent. The critical mass for core collapse, i.e., the process by which neutron stars are created, is modified in the presence of mirror baryons. We calculate the modified Chandrasekhar mass and fit it with a polynomial. A few percent mirror baryons is sufficient to lower the critical mass for core collapse by ∼ 0.1 M ⊙ . This could allow for the formation of extraordinary compact neutron stars with low mass.

HYPERACCRETING NEUTRON STAR DISKS AND NEUTRINO ANNIHILATION

Newborn neutron stars surrounded by hyperaccreting and neutrino-cooled disks may exist in some gamma-ray bursts and/or supernovae. In this paper, we further study the structure of such a neutron star disk based on the two-region (i.e., inner and outer) disk scenario following our previous work, and calculate the neutrino annihilation luminosity from the disk in various cases. We investigate the effects of the viscosity parameter {alpha}, energy parameter {epsilon} (measuring the neutrino cooling efficiency of the inner disk), and outflow strength on the structure of the entire disk as well as the effect of emission from the neutron star surface boundary emission on the total neutrino annihilation rate. The inner disk satisfies the entropy-conservation self-similar structure for the energy parameter {epsilon} {approx_equal} 1 and the advection-dominated structure for {epsilon} < 1. An outflow from the disk decreases the density and pressure but increases the thickness of the disk. Moreover, compared with the black hole disk, the neutrino annihilation luminosity above the neutron star disk is higher, and the neutrino emission from the boundary layer could increase the neutrino annihilation luminosity by about one order of magnitude higher than the disk without boundary emission. The neutron star disk with the advection-dominatedmore » inner disk could produce the highest neutrino luminosity while the disk with an outflow has the lowest. Although a heavily mass-loaded outflow from the neutron star surface at early times of neutron star formation prevents the outflow material from being accelerated to a high bulk Lorentz factor, an energetic ultrarelativistic jet via neutrino annihilation can be produced above the stellar polar region at late times if the disk accretion rate and the neutrino emission luminosity from the surface boundary layer are sufficiently high.« less

Nucleosynthesis of the Elements in Faint Supernovae and Hypernovae

AbstractWe review the properties of supernovae (SNe) as a function of the progenitor's massM. (1)Mup- 10M⊙stars are super-AGB stars and resultant electron capture SNe may be Faint supernovae like Type IIn SN 2008S. (2) 10 - 12M⊙stars undergo Fe-core collapse to form neutron stars (NSs) and Faint supernovae. (3) 12M⊙-MBNstars undergo Fe-core collapse to form NSs and normal core-collapse supernovae. (4)MBN- 90M⊙stars undergo Fe-core collapse to form Black Holes. Resultant supernovae are bifurcate into Hypernovae and Faint supernovae. The observed properties of SN 2008ha can be explained with this type of Faint supernovae. (5) 90 - 140M⊙stars produce Luminous SNe, like SNe 2007bi and 2006gy. (6) 140 - 300M⊙stars become pair-instability supernovae which could be Luminous supernovae (SNe 2007bi and 2006gy). (7) Very massive stars withM≳ 300M⊙undergo core-collapse to form intermediate mass black holes. Some SNe could be more Luminous supernovae (like SN 2006gy).

PROBING ELECTRON-CAPTURE SUPERNOVAE: X-RAY BINARIES IN STARBURSTS

We develop population models of high-mass X-ray binaries (HMXBs) formed after bursts of star formation and we investigate the effect of electron-capture supernovae (ECS) of massive ONeMg white dwarfs and the hypothesis that ECS events are associated with typically low supernova kicks imparted to the nascent neutron stars. We identify an interesting ECS bump in the time evolution of HMXB numbers; this bump is caused by significantly increased production of wind-fed HMXBs 20–60 Myr post-starburst. The amplitude and age extent of the ECS bump depend on the strength of ECS kicks and the mass range of ECS progenitors. We also find that ECS-HMXBs form through a specific evolutionary channel that is expected to lead to binaries with Be donors in wide orbits. These characteristics, along with their sensitivity to ECS properties, provide us with an intriguing opportunity to probe ECS physics and progenitors through studies of starbursts of different ages. Specifically, the case of the Small Magellanic Cloud, with a significant observed population of Be-HMXBs and starburst activity 30–60 Myr ago, arises as a promising laboratory for understanding the role of ECS in neutron star formation.

Deep optical observations of the central X-ray source in the Puppis A supernova remnant

X-ray observations reveiled a group of radio-silent isolated neutron stars (INSs) at the centre of young supernova remnants (SNRs), dubbed central compact objects or CCOs, with properties different from those of classical rotation-powered pulsars. In at least three cases, evidence points towards CCOs being low-magnetized INSs, born with slow rotation periods, and possibly accreting from a debris disc of material formed out of the supernova event. Understanding the origin of the diversity of the CCOs can shed light on supernova explosion and neutron star formation models. Optical/infrared (IR) observations are crucial to test different CCO interpretations. The aim of our work is to perform a deep optical investigation of the CCO RX J0822.0-4300 in the Puppis A SNR, one of the most poorly understood in the CCO family. By using as a reference the Chandra X-ray coordinates of RX J0822.0-4300, we performed deep optical observations in the B, V and I bands with the Very Large Telescope (VLT). We found no candidate optical counterpart within 3 sigma of the computed Chandra X-ray position down to 5 sigma limits of B~27.2, V~26.9, and I~25.6, the deepest obtained in the optical band for this source. These limits confirm the non-detection of a companion brighter than an M5 dwarf. At the same time, they do not constrain optical emission from the neutron star surface, while emission from the magnetosphere would require a spectral break in the optical/IR.

DELAYED NEUTRINO-DRIVEN SUPERNOVA EXPLOSIONS AIDED BY THE STANDING ACCRETION-SHOCK INSTABILITY

We present results of 2D hydrodynamic simulations of stellar core collapse, which confirm that the neutrino-heating mechanism remains viable for the explosion of a wider mass range of supernova progenitors with iron cores. We used an energy-dependent treatment of the neutrino transport based on the ray-by-ray plus approximation, in which the number, energy, and momentum equations are closed with a variable Eddington factor obtained by iteratively solving a model Boltzmann equation. We focus on the evolution of a 15 Msun progenitor and show that shock revival and the explosion are initiated at about 600 ms post bounce, powered by neutrino energy deposition. Similar to previous findings for an 11.2 Msun star, but significantly later, the onset of the explosion is fostered by the standing accretion shock instability (SASI). This instability exhibits highest growth rates for the dipole and quadrupole modes, which lead to large-amplitude bipolar shock oscillations and push the shock to larger radii, thus increasing the time accreted matter is exposed to neutrino heating in the gain layer. Therefore also convective overturn behind the shock is strengthened. A soft nuclear equation of state that causes a rapid contraction and a smaller radius of the forming neutron star and thus a fast release of gravitational binding energy, seems to be more favorable for an explosion. Rotation has the opposite effect because it leads to a more extended and cooler neutron star and thus lower neutrino luminosities and mean energies and overall less neutrino heating. Neutron star g-mode oscillations and the acoustic mechanism play no important role in our simulations. (abridged)

TYPE Ib SUPERNOVA 2008D ASSOCIATED WITH THE LUMINOUS X-RAY TRANSIENT 080109: AN ENERGETIC EXPLOSION OF A MASSIVE HELIUM STAR

We present a theoretical model for supernova (SN) 2008D associated with the luminous X-ray transient 080109. The bolometric light curve and optical spectra of the SN are modelled based on the progenitor models and the explosion models obtained from hydrodynamic/nucleosynthetic calculations. We find that SN 2008D is a more energetic explosion than normal core-collapse supernovae, with an ejecta mass of Mej = 5.3 +- 1.0 Msun and a kinetic energy of E = 6.0 +- 2.5 x 10^{51} erg. The progenitor star of the SN has a 6-8 Msun He core with essentially no H envelope (< 5 x 10^{-4} Msun) prior to the explosion. The main-sequence mass of the progenitor is estimated to be Mms =20-25 Msun, with additional systematic uncertainties due to convection, mass loss, rotation, and binary effects. These properties are intermediate between those of normal SNe and hypernovae associated with gamma-ray bursts. The mass of the central remnant is estimated as 1.6 - 1.8 Msun, which is near the boundary between neutron star and black hole formation.

SGRs and AXPs: Evidence for Delayed Amplification of Magnetic Field after Neutron Star Formation?

We present new analysis of the birth rate of AXPs and SGRS and their associated SNRs. Using Kolmogorov-Smirnov statistics together with parametric fits based on a robust estimator, we find a birth rate of ∼1/(1000 years) for AXPs/SGRs and their associated SNRs. These high rates suggest that all massive stars (greater than ∼(23–32)M⊙) give rise to remnants with magnetar-like fields. Observations indicate a limited fraction of high magnetic fields in these progenitors; thus our study is suggestive of magnetic field amplification. Dynamo mechanisms during the birth of the neutron stars require spin rates much faster than either observations or theory indicate. We propose that massive stars produce neutron stars with normal (∼1012 G) magnetic fields, which are then amplified to1014-1015 G after a delay of hundreds of years. The amplification is speculated to be a consequence of color ferromagnetism and to occur with a delay after the neutron star core reaches quark deconfinement density (i.e., the quark-nova scenario). The delayed amplification allows one to interpret simultaneously the high birth rate and high magnetic fields of AXPs/SGRs and their link to massive stars.

Chemical Evolution of the Juvenile Universe

AbstractModels of average Galactic chemical abundances are in good general agreement with observations for [Fe/H] &gt; –1.5, but there are gross discrepancies at lower metallicities. Only massive stars contribute to the chemical evolution of the ‘juvenile universe’ corresponding to [Fe/H] ≲ –1.5. If Type II supernovae (SNe II) are the only relevant sources, then the abundances in the interstellar medium of the juvenile epoch are simply the sum of different SN II contributions. Both low-mass (∼8–11 M⊙) and normal (∼12–25 M⊙) SNe II produce neutron stars, which have intense neutrino-driven winds in their nascent stages. These winds produce elements such as Sr, Y and Zr through charged-particle reactions (CPR). Such elements are often called the ‘light r-process elements’, but are considered here as products of CPR and not the r process. The observed absence of production of the low-A elements (Na through Zn including Fe) when the true r-process elements (Ba and above) are produced requires that only low-mass SNe II be the site if the r process occurs in SNe II. Normal SNe II produce the CPR elements in addition to the low-A elements. This results in a two-component model that is quantitatively successful in explaining the abundances of all elements relative to hydrogen for –3 ≲ [Fe/H] ≲ –1.5. This model explicitly predicts that [Sr/Fe] ≥ –0.32. Recent observations show that there are stars with [Sr/Fe] ≲ –2 and [Fe/H] &lt; –3. This proves that the two-component model is not correct and that a third component is necessary to explain the observations. The production of CPR elements associated with the formation of neutron stars requires that the third component must be massive stars ending as black holes. It is concluded that stars of ∼25–50 M⊙ (possibly up to ∼100 M⊙) are the appropriate candidates. These produce hypernovae (HNe) that have very high Fe yields and are observed today. Stars of ∼140–260 M⊙ are completely disrupted upon explosion. However, they produce an abundance pattern greatly deficient in elements of odd atomic numbers, which is not observed, and therefore they are not considered as a source here. Using a Salpeter initial mass function, it is shown that HNe are a source of Fe that far outweighs normal SNe II, with the former and the latter contributing ∼24% and ∼9% of the solar Fe abundance, respectively. It follows that the usual assignment of ∼⅓ of the solar Fe abundance to normal SNe II is not correct. This leads to a simple three-component model including low-mass and normal SNe II and HNe, which gives a good description of essentially all the data for stars with [Fe/H] ≲ –1.5. We conclude that HNe are more important than normal SNe II in the chemical evolution of the low-A elements from Na through Zn (including Fe), in sharp distinction to earlier models.

Core-collapse supernovae: Magnetorotational explosions and jet formation

Core-collapse supernovae are accompanied by formation of neutron stars. The gravitational energy is transformed into the energy of the explosion, observed as SN II, SN Ib,c type supernovae. We present results of 2D MHD simulations, where the source of energy is rotation and the magnetic field serves as a “transition belt” for the transformation of the rotation energy into the energy of the explosion. The toroidal part of the magnetic energy initially grows linearly with time due to differential rotation. When the twisted toroidal component strongly exceeds the poloidal field, magneto-rotational instability develops, leading to a drastic acceleration in the growth of magnetic energy. Finally, a fast MHD shock is formed, producing a supernova explosion. A mildly collimated jet is produced for the dipolelike type of the initial field.