Researcher profile

Troels Haugboelle

Troels Haugboelle contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
12works
0followers
11topics
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

12 published item(s)

preprint2016arXiv

Supernova Driving. I. The Origin of Molecular Cloud Turbulence

Turbulence is ubiquitous in molecular clouds (MCs), but its origin is still unclear because MCs are usually assumed to live longer than the turbulence dissipation time. Interstellar medium (ISM) turbulence is likely driven by SN explosions, but it has never been demonstrated that SN explosions can establish and maintain a turbulent cascade inside MCs consistent with the observations. In this work, we carry out a simulation of SN-driven turbulence in a volume of (250 pc)$^3$, specifically designed to test if SN driving alone can be responsible for the observed turbulence inside MCs. We find that SN driving establishes a velocity scaling consistent with the usual scaling laws of supersonic turbulence, suggesting that previous idealized simulations of MC turbulence, driven with a random, large-scale volume force, were correctly adopted as appropriate models for MC turbulence, despite the artificial driving. We also find that the same scaling laws extend to the interior of MCs, and that the velocity-size relation of the MCs selected from our simulation is consistent with that of MCs from the Outer-Galaxy Survey, the largest MC sample available. The mass-size relation and the mass and size probability distributions also compare successfully with those of the Outer Galaxy Survey. Finally, we show that MC turbulence is super-Alfvénic with respect to both the mean and rms magnetic-field strength. We conclude that MC structure and dynamics are the natural result of SN-driven turbulence.

preprint2016arXiv

Tracking the distribution of $^{26}$Al and $^{60}$Fe during the early phases of star and disk evolution

The short-lived $^{26}$Al and $^{60}$Fe radionuclides are synthesized and expelled in the interstellar medium by core-collapse supernova events. The solar system's first solids, calcium-aluminium refractory inclusions (CAIs), contain evidence for the former presence of the $^{26}$ Al nuclide defining the canonical $^{26}$Al/$^{27}$ Al ratio of $\sim5 \times10^{-5}$. A different class of objects temporally related to canonical CAIs are CAIs with fractionation and unidentified nuclear effects (FUN CAIs), which record a low initial $^{26}$Al/$^{27}$Al of $10^{-6}$. The contrasting level of $^{26}$Al between these objects is often interpreted as reflecting the admixing of the $^{26}$Al nuclide during the early formative phase of the Sun. We use giant molecular cloud (GMC) scale adaptive mesh-refinement numerical simulations to trace the abundance of $^{26}$Al and $^{60}$Fe in star-forming gas during the early stages of accretion of individual low mass protostars. We find that the $^{26}$Al/$^{27}$Al and $^{60}$Fe/$^{56}$Fe ratios of accreting gas within a vicinity of 1000 AU of the stars follow the predicted decay curves of the initial abundances at time of star formation without evidence of spatial or temporal heterogeneities for the first 100 kyr of star formation. Therefore, the observed differences in $^{26}$Al/$^{27}$Al ratios between FUN and canonical CAIs are likely not caused by admixing of supernova material during the early evolution of the proto-Sun. Selective thermal processing of dust grains is a more viable scenario to account for the heterogeneity in $^{26}$Al/$^{27}$Al ratios at the time of solar system formation.

preprint2012arXiv

A Simple Law of Star Formation

We show that supersonic MHD turbulence yields a star formation rate (SFR) as low as observed in molecular clouds (MCs), for characteristic values of the free-fall time divided by the dynamical time, $t_{\rm ff}/t_{\rm dyn}$, the alfvénic Mach number, ${\cal M}_{\rm a}$, and the sonic Mach number, ${\cal M}_{\rm s}$. Using a very large set of deep adaptive-mesh-refinement simulations, we quantify the dependence of the SFR per free-fall time, $ε_{\rm ff}$, on the above parameters. Our main results are: i) $ε_{\rm ff}$ decreases exponentially with increasing $t_{\rm ff}/t_{\rm dyn}$, but is insensitive to changes in ${\cal M}_{\rm s}$, for constant values of $t_{\rm ff}/t_{\rm dyn}$ and ${\cal M}_{\rm a}$. ii) Decreasing values of ${\cal M}_{\rm a}$ (stronger magnetic fields) reduce $ε_{\rm ff}$, but only to a point, beyond which $ε_{\rm ff}$ increases with a further decrease of ${\cal M}_{\rm a}$. iii) For values of ${\cal M}_{\rm a}$ characteristic of star-forming regions, $ε_{\rm ff}$ varies with ${\cal M}_{\rm a}$ by less than a factor of two. We propose a simple star-formation law, based on the empirical fit to the minimum $ε_{\rm ff}$, and depending only on $t_{\rm ff}/t_{\rm dyn}$: $ε_{\rm ff} \approx ε_{\rm wind} \exp(-1.6 \,t_{\rm ff}/t_{\rm dyn})$. Because it only depends on the mean gas density and rms velocity, this law is straightforward to implement in simulations and analytical models of galaxy formation and evolution.

