Central Compact Object Research Articles

Dust scattering halo around the CCO in HESS J1731–347: A detailed analysis

The supernova remnant (SNR) HESS J1731−347 is one of the few objects exhibiting emission up to the TeV energy band and it stands as a prime target for the study of cosmic ray acceleration in SNRs. It also hosts a central compact object (CCO), which is of interest in the context of the ultra-dense matter equation of state in neutron stars. For both types of studies, however, the parameters of the respective models depend crucially on the assumed distance to the object and are affected to a certain extent by the assumed interstellar medium (ISM) properties around the SNR. Here, we report on the first quantitative analysis of the properties of the compact X-ray dust scattering halo that is assumed to be present around the CCO based on Chandra observations of the source. Our findings unambiguously confirm the presence of a compact halo around the CCO, and we show that the observed halo properties are consistent with expectations from independent measurements of the dust distribution along the line of sight and the distance to the source. Although we were not able to significantly improve those constraints, our results are important for future studies of the CCO itself. Indeed, the halo contribution is expected to affect the X-ray spectrum and the derived parameters of the neutron star when observed with moderate angular resolutions. Our results, which offer a quantitative characterization of the halo properties, will be useful in accounting for this effect.

Measuring the Nonaxially Symmetric Surface Temperature Distribution of the Central Compact Object in Puppis A

Abstract The surface temperature distributions of central compact objects (CCOs) are powerful probes of their crustal magnetic field strengths and geometries. Here we model the surface temperature distribution of RX J0822−4300, the CCO in the Puppis A supernova remnant, using 471 ks of XMM-Newton data. We compute the energy-dependent pulse profiles in 16 energy bands, fully including the general relativistic effects of gravitational redshift and light bending, to accurately model the two heated surface regions of different temperatures and areas, in addition to constraining the viewing geometry. This results in precise measurements of the two temperatures: kT warm = ( 1 + z ) × 0.222 − 0.019 + 0.018 keV and kT hot = (1 + z) × 0.411 ± 0.011 keV. The two heated surface regions are likely located very close to the rotational poles, with the most probable position of the hotter component ≈ 6° from the rotational pole. For the first time, we are able to measure a deviation from a pure antipodal hot-spot geometry, with a longitudinal offset δ γ = 11 .° 7 − 2 .° 5 + 2 .° 6 . The discovery of this asymmetry, along with the factor of ≈2 temperature difference between the two emitting regions, may indicate that RX J0822−4300 was born with a strong, tangled crustal magnetic field.

GRMHD Simulations and Modeling for Jet Formation and Acceleration Region in AGNs

Relativistic jets are collimated plasma outflows with relativistic speeds. Astrophysical objects involving relativistic jets are a system comprising a compact object such as a black hole, surrounded by rotating accretion flows, with the relativistic jets produced near the central compact object. The most accepted models explaining the origin of relativistic jets involve magnetohydrodynamic (MHD) processes. Over the past few decades, many general relativistic MHD (GRMHD) codes have been developed and applied to model relativistic jet formation in various conditions. This short review provides an overview of the recent progress of GRMHD simulations in generating relativistic jets and their modeling for observations.

What causes the absence of pulsations in Central Compact Objects in supernova remnants?

Most young neutron stars belonging to the class of Central Compact Objects (CCOs) in supernova remnants do not have known periodicities. We investigated seven such CCOs to understand the common reasons for the absence of detected pulsations. Making use of XMM-Newton, Chandra, and NICER observations, we perform a systematic timing and spectral analysis to derive updated sensitivity limits for both periodic signals and multi-temperature spectral components that could be associated with radiation from hotspots on the neutron star surface. Based on these limits, we then investigated for each target the allowed viewing geometry that could explain the lack of pulsations. We find that it is unlikely (< 10−6) to attribute that we do not see pulsations to an unfavorable viewing geometry for five considered sources. Alternatively, the carbon atmosphere model, which assumes homogeneous temperature distribution on the surface, describes the spectra equally well and provides a reasonable interpretation for the absence of detected periodicities within current limits. The unusual properties of CCOs with respect to other young neutron stars could suggest a different evolutionary path, as that proposed for sources experiencing episodes of significant fallback accretion after the supernova event.

