SIROCCO: a publicly available Monte Carlo ionization and radiative transfer code for astrophysical outflows

Monthly Notices of the Royal Astronomical Society Oxford University Press 536:1 (2024) 879-904

Authors:

James H Matthews, Knox S Long, Christian Knigge, Stuart A Sim, Edward J Parkinson, Nick Higginbottom, Samuel W Mangham, Nicolas Scepi, Austen Wallis, Henrietta A Hewitt, Amin Mosallanezhad

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.

Simultaneous optical and X-ray detection of a Thermonuclear Burst in the 2024 outburst of EXO 0748–676

Monthly Notices of the Royal Astronomical Society: Letters Oxford University Press (OUP) (2024) slae103

Authors:

Amy H Knight, Lauren Rhodes, Douglas JK Buisson, James H Matthews, Noel Castro Segura, Adam Ingram, Matthew Middleton, Timothy P Roberts

Distance estimation of gamma-ray-emitting BL Lac objects from imaging observations

Astronomy & Astrophysics EDP Sciences 691 (2024) a154

Authors:

K Nilsson, V Fallah Ramazani, E Lindfors, P Goldoni, J Becerra González, JA Acosta Pulido, R Clavero, J Otero-Santos, T Pursimo, S Pita, PM Kouch, C Boisson, M Backes, G Cotter, F D’Ammando, E Kasai

Hidden by a star: The redshift and the offset broad line of the flat-spectrum radio quasar PKS 0903–57

Astronomy & Astrophysics EDP Sciences 691 (2024) l5

Authors:

P Goldoni, C Boisson, S Pita, F D’Ammando, E Kasai, W Max-Moerbeck, M Backes, G Cotter

The Early Radio Afterglow of Short GRB 230217A

The Astrophysical Journal Letters American Astronomical Society 975:1 (2024) L13

Authors:

GE Anderson, G Schroeder, AJ van der Horst, L Rhodes, A Rowlinson, A Bahramian, SI Chastain, BP Gompertz, PJ Hancock, T Laskar, JK Leung, RAMJ Wijers

Abstract:

We present the radio afterglow of short gamma-ray burst (GRB) 230217A, which was detected less than 1 day after the gamma-ray prompt emission with the Australia Telescope Compact Array (ATCA) and the Karl G. Jansky Very Large Array. The ATCA rapid-response system automatically triggered an observation of GRB 230217A following its detection by the Neil Gehrels Swift Observatory and began observing the event just 32 minutes postburst at 5.5 and 9 GHz for 7 hr. Dividing the 7 hr observation into three time-binned images allowed us to obtain radio detections with logarithmic central times of 1, 2.8, and 5.2 hr postburst, the first of which represents the earliest radio detection of any GRB to date. The decline of the light curve is consistent with reverse shock emission if the observing bands are below the spectral peak and not affected by synchrotron self-absorption. This makes GRB 230217A the fifth short GRB (SGRB) with radio detections attributed to a reverse shock at early times (<1 day postburst). Following brightness temperature arguments, we have used our early radio detections to place the highest minimum Lorentz factor ( Γmin>50 at ∼1 hr) constraints on a GRB in the radio band. Our results demonstrate the importance of rapid radio follow-up observations with long integrations and good sensitivity for detecting the fast-evolving radio emission from SGRBs and probing their reverse shocks.