Astrophysics

We present Hubble Space Telescope imaging of a pre-explosion counterpart to SN 2019yvr obtained 2.6 years before its explosion as a type Ib supernova (SN Ib). Aligning to a post-explosion Gemini-S/GSAOI image, we demonstrate that there is a single source consistent with being the SN 2019yvr progenitor system, the second SN Ib progenitor candidate after iPTF13bvn. We also analyzed pre-explosion Spitzer/IRAC imaging, but we do not detect any counterparts at the SN location. SN 2019yvr was highly reddened, and comparing its spectra and photometry to those of other, less extinguished SNe Ib we derive E(B?�V)=0.51 +0.27 ??.16 mag for SN 2019yvr. Correcting photometry of the pre-explosion source for dust reddening, we determine that this source is consistent with a log(L/ L ??)=5.3±0.2 and T eff =6800 +400 ??00 K star. This relatively cool photospheric temperature implies a radius of 320 +30 ??0 R ??, much larger than expectations for SN Ib progenitor stars with trace amounts of hydrogen but in agreement with previously identified SN IIb progenitor systems. The photometry of the system is also consistent with binary star models that undergo common envelope evolution, leading to a primary star hydrogen envelope mass that is mostly depleted but seemingly in conflict with the SN Ib classification of SN 2019yvr. SN 2019yvr had signatures of strong circumstellar interaction in late-time ( > 150 day) spectra and imaging, and so we consider eruptive mass loss and common envelope evolution scenarios that explain the SN Ib spectroscopic class, pre-explosion counterpart, and dense circumstellar material. We also hypothesize that the apparent inflation could be caused by a quasi-photosphere formed in an extended, low-density envelope or circumstellar matter around the primary star.

Neutrino-matter interactions play an important role in core-collapse supernova (CCSN) explosions as they contribute to both lepton number and/or four-momentum exchange between neutrinos and matter, and thus act as the agent for neutrino-driven explosions. Due to the multiscale nature of neutrino transport in CCSN simulations, an implicit treatment of neutrino-matter interactions is desired, which requires solutions of coupled nonlinear systems in each step of the time integration scheme. In this paper we design and compare nonlinear iterative solvers for implicit systems with energy coupling neutrino-matter interactions commonly used in CCSN simulations. Specifically, we consider electron neutrinos and antineutrinos, which interact with static matter configurations through the Bruenn~85 opacity set. The implicit systems arise from the discretization of a non-relativistic two-moment model for neutrino transport, which employs the discontinuous Galerkin (DG) method for phase-space discretization and an implicit-explicit (IMEX) time integration scheme. In the context of this DG-IMEX scheme, we propose two approaches to formulate the nonlinear systems -- a coupled approach and a nested approach. For each approach, the resulting systems are solved with Anderson-accelerated fixed-point iteration and Newton's method. The performance of these four iterative solvers has been compared on relaxation problems with various degree of collisionality, as well as proto-neutron star deleptonization problems with several matter profiles adopted from spherically symmetric CCSN simulations. Numerical results suggest that the nested Anderson-accelerated fixed-point solver is more efficient than other tested solvers for solving implicit nonlinear systems with energy coupling neutrino-matter interactions.

The Cygnus Superbubble (CSB) is a region of soft X-ray emission approximately 13 degrees wide in the direction of the local spiral arm. Such a large region might be the result of strong stellar winds and supernovae from nearby stellar nurseries, or it could be the result of a single event - a hypernova. HaloSat observed 4 non-overlapping 10 degree diameter fields in the CSB region over the 0.4-7 keV band. The CSB absorption and temperature was found to be consistent over all 4 fields, with a weighted average of 6.1x10^21 cm^-2 and 0.190 keV, respectively. These observations suggest that the CSB is a cohesive object with a singular origin. The total thermal energy for the CSB is estimated at 4x10^52 erg, based upon a shell-like physical model of the CSB. Absorption and distance estimates to Cyg OB associations are examined. The CSB absorption is found to be most consistent with the absorption seen in Cyg OB1, implying that the CSB lies at a similar distance of 1.1-1.4 kpc.

