Researcher profile

Dominique Aubert

Dominique Aubert contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
2topics
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)

preprint2026arXiv

Higher order methods for Radiative Transfer in Astrophysical simulations: Pn vs M1

In current cosmological simulations, the radiative transfer modules generally rely on the M1 approximation, which has some glaring flaws related to its fluid-like behaviour, such as spurious pseudo-sources and loss of directionality when radiation fronts from different directions collide. Pn, another moment-based model used in other fields of physics, may correct these issues. We aim at testing out Pn in an astrophysical setting and compare it to M1, in order to see if it can indeed correct M1's flaws. Also, we want to use Pn's solutions to better pinpoint M1 errors. We implement a Pn radiation transport method and couple it to a photo-thermo-chemistry module to account for the interaction of ionising radiation with the Hydrogen gas, and benchmark it using tests for radiative transfer models comparison in astrophysics as defined in arXiv:astro-ph/0603199. We find that high order P_n (e.g. P9) indeed correct M1's flaws, while faring as well or even better in some aspects in the tests, in particular when directionality is important or colliding radiation fronts occur. By comparing P9 and M1 radiation fields in an idealised and cosmological test case, we highlight a new, thus far unreported artefact of M_1, the 'dark sombrero'. A dark sombrero appears as a spherical photon-deficit shell around the source. The photon density in dark sombreros can be underestimated by a factor up to 2-3. They occur in regions where a source's radiation field connects with that of another source or group of sources. These basic properties (position and amplitude) of the dark sombreros may depend on the sources' relative intensities, positions, spatial resolution, although we have not been able to test this in detail in this study.

preprint2023arXiv

DUSTiER (DUST in the Epoch of Reionization): dusty galaxies in cosmological radiation-hydrodynamical simulations of the Epoch of Reionization with RAMSES-CUDATON

In recent years, interstellar dust has become a crucial topic in the study of the high and very high redshift Universe. Evidence points to the existence of high dust masses in massive star forming galaxies already during the Epoch of Reionization, potentially affecting the escape of ionising photons into the intergalactic medium. Moreover, correctly estimating dust extinction at UV wavelengths is essential for precise ultra-violet luminosity function (UVLF) prediction and interpretation. In this paper, we investigate the impact of dust on the observed properties of high redshift galaxies, and cosmic reionization. To this end, we couple a physical model for dust production to the fully coupled radiation-hydrodynamics cosmological simulation code RAMSES-CUDATON, and perform a $16^3$, $2048^3$, simulation, that we call DUSTiER for DUST in the Epoch of Reionization. It yields galaxies with dust masses and UV slopes compatible with constraints at z $\geq 5$. We find that extinction has a dramatic impact on the bright end of the UVLF, even as early as $\rm z=8$, and our dusty UVLFs are in better agreement with observations than dust-less UVLFs. The fraction of obscured star formation rises up to 45% at $\rm z=5$, in agreement with some of the latest results from ALMA. Finally, we find that dust reduces the escape of ionising photons from galaxies more massive than $10^{10} M_\odot$ (brighter than $\approx -18$ MAB1600) by >10%, and possibly up to 80-90% for our most massive galaxies. Nevertheless, we find that the ionising escape fraction is first and foremost set by neutral Hydrogen in galaxies, as the latter produces transmissions up to 100 times smaller than through dust alone.

preprint2022arXiv

Reionization time of the Local Group and Local-Group-like halo pairs

Patchy cosmic reionization resulted in the ionizing UV background asynchronous rise across the Universe. The latter might have left imprints visible in present day observations. Several numerical simulation-based studies show correlations between reionization time and overdensities and object masses today. To remove the mass from the study, as it may not be the sole important parameter, this paper focuses solely on the properties of paired halos within the same mass range as the Milky Way. For this purpose, it uses CoDaII, a fully-coupled radiation hydrodynamics reionization simulation of the local Universe. This simulation holds a halo pair representing the Local Group, in addition to other pairs, sharing similar mass, mass ratio, distance separation and isolation criteria but in other environments, alongside isolated halos within the same mass range. Investigations of the paired halo reionization histories reveal a wide diversity although always inside-out given our reionization model. Within this model, halos in a close pair tend to be reionized at the same time but being in a pair does not bring to an earlier time their mean reionization. The only significant trend is found between the total energy at z = 0 of the pairs and their mean reionization time: pairs with the smallest total energy (bound) are reionized up to 50 Myr earlier than others (unbound). Above all, this study reveals the variety of reionization histories undergone by halo pairs similar to the Local Group, that of the Local Group being far from an average one. In our model, its reionization time is ~625 Myr against 660+/-4 Myr (z~8.25 against 7.87+/-0.02) on average.

