Researcher profile

Joachim Stadel

Joachim Stadel contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

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

10 published item(s)

preprint2022arXiv

From EMBER to FIRE: predicting high resolution baryon fields from dark matter simulations with Deep Learning

Hydrodynamic simulations provide a powerful, but computationally expensive, approach to study the interplay of dark matter and baryons in cosmological structure formation. Here we introduce the EMulating Baryonic EnRichment (EMBER) Deep Learning framework to predict baryon fields based on dark-matter-only simulations thereby reducing computational cost. EMBER comprises two network architectures, U-Net and Wasserstein Generative Adversarial Networks (WGANs), to predict two-dimensional gas and HI densities from dark matter fields. We design the conditional WGANs as stochastic emulators, such that multiple target fields can be sampled from the same dark matter input. For training we combine cosmological volume and zoom-in hydrodynamical simulations from the Feedback in Realistic Environments (FIRE) project to represent a large range of scales. Our fiducial WGAN model reproduces the gas and HI power spectra within 10% accuracy down to ~10 kpc scales. Furthermore, we investigate the capability of EMBER to predict high resolution baryon fields from low resolution dark matter inputs through upsampling techniques. As a practical application, we use this methodology to emulate high-resolution HI maps for a dark matter simulation of a L=100 Mpc/h comoving cosmological box. The gas content of dark matter haloes and the HI column density distributions predicted by EMBER agree well with results of large volume cosmological simulations and abundance matching models. Our method provides a computationally efficient, stochastic emulator for augmenting dark matter only simulations with physically consistent maps of baryon fields.

preprint2021arXiv

Growing Mars fast: High-resolution GPU simulations of embryo formation

Recent high precision meteoritic data improve constraints on the formation timescale and bulk composition of the terrestrial planets. High resolution N-body simulations allow direct comparison of embryo growth timescale and accretion zones to these constraints. In this paper, we present results of high resolution simulations for embryo formation from a disc of up to 41,000 fully-self gravitating planetesimals with the GPU-based N-body code GENGA. Our results indicate that the growth of embryos are highly dependent on the initial conditions. More massive initial planetesimals, a shorter gas disc decay timescale and initially eccentric Jupiter and Saturn (EJS) all lead to faster growth of embryos. Asteroid belt material can thereby be implanted into the terrestrial planet region via sweeping secular resonances. This could possibly explain the rapid growth of Mars within 10 Myr inferred from its Hf-W chronology. The sweeping secular resonance almost completely clears the asteroid belt and deposits this material in the Mercury-Venus region, altering the composition of embryos there. This could result in embryos in the Mercury-Venus region accreting an unexpectedly high mass fraction from beyond 2 AU. Changing the initial orbits of Jupiter and Saturn to more circular (CJS) or assuming embryos formed in a gas free environment removes the sweeping secular resonance effect and thus greatly decreases material accreted from beyond 2 AU for Mercury-Venus region embryos. We therefore propose that rock samples from Mercury and Venus could aid greatly in deducing the condition and lifetime of the initial protoplanetary gas disc during planetesimal and embryo formation, as well as the initial architecture of the giant planets.

preprint2021arXiv

Parameter inference with non-linear galaxy clustering: accounting for theoretical uncertainties

We implement EuclidEmulator (version 1), an emulator for the non-linear correction of the matter power spectrum, into the MCMC forecasting code MontePython. We compare the performance of Halofit, HMCode, and EuclidEmulator1, both at the level of power spectrum prediction and at the level of posterior probability distributions of the cosmological parameters, for different cosmological models and different galaxy power spectrum wave number cut-offs. We confirm that the choice of the power spectrum predictor has a non-negligible effect on the computed sensitivities when doing cosmological parameter forecasting, even for a conservative wave number cut-off of $0.2\,h\,{\rm Mpc}^{-1}$. We find that EuclidEmulator1 is on average up to $17\%$ more sensitive to the cosmological parameters than the other two codes, with the most significant improvements being for the Hubble parameter of up to $42\%$ and the equation of state of dark energy of up to $26\%$, depending on the case. In addition, we point out that the choice of the power spectrum predictor contributes to the risk of computing a significantly biased mean cosmology when doing parameter estimations. For the four tested scenarios we find biases, averaged over the cosmological parameters, of between 0.5 and 2$σ$ (from below $1σ$ up to $6σ$ for individual parameters). This paper provides a proof of concept that this risk can be mitigated by taking a well-tailored theoretical uncertainty into account as this allows to reduce the bias by a factor of 2 to 5, depending on the case under consideration, while keeping posterior credibility contours small: the standard deviations are amplified by a factor of $\leq1.4$ in all cases.