Ultra-delayed Neutrino-driven Explosion of Rotating Massive-star Collapse

Abstract Long-term neutrino-radiation hydrodynamics simulations in full general relativity are performed for the collapse of rotating massive stars that are evolved from He-stars with initial masses of 20 and 32 M ⊙. It is shown that if the collapsing stellar core has sufficient angular momentum, the rotationally supported proto-neutron star (PNS) survives for seconds accompanying the formation of a massive torus of mass larger than 1 M ⊙. Subsequent mass accretion onto the central region produces a massive and compact central object, and eventually enhances the neutrino luminosity beyond 1053 erg s−1, resulting in a very delayed neutrino-driven explosion, in particular toward the polar direction. The kinetic energy of the explosion can be appreciably higher than 1052 erg for a massive progenitor star and compatible with that of energetic supernovae like broad-line type-Ic supernovae. By the subsequent accretion, the massive PNS collapses eventually into a rapidly spinning black hole, which could be a central engine for gamma-ray bursts if a massive torus surrounds it.

A Necessary Condition for Supernova Fallback Invading Newborn Neutron-star Magnetosphere

Abstract We numerically investigate the dynamics of a supernova fallback accretion confronting with a relativistic wind from a newborn neutron star (NS). The time evolution of the accretion shock in the radial direction is basically characterized by the encounter radius of the flow r enc and a dimensionless parameter ζ ≡ L / M ̇ fb c 2 , where L is the NS wind luminosity and M ̇ fb is the fallback mass accretion rate. We find that the critical condition for the fallback matter to reach near the NS surface can be simply described as ζ &lt; ζ min ≡ GM * / c 2 r enc or r enc L / GM * M ̇ fb &lt; 1 independent of the wind Lorentz factor, where M * is the NS mass. With combining the condition for the fallback matter to bury the surface magnetic field under the NS crust, we discuss the possibility that the trifurcation of NSs into rotation-powered pulsars, central compact objects, and magnetars can be induced by supernova fallback.

Maximum accretion rate of supermassive stars

Context. The formation of the most massive quasars observed at high redshifts requires extreme inflows of gas down to the length scales of the central compact object. Aims. Here we estimate the maximum inflow rate allowed by gravity down to the surface of supermassive stars, the possible progenitors of these supermassive black holes. Methods. We use the continuity equation and the assumption of spherical symmetry and free fall to derive the maximum allowed inflow rates for various density profiles. We apply our approach to the mass–radius relation of rapidly accreting supermassive stars to estimate an upper limit to the accretion rates allowed during the formation of these objects. Results. We find that, as long as the density of the accreted gas is smaller than or equal to the average density of the accretor, the maximum allowed rate, Ṁmax, is given uniquely by the compactness of the accretor. We argue that a density inversion between accreting matter and the accretor is inconsistent with gravitational collapse. For the compactness of rapidly accreting supermassive stars, Ṁmax is related to the stellar mass, M, by a power law, Ṁmax ∝ M3/4. The rates of atomically cooled halos (0.1−10 M⊙ yr−1) are allowed as soon as M ≳ 1 M⊙. The largest rates expected in galaxy mergers (104 − 105 M⊙ yr−1) become accessible once the accretor is supermassive (M ≳ 104 M⊙). Conclusions. These results suggest that supermassive stars can accrete up to masses &gt; 106 M⊙ before they collapse via the general-relativistic instability. At such masses, the collapse is expected to lead to the direct formation of a supermassive black hole, even within metal-rich gas, resulting in a black hole seed that is significantly heavier than in conventional direct collapse models for atomic cooling halos.

Fallback Accretion Model for the Years-to-decades X-Ray Counterpart to GW170817

