Orbital phase-resolved spectroscopy of 4U 1538−52 with MAXI

Context. 4U 1538−52, an absorbed high mass X-ray binary with an orbital period of ~3.73 days, shows moderate orbital intensity modulations with a low level of counts during the eclipse. Several models have been proposed to explain the accretion at different orbital phases by a spherically symmetric stellar wind from the companion. Aims. The aim of this work is to study both the light curve and orbital phase spectroscopy of this source in the long term. In particular, we study the folded light curve and the changes in the spectral parameters with orbital phase to analyse the stellar wind of QV Nor, the mass donor of this binary system. Methods. We used all the observations made from the Gas Slit Camera on board MAXI of 4U 1538−52 covering many orbits continuously. We obtained the good interval times for all orbital phase ranges, which were the input for extracting our data. We estimated the orbital period of the system and then folded the light curves, and we fitted the X-ray spectra with the same model for every orbital phase spectrum. We also extracted the averaged spectrum of all the MAXI data available. Results. The MAXI spectra in the 2–20 keV energy range were fitted with an absorbed Comptonisation of cool photons on hot electrons. We found a strong orbital dependence of the absorption column density but neither the fluorescence iron emission line nor low energy excess were needed to fit the MAXI spectra. The variation in the spectral parameters over the binary orbit were used to examine the mode of accretion onto the neutron star in 4U 1538−52. We deduce a best value of Ṁ/v∞ = 0.65 × 10-9M⊙ yr-1/ (km s-1) for QV Nor.

Magnetars: Properties, Origin and Evolution

Magnetars are neutron stars in which a strong magnetic field is the main energy source. About two dozens of magnetars, plus several candidates, are currently known in our Galaxy and in the Magellanic Clouds. They appear as highly variable X-ray sources and, in some cases, also as radio and/or optical pulsars. Their spin periods (2–12 s) and spin-down rates (∼10−13–10−10 s s−1) indicate external dipole fields of ∼1013−15 G, and there is evidence that even stronger magnetic fields are present inside the star and in non-dipolar magnetospheric components. Here we review the observed properties of the persistent emission from magnetars, discuss the main models proposed to explain the origin of their magnetic field and present recent developments in the study of their evolution and connection with other classes of neutron stars.

Modelling of the surface emission of the low magnetic field magnetar SGR 0418+5729

We perform a detailed modelling of the post-outburst surface emission of the low magnetic field magnetar SGR 0418+5729. The dipolar magnetic field of this source, B=6×1012G estimated from its spin-down rate, is in the observed range of magnetic fields for normal pulsars. The source is further characterized by a high pulse fraction and a single-peak profile. Using synthetic temperature distribution profiles, and fully accounting for the general-relativistic effects of light deflection and gravitational redshift, we generate synthetic X-ray spectra and pulse profiles that we fit to the observations. We find that asymmetric and symmetric surface temperature distributions can reproduce equally well the observed pulse profiles and spectra of SGR 0418. None the less, the modelling allows us to place constraints on the system geometry (i.e. the angles ψ and ξ that the rotation axis makes with the line of sight and the dipolar axis, respectively), as well as on the spot size and temperature contrast on the neutron star surface. After performing an analysis iterating between the pulse profile and spectra, as done in similar previous works, we further employed, for the first time in this context, a Markov-Chain Monte Carlo approach to extract constraints on the model parameters from the pulse profiles and spectra, simultaneously. We find that, to reproduce the observed spectrum and flux modulation: (a) the angles must be restricted to 65° ≲ ψ + ξ ≲ 125° or 235° ≲ ψ + ξ ≲ 295°; (b) the temperature contrast between the poles and the equator must be at least a factor of ∼6, and (c) the size of the hottest region ranges between 0.2 and 0.7 km (including uncertainties on the source distance). Lastly, we interpret our findings within the context of internal and external heating models.

Constraining the GRB-Magnetar Model by Means of the Galactic Pulsar Population

