Search

Recent studies of nearby globular clusters have discovered excess dark mass in their cores, apparently in an extended distribution, and simulations indicate that this mass is composed mostly of white dwarfs (respectively stellar-mass black holes) in clusters that are core-collapsed (respectively with a flatter core). We perform mass-anisotropy modelling of the closest globular cluster, M4, with intermediate slope for the inner stellar density. We use proper-motion data from Gaia EDR3 and from observations by the Hubble Space Telescope. We extract the mass profile employing Bayesian Jeans modelling, and check our fits with realistic mock data. Our analyses return isotropic motions in the cluster core and tangential motions ($\beta\approx -0.4$$\pm$$0.1$) in the outskirts. We also robustly measure a dark central mass of roughly $800\pm300 \,$M$_{\odot}$, but it is not possible to distinguish between a point-like source, such as an intermediate-mass black hole (IMBH), or a dark population of stellar remnants of extent $\approx 0.016\,\rm pc \simeq 3300\,AU$. However, when removing a high-velocity star from the cluster centre, the same mass excess is found, but more extended ($\sim 0.034\, \rm{pc} \approx 7000\,\rm AU$). We use Monte Carlo $N$-body models of M4 to interpret the second outcome, and find that our excess mass is not sufficiently extended to be confidently associated with a dark population of remnants. Finally, we discuss the feasibility of these two scenarios (i.e., IMBH vs. remnants), and propose new observations that could help to better grasp the complex dynamics in M4's core., Comment: 19 page, 15 figures, 3 tables. Accepted for publication in MNRAS

Gravitational wave observations of the inspiral of stellar-mass compact objects into massive black holes (MBHs), extreme mass ratio inspirals (EMRIs), enable precision measurements of parameters such as the MBH mass and spin. The Laser Interferometer Space Antenna is expected to detect sufficient EMRIs to probe the underlying source population, testing theories of the formation and evolution of MBHs and their environments. Population studies are subject to selection effects that vary across the EMRI parameter space, which bias inference results if unaccounted for. This bias can be corrected, but evaluating the detectability of many EMRI signals is computationally expensive. We mitigate this cost by (i) constructing a rapid and accurate neural network interpolator capable of predicting the signal-to-noise ratio of an EMRI from its parameters, and (ii) further accelerating detectability estimation with a neural network that learns the selection function, leveraging our first neural network for data generation. The resulting framework rapidly estimates the selection function, enabling a full treatment of EMRI detectability in population inference analyses. We apply our method to an astrophysically motivated EMRI population model, demonstrating the potential selection biases and subsequently correcting for them. Accounting for selection effects, we predict that LISA will measure the MBH mass function slope to a precision of 8.8%, the CO mass function slope to a precision of 4.6%, the width of the MBH spin magnitude distribution to a precision of 10% and the event rate to a precision of 12% with EMRIs at redshifts below z=6., 12 pages, 4 figures

The Granger causality is an econometric test for determining whether one time series is useful for forecasting another one with a certain Granger lag. Here, the light curves in the 0.3-1 keV (reflection dominated, soft) and 1.2-5 keV (continuum dominated, hard) bands of Active Galactic Nuclei (AGNs) are produced, and the Granger lags are estimated and compared to the traditional lag-frequency spectra. We find that the light curves in the hard band Granger-cause (lead) those in the soft band, whereas the obtained Granger lags could be interpreted as the intrinsic reverberation lags. Then, we extract the Granger-lag profiles from 14 XMM-Newton observations of IRAS 13224-3809, and find that the lags are significant in 12 observations. The majority of the obtained Granger (intrinsic) lags are ~ 200-500 s. With the IRAS 13224-3809 mass of 2 $\times 10^6$ $M_{\odot}$, these lags correspond to the true light-travel distance of ~ 20-50 $r_g$. Hence, by assuming a lamp-post geometry and a face-on disc, this places the corona at ~ 10-25 $r_{g}$ above the central black hole. Moreover, multiple Granger lags consisting of the small and large lags of < 500 s and > 1000 s are detected in 4 observations. This means that the corona height can significantly evolve from ~ 10-25 $r_{g}$ to ~ 55 $r_{g}$, or vice versa, towards the end of the observation. Granger lags are a promising way to measure the intrinsic lags, and provide evidence of coronal height variability within each individual observation., 13 pages, 10 figures, 1 Table, accepted for publication in MNRAS

We measure the black hole mass and investigate the accretion flow around the local ($z=0.0502$) quasar PG 1119+120. Spectroscopic monitoring with Calar Alto provides H$\beta$ lags and linewidths from which we estimate a black hole mass of $\log \left(M_{\bullet}/\mathrm{M}_{\odot} \right) = 7.0$, uncertain by $\sim0.4$ dex. High cadence photometric monitoring over two years with the Las Cumbres Observatory provides lightcurves in 7 optical bands suitable for intensive continuum reverberation mapping. We identify variability on two timescales. Slower variations on a 100-day timescale exhibit excess flux and increased lag in the $u'$ band and are thus attributable to diffuse bound-free continuum emission from the broad line region. Faster variations that we attribute to accretion disc reprocessing lack a $u'$-band excess and have flux and delay spectra consistent with either $\tau \propto \lambda^{4/3}$, as expected for a temperature structure of $T(R) \propto R^{-3/4}$ for a thin accretion disc, or $\tau \propto \lambda^{2}$ expected for a slim disc. Decomposing the flux into variable (disc) and constant (host galaxy) components, we find the disc SED to be flatter than expected with $f_{\nu} \sim \rm{const}$. Modelling the SED predicts an Eddington ratio of $\lambda_{\rm Edd} > 1$, where the flat spectrum can be reproduced by a slim disc with little dust extinction or a thin disc which requires more dust extinction. While this accretion is super-Eddington, the geometry is still unclear, however a slim disc is expected due to the high radiation pressure at these accretion rates, and is entirely consistent with our observations., Comment: 23 pages, 16 figures, 10 tables. Accepted for publication in MNRAS

The effects of the interaction between Type Ia supernova ejecta and their circumstellar wind on the photometric properties of Type Ia supernovae are investigated. We assume that a hydrogen-rich, dense, and extended circumstellar matter (CSM) is formed by the steady mass loss of their progenitor systems. The CSM density is assumed to be proportional to r^{-2}. When the mass-loss rate is above 1e-4 Msun/yr with a wind velocity of 100 km/s, CSM interaction results in an early flux excess in optical light-curves within 4 days of explosion. In these cases, the optical colour quickly evolves to the blue. The ultraviolet flux below 3000 A is found to have a persistent flux excess compared to Type Ia supernovae as long as CSM interaction continues. Type Ia supernovae with progenitor mass-loss rates between 1e-4 and 1e-3 Msun/yr may not have a CSM that is dense enough to affect spectra to make them Type Ia-CSM, but they may still result in Type Ia supernovae with an early optical flux excess. Because they have a persistent ultraviolet flux excess, ultraviolet light curves around the luminosity peak would be significantly different from those with a low-density CSM., Comment: 8 pages, 8 figures, accepted by Monthly Notices of the Royal Astronomical Society

We study the long-term (t >> 10 s) evolution of the accretion disk after a neutron star(NS)-NS or NS-black hole merger, taking into account the radioactive heating by r-process nuclei formed in the first few seconds. We find that the cumulative heating eventually exceeds the disk's binding energy at t ~ 10^2 s (\alpha/0.1)^{-1.8} (M/2.6 Msun)^{1.8} after the merger, where \alpha is the Shakura-Sunyaev viscosity parameter and M is the mass of the remnant object. This causes the disk to evaporate rapidly and the jet power to shut off. We propose that this is the cause of the steep flux decline at the end of the extended emission (EE) or X-ray plateau seen in many short gamma-ray bursts (GRBs). The shallow flux evolution before the steep decline is consistent with a plausible scenario where the jet power scales linearly with the disk mass. We suggest that the jets from NS mergers have two components -- a short-duration narrow one corresponding to the prompt gamma-ray emission and a long-lasting wide component producing the EE. This leads to a prediction that "orphan EE" (without the prompt gamma-rays) may be a promising electromagnetic counterpart for NS mergers observable by future wide-field X-ray surveys. The long-lived disk produces a slow ejecta component that can efficiently thermalize the energy carried by beta-decay electrons up to t ~ 100 d and contributes 10% of the kilonova's bolometric luminosity at these late epochs. We predict that future ground-based and JWST near-IR spectroscopy of nearby (< 100 Mpc) NS mergers will detect narrow (~0.01c) line features a few weeks after the merger, which provides a powerful probe of the atomic species formed in these events., Comment: 13 pages, 6 figures, plus appendix. MNRAS submitted. Comments welcome!