GW190814 was a compact object binary coalescence detected in gravitational waves by Advanced LIGO and Advanced Virgo that garnered exceptional community interest due to its excellent localization and the uncertain nature of the binary's lighter-mass component (either the heaviest known neutron star, or the lightest known black hole). Despite extensive follow up observations, no electromagnetic counterpart has been identified. Here we present new radio observations of 75 galaxies within the localization volume at ?t??5??66 days post-merger. Our observations cover ??2 % of the total stellar luminosity in the final localization volume and extend to later timescales than previously-reported searches, allowing us to place the deepest constraints to date on the existence of a radio afterglow from a highly off-axis relativistic jet launched during the merger (assuming that the merger occurred within the observed area). For a viewing angle of ??46 ??(the best-fit binary inclination derived from the gravitational wave signal) and assumed electron and magnetic field energy fractions of ϵ e =0.1 and ϵ B =0.01 , we can rule out a typical short gamma-ray burst-like Gaussian jet with isotropic-equivalent kinetic energy 2? 10 51 erg propagating into a constant density medium n??.01 cm ?? . These are the first limits resulting from a galaxy-targeted search for a radio counterpart to a gravitational wave event, and we discuss the challenges, and possible advantages, of applying similar search strategies to future events using current and upcoming radio facilities.

In this paper, we employ Machine Learning algorithms on multi-mission observations for the classification of accretion states of outbursting Black hole X-ray binaries for the first time. Archival data from RXTE, Swift, MAXI and AstroSat observatories are used to generate the hardness intensity diagrams (HIDs) for outbursts of the sources XTE J1859+226 (1999 outburst), GX 339-4 (2002, 2004, 2007 and 2010 outbursts), IGR J17091-3624 (2016 outburst), and MAXI J1535-571 (2017 outburst). Based on variation of X-ray flux, hardness ratios, presence of various types of Quasi-periodic Oscillations (QPOs), photon indices and disk temperature, we apply clustering algorithms like K-Means clustering and Hierarchical clustering to classify the accretion states (clusters) of each outburst. As multiple parameters are involved in the classification process, we show that clustering algorithms club together the observations of similar characteristics more efficiently than the `standard' method of classification. We also infer that K-Means clustering provides more reliable results than Hierarchical clustering. We demonstrate the importance of the classification based on machine learning by comparing it with results from `standard' classification.

Similar to black-hole X-ray binary (BHXRB) transients, hysteresis-like state transitions are also seen in some neutron-star X-ray binaries (NSXRBs). Using a method based on wavelets and lightcurves constructed from archival RXTE observations, we extract a minimal timescale over the complete range of transitions for 4U 1608-52 during the 2002 and 2007 outbursts and the 1999 and 2000 outbursts for Aql X-1. We present evidence for a strong positive correlation between this minimal timescale and a similar timescale extracted from the corresponding power spectra of these sources.

This paper describes an analysis of the NuSTAR data of the fastest-rotating magnetar 1E 1547 ??5408, acquired in 2016 April for a time lapse of 151 ks. The source was detected with a 1-60 keV flux of 1.7? 10 ??1 ergs s ?? cm ?? , and its pulsation at a period of 2.086710(5) sec. In 8-25 keV, the pulses were phase-modulated with a period of T=36.0±2.3 ks, and an amplitude of ??.2 sec. This reconfirms the Suzaku discovery of the same effect at T= 36.0 +4.5 ??.5 ks, made in the 2009 outburst. These results strengthen the view derived from the Suzaku data, that this magnetar performs free precession as a result of its axial deformation by ??.6? 10 ?? , possibly caused by internal toroidal magnetic fields reaching ??10 16 G. Like in the Suzaku case, the modulation was not detected in energies below ?? keV. Above 10 keV, the pulse-phase behaviour, including the 36 ks modulation parameters, exhibited complex energy dependences: at ??2 keV, the modulation amplitude increased to ??.5 sec, and the modulation phase changed by ??65 ??over 10--27 keV, followed by a phase reversal. Although the pulse significance and pulsed fraction were originally very low in >10 keV, they both increased noticeably, when the arrival times of individual photons were corrected for these systematic pulse-phase variations. Possible origins of these complex phenomena are discussed, in terms of several physical processes that are specific to ultra-strong magnetic fields.