A large fraction of Gamma-ray bursts (GRBs) displays an X-ray plateau phase within <105 s from the prompt emission, proposed to be powered by the spin-down energy of a rapidly spinning newly born magnetar. In this work we use the properties of the Galactic neutron star population to constrain the GRB-magnetar scenario. We re-analyze the X-ray plateaus of all Swift GRBs with known redshift, between 2005 January and 2014 August. From the derived initial magnetic field distribution for the possible magnetars left behind by the GRBs, we study the evolution and properties of a simulated GRB-magnetar population using numerical simulations of magnetic field evolution, coupled with Monte Carlo simulations of Pulsar Population Synthesis in our Galaxy. We find that if the GRB X-ray plateaus are powered by the rotational energy of a newly formed magnetar, the current observational properties of the Galactic magnetar population are not compatible with being formed within the GRB scenario (regardless of the GRB type or rate at z = 0). Direct consequences would be that we should allow the existence of magnetars and "super-magnetars" having different progenitors, and that Type Ib/c SNe related to Long GRBs form systematically neutron stars with higher initial magnetic fields. We put an upper limit of ≤16 "super-magnetars" formed by a GRB in our Galaxy in the past Myr (at 99% c.l.). This limit is somewhat smaller than what is roughly expected from Long GRB rates, although the very large uncertainties do not allow us to draw strong conclusion in this respect.

The outburst decay of the low magnetic field magnetar SWIFT J1822.3−1606: phase-resolved analysis and evidence for a variable cyclotron feature

We study the timing and spectral properties of the low-magnetic field, transient magnetar SWIFT J1822.3−1606 as it approached quiescence. We coherently phase-connect the observations over a time-span of ∼500 d since the discovery of SWIFT J1822.3−1606 following the Swift-Burst Alert Telescope (BAT) trigger on 2011 July 14, and carried out a detailed pulse phase spectroscopy along the outburst decay. We follow the spectral evolution of different pulse phase intervals and find a phase and energy-variable spectral feature, which we interpret as proton cyclotron resonant scattering of soft photon from currents circulating in a strong (≳1014 G) small-scale component of the magnetic field near the neutron star surface, superimposed to the much weaker (∼3 × 1013 G) magnetic field. We discuss also the implications of the pulse-resolved spectral analysis for the emission regions on the surface of the cooling magnetar.

Are pulsars born with a hidden magnetic field?

The observation of several neutron stars in the centre of supernova remnants and with significantly lower values of the dipolar magnetic field than the average radio-pulsar population has motivated a lively debate about their formation and origin, with controversial interpretations. A possible explanation requires the slow rotation of the protoneutron star at birth, which is unable to amplify its magnetic field to typical pulsar levels. An alternative possibility, the hidden magnetic field scenario, considers the accretion of the fallback of the supernova debris on to the neutron star as responsible for the submergence (or screening) of the field and its apparently low value. In this paper, we study under which conditions the magnetic field of a neutron star can be buried into the crust due to an accreting, conducting fluid. For this purpose, we consider a spherically symmetric calculation in general relativity to estimate the balance between the incoming accretion flow and the magnetosphere. Our study analyses several models with different specific entropy, composition, and neutron star masses. The main conclusion of our work is that typical magnetic fields of a few times 1012 G can be buried by accreting only 10−3–10−2 M⊙, a relatively modest amount of mass. In view of this result, the central compact object scenario should not be considered unusual, and we predict that anomalously weak magnetic fields should be common in very young (< few kyr) neutron stars.

Measuring the stellar wind parameters in IGR J17544-2619 and Vela X-1 constrains the accretion physics in supergiant fast X-ray transient and classical supergiant X-ray binaries

