Researcher profile

Pierre Ocvirk

Pierre Ocvirk contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
3topics
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

12 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

Best Practices for Data Publication in the Astronomical Literature

We present an overview of best practices for publishing data in astronomy and astrophysics journals. These recommendations are intended as a reference for authors to help prepare and publish data in a way that will better represent and support science results, enable better data sharing, improve reproducibility, and enhance the reusability of data. Observance of these guidelines will also help to streamline the extraction, preservation, integration and cross-linking of valuable data from astrophysics literature into major astronomical databases, and consequently facilitate new modes of science discovery that will better exploit the vast quantities of panchromatic and multi-dimensional data associated with the literature. We encourage authors, journal editors, referees, and publishers to implement the best practices reviewed here, as well as related recommendations from international astronomical organizations such as the International Astronomical Union (IAU) for publication of nomenclature, data, and metadata. A convenient Checklist of Recommendations for Publishing Data in the Literature is included for authors to consult before the submission of the final version of their journal articles and associated data files. We recommend that publishers of journals in astronomy and astrophysics incorporate a link to this document in their Instructions to Authors.

preprint2022arXiv

LyC escape from SPHINX galaxies in the Epoch of Reionization

We measure escape fractions, $f_{\rm esc}$, of ionizing radiation from galaxies in the SPHINX suite of cosmological radiation-hydrodynamical simulations of reionization, resolving halos with $M_{\rm vir} \gtrapprox 7.5 \times 10^7 \ M_{\odot}$ with a minimum cell width of $\approx 10$ pc. Our new and largest $20$ co-moving Mpc wide volume contains tens of thousands of star-forming galaxies with halo masses up to a few times $10^{11} \ M_{\odot}$. The simulated galaxies agree well with observational constraints of the UV luminosity function in the Epoch of Reionization. The escape fraction fluctuates strongly in individual galaxies over timescales of a few Myrs, due to its regulation by supernova and radiation feedback, and at any given time a tiny fraction of star-forming galaxies emits a large fraction of the ionizing radiation escaping into the inter-galactic medium. Statistically, $f_{\rm esc}$ peaks in intermediate-mass, intermediate-brightness, and low-metallicity galaxies ($M_{*} \approx 10^7 \ M_{\odot}$, $M_{1500} \approx -17$, $Z\lesssim 5 \times 10^{-3} \ Z_{\odot}$), dropping strongly for lower and higher masses, brighter and dimmer galaxies, and more metal-rich galaxies. The escape fraction correlates positively with both the short-term and long-term specific star formation rate. According to SPHINX, galaxies too dim to be yet observed, with $M_{1500} \gtrapprox -17$, provide about $55$ percent of the photons contributing to reionization. The global averaged $f_{\rm esc}$ naturally decreases with decreasing redshift, as predicted by UV background models and low-redshift observations. This evolution is driven by decreasing specific star formation rates over cosmic time.

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

Scattering of Lyα Photons through the Reionizing Intergalactic Medium: I. Spectral Energy Distribution

During reionization, a fraction of galactic Ly$α$ emission is scattered in the intergalactic medium (IGM) and appears as a diffuse light extending megaparsecs from the source. We investigate how to probe the properties of the early galaxies and their surrounding IGM using this scattered light. We create a Monte Carlo algorithm to track individual photons and reproduce several test cases from previous literature. Then, we run our code on the simulated IGM of the CoDaII simulation. We find that the scattered light can leave an observable imprint on the emergent spectrum if collected over several square arcminutes. Scattering can redden the emission by increasing the path lengths of photons, but it can also make the photons bluer by upscattering them according to the peculiar motion of the scatterer. The photons emitted on the far blue side of the resonance appear more extended in both frequency and space compared to those emitted near the resonance. This provides a discriminating feature for the blueward emission, which cannot be constrained from the unscattered light coming directly from the source. The ionization state of the IGM also affects the scattered light spectrum. When the source is in a small HII region, the emission goes through more scatterings in the surrounding HI region regardless of the initial frequency and ends up more redshifted and spatially extended. This can result in a weakening of the scattered light toward high $z$ during reionization. Our results provide a framework for interpreting the scattered light to be measured by high-$z$ integral-field-unit surveys.

preprint2022arXiv

The Nature of High ${\rm [OIII]}_{\rm 88μm}$/${\rm [CII]}_{\rm 158μm}$ Galaxies in the Epoch of Reionization: Low Carbon Abundance and a Top-Heavy IMF?

ALMA observations of $z>6$ galaxies have revealed abnormally high [OIII]$_{\rm 88μm}$/[CII]$_{\rm 158μm}$ ratios and [CII]$_{\rm 158μm}$ deficits compared to local galaxies. Numerous solutions have been proposed including differences in C and O abundance ratios, observational bias, and differences in ISM properties. In order to elucidate the underlying physics that drives this high-redshift phenomenon, we employ SPHINX$^{20}$, a state-of-the-art, cosmological radiation-hydrodynamics simulation, that resolves detailed ISM properties of thousands of galaxies in the epoch of reionization. We find that the observed $z>6$ [OIII]-SFR and [CII]-SFR relations can only be reproduced when the C/O abundance ratio is $\sim8\times$ lower than Solar and the total metal production is $\sim4\times$ higher than that of a Kroupa IMF. This implies that high-redshift galaxies are potentially primarily enriched by low-metallicity core-collapse supernovae with a more top-heavy IMF. As AGB stars and type-Ia supernova begin to contribute to the galaxy metallicity, both the [CII]-SFR and [CII] luminosity functions are predicted to converge to observed values at $z\sim4.5$. While we demonstrate that ionisation parameter, LyC escape fraction, ISM gas density, and CMB attenuation all drive galaxies towards higher [OIII]/[CII], observed values at $z>6$ can only be reproduced with substantially lower C/O abundances compared to Solar. The combination of [CII] and [OIII] can be used to predict the values of ionisation parameter, ISM gas density, and LyC escape fraction and we provide estimates of these quantities for nine observed $z>6$ galaxies. Finally, we demonstrate that [OI]$_{\rm 63μm}$ can be used as a replacement for [CII] in high-redshift galaxies where [CII] is unobserved and argue that more observation time should be used to target [OI] at $z>6$. (Abridged)

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

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.

preprint2011arXiv

A signature of the internal reionisation of the Milky Way?

We present a new semi-analytical model of the population of satellite galaxies of the Milky Way, aimed at estimating the effect of the spatial structure of reionisation at galaxy scale on the properties of the satellites. In this model reionisation can be either: (A) externally-driven and uniform, or (B) internally-driven, by the most massive progenitor of the Milky-Way. In the latter scenario the propagation of the ionisation front and photon dilution introduce a delay in the photo-evaporation of the outer satellites&#39; gas with respect to the inner satellites. As a consequence, outer satellites experience a longer period of star formation than those in the inner halo. We use simple models to account for star formation, the propagation of the ionisation front, photo-evaporation and observational biases. Both scenarios yield a model satellite population at z=0 that matches the observed luminosity function and mass-to-light ratios. However, the predicted population for scenario (B) is significantly more extended spatially than for scenario (A), by about 0.3 dex in distance, resulting in a much better match to the observations. The survival of this structural signature imprinted by the local UV field during reionisation on the radial distribution of satellites makes it a promising tool for studying the reionisation epoch at galaxy scale in the Milky Way and nearby galaxies resolved in stars with forthcoming large surveys. However, more work is needed to determine how the effect reported here can be disentangled from that of cosmic variance between different realisations of Milky Way haloes.