preprint2022arXiv

The short ionizing photon mean free path at z=6 in Cosmic Dawn III, a new fully-coupled radiation-hydrodynamical simulation of the Epoch of Reionization

Recent determinations of the mean free path of ionising photons (mfp) in the intergalactic medium (IGM) at $\rm z=6$ are lower than many theoretical predictions. To gain insight into this issue, we investigate the evolution of the mfp in our new massive fully coupled radiation hydrodynamics cosmological simulation of reionization: Cosmic Dawn III (CoDaIII). CoDaIII&#39;s scale ($\rm 94^3 \, cMpc^3$) and resolution ($\rm 8192^3$ grid) make it particularly suitable to study the evolution of the IGM during Reionization. The simulation was performed with RAMSES-CUDATON on Summit, and used 131072 processors coupled to 24576 GPUs, making it the largest Reionization simulation, and largest RAMSES simulation ever performed. A superior agreement with global constraints on Reionization is obtained in CoDaIII over CoDaII especially for the evolution of the neutral hydrogen fraction and the cosmic photo-ionization rate, thanks to an improved calibration, later end of reionization ($\rm z=5.6$), and higher spatial resolution. Analyzing the mfp, we find that CoDaIII reproduces the most recent observations very well, from $\rm z=6$ to $\rm z=4.6$. We show that the distribution of the mfp in CoDaIII is bimodal, with short (neutral) and long (ionized) mfp modes, respectively, due to the patchiness of reionization and the co-existence of neutral versus ionized regions during Reionization. The neutral mode peaks at sub-kpc to kpc scales of mfp, while the ionized mode peak evolves from $\rm 0.1 Mpc/h$ at $\rm z=7$ to $\sim 10$ Mpc/h at $\rm z=5.2$. Computing the mfp as the average of the ionized mode provides the best match to the recent observational determinations. The distribution reduces to a single neutral (ionized) mode at $\rm z>13$ ($\rm z<5$).

preprint2021arXiv

21cm signal predictions at Cosmic Dawn and Reionization with coupled radiative-hydrodynamics

The process of heating and reionization of the Universe at high redshift links small scale structures/galaxy formation and large scale inter-galactic medium properties. Even if the first one is difficult to observe, an observation window is opening on the second one, with the promising development of current and future radio telescopes. They will permit to observe the 21cm brightness temperature global signal and fluctuations. The need of large scale simulations is therefore strong to understand the properties of the IGM that will be observed. But at the same time the urge to resolve the structures responsible of those process is important. We introduce in this study, a coupled hydro-radiative transfer simulations of the Cosmic Dawn and Reionization with a simple sub-grid star formation process developed and calibrated on the state of the art simulation CoDaII. This scheme permits to follow consistently dark matter, hydrodynamics and radiative transfer evolution&#39;s on large scales, while the sub-grid models bridges to the galaxy formation scale. We process the simulation to produce 21cm signal as close as possible to the observations.

preprint2021arXiv

A first look at the topology of reionization redshifts in models of the Epoch of the Reionization

During the EoR, the first stars and galaxies appear while creating ionized bubbles that will eventually percolate near z=6. These ionized bubbles and percolation process are nowadays under a lot of scrutiny since observations of the HI gas will be carried on in the next decade with e.g. the SKA radiotelescope. Studies of the EoR are performed on semi-analytical and fully numerical cosmological simulations to investigate e.g. the topology of the process. We analyse the topology of EoR models through regions that are under the radiative influence of ionization sources. They are associated with peak patches of reionization redshift (zreion) field, for which we measure their general properties (e.g. number, size, shape, orientation). We aim at gaining insights on the geometry of the reionization process and how it relates to the matter distribution for example. We also assess how such measurements can be used to quantify the influence of physical parameters on the reionization models or the differences between fully numerical simulations and semi-analytical models. We use the DisPerSE framework (which applies the Morse theory and the persistent homology) on different EoR scenarios through gas density and zreion maps. We find that we can distinguish between EoR models with different sources using simple analyses on the number, shape and size distributions of the peak patches. For every model, we statistically show that they are rather prolate and aligned with the gas filaments. We briefly highlight that the percolation process can be followed studying zreion fields with different persistence thresholds. We show that fully numerical EMMA simulations can be made consistent with 21cmFAST models in this topological framework as long as the source distribution is diffuse enough.

preprint2020arXiv

Constraining the second half of reionization with the Lyman-$β$ forest