Abstract A new component was reported in the X-ray counterpart to the binary neutron star merger and gravitational-wave event GW170817, exceeding the afterglow emission from an off-axis structured jet. The afterglow emission from the kilonova/macronova ejecta may explain the X-ray excess but exceeds the radio observations if the spectrum is the same. We propose a fallback accretion model that a part of ejecta from the neutron star merger falls back and forms a disk around the central compact object. In the super-Eddington accretion phase, the X-ray luminosity stays near the Eddington limit of a few solar masses and the radio is weak, as observed. This will be followed by a power-law decay. The duration of the constant luminosity phase conveys the initial fallback timescale t 0 in the past. The current multiyear duration requires t 0 &gt; 3–30 s, suggesting that the disk wind rather than the dynamical ejecta falls back after the jet launch. Future observations in the next decades will probe the timescale of t 0 ∼ 10–104 s, around the time of extended emission in short gamma-ray bursts. The fallback accretion has not been halted by the r-process heating, implying that fission is weak on the year scale. We predict that the X-ray counterpart will disappear in a few decades due to the r-process halting or the depletion of fallback matter.

X-ray bounds on cooling, composition, and magnetic field of the Cassiopeia A neutron star and young central compact objects

ABSTRACT We present analysis of multiple Chandra and XMM–Newton spectra, separated by 9–19 yr, of four of the youngest central compact objects (CCOs) with ages &amp;lt;2500 yr: CXOU J232327.9+584842 (Cassiopeia A), CXOU J160103.1−513353 (G330.2+1.0), 1WGA J1713.4−3949 (G347.3−0.5), and XMMU J172054.5−372652 (G350.1−0.3). By fitting these spectra with thermal models, we attempt to constrain each CCO’s long-term cooling rate, composition, and magnetic field. For the CCO in Cassiopeia A, 14 measurements over 19 yr indicate a decreasing temperature at a 10-yr rate of 2.2 ± 0.2 or 2.8 ± 0.3 per cent (1σ error) for a constant or changing X-ray absorption, respectively. We obtain cooling rate upper limits of 17 per cent for CXOU J160103.1−513353 and 6 per cent for XMMU J172054.5−372652. For the oldest CCO, 1WGA J1713.4−3949, its temperature seems to have increased by 4 ± 2 per cent over a 10-yr period. Assuming each CCO’s preferred distance and an emission area that is a large fraction of the total stellar surface, a non-magnetic carbon atmosphere spectrum is a good fit to spectra of all four CCOs. If distances are larger and emission areas are somewhat smaller, then equally good spectral fits are obtained using a hydrogen atmosphere with B ≤ 7 × 1010 G or B ≥ 1012 G for CXOU J160103.1−513353 and B ≤ 1010 G or B ≥ 1012 G for XMMU J172054.5−372652 and non-magnetic hydrogen atmosphere for 1WGA J1713.4−3949. In a unified picture of CCO evolution, our results suggest most CCOs, and hence a sizable fraction of young neutron stars, have a surface magnetic field that is low early in their life but builds up over several thousand years.

A kinematic study of central compact objects and their host supernova remnants

Context. Central compact objects (CCOs) are a peculiar class of neutron stars, primarily encountered close to the center of young supernova remnants (SNRs) and characterized by thermal X-ray emission. Measurements of their proper motion and the expansion of the parent SNR are powerful tools for constraining explosion kinematics and the age of the system. Aims. Our goal is to perform a systematic study of the proper motion of all known CCOs with appropriate data available. From this, we hope to obtain constraints on the violent kick acting on the neutron star during the supernova explosion and infer the exact site of the explosion within the SNR. In addition, we aim to measure the expansion of three SNRs within our sample to obtain a direct handle on their kinematics and age. Methods. We analyze multiple archival Chandra data sets that consist of HRC and ACIS observations separated by temporal baselines of between 8 and 15 years. We achieve accurate source positions by fitting the imaging data with ray-tracing models of the Chandra point spread function. In order to correct for Chandra’s systematic astrometric uncertainties, we establish a reference frame using X-ray-detected sources in Gaia DR2 to provide accurate proper motion estimates for our target CCOs. Complementarily, we use our coaligned data sets to trace the expansion of three SNRs by directly measuring the spatial offset of various filaments and ejecta clumps between different epochs. Results. In total, we present new proper motion measurements for six CCOs. Within our sample, we do not find any indication of a hypervelocity object, and we determine comparatively tight upper limits (&lt; 230 km s−1) on the transverse velocities of the CCOs in G330.2+1.0 and RX J1713.7−3946. We tentatively identify direct signatures of expansion for the SNRs G15.9+0.2 and Kes 79, at an estimated significance of 2.5σ and 2σ, respectively. Moreover, we confirm recently reported results, measuring the rapid expansion of G350.1−0.3 at almost 6000 km s−1, which places its maximal age at 600 − 700 years, making this object one of the youngest Galactic core-collapse SNRs. The observed expansion, combined with the proper motion of its CCO, which is much slower than previously predicted, implies the need for a very inhomogeneous circumstellar medium to explain the highly asymmetric appearance of the SNR. Finally, for the SNR RX J1713.7−3946, we combine previously published expansion measurements with our measurement of the CCO’s proper motion to obtain a constraining upper limit of 1700 years on the system’s age.

