Radiative transfer is a key component in almost all astrophysical and cosmological simulations. We present Magritte: a modern open-source software library for 3D radiative transfer. It uses a deterministic ray-tracer and formal solver, i.e. it computes the radiation field by tracing rays through the model and solving the radiative transfer equation in its second-order form along a fixed set of rays originating from each point. Magritte can handle structured and unstructured input meshes, as well as smoothed-particle hydrodynamics (SPH) particle data. In this first paper, we describe the numerical implementation, semi-analytic tests and cross-code benchmarks for the non-LTE line radiative transfer module of Magritte. This module uses the radiative transfer solver to self-consistently determine the populations of the quantised energy levels of atoms and molecules using an accelerated Lambda iteration (ALI) scheme. We compare Magritte with the established radiative transfer solvers Ratran (1D) and Lime (3D) on the van Zadelhoff benchmark and present a first application to a simple Keplerian disc model. Comparing with Lime, we conclude that Magritte produces more accurate and more precise results, especially at high optical depth, and that it is faster.

We aim to determine the distributions of molecular SiS and CS in the circumstellar envelopes of oxygen-rich asymptotic giant branch stars and how these distributions differ between stars that lose mass at different rates. In this study we analyse ALMA observations of SiS and CS emission lines for three oxygen-rich galactic AGB stars: IK Tau, with a moderately high mass-loss rate of $5\times10^{-6}$M$_\odot$ yr$^{-1}$, and W Hya and R Dor with low mass loss rates of $\sim1\times10^{-7}$M$_\odot$ yr$^{-1}$. These molecules are usually more abundant in carbon stars but the high sensitivity of ALMA allows us to detect their faint emission in the low mass-loss rate AGB stars. The high spatial resolution of ALMA also allows us to precisely determine the spatial distribution of these molecules in the circumstellar envelopes. We run radiative transfer models to calculate the molecular abundances and abundance distributions for each star. We find a spread of peak SiS abundances with $\sim10^{-8}$ for R Dor, $\sim10^{-7}$ for W Hya, and $\sim3\times10^{-6}$ for IK Tau relative to H$_2$. We find lower peak CS abundances of $\sim7\times10^{-9}$ for R Dor, $\sim7\times10^{-8}$ for W Hya and $\sim4\times10^{-7}$ for IK Tau, with some stratifications in the abundance distributions. For IK Tau we also calculate abundances for the detected isotopologues: C$^{34}$S, $^{29}$SiS, $^{30}$SiS, Si$^{33}$S, Si$^{34}$S, $^{29}$Si$^{34}$S, and $^{30}$Si$^{34}$S. Overall the isotopic ratios we derive for IK Tau suggest a lower metallicity than solar.

(abridged) In the inner regions of AGB outflows, several molecules have been detected with abundances much higher than those predicted from thermodynamic equilibrium (TE) chemical models. The presence of the majority of these species can be explained by shock-induced non-TE chemical models, where shocks caused by the pulsating star take the chemistry out of TE in the inner region. Moreover, a non-uniform density structure has been detected in several AGB outflows. A detailed parameter study on the quantitative effects of a non-homogeneous outflow has so far not been performed. We implement a porosity formalism for treating the increased leakage of light associated with radiation transport through a clumpy, porous medium. The effects from the altered UV radiation field penetration on the chemistry, accounting also for the increased reaction rates of two-body processes in the overdense clumps, are examined. We present a parameter study of the effect of clumping and porosity on the chemistry throughout the outflow. Both the higher density within the clumps and the increased UV radiation field penetration have an important impact on the chemistry, as they both alter the chemical pathways. The increased amount of UV radiation in the inner region leads to photodissociation of parent species, releasing the otherwise deficient elements. We find an increased abundance in the inner region of all species not expected to be present assuming TE chemistry, such as HCN in O-rich outflows, H$_2$O in C-rich outflows, and NH$_3$ in both. Outflows whose clumps have a large overdensity and that are very porous to the interstellar UV radiation field yield abundances comparable to those observed in O- and C-rich outflows for most of the unexpected species investigated. The inner wind abundances of H$_2$O in C-rich outflows and of NH$_3$ in O- and C-rich outflows are however underpredicted.