A striking feature of the Vela pulsar (PSR J0835$-$4510) is that it undergoes sudden increases in its spin frequency, known as glitches, with a fractional amplitude on the order of $10^{-6}$ approximately every 900 days. Glitches of smaller magnitudes are also known to occur in Vela. Their distribution in both time and amplitude is less well constrained but equally important for understanding the physical process underpinning these events. In order to better understand these small glitches in Vela, an analysis of high-cadence observations from the Mount Pleasant Observatory is presented. A hidden Markov model (HMM) is used to search for small, previously undetected glitches across 24 years of observations covering MJD 44929 to MJD 53647. One previously unknown glitch is detected around MJD 48636 (Jan 15 1992), with fractional frequency jump $\Delta f/f = (8.19 \pm 0.04) \times 10^{-10}$ and frequency derivative jump $\Delta\dot{f}/\dot{f} = (2.98 \pm 0.01) \times 10^{-4}$. Two previously reported small glitches are also confidently re-detected, and independent estimates of their parameters are reported. Excluding these events, 90% confidence frequentist upper limits on the sizes of missed glitches are also set, with a median upper limit of $\Delta f^{90\%}/f = 1.35 \times 10^{-9}$. Upper limits of this kind are enabled by the semi-automated and computationally efficient nature of the HMM, and are crucial to informing studies which are sensitive to the lower end of the glitch size distribution., Comment: 10 pages, 7 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

We report a detailed analysis of the XMM-Newton spectra of six Narrow-Line Seyfert 1 (NLS1) galaxies at redshift z = 0.35--0.92. Compared with the NLS1s at lower redshift in the previously most-studied sample, these NLS1s have larger black hole (BH) masses ($\log\,M_\text{BH}>7.5$) with similar or even lower Eddington ratios. Our extended XMM-Newton sample of NLS1s shows strong soft X-ray excess emission below 2 keV. The quantified soft excess strength does not show an obvious discrepancy from previous studies of the lower-redshift NLS1s. The systematic effect in the measurement of the Eddington ratio mainly lies in the bolometric correction factor. We also tentatively fit the spectra assuming two more physical models for the soft excess: warm Comptonization and relativistic reflection from the inner accretion disk. In the first scenario, we confirm the ubiquity of a warm and optically thick corona. The behavior of a single source can be better explained by relativistic reflection, although we cannot distinguish which model is a more favorable explanation for the soft excess based on the best-fit statistics., 13 pages, 10 figures

Fast radio bursts (FRBs) are mysterious astronomical phenomena, and it is still uncertain whether they consist of multiple types. In this study we use two nonlinear dimensionality reduction algorithms - Uniform Manifold Approximation and Projection (UMAP) and t-distributed stochastic neighbour embedding (t-SNE) - to differentiate repeaters from apparently non-repeaters in FRBs. Based on the first Canadian Hydrogen Intensity Mapping Experiment (CHIME) FRB catalogue, these two methods are applied to standardized parameter data and image data from a sample of 594 sub-bursts and 535 FRBs, respectively. Both methods are able to differentiate repeaters from apparently non-repeaters. The UMAP algorithm using image data produces more accurate results and is a more model-independent method. Our result shows that in general repeater clusters tend to be narrowband, which implies a difference in burst morphology between repeaters and apparently non-repeaters. We also compared our UMAP predictions with the CHIME/FRB discovery of 6 new repeaters, the performance was generally good except for one outlier. Finally, we highlight the need for a larger and more complete sample of FRBs., 10 pages, 9 figures, 4 tables, accepted by MNRAS

An open question in SN Ia research is where the boundary lies between 'normal' Type Ia supernovae (SNe Ia) that are used in cosmological measurements and those that sit off the Phillips relation. We present the spectroscopic modelling of one such '86G-like' transitional SN Ia, SN 2021rhu, that has recently been employed as a local Hubble Constant calibrator using a tip of the red-giant branch measurement. We detail its modelling from -12 d until maximum brightness using the radiative-transfer spectral-synthesis code tardis. We base our modelling on literature delayed-detonation and deflagration models of Chandrasekhar mass white dwarfs, as well as the double-detonation models of sub-Chandrasekhar mass white dwarfs. We present a new method for 'projecting' abundance profiles to different density profiles for ease of computation. Due to the small velocity extent and low outer densities of the W7 profile, we find it inadequate to reproduce the evolution of SN 2021rhu as it fails to match the high-velocity calcium components. The host extinction of SN 2021rhu is uncertain but we use modelling with and without an extinction correction to set lower and upper limits on the abundances of individual species. Comparing these limits to literature models we conclude that the spectral evolution of SN 2021rhu is also incompatible with double-detonation scenarios, lying more in line with those resulting from the delayed detonation mechanism (although there are some discrepancies, in particular a larger titanium abundance in SN 2021rhu compared to the literature). This suggests that SN 2021rhu is likely a lower luminosity, and hence lower temperature, version of a normal SN Ia., 25 pages, 22 figures, accepted for publication in MNRAS

Natal supernova kicks, the linear momentum compact remnants receive during their formation, are an essential part of binary population synthesis (BPS) models. Although these kicks are well-supported by evidence, their underlying distributions and incorporation into BPS models is uncertain. In this work, we investigate the nature of natal kicks using a previously proposed analytical prescription where the strength of the kick is given by $v_\text{k}=\alpha\frac{m_\text{ejecta}}{m_\text{remnant}}+\beta~\text{km s}^{-1}$, for free parameters $\alpha$ and $\beta$. We vary the free parameters over large ranges of possible values, comparing these synthetic populations simultaneously against four constraints: the merger rate of compact binary neutron star (BNS) systems, the period-eccentricity distribution of galactic BNSs, the velocity distribution of single-star pulsars, and the likelihood for low-ejecta mass supernovae to produce low-velocity kicks. We find that different samples of the parameter space satisfy each tests, and only 1 per cent of the models satisfy all four constraints simultaneously. Although we cannot identify a single best kick model, we report $\alpha=115^{+40}_{-55}~\text{km s}^{-1}, \beta=15^{+10}_{-15}~\text{km s}^{-1}$ as the center of the region of the parameter space that fulfils all of our constraints, and expect $\beta\geq0~\text{km s}^{-1}$ as a further constraint. We also suggest further observations that will enable future refinement of the kick model. A sensitive test for the kick model will be the redshift evolution of the BNS merger rate since this is effectively a direct measure of the delay-time distribution for mergers. For our best fitting values, we find that the peak of the BNS merger rate is the present-day., Comment: Accepted by MNRAS

We report detailed X-ray observations of the unique binary system $\epsilon$ Lupi, the only known short-period binary consisting of two magnetic early-type stars. The components have comparably strong, but anti-aligned magnetic fields. The orbital and magnetic properties of the system imply that the magnetospheres overlap at all orbital phases, suggesting the possibility of variable inter-star magnetospheric interaction due to the non-negligible eccentricity of the orbit. To investigate this effect, we observed the X-ray emission from $\epsilon$ Lupi both near and away from periastron passage, using the Neutron Star Interior Composition Explorer mission (NICER) X-ray Telescope. We find that the system produces excess X-ray emission at the periastron phase, suggesting the presence of variable inter-star magnetospheric interaction. We also discover that the enhancement at periastron is confined to a very narrow orbital phase range ($\approx 5\%$ of the orbital period), but the X-ray properties close to periastron phase are similar to those observed away from periastron. From these observations, we infer that the underlying cause is magnetic reconnection heating the stellar wind plasma, rather than shocks produced by wind-wind collision. Finally, by comparing the behavior of $\epsilon$ Lupi with that observed for cooler magnetic binary systems, we propose that elevated X-ray flux at periastron phase is likely a general characteristic of interacting magnetospheres irrespective of the spectral types of the constituent stars., Comment: Accepted for publication in MNRAS (20 pages, 17 figures)

