SIROCCO: A Publicly Available Monte Carlo Ionization and Radiative Transfer Code for Astrophysical Outflows
arXiv:2410.19908 · doi:10.1093/mnras/stae2677
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.
26 pages, 16 figures. Published in MNRAS. This is the release paper for the SIROCCO code, which can be found at https://github.com/sirocco-rt/sirocco with links to documentation. Underlying data and figure scripts available at https://github.com/sirocco-rt/release-models -- comments on the code, paper or documentation are welcomed
References in corpus (112)
- Astropy: A Community Python Package for Astronomy
- The Astropy Project: Building an inclusive, open-science project and status of the v2.0 core package
- Observational Evidence of AGN Feedback
- PLUTO: a Numerical Code for Computational Astrophysics
- Black Holes, Galaxy Formation, and the M_BH-sigma Relation
- A Structure for Quasars
- Mass loss from hot massive stars
- Intrinsic disc emission and the Soft X-ray Excess in AGN
- Optical Flares from the Tidal Disruption of Stars by Massive Black Holes
- Radiative Transfer Simulations for Neutron Star Merger Ejecta
- A high velocity ionised outflow and XUV photosphere in the narrow emission line quasar PG1211+143
- Ubiquitous equatorial accretion disc winds in black hole soft states
- Unification of Luminous Type 1 Quasars through CIV Emission
- A Catalog of Broad Absorption Line Quasars from the Sloan Digital Sky Survey Third Data Release
- CHIANTI -- an atomic database for emission lines -- Paper XVI: Version 10, further extensions
- Seventeen Tidal Disruption Events from the First Half of ZTF Survey Observations: Entering a New Era of Population Studies
- AGN outflows and feedback twenty years on
- Time Dependent Monte Carlo Radiative Transfer Calculations For 3-Dimensional Supernova Spectra, Lightcurves, and Polarization
- A Massive X-Ray Outflow from the Quasar PDS 456
- SKIRT: an Advanced Dust Radiative Transfer Code with a User-Friendly Architecture
- The magnetic nature of disk accretion onto black holes
- The Frequency and Radio Properties of Broad Absorption Line Quasars
- The Suzaku view of highly-ionised outflows in AGN: I - Statistical detection and global absorber properties
- Continuum and Emission-Line Properties of Broad Absorption Line Quasars
- A Characteristic Scale for Cold Gas
- Mocassin: A fully three-dimensional Monte Carlo photoionization code
- A strong redshift dependence of the broad absorption line quasar fraction
- HYPERION: An open-source parallelized three-dimensional dust continuum radiative transfer code
- Radiative recombination data for modelling dynamic finite-density plasmas
- Time-dependent 3D spectrum synthesis for type Ia supernovae
- High resolution Chandra HETG and RXTE observations of GRS 1915+105 : A hot disk atmosphere & cold gas enriched in Iron and Silicon
- The Chandra HETGS and RXTE view of GRS 1915+105
- The dusty MOCASSIN: fully self-consistent 3D photoionisation and dust radiative transfer models
- The Intrinsic Fraction of Broad Absorption Line Quasars
- A spectral synthesis code for rapid modelling of supernovae
- The Many Routes to AGN Feedback
- Simultaneous Chandra and RXTE Spectroscopy of the Microquasar H~1743-322: Clues to Disk Wind and Jet Formation from a Variable Ionized Outflow
- Chandra/HETGS Spectroscopy of the Galactic Black Hole GX 339-4: A Relativistic Iron Line and Evidence for a Seyfert-like Warm Absorber
- A global view of the inner accretion and ejection flow around super massive black holes: radiation driven accretion disk winds in a physical context
- On the consistent treatment of the quasi-hydrostatic layers in hot star atmospheres
- On the role of the UV and X-ray radiation in driving a disk wind in X-ray binaries
- Geometric and Dynamical Models of Reverberation Mapping Data
- Multi-dimensional modelling of X-ray spectra for AGN accretion-disk outflows III: application to a hydrodynamical simulation
- Regulation of black-hole accretion by a disk wind during a violent outburst of V404 Cygni
- 2MASS Reveals a Large Intrinsic Fraction of BALQSOs
- Accretion disc atmospheres and winds in low-mass X-ray binaries
- grmonty: a Monte Carlo Code for Relativistic Radiative Transport
- Modeling reverberation mapping data I: improved geometric and dynamical models and comparison with cross-correlation results
- Line-driven Disk Winds in Active Galactic Nuclei: The Critical Importance of Ionization and Radiative Transfer
- Modeling the Spectral Signatures of Accretion Disk Winds: A New Monte Carlo Approach
- Broad absorption features in wind-dominated ULXs?
- Supercritical Accretion Discs in Ultraluminous X-ray Sources and SS 433
- Multidimensional models of hydrogen and helium emission line profiles for classical T Tauri Stars: method, tests and examples
- The Broad Absorption Line Tidal Disruption Event iPTF15af: Optical and Ultraviolet Evolution
- Monte Carlo transition probabilities
- Monte Carlo techniques for time-dependent radiative transfer in 3-D supernovae
- Monte Carlo transition probabilities. II
- Multi-dimensional modelling of X-ray spectra for AGN accretion-disk outflows
- BAL and non-BAL quasars: continuum, emission and absorption properties establish a common parent sample
- Monte Carlo Radiative Transfer
- Modelling the spectra of the kilonova AT2017gfo -- I: The photospheric epochs
- What Sets the Line Profiles in Tidal Disruption Events?
- Hard-state accretion disk winds from black holes: the revealing case of MAXI J1820+070
- The Impact of Accretion Disk Winds on the Optical Spectra of Cataclysmic Variables
- The origin of ultra-fast outflows in AGN: Monte-Carlo simulations of the wind in PDS 456
- The Ultraviolet Spectroscopic Evolution of the Low-Luminosity Tidal Disruption Event iPTF16fnl
- Testing Quasar Unification: Radiative Transfer in Clumpy Winds
- Polarisation spectral synthesis for Type Ia supernova explosion models
- Discovery of Highly Blueshifted Broad Balmer and Metastable Helium Absorption Lines in a Tidal Disruption Event
- A Simple Disk Wind Model for Broad Absorption Line Quasars
- Three-Dimensional Simulations of Dynamics of Accretion Flows Irradiated by a Quasar
- Absorption lines from magnetically-driven winds in X-ray binaries
- Super-Eddington Accretion in Tidal Disruption Events: the Impact of Realistic Fallback Rates on Accretion Rates
- Magnetically-driven jets and winds from weakly magnetized accretion disks
- Monte Carlo radiative transfer for the nebular phase of Type Ia supernovae
- Stratified disc wind models for the AGN broad-line region: ultraviolet, optical and X-ray properties
- The low-luminosity accretion disc wind of the black hole transient V4641 Sagittarii
- Testing AGN outflow and accretion models with CIV and HeII emission line demographics in z=2 quasars
- A persistent ultraviolet outflow from an accreting neutron star binary transient
- Thermal Disk Winds in X-ray Binaries: Realistic Heating and Cooling Rates Give Rise to Slow, but Massive Outflows
- Testing the line-driven disk wind model: time-resolved UV spectroscopy of IX Vel and V3885 Sgr
- A Balmer-line Broad Absorption Line Quasar
- 2-D Monte Carlo simulations of H I line formation in massive YSO disk winds
- The Geometry and Ionization Structure of the Wind in the Eclipsing Nova-like Variables RW Tri and UX UMa
- Photoionization calculations of the radiation force due to spectral lines in AGNs
- The thermal-radiative wind in low mass X-ray binary H 1743-322; Radiation hydrodynamic simulations
- Polarization in Monte Carlo radiative transfer and dust scattering polarization signatures of spiral galaxies
- Thermal and radiation driving can produce observable disk winds in hard-state X-ray binaries
- An equatorial outflow in the black hole optical dipper Swift J1357.2-0933
- Radiation-hydrodynamic simulations of thermally-driven disc winds in X-ray binaries: A direct comparison to GRO J1655-40
- Hot, dense He II outflows during the 2017 outburst of the X-ray transient Swift J1357.2-0933
- The luminosity dependence of thermally-driven disc winds in low-mass X-ray binaries
- What powers the wind from the black hole accretion disc in GRO J1655-40?
- Cloud coalescence: a dynamical instability affecting multiphase environments
- Synthetic X-ray spectra for simulations of the dynamics of an accretion flow irradiated by a quasar
- Do Reverberation Mapping Analyses Provide an Accurate Picture of the Broad Line Region?
- Magnetic wind-driven accretion in dwarf novae
- Modeling Magnetic Disk-Wind State Transitions in Black Hole X-ray Binaries
- Angular Momentum in Disk Wind Revealed in the Young Star MWC349A
- Self Regulated Shocks in Massive Star Binary Systems
- The reverberation signatures of rotating disc winds in active galactic nuclei
- Monte-Carlo methods for NLTE spectral synthesis of supernovae
- A disc wind model for blueshifts in quasar broad emission lines
- Disc wind models for FU Ori objects
- Quasar emission lines as probes of orientation: implications for disc wind geometries and unification
- Modelling the X-ray spectra of high velocity outflows from quasars
- Accretion Disc Winds in Tidal Disruption Events: Ultraviolet Spectral Lines as Orientation Indicators
- Thermally driven winds in ULXs
- The origin of optical emission lines in the soft state of X-ray binary outbursts: the case of MAXI J1820+070
- Radiative hydrodynamical simulations of super-Eddington accretion flow in tidal disruption event: the accretion flow and wind
- Unveiling optical signatures of outflows in accreting white dwarfs
- A disc wind origin for the optical spectra of dwarf novae in outburst
Cited by in corpus (5)
- A multi-dimensional view of a unified model for TDEs
- Accretion disc winds in X-ray binaries
- A Diagnostic Kit for Optical Emission Lines Shaped by Accretion Disc Winds
- Flux-Averaged Force Multipliers
- Unveiling BLR Structure in AGN with High Resolution X-Ray Spectra: An Analytic Approach to Wind Emission Line Profiles