Pulsars, transients and relativistic astrophysics

Abstract:

Outflows are critical components of many astrophysical systems, including accreting compact binaries and active galactic nuclei (AGN). These outflows can significantly affect a system’s evolution and alter its observational appearance by reprocessing the radiation produced by the central engine. sirocco (Simulating Ionization and Radiation in Outflows Created by Compact Objects – or ‘the code formerly known as python’) is a Sobolev-based Monte Carlo ionization and radiative transfer code. It is designed to simulate the spectra produced by any system with an azimuthally symmetric outflow, from spherical stellar winds to rotating, biconical accretion disc winds. Wind models can either be parametrized or imported, e.g. from hydrodynamical simulations. The radiation sources include an optically thick accretion disc and various central sources with flexible spectra and geometries. The code tracks the ‘photon packets’ produced by the sources in any given simulation as they traverse and interact with the wind. The code assumes radiative near-equilibrium, so the thermal and ionization state can be determined iteratively from these interactions. Once the physical properties in the wind have converged, sirocco can be used to generate synthetic spectra at a series of observer sightlines. Here, we describe the physical assumptions, operation, performance and limitations of the code. We validate it against tardis, cmfgen, and cloudy, finding good agreement, and present illustrative synthetic spectra from disc winds in cataclysmic variables, tidal disruption events, AGN, and X-ray binaries. sirocco is publicly available on GitHub, alongside its associated data, documentation and sample input files covering a wide range of astrophysical applications.

Abstract:

While dynamical tides only become relevant during the last couple of orbits for circular inspirals, orbital eccentricity can increase their impact during earlier phases of the inspiral by exciting tidal oscillations at each close encounter. We investigate the effect of dynamical tides on the orbital evolution of eccentric neutron star binaries using post-Newtonian numerical simulations and construct an analytic stochastic model that reproduces the numerical results. Our study reveals a strong dependence of dynamical tides on the pericenter distance, with the fractional energy transferred to dynamical tides over that dissipated in gravitational waves (GWs) exceeding ∼1% at separations rp≲50 km for large eccentricities. We demonstrate that the effect of dynamical tides on orbital evolution can manifest as a phase shift in the GW signal. We show that the signal-to-noise ratio of the GW phase shift can reach the detectability threshold of 8 with a single advanced Laser Interferometer Gravitational Wave Observatory detector at design sensitivity for eccentric neutron star binaries at a distance of 40 Mpc. This requires a pericenter distance of rp0≲68 km (rp0≲76 km) at binary formation with eccentricity close to 1 for a reasonable tidal deformability and f-mode frequency of 500 and 1.73 kHz (700 and 1.61 kHz), respectively. The observation of the phase shift will enable measuring the f-mode frequency of neutron stars independently from their tidal deformability, providing significant insights into neutron star seismology and the properties of the equation of state. We also explore the potential of distinguishing between equal-radius and twin-star binaries, which could provide an opportunity to reveal strong first-order phase transitions in the nuclear equation of state.

Abstract:

We present extensive observations of the Type II supernova (SN II) SN 2023ufx, which is likely the most metal-poor SN II observed to date. It exploded in the outskirts of a low-metallicity (Z host ∼ 0.1 Z ⊙) dwarf (M g = −13.39 ± 0.16 mag, r proj ∼ 1 kpc) galaxy. The explosion is luminous, peaking at M g ≈ −18.5 mag, and shows rapid evolution. The r-band (pseudobolometric) light curve has a shock-cooling phase lasting 20 (17) days followed by a 19 (23) day plateau. The entire optically thick phase lasts only ≈55 days following explosion, indicating that the red supergiant progenitor had a thinned H envelope prior to explosion. The early spectra obtained during the shock-cooling phase show no evidence for narrow emission features and limit the preexplosion mass-loss rate to Ṁ≲10−3 M ⊙ yr−1. The photospheric-phase spectra are devoid of prominent metal absorption features, indicating a progenitor metallicity of ≲0.1 Z ⊙. The seminebular (∼60–130 days) spectra reveal weak Fe ii, but other metal species typically observed at these phases (Ti ii, Sc ii, and Ba ii) are conspicuously absent. The late-phase optical and near-infrared spectra also reveal broad (≈104 km s−1) double-peaked Hα, Pβ, and Pγ emission profiles suggestive of a fast outflow launched during the explosion. Outflows are typically attributed to rapidly rotating progenitors, which also prefer metal-poor environments. This is only the second SN II with ≲0.1 Z ⊙ and both exhibit peculiar evolution, suggesting a sizable fraction of metal-poor SNe II have distinct properties compared to nearby metal-enriched SNe II. These observations lay the groundwork for modeling the metal-poor SNe II expected in the early Universe.