Source author record

J. Rosdahl

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

8works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

8 published item(s)

preprint2022arXiv

Sub-parsec resolution cosmological simulations of star-forming clumps at high redshift with feedback of individual stars

We introduce a new set of zoom-in cosmological simulations with sub-pc resolution, intended to model extremely faint, highly magnified star-forming stellar clumps, detected at z=6.14 thanks to gravitational lensing. The simulations include feedback from individual massive stars (in both the pre-supernova and supernova phases), generated via stochastic, direct sampling of the stellar initial mass function. We adopt a modified 'delayed cooling' feedback scheme, specifically created to prevent artificial radiative loss of the energy injected by individual stars in very dense gas (n~10^3-10^5 cm^{-3}). The sites where star formation ignites are characterised by maximum densities of the order of 10^5 cm^{-3} and gravitational pressures P/k>10^7 K/cm^3, corresponding to the values of the local, turbulent regions where the densest stellar aggregates form. The total stellar mass at z=6.14 is 3.4x10^7 M_sun, in satisfactory agreement with the observed stellar mass of the observed systems. The most massive clumps have masses of ~10^6 M_sun and half-mass sizes of ~100 pc. These sizes are larger than the observed ones, including also other samples of lensed high-redshift clumps, and imply an average density one order of magnitude lower than the observed one. In the size-mass plane, our clumps populate a sequence that is intermediate between the ones of observed high-redshift clumps and local dSph galaxies.

preprint2021arXiv

Second generation star formation in globular clusters of different masses

By means of three-dimensional hydrodynamical simulations, we investigate the formation of second generation (SG) stars in young globular clusters of different masses. We consider clusters with a first generation of asymptotic giant branch (AGB) stars with mass 10^5 and 10^6 Msun moving at constant velocity through a uniform gas with density 10^(-24) and 10^(-23) g cm^(-3). Our setup is designed to reproduce the encounter of a young cluster with a reservoir of dense gas, e. g. during its orbital motion in the host galaxy. In the low-density models, as a result of the cooling AGB ejecta which collect in the centre, weakly perturbed by the external ram pressure, a compact central He-rich SG stellar component is formed on a timescale which decreases with increasing initial cluster mass. Our high-density models are subject to stronger ram pressure, which prevents the accumulation of the most He-rich AGB ejecta in the cluster centre. As a result, the SG is more extended and less He-enhanced than in the low-density models. By combining our results with previous simulations, we are able to study relevant, cluster-related scaling relations across a dynamical range of two orders of magnitude in mass (from 10^5 Msun to 10^7 Msun). In agreement with current observationally-based estimates, we find positive correlations between the SG-to-total number ratio and maximum He enhancement in SG stars as a function of the initial cluster mass.

preprint2020arXiv

A new hybrid radiative transfer method for massive star formation

Frequency-dependent/hybrid approaches for stellar irradiation are of primary importance in numerical simulations of massive star formation. We seek to compare outflow and accretion mechanisms in star formation simulations. We investigate the accuracy of a hybrid radiative transfer method using the gray M1 closure relation for proto-stellar irradiation and gray flux-limited diffusion (FLD) for photons emitted everywhere else. We have coupled the FLD module of the adaptive-mesh refinement code Ramses with Ramses-RT, which is based on the M1 closure relation. Our hybrid (M1+FLD) method takes an average opacity at the stellar temperature for the M1 module, instead of the local environmental radiation field. We have tested this approach in radiative transfer tests of disks irradiated by a star for three levels of optical thickness and compared the temperature structure with RADMC-3D and MCFOST. We applied it to a radiation-hydrodynamical simulation of massive star formation. Our tests validate our hybrid approach for determining the temperature structure of an irradiated disk in the optically-thin and moderately optically-thick regimes and the most optically-thick test shows the limitation of our approach. The optically-thick setups highlight the ability of the hybrid method to partially capture the self-shielding in the disk while the FLD cannot. The radiative acceleration is 100 times greater with the hybrid method. It consistently leads to about +50% more extended and wider-angle radiative outflows in the massive star formation simulation. We obtain a $17.6 M_\odot$ at $t{\simeq}0.7 τ_\mathrm{ff}$, while the accretion phase is ongoing. Finally, despite the use of refinement to resolve the radiative cavities, no Rayleigh-Taylor instability appears in our simulations, and we justify their absence by physical arguments based on the entropy gradient. (abridged)

preprint2020arXiv

UV absorption lines and their potential for tracing the Lyman continuum escape fraction

The neutral intergalactic medium above redshift 6 is opaque to ionizing radiation, therefore one needs indirect measurements of the escape fraction of ionizing photons from galaxies of this epoch. Low-ionization state absorption lines are a common feature in the spectrum of galaxies, showing a diversity of strengths and shapes. Since these lines indicate the presence of neutral gas in front of the stars, they have been proposed to carry information on the escape of ionizing radiation from galaxies. We study which processes are responsible for the shape of the absorption lines, to better understand their origin. We then explore whether the absorption lines can be used to predict the escape fractions. Using a radiation-hydrodynamical zoom-in simulation and the radiative transfer code RASCAS, we generate mock CII 1334 and LyB lines of a virtual galaxy at redshift 3 as seen from many directions of observation. We also compute the escape fraction of ionizing photons in those directions and look for correlations between the lines and the escape fractions. We find that the resulting mock absorption lines are comparable to observations and that the lines and the escape fractions vary strongly depending on the direction of observation. Gas velocity and dust always affect the absorption profile significantly. We find no strong correlations between observable LyB or CII 1334 and the escape fraction. After correcting the continuum for attenuation by dust to recover the intrinsic continuum, the residual flux of CII 1334 correlates well with the escape fraction for directions with a dust corrected residual flux larger than 30%. For other directions, the relations have a strong dispersion, and the residual flux overestimates the escape fraction for most cases. Concerning LyB, the residual flux after dust correction does not correlate with the escape fraction but can be used as a lower limit. (abridged)