Short gamma-ray burst jet propagation in binary neutron star merger environments

ABSTRACT The multimessenger event GW170817/GRB 170817A confirmed that binary neutron star (BNS) mergers can produce short gamma-ray burst (SGRB) jets. This evidence promoted new investigations on the mechanisms through which a BNS merger remnant can launch such a powerful relativistic outflow and on the propagation of the latter across the surrounding post-merger environment. In particular, great strides have been made in jet propagation models, establishing connections between the initial jet launching conditions, including the incipient jet launching time (with respect to merger) and the injection parameters, and the observable SGRB prompt and afterglow emission. However, present semi-analytical models and numerical simulations (with one notable exception) adopt simple handmade prescriptions to account for the post-merger environment, lacking a direct association with any specific merging BNS system. Here, we present the first three-dimensional relativistic hydrodynamics simulations of incipient SGRB jets propagating through a post-merger environment that is directly imported from the outcome of a previous general relativistic BNS merger simulation. Our results show that the evolution and final properties of the jet can be largely affected by the anisotropies and the deviations from axisymmetry and homologous expansion characterizing more realistic BNS merger environments. In addition, we find that the inclusion of the gravitational pull from the central compact object, often overlooked, can have a major impact. Finally, we consider different jet launching times referred to the same BNS merger model and discuss the consequences for the ultimate jet properties.

Emission lines from X-ray illuminated accretion disc in black hole binaries

ABSTRACT X-ray flux from the inner hot region around central compact object in a binary system illuminates the upper surface of an accretion disc and it behaves like a corona. This region can be photoionized by the illuminating radiation, and thus can emit different emission lines. We study those line spectra in black hole X-ray binaries for different accretion flow parameters including its geometry. The varying range of model parameters captures maximum possible observational features. We also put light on the routinely observed Fe-line emission properties based on different model parameters, ionization rate, and Fe abundances. We find that the Fe-line equivalent width WE decreases with increasing disc accretion rate and increases with the column density of the illuminated gas. Our estimated line properties are in agreement with observational signatures.

The Center of Expansion and Age of the Oxygen-rich Supernova Remnant 1E 0102.2-7219

Abstract We present new proper-motion measurements of optically emitting oxygen-rich knots of supernova remnant 1E 0102.2-7219 (E0102), which are used to estimate the remnant’s center of expansion and age. Four epochs of high-resolution Hubble Space Telescope images spanning 19 yr were retrieved and analyzed. We found a robust center of expansion of α = 1h04m02.ˢ48 and δ = −72°01′53.″92 (J2000) with 1σ uncertainty of 1.″77 using 45 knots from images obtained with the Advanced Camera for Surveys using the F475W filter in 2003 and 2013 having the highest signal-to-noise ratio. We also estimate an upper limit explosion age of 1738 ± 175 yr by selecting the knots with the highest proper motions and these knots are assumed to be the least decelerated. We find evidence of an asymmetry in the proper motions of the knots as a function of position angle. We conclude that these asymmetries were most likely caused by interaction between E0102's original supernova blast wave and an inhomogeneous surrounding environment, as opposed to intrinsic explosion asymmetry. The observed nonhomologous expansion suggests that the use of a free expansion model inaccurately offsets the center of expansion and leads to an overestimated explosion age. We discuss our findings as they compare to previous age and center of expansion estimates of E0102, as well as their relevance to a recently identified candidate central compact object.

Super-Eddington accretion on to a stellar mass ultraluminous X-ray source NGC 4190 ULX1

