Source author record

Tuomas Lunttila

Tuomas Lunttila appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.GA astro-ph.IM astro-ph

Catalog footprint

What is connected

5works

3topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2016arXiv

Composite biasing in Monte Carlo radiative transfer

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.

preprint2015arXiv

Benchmarking the Calculation of Stochastic Heating and Emissivity of Dust Grains in the Context of Radiative Transfer Simulations

We define an appropriate problem for benchmarking dust emissivity calculations in the context of radiative transfer (RT) simulations, specifically including the emission from stochastically heated dust grains. Our aim is to provide a self-contained guide for implementors of such functionality, and to offer insights in the effects of the various approximations and heuristics implemented by the participating codes to accelerate the calculations. 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; a series of analytically defined radiation fields to which the dust population is to be exposed; and instructions for the desired output. We process this problem using six RT codes participating in this benchmark effort, and compare the results to a reference solution computed with the publicly available dust emission code DustEM. 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 from 3 to 1000 micron. We conclude that the relevant modules in RT codes can and do produce fairly consistent results for the emissivity spectra of stochastically heated dust grains.

preprint2011arXiv

Mapping the prestellar core Ophiuchus D (L1696A) in ammonia

The gas kinetic temperature in the centres of starless, high-density cores is predicted to fall as low as 5-6 K. The aim of this study was to determine the kinetic temperature distribution in the low-mass prestellar core Oph D where previous observations suggest a very low central temperature. The densest part of the Oph D core was mapped in the NH3(1,1) and (2,2) inversion lines using the Very Large Array (VLA). The physical quantities were derived from the observed spectra by fitting the hyperfine structure of the lines, and subsequently the temperature distribution of Oph D was calculated using the standard rotational temperature techniques. A physical model of the cores was constructed, and the simulated spectra after radiative transfer calculations with a 3D Monte Carlo code were compared with the observed spectra. Temperature, density, and ammonia abundance of the core model were tuned until a satisfactory match with the observation was obtained. The high resolution of the interferometric data reveals that the southern part of Oph D comprises of two small cores. The gas kinetic temperatures, as derived from ammonia towards the centres of the southern and northern core are 7.4 and 8.9 K, respectively. The observed masses of the cores are only 0.2 M_Sun. Their potential collapse could lead to formation of brown dwarfs or low-mass stars.

preprint2009arXiv

The Super-Alfvénic Model of Molecular Clouds: Predictions for Mass-to-Flux and Turbulent-to-Magnetic Energy Ratios

Recent measurements of the Zeeman effect in dark-cloud cores provide important tests for theories of cloud dynamics and prestellar core formation. In this Letter we report results of simulated Zeeman measurements, based on radiative transfer calculations through a snapshot of a simulation of supersonic and super-Alfvénic turbulence. We have previously shown that the same simulation yields a relative mass-to-flux ratio (core versus envelope) in agreement with the observations (and in contradiction with the ambipolar-drift model of core formation). Here we show that the mass-to-flux and turbulent-to-magnetic-energy ratios in the simulated cores agree with observed values as well. The mean magnetic field strength in the simulation is very low, \bar{B}=0.34 \muG, presumably lower than the mean field in molecular clouds. Nonetheless, high magnetic field values are found in dense cores, in agreement with the observations (the rms field, amplified by the turbulence, is B_{rms}=3.05 \muG). We conclude that a strong large-scale mean magnetic field is not required by Zeeman effect measurements to date, although it is not ruled out by this work.

preprint2008arXiv

The Super-Alfvénic Model of Molecular Clouds: Predictions for Zeeman Splitting Measurements

We present synthetic OH Zeeman splitting measurements of a super-Alfvenic model of molecular clouds. We select dense cores from synthetic 13CO maps computed from the largest simulation to date of supersonic and super-Alfvenic turbulence. The synthetic Zeeman splitting measurements in the cores yield a relation between the magnetic field strength, B, and the column density, N, in good agreement with the observations. The large scatter in B at a fixed value of N is partly due to intrinsic variations in the magnetic field strength from core to core. We also compute the relative mass-to-flux ratio between the center of the cores and their envelopes, ${\cal R}_μ$, and show that super-Alfvenic turbulence produces a significant scatter also in ${\cal R}_μ$, including negative values (field reversal between core center and envelope). We find ${\cal R}_μ<1$ for 70% of the cores, and ${\cal R}_μ<0$ for 12%. Of the cores with $|B_{\rm LOS}|>10$ \muG, 81% have ${\cal R}_μ<1$. These predictions of the super-Alfvenic model are in stark contrast to the ambipolar drift model of core formation, where only ${\cal R}_μ>1$ is allowed.