Researcher profile

Troels Haugbølle

Troels Haugbølle contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
11works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

11 published item(s)

preprint2022arXiv

$ν\scriptstyle{\rm CO}N{\rm CEPT}$: Cosmological neutrino simulations from the non-linear Boltzmann hierarchy

In this paper the non-linear effect of massive neutrinos on cosmological structures is studied in a conceptually new way. We have solved the non-linear continuity and Euler equations for the neutrinos on a grid in real space in $N$-body simulations, and closed the Boltzmann hierarchy at the non-linear Euler equation using the stress and pressure perturbations from linear theory. By comparing with state-of-the art cosmological neutrino simulations, we are able to simulate the non-linear neutrino power spectrum very accurately. This translates into a negligible error in the matter power spectrum, and so our $\scriptstyle{\rm CO}N{\rm CEPT}$ code is ideally suited for extracting the neutrino mass from future high precision non-linear observational probes such as Euclid.

preprint2022arXiv

Natural separation of two primordial planetary reservoirs in an expanding solar protoplanetary disk

Meteorites display an isotopic composition dichotomy between non-carbonaceous (NC) and carbonaceous (CC) groups, indicating that planetesimal formation in the solar protoplanetary disk occurred in two distinct reservoirs. The prevailing view is that a rapidly formed Jupiter acted as a barrier between these reservoirs. We show a fundamental inconsistency in this model: if Jupiter is an efficient blocker of drifting pebbles, then the interior NC reservoir is depleted by radial drift within a few hundred thousand years. If Jupiter lets material pass it, then the two reservoirs will be mixed. Instead, we demonstrate that the arrival of the CC pebbles in the inner disk is delayed for several million years by the viscous expansion of the protoplanetary disk. Our results support that Jupiter formed in the outer disk (>10 AU) and allowed a considerable amount of CC material to pass it and become accreted by the terrestrial planets.

preprint2020arXiv

The Origin of Massive Stars: The Inertial-Inflow Model

We address the problem of the origin of massive stars, namely the origin, path and timescale of the mass flows that create them. Based on extensive numerical simulations, we propose a scenario where massive stars are assembled by large-scale, converging, inertial flows that naturally occur in supersonic turbulence. We refer to this scenario of massive-star formation as the "Inertial-Inflow Model". This model stems directly from the idea that the mass distribution of stars is primarily the result of turbulent fragmentation. Under this hypothesis, the statistical properties of the turbulence determine the formation timescale and mass of prestellar cores, posing definite constraints on the formation mechanism of massive stars. We quantify such constraints by the analysis of a simulation of supernova-driven turbulence in a 250-pc region of the interstellar medium, describing the formation of hundreds of massive stars over a time of approximately 30 Myr. Due to the large size of our statistical sample, we can say with full confidence that massive stars in general do not form from the collapse of massive cores, nor from competitive accretion, as both models are incompatible with the numerical results. We also compute synthetic continuum observables in Herschel and ALMA bands. We find that, depending on the distance of the observed regions, estimates of core mass based on commonly-used methods may exceed the actual core masses by up to two orders of magnitude, and that there is essentially no correlation between estimated and real core masses.

preprint2016arXiv

Supernova Driving. III. Synthetic Molecular Cloud Observations

We present a comparison of molecular clouds (MCs) from a simulation of supernova-driven interstellar medium (ISM) turbulence with real MCs from the Outer Galaxy Survey. The radiative transfer calculations to compute synthetic CO spectra are carried out assuming the CO relative abundance depends only on gas density, according to four different models. Synthetic MCs are selected above a threshold brightness temperature value, $T_{\rm B,min}=1.4$ K, of the $J=1-0$ $^{12}$CO line, generating 16 synthetic catalogs (four different spatial resolutions and four CO abundance models), each containing up to several thousands MCs. The comparison with the observations focuses on the mass and size distributions and on the velocity-size and mass-size Larson relations. The mass and size distributions are found to be consistent with the observations, with no significant variations with spatial resolution or chemical model, except in the case of the unrealistic model with constant CO abundance. The velocity-size relation is slightly too steep for some of the models, while the mass-size relation is a bit too shallow for all models only at a spatial resolution $dx\approx 1$ pc. The normalizations of the Larson relations show a clear dependence on spatial resolution, for both the synthetic and the real MCs. The comparison of the velocity-size normalization suggests that the SN rate in the Perseus arm is approximately 70\% or less of the rate adopted in the simulation. Overall, the realistic properties of the synthetic clouds confirm that supernova-driven turbulence can explain the origin and dynamics of MCs.