ABSTRACT We present the results of high-quality XMM-NEWTON observations of a ultraluminous X-ray source (ULX) in the galaxy NGC 4190. The detection of spectral cutoff in NGC 4190 ULX1 spectra rules out the interpretation of the ULX to be in a standard low/hard canonical accretion state. We report that the high quality EPIC spectra can be better described by broad thermal component, such as a slim disc. In addition we found long-term spectral and flux variability in the source using several XMM-NEWTON and Swift data. A clear anticorrelation between flux and power-law photon index is found which further confirms the unusual spectral state evolution of the ULX. Spectral properties of the ULX suggest that the source is in a broadened disc state with luminosities [≈(3 − 10) × 1039 erg s−1] falling in the ultraluminous regime. The positive luminosity–temperature relation further suggests that the multicolour disc model follows the L ∝ T4 relation that is expected for a blackbody disc emission from a constant area and the slim disc model seems to favour L ∝ T2 relation consistent with an advection-dominated disc emission. From the broadened disc-like spectral feature at such luminosity, we estimated the upper limit of the mass of the central compact object from the inner disc radius and found that the ULX hosts a stellar mass black hole.

Time-scales to reach chemical equilibrium in ices at snowline distance around compact objects: the influence of accretion mass in the central object

ABSTRACT In this work, we analyse soft X-ray emission due to mass accretion on to compact stars and its effects on the time-scale to reach chemical equilibrium of eventual surrounding astrophysical ices exposed to that radiation. Reaction time-scales due to soft X-ray in water-rich and pure ices of methanol, acetone, acetonitrile, formic acid, and acetic acid were determined. For accretion rates in the range $\dot{m}=10^{-12}\!-\!10^{-8}\,{\rm M}_\odot$ yr−1 and distances in the range 1–3 LY from the central compact objects, the time-scales lie in the range 10–108 yr, with shorter time-scales corresponding to higher accretion rates. Obtained time-scales for ices at snow-line distances can be small when compared to the lifetime (or age) of the compact stars, showing that chemical equilibrium could have been achieved. Time-scales for ices to reach chemical equilibrium depend on X-ray flux and, hence, on accretion rate, which indicates that systems with low accretion rates may not have reached chemical equilibrium.

3D Magnetothermal Simulations of Tangled Crustal Magnetic Field in Central Compact Objects

Abstract Central compact objects (CCOs) are young neutron stars emitting thermal X-rays with bolometric luminosities L X in the range of 1032–1034 erg s−1. Gourgouliatos, Hollerbach, and Igoshev recently suggested that peculiar emission properties of CCOs can be explained by tangled magnetic field configurations formed in a stochastic dynamo during the proto–neutron star stage. In this case the magnetic field consists of multiple small-scale components with negligible contribution of global dipolar field. We study numerically three-dimensional magnetothermal evolution of tangled crustal magnetic fields in neutron stars. We find that all configurations produce complicated surface thermal patterns that consist of multiple small hot regions located at significant separations from each other. The configurations with initial magnetic energy of (2.5–10) × 1047 erg have temperatures of hot regions that reach ≈ 0.2 keV, to be compared with the bulk temperature of ≈ 0.1 keV in our simulations with no cooling. A factor of two in temperature is also seen in observations of CCOs. The hot spots produce periodic modulations in light curve with typical amplitudes of ≤9%–11%. Therefore, the tangled magnetic field configuration can explain thermal emission properties of some CCOs.

Fast radio burst repeaters produced via Kozai-Lidov feeding of neutron stars in binary systems

Neutron stars are likely surrounded by gas, debris, and asteroid belts. Kozai-Lidov perturbations, induced by a distant, but gravitationally bound companion, can trigger the infall of such orbiting bodies onto a central compact object. These effects could lead to the emission of fast radio bursts (FRBs), for example by asteroid-induced magnetic wake fields in the wind of the compact object. A few percent of binary neutron star systems in the Universe, such as neutron star-main sequence star, neutron star-white dwarf, double neutron star, and neutron star-black hole systems, can account for the observed non-repeating FRB rates. More remarkably, we find that wide and close companion orbits lead to non-repeating and repeating sources, respectively, and they allow for one to compute a ratio between repeating and non-repeating sources of a few percent, which is in close agreement with the observations. Three major predictions can be made from our scenario, which can be tested in the coming years: (1) most repeaters should stop repeating after a period between 10 years to a few decades, as their asteroid belts become depleted; (2) some non-repeaters could occasionally repeat, if we hit the short period tail of the FRB period distribution; and (3) series of sub-Jansky level short radio bursts could be observed as electromagnetic counterparts of the mergers of binary neutron star systems.

