LOC program for line radiative transfer

Abstract

Radiative transfer modelling is part of many astrophysical simulations and is used to make synthetic observations and to assist analysis of observations. We concentrate on the modelling of the radio lines emitted by the interstellar medium. In connection with high-resolution models, this can be significant computationally challenge. Our goal is a line radiative transfer (RT) program that makes good use of multi-core CPUs and GPUs. Parallelisation is essential to speed up computations and to enable the tackling of large modelling tasks with personal computers. The program LOC is based on ray-tracing and uses standard accelerated lambda iteration (ALI) methods for faster convergence. The program works on 1D and 3D grids. The 1D version makes use of symmetries to speed up the RT calculations. The 3D version works with octree grids and, to enable calculations with large models, is optimised for low memory usage. Tests show that LOC gives results that are in agreement with other RT codes to within ~2%. This is typical of code-to-code differences, which often are related to different interpretations of the model set-up. LOC run times compare favourably with those of Monte Carlo codes. In 1D tests, LOC runs were by up to a factor ~20 faster on a GPU than on a single CPU core. In spite of the complex path calculations, up to ~10 speed-up was observed also for 3D models using octree discretisation. GPUs enable calculations of models with hundreds of millions of cells, as encountered in the context of large-scale simulations of interstellar clouds. LOC shows good performance and accuracy and and is able to handle many RT modelling tasks on personal computers. Being written in Python, with the computing-intensive parts implemented as compiled OpenCL kernels, it can also a serve as a platform for further experimentation with alternative RT implementations.

accepted for publication in A&A

References in corpus (34)

• HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
• Dust Grain Size Distributions and Extinction in the Milky Way, LMC, and SMC
• A computer program for fast non-LTE analysis of interstellar line spectra
• An atomic and molecular database for analysis of submillimetre line observations
• Cores, filaments, and bundles: hierarchical core formation in the L1495/B213 Taurus region
• The SILCC (SImulating the LifeCycle of molecular Clouds) project: I. Chemical evolution of the supernova-driven ISM
• LIME - a flexible, non-LTE line excitation and radiation transfer method for millimeter and far-infrared wavelengths
• An ALMA study of the Orion Integral Filament: I. Evidence for narrow fibers in a massive cloud
• HELIOS: An Open-source, GPU-accelerated Radiative Transfer Code For Self-consistent Exoplanetary Atmospheres
• Formation and evolution of interstellar filaments; Hints from velocity dispersion measurements
• Sulphur chemistry in the L1544 pre-stellar core
• Numerical methods for non-LTE line radiative transfer: Performance and convergence characteristics
• Gas phase Elemental abundances in Molecular cloudS (GEMS) I. The prototypical dark cloud TMC 1
• Cosmological Radiative Transfer Comparison Project II: The Radiation-Hydrodynamic Tests
• The Cloud Factory I: Generating resolved filamentary molecular clouds from galactic-scale forces
• Hybrid Characteristics: 3D radiative transfer for parallel adaptive mesh refinement hydrodynamics
• Modeling molecular hyperfine line emission
• Why does ammonia not freeze out in the center of pre-stellar cores?
• Measuring the Alfvenic Nature of the Interstellar Medium: Velocity Anisotropy Revisited
• Self-gravitating Filament Formation from Shocked Flows: Velocity Gradients across Filaments
• Fully Threaded Transport Engine: New Method for Multi-Scale Radiative Transfer
• SOC program for dust continuum radiative transfer
• TurbuStat: Turbulence Statistics in Python
• Velocity Gradient in the Presence of Self-Gravity: Identifying Gravity-induced Inflow and Determining Collapsing Stage
• GPU-based Monte Carlo dust radiative transfer scheme applied to AGN
• Magritte, a modern software library for 3D radiative transfer: I. Non-LTE atomic and molecular line modelling
• Supernova Driving. III. Synthetic Molecular Cloud Observations
• The hierarchical fragmentation of filaments and the role of sub-filaments
• Radiative transfer on hierarchial grids
• Radiation Hydrodynamics using Characteristics on Adaptive Decomposed Domains for Massively Parallel Star Formation Simulations
• Accretion driven turbulence in filaments II: Effects of self-gravity
• ART${}^2$: A 3D Parallel Multi-wavelength Radiative Transfer Code for Continuum and Atomic and Molecular Lines
• Multi-Resolution Radiative Transfer for Line Emission
• ARC: Adaptive Ray-tracing with CUDA, a New Ray Tracing Code for Parallel GPUs

Cited by in corpus (18)

• The Central 1000 au of a Pre-stellar Core Revealed with ALMA. II. Almost Complete Freeze-out
• An Interferometric View of H-MM1. I. Direct Observation of NH3 Depletion
• Self-consistent dust and non-LTE line radiative transfer with SKIRT
• Multi-line observations of CH${}_3$OH, c-C${}_3$H${}_2$ and HNCO towards L1544: Dissecting the core structure with chemical differentiation
• Cyanopolyyne chemistry in the L1544 prestellar core: new insights from GBT observations
• Chemistry and dynamics of the prestellar core L1544
• Similar levels of deuteration in the pre-stellar core L1544 and the protostellar core HH211
• Hunting pre-stellar cores with APEX: IRAS16293E (Oph464)
• Tracing the contraction of the pre-stellar core L1544 with HC${}^{17}$O${}^{+}$ $J$ = 1-0 emission
• Mapping Synthetic Observations to Prestellar Core Models: An Interpretable Machine Learning Approach
• Hunting pre-stellar cores with APEX: Corona Australis 151, the densest pre-stellar core or the youngest protostar?
• Synthetic ngVLA line observations of a massive star-forming cloud
• A low cosmic-ray ionisation rate in the prestellar core Ophiuchus/H-MM1. Mapping of the molecular ions ortho-H2D+, N2H+, and DCO+
• Beetroots: spatially-regularized Bayesian inference of physical parameter maps -- Application to Orion
• Role of NH3 Binding Energy in the Early Evolution of Protostellar Cores
• Fast spectral line calculations with the escape probability method and tests with synthetic observations of interstellar clouds
• Properties of Contracting Massive Protostellar Cores
• The Dynamical State of Massive Clumps