Context. Classical supergiant X-ray binaries (SGXBs) and supergiant fast X-ray transients (SFXTs) are two types of high-mass X-ray binaries (HMXBs) that present similar donors but, at the same time, show very different behavior in the X-rays. The reason for this dichotomy of wind-fed HMXBs is still a matter of debate. Among the several explanations that have been proposed, some of them invoke specific stellar wind properties of the donor stars. Only dedicated empiric analysis of the donors’ stellar wind can provide the required information to accomplish an adequate test of these theories. However, such analyses are scarce. Aims. To close this gap, we perform a comparative analysis of the optical companion in two important systems: IGR J17544-2619 (SFXT) and Vela X-1 (SGXB). We analyze the spectra of each star in detail and derive their stellar and wind properties. As a next step, we compare the wind parameters, giving us an excellent chance of recognizing key differences between donor winds in SFXTs and SGXBs. Methods. We use archival infrared, optical and ultraviolet observations, and analyze them with the non-local thermodynamic equilibrium (NLTE) Potsdam Wolf-Rayet model atmosphere code. We derive the physical properties of the stars and their stellar winds, accounting for the influence of X-rays on the stellar winds. Results. We find that the stellar parameters derived from the analysis generally agree well with the spectral types of the two donors: O9I (IGR J17544-2619) and B0.5Iae (Vela X-1). The distance to the sources have been revised and also agree well with the estimations already available in the literature. In IGR J17544-2619 we are able to narrow the uncertainty to d = 3.0 ± 0.2 kpc. From the stellar radius of the donor and its X-ray behavior, the eccentricity of IGR J17544-2619 is constrained to e< 0.25. The derived chemical abundances point to certain mixing during the lifetime of the donors. An important difference between the stellar winds of the two stars is their terminal velocities (ν∞ = 1500 km s-1 in IGR J17544-2619 and ν∞ = 700 km s-1 in Vela X-1), which have important consequences on the X-ray luminosity of these sources. Conclusions. The donors of IGR J17544-2619 and Vela X-1 have similar spectral types as well as similar parameters that physically characterize them and their spectra. In addition, the orbital parameters of the systems are similar too, with a nearly circular orbit and short orbital period. However, they show moderate differences in their stellar wind velocity and the spin period of their neutron star which has a strong impact on the X-ray luminosity of the sources. This specific combination of wind speed and pulsar spin favors an accretion regime with a persistently high luminosity in Vela X-1, while it favors an inhibiting accretion mechanism in IGR J17544-2619. Our study demonstrates that the relative wind velocity is critical in class determination for the HMXBs hosting a supergiant donor, given that it may shift the accretion mechanism from direct accretion to propeller regimes when combined with other parameters.

A simple approach to the supernova progenitor–explosion connection

We present a new approach to understand the landscape of supernova explosion energies, ejected nickel masses, and neutron star birth masses. In contrast to other recent parametric approaches, our model predicts the properties of neutrino-driven explosions based on the pre-collapse stellar structure without the need for hydrodynamic simulations. The model is based on physically motivated scaling laws and simple differential equations describing the shock propagation, the contraction of the neutron star, the neutrino emission, the heating conditions, and the explosion energetics. Using model parameters compatible with multi-D simulations and a fine grid of thousands of supernova progenitors, we obtain a variegated landscape of neutron star and black hole formation similar to other parametrized approaches and find good agreement with semi-empirical measures for the ‘explodability’ of massive stars. Our predicted explosion properties largely conform to observed correlations between the nickel mass and explosion energy. Accounting for the coexistence of outflows and downflows during the explosion phase, we naturally obtain a positive correlation between explosion energy and ejecta mass. These correlations are relatively robust against parameter variations, but our results suggest that there is considerable leeway in parametric models to widen or narrow the mass ranges for black hole and neutron star formation and to scale explosion energies up or down. Our model is currently limited to an all-or-nothing treatment of fallback and there remain some minor discrepancies between model predictions and observational constraints.

Progenitor-Dependent Explosion Dynamics in Self-Consistent, Axisymmetric Simulations of Neutrino-Driven Core-Collapse Supernovae

