High Energy Astrophysical Phenomena (astro-ph.HE)

Abstract: The Galactic black hole candidate MAXI J0637-430 was first discovered by $\textit{MAXI/GSC}$ on 2019 November 02. We study the spectral properties of MAXI J0637-430 by using the archived $\textit{NuSTAR}$ data and $\textit{Swift}$/XRT data. After fitting the eight spectra by using a disk component and a powerlaw component model with absorption, we select the spectra with relatively strong reflection components for detailed X-ray reflection spectroscopy. Using the most state-of-art reflection model $\tt{relxillCp}$, the spectral fitting measures a black hole spin $\textit{a}_{\rm{*}} > 0.72$ and the inclination angle of the accretion disk $i$ = $46.1_{-5.3}^{+4.0}$ degrees, at 90 per cent confidence level. In addition, the fitting results show an extreme supersolar iron abundance. Combined with the fitting results of the reflection model $\tt{reflionx\_hd}$, we consider that this unphysical iron abundance may be caused by a very high density accretion disk ( $n_{\rm{e}} > 2.34 \times 10^{21}$ $\rm{cm}^{-3}$ ) or a strong Fe K$\alpha$ emission line. The soft excess is found in the soft state spectral fitting results, which may be an extra free-free heating effect caused by high density of the accretion disk. Finally, we discuss the robustness of black hole spin obtained by X-ray reflection spectroscopy. The result of relatively high spin is self-consistent with broadened Fe K$\alpha$ line. Iron abundance and disk density have no effect on the spin results.

Abstract: Pulsar timing arrays (PTAs) and the Laser Interferometer Space Antenna (LISA) will open complementary observational windows on massive black-hole binaries (MBHBs), i.e., with masses in the range $\sim 10^6 - 10^{10}\,$ M$_{\odot}$. While PTAs may detect a stochastic gravitational-wave background from a population of MBHBs, during operation LISA will detect individual merging MBHBs. To demonstrate the profound interplay between LISA and PTAs, we estimate the number of MBHB mergers that one can expect to observe with LISA by extrapolating direct observational constraints on the MBHB merger rate inferred from PTA data. For this, we postulate that the common noise currently detected in PTAs is an astrophysical background sourced by a single MBHB population. We then constrain the LISA detection rate, $\mathcal{R}$, in the mass-redshift space by combining our Bayesian-inferred merger rate with LISA's sensitivity to spin-aligned, inspiral-merger-ringdown waveforms. Using an astrophysically-informed formation model, we predict a 95$\%$ upper limit on the detection rate of $\mathcal{R} < 134\,{\rm yr}^{-1}$ for binaries with total masses in the range $10^7 - 10^8\,$ M$_{\odot}$. For higher masses, i.e., $>10^8\,$ M$_{\odot}$, we find $\mathcal{R} < 2\,(1)\,\mathrm{yr}^{-1}$ using an astrophysically-informed (agnostic) formation model, rising to $11\,(6)\,\mathrm{yr}^{-1}$ if the LISA sensitivity bandwidth extends down to $10^{-5}$ Hz. Forecasts of LISA science potential with PTA background measurements should improve as PTAs continue their search.

Abstract: Possible strong first-order hadron-quark phase transitions in neutron star interiors leave an imprint on gravitational waves, which could be detected with planned third-generation interferometers. Given a signal from the late inspiral of a binary neutron star (BNS) coalescence, %the possibility of assessing the presence of such a phase transition depends on the precision that can be attained in the determination of the tidal deformability parameter, as well as on the model used to describe the hybrid star equation of state. For the latter, we employ here a phenomenological meta-modelling of the equation of state that largely spans the parameter space associated with both the low density phase and the quark high density compatible with current constraints. We show that with a network of third-generation detectors, a single loud BNS event might be sufficient to infer the presence of a phase transition at low baryon densities with an average Bayes factor $B\approx 100$, up to a luminosity distance ($\mathcal{D}_L \lesssim$ 300 Mpc).

Abstract: In MeerKAT observations pointed at a Galactic X-ray binary located on the Galactic plane we serendipitously discovered a radio nebula with cometary-like morphology. The feature, which we named `the Mini Mouse' based on its similarity with the previously discovered `Mouse' nebula, points back towards the previously unidentified candidate supernova remnant G45.24$+$0.18. We observed the location of the Mini Mouse with MeerKAT in two different observations, and we localised with arcsecond precision the 138 ms radio pulsar PSR J1914+1054g, recently discovered by the FAST telescope, to a position consistent with the head of the nebula. We confirm a dispersion measure of about 418 pc cm$^{-3}$ corresponding to a distance between 7.8 and 8.8 kpc based on models of the electron distribution. Using our accurate localisation and 2 period measurements spaced 90 days apart we calculate a period derivative of (2.7 $\pm$ 0.3) $\times$ 10 $^{-14}$ s s$^{-1}$. We derive a characteristic age of approximately 82 kyr and a spin down luminosity of 4$\times$10$^{35}$ erg s$^{-1}$, respectively. For a pulsar age comparable with the characteristic age, we find that the projected velocity of the neutron star is between 320 and 360 km/s if it was born at the location of the supernova remnant. The size of the proposed remnant appears small if compared with the pulsar characteristic age, however the relatively high density of the environment near the Galactic plane could explain a suppressed expansion rate and thus a smaller remnant.