preprint2015arXiv

Large-scale numerical simulations of star formation put to the test: Comparing synthetic images and actual observations for statistical samples of protostars

(abridged) Context: Both observations and simulations of embedded protostars have progressed rapidly in recent years. Bringing them together is an important step in advancing our knowledge about the earliest phases of star formation. Aims: To compare synthetic continuum images and SEDs, calculated from large-scale numerical simulations, to observational studies, thereby aiding in both the interpretation of the observations and in testing the fidelity of the simulations. Methods: The radiative transfer code RADMC-3D is used to create synthetic continuum images and SEDs of protostellar systems in a large numerical simulation of a molecular cloud. More than 13000 unique radiative transfer models are produced of a variety of different protostellar systems. Results: Over the course of 0.76 Myr the simulation forms more than 500 protostars, primarily within two sub-clusters. Synthetic SEDs are used to the calculate evolutionary tracers Tbol and Lsmm/Lbol. It is shown that, while the observed distributions of the tracers are well matched by the simulation, they generally do a poor job of tracking the protostellar ages. Disks form early in the simulation, with 40 % of the Class 0 protostars being encircled by one. The flux emission from the simulated disks is found to be, on average, a factor of 6 too low relative to real observations. The distribution of protostellar luminosities spans more than three order of magnitudes, similar to the observed distribution. Cores and protostars are found to be closely associated with one another, with the distance distribution between them being in excellent agreement with observations. Conclusions: The analysis and statistical comparison of synthetic observations to real ones is established as a powerful tool in the interpretation of observational results.

preprint2015arXiv

On the local variation of the Hubble constant

We have carefully studied how local measurements of the Hubble constant, $H_0$, can be influenced by a variety of different parameters related to survey geometry, depth, and size, as well as observer position in space. Our study is based on N-body simulations of structure in the standard $Λ$CDM model and our conclusion is that the expected variance in measurements of $H_0$ is far too small to explain the current discrepancy between the low value of $H_0$ inferred from measurements of the cosmic microwave background (CMB) by the Planck collaboration and the value measured directly in the local universe by use of Type Ia supernovae. This conclusion is very robust and does not change with different assumptions about survey geometry and effective sky coverage or observer position in space.

preprint2015arXiv

Protostellar accretion traced with chemistry: Comparing synthetic C18O maps of embedded protostars to real observations

Context: Understanding how protostars accrete their mass is a central question of star formation. One aspect of this is trying to understand whether the time evolution of accretion rates in deeply embedded objects is best characterised by a smooth decline from early to late stages or by intermittent bursts of high accretion. Aims: We create synthetic observations of deeply embedded protostars in a large numerical simulation of a molecular cloud, which are compared directly to real observations. The goal is to compare episodic accretion events in the simulation to observations and to test the methodology used for analysing the observations. Methods: Simple freeze-out and sublimation chemistry is added to the simulation, and synthetic C$^{18}$O line cubes are created for a large number of simulated protostars. The spatial extent of C$^{18}$O is measured for the simulated protostars and compared directly to a sample of 16 deeply embedded protostars observed with the Submillimeter Array. If CO is distributed over a larger area than predicted based on the protostellar luminosity, it may indicate that the luminosity has been higher in the past and that CO is still in the process of refreezing. Results: Approximately 1% of the protostars in the simulation show extended C$^{18}$O emission, as opposed to approximately 50% in the observations, indicating that the magnitude and frequency of episodic accretion events in the simulation is too low relative to observations. The protostellar accretion rates in the simulation are primarily modulated by infall from the larger scales of the molecular cloud, and do not include any disk physics. The discrepancy between simulation and observations is taken as support for the necessity of disks, even in deeply embedded objects, to produce episodic accretion events of sufficient frequency and amplitude.

preprint2014arXiv

ALMA observations of the kinematics and chemistry of disc formation