Expansion and Age of the Supernova Remnant G350.1–0.3: High-velocity Iron Ejecta from a Core-collapse Event

Abstract We report Chandra observations of the highly asymmetric core-collapse supernova remnant G350.1−0.3. We document expansion over 9 yr away from the roughly stationary central compact object, with sky-plane velocities up to 5000 d 4.5 km s−1 (d 4.5 is the distance in units of 4.5 kpc), redshifts ranging from 900 to 2600 km s−1, and three-dimensional space velocities approaching 6000 km s−1. Most of the bright emission comes from heavy-element ejecta particularly strong in iron. Iron-enhanced ejecta are seen at 4000–6000 km s−1, strongly suggesting that the supernova was not a common Type IIP event. While some fainter regions have roughly solar abundances, we cannot identify clear blast-wave features. Our expansion proper motions indicate that G350.1−0.3 is no more than about 600 yr old, independent of distance: the third youngest known core-collapse supernova in the Galaxy, and one of the most asymmetric.

Einstein’s Geometrical versus Feynman’s Quantum-Field Approaches to Gravity Physics: Testing by Modern Multimessenger Astronomy

Modern multimessenger astronomy delivers unique opportunity for performing crucial observations that allow for testing the physics of the gravitational interaction. These tests include detection of gravitational waves by advanced LIGO-Virgo antennas, Event Horizon Telescope observations of central relativistic compact objects (RCO) in active galactic nuclei (AGN), X-ray spectroscopic observations of Fe Kα line in AGN, Galactic X-ray sources measurement of masses and radiuses of neutron stars, quark stars, and other RCO. A very important task of observational cosmology is to perform large surveys of galactic distances independent on cosmological redshifts for testing the nature of the Hubble law and peculiar velocities. Forthcoming multimessenger astronomy, while using such facilities as advanced LIGO-Virgo, Event Horizon Telescope (EHT), ALMA, WALLABY, JWST, EUCLID, and THESEUS, can elucidate the relation between Einstein’s geometrical and Feynman’s quantum-field approaches to gravity physics and deliver a new possibilities for unification of gravitation with other fundamental quantum physical interactions.

On the X-Ray Properties of the Putative Central Compact Object in 1E 0102.2−7219

Abstract We have analyzed the archival Chandra X-ray Observatory observations of the compact feature in the Small Magellanic Cloud supernova remnant 1E 0102.2−7219, which has recently been suggested to be the central compact object remaining after the supernova explosion. In our analysis, we have used appropriate, time-dependent responses for each of the archival observations, have modeled the background instead of subtracting it, and have fit unbinned spectra to preserve the maximal spectral information. The spectrum of this feature is similar to the spectrum of the surrounding regions that have significantly enhanced abundances of O, Ne, and Mg. We find that the previously suggested blackbody model is inconsistent with the data, as Monte Carlo simulations indicate that more than 99% of the simulated data sets have a test statistic value lower than that of the data. The spectrum is described adequately by a nonequilibrium ionization thermal model with two classes of models that fit the data equally well. One class of models has a temperature of kT ∼ 0.79 keV, an ionization timescale of ∼3 × 1011 cm−3 s, and marginal evidence for enhanced abundances of O and Ne, and the other has a temperature of kT ∼ 0.91 keV, an ionization timescale of ∼7 × 1010 cm−3 s, and abundances consistent with local interstellar medium values. We also performed an image analysis and find that the spatial distribution of the counts is not consistent with that of a point source. The hypothesis of a point-source distribution can be rejected at the 99.9% confidence level. Therefore, this compact feature is most likely a knot of O- and Ne-rich ejecta associated with the reverse shock.