Stars and stellar remnants orbiting a supermassive black hole (SMBH) can interact with an active galactic nucleus (AGN) disc. Over time, prograde orbiters (inclination $ii_{\rm ret}$ sBH decay towards embedded retrograde orbits $(i \to 180^{\circ})$. sBH near polar orbits $(i \sim 90^{\circ})$ and stars on nearly embedded retrograde orbits $(i \sim 180^{\circ})$ show the greatest decreases in $a$. Whether a star is captured by the disc within an AGN lifetime depends primarily on disc density, and secondarily on stellar type and initial $a$. For sBH, disc capture-time is longest for polar orbits, low mass sBH and lower density discs. Larger mass sBH should typically spend more time in AGN discs, with implications for the embedded sBH spin distribution., 9 pages, 7 figures, 1 table

Double white dwarfs (DWDs) will be the most numerous gravitational-wave (GW) sources for the Laser Interferometer Space Antenna (LISA). Most of the Galactic DWDs will be unresolved and will superpose to form a confusion noise foreground, the dominant LISA noise source around $\sim 0.5\mathrm{-}3\,\mathrm{mHz}$. A small fraction of these sources will stand out from the background and be individually detectable. Uniquely among GW sources, a handful of these binaries will be known in advance from electromagnetic (EM) observations and will be guaranteed sources of detectable GWs in the LISA band; these are known as verification binaries (VBs). High-cadence photometric surveys are continuously discovering new VB systems, and their number will continue to grow ahead of the launch of LISA. We analyse, in a fully Bayesian framework, all the currently known VB candidates with the latest design requirements for the LISA mission and find that 25 of the considered sources can be detected within a $4\,\mathrm{yr}$ observation time. We explore what can be expected from GW observations, both alone and in combination with EM observations, and estimate the VB's time to detection in the early months of LISA operations. We also show how VBs can be analysed in the case where their GW signals compete with many other unknown binary signals (both resolved and unresolved) from a realistic Galactic population of DWDs., Comment: 13 pages, 8 figures, plus appendices (accepted in Monthly Notices of the Royal Astronomical Society on 26 April 2023)

MAXI J1348-630 is a low-mass X-ray black hole binary located in the Galaxy and undergone the X-ray outburst in 2019. We analyzed the observation data in very soft state during the outburst between MJD 58588 and MJD 58596 based on the Insight-HXMT observations from 2 -- 20 keV via the continuum fitting method to measure the spin of the stellar-mass black hole in MAXI J1348-630. The inner disk temperature and the apparent inner disk radius were found to be $0.47\pm 0.01 \rm keV$ and $5.33\pm 0.10 \ R_{g}$ from the observation data modeled by the multicolor disc blackbody model. Assuming the distance of the source $D\sim 3.4 \rm kpc$, the mass of the black hole $M\sim 11 \ M_{\odot}$, and the inclination of the system $i\sim 29.2^{\circ}$, the spin is determined to be $a_{\star}=0.41\pm 0.03$ for fixing hardening factor at 1.6 and $n_{H}=8.6\times 10^{21} \rm cm^{-2}$. Besides, considering the uncertainty of the parameters $D, M, i$ of this system, with the Monte Carlo analysis, we still confirm the moderate spin of the black hole as $a_{\star}=0.42^{+0.13}_{-0.50}$. Some spectral parameters (e.g., column density and hardening factor) which could affect the measurements of the BH spin are also briefly discussed., Comment: 10 pages, 14 figures, 5 tables, accept for publication in MNRAS

Quasiperiodic eruptions (QPEs) from low-mass galaxy centres may result from accretion from a white dwarf in a very eccentric orbit about the central massive black hole. Evolution under gravitational radiation losses reduces the separation and eccentricity. I note that below a critical eccentricity $e_{\rm crit} \simeq 0.97$, the accretion disc's viscous timescale at pericentre passage is probably longer than the orbital period $P$, and periodic eruptive behaviour is no longer possible. These QPE descendant systems (QPEDs) are then likely to produce quasiperiodic oscillations (QPOs) rather than eruptions, varying more smoothly over the orbital cycle, with duty cycles $\sim1$. I identify 2XMM J123103.2+110648 ($P \simeq 3.9$~hr) and (more tentatively) RE J1034+396 ($P \simeq 1$~hr) as candidate systems of this type, and find agreement with their deduced eccentricities $e < e_{\rm crit}$. The absence of eruptions and the lower accretion luminosities resulting from the smaller gravitational radiation losses may make QPED systems harder to discover. Ultimately they must evolve to have viscous times much longer than the orbital period, and either remain steady, or possibly have infrequent but large outbursts. The latter systems would be massive analogues of the soft X-ray transients produced by low stellar-mass X-ray binaries., accepted for publication in MNRAS

Measuring the proper motion of the emission component in radio-quiet quasars (RQQs) could help to distinguish between the origins of the radio emission and to understand whether the jet production mechanism is the same in radio-loud quasars (RLQs) and RQQs. PG 1351+640 is one of the few RQQs suitable for proper motion studies: it has two compact components on milli-arcsecond scales, a flat-spectrum core and a steep-spectrum jet; both components are >2 mJy at 5 GHz and are well suited for Very Long Baseline Array (VLBA) observations. We compare recent VLBA observations with that made seventeen years ago and find no significant change in the core-jet separation between 2005 and 2015 (a proper motion of 0.003 mas yr-1). However, the core-jet separation increased significantly between 2015 and 2022, inferring a jet proper motion velocity of 0.063 mas yr-1, which corresponds to an apparent transverse velocity of 0.37c. The result suggests that the jet of the RQQ PG 1351+640 is mildly relativistic and oriented at a relatively small viewing angle., The article has been published by Oxford University Press: https://academic.oup.com/mnrasl/advance-article-abstract/doi/10.1093/mnrasl/slad051/7146829?utm_source=advanceaccess&utm_campaign=mnrasl&utm_medium=email

A comparative analysis of the individual bursts between FRB 20190520B and FRB 20121102A is presented by compiling a sample of bursts in multiple wavelengths. It is found that the peak frequency ($\nu_p$) distribution of the bursts of FRB 20190520B illustrates four discrete peaks in $\sim1-6$ GHz and their spectral width distribution can be fitted with a log-normal function peaking at 0.35 GHz. The discrete $\nu_p$ distribution and the narrow-banded spectral feature are analogous to FRB 20121102A. The burst duration of FRB 20190520B in the rest frame averages 10.72 ms, longer than that of FRB 20121102A by a factor 3. The specific energy ($E_{\rm \mu_{\rm c}}$) at 1.25 GHz of FRB 20190520B observed with the FAST telescope narrowly ranges in $[0.4, 1]\times 10^{38}$ erg, different from the bimodal $E_{\rm \mu_{\rm c}}$ distribution of FRB 20121102A. Assuming a Gaussian spectral profile of the bursts, our Monte Carlo simulation analysis suggests that a power-law (PL) or a cutoff power-law (CPL) energy function can comparably reproduce the $E_{\rm \mu_{\rm c}}$ distribution of FRB 20190520B. The derived energy function index of the PL model is $4.46\pm 0.17$, much steeper than that of FRB 20121102A ($1.82^{+0.10}_{-0.30}$). For the CPL model, we obtain an index of $0.47$ and a cutoff energy of $7.4\times 10^{37}$ erg. Regarding the predicted $\nu_p$ distribution in 1-2 GHz, the CPL model is more preferred than the PL model. These results indicate that FRB 20190520B and FRB 20121102A shares similar spectral properties, but their energy functions are intrinsically different., Comment: 2023,MNRAS in press, 12 pages, 7 figures

The nearby radio galaxy M87 is a prime target for studying black hole accretion and jet formation^{1,2}. Event Horizon Telescope observations of M87 in 2017, at a wavelength of 1.3 mm, revealed a ring-like structure, which was interpreted as gravitationally lensed emission around a central black hole^3. Here we report images of M87 obtained in 2018, at a wavelength of 3.5 mm, showing that the compact radio core is spatially resolved. High-resolution imaging shows a ring-like structure of 8.4_{-1.1}^{+0.5} Schwarzschild radii in diameter, approximately 50% larger than that seen at 1.3 mm. The outer edge at 3.5 mm is also larger than that at 1.3 mm. This larger and thicker ring indicates a substantial contribution from the accretion flow with absorption effects in addition to the gravitationally lensed ring-like emission. The images show that the edge-brightened jet connects to the accretion flow of the black hole. Close to the black hole, the emission profile of the jet-launching region is wider than the expected profile of a black-hole-driven jet, suggesting the possible presence of a wind associated with the accretion flow., 50 pages, 18 figures, 3 tables, author's version of the paper published in Nature

We present ultraviolet (UV) to near-infrared (NIR) observations and analysis of the nearby Type Ia supernova SN 2021fxy. Our observations include UV photometry from Swift/UVOT, UV spectroscopy from HST/STIS, and high-cadence optical photometry with the Swope 1-m telescope capturing intra-night rises during the early light curve. Early $B-V$ colours show SN 2021fxy is the first "shallow-silicon" (SS) SN Ia to follow a red-to-blue evolution, compared to other SS objects which show blue colours from the earliest observations. Comparisons to other spectroscopically normal SNe Ia with HST UV spectra reveal SN 2021fxy is one of several SNe Ia with flux suppression in the mid-UV. These SNe also show blue-shifted mid-UV spectral features and strong high-velocity Ca II features. One possible origin of this mid-UV suppression is the increased effective opacity in the UV due to increased line blanketing from high velocity material, but differences in the explosion mechanism cannot be ruled out. Among SNe Ia with mid-UV suppression, SNe 2021fxy and 2017erp show substantial similarities in their optical properties despite belonging to different Branch subgroups, and UV flux differences of the same order as those found between SNe 2011fe and 2011by. Differential comparisons to multiple sets of synthetic SN Ia UV spectra reveal this UV flux difference likely originates from a luminosity difference between SNe 2021fxy and 2017erp, and not differing progenitor metallicities as suggested for SNe 2011by and 2011fe. These comparisons illustrate the complicated nature of UV spectral formation, and the need for more UV spectra to determine the physical source of SNe Ia UV diversity., 26 pages, 19 figures, 9 tables; submitted to MNRAS, posted after receiving referee comments

The transient X-ray source MAXI J1810-222 was discovered in 2018 and has been active ever since. A long combined radio and X-ray monitoring campaign was performed with ATCA and Swift respectively. It has been proposed that MAXI J1810-222 is a relatively distant black hole X-ray binary, albeit showing a very peculiar outburst behaviour. Here, we report on the spectral study of this source making use of a large sample of NICER observations performed between 2019 February and 2020 September. We detected a strong spectral absorption feature at $\sim$1 keV, which we have characterised with a physical photoionisation model. Via a deep scan of the parameters space, we obtained evidence for a spectral-state dependent outflow, with mildly relativistic speeds. In particular, the soft and intermediate states point to a hot plasma outflowing at 0.05-0.15 $c$. This speeds rule-out thermal winds and, hence, they suggest that such outflows could be radiation pressure or (most likely) magnetically-driven winds. Our results are crucial to test current theoretical models of wind formation in X-ray binaries., 6 pages, 10 figures, accepted for publication in MNRAS Letter

We present polarization observations of the young supernova remnant (SNR) Cas A using the High-resolution Airborne Wideband Camera-Plus (HAWC+) instrument onboard the Stratospheric Observatory for Infrared Astronomy (SOFIA). The polarization map at 154 microns reveals dust grains with strong polarization fractions (5 - 30 percent), supporting previous measurements made over a smaller region of the remnant at 850 microns. The 154 microns emission and the polarization signal is coincident with a region of cold dust observed in the southeastern shell and in the unshocked central ejecta. The highly polarized far-IR emission implies the grains are large (greater than 0.14 microns) and silicate-dominated. The polarization level varies across the SNR, with an inverse correlation between the polarization degree and the intensity and smaller polarization angle dispersion for brighter SNR emission. Stronger polarization is detected between the bright structures. This may result from a higher collision rate between the gas and dust producing a lower grain alignment efficiency where the gas density is higher. We use the dust emission to provide an estimate of the magnetic field strength in Cas A using the Davis-Chandrasekhar-Fermi method. The high polarization level is direct evidence that grains are highly elongated and strongly aligned with the magnetic field of the SNR. The dust mass from the polarized region is 0.14+-0.04 Msun, a lower limit of the amount of dust present within the ejecta of Cas A. This result strengthens the hypothesis that core-collapse SNe are an important contributor to the dust mass in high redshift galaxies., Comment: MNRAS, accepted (18 pages with 14 figures)

With tantalizing evidence of the recent e-Rosita mission, re-discovering very soft X-rays and EUV radiation from a cluster of galaxies or its environment, the question of the origin of cluster EUV excess is revisited in this work. It will be shown that the gas temperature, density, and frozen-in magnetic field of the intracluster medium, collectively support the emission and propagation of coherent \u{C}erenkov radiation, which is low frequency and large amplitude radiation capable of accelerating charged particles to relativistic speeds. Owing to the spectrum of \u{C}erenkov radiation, most of the incipient relativistic electrons undergo inverse-Compton scattering with the cosmic microwave background. It turns out the scattered radiation has observable ramifications only in the EUV band, of photon energy $70 -- 100$~eV, having a luminosity $\approx 10^{44}$~ergs~s$^{-1}$. This luminosity is on par with the EUV excess level detected from Abell 1795 and the Coma cluster. It should be stressed, as {\it caveat emptor}, that although the main subject is the putative large amplitude coherent \u{C}erenkov modes which are highly nonlinear, the results presented were derived using a quasi-linear approach to highlight the observable features of the phenomenon, namely the EUV emission., Comment: MNRAS in press

AT 2022cmc is a luminous optical transient ($\nu L_{\nu} \gtrsim 10^{45}$ erg s$^{-1}$) accompanied by decaying non-thermal X-rays (peak duration $t_{\rm X} \lesssim$ days and isotropic energy $E_{\rm X,iso} \gtrsim 10^{53}$ erg) and a long-lived radio/mm synchrotron afterglow, which has been interpreted as a jetted tidal disruption event (TDE). Both an equipartition analysis and a detailed afterglow model reveals the radio/mm emitting plasma to be expanding mildly relativistically (Lorentz factor $\Gamma \gtrsim\,few$) with an opening angle $\theta_{\rm j}\simeq0.1$ and roughly fixed energy $E_{\rm j,iso} \gtrsim few \times 10^{53}$ erg into an external medium of density profile $n \propto R^{-k}$ with $k \simeq 1.5-2$, broadly similar to that of the first jetted TDE candidate Swift J1644+57 and consistent with Bondi accretion at a rate $\sim 10^{-3}\dot{M}_{\rm Edd}$ onto a $10^{6}M_{\odot}$ black hole before the outburst. The rapidly decaying optical emission over the first days is consistent with fast-cooling synchrotron radiation from the same forward shock as the radio/mm emission, while the bluer slowly decaying phase to follow likely represents a separate thermal emission component. Emission from the reverse shock may have peaked during the first days, but whose non-detection in the optical band places an upper bound $\Gamma_{\rm j} \lesssim 100$ on the Lorentz factor of the unshocked jet. Although a TDE origin for AT 2022cmc is indeed supported by some observations, the vast difference between the short-lived jet activity phase $t_{\rm X} \lesssim$ days relative to the months-long thermal optical emission, also challenges this scenario. A stellar core-collapse event giving birth to a magnetar or black hole engine of peak duration $\sim 1$ day offers an alternative model also consistent with the circumburst environment, if interpreted as a massive-star wind., Comment: 11 pages, 9 figures, 2 tables, accepted for publication in MNRAS

In this work, we report the presence of rapid intra-night optical variations in both -- flux and polarization of the blazar BL Lacertae during its unprecedented 2020--2021 high state of brightness. The object showed significant flux variability and some color changes, but no firmly detectable time delays between the optical bands. The linear polarization was also highly variable in both -- polarization degree and angle (EVPA). The object was observed from several observatories throughout the world, covering in a total of almost 300 hours during 66 nights. Based on our results, we suggest, that the changing Doppler factor of an ensemble of independent emitting regions, travelling along a curved jet that at some point happens to be closely aligned with the line of sight can successfully reproduce our observations during this outburst. This is one of the most extensive variability studies of the optical polarization of a blazar on intra-night timescales., 23 pages,7 figures, 5 Tables (2 as appendix). Accepted for publication in MNRAS

We report on the case of optical periodic modulation discovered in two Active Galactic Nuclei (AGN) and one candidate AGN. Analyzing the archival optical data obtained from large transient surveys, namely the Catalina Real-Transient Survey (CRTS) and the Zwicky Transient Facility (ZTF), we find periodicities of 2169.7, 2103.1, and 1462.6\,day in sources J0122+1032, J1007+1248 (or PG~1004+1248), and J2131$-$1127, respectively. The optical spectra of the first two indicate that the first is likely a blazar and the second a type 1 Seyfert galaxy, and while no spectroscopic information is available for the third one, its overall properties suggest that it is likely an AGN. In addition, mid-infrared (MIR) light curve data of the three sources, taken by the Wide-field InfraredSurvey Explorer (WISE), are also analyzed. The light curves show significant variations, but not appearing related to the optical periodicities. Based on the widely-discussed super-massive black-hole binary (SMBHB) scenario, we discuss the origin of the optical modulation. Two possible interesting features, an additional 162-day short optical periodicity in J2131$-$1127 and the consistency the X-ray flux variations of J1007+1248 with its optical periodicity, are also discussed within the SMBHB scenario., Comment: 10 pages, 6 figures, 5 tables, accepted for publication in MNRAS

We present the to-date largest parameter space exploration of binaries in circumbinary disks (CBDs), deriving orbital evolution prescriptions for eccentric, unequal mass binaries from our suite of hydrodynamic simulations. In all cases, binary eccentricities evolve towards steady state values that increase with mass ratio, and saturate at an equilibrium eccentricity $e_{\rm b, eq} \sim 0.5$ in the large mass ratio regime, in line with resonant theory. For binaries accreting at their combined Eddington limit, a steady state eccentricity can be achieved within a few Megayears. Once at their steady state eccentricities, binaries with $q_{\rm b} \gtrsim 0.3$ evolve towards coalescence, while lower mass ratio systems expand due to CBD torques. We discuss implications for population studies of massive black hole binaries, protostars in binary systems, and post-common envelope binaries observed by ground-based gravitational wave detectors., 12 pages, 8 figures. Submitted to MNRAS

We have found evidence of bulk velocities following active galactic nucleus (AGN) bubbles in the Virgo cluster and galaxy motions in the Centaurus cluster. In order to increase the sample and improve our understanding of the intracluster medium (ICM), we present the results of a detailed mapping of the Ophiuchus cluster with {\it XMM-Newton} to measure bulk flows through very accurate Fe~K measurements. To measure the gas velocities we use a novel EPIC-pn energy scale calibration, which uses the Cu K$\alpha$ instrumental line as reference for the line emission. We created 2D spectral maps for the velocity, metallicity, temperature, density, entropy and pressure with a spatial resolution of 0.25$'$ ($\sim 26$~kpc). The ICM velocities in the central regions where AGN feedback is most important are similar to the velocity of the brightest cluster galaxy (BCG). We have found a large interface region where the velocity changes abruptly from blueshifted to redshifted gas which follows a sharp surface brightness discontinuity. We also found that the metallicities and temperatures do not change as we move outwards from the giant radio fossil previously identified in radio observations of the cluster. Finally, we have found a contribution from the kinetic component of $, Comment: 15 pages, 20 figures. This paper is part of a series on the ICM velocity structure using XMM-Newton observations. Related series papers: arXiv:2109.06213, arXiv:2203.12635

The nonlinear interaction between electromagnetic waves and plasmas attracts significant attention in astrophysics because it can affect the propagation of Fast Radio Bursts (FRBs) -- luminous millisecond-duration pulses detected at radio frequency. The filamentation instability (FI) -- a type of nonlinear wave-plasma interaction -- is considered to be dominant near FRB sources, and its nonlinear development may also affect the inferred dispersion measure of FRBs. In this paper, we carry out fully kinetic particle-in-cell simulations of the FI in unmagnetized pair plasmas. Our simulations show that the FI generates transverse density filaments, and that the electromagnetic wave propagates in near vacuum between them, as in a waveguide. The density filaments keep merging until force balance between the wave ponderomotive force and the plasma pressure gradient is established. We estimate the merging timescale and discuss the implications of filament merging for FRB observations., Submitted to MNRAS

High mass accretion rate onto strongly magnetised neutron stars results in the appearance of accretion columns supported by the radiation pressure and confined by the strong magnetic field of a star. At mass accretion rates above $\sim 10^{19}\,{\rm g\,s^{-1}}$, accretion columns are expected to be advective. Under such conditions, a noticeable part of the total energy release can be carried away by neutrinos of a MeV energy range. Relying on a simple model of the neutrino luminosity of accreting strongly magnetised neutron stars, we estimate the neutrino energy fluxes expected from six ULX pulsars known up to date and three brightest Be X-ray transits hosting magnetised neutron stars. Despite the large neutrino luminosity expected in ULX pulsars, the neutrino energy flux from the Be X-ray transients of our Galaxy, SMC and LMC is dominant. However, the neutrino flux from the brightest X-ray transients is estimated to be below the isotropic background by two orders of magnitude at least, which makes impossible direct registration of neutrino emission from accreting strongly magnetised neutron stars nowadays., Comment: 7 pages, 6 figures, submitted to MNRAS

The ground-based gravitational wave (GW) observatories discover a population of merging stellar binary black holes (BBHs), which are promising targets for multiband observations by the low-, middle-, and high-frequency GW detectors. In this paper, we investigate the multiband GW detections of BBHs and demonstrate the advantages of such observations in improving the localization and parameter estimates of the sources. We generate mock samples of BBHs by considering different formation models as well as the merger rate density constrained by the current observations (GWTC-3). We specifically consider the astrodynamical middle-frequency interferometer GW observatory (AMIGO) in the middle-frequency band and estimate that it may detect $21$-$91$ BBHs with signal-to-noise ratio $\varrho\geq8$ in a $4$-yr observation period. The multiband observations by the low-frequency detectors [Laser Interferometer Space Antenna (LISA) and Taiji] and AMIGO may detect $5$-$33$ BBHs with $\varrho_{\rm LT}\geq5$ and $\varrho_{\rm AMI}\geq5$, which can evolve to the high-frequency band within $4$ yr and can be detected by the Cosmic Explorer (CE) and Einstein Telescope (ET). The joint observations of LISA-Taiji-AMIGO-ET-CE may localize the majority of the detectable BBHs in sky areas of $7\times10^{-7}$ to $2\times10^{-3}$ deg$^2$, which is improved by a factor of $\sim120$, $\sim2.4\times10^{5}$, $\sim1.8\times10^{4}$, or $\sim1.2\times10^{4}$, comparing with those by only adopting CE-ET, AMIGO, LISA-Taiji, or LISA-Taiji-AMIGO. These joint observations can also lead to an improvement of the measurement precision of the chirp mass (symmetric mass ratio) by a factor of $\sim5.5\times10^{4}$ ($33$), $\sim16$ ($8$), $\sim120$ ($90$), or $\sim5$ ($5$), comparing with those by CE-ET, AMIGO, LISA-Taiji, or LISA-Taiji-AMIGO., 16 pages, 12 figures, MNRAS published

Hydrodynamic studies of stellar-mass compact objects (COs) in a common envelope (CE)have shown that the accretion rate onto the CO is a few orders of magnitude below the Bondi-Hoyle-Lyttleton (BHL) estimate. This is several orders of magnitude above the Eddington limit and above the limit for neutrino-cooled accretion (i.e., hypercritical accretion, or HCA). Considering that a binary system inside the CE of a third star accretes material at nearly the same rate as a single object of the same total mass, we propose stellar-evolution channels which form binary black hole (BBH) systems with its component masses within the pair-instability supernova (PISN) mass gap. Our model is based on HCA onto the BBH system engulfed into the CE of a massive tertiary star. Furthermore, we propose a mass transfer mode which allows to store mass lost by the binary onto a third star. Through the use of population synthesis simulations for the evolution of BBHs and standard binary-evolution principles for the interaction with a tertiary star, we are able to produce BBHs masses consistent with those estimated for GW190521. We also discuss the massive binary system Mk34 as a possible progenitor of BBHs in the PISN gap, as well as the spin distribution of the observed mergers in the gravitational-wave catalog., 10 pages, 3 figures, 2 tables, submitted. Comments welcome

We survey the spectro-temporal properties of fast radio bursts from FRB 20121102A observed by earlier studies across a wide range of frequencies. We investigate 167 bursts from FRB 20121102A spanning frequencies 1--7.5GHz, durations of less than 1 ms to approximately 10 ms, with low and high energies, and with different wait-times. We find from this sample of bursts a strong agreement with the inverse relationship between sub-burst slope and duration and with other predictions made by the triggered relativistic dynamical model (TRDM). Earlier results found agreement with those predictions across three different repeating FRB sources. For this sample of bursts, we find that the sub-burst slope as well as the 'sad trombone' drift rate are consistent with being in a quadratic relationship with frequency and that both these quantities are inversely proportional to the duration. We also find that the duration decreases with increasing frequency as well as a statistically significant correlation between the sub-burst duration and bandwidth (proportional to $t^{-1/2}$) that is unexpected. No distinct group of bursts in this sample deviated from these relationships, however significant scatter can be seen in measurements. This study demonstrates the consistent existence of relationships between the spectro-temporal properties of bursts from a repeating FRB source. A simple explanation for the inverse relation between the sub-burst slope and duration is an inherently narrowband emission process. We make all measurements available as well as a graphical user interface called Frbgui developed and used to perform measurements of burst waterfalls., Comment: Accepted at MNRAS. 13 pages, 7 figures

We examine the time delay between radio and X-ray and between narrow radio frequency flares in Sagittarius A* (Sgr A*), from analyses of the synchrotron, bremsstrahlung and monochromatic luminosity curves. Using the results of 2D relativistic radiation magnetohydrodynamic (MHD) simulations based on the shock oscillation model, we find three types of time delay between the synchrotron and bremsstrahlung emissions: Type A with a time delay of 2 -- 3 h on the shock descending branch, Type B with no time delay and Type C with an inverse time delay of 0.5 -- 1 h on the shock ascending branch. The time delays in Types A and C are interpreted as a transit time of Alfv\'{e}n and acoustic waves between both emission dominant regions, respectively. The delay times between 22 and 43 GHz flares and between 8 and 10 GHz flares are $\sim$ 13 -- 26 min and 13 min, respectively, while the inverse delay also occurs dependently on the shock location branch. These time delays between the narrow radio bands are interpreted as the transit time of the acoustic wave between the frequency-dependent effective radii $R_{\tau_{\rm \nu=1}}$, at which the optical depth $\tau_{\rm \nu}$ at the accretion disc surface becomes $\sim$ unity. The shock oscillation model explains well the observed delay times of 0.5 -- 5 h between radio and X-ray, 20 -- 30 min between 22 and 43 GHz and $\sim$ 18 min between 8 and 10 GHz in Sgr A*., Comment: 14 pages, 14 figures, (accepted for publication in MNRAS)

Despite advances in our understanding of low luminosity active galactic nuclei (LLAGNs), the fundamental details about the mechanisms of radiation and flare/outburst in hot accretion flow are still largely missing. We have systematically analyzed the archival Chandra and NuSTAR X-ray data of the nearby LLAGN M81*, whose $L_{\rm bol}\sim 10^{-5} L_{\rm Edd}$. Through a detailed study of X-ray light curve and spectral properties, we find that the X-ray continuum emission of the power-law shape more likely originates from inverse Compton scattering within the hot accretion flow. In contrast to Sgr A*, flares are rare in M81*. Low-amplitude variation can only be observed in soft X-ray band (amplitude usually $\lesssim 2$). Several simple models are tested, including sinusoidal-like and quasi-periodical. Based on a comparison of the dramatic differences of flare properties among Sgr A*, M31* and M81*, we find that, when the differences in both the accretion rate and the black hole mass are considered, the flares in LLAGNs can be understood universally in a magneto-hydrodynamical model., 11 pages, 8 figures, and 4 tables. Accepted to MNRAS

We present observations from X-ray to mid-infrared wavelengths of the most energetic non-quasar transient ever observed, AT2021lwx. Our data show a single optical brightening by a factor $>100$ to a luminosity of $7\times10^{45}$ erg s$^{-1}$, and a total radiated energy of $1.5\times10^{53}$ erg, both greater than any known optical transient. The decline is smooth and exponential and the ultra-violet - optical spectral energy distribution resembles a black body with temperature $1.2\times10^4$ K. Tentative X-ray detections indicate a secondary mode of emission, while a delayed mid-infrared flare points to the presence of dust surrounding the transient. The spectra are similar to recently discovered optical flares in known active galactic nuclei but lack some characteristic features. The lack of emission for the previous seven years is inconsistent with the short-term, stochastic variability observed in quasars, while the extreme luminosity and long timescale of the transient disfavour the disruption of a single solar-mass star. The luminosity could be generated by the disruption of a much more massive star, but the likelihood of such an event occurring is small. A plausible scenario is the accretion of a giant molecular cloud by a dormant black hole of $10^8 - 10^9$ solar masses. AT2021lwx thus represents an extreme extension of the known scenarios of black hole accretion., 11 pages, 5 figures, Accepted for publication in MNRAS

We present the X-ray spectral and temporal analysis of the black hole X-ray transient Swift J1658.2--4242 observed by {\it AstroSat}. Three epochs of data have been analysed using the JeTCAF model to estimate the mass accretion rates and to understand the geometry of the flow. The best-fit disc mass accretion rate ($\dot m_d$) varies between $0.90^{+0.02}_{-0.01}$ to $1.09^{+0.04}_{-0.03}$ $\dot M_{\rm Edd}$ in these observations, while the halo mass accretion rate changes from $0.15^{+0.01}_{-0.01}$ to $0.25^{+0.02}_{-0.01}$ $\dot M_{\rm Edd}$. We estimate the size of the dynamic corona, that varies substantially from $64.9^{+3.9}_{-3.1}$ to $34.5^{+2.0}_{-1.5}$ $r_g$ and a moderately high jet/outflow collimation factor stipulates isotropic outflow. The inferred high disc mass accretion rate and bigger corona size indicate that the source might be in the intermediate to soft spectral state of black hole X-ray binaries. The mass of the black hole estimated from different model combinations is $\sim 14 M_\odot$. In addition, we compute the quasi-periodic oscillation (QPO) frequencies from the model-fitted parameters, which match the observed QPOs. We further calculate the binary parameters of the system from the decay profile of the light curve and the spectral parameters. The estimated orbital period of the system is $4.0\pm0.4$ hr by assuming the companion as a mid or late K-type star. Our analysis using the JeTCAF model sheds light on the physical origin of the spectro-temporal behaviour of the source, and the observed properties are mainly due to the change in both the mass accretion rates and absorbing column density., 10 pages, 5 figures, 3 tables, accepted (06/03/2023) for publication in MNRAS

The spin evolution of main sequence stars has long been of interest for basic stellar evolution, stellar aging, stellar activity, and consequent influence on companion planets. Observations of older than solar late-type main-sequence stars have been interpreted to imply that a change from a dipole-dominated magnetic field to one with more prominent higher multipoles might be necessary to account for the data. The spin-down models that lead to this inference are essentially tuned to the sun. Here we take a different approach which considers individual stars as fixed points rather than just the Sun. We use a time-dependent theoretical model to solve for the spin evolution of low-mass main-sequence stars that includes a Parker-type wind and a time-evolving magnetic field coupled to the spin. Because the wind is exponentially sensitive to the stellar mass over radius and the coronal base temperature, the use of each observed star as a separate fixed point is more appropriate and, in turn, produces a set of solution curves that produces a solution envelope rather than a simple line. This envelope of solution curves, unlike a single line fit, is consistent with the data and does not unambiguously require a modal transition in the magnetic field to explain it., 8 pages, 4 figures; Accepted for publication in MNRAS

We searched for thermal gyro-synchrotron radio emission from a sample of five radio-loud stars whose X-ray coronae contain a hot ($T_e>10^7$ K) thermal component. We used the JVLA to measure Stokes I and V/I spectral energy distributions (SEDs) over the frequency range 15--45 GHz, determining the best-fitting model parameters using power-law and thermal gyro-synchrotron emission models. The SEDs of the three chromospherically active binaries (Algol, UX Arietis, HR 1099) were well-fit by a power-law gyro-synchrotron model, with no evidence for a thermal component. However, the SEDs of the two weak-lined T Tauri stars (V410 Tau, HD 283572) had a circularly polarized enhancement above 30 GHz that was inconsistent with a pure power-law distribution. These spectra were well-fit by summing the emission from an extended coronal volume of power-law gyro-synchrotron emission and a smaller region with thermal plasma and a much stronger magnetic field emitting thermal gyro-synchrotron radiation. We used Bayesian inference to estimate the physical plasma parameters of the emission regions (characteristic size, electron density, temperature, power-law index, and magnetic field strength and direction) using independently measured radio sizes, X-ray luminosities, and magnetic field strengths as priors, where available. The derived parameters were well-constrained but somewhat degenerate. The power-law and thermal volumes in the pre-main-sequence stars are probably not co-spatial, and we speculate they may arise from two distinct regions: a tangled-field magnetosphere where reconnection occurs and a recently discovered axisymmetric toroidal magnetic field, respectively., Submitted to MNRAS; 18 pages, 10 figures

We study the effects of the time-variable properties of thermonuclear X-ray bursts on modeling their millisecond-period burst oscillations. We apply the pulse profile modeling technique that is being used in the analysis of rotation-powered millisecond pulsars by the Neutron Star Interior Composition Explorer (NICER) to infer masses, radii, and geometric parameters of neutron stars. By simulating and analyzing a large set of models, we show that overlooking burst time-scale variability in temperatures and sizes of the hot emitting regions can result in substantial bias in the inferred mass and radius. To adequately infer neutron star properties, it is essential to develop a model for the time variable properties or invest a substantial amount of computational time in segmenting the data into non-varying pieces. We discuss prospects for constraints from proposed future X-ray telescopes., Comment: Accepted for publication in MNRAS

We consider the gas dynamics in an accreting binary system of degenerate stars within the framework of the Newtonian approximation. In such a system, the accretion stream can impact the surface of a white dwarf (WD) or neutron star (NS) as a result of the very compact orbit. This causes a loss of angular momentumfrom the orbit and spin-up of the accretor. We construct approximations for the specific angular momentum of the accreting matter which goes to spin up the accretor and some other parameters of the system. It is shown that the obtained approximation of the specific momentum is qualitatively different from the widely used Keplerian formula. It should affect the boundary between scenarios of immediate tidal disruption and slow mass loss of the donor in WD-WD and NS-NS binaries, as well as the time of stable mass transfer in the stripping scenario., 9 pages, 9 figs, 6 tables. Published in MNRAS

The Search for Extra-Terrestrial Intelligence (SETI) has been conducted for over sixty years, yet no technosignatures have been identified. Previous studies have focused on stars in our galaxy, with few searches in the extragalactic Universe despite a larger volume being available. Civilizations capable of harvesting energy from a star or a galaxy are classified as KII or KIII on the Kardashev scale, respectively. Technosignatures from such advanced civilizations would be extremely luminous and detectable by current radio telescopes, even from distant galaxies. To explore the frontier of extragalactic SETI, we investigate the likely prevalence of extragalactic civilizations possessing a radio transmitter, known as the transmitter rate, based on observational results from the Breakthrough Listen (BL) observations. We calculated the transmitter rate by considering the background galaxies in the field of view of target stars in BL observations. We used a statistical method to derive the total mass of stars in those background galaxies from a galaxy stellar mass function. Our statistical method suggests that less than one in hundreds of trillions of extragalactic civilizations within 969 Mpc possess a radio transmitter above 7.7$\times$10$^{26}$ W of power, assuming one civilization per one-solar-mass stellar system. Additionally, we cross-matched the BL survey fields with the WISE$\times$SuperCOSMOS Photometric Redshift Catalogue and compared with the statistical method. Our result sets the strictest limits to date on the transmitter rate at such high power levels, emphasizing the high efficiency of searching for radio transmitters in galaxies and the rarity of technologically advanced civilizations in our Universe., 5 pages, 1 figure, accepted for publication in MNRAS. A summary video is available: https://www.youtube.com/watch?v=_xaRCnLMjsY&t=21s&ab_channel=NCHUAstronomy

In the paper we develop multi-class classification of Fermi-LAT gamma-ray sources using machine learning with hierarchical determination of classes. One of the main challenges in the multi-class classification of the Fermi-LAT sources is that the size of some of the classes is relatively small, for example with less than 10 associated sources belonging to a class. In the paper we propose a hierarchical structure for the determination of the classes. This enables us to have control over the size of classes and to compare the performance of the classification for different numbers of classes. In particular, the class probabilities in the two-class case can be computed either directly by the two-class classification or by summing probabilities of children classes in multi-class classification. We find that the classifications with few large classes have comparable performance with classifications with many smaller classes. Thus, on the one hand, the few-class classification can be recovered by summing probabilities of classification with more classes while, on the other hand, the classification with many classes gives a more detailed information about the physical nature of the sources. As a result of this work, we construct three probabilistic catalogs. Two catalogs are based on Gaussian mixture model for the determination of the groups of classes: one with the requirement for the minimal number of sources in a group to be larger than 100, which results in six final groups of physical classes, while the other one has the requirement that the minimal number of sources in a group is larger than 15, which results in nine groups of physical classes. The third catalog is based on a random forest determination of the groups of classes with the requirement that the minimal number of sources in a group is larger than 100, which results in six groups of classes. (abridged), Comment: 16 pages, 18 figures, results are available online at https://zenodo.org/record/7538664, v2: minor changes to Figure 1 and related discussion

We report the discovery of millihertz quasi-periodic oscillations (mHz QPOs) from the neutron star (NS) low-mass X-ray binary 4U 1730--22 using the Neutron Star Interior Composition Explorer (NICER). After being inactive for almost 50 years, 4U 1730--22 went into outburst twice between June and August 2021, and between February and July 2022. We analyse all the NICER observations of this source, and detect mHz QPOs with a significance > $4\sigma$ in 35 observations. The QPO frequency of the full data set ranged between ~4.5 and ~8.1 mHz with an average fractional rms amplitude of the order of ~2%. The X-ray colour analysis strongly suggests that 4U 1730--22 was in a soft spectral state during the QPO detections. Our findings are consistent with those reported for other sources where the mHz QPOs have been interpreted as the result of a special mode of He burning on the NS surface called marginally stable nuclear burning (MSNB). We conclude that the mHz QPOs reported in this work are also associated with the MSNB, making 4U 1730--22 the eighth source that shows this phenomenology. We discuss our findings in the context of the heat flux from the NS crust., Comment: 9 pages, 6 figures. Accepted for publication in MNRAS

We analyse all the available Chandra observations of the Cartwheel Galaxy and its compact group, taken between 2001 and 2008, with the main aim of addressing the variability in the X-ray band for this spectacular collisional ring galaxy. We focus on the study of point-like sources, in particular we are interested in Ultraluminous X-ray sources (ULXs, Lx >= 10^39 erg/s), that we treat as a class. We exploit archival XMM-Newton data to enrich the study of the long-term variability, on timescales of months to years. We find a total of 44 sources in the group area, of which 37 in total are ULXs positionally linked with the galaxies and of which we can study variability. They are 29 in the Cartwheel itself, 7 in G1 and 1 in G3. About one third of these 37 sources show long-term variability, while no variability is detected within the single observations. Of those, 5 ULXs have a transient behaviour with a maximum range of variability (Lmax/Lmin) of about one order of magnitude and are the best candidate neutron stars. The X-ray Luminosity Function (XLF) of the point-like sources remains consistent in shape between the Chandra observations both for the Cartwheel galaxy itself and for G1, suggesting that flux variability does not strongly influence the average properties of the population on the observation timescales., 17 pages, 5 figures, Accepted for publication in MNRAS

We propose that the strength of angular momentum transport in accretion discs threaded by net vertical magnetic field is determined by a self-regulation mechanism: the magnetorotational instability (MRI) grows until its own turbulent resistivity damps the fastest growing mode on the scale of the disc thickness. Given weak assumptions as to the structure of MRI-derived turbulence, supported by prior simulation evidence, the proposed mechanism reproduces the known scaling of the viscous $\alpha$-parameter, $\alpha \propto \beta_z^{-1/2}$. Here, $\beta_z = 8\pi p_g/B_{z0}^2$ is the initial plasma $\beta$-parameter on the disc midplane, $B_{z0}$ is the net field, and $p_g $ is the midplane gas pressure. We generalize the argument to discs with strong suprathermal toroidal magnetic fields, where the MRI growth rate is modified from the weak-field limit. Additional sources of turbulence are required if such discs are to become magnetically elevated, with the increased scale heights near the midplane that are seen in simulations. We speculate that tearing modes, associated with current sheets broadened by the effective resistivity, are a possible source of enhanced turbulence in elevated discs., Comment: 9 pages, 2 figures, to be published in Monthly Notices of the Royal Astronomical Society. Replaces arXiv:2211.09261, eliminating section 2.2 and with moderate revisions and clarifications in sections 3.1, 3.3, 4.2 and 4.3

We examine the binding energies of massive stripped-envelope core collapse supernova (SECCSN) progenitors with the stellar evolution code MESA, and find that the jittering jets explosion mechanism is preferred for explosions where carbon-oxygen cores with masses of $>20 M_\odot$ collapse to leave a neutron star (NS) remnant. We calculate the binding energy at core collapse under the assumption that the remnant is a NS. Namely, stellar gas above mass coordinate of $~1.5-2.5 M_\odot$ is ejected in the explosion. We find that the typical binding energy of the ejecta of stripped-envelope progenitors with carbon-oxygen core masses of $M_{CO} > 20 M_\odot$ is $E_{bind}>2 \times 10^{51} erg$. We claim that jets are most likely to explode such cores as jet-driven explosion mechanisms can supply high energies to the explosion. We apply our results to SN 2020qlb, which is a SECCSN with a claimed core mass of $~30-50 M_\odot$, and conclude that the jittering jets explosion mechanism best accounts for such an explosion that leaves a NS., Accepted for publication in MNRAS

The Zwicky Transient Transient Facility (ZTF) performs a systematic neutrino follow-up program, searching for optical counterparts to high-energy neutrinos with dedicated Target-of-Opportunity (ToO) observations. Since first light in March 2018, ZTF has taken prompt observations for 24 high-quality neutrino alerts from the IceCube Neutrino Observatory, with a median latency of 12.2 hours from initial neutrino detection. From two of these campaigns, we have already reported tidal disruption event (TDE) AT2019dsg and likely TDE AT2019fdr as probable counterparts, suggesting that TDEs contribute >7.8% of the astrophysical neutrino flux. We here present the full results of our program through to December 2021. No additional candidate neutrino sources were identified by our program, allowing us to place the first constraints on the underlying optical luminosity function of astrophysical neutrino sources. Transients with optical absolutes magnitudes brighter that -21 can contribute no more than 87% of the total, while transients brighter than -22 can contribute no more than 58% of the total, neglecting the effect of extinction. These are the the first observational constraints on the neutrino emission of bright populations such as superluminous supernovae. None of the neutrinos were coincident with bright optical AGN flares comparable to that observed for TXS 0506+056/IC170922A, suggesting that most astrophysical neutrinos are not produced during such optical flares. We highlight the outlook for electromagnetic neutrino follow-up programs, including the expected potential for the Rubin Observatory., Submitted to MNRAS, comments welcome!

We present a systematic study of periodic X-ray sources in the massive globular cluster 47 Tuc, utilizing deep archival Chandra observations that resolve the cluster core and recently available eROSITA observations that cover the cluster outskirt. By applying the Gregory-Loredo algorithm, we detect 20 periodic signals among 18 X-ray sources, ranging between 205-95731 second. Fourteen periods are newly discovered in the X-ray band. We classify these periodic sources into four quiescent low-mass X-ray binaries, one milli-second pulsar, two coronally-active binaries and eleven cataclysmic variables (CVs), based on their X-ray temporal and spectral properties, as well as multi-band information. Despite a small sample subject to potential selection bias against faint and non-magnetic CVs, the 11 CVs together define an orbital period distribution significantly different from that of the CVs previously found in the solar neighborhood and the Galactic bulge. In particular, there exists in 47 Tuc an apparent paucity of short-period CVs below the period gap, which might be attributed to a high occupation fraction of non-magnetic CVs. Also characteristic of the 47 Tuc CVs are an overabundance of long-period CVs with a subgiant donor, a substantial fraction of CVs within the period gap, and a steep radial surface density profile. These are best understood as a group of CVs having recently formed via dynamical interactions in the dense cluster core. Despite sufficient sensitivity of the X-ray data, only one periodic source is found between one-third of the half-light radius and the tidal radius, the nature of which is unclear., Comment: 20 pages, 10 figures, Accepted for publication in MNRAS

The determination of distance is fundamental in astrophysics. Gamma-ray sources are poorly characterized in this sense, as the limited angular resolution and poor photon-count statistics in gamma-ray astronomy makes it difficult to associate them to a multiwavelength object with known redshift. Taking the 1794 active galactic nuclei (AGNs) with known redshift from the Fermi-LAT latest AGN catalog, 4LAC-DR3, we employ machine learning techniques to predict the distance of the rest of AGNs based on their spectral and spatial properties. The state-of-the-art CatBoost algorithm reaches an average 0.56 R2 score with 0.46 root-mean-squared error (RMSE), predicting an average redshift value of $z_{avg}=0.63$, with a maximum $z_{max}=1.97$. We use the SHAP explainer package to gain insights into the variables influence on the outcome, and also study the extragalactic bakground light (EBL) implications. In a second part, we use this regression model to predict the redshift of the unassociated sample of the latest LAT point-source catalog, 4FGL-DR3, using the results of a previous paper to determine the possible AGNs within them., 7 pages, 7 figures. Matches the accepted MNRAS version