Abstract: Core-collapse supernova explosions play a wide role in astrophysics by producing compact remnants (neutron stars, black holes) and the synthesis and injection of many heavy elements into their host Galaxy. Because they are produced in some of the most extreme conditions in the universe, they can also probe physics in extreme conditions (matter at nuclear densities and extreme temperatures and magnetic fields). To quantify the impact of supernovae on both fundamental physics and our understanding of the Universe, we must leverage a broad set of observables of this engine. In this paper, we study a subset of these probes using a suite of 1-dimensional, parameterized mixing models: ejecta remnants from supernovae, ultraviolet, optical and infra-red lightcurves, and transient gamma-ray emission. We review the other diagnostics and show how the different probes tie together to provide a more clear picture of the supernova engine.

Abstract: We conduct a three-dimensional hydrodynamical simulation of a common envelope evolution (CEE) where a neutron star (NS) spirals-in inside the envelope of a red supergiant (RSG) star in a predetermined orbit. We find that the jets shed pairs of vortices in an expanding spiral pattern, inflate two expanding spirally-shaped low-density bubbles, one above and one below the equatorial plane, and deposit angular momentum to the envelope. In the simulation we do not include the gravity of the NS such that all effects we find are solely due to the jets that the spiralling-in NS launches. The angular momentum that the jets deposit to the envelope is of the same order of magnitude as the orbital angular momentum and has the same direction. The turbulence that the jets induce in the common envelope might play a role in transporting energy and angular momentum. The jet-deposited energy that is radiated away (a process not studied here) leads to a transient event that is termed common envelope jets supernova (CEJSN) and might mimic an energetic core collapse supernova. The turbulence and the spiral pattern that we explore here might lead to bumps in the late light curve of the CEJSN when different segments of the ejected envelope collide with each other. This study emphasises the roles that jets can play in CEE (including jets launched by black hole companions) and adds to the rich variety of processes in CEJSN events.

Abstract: We present a novel fully Bayesian analysis to constrain short gamma-ray burst jet structures associated with cocoon, wide-angle and simple top-hat jet models, as well as the binary neutron star merger rate. These constraints are made given the distance and inclination information from GW170817, observed flux of GRB170817A, observed rate of short gamma-ray bursts detected by Swift, and the neutron star merger rate inferred from LIGO's first and second observing runs. A separate analysis is conducted where a fitted short gamma-ray burst luminosity function is included to provide further constraints. The jet structure models are further constrained using the observation of GW190425 and we find that the assumption that it produced a GRB170817-like short gamma-ray burst that went undetected due to the jet geometry is consistent with previous observations. We find and quantify evidence for low luminosity and wide-angled jet structuring in the short gamma-ray burst population, independently from afterglow observations, with log Bayes factors of $0.45{-}0.55$ for such models when compared to a classical top-hat jet. Slight evidence is found for a Gaussian jet structure model over all others when the fitted luminosity function is provided, producing log Bayes factors of $0.25{-}0.9\pm0.05$ when compared to the other models. However without considering GW190425 or the fitted luminosity function, the evidence favours a cocoon-like model with log Bayes factors of $0.14\pm0.05$ over the Gaussian jet structure. We provide new constraints to the binary neutron star merger rates of $1{-}1300$Gpc$^{-3}$yr$^{-1}$ or $2{-}680$Gpc$^{-3}$yr$^{-1}$ when a fitted luminosity function is assumed.

Abstract: We present the first 3D fully coupled magneto-thermal simulations of neutron stars (including the most realistic background structure and microphysical ingredients so far) applied to a very complex initial magnetic field topology in the crust, similar to what recently obtained by proto-neutron star dynamo simulations. In such configurations, most of the energy is stored in the toroidal field, while the dipolar component is a few percent of the mean magnetic field. This initial feature is maintained during the long-term evolution (1e6 yr), since the Hall term favours a direct cascade (compensating for Ohmic dissipation) rather than a strong inverse cascade, for such an initial field topology. The surface dipolar component, responsible for the dominant electromagnetic spin-down torque, does not show any increase in time, when starting from this complex initial topology. This is at contrast with the timing properties of young pulsars and magnetars which point to higher values of the surface dipolar fields. A possibility is that the deep-seated magnetic field (currents in the core) is able to self-organize in large scales (during the collapse or in the early life of a neutron star). Alternatively, the dipolar field might be lower than is usually thought, with magnetosphere substantially contributing to the observed high spin-down, via e.g., strong winds or strong coronal magnetic loops, which can also provide a natural explanation to the tiny surface hotspots inferred from X-ray data.