| 1. |
- Baes, Maarten, et al.
(författare)
-
Composite biasing in Monte Carlo radiative transfer
- 2016
-
Ingår i: Astronomy and Astrophysics. - : EDP Sciences. - 0004-6361 .- 1432-0746. ; 590, s. Art. no. A55-
-
Tidskriftsartikel (refereegranskat)abstract
- Biasing or importance sampling is a powerful technique in Monte Carlo radiative transfer, and can be applied in different forms to increase the accuracy and efficiency of simulations. One of the drawbacks of the use of biasing is the potential introduction of large weight factors. We discuss a general strategy, composite biasing, to suppress the appearance of large weight factors. We use this composite biasing approach for two different problems faced by current state-of-the-art Monte Carlo radiative transfer codes: the generation of photon packages from multiple components, and the penetration of radiation through high optical depth barriers. In both cases, the implementation of the relevant algorithms is trivial and does not interfere with any other optimisation techniques. Through simple test models, we demonstrate the general applicability, accuracy and efficiency of the composite biasing approach. In particular, for the penetration of high optical depths, the gain in efficiency is spectacular for the specific problems that we consider: in simulations with composite path length stretching, high accuracy results are obtained even for simulations with modest numbers of photon packages, while simulations without biasing cannot reach convergence, even with a huge number of photon packages.
|
|
| 2. |
- Camps, Peter, et al.
(författare)
-
Benchmarking the calculation of stochastic heating and emissivity of dust grains in the context of radiative transfer simulations
- 2015
-
Ingår i: Astronomy and Astrophysics. - : EDP Sciences. - 0004-6361 .- 1432-0746. ; 580
-
Tidskriftsartikel (refereegranskat)abstract
- Context. Thermal emission by stochastically heated dust grains (SHGs) plays an important role in the radiative transfer (RT) problem for a dusty medium. It is therefore essential to verify that RT codes properly calculate the dust emission before studying the effects of spatial distribution and other model parameters on the simulated observables. Aims. We define an appropriate problem for benchmarking dust emissivity calculations in the context of RT simulations, specifically including the emission from SHGs. Our aim is to provide a self-contained guide for implementors of such functionality and to offer insight into the effects of the various approximations and heuristics implemented by the participating codes to accelerate the calculations. Methods. The benchmark problem definition includes the optical and calorimetric material properties and the grain size distributions for a typical astronomical dust mixture with silicate, graphite, and PAH components. It also includes a series of analytically defined radiation fields to which the dust population is to be exposed and instructions for the desired output. We processed this problem using six RT codes participating in this benchmark effort and compared the results to a reference solution computed with the publicly available dust emission code DustEM. Results. The participating codes implement different heuristics to keep the calculation time at an acceptable level.We study the effects of these mechanisms on the calculated solutions and report on the level of (dis)agreement between the participating codes. For all but the most extreme input fields, we find agreement within 10% across the important wavelength range 3 μm ≤1000 μm. Conclusions. We conclude that the relevant modules in RT codes can and do produce fairly consistent results for the emissivity spectra of SHGs. This work can serve as a reference for implementors of dust RT codes, and it will pave the way for a more extensive benchmark effort focusing on the RT aspects of the various codes.
|
|
| 3. |
- Feher, O., et al.
(författare)
-
A CO survey on a sample of Herschel cold clumps
- 2017
-
Ingår i: Astronomy and Astrophysics. - : EDP Sciences. - 0004-6361 .- 1432-0746. ; 606
-
Tidskriftsartikel (refereegranskat)abstract
- Context. The physical state of cold cloud clumps has a great impact on the process and efficiency of star formation and the masses of the forming stars inside these objects. The sub-millimetre survey of the Planck space observatory and the far-infrared follow-up mapping of the Herschel space telescope provide an unbiased, large sample of these cold objects. Aims. We have observed (CO)-C-12(1-0) and (CO)-C-13(1-0) emission in 35 high-density clumps in 26 Herschel fields sampling different environments in the Galaxy. Here, we aim to derive the physical properties of the objects and estimate their gravitational stability. Methods. The densities and temperatures of the clumps were calculated from both the dust continuum and the molecular line data. Kinematic distances were derived using (CO)-C-13(1-0) line velocities to verify previous distance estimates and the sizes and masses of the objects were calculated by fitting 2D Gaussian functions to their optical depth distribution maps on 250 mu m. The masses and virial masses were estimated assuming an upper and lower limit on the kinetic temperatures and considering uncertainties due to distance limitations. Results. The derived excitation temperatures are between 8.5-19.5 K, and for most clumps between 10 15 K, while the Herschel-derived dust colour temperatures are more uniform, between 12 16 K. The sizes (0.1-3 pc), (CO)-C-13 column densities (0.5-44 x 10(15) cm(-2)) and masses (from less than 0.1 M-circle dot to more than 1500 M-circle dot) of the objects all span broad ranges. We provide new kinematic distance estimates, identify gravitationally bound or unbound structures and discuss their nature. Conclusions. The sample contains objects on a wide scale of temperatures, densities and sizes. Eleven gravitationally unbound clumps were found, many of them smaller than 0.3 pc, but large, parsec-scale clouds with a few hundred solar masses appear as well. Colder clumps have generally high column densities but warmer objects appear at both low and higher column densities. The clump column densities derived from the line and dust observations correlate well, but are heavily affected by uncertainties of the dust properties, varying molecular abundances and optical depth effects.
|
|
| 4. |
- Gordon, K. D., et al.
(författare)
-
I. A 3D externally illuminated slab benchmark for dust radiative transfer
- 2017
-
Ingår i: Astronomy and Astrophysics. - : EDP Sciences. - 0004-6361 .- 1432-0746. ; 603, s. A114-
-
Tidskriftsartikel (refereegranskat)abstract
- Context. The radiative transport of photons through arbitrary three-dimensional (3D) structures of dust is a challenging problem due to the anisotropic scattering of dust grains and strong coupling between different spatial regions. The radiative transfer problem in 3D is solved using Monte Carlo or Ray Tracing techniques as no full analytic solution exists for the true 3D structures. Aims. We provide the first 3D dust radiative transfer benchmark composed of a slab of dust with uniform density externally illuminated by a star. This simple 3D benchmark is explicitly formulated to provide tests of the different components of the radiative transfer problem including dust absorption, scattering, and emission. Methods. The details of the external star, the slab itself, and the dust properties are provided. This benchmark includes models with a range of dust optical depths fully probing cases that are optically thin at all wavelengths to optically thick at most wavelengths. The dust properties adopted are characteristic of the diffuse Milky Way interstellar medium. This benchmark includes solutions for the full dust emission including single photon (stochastic) heating as well as two simplifying approximations: One where all grains are considered in equilibrium with the radiation field and one where the emission is from a single effective grain with size-distribution-averaged properties. A total of six Monte Carlo codes and one Ray Tracing code provide solutions to this benchmark. Results. The solution to this benchmark is given as global spectral energy distributions (SEDs) and images at select diagnostic wavelengths from the ultraviolet through the infrared. Comparison of the results revealed that the global SEDs are consistent on average to a few percent for all but the scattered stellar flux at very high optical depths. The image results are consistent within 10%, again except for the stellar scattered flux at very high optical depths. The lack of agreement between different codes of the scattered flux at high optical depths is quantified for the first time. Convergence tests using one of the Monte Carlo codes illustrate the sensitivity of the solutions to various model parameters. Conclusions. We provide the first 3D dust radiative transfer benchmark and validate the accuracy of this benchmark through comparisons between multiple independent codes and detailed convergence tests.
|
|
| 5. |
- Liseau, René, 1949, et al.
(författare)
-
Gas and dust in the star-forming region ρ Oph A ∗, ∗∗, ∗∗∗: The dust opacity exponent β and the gas-to-dust mass ratio g2d
- 2015
-
Ingår i: Astronomy and Astrophysics. - : EDP Sciences. - 0004-6361 .- 1432-0746. ; 578
-
Tidskriftsartikel (refereegranskat)abstract
- © ESO, 2015. Aims. We aim at determining the spatial distribution of the gas and dust in star-forming regions and address their relative abundances in quantitative terms. We also examine the dust opacity exponent β for spatial and/or temporal variations. Methods. Using mapping observations of the very dense ρ Oph A core, we examined standard 1D and non-standard 3D methods to analyse data of far-infrared and submillimetre (submm) continuum radiation. The resulting dust surface density distribution can be compared to that of the gas. The latter was derived from the analysis of accompanying molecular line emission, observed with Herschel from space and with APEX from the ground. As a gas tracer we used N2H+, which is believed to be much less sensitive to freeze-out than CO and its isotopologues. Radiative transfer modelling of the N2H+ (J = 3-2) and (J = 6-5) lines with their hyperfine structure explicitly taken into account provides solutions for the spatial distribution of the column density N(H2), hence the surface density distribution of the gas. Results. The gas-to-dust mass ratio is varying across the map, with very low values in the central regions around the core SM 1. The global average, = 88, is not far from the canonical value of 100, however. In ρ Oph A, the exponent β of the power-law description for the dust opacity exhibits a clear dependence on time, with high values of 2 for the envelope-dominated emission in starless Class -1 sources to low values close to 0 for the disk-dominated emission in Class III objects. β assumes intermediate values for evolutionary classes in between. Conclusions. Since β is primarily controlled by grain size, grain growth mostly occurs in circumstellar disks. The spatial segregation of gas and dust, seen in projection toward the core centre, probably implies that, like C18O, also N2H+ is frozen onto the grains.
|
|
