Solving radiative transfer with line overlaps using Gauss Seidel algorithms

Abstract

The improvement in observational facilities requires refining the modelling of the geometrical structures of astrophysical objects. Nevertheless, for complex problems such as line overlap in molecules showing hyperfine structure, a detailed analysis still requires a large amount of computing time and thus, misinterpretation cannot be dismissed due to an undersampling of the whole space of parameters. We extend the discussion of the implementation of the Gauss--Seidel algorithm in spherical geometry and include the case of hyperfine line overlap. We first review the basics of the short characteristics method that is used to solve the radiative transfer equations. Details are given on the determination of the Lambda operator in spherical geometry. The Gauss--Seidel algorithm is then described and, by analogy to the plan--parallel case, we see how to introduce it in spherical geometry. Doing so requires some approximations in order to keep the algorithm competitive. Finally, line overlap effects are included. The convergence speed of the algorithm is compared to the usual Jacobi iterative schemes. The gain in the number of iterations is typically factors of 2 and 4 for the two implementations made of the Gauss--Seidel algorithm. This is obtained despite the introduction of approximations in the algorithm. A comparison of results obtained with and without line overlaps for N2H+, HCN, and HNC shows that the J=3-2 line intensities are significantly underestimated in models where line overlap is neglected.

11 pages, 13 figures, accepted in A&A

References in corpus (5)

• Numerical methods for non-LTE line radiative transfer: Performance and convergence characteristics
• Water Vapor around Sgr B2
• Fast multilevel radiative transfer
• Fast 2D non-LTE radiative modelling of prominences I. Numerical methods and benchmark results
• Very Efficient Methods for Multilevel Radiative Transfer in Atomic and Molecular Lines

Cited by in corpus (21)

• Molecular abundances in the inner layers of IRC +10216
• Nitrogen isotopic ratios in Barnard 1: a consistent study of the N2H+, NH3, CN, HCN and HNC isotopologues
• A survey of C2H, HCN, and C18O in protoplanetary disks
• Unveiling the dust nucleation zone of IRC+10216 with ALMA
• Observation of HCN hyperfine line anomalies towards low- and high-mass star-forming cores
• Hyperfine excitation of N${}_2$H${}^{+}$ by H${}_2$: Toward a revision of N${}_2$H${}^{+}$ abundance in cold molecular clouds
• Multiple nitrogen reservoirs in a protoplanetary disk at the epoch of comet and giant planet formation
• Abundance of HCN and its C and N isotopologues in L1498
• The complex dust formation zone of the AGB star IRC+10216 probed with CARMA 0.25 arcsec angular resolution molecular observations
• Tentative Detection of the Nitrosylium Ion in Space
• N${}_2$H${}^{+}$ and N${}^{15}$NH${}^{+}$ towards the prestellar core 16293E in L1689N
• The 35Cl/37Cl isotopic ratio in dense molecular clouds: HIFI observations of hydrogen chloride towards W3A
• Anomalous HCN emission from warm giant molecular clouds
• Collisional excitation of doubly and triply deuterated ammonia ND${}_2$H and ND${}_3$ by H${}_2$
• Stratified NH and ND emission in the prestellar core 16293E in L1689N
• Herschel/HIFI observation of highly excited rotational lines of HNC toward IRC +10 216
• HCN emission from translucent gas and UV-illuminated cloud edges revealed by wide-field IRAM 30m maps of Orion B GMC: Revisiting its role as tracer of the dense gas reservoir for star formation
• Non--local radiative transfer in strongly inverted masers
• Preconditioned Bi-Conjugate Gradient Method for Radiative Transfer in Spherical Media
• On the influence of collisional rate coefficients on the water vapour excitation
• The maser emitting structure and time variability of the SiS lines J=14-13 and 15-14 in IRC+10216