The evolution and masses of the neutron star and donor star in the high mass X-ray binary OAO 1657−415

We report near-infrared radial velocity (RV) measurements of the recently identified donor star in the high mass X-ray binary (HMXB) system OAO 1657−415 obtained in the H band using ISAAC on the Very Large Telescope. Cross-correlation methods were employed to construct a RV curve with a semi-amplitude of 22.1 ± 3.5 km s−1. Combined with other measured parameters of this system it provides a dynamically determined neutron star (NS) mass of 1.42 ± 0.26 M⊙ and a mass of 14.3 ± 0.8 M⊙ for the Ofpe/WN9 highly evolved donor star. OAO 1657−415 is an eclipsing HMXB pulsar with the largest eccentricity and orbital period of any within its class. Of the 10 known eclipsing X-ray binary pulsars OAO 1657−415 becomes the ninth with a dynamically determined NS mass solution and only the second in an eccentric system. Furthermore, the donor star in OAO 1657−415 is much more highly evolved than the majority of the supergiant donors in other HMXBs, joining a small but growing list of HMXBs donors with extensive hydrogen depleted atmospheres. Considering the evolutionary development of OAO 1657−415, we have estimated the binding energy of the envelope of the mass donor and find that there is insufficient energy for the removal of the donor’s envelope via spiral-in, ruling out a common envelope evolutionary scenario. With its non-zero eccentricity and relatively large orbital period the identification of a definitive evolutionary pathway for OAO 1657−415 remains problematic, we conclude by proposing two scenarios which may account for OAO 1657−415 current orbital configuration.

Simulated magnetic field expulsion in neutron star cores

The study of long-term evolution of neutron star (NS) magnetic fields is key to understanding the rich diversity of NS observations, and to unifying their nature despite the different emission mechanisms and observed properties. Such studies in principle permit a deeper understanding of the most important parameters driving their apparent variety, e.g. radio pulsars, magnetars, X-ray dim isolated NSs, gamma-ray pulsars. We describe, for the first time, the results from self-consistent magnetothermal simulations considering not only the effects of the Hall-driven field dissipation in the crust, but also adding a complete set of proposed driving forces in a superconducting core. We emphasize how each of these core-field processes drive magnetic evolution and affect observables, and show that when all forces are considered together in vectorial form, the net expulsion of core magnetic flux is negligible, and will have no observable effect in the crust (consequently in the observed surface emission) on megayear time-scales. Our new simulations suggest that strong magnetic fields in NS cores (and the signatures on the NS surface) will persist long after the crustal magnetic field has evolved and decayed, due to the weak combined effects of dissipation and expulsion in the stellar core.

Dynamical analysis on heavy-ion fusion reactions near Coulomb barrier

The shell correction is proposed in the improved isospin dependent quantum molecular dynamics (Im-IQMD) model, which plays an important role in heavy-ion fusion reactions near Coulomb barrier. By using the ImIQMD model, the static and dynamical fusion barriers, dynamical barrier distribution in the fusion reactions are analyzed systematically. The fusion and capture excitation functions for a series of reaction systems are calculated and compared with experimental data. It is found that the fusion cross sections for neutron-rich systems increase obviously, and the strong shell effects of two colliding nuclei result in a decrease of the fusion cross sections at the sub-barrier energies. The lowering of the dynamical fusion barriers favors the enhancement of the sub-barrier fusion cross sections, which is related to the nucleon transfer and the neck formation in the fusion reactions.

Dynamics of pion production in heavy-ion collisions around 1A GeV energies

Within the framework of the improved isospin-dependent quantum molecular dynamics model, the dynamics of pion emission in heavy-ion collisions in the region of 1A GeV energies as a probe of nuclear symmetry energy at suprasaturation densities is investigated systematically. The total pion multiplicities and the pi(-)/pi(+) yields are calculated for selected Skyrme parameters SkP, SLy6, Ska, and SIII and also for the cases of different stiffness of symmetry energy with the parameter SLy6. The influence of Coulomb potential, symmetry energy, and in-medium pion potential on the pion production is investigated and compared to each other by analyzing the distributions of transverse momentum and longitudinal rapidity and also the excitation functions of the total pion and the pi(-)/pi(+) ratio. The directed flow, elliptic flow, and polar-angle distributions are calculated for the cases of different collision centralities and also the various stiffnesses of the symmetry energies. A comparison of the calculations with the available experimental data is performed.

The masses of the neutron and donor star in the high-mass X-ray binary IGR J18027-2016