We present self-consistent, axisymmetric core-collapse supernova simulations performed with the Prometheus-Vertex code for 18 pre-supernova models in the range of 11–28 M ⊙, including progenitors recently investigated by other groups. All models develop explosions, but depending on the progenitor structure, they can be divided into two classes. With a steep density decline at the Si/Si–O interface, the arrival of this interface at the shock front leads to a sudden drop of the mass-accretion rate, triggering a rapid approach to explosion. With a more gradually decreasing accretion rate, it takes longer for the neutrino heating to overcome the accretion ram pressure and explosions set in later. Early explosions are facilitated by high mass-accretion rates after bounce and correspondingly high neutrino luminosities combined with a pronounced drop of the accretion rate and ram pressure at the Si/Si–O interface. Because of rapidly shrinking neutron star radii and receding shock fronts after the passage through their maxima, our models exhibit short advection timescales, which favor the efficient growth of the standing accretion-shock instability. The latter plays a supportive role at least for the initiation of the re-expansion of the stalled shock before runaway. Taking into account the effects of turbulent pressure in the gain layer, we derive a generalized condition for the critical neutrino luminosity that captures the explosion behavior of all models very well. We validate the robustness of our findings by testing the influence of stochasticity, numerical resolution, and approximations in some aspects of the microphysics.

The dynamics of neutrino-driven supernova explosions after shock revival in 2D and 3D

We study the growth of the explosion energy after shock revival in neutrino-driven explosions in two and three dimensions (2D/3D) using multi-group neutrino hydrodynamics simulations of an 11.2 M⊙ star. The 3D model shows a faster and steadier growth of the explosion energy and already shows signs of subsiding accretion after one second. By contrast, the growth of the explosion energy in 2D is unsteady, and accretion lasts for several seconds as confirmed by additional long-time simulations of stars of similar masses. Appreciable explosion energies can still be reached, albeit at the expense of rather high neutron star masses. In 2D, the binding energy at the gain radius is larger because the strong excitation of downward-propagating g modes removes energy from the freshly accreted material in the downflows. Consequently, the mass outflow rate is considerably lower in 2D than in 3D. This is only partially compensated by additional heating by outward-propagating acoustic waves in 2D. Moreover, the mass outflow rate in 2D is reduced because much of the neutrino energy deposition occurs in downflows or bubbles confined by secondary shocks without driving outflows. Episodic constriction of outflows and vertical mixing of colder shocked material and hot, neutrino-heated ejecta due to Rayleigh–Taylor instability further hamper the growth of the explosion energy in 2D. Further simulations will be necessary to determine whether these effects are generic over a wider range of supernova progenitors.

Outer crust of a cold non-accreting magnetar

The outer-crust structure and composition of a cold, non-accreting magnetar are studied. We model the outer crust to be made of fully equilibrated matter where ionized nuclei form a Coulomb crystal embedded in an electron gas. The main effects of the strong magnetic field are those of quantizing the electron motion in Landau levels and of modifying the nuclear single-particle levels producing, on average, an increased binding of nucleons in nuclei present in the Coulomb lattice. The effect of a homogeneous and constant magnetic field on nuclear masses has been predicted by using a covariant density functional in which induced currents and axial deformation due to the presence of a magnetic field that breaks time-reversal symmetry have been included self-consistently in the nucleon and meson equations of motion. Although not yet observed, for Ba 1016 G both effects contribute to produce different compositions - odd-mass nuclei are frequently predicted - and to increase the neutron-drip pressure as compared to a typical neutron star. Specifically, in such a regime, the magnetic-field effects on nuclei favor the appearance of heavier nuclei at low pressures. As B increases, such heavier nuclei are also preferred up to larger pressures. For the most extreme magnetic field considered, B=1018 G, and for the models studied, almost the whole outer crust is made of 4092Zr52.

Exploring the white dwarf cooling sequence in Globular Clusters