We report the detection of a probable γ -ray quasi-periodic oscillation (QPO) of around 314 days in the monthly binned 0.1 -- 300 GeV γ -ray {\it Fermi}-LAT light curve of the well known BL Lac blazar OJ 287. To identify and quantify the QPO nature of the γ -ray light curve of OJ 287, we used the Lomb-Scargle periodogram (LSP), REDFIT, and weighted wavelet z-transform (WWZ) analyses. We briefly discuss possible emission models for radio-loud active galactic nuclei (AGN) that can explain a γ -ray QPO of such a period in a blazar. Reports of changes in the position of quasi-stationary radio knots over a yearly timescale as well as a strong correlation between gamma-ray and mm-radio emission in previous studies indicate that the signal is probably associated with these knots.

Approximately 8% of the ??2800 known pulsars exhibit "nulling," a temporary broadband cessation of normal pulsar emission. Nulling behaviour can be coarsely quantified by the nulling fraction, which describes the percentage of time a given pulsar will be found in a null state. In this paper, we perform the most thorough statistical analysis thus far of the properties of 141 known nulling pulsars. We find weak, non-linear correlations between nulling fraction and pulse width, as well as nulling fraction and spin period which could be attributed to selection effects. We also further investigate a recently-hypothesized gap at 40% nulling fraction. While a local minimum does exist in the distribution, we cannot confirm a consistent and unique break in the distribution when we investigate with univariate and multivariate clustering methods, nor can we prove the existence of two statistically distinct populations about this minimum. Using the same methods, we find that nulling pulsars are a statistically different population from normal, radio, non-nulling pulsars, which has never been quantitatively verified. In addition, we summarize the findings of the prior nulling pulsar statistics literature, which are notoriously contradictory. This study, in context, furthers the idea that nulling fraction alone does not contain enough information to describe the behaviour of a nulling pulsar and that other parameters such as null lengths and null randomness, in addition to a better understanding of selection effects, are required to fully understand this phenomenon.

A sample of narrow-line (NLS1) and broad-line Seyfert 1 (BLS1) galaxies observed with Suzaku is presented. The final sample consists of 22 NLS1s and 47 BLS1s, for a total of 69 AGN that are all at low redshift (z<0.5) and exhibit low host galaxy column densities (<10^22 cm^-2). The average spectrum for each object is fit with a toy model to characterise important parameters, including the photon index, soft excess, Compton hump (or hard excess), narrow iron line strength, luminosity and X-ray Eddington ratio (L_x/L_Edd). We confirm previous findings that NLS1s have steeper power laws and higher X-ray Eddington ratios, but also find that NLS1 galaxies have stronger soft and hard excesses than their BLS1 counterparts. Studying the correlations between parameters shows that the soft and hard excesses are correlated for NLS1 galaxies, while no such correlation is observed for BLS1s. Performing a principal component analysis (PCA) on the measured X-ray parameters shows that while the X-ray Eddington ratio is the main source of variations within our sample (PC1), variations in the soft and hard excesses form the second principal component (PC2) and it is dominated by the NLS1s. The correlation between the soft and hard excess in NLS1 galaxies may suggest a common origin for the two components, such as a blurred reflection model. The presented Suzaku sample of Seyfert 1 galaxies is a useful tool for analysis of the X-ray properties of AGN, and for the study of the soft and hard excesses observed in AGN.