preprint2015arXiv

A scheme for radiation pressure and photon diffusion with the M1 closure in RAMSES-RT

We describe and test an updated version of radiation-hydrodynamics (RHD) in the RAMSES code, that includes three new features: i) radiation pressure on gas, ii) accurate treatment of radiation diffusion in an unresolved optically thick medium, and iii) relativistic corrections that account for Doppler effects and work done by the radiation to first order in v/c. We validate the implementation in a series of tests, which include a morphological assessment of the M1 closure for the Eddington tensor in an astronomically relevant setting, dust absorption in a optically semi-thick medium, direct pressure on gas from ionising radiation, convergence of our radiation diffusion scheme towards resolved optical depths, correct diffusion of a radiation flash and a constant luminosity radiation, and finally, an experiment from Davis et al. of the competition between gravity and radiation pressure in a dusty atmosphere, and the formation of radiative Rayleigh-Taylor instabilities. With the new features, RAMSES-RT can be used for state-of-the-art simulations of radiation feedback from first principles, on galactic and cosmological scales, including not only direct radiation pressure from ionising photons, but also indirect pressure via dust from multi-scattered IR photons reprocessed from higher-energy radiation, both in the optically thin and thick limits.

preprint2015arXiv

Modeling small galaxies during the Epoch of Reionisation

Small galaxies are thought to be the main contributors to the ionising budget of the Universe before reionisation was complete. There have been a number of numerical studies trying to quantify their ionising efficiency through the escape fraction $f_{esc}$. While there is a clear trend that $f_{esc}$ is higher for smaller haloes, there is a large scatter in the distribution of $f_{esc}$ for a single halo mass. We propose that this is due to the intrinsic burstiness of star formation in low mass galaxies. We performed high resolution radiative hydrodynamics simulations with Ramses-RT to model the evolution of three galaxies and their ionising efficiency. We found that the variability of $f_{esc}$ follows that of the star formation rate. We then discuss the consequences of this variability on the observability of such galaxies by JWST.

preprint2015arXiv

StarBench: The D-type expansion of an HII region

StarBench is a project focused on benchmarking and validating different star-formation and stellar feedback codes. In this first StarBench paper we perform a comparison study of the D-type expansion of an HII region. The aim of this work is to understand the differences observed between the twelve participating numerical codes against the various analytical expressions examining the D-type phase of HII region expansion. To do this, we propose two well-defined tests which are tackled by 1D and 3D grid- and SPH- based codes. The first test examines the `early phase' D-type scenario during which the mechanical pressure driving the expansion is significantly larger than the thermal pressure of the neutral medium. The second test examines the `late phase' D-type scenario during which the system relaxes to pressure equilibrium with the external medium. Although they are mutually in excellent agreement, all twelve participating codes follow a modified expansion law that deviates significantly from the classical Spitzer solution in both scenarios. We present a semi-empirical formula combining the two different solutions appropriate to both early and late phases that agrees with high-resolution simulations to $\lesssim2\%$. This formula provides a much better benchmark solution for code validation than the Spitzer solution. The present comparison has validated the participating codes and through this project we provide a dataset for calibrating the treatment of ionizing radiation hydrodynamics codes.

preprint2012arXiv

Extended Lyman-alpha emission from cold accretion streams

{Abridged} We investigate the observability of cold accretion streams at redshift 3 via Lyman-alpha (Lya) emission and the feasibility of cold accretion as the main driver of Lya blobs (LABs). We run cosmological zoom simulations focusing on 3 halos spanning two orders of magnitude in mass, roughly from 10^11 to 10^13 solar masses. We use a version of the Ramses code that includes radiative transfer of UV photons, and we employ a refinement strategy that allows us to resolve accretion streams in their natural environment to an unprecedented level. For the first time, we self-consistently model self-shielding in the cold streams from the cosmological UV background, which enables us to predict their temperatures, ionization states and Lya luminosities with improved accuracy. We find the efficiency of gravitational heating in cold streams in a ~10^11 solar mass halo is around 10-20% throughout most of the halo but reaching much higher values close to the center. As a result most of the Lya luminosity comes from gas which is concentrated at the central 20% of the halo radius, leading to Lya emission which is not extended. In more massive halos, of >10^12 solar masses, cold accretion is complex and disrupted, and gravitational heating does not happen as a steady process. Ignoring the factors of Lya scattering, local UV enhancement, and SNe feedback, cold accretion alone in these massive halos can produce LABs that largely agree with observations in terms of morphology, extent, and luminosity. Our simulations slightly and systematically over-predict LAB abundances, perhaps hinting that the interplay of these ignored factors may have a negative net effect on extent and luminosity. We predict that a factor of a few increase in sensitivity from current observational limits should unambiguously reveal continuum-free accretion streams around massive galaxies at z=3.