preprint2021arXiv

The EOS/Resolution Conspiracy: Convergence in Proto-Planetary Collision Simulations

We investigate how the choice of equation of state (EOS) and resolution conspire to affect the outcomes of giant impact (GI) simulations. We focus on the simple case of equal mass collisions of two Earth-like $0.5\,M_\oplus$ proto-planets showing that the choice of EOS has a profound impact on the outcome of such collisions as well as on the numerical convergence with resolution. In simulations where the Tillotson EOS is used, impacts generate an excess amount of vapour due to the lack of a thermodynamically consistent treatment of phase transitions and mixtures. In oblique collisions this enhances the artificial angular momentum (AM) transport from the planet to the circum-planetary disc reducing the planet's rotation period over time. Even at a resolution of $1.3 \times 10^6$ particles the result is not converged. In head-on collisions the lack of a proper treatment of the solid/liquid-vapour phase transition allows the bound material to expand to very low densities which in turn results in very slow numerical convergence of the critical specific impact energy for catastrophic disruption $Q_{RD}^*$ with increasing resolution as reported in prior work. The simulations where ANEOS is used for oblique impacts are already converged at a modest resolution of $10^5$ particles, while head-on collisions converge when they evidence the post-shock formation of a dense iron-rich ring, which promotes gravitational re-accumulation of material. Once sufficient resolution is reached to resolve the liquid-vapour phase transition of iron in the ANEOS case, and this ring is resolved, the value of $Q_{RD}^*$ has then converged.

preprint2020arXiv

Baryonic effects for weak lensing. Part I. Power spectrum and covariance matrix

Baryonic feedback effects lead to a suppression of the weak lensing angular power spectrum on small scales. The poorly constrained shape and amplitude of this suppression is an important source of uncertainties for upcoming cosmological weak lensing surveys such as Euclid or LSST. In this first paper in a series of two, we use simulations to build a Euclid-like tomographic mock data-set for the cosmic shear power spectrum and the corresponding covariance matrix, which are both corrected for baryonic effects following the baryonification method of Schneider et al. (2019). In addition, we develop an emulator to obtain fast predictions of the baryonic power suppression, allowing us to perform a likelihood inference analysis for a standard $Λ$CDM cosmology with both cosmological and astrophysical parameters. Our main findings are the following: (i) ignoring baryonic effects leads to a greater than 5$σ$ bias on the cosmological parameters $Ω_m$ and $σ_8$; (ii) restricting the analysis to the largest scales, that are mostly unaffected by baryons, makes the bias disappear, but results in a blow-up of the $Ω_m$-$σ_8$ contour area by more than a factor of 10; (iii) ignoring baryonic effects on the covariance matrix does not significantly affect cosmological parameter estimates; (iv) while the baryonic suppression is mildly cosmology dependent, this effect does not noticeably modify the posterior contours. Overall, we conclude that including baryonic uncertainties in terms of nuisance parameters results in unbiased and surprisingly tight constraints on cosmology.

preprint2020arXiv

Baryonic effects for weak lensing. Part II. Combination with X-ray data and extended cosmologies

An accurate modelling of baryonic feedback effects is required to exploit the full potential of future weak-lensing surveys such as Euclid or LSST. In this second paper in a series of two, we combine Euclid-like mock data of the cosmic shear power spectrum with an eROSITA X-ray mock of the cluster gas fraction to run a combined likelihood analysis including both cosmological and baryonic parameters. Following the first paper of this series, the baryonic effects (based on the baryonic correction model of Schneider et al. 2019) are included in both the tomographic power spectrum and the covariance matrix. However, this time we assume the more realistic case of a $Λ$CDM cosmology with massive neutrinos, and we consider several extensions of the currently favoured cosmological model. For the standard $Λ$CDM case, we show that including X-ray data reduces the uncertainties on the sum of the neutrino mass by $\sim30$ percent, while there is only a mild improvement on other parameters such as $Ω_m$ and $σ_8$. As extensions of $Λ$CDM, we consider the cases of a dynamical dark energy model (wCDM), a $f(R)$ gravity model (fRCDM), and a mixed dark matter model ($Λ$MDM) with both a cold and a warm/hot dark matter component. We find that combining weak lensing with X-ray data only leads to a mild improvement of the constraints on the additional parameters of wCDM, while the improvement is more substantial for both fRCDM and $Λ$MDM. Ignoring baryonic effects in the analysis pipeline leads to significant false-detections of either phantom dark energy or a light subdominant dark matter component. Overall we conclude that for all cosmologies considered, a general parametrisation of baryonic effects is both necessary and sufficient to obtain tight constraints on cosmological parameters.