Context: The R CrA cloud hosts a handful of Class 0/I low-mass young stellar objects. The chemistry and physics at scales $>500$ AU in this cloud are dominated by the irradiation from the nearby Herbig Be star R CrA. The luminous large-scale emission makes it necessary to use high-resolution spectral imaging to study the chemistry and dynamics of the inner envelopes and discs of the protostars. Aims: We aim to better understand the structure of the inner regions of these protostars and, in particular, the interplay between the chemistry and the presence of discs. Methods: Using Atacama Large Millimeter/submillimeter Array (ALMA) high-resolution spectral imaging interferometry observations, we study the molecular line and dust continuum emission at submillimetre wavelengths. Results: We detect dust continuum emission from four circumstellar discs around Class 0/I objects within the R CrA cloud. Towards IRS7B we detect C$^{17}$O emission showing a rotation curve consistent with a Keplerian disc with a well-defined edge that gives a good estimate for the disc radius at 50 AU. We derive the central object mass to $2.3M_{\odot}$ and the disc mass to $0.024M_{\odot}$. The observations are also consistent with a model of material infalling under conservation of angular momentum; however, this model provides a worse fit to the data. We also report a likely detection of faint CH$_3$OH emission towards this point source, as well as more luminous CH$_3$OH emission in an outflow orthogonal to the major axis of the C$^{17}$O emission. Conclusions: The faint CH$_3$OH emission seen towards IRS7B can be explained by a flat density profile of the inner envelope caused by the disc with a radius $\lesssim50$ AU. We propose that the regions of the envelopes where complex organic molecules are present in Class 0/I young stellar objects can become quenched as the disc grows.

preprint2014arXiv

Infall-Driven Protostellar Accretion and the Solution to the Luminosity Problem

We investigate the role of mass infall in the formation and evolution of protostars. To avoid ad hoc initial and boundary conditions, we consider the infall resulting self-consistently from modeling the formation of stellar clusters in turbulent molecular clouds. We show that infall rates in turbulent clouds are comparable to accretion rates inferred from protostellar luminosities or measured in pre-main-sequence stars. They should not be neglected in modeling the luminosity of protostars and the evolution of disks, even after the embedded protostellar phase. We find large variations of infall rates from protostar to protostar, and large fluctuations during the evolution of individuals protostars. In most cases, the infall rate is initially of order 10$^{-5}$\msun\ yr$^{-1}$, and may either decay rapidly in the formation of low-mass stars, or remain relatively large when more massive stars are formed. The simulation reproduces well the observed characteristic values and scatter of protostellar luminosities and matches the observed protostellar luminosity function. The luminosity problem is therefore solved once realistic protostellar infall histories are accounted for, with no need for extreme accretion episodes. These results are based on a simulation of randomly-driven magneto-hydrodynamic turbulence on a scale of 4pc, including self-gravity, adaptive-mesh refinement to a resolution of 50AU, and accreting sink particles. The simulation yields a low star formation rate, consistent with the observations, and a mass distribution of sink particles consistent with the observed stellar initial mass function during the whole duration of the simulation, forming nearly 1,300 sink particles over 3.2 Myr.

preprint2011arXiv

Neutrinos in Non-linear Structure Formation - a Simple SPH Approach

We present a novel method for implementing massive neutrinos in N-body simulations. Instead of sampling the neutrino velocity distribution by individual point particles we take neutrino free-streaming into account by treating it as an effective redshift dependent sound speed in a perfect isothermal fluid, and assume a relation between the sound speed and velocity dispersion of the neutrinos. Although the method fails to accurately model the true neutrino power spectrum, it is able to calculate the total matter power spectrum to the same accuracy as more complex hybrid neutrino methods, except on very small scales. We also present an easy way to update the publicly available Gadget-2 version with this neutrino approximation.

preprint2011arXiv

The effect of peculiar velocities on supernova cosmology

We present an analysis of peculiar velocities and their effect on supernova cosmology. In particular, we study (a) the corrections due to our own motion, (b) the effects of correlations in peculiar velocities induced by large-scale structure, and (c) uncertainties arising from a possible local under- or over-density. For all of these effects we present a case study of their impact on the cosmology derived by the Sloan Digital Sky Survey-II Supernova Survey (SDSS-II SN Survey). Correcting supernova redshifts for the CMB dipole slightly over-corrects nearby supernovae that share some of our local motion. We show that while neglecting the CMB dipole would cause a shift in the derived equation of state of Delta w ~ 0.04 (at fixed matter density) the additional local-motion correction is currently negligible (Delta w<0.01). We use a covariance-matrix approach to statistically account for correlated peculiar velocities. This down-weights nearby supernovae and effectively acts as a graduated version of the usual sharp low-redshift cut. Neglecting coherent velocities in the current sample causes a systematic shift of ~2% in the preferred value of w and will therefore have to be considered carefully when future surveys aim for percent-level accuracy. Finally, we perform n-body simulations to estimate the likely magnitude of any local density fluctuation (monopole) and estimate the impact as a function of the low-redshift cutoff. We see that for this aspect the low-z cutoff of z=0.02 is well-justified theoretically, but that living in a putative local density fluctuation leaves an indelible imprint on the magnitude-redshift relation.