We present an analysis of the evolution of the Lyman-series forest into the epoch of reionization using cosmological radiative transfer simulations in a scenario where reionization ends late. We explore models with different midpoints of reionization and gas temperatures. We find that once the simulations have been calibrated to match the mean flux of the observed Lyman-$α$ forest at $4 < z < 6$, they also naturally reproduce the distribution of effective optical depths of the Lyman-$β$ forest in this redshift range. We note that the tail of the largest optical depths that is most challenging to match corresponds to the long absorption trough of ULAS J0148+0600, which we have previously shown to be rare in our simulations. We consider the evolution of the Lyman-series forest out to higher redshifts, and show that future observations of the Lyman-$β$ forest at $z>6$ will discriminate between different reionization histories. The evolution of the Lyman-$α$ and Lyman-$γ$ forests are less promising as a tool for pushing studies of reionization to higher redshifts due to the stronger saturation and foreground contamination, respectively.

preprint2020arXiv

Galactic ionising photon budget during the Epoch of Reionisation in the Cosmic Dawn II simulation

Cosmic Dawn (&#34;CoDa&#34;) II yields the first statistically-meaningful determination of the relative contribution to reionization by galaxies of different halo mass, from a fully-coupled radiation-hydrodynamics simulation of the epoch of reionization large enough ($\sim$ 100 Mpc) to model global reionization while resolving the formation of all galactic halos above $\sim 10^8 M_\odot$. Cell transmission inside high-mass haloes is bi-modal -- ionized cells are transparent, while neutral cells absorb the photons their stars produce - and the halo escape fraction $f_{esc}$ reflects the balance of star formation rate (&#34;SFR&#34;) between these modes. The latter is increasingly prevalent at higher halo mass, driving down $f_{esc}$ (we provide analytical fits to our results), whereas halo escape luminosity, proportional to $f_{esc} \times$SFR, increases with mass. Haloes with dark matter masses within $6.10^{8} M_\odot < M_h < 3.10^{10} M_\odot$ produce $\sim 80$% of the escaping photons at z=7, when the Universe is 50% ionized, making them the main drivers of cosmic reionization. Less massive haloes, though more numerous, have low SFRs and contribute less than 10% of the photon budget then, despite their high $f_{esc}$. High mass haloes are too few and too opaque, contributing $<10$% despite their high SFRs. The dominant mass range is lower (higher) at higher (lower) redshift, as mass function and reionization advance together (e.g. at z$=8.5$, x$_{\rm HI}=0.9$, $M_h < 5.10^9 M_\odot$ haloes contributed $\sim$80%). Galaxies with UV magnitudes $M_{AB1600}$ between $-12$ and $-19$ dominated reionization between z$=6$ and 8.

preprint2020arXiv

Improved constraints on reionisation from CMB observations: A parameterisation of the kSZ effect

We show that, in the context of patchy reionisation, an accurate description of the angular power spectrum of the kinetic Sunyaev-Zel&#39;dovich (kSZ) effect is not possible with simple scaling relations between the amplitude of the spectrum and global parameters, such as the reionisation midpoint and its duration. We introduce a new parameterisation of this spectrum, based on a novel description of the power spectrum of the free electrons density contrast Pee (k,z) in terms of the reionisation global history and morphology. We directly relate features of the spectrum to the typical ionised bubble size at different stages in the process and, subsequently, to the angular scale at which the patchy kSZ power spectrum reaches its maximum. We successfully calibrated our results on a custom set of advanced radiative hydrodynamical simulations and later found our parameterisation to be a valid description of a wide range of other simulations and, therefore, reionisation physics. In the end, and as long as the global reionisation history is known, two parameters are sufficient to derive the angular power spectrum. Such an innovative framework applied to cosmic microwave background data and combined with 21cm intensity mapping will allow a first consistent detection of the amplitude and shape of the patchy kSZ signal, giving in turn access to the physics of early light sources.

preprint2019arXiv

Long troughs in the Lyman-$α$ forest below redshift 6 due to islands of neutral hydrogen

A long (110 cMpc/$h$) and deep absorption trough in the Ly$α$ forest has been observed extending down to redshift 5.5 in the spectrum of ULAS J0148+0600. Although no Ly$α$ transmission is detected, Ly$β$ spikes are present which has led to claims that the gas along this trough must be ionized. Using high resolution cosmological radiative transfer simulations in large volumes, we show that in a scenario where reionization ends late ($z \sim 5.2$), our simulations can reproduce troughs as long as observed. In this model, we find that the troughs are caused by islands of neutral hydrogen. Small ionized holes within the neutral islands allow for the transmission of Ly$β$. We have also modelled the Ly$α$ emitter population around the simulated troughs, and show that there is a deficit of Ly$α$ emitters close to the trough as is observed.