preprint2012arXiv

Halo abundances and shear in void models

We study the non-linear gravitational collapse of dark matter into halos through numerical N-body simulations of Lemaitre-Tolman-Bondi void models. We extend the halo mass function formalism to these models in a consistent way. This extension not only compares well with the simulated data at all times and radii, but it also gives interesting clues about the impact of the background shear on the growth of perturbations. Our results give hints about the possibility of constraining the background shear via cluster number counts, which could then give rise to strong constraints on general inhomogeneous models, of any scale.

preprint2012arXiv

Photon-Plasma: a modern high-order particle-in-cell code

We present the Photon-Plasma code, a modern high order charge conserving particle-in-cell code for simulating relativistic plasmas. The code is using a high order implicit field solver and a novel high order charge conserving interpolation scheme for particle-to-cell interpolation and charge deposition. It includes powerful diagnostics tools with on-the-fly particle tracking, synthetic spectra integration, 2D volume slicing, and a new method to correctly account for radiative cooling in the simulations. A robust technique for imposing (time-dependent) particle and field fluxes on the boundaries is also presented. Using a hybrid OpenMP and MPI approach the code scales efficiently from 8 to more than 250.000 cores with almost linear weak scaling on a range of architectures. The code is tested with the classical benchmarks particle heating, cold beam instability, and two-stream instability. We also present particle-in-cell simulations of the Kelvin-Helmholtz instability, and new results on radiative collisionless shocks.

preprint2011arXiv

Large scale structure simulations of inhomogeneous LTB void models

We perform numerical simulations of large scale structure evolution in an inhomogeneous Lemaitre-Tolman-Bondi (LTB) model of the Universe. We follow the gravitational collapse of a large underdense region (a void) in an otherwise flat matter-dominated Einstein-deSitter model. We observe how the (background) density contrast at the centre of the void grows to be of order one, and show that the density and velocity profiles follow the exact non-linear LTB solution to the full Einstein equations for all but the most extreme voids. This result seems to contradict previous claims that fully relativistic codes are needed to properly handle the non-linear evolution of large scale structures, and that local Newtonian dynamics with an explicit expansion term is not adequate. We also find that the (local) matter density contrast grows with the scale factor in a way analogous to that of an open universe with a value of the matter density OmegaM(r) corresponding to the appropriate location within the void.

preprint2011arXiv

Radiation from sub-Larmor scale magnetic fields

Spontaneous rapid growth of strong magnetic fields is rather ubiquitous in high-energy density environments ranging from astrophysical sources (e.g., gamma-ray bursts and relativistic shocks), to reconnection, to laser-plasma interaction laboratory experiments, where they are produced by kinetic streaming instabilities of the Weibel type. Relativistic electrons propagating through these sub-Larmor-scale magnetic fields radiate in the jitter regime, in which the anisotropy of the magnetic fields and the particle distribution have a strong effect on the produced radiation. Here we develop the general theory of jitter radiation, which includes (i) anisotropic magnetic fields and electron velocity distributions, (ii) the effects of trapped electrons and (iii) extends the description to large deflection angles of radiating particles thus establishing a cross-over between the classical jitter and synchrotron regimes. Our results are in remarkable agreement with the radiation spectra obtained from particle-in-cell simulations of the classical Weibel instability. Particularly interesting is the onset of the field growth, when the transient hard synchrotron-violating spectra are common as a result of the dominant role of the trapped population. This effect can serve as a distinct observational signature of the violent field growth in astrophysical sources and lab experiments. It is also interesting that a system with small-scale fields tends to evolve toward the small-angle jitter regime, which can, under certain conditions, dominate the overall emission of a source.

preprint2010arXiv

3D modeling of collisionless shocks

Two dimensional modeling of collisionless shocks has been of tremendous importance in understanding the physics of the non-linear evolution, momentum transfer and particle acceleration, but current computer capacities have now reached a point where three dimensional modeling is becoming feasible. We present the first three dimensional model of a fully developed and relaxed relativistic ion-electron shock, and analyze and compare it to similar 2D models. Quantitative and qualitative differences are found with respect to the two-dimensional models. The shock jump conditions are of course different, because of the extra degree of freedom, but in addition it is found that strong parallel electric fields develop at the shock interface, the level of magnetic field energy is lower, and the non-thermal particle distribution is shallower with a powerlaw index of ~2.2.