White dwarfs (WDs) are electron-degenerate structures that are commonly assumed to evolve via a pure cooling process, with no stable thermonuclear activity at work. Their cooling rate is adopted as a cosmic chronometer to constrain the age of several Galactic populations, including the disk, Globular Clusters (GCs) and open clusters. This thesis work is aimed at the study of the WD populations in globular clusters and is articulated in two branches. The first was focused on the study of the bright portion of the WD cooling sequence. By analyzing high resolution UV data acquired with the Hubble Space Telescope (HST), we compared the WD luminosity functions (LFs) in four Galactic GCs (namely M13, M3, NGC6752, and M5) finding an unexpected over-abundance of WDs in M13 and NGC6752 with respect to M3 and M5. Theoretical models suggest that, consistently with the blue-tail horizontal branch (HB) morphology of M13 and NGC6752, this overabundance is due to a population of slowly cooling WDs, i.e., WDs fading more slowly than in a pure cooling process thanks to an extra-energy source provided by stable thermonuclear burning in their residual hydrogen-rich envelope. This is the first empirical evidence of WDs fading at a slower rate than usually assumed, and has a crucial impact on the use of the cooling sequence as a cosmic chronometer. The second branch was focused on the search for the companion star to binary millisecond Pulsars (MSP) in the globular clusters M13 and NGC 6652: the identified companions turned out to be helium-core WDs, and provided a invaluable constraints on the mass of the neutron star and the epoch of the MSP formation.

Equation of state for the magnetic-color-flavor-locked phase and its implications for compact star models

Using the solutions of the gap equations of the magnetic-color-flavor-locked (MCFL) phase of paired quark matter in a magnetic field, and taking into consideration the separation between the longitudinal and transverse pressures due to the field-induced breaking of the spatial rotational symmetry, the equation of state of the MCFL phase is self-consistently determined. This result is then used to investigate the possibility of absolute stability, which turns out to require a field-dependent ""bag constant"" to hold. That is, only if the bag constant varies with the magnetic field, there exists a window in the magnetic field vs bag constant plane for absolute stability of strange matter. Implications for stellar models of magnetized (self-bound) strange stars and hybrid (MCFL core) stars are calculated and discussed.

A Be-type star with a black-hole companion

Stellar-mass black holes have all been discovered through X-ray emission, which arises from the accretion of gas from their binary companions (this gas is either stripped from low-mass stars or supplied as winds from massive ones). Binary evolution models also predict the existence of black holes accreting from the equatorial envelope of rapidly spinning Be-type stars (stars of the Be type are hot blue irregular variables showing characteristic spectral emission lines of hydrogen). Of the ~80 Be X-ray binaries known in the Galaxy, however, only pulsating neutron stars have been found as companions. A black hole was formally allowed as a solution for the companion to the Be star MWC 656 (also known as HD 215227), although that was based on a single radial velocity curve of the Be star, a mistaken spectral classification and rough estimates of the inclination angle. Here we report observations of an accretion disk line mirroring the orbit of the Be star. This, together with an improved radial velocity curve of the Be star through fitting sharp Fe II profiles from the equatorial disk, and a refined Be classification (to that of a B1.5-B2 III star), reveals a black hole of 3.8 to 6.9 solar masses orbiting MWC 656, the candidate counterpart of the gamma-ray source AGL J2241+4454. The black hole is X-ray quiescent and fed by a radiatively inefficient accretion flow giving a luminosity less than 1.6 x 10-7 times the Eddington luminosity. This implies that Be binaries with black-hole companions are difficult to detect by conventional X-ray surveys.

Is SAX J1808.4-3658 a Strange Star?

The possibility of strange stars is one of the most important issues in the study of compact objects. Here we use the observations of the newly discovered millisecond x-ray pulsar SAX J1808.4-3658 to constrain the radius of the compact star. Comparing the mass-radius relation of SAX J1808.4-3658 with theoretical models for both neutron stars and strange stars, we argue that a strange star model could be more consistent with SAX J1808.4-3658, and suggest that it is a likely strange star candidate. Our results are useful in constraining microscopic chiral symmetry restoration parameters in the quantum chromodynamics (QCD) modeling of strange matter.