Aims. We report near-infrared observations of the supergiant donor to the eclipsing high mass X-ray binary pulsar IGR J18027-2016. We aim to determine its spectral type and measure its radial velocity curve and hence determine the stellar masses of the components. Methods. ESO/VLT observations of the donor utilising the NIR spectrograph ISAAC were obtained in the H and K bands. The multi-epoch H band spectra were cross-correlated with RV templates in order to determine a radial solution for the system. Results. The spectral type of the donor was confirmed as B0-1 I. The radial velocity curve constructed has a semi-amplitude of 23.8 ± 3.1 km s-1. Combined with other measured system parameters, a dynamically determined neutron star mass of 1.4  ±  0.2–1.6  ±  0.3 M⊙ is found. The mass range of the B0-B1 I donor was 18.6  ±  0.8–21.8  ±  2.4 M⊙. These lower and upper limits were obtained under the assumption that the system is viewed edge-on (i = 90° with β = 0.89) for the lower limit and the donor fills its Roche lobe (β = 1 with i = 73.1°) for the upper limit respectively.

Dynamical stability of nascent neutron stars

This thesis presents a study of the dynamical stability of nascent neutron stars resulting from the accretion induced collapse of rapidly rotating white dwarfs.

Chapter 2 and part of Chapter 3 study the equilibrium models for these neutron stars. They are constructed by assuming that the neutron stars have the same masses, angular momenta, and specific angular momentum distributions as the pre-collapse white dwarfs. If the pre-collapse white dwarf is rapidly rotating, the collapsed object will contain a high density central core of size about 20 km, surrounded by a massive accretion torus extending to hundreds of kilometers from the rotation axis. The ratio of the rotational kinetic energy to gravitational binding energy, β, of these neutron stars is all found to be less than 0.27.

Chapter 3 studies the dynamical stability of these neutron stars by numerically evolving the linearized hydrodynamical equations. A dynamical bar-mode instability is observed when the β of the star is greater than the critical value β_d ≈ 0.25. It is expected that the unstable mode will persist until a substantial amount of angular momentum is carried away by gravitational radiation. The detectability of these sources is studied and it is estimated that LIGO II is unlikely to detect them unless the event rate is greater than 10^-6/year/galaxy.

All the calculations on the structure and stability of the neutron stars in Chapters 2 and 3 are carried out using Newtonian hydrodynamics and gravity. Chapter 4 studies the relativistic effects on the structure of these neutron stars. New techniques are developed and used to construct neutron star models to the first post-Newtonian (1PN) order. The structures of the 1PN models are qualitatively similar to the corresponding Newtonian models, but the values of β are somewhat smaller. The maximum β for these 1PN neutron stars is found to be 0.24, which is 8% smaller than the Newtonian result (0.26). However, relativistic effects will also change the critical value β_d. A detailed post-Newtonian stability analysis has yet to be carried out to study the relativistic effects on the dynamical stability of these neutron stars.

(PF2)-P-3 neutron superfluidity in neutron stars and three-body force effect

We investigate the (PF2)-P-3 neutron superfluidity in beta-stable neutron star matter and neutron stars by using the BCS theory and the Brueckner-Hartree-Fock approach. We adopt the Argonne V-18 potential supplemented with a microscopic three-body force as the realistic nucleon-nucleon interaction. We have concentrated on studying the three-body force effect on the (PF2)-P-3 neutron pairing gap. It is found that the three-body force effect is to enhance remarkably the (PF2)-P-3 neutron superfluidity in neutron star matter and neutron stars.

Hyperon coupling dependence of hadron matter properties in relativistic mean field model

The properties of hadronic matter at beta equilibrium in a wide range of densities are described by appropriate equations of state in the framework of the relativistic mean field model. Strange meson fields, namely the scalar meson field sigma*(975) and the vector meson field sigma*(1020), are included in the present work. We discuss and compare the results of the equation of state, nucleon effective mass, and strangeness fraction obtained by adopting the TM1, TMA, and GL parameter sets for nuclear sector and three different choices for the hyperon couplings. We find that the parameter set TM1 favours the onset of hyperons most, while at high densities the GL parameter set leads to the most hyperon-rich matter. For a certain parameter set (e.g. TM1), the most hyperon-rich matter is obtained for the hyperon potential model. The influence of the hyperon couplings on the effective mass of nucleon, is much weaker than that on the nucleon parameter set. The nonstrange mesons dominate essentially the global properties of dense hyperon matter. The hyperon potential model predicts the lowest value of the neutron star maximum mass of about 1.45 M-sun to be 0.4-0.5 M-sun lower than the prediction by using the other choices for hyperon couplings.

Hyperon-hyperon interaction in relativistic mean field model

By using the new experimental data of Lambda Lambda potential, this paper has performed a full calculation for strange hadronic matter with different strangeness contents as well as its consequences on the global properties of neutron star matter in relativistic mean field model. It finds that the new weak hyperon - hyperon interaction makes the equations of state much stiffer than the result of the previous strong hyperon-hyperon interaction, and even stiffer than the result without consideration of hyperon -hyperon interaction. This new hyperon -hyperon interaction results in a maximum mass of 1.75M(circle dot) ( where M-circle dot stands for the mass of the Sun), about 0.2-0.5M(circle dot) larger than the previous prediction with the presence of hyperons. After examining carefully the onset densities of kaon condensation it finds that this new weak version of hyperon -hyperon interaction favours the occurrence of kaons in comparison with the strong one.