preprint2010arXiv

Neutrinos in Non-linear Structure Formation - The Effect on Halo Properties

We use N-body simulations to find the effect of neutrino masses on halo properties, and investigate how the density profiles of both the neutrino and the dark matter components change as a function of the neutrino mass. We compare our neutrino density profiles with results from the N-one-body method and find good agreement. We also show and explain why the Tremaine-Gunn bound for the neutrinos is not saturated. Finally we study how the halo mass function changes as a function of the neutrino mass and compare our results with the Sheth-Tormen semi-analytic formulae. Our results are important for surveys which aim at probing cosmological parameters using clusters, as well as future experiments aiming at measuring the cosmic neutrino background directly.

preprint2010arXiv

Non-Gaussianity from Axion Monodromy Inflation

We study the primordial non-Gaussinity predicted from simple models of inflation with a linear potential and superimposed oscillations. This generic form of the potential is predicted by the axion monodromy inflation model, that has recently been proposed as a possible realization of chaotic inflation in string theory, where the monodromy from wrapped branes extends the range of the closed string axions to beyond the Planck scale. The superimposed oscillations in the potential can lead to new signatures in the CMB spectrum and bispectrum. In particular the bispectrum will have a new distinct shape. We calculate the power spectrum and bispectrum of curvature perturbations in the model, as well as make analytic estimates in various limiting cases. From the numerical analysis we find that for a wide range of allowed parameters the model produces a feature in the bispectrum with fnl ~ 50 or larger while the power spectrum is almost featureless. This model is therefore an example of a string-inspired inflationary model which is testable mainly through its non-Gaussian features. Finally we provide a simple analytic fitting formula for the bispectrum which is accurate to approximately 5% in all cases, and easily implementable in codes designed to provide non-Gaussian templates for CMB analyses.

preprint2010arXiv

Radiation Spectral Synthesis of Relativistic Filamentation

Radiation from many astrophysical sources, e.g. gamma-ray bursts and active galactic nuclei, is believed to arise from relativistically shocked collisionless plasmas. Such sources often exhibit highly transient spectra evolving rapidly, compared with source lifetimes. Radiation emitted from these sources is typically associated with non-linear plasma physics, complex field topologies and non-thermal particle distributions. In such circumstances a standard synchrotron paradigm may fail to produce accurate conclusions regarding the underlying physics. Simulating spectral emission and spectral evolution numerically in various relativistic shock scenarios is then the only viable method to determine the detailed physical origin of the emitted spectra. In this Letter we present synthetic radiation spectra representing the early stage development of the filamentation (streaming) instability of an initially unmagnetized plasma, which is relevant for both collisionless shock formation and reconnection dynamics in relativistic astrophysical outflows, as well as for laboratory astrophysics experiments. Results were obtained using a highly efficient "in situ" diagnostics method, based on detailed particle-in-cell modeling of collisionless plasmas. The synthetic spectra obtained here are compared with those predicted by a semi-analytical model for jitter radiation from the filamentation instability, the latter including self-consistent generated field topologies and particle distributions obtained from the simulations reported upon here. Spectra exhibit dependence on the presence - or absence - of an inert plasma constituent, when comparing baryonic plasmas (i.e. containing protons) with pair plasmas. The results also illustrate that considerable care should be taken when using lower-dimensional models to obtain information about the astrophysical phenomena generating observed spectra.

preprint2008arXiv

Clustering properties of ultrahigh energy cosmic rays and the search for their astrophysical sources

The arrival directions of ultrahigh energy cosmic rays (UHECRs) may show anisotropies on all scales, from just above the experimental angular resolution up to medium scales and dipole anisotropies. We find that a global comparison of the two-point auto-correlation function of the data with the one of catalogues of potential sources is a powerful diagnostic tool. In particular, this method is far less sensitive to unknown deflections in magnetic fields than cross-correlation studies while keeping a strong discrimination power among source candidates. We illustrate these advantages by considering ordinary galaxies, gamma ray bursts and active galactic nuclei as possible sources. Already the sparse publicly available data suggest that the sources of UHECRs may be a strongly clustered sub-sample of galaxies or of active galactic nuclei. We present forecasts for various cases of source distributions which can be checked soon by the Pierre Auger Observatory.