preprint2020arXiv

Could Uranus and Neptune form by collisions of planetary embryos?

The origin of Uranus and Neptune remains a challenge for planet formation models. A potential explanation is that the planets formed from a population of a few planetary embryos with masses of a few Earth masses which formed beyond Saturn's orbit and migrated inwards. These embryos can collide and merge to form Uranus and Neptune. In this work we revisit this formation scenario and study the outcomes of such collisions using 3D hydrodynamical simulations. We investigate under what conditions the perfect-merging assumption is appropriate, and infer the planets' final masses, obliquities and rotation periods, as well as the presence of proto-satellite disks. We find that the total bound mass and obliquities of the planets formed in our simulations generally agree with N-body simulations therefore validating the perfect-merging assumption. The inferred obliquities, however, are typically different from those of Uranus and Neptune, and can be roughly matched only in a few cases. In addition, we find that in most cases the planets formed in this scenario rotate faster than Uranus and Neptune, close to break-up speed, and have massive disks. We therefore conclude that forming Uranus and Neptune in this scenario is challenging, and further research is required. We suggest that future planet formation models should aim to explain the various physical properties of the planets such as their masses, compositions, obliquities, rotation rates and satellite systems.

preprint2020arXiv

Report from the Tri-Agency Cosmological Simulation Task Force

The Tri-Agency Cosmological Simulations (TACS) Task Force was formed when Program Managers from the Department of Energy (DOE), the National Aeronautics and Space Administration (NASA), and the National Science Foundation (NSF) expressed an interest in receiving input into the cosmological simulations landscape related to the upcoming DOE/NSF Vera Rubin Observatory (Rubin), NASA/ESA's Euclid, and NASA's Wide Field Infrared Survey Telescope (WFIRST). The Co-Chairs of TACS, Katrin Heitmann and Alina Kiessling, invited community scientists from the USA and Europe who are each subject matter experts and are also members of one or more of the surveys to contribute. The following report represents the input from TACS that was delivered to the Agencies in December 2018.

preprint2019arXiv

Bifurcation in the history of Uranus and Neptune: the role of giant impacts

Despite many similarities, there are significant observed differences between Uranus and Neptune: while Uranus is tilted and has a regular set of satellites, suggesting their accretion from a disk, Neptune's moons are irregular and are captured objects. In addition, Neptune seems to have an internal heat source, while Uranus is in equilibrium with solar insulation. Finally, structure models based on gravity data suggest that Uranus is more centrally condensed than Neptune. We perform a large suite of high resolution SPH simulations to investigate whether these differences can be explained by giant impacts. For Uranus, we find that an oblique impact can tilt its spin axis and eject enough material to create a disk where the regular satellites are formed. Some of the disks are massive and extended enough, and consist of enough rocky material to explain the formation of Uranus' regular satellites. For Neptune, we investigate whether a head-on collision could mix the interior, and lead to an adiabatic temperature profile, which may explain its larger flux and higher moment of inertia value. We find that massive and dense projectiles can penetrate towards the centre and deposit mass and energy in the deep interior, leading to a less centrally concentrated interior for Neptune. We conclude that the dichotomy between the ice giants can be explained by violent impacts after their formation.

preprint2019arXiv

Primordial Earth mantle heterogeneity caused by the Moon-forming giant impact

The giant impact hypothesis for Moon formation successfully explains the dynamic properties of the Earth-Moon system but remains challenged by the similarity of isotopic fingerprints of the terrestrial and lunar mantles. Moreover, recent geochemical evidence suggests that the Earth's mantle preserves ancient (or "primordial") heterogeneity that predates the Moon-forming giant impact. Using a new hydrodynamical method, we here show that Moon-forming giant impacts lead to a stratified starting condition for the evolution of the terrestrial mantle. The upper layer of the Earth is compositionally similar to the disk, out of which the Moon evolves, whereas the lower layer preserves proto-Earth characteristics. As long as this predicted compositional stratification can at least partially be preserved over the subsequent billions of years of Earth mantle convection, the compositional similarity between the Moon and the accessible Earth's mantle is a natural outcome of realistic and high-probability Moon-forming impact scenarios. The preservation of primordial heterogeneity in the modern Earth not only reconciles geochemical constraints but is also consistent with recent geophysical observations. Furthermore, for significant preservation of a proto-Earth reservoir, the bulk composition of the Earth-Moon system may be systematically shifted towards chondritic values.