The spotty donor star in the X-ray transient Cen X-4

We accurately determine the fundamental system parameters of the neutron star X-ray transient Cen X-4 solely using phase-resolved high-resolution UV-Visual Echelle Spectrograph spectroscopy. We first determine the radial-velocity curve of the secondary star and then model the shape of the phase-resolved absorption line profiles using an X-ray binary model. The model computes the exact rotationally broadened, phase-resolved spectrum and does not depend on assumptions about the rotation profile, limb-darkening coefficients and the effects of contamination from an accretion disc. We determine the secondary star-to-neutron star binary mass ratio to be 0.1755 ± 0.0025, which is an order of magnitude more accurate than previous estimates. We also constrain the inclination angle to be 32^{+8}_{-2} degrees. Combining these values with the results of the radial-velocity study gives a neutron star mass of 1.94^{+0.37}_{-0.85}M⊙ consistent with previous estimates. Finally, we perform the first Roche tomography reconstruction of the secondary star in an X-ray binary. The tomogram reveals surface inhomogeneities that are due to the presence of cool starspots. A large cool polar spot, similar to that seen in Doppler images of rapidly rotating isolated stars, is present on the Northern hemisphere of the K7 secondary star and we estimate that ~4 percent of the total surface area of the donor star is covered with spots.This evidence for starspots supports the idea that magnetic braking plays an important role in the evolution of low-mass X-ray binaries.

Astrophysical S factors for fusion reactions involving C, O, Ne, and Mg isotopes

Using the Sao Paulo potential and the barrier penetration formalism we have calculated the astrophysical factor S(E) for 946 fusion reactions involving stable and neutron-rich isotopes of C, O, Ne, and Mg for center-of-mass energies E varying from 2 to approximate to 18-30 MeV (covering the range below and above the Coulomb barrier). We have parameterized the energy dependence, S(E), by an accurate universal 9-parameter analytic expression and present tables of fit parameters for all the reactions. We also discuss the reduced 3-parameter version of our fit which is highly accurate at energies below the Coulomb barrier, and outline the procedure for calculating the reaction rates. The results can be easily converted to thermonuclear or pycnonuclear reaction rates to simulate various nuclear burning phenomena, in particular, stellar burning at high temperatures and nucleosynthesis in high density environments. (C) 2010 Elsevier Inc. All rights reserved

Evolution of magnetic field and spin period in accreting neutron stars

Based on the accretion-induced magnetic field decay model, in which a frozen field and an incompressible fluid are assumed, we obtain the following results: (1) an analytic relation between the magnetic field and spin period, if the fastness parameter of the accretion disk is neglected (The evolutionary tracks of accreting neutron stars in the P-B diagram in our model are different from the equilibrium period lines when the influence of the fastness parameter is taken into account.); (2) the theoretical minimum spin period of an accreting neutron star is max(1.1ms (DeltaM/M(circle dot))(-1)R(6)(-5/14) I(45)(M/M(circle dot))(-1/2),1.1ms (M/M(circle dot))(-1/2) R(6)(17/14)), independent of the accretion rate (X-ray luminosity) but dependent on the total accretion mass, DeltaM; however, the minimum magnetic field depends on the accretion rate; (3) the magnetic field strength decreases faster with time than does the period.

The saturated magnetic field of the neutron stars

Considering the ferromagnetic screening for the decay of the X-ray neutron star magnetic field in the binary accretion phase, the phase transition of ferromagnetic materials in the crust of neutron star induces the ferromagnetic screening saturation of the accreted crust, which results in the minimum surface magnetic field of the accreting neutron star, about 108 G, if the accreted matter has completely replaced the crust mass of the neutron star. The magnetic field evolution versus accreted mass is given as Bs ∝ ΔM-0.9, and the obtained magnetic field versus spin period relation is consistent with the distribution of the binary X-ray sources and recycled pulsars. The further thermal effect on the magnetic evolution is also studied.

Self-bound interacting QCD matter in compact stars

The quark gluon plasma (QGP) at zero temperature and high baryon number is a system that may be present inside compact stars. It is quite possible that this cold QGP shares some relevant features with the hot QGP observed in heavy ion collisions, being also a strongly interacting system. In a previous work we have derived from the QCD Lagrangian an equation of state (EOS) for the cold QGP, which can be considered an improved version of the MIT bag-model EOS. Compared to the latter, our EOS reaches higher values of the pressure at comparable baryon densities. This feature is due to perturbative corrections and also to nonperturbative effects. Here we apply this EOS to the study of neutron stars, discussing the absolute stability of quark matter and computing the mass-radius relation for self-bound (strange) stars. The maximum masses of the sequences exceed two solar masses, in agreement with the recently measured values of the mass of the pulsar PSR J1614-2230, and the corresponding radii of around 10-11 km.