Showing 1–50 of 57 results for author: Semelin, B

YOLO-CIANNA: Galaxy detection with deep learning in radio data. I. A new YOLO-inspired source detection method applied to the SKAO SDC1

Authors: D. Cornu, P. Salomé, B. Semelin, A. Marchal, J. Freundlich, S. Aicardi, X. Lu, G. Sainton, F. Mertens, F. Combes, C. Tasse

Abstract: The upcoming Square Kilometer Array (SKA) will set a new standard regarding data volume generated by an astronomical instrument, which is likely to challenge widely adopted data analysis tools that scale inadequately with the data size. This study aims to develop a new source detection and characterization method for massive radio astronomical datasets by adapting modern deep-learning object detec… ▽ More The upcoming Square Kilometer Array (SKA) will set a new standard regarding data volume generated by an astronomical instrument, which is likely to challenge widely adopted data analysis tools that scale inadequately with the data size. This study aims to develop a new source detection and characterization method for massive radio astronomical datasets by adapting modern deep-learning object detection techniques. These approaches have proved their efficiency on complex computer vision tasks, and we seek to identify their specific strengths and weaknesses when applied to astronomical data. We introduce YOLO-CIANNA, a highly customized deep-learning object detector designed specifically for astronomical datasets. This paper presents the method and describes all the low-level adaptations required to address the specific challenges of radio-astronomical images. We demonstrate this method's capabilities using simulated 2D continuum images from the SKAO SDC1 dataset. Our method outperforms every other published result on the specific SDC1 dataset. Using the SDC1 metric, we improve the challenge-winning score by +139\% and the score of the only other post-challenge participation by +61\%. Our catalog has a detection purity of 94\% while detecting 40 to 60 \% more sources than previous top-score results. The trained model can also be forced to reach 99\% purity in post-process and still detect 10 to 30\% more sources than the other top-score methods. It is also capable of real-time detection, with a peak prediction speed of 500 images of 512x512 pixels per second on a single GPU. YOLO-CIANNA achieves state-of-the-art detection and characterization results on the simulated SDC1 dataset. The method is open source and included in the wider CIANNA framework. We provide scripts to train and apply this method to the SDC1 dataset in the CIANNA repository. △ Less

Submitted 8 February, 2024; originally announced February 2024.

Comments: 40 pages

First upper limits on the 21 cm signal power spectrum from cosmic dawn from one night of observations with NenuFAR

Authors: S. Munshi, F. G. Mertens, L. V. E. Koopmans, A. R. Offringa, B. Semelin, D. Aubert, R. Barkana, A. Bracco, S. A. Brackenhoff, B. Cecconi, E. Ceccotti, S. Corbel, A. Fialkov, B. K. Gehlot, R. Ghara, J. N. Girard, J. M. Grießmeier, C. Höfer, I. Hothi, R. Mériot, M. Mevius, P. Ocvirk, A. K. Shaw, G. Theureau, S. Yatawatta , et al. (2 additional authors not shown)

Abstract: The redshifted 21 cm signal from neutral hydrogen is a direct probe of the physics of the early universe and has been an important science driver of many present and upcoming radio interferometers. In this study we use a single night of observations with the New Extension in Nançay Upgrading LOFAR (NenuFAR) to place upper limits on the 21 cm power spectrum from cosmic dawn at a redshift of $z$ = 2… ▽ More The redshifted 21 cm signal from neutral hydrogen is a direct probe of the physics of the early universe and has been an important science driver of many present and upcoming radio interferometers. In this study we use a single night of observations with the New Extension in Nançay Upgrading LOFAR (NenuFAR) to place upper limits on the 21 cm power spectrum from cosmic dawn at a redshift of $z$ = 20.3. NenuFAR is a new low-frequency radio interferometer, operating in the 10-85 MHz frequency range, currently under construction at the Nançay Radio Observatory in France. It is a phased array instrument with a very dense uv coverage at short baselines, making it one of the most sensitive instruments for 21 cm cosmology analyses at these frequencies. Our analysis adopts the foreground subtraction approach, in which sky sources are modeled and subtracted through calibration and residual foregrounds are subsequently removed using Gaussian process regression. The final power spectra are constructed from the gridded residual data cubes in the uv plane. Signal injection tests are performed at each step of the analysis pipeline, the relevant pipeline settings are optimized to ensure minimal signal loss, and any signal suppression is accounted for through a bias correction on our final upper limits. We obtain a best 2$σ$ upper limit of $2.4\times 10^7$ $\text{mK}^{2}$ at $z$ = 20.3 and $k$ = 0.041 $h\,\text{cMpc}^{-1}$. We see a strong excess power in the data, making our upper limits two orders of magnitude higher than the thermal noise limit. We investigate the origin and nature of this excess power and discuss further improvements to the analysis pipeline that can potentially mitigate it and consequently allow us to reach thermal noise sensitivity when multiple nights of observations are processed in the future. △ Less

Submitted 30 April, 2024; v1 submitted 9 November, 2023; originally announced November 2023.

Comments: 27 pages, 21 figures, and 6 tables; accepted for publication in Astronomy and Astrophysics (A&A); language edits implemented; typos corrected

Journal ref: A&A 681, A62 (2024)

Wavelet Based Statistics for Enhanced 21cm EoR Parameter Constraints

Authors: Ian Hothi, Erwan Allys, Benoit Semelin, Francois Boulanger

Abstract: We propose a new approach to improve the precision of astrophysical parameter constraints for the 21cm signal from the epoch of reionisation (EoR). Our method introduces new sets of summary statistics, hereafter `evolution compressed' statistics, which quantify the spectral evolution of the 2D spatial statistics computed at fixed redshift. We defined such compressed statistics for power spectrum (… ▽ More We propose a new approach to improve the precision of astrophysical parameter constraints for the 21cm signal from the epoch of reionisation (EoR). Our method introduces new sets of summary statistics, hereafter `evolution compressed' statistics, which quantify the spectral evolution of the 2D spatial statistics computed at fixed redshift. We defined such compressed statistics for power spectrum (PS), wavelet scattering transforms (WST), and wavelet moments (WM), which also characterise non-Gaussian features. To compare these different statistics, along with the 3D power spectrum, we estimated their Fisher information on three cosmological parameters from an ensemble of simulations of 21cm EoR data, both in noiseless and noisy scenarios using Square Kilometre Array (SKA) noise levels equivalent to 100 and 1000 hours of observations. We also compare wavelet statistics, in particular WST, built from standard directional Morlet wavelets, as well as from a set of isotropic wavelets derived from the binning window function of the 2D power spectrum. For the noiseless case, the compressed wavelet statistics give constraints that are up to five times more precise than those obtained from the 3D isotropic power spectrum. At the same time, for 100h SKA noise, from which it is difficult to extract non-Gaussian features, compressed wavelet statistics still give over 30\% tighter constraints. We find that the wavelet statistics with wavelets derived from the power-spectrum binning window function provide the tightest constraints of all the statistics, with the WSTs seemingly performing better than the WMs, in particular when working with noisy data. The findings of this study demonstrate that evolution-compressed statistics extract more information than usual 3D isotropic power-spectra approaches and that our wavelet-based statistics can consistently outmatch power-spectrum-based statistics. △ Less

Submitted 5 April, 2024; v1 submitted 31 October, 2023; originally announced November 2023.

Comments: 17 pages, 8 figures. Accepted to A&A. Difference to version 1 includes the addition of the same wavelet used for the Wavelet Moments application within the Wavelet Scattering Transform formalism, as well as better normalising of the scattering coefficients

Journal ref: A&A 686, A212 (2024)

The LoReLi database: 21 cm signal inference with 3D radiative hydrodynamics simulations

Authors: Romain Meriot, Benoit Semelin

Abstract: The Square Kilometer array is expected to measure the 21cm signal from the Epoch of Reionization (EoR) in the coming decade, and its pathfinders may provide a statistical detection even earlier. The currently reported upper limits provide tentative constraints on the astrophysical parameters of the models of the EoR. In order to interpret such data with 3D radiative hydrodynamics simulations usi… ▽ More The Square Kilometer array is expected to measure the 21cm signal from the Epoch of Reionization (EoR) in the coming decade, and its pathfinders may provide a statistical detection even earlier. The currently reported upper limits provide tentative constraints on the astrophysical parameters of the models of the EoR. In order to interpret such data with 3D radiative hydrodynamics simulations using Bayesian inference, we present the latest developments of the \textsc{Licorice} code. Relying on an implementation of the halo conditional mass function to account for unresolved star formation, this code now allows accurate simulations of the EoR at $256^3$ resolution. We use this version of \textsc{Licorice} to produce the first iteration of \textsc{LoReLi}, a public dataset now containing hundreds of 21cm signals computed from radiative hydrodynamics simulations. We train a neural network on \textsc{LoReLi} to provide a fast emulator of the \textsc{Licorice} power spectra, \textsc{LorEMU}, which has $\sim 5\%$ rms error relative to the simulated signals. \textsc{LorEMU} is used in a Markov Chain Monte Carlo framework to perform Bayesian inference, first on a mock observation composed of a simulated signal and thermal noise corresponding to 100h observations with the SKA. We then apply our inference pipeline to the latest measurements from the HERA interferometer. We report constraints on the X-ray emissivity, and confirm that cold reionization scenarios are unlikely to accurately represent our Universe. △ Less

Submitted 1 December, 2023; v1 submitted 4 October, 2023; originally announced October 2023.

Comments: 15 figures, 15 pages. Accepted by A&A

SKA Science Data Challenge 2: analysis and results

Authors: P. Hartley, A. Bonaldi, R. Braun, J. N. H. S. Aditya, S. Aicardi, L. Alegre, A. Chakraborty, X. Chen, S. Choudhuri, A. O. Clarke, J. Coles, J. S. Collinson, D. Cornu, L. Darriba, M. Delli Veneri, J. Forbrich, B. Fraga, A. Galan, J. Garrido, F. Gubanov, H. Håkansson, M. J. Hardcastle, C. Heneka, D. Herranz, K. M. Hess , et al. (83 additional authors not shown)

Abstract: The Square Kilometre Array Observatory (SKAO) will explore the radio sky to new depths in order to conduct transformational science. SKAO data products made available to astronomers will be correspondingly large and complex, requiring the application of advanced analysis techniques to extract key science findings. To this end, SKAO is conducting a series of Science Data Challenges, each designed t… ▽ More The Square Kilometre Array Observatory (SKAO) will explore the radio sky to new depths in order to conduct transformational science. SKAO data products made available to astronomers will be correspondingly large and complex, requiring the application of advanced analysis techniques to extract key science findings. To this end, SKAO is conducting a series of Science Data Challenges, each designed to familiarise the scientific community with SKAO data and to drive the development of new analysis techniques. We present the results from Science Data Challenge 2 (SDC2), which invited participants to find and characterise 233245 neutral hydrogen (Hi) sources in a simulated data product representing a 2000~h SKA MID spectral line observation from redshifts 0.25 to 0.5. Through the generous support of eight international supercomputing facilities, participants were able to undertake the Challenge using dedicated computational resources. Alongside the main challenge, `reproducibility awards' were made in recognition of those pipelines which demonstrated Open Science best practice. The Challenge saw over 100 participants develop a range of new and existing techniques, with results that highlight the strengths of multidisciplinary and collaborative effort. The winning strategy -- which combined predictions from two independent machine learning techniques to yield a 20 percent improvement in overall performance -- underscores one of the main Challenge outcomes: that of method complementarity. It is likely that the combination of methods in a so-called ensemble approach will be key to exploiting very large astronomical datasets. △ Less

Submitted 14 March, 2023; originally announced March 2023.

Comments: Under review by MNRAS; 28 pages, 16 figures

Studying the Multi-frequency Angular Power Spectrum of the Cosmic Dawn 21-cm Signal

Authors: Abinash Kumar Shaw, Raghunath Ghara, Saleem Zaroubi, Rajesh Mondal, Garrelt Mellema, Florent Mertens, Léon V. E. Koopmans, Benoît Semelin

Abstract: The light-cone (LC) anisotropy arises due to cosmic evolution of the cosmic dawn 21-cm signal along the line-of-sight (LoS) axis of the observation volume. The LC effect makes the signal statistically non-ergodic along the LoS axis. The multi-frequency angular power spectrum (MAPS) provides an unbiased alternative to the popular 3D power spectrum as it does not assume statistical ergodicity along… ▽ More The light-cone (LC) anisotropy arises due to cosmic evolution of the cosmic dawn 21-cm signal along the line-of-sight (LoS) axis of the observation volume. The LC effect makes the signal statistically non-ergodic along the LoS axis. The multi-frequency angular power spectrum (MAPS) provides an unbiased alternative to the popular 3D power spectrum as it does not assume statistical ergodicity along every direction in the signal volume. Unlike the 3D power spectrum which mixes the cosmic evolution of the 21-cm signal along the LoS $k$ modes, MAPS keeps the evolution information disentangled. Here we first study the impact of different underlying physical processes during cosmic dawn on the behaviour of the 21-cm MAPS using simulations of various different scenarios and models. We also make error predictions in 21-cm MAPS measurements considering only the system noise and cosmic variance for mock observations of HERA, NenuFAR and SKA-Low. We find that $100~{\rm h}$ of HERA observations will be able to measure 21-cm MAPS at $\geq 3σ$ for $\ell \lesssim 1000$ with $0.1\,{\rm MHz}$ channel-width. The better sensitivity of SKA-Low allows reaching this sensitivity up to $\ell \lesssim 3000$. Note that due to the difference in the frequency coverage of the various experiements, the CD-EoR model considered for NenuFAR is different than those used for the HERA and SKA-Low predictions. Considering NenuFAR with the new model, measurements $\geq 2σ$ are possible only for $\ell \lesssim 600$ with $0.2\,{\rm MHz}$ channel-width and for a ten times longer observation time of $t_{\rm obs} = 1000~{\rm h}$. However, for the range $300 \lesssim \ell \lesssim 600$ and $t_{\rm obs}=1000~{\rm h}$ more than $3σ$ measurements are still possible for NenuFAR when combining consecutive frequency channels within a $5 ~{\rm MHz}$ band. △ Less

Submitted 13 April, 2023; v1 submitted 2 February, 2023; originally announced February 2023.

Comments: 21 pages, 12 figures, 3 tables, accepted for publication in MNRAS

Accurate modelling of the Lyman-$α$ coupling for the 21-cm signal, observability with NenuFAR and SKA

Authors: Benoit Semelin, Romain Mériot, Florent Mertens, Léon V. E. Koopmans, Dominique Aubert, Rennan Barkana, Anastasia Fialkov, Satyapan Munshi, Pierre Ocvirk

Abstract: The measurement of the $21$ cm signal from the Cosmic Dawn is a major goal for several existing and upcoming radio interferometers such as NenuFAR and the SKA. During this era before the beginning of the Epoch of Reionization, the signal is more difficult to observe due to brighter foregrounds but reveals additional information on the underlying astrophysical processes encoded in the spatial fluct… ▽ More The measurement of the $21$ cm signal from the Cosmic Dawn is a major goal for several existing and upcoming radio interferometers such as NenuFAR and the SKA. During this era before the beginning of the Epoch of Reionization, the signal is more difficult to observe due to brighter foregrounds but reveals additional information on the underlying astrophysical processes encoded in the spatial fluctuations of the spin temperature of hydrogen. To interpret future measurements, controlling the level of accuracy of the Lyman-$α$ flux modelling is mandatory. In this work, we evaluate the impact of various approximations that exist in the main fast modelling approach compared to the results of a costly full radiative transfer simulation. The fast SPINTER code, presented in this work, computes the Lyman-$α$ flux including the effect of wing scatterings for an inhomogeneous emissivity field, but assuming an otherwise homogeneous expanding universe. The LICORICE code computes the full radiative transfer in the Lyman-$α$ line without any substantial approximation. We find that the difference between homogeneous and inhomogeneous gas density and temperature is very small for the computed flux. On the contrary, neglecting the effect of gas velocities produces a significant change in the computed flux. We identify the causes (mainly Doppler shifts due to velocity gradients) and quantify the magnitude of the effect in both an idealised setup and a realistic cosmological situation. We find that the amplitude of the effect, up to a factor of $\sim 2$ on the $21$ cm signal power spectrum on some scales (depending on both other model parameters and the redshift), can be easily discriminated with an SKA-like survey and already be approached, particularly for exotic signals, by the ongoing NenuFAR Cosmic Dawn Key Science Program. △ Less

Submitted 30 January, 2023; originally announced January 2023.

Comments: 14 pages, accepted for publication in A&A

Journal ref: A&A 672, A162 (2023)

A bubble size distribution model for the Epoch of Reionization

Authors: Aristide Doussot, Benoit Semelin

Abstract: The bubble size distribution is a summary statistics that can be computed from the observed 21-cm signal from the Epoch of Reionization. As it depends only on the ionization field and is not limited to gaussian information, it is an interesting probe, complementary to the power spectrum of the full 21-cm signal. Devising a flexible and reliable theoretical model for the bubble size distribution pa… ▽ More The bubble size distribution is a summary statistics that can be computed from the observed 21-cm signal from the Epoch of Reionization. As it depends only on the ionization field and is not limited to gaussian information, it is an interesting probe, complementary to the power spectrum of the full 21-cm signal. Devising a flexible and reliable theoretical model for the bubble size distribution paves the way for using it for astrophysical parameters inference. The proposed model is built from the excursion set theory and a functional relation between the bubble volume and the collapsed mass in the bubble. Unlike previous models it accommodates any functional relation or distributions. Using parameterized relations allows us to test the predictive power of the model by performing a minimization best-fit to the bubble size distribution obtained from a high resolution, fully coupled radiative hydrodynamics simulations, HIRRAH-21. Our model is able to provide a better fit to the numerical bubble size distribution at ionization fraction of $x_{\text{H}_{\text{II}}} \sim 1\%$ and $3\%$ than other existing models. Moreover, the bubble volume to collapsed mass relation corresponding to the best-fit parameters, which is not an observable, is compared to numerical simulation data. A good match is obtained, confirming the possibility to infer this relation from an observed bubble size distribution using our model. Finally we present a simple algorithm that empirically implements the process of percolation. We show that it extends the usability of our bubble size distribution model up to $x_{\text{H}_{\text{II}}} \sim 30\%$. △ Less

Submitted 30 August, 2022; originally announced August 2022.

Comments: 15 pages, 11 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 667, A118 (2022)

Exploring the Cosmic Dawn with NenuFAR

Authors: F. G. Mertens, B. Semelin, L. V. E. Koopmans

Abstract: The exploration of the Cosmic Dawn, the period of the Universe during which the first stars and galaxies were formed, is one of the last frontiers of modern astronomy and cosmology. The redshifted 21-cm line emission from neutral hydrogen is a unique probe and can open this era for astrophysical and cosmological studies. The tentative detection of the 21-cm global signal by the EDGES team at… ▽ More The exploration of the Cosmic Dawn, the period of the Universe during which the first stars and galaxies were formed, is one of the last frontiers of modern astronomy and cosmology. The redshifted 21-cm line emission from neutral hydrogen is a unique probe and can open this era for astrophysical and cosmological studies. The tentative detection of the 21-cm global signal by the EDGES team at $z \sim 17$ underlines the need for an interferometric detection of this signal to discriminate between the numerous models trying to explain this unexpected discovery. The NenuFAR Cosmic-Dawn Key-Science Program aims to perform this detection in the redshift range $z \sim 15 - 31$ with a novel SKA precursor, the NenuFAR radio telescope located at the Station de Radioastronomie de Nançay and that started operating in 2019. Due to its compactness it is particularly sensitive to the large scale of the 21-cm signal. Only 100 hours of observation are needed to reach the level of the most extreme models, while 1000 hours are needed for the more standard models. Observations have already started, accumulating to almost 500 hours on the North Celestial Pole field. In this contribution, we introduce the project, our first results and the developments in calibration and RFI mitigation specific to this new instrument. △ Less

Submitted 21 September, 2021; originally announced September 2021.

Comments: 4 pages, SF2A 2021 proceeding

Peering into the Dark (Ages) with Low-Frequency Space Interferometers

Authors: Leon Koopmans, Rennan Barkana, Mark Bentum, Gianni Bernardi, Albert-Jan Boonstra, Judd Bowman, Jack Burns, Xuelei Chen, Abhirup Datta, Heino Falcke, Anastasia Fialkov, Bharat Gehlot, Leonid Gurvits, Vibor Jelić, Marc Klein-Wolt, Léon Koopmans, Joseph Lazio, Daan Meerburg, Garrelt Mellema, Florent Mertens, Andrei Mesinger, André Offringa, Jonathan Pritchard, Benoit Semelin, Ravi Subrahmanyan , et al. (6 additional authors not shown)

Abstract: Neutral hydrogen pervades the infant Universe, and its redshifted 21-cm signal allows one to chart the Universe. This signal allows one to probe astrophysical processes such as the formation of the first stars, galaxies, (super)massive black holes and enrichment of the pristine gas from z~6 to z~30, as well as fundamental physics related to gravity, dark matter, dark energy and particle physics at… ▽ More Neutral hydrogen pervades the infant Universe, and its redshifted 21-cm signal allows one to chart the Universe. This signal allows one to probe astrophysical processes such as the formation of the first stars, galaxies, (super)massive black holes and enrichment of the pristine gas from z~6 to z~30, as well as fundamental physics related to gravity, dark matter, dark energy and particle physics at redshifts beyond that. As one enters the Dark Ages (z>30), the Universe becomes pristine. Ground-based low-frequency radio telescopes aim to detect the spatial fluctuations of the 21-cm signal. Complementary, global 21-cm experiments aim to measure the sky-averaged 21-cm signal. Escaping RFI and the ionosphere has motivated space-based missions, such as the Dutch-Chinese NCLE instrument (currently in lunar L2), the proposed US-driven lunar or space-based instruments DAPPER and FARSIDE, the lunar-orbit interferometer DSL (China), and PRATUSH (India). To push beyond the current z~25 frontier, though, and measure both the global and spatial fluctuations (power-spectra/tomography) of the 21-cm signal, low-frequency (1-100MHz; BW~50MHz; z>13) space-based interferometers with vast scalable collecting areas (1-10-100 km2), large filling factors (~1) and large fields-of-view (4pi sr.) are needed over a mission lifetime of >5 years. In this ESA White Paper, we argue for the development of new technologies enabling interferometers to be deployed, in space (e.g. Earth-Sun L2) or in the lunar vicinity (e.g. surface, orbit or Earth-Moon L2), to target this 21-cm signal. This places them in a stable environment beyond the reach of most RFI from Earth and its ionospheric corruptions, enabling them to probe the Dark Ages as well as the Cosmic Dawn, and allowing one to investigate new (astro)physics that is inaccessible in any other way in the coming decades. [Abridged] △ Less

Submitted 12 August, 2019; originally announced August 2019.

Comments: 24 pages, 8 figures, 3 tables; White paper submitted to ESA Voyage 2050

Improved supervised learning methods for EoR parameters reconstruction

Authors: Aristide Doussot, Evan Eames, Benoit Semelin

Abstract: Within the next few years, the Square Kilometer Array (SKA) or one of its pathfinders will hopefully provide a detection of the 21-cm signal fluctuations from the Epoch of Reionization (EoR). Then, the main goal will be to accurately constrain the underlying astrophysical parameters. Currently, this is mainly done with Bayesian inference using Markov Chain Monte Carlo sampling. Recently, studies u… ▽ More Within the next few years, the Square Kilometer Array (SKA) or one of its pathfinders will hopefully provide a detection of the 21-cm signal fluctuations from the Epoch of Reionization (EoR). Then, the main goal will be to accurately constrain the underlying astrophysical parameters. Currently, this is mainly done with Bayesian inference using Markov Chain Monte Carlo sampling. Recently, studies using neural networks trained to performed inverse modelling have shown interesting results. We build on these by improving the accuracy of the predictions using neural network and exploring other supervised learning methods: the kernel and ridge regressions. Based on a large training set of 21-cm power spectra, we compare the performances of these supervised learning methods. When using an un-noised signal as input, we improve on previous neural network accuracy by one order of magnitude and, using local ridge kernel regression, we gain another factor of a few. We then reach a rms prediction error of a few percents of the 1-sigma confidence level due to SKA thermal noise (as estimated with Bayesian inference). This last performance level requires optimizing the hyper-parameters of the method: how to do that perfectly in the case of an unknown signal remains an open question. For an input signal altered by a SKA-type thermal noise, our neural network recovers the astrophysical parameter values with an error within half of the 1$σ$ confidence level due to the SKA thermal noise. This accuracy improves to 10$\%$ of the 1$σ$ level when using the local ridge kernel regression (with optimized hyper-parameters). We are thus reaching a performance level where supervised learning methods are a viable alternative to determine the best-fit parameters values. △ Less

Submitted 8 April, 2019; originally announced April 2019.

Comments: 14 pages, 6 figures, 4 tables, submitted to MNRAS

Working towards an optimal sampling of the 21 cm signal parameter space

Authors: Evan Eames, Aristide Doussot, Benoît Semelin

Abstract: With a statistical detection of the 21 cm signal fluctuations from the Epoch of Reionization (EoR) expected in the next few years, there is an interest in developing robust and precise techniques to constrain the underlying astrophysical parameters. Bayesian inference with Markov Chain Monte Carlo, or different types of supervised learning for backward modelling (from signal to parameters) are exa… ▽ More With a statistical detection of the 21 cm signal fluctuations from the Epoch of Reionization (EoR) expected in the next few years, there is an interest in developing robust and precise techniques to constrain the underlying astrophysical parameters. Bayesian inference with Markov Chain Monte Carlo, or different types of supervised learning for backward modelling (from signal to parameters) are examples of such techniques. They usually require many instances of forward modelling (from parameters to signal) in sampling the parameters space, either when performing the steps of the Markov Chain or when building a training sample for supervised learning. As forward modelling can be costly (if performed with numerical simulations for example), we should attempt to perform an optimal sampling according to some principle. With this goal in mind, we present an approach based on defining a metric on the space of observables, induced by the manner through which the modelling creates a mapping from the parameter space onto the space of observables. This metric bears a close connection to Jeffreys' prior from information theory. It is used to generate a homogeneous and isotropic sampling of the signal space with two different methods. We show that when the resulting optimized samplings, created with 21cmFAST, are used to train a neural network we obtain a modest reduction of the error on parameter reconstruction of ~10% (compared to a naïve sampling of the same size). Excluding the borders of the parameter space region, the improvement is more substantial, on the order of 30-40%. △ Less

Submitted 29 November, 2018; originally announced November 2018.

Comments: Submitted to MNRAS

Imprints of quasar duty cycle on the 21cm signal from the Epoch of Reionization

Authors: Florian Bolgar, Evan Eames, Clément Hottier, Benoit Semelin

Abstract: Quasars contribute to the 21-cm signal from the Epoch of Reionization (EoR) primarily through their ionizing UV and X-ray emission. However, their radio continuum and Lyman-band emission also regulates the 21-cm signal in their direct environment, potentially leaving the imprint of their duty cycle. We develop a model for the radio and UV luminosity functions of quasars from the EoR, and constra… ▽ More Quasars contribute to the 21-cm signal from the Epoch of Reionization (EoR) primarily through their ionizing UV and X-ray emission. However, their radio continuum and Lyman-band emission also regulates the 21-cm signal in their direct environment, potentially leaving the imprint of their duty cycle. We develop a model for the radio and UV luminosity functions of quasars from the EoR, and constrain it using recent observations. Our model is consistent with the z~7.5 quasar from Banados et al 2017, and also predicts only a few quasars suitable for 21-cm forest observations (10mJy) in the sky. We exhibit a new effect on the 21-cm signal observed against the CMB: a radio-loud quasar can leave the imprint of its duty cycle on the 21-cm tomography. We apply this effect in a cosmological simulation and conclude that the effect of typical radio-loud quasars is most likely negligible in an SKA field of view. For a 1-10mJy quasar the effect is stronger though hardly observable at SKA resolution. Then we study the contribution of the lyman band Ly-alpha to Ly-beta) emission of quasars to the Wouthuisen-Field coupling. The collective effect of quasars on the 21-cm power spectrum is larger than the thermal noise at low k, though featureless. However, a distinctive pattern around the brightest quasars in an SKA field of view may be observable in the tomography, encoding the duration of their duty cycle. This pattern has a high signal-to-noise ratio for the brightest quasar in a typical SKA shallow survey. △ Less

Submitted 16 February, 2018; originally announced February 2018.

Comments: 14 pages, 10 figures

SEURAT: SPH scheme extended with ultraviolet line radiative transfer

Authors: Makito Abe, Hiroyuki Suzuki, Kenji Hasegawa, Benoit Semelin, Hidenobu Yajima, Masayuki Umemura

Abstract: We present a novel Lyman alpha (Ly$α$) radiative transfer code, SEURAT, where line scatterings are solved adaptively with the resolution of the smoothed particle hydrodynamics (SPH). The radiative transfer method implemented in SEURAT is based on a Monte Carlo algorithm in which the scattering and absorption by dust are also incorporated. We perform standard test calculations to verify the validit… ▽ More We present a novel Lyman alpha (Ly$α$) radiative transfer code, SEURAT, where line scatterings are solved adaptively with the resolution of the smoothed particle hydrodynamics (SPH). The radiative transfer method implemented in SEURAT is based on a Monte Carlo algorithm in which the scattering and absorption by dust are also incorporated. We perform standard test calculations to verify the validity of the code; (i) emergent spectra from a static uniform sphere, (ii) emergent spectra from an expanding uniform sphere, and (iii) escape fraction from a dusty slab. Thereby we demonstrate that our code solves the Ly$α$ radiative transfer with sufficient accuracy. We emphasise that SEURAT can treat the transfer of Ly$α$ photons even in highly complex systems that have significantly inhomogeneous density fields. The high adaptivity of SEURAT is desirable to solve the propagation of Ly$α$ photons in the interstellar medium of young star-forming galaxies like Ly$α$ emitters (LAEs). Thus, SEURAT provides a powerful tool to model the emergent spectra of Ly$α$ emission, which can be compared to the observations of LAEs. △ Less

Submitted 28 January, 2018; originally announced January 2018.

Comments: 11 pages, 10 figures, accepted for publication in MNRAS

French SKA White Book - The French Community towards the Square Kilometre Array

Authors: F. Acero, J. -T. Acquaviva, R. Adam, N. Aghanim, M. Allen, M. Alves, R. Ammanouil, R. Ansari, A. Araudo, E. Armengaud, B. Ascaso, E. Athanassoula, D. Aubert, S. Babak, A. Bacmann, A. Banday, K. Barriere, F. Bellossi, J. -P. Bernard, M. G. Bernardini, M. Béthermin, E. Blanc, L. Blanchet, J. Bobin, S. Boissier , et al. (153 additional authors not shown)

Abstract: The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive p… ▽ More The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration. △ Less

Submitted 28 March, 2018; v1 submitted 19 December, 2017; originally announced December 2017.

Comments: Editor in chief: C. Ferrari; Editors: M. Alves, S. Bosse, S. Corbel, A. Ferrari, K. Ferrière, S. Gauffre, E. Josselin, G. Lagache, S. Lambert, G. Marquette, J.-M. Martin, M.-A. Miville-Deschênes, L. Montier, B. Semelin, G. Theureau, S. Vergani, N. Vilmer, P. Zarka; Original file with high resolution figures at SKA-France link: https://ska-france.oca.eu/images/SKA-France-Media/FWB_051017.pdf

What the Milky Way bulge reveals about the initial metallicity gradients in the disc

Authors: F. Fragkoudi, P. Di Matteo, M. Haywood, S. Khoperskov, A. Gómez, M. Schultheis, F. Combes, B. Semelin

Abstract: We examine the metallicity trends in the Milky Way (MW) bulge - using APOGEE DR13 data - and explore their origin by comparing two N-body models of isolated galaxies which develop a bar and a boxy/peanut (b/p) bulge. Both models have been proposed as scenarios for reconciling a disc origin of the MW bulge with a negative vertical metallicity gradient. The first model is a superposition of co-spati… ▽ More We examine the metallicity trends in the Milky Way (MW) bulge - using APOGEE DR13 data - and explore their origin by comparing two N-body models of isolated galaxies which develop a bar and a boxy/peanut (b/p) bulge. Both models have been proposed as scenarios for reconciling a disc origin of the MW bulge with a negative vertical metallicity gradient. The first model is a superposition of co-spatial, i.e. overlapping, disc populations with different scaleheights, kinematics and metallicities. In this model the thick, metal-poor, and centrally concentrated disc populations contribute significantly to the stellar mass budget in the inner galaxy. The second model is a single disc with an initial steep radial metallicity gradient, which is mapped by the bar into the b/p bulge in such a way that the vertical metallicity gradient of the MW bulge is reproduced -- as shown already in previous works in the literature. However, as we show here, the latter model does not reproduce the positive longitudinal metallicity gradient of the inner disc, nor the metal-poor innermost regions of the Bulge seen in the data. On the other hand, the model with co-spatial thin and thick disc populations reproduces all the aforementioned trends. We therefore see that it is possible to reconcile a (primarily) disc origin for the MW bulge with the observed trends in metallicity by mapping the inner thin and thick discs of the MW into a b/p. For this scenario to reproduce the observations, the $α$-enhanced, metal-poor, thick disc populations must have a significant mass contribution in the inner regions -- as has been suggested for the Milky Way. △ Less

Submitted 31 October, 2017; v1 submitted 18 October, 2017; originally announced October 2017.

Comments: 6 pages, accepted for publication in A&A Letters

Journal ref: A&A 607, L4 (2017)

21SSD: a public database of simulated 21-cm signals from the epoch of reionization

Authors: Benoit Semelin, Evan Eames, Florian Bolgar, Michel Caillat

Abstract: The 21-cm signal from the Epoch of Reionization (EoR) is expected to be detected in the next few years, either with existing instruments or by the upcoming SKA and HERA projects. In this context there is a pressing need for publicly available high-quality templates covering a wide range of possible signals. These are needed both for end-to-end simulations of the up-coming instruments, as well as t… ▽ More The 21-cm signal from the Epoch of Reionization (EoR) is expected to be detected in the next few years, either with existing instruments or by the upcoming SKA and HERA projects. In this context there is a pressing need for publicly available high-quality templates covering a wide range of possible signals. These are needed both for end-to-end simulations of the up-coming instruments, as well as to develop signal analysis methods. In this work we present such a set of templates, publicly available, for download at https://21ssd.obspm.fr/. The database contains 21-cm brightness temperature lightcones at high and low resolution, and several derived statistical quantities for 45 models spanning our choice of 3D parameter space. These data are the result of fully coupled radiative hydrodynamic high resolution ($1024^3$) simulations performed with the LICORICE code. Both X-ray and Lyman line transfer is performed to account for heating and Wouthuysen-Field coupling fluctuations. We also present a first exploitation of the data using the power spectrum and the Pixel Distribution Function (PDF) as functions of redshifts, computed from lightcone data. We analyse how these two quantities behave when varying the model parameters while taking into account the thermal noise expected of a typical SKA survey. Finally, we show that the power spectrum and the PDF have different -- and to some extent complementary -- abilities to distinguish between different models. This opens the door to formulating an optimal sampling of the parameter space, dependant on the chosen diagnostics. △ Less

Submitted 7 July, 2017; originally announced July 2017.

Comments: 14 pages,12 figures. Submitted to MNRAS

Bars & boxy/peanut bulges in thin & thick discs: I. Morphology and line-of-sight velocities of a fiducial model

Authors: F. Fragkoudi, P. Di Matteo, M. Haywood, A. Gómez, F. Combes, D. Katz, B. Semelin

Abstract: We explore trends in the morphology and line-of-sight (los) velocity of stellar populations in the inner regions of disc galaxies, using N-body simulations with both a thin (kinematically cold) and a thick (kinematically hot) disc which form a bar and boxy/peanut (b/p) bulge. The bar in the thin disc component is $\sim$50\% stronger than the thick disc bar and is more elongated, with an axis ratio… ▽ More We explore trends in the morphology and line-of-sight (los) velocity of stellar populations in the inner regions of disc galaxies, using N-body simulations with both a thin (kinematically cold) and a thick (kinematically hot) disc which form a bar and boxy/peanut (b/p) bulge. The bar in the thin disc component is $\sim$50\% stronger than the thick disc bar and is more elongated, with an axis ratio almost half that of the thick disc bar. The thin disc b/p bulge has a pronounced X-shape, while the thick disc b/p is weaker with a rather boxy shape. This leads to the signature of the b/p bulge in the thick disc to be weaker and further away from the plane than in the thin disc. Regarding the kinematics, we find that the los velocity of thick disc stars in the outer parts of the b/p bulge can be \emph{larger} than that of thin disc stars, by up to 40\% and 20\% for side-on and Milky Way-like orientations of the bar respectively. This is due to the different orbits followed by thin and thick disc stars in the bar-b/p region, which are affected by the fact that: i) thin disc stars are trapped more efficiently in the bar - b/p instability and thus lose more angular momentum than their thick disc counterparts and ii) thick disc stars have large radial excursions and therefore stars from large radii with high angular momenta can be found in the bar region. We also find that the difference between the los velocities of the thin and thick disc in the b/p bulge ($Δv_{los}$) correlates with the initial difference between the radial velocity dispersions of the two discs ($Δσ$) . We therefore conclude that stars in the bar - b/p bulge will have considerably different morphologies and kinematics depending on the kinematic properties of the disc population they originate from. △ Less

Submitted 30 May, 2017; v1 submitted 3 April, 2017; originally announced April 2017.

Comments: Accepted for publication in A&A. 15 pages (2 page appendix). 16 figures

Journal ref: A&A 606, A47 (2017)

Analysing the 21cm signal from the Epoch of Reionization with artificial neural networks

Authors: Hayato Shimabukuro, Benoit Semelin

Abstract: The 21 cm signal from the Epoch of Reionization should be observed within the next decade. While a simple statistical detection is expected with SKA pathfinders, the SKA will hopefully produce a full 3D mapping of the signal. To extract from the observed data constraints on the parameters describing the underlying astrophysical processes, inversion methods must be developed. For example, the Marko… ▽ More The 21 cm signal from the Epoch of Reionization should be observed within the next decade. While a simple statistical detection is expected with SKA pathfinders, the SKA will hopefully produce a full 3D mapping of the signal. To extract from the observed data constraints on the parameters describing the underlying astrophysical processes, inversion methods must be developed. For example, the Markov Chain Monte Carlo method has been successfully applied. Here we test another possible inversion method: artificial neural networks (ANN). We produce a training set which consists of 70 individual sample. Each sample is made of the 21 cm power spectrum at different redshifts produced with the 21cmFast code plus the value of three parameters used in the semi-numerical simulations that describe astrophysical processes. Using this set we train the network to minimize the error between the parameter values it produces as an output and the true values. We explore the impact of the architecture of the network on the quality of the training. Then we test the trained network on the new set of 54 test samples with different values of the parameters. We find that the quality of the parameter reconstruction depends on the sensitivity of the power spectrum to the different parameters at a given redshift, that including thermal noise and sample variance decreases the quality of the reconstruction and that using the power spectrum at several redshifts as an input to the ANN improves the quality of the reconstruction. We conclude that ANNs are a viable inversion method whose main strength is that they require a sparse exploration of the parameter space and thus should be usable with full numerical simulations. △ Less

Submitted 22 March, 2017; v1 submitted 24 January, 2017; originally announced January 2017.

Comments: 8 figures, 2 tables, accepted in MNRAS

Journal ref: Mon Not R Astron Soc (2017) 468 (4): 3869-3877

On the kinematic detection of accreted streams in the Gaia era: a cautionary tale

Authors: I. Jean-Baptiste, P. Di Matteo, M. Haywood, A. Gomez, M. Montuori, F. Combes, B. Semelin

Abstract: The $Λ$CDM cosmological scenario predicts that our Galaxy should contain hundreds of stellar streams at the solar vicinity, fossil relics of the merging history of the Milky Way and more generally of the hierarchical growth of galaxies. Because of the mixing time scales in the inner Galaxy, it has been claimed that these streams should be difficult to detect in configuration space but can still be… ▽ More The $Λ$CDM cosmological scenario predicts that our Galaxy should contain hundreds of stellar streams at the solar vicinity, fossil relics of the merging history of the Milky Way and more generally of the hierarchical growth of galaxies. Because of the mixing time scales in the inner Galaxy, it has been claimed that these streams should be difficult to detect in configuration space but can still be identifiable in kinematic-related spaces like the energy/angular momenta spaces, E-Lz and Lperp-Lz, or spaces of orbital/velocity parameters. By means of high-resolution, dissipationless N-body simulations, containing between 25$\times10^6$ and 35$\times10^6$ particles, we model the accretion of a series of up to four 1:10 mass ratio satellites then up to eight 1:100 satellites and we search systematically for the signature of these accretions in these spaces. In all spaces considered (1) each satellite gives origin to several independent overdensities; (2) overdensities of multiple satellites overlap; (3) satellites of different masses can produce similar substructures; (4) the overlap between the in-situ and the accreted population is considerable everywhere; (5) in-situ stars also form substructures in response to the satellite(s) accretion. These points are valid even if the search is restricted to kinematically-selected halo stars only. As we are now entering the 'Gaia era', our results warn that an extreme caution must be employed before interpreting overdensities in any of those spaces as evidence of relics of accreted satellites. Reconstructing the accretion history of our Galaxy will require a substantial amount of accurate spectroscopic data, that, complemented by the kinematic information, will possibly allow us to (chemically) identify accreted streams and measure their orbital properties. (abridged) △ Less

Submitted 22 November, 2016; originally announced November 2016.

Comments: Accepted on A&A. A high-resolution version of the paper is available at http://aramis.obspm.fr/~paola/ELZ/Elz.pdf

Journal ref: A&A 604, A106 (2017)

Detailed modelling of the 21-cm Forest

Authors: Benoit Semelin

Abstract: The 21-cm forest is a promising probe of the Epoch of Reionization. The local state of the intergalactic medium (IGM) is encoded in the spectrum of a background source (radio-loud quasars or gamma ray burst afterglow) by absorption at the local 21-cm wavelength, resulting in a continuous and fluctuating absorption level. Small-scale structures (filaments and minihaloes) in the IGM are responsible… ▽ More The 21-cm forest is a promising probe of the Epoch of Reionization. The local state of the intergalactic medium (IGM) is encoded in the spectrum of a background source (radio-loud quasars or gamma ray burst afterglow) by absorption at the local 21-cm wavelength, resulting in a continuous and fluctuating absorption level. Small-scale structures (filaments and minihaloes) in the IGM are responsible for the strongest absorption features. The absorption can also be modulated on large scales by inhomogeneous heating and Wouthuysen-Field coupling. We present the results from a simulation that attempts to preserve the cosmological environment while resolving some of the small-scale structures (a few kpc resolution in a 50 Mpc/h box). The simulation couples the dynamics and the ionizing radiative transfer and includes X-ray and Lyman lines radiative transfer for a detailed physical modelling. As a result we find that soft X-ray self-shielding, Lyman-alpha self-shielding and shock heating all have an impact on the predicted values of the 21-cm optical depth of moderately overdense structures like filaments. An correct treatment of the peculiar velocities is also critical. Modelling these processes seems necessary for accurate predictions and can be done only at high enough resolution. As a result, based on our fiducial model, we estimate that LOFAR should be able to detect a few (strong) absorptions features in a frequency range of a few tens of MHz for a 20 mJy source located at z=10, while the SKA would extract a large fraction of the absorption information for the same source. △ Less

Submitted 8 October, 2015; originally announced October 2015.

Comments: 13 pages, accepted for publication in MNRAS

The Cosmic Dawn and Epoch of Reionization with the Square Kilometre Array

Authors: L. V. E. Koopmans, J. Pritchard, G. Mellema, F. Abdalla, J. Aguirre, K. Ahn, R. Barkana, I. van Bemmel, G. Bernardi, A. Bonaldi, F. Briggs, A. G. de Bruyn, T. C. Chang, E. Chapman, X. Chen, B. Ciardi, K. K. Datta, P. Dayal, A. Ferrara, A. Fialkov, F. Fiore, K. Ichiki, I. T. Illiev, S. Inoue, V. Jelić , et al. (20 additional authors not shown)

Abstract: Concerted effort is currently ongoing to open up the Epoch of Reionization (EoR) ($z\sim$15-6) for studies with IR and radio telescopes. Whereas IR detections have been made of sources (Lyman-$α$ emitters, quasars and drop-outs) in this redshift regime in relatively small fields of view, no direct detection of neutral hydrogen, via the redshifted 21-cm line, has yet been established. Such a direct… ▽ More Concerted effort is currently ongoing to open up the Epoch of Reionization (EoR) ($z\sim$15-6) for studies with IR and radio telescopes. Whereas IR detections have been made of sources (Lyman-$α$ emitters, quasars and drop-outs) in this redshift regime in relatively small fields of view, no direct detection of neutral hydrogen, via the redshifted 21-cm line, has yet been established. Such a direct detection is expected in the coming years, with ongoing surveys, and could open up the entire universe from $z\sim$6-200 for astrophysical and cosmological studies, opening not only the EoR, but also its preceding Cosmic Dawn ($z\sim$30-15) and possibly even the later phases of the Dark Ages ($z\sim$200-30). All currently ongoing experiments attempt statistical detections of the 21-cm signal during the EoR, with limited signal-to-noise. Direct imaging, except maybe on the largest (degree) scales at lower redshifts, as well as higher redshifts will remain out of reach. The Square Kilometre Array(SKA) will revolutionize the field, allowing direct imaging of neutral hydrogen from scales of arc-minutes to degrees over most of the redshift range $z\sim$6-28 with SKA1-LOW, and possibly even higher redshifts with the SKA2-LOW. In this SKA will be unique, and in parallel provide enormous potential of synergy with other upcoming facilities (e.g. JWST). In this chapter we summarize the physics of 21-cm emission, the different phases the universe is thought to go through, and the observables that the SKA can probe, referring where needed to detailed chapters in this volume (Abridged). △ Less

Submitted 28 May, 2015; originally announced May 2015.

Comments: Accepted for publication in the SKA Science Book 'Advancing Astrophysics with the Square Kilometre Array', to appear in 2015. PoS(AASKA14)001

The physics of Reionization: processes relevant for SKA observations

Authors: Benoit Semelin, Ilian Iliev

Abstract: The local intensity of the 21 cm signal emitted during the Epoch of Reionization that will be mapped by the SKA is modulated by the amount of neutral hydrogen. Consequently, understanding the process of reionization of the intergalactic medium (IGM) is crucial for predicting and interpreting the upcoming observations. After presenting the basic physics and most meaningful quantities pertaining to… ▽ More The local intensity of the 21 cm signal emitted during the Epoch of Reionization that will be mapped by the SKA is modulated by the amount of neutral hydrogen. Consequently, understanding the process of reionization of the intergalactic medium (IGM) is crucial for predicting and interpreting the upcoming observations. After presenting the basic physics and most meaningful quantities pertaining to the process of reionization, we will review recent progress in our understanding of the production and escape of ionizing photons in primordial galaxies and of their absorption in the IGM especially in so-called minihalos and Lyman Limit Systems. △ Less

Submitted 19 January, 2015; originally announced January 2015.

Comments: 10 pages, 1 figure. Accepted for publication in the SKA Science Book 'Advancing Astrophysics with the Square Kilometre Array', to appear in 2015

Why the Milky Way's bulge is not only a bar formed from a cold thin disk

Authors: P. Di Matteo, A. Gomez, M. Haywood, F. Combes, M. D. Lehnert, M. Ness, O. N. Snaith, D. Katz, B. Semelin

Abstract: By analyzing a N-body simulation of a bulge formed simply via a bar instability mechanism operating on a kinematically cold stellar disk, and by comparing the results of this analysis with the structural and kinematic properties of the main stellar populations of the Milky Way bulge, we conclude that the bulge of our Galaxy is not a pure stellar bar formed from a pre-existing thin stellar disk, as… ▽ More By analyzing a N-body simulation of a bulge formed simply via a bar instability mechanism operating on a kinematically cold stellar disk, and by comparing the results of this analysis with the structural and kinematic properties of the main stellar populations of the Milky Way bulge, we conclude that the bulge of our Galaxy is not a pure stellar bar formed from a pre-existing thin stellar disk, as some studies have recently suggested. On the basis of several arguments emphasized in this paper, we propose that the bulge population which, in the Milky Way, is observed not to be part of the peanut structure corresponds to the old galactic thick disk, thus implying that the Milky Way is a pure thin+thick disk galaxy, with only a possible limited contribution of a classical bulge. △ Less

Submitted 5 November, 2014; originally announced November 2014.

Comments: 9 pages, 4 figures, accepted for publication on A&A

Journal ref: A&A 577, A1 (2015)

Mapping a stellar disk into a boxy bulge: The outside-in part of the Milky Way bulge formation

Authors: P. Di Matteo, M. Haywood, A. Gomez, L. van Damme, F. Combes, A. Halle, B. Semelin, M. D. Lehnert, D. Katz

Abstract: By means of idealized, dissipationless N-body simulations which follow the formation and subsequent buckling of a stellar bar, we study the characteristics of boxy/peanut-shaped bulges and compare them with the properties of the stellar populations in the Milky Way bulge. The main results of our modeling, valid for the general family of boxy/peanut shaped bulges, are the following: (i) because of… ▽ More By means of idealized, dissipationless N-body simulations which follow the formation and subsequent buckling of a stellar bar, we study the characteristics of boxy/peanut-shaped bulges and compare them with the properties of the stellar populations in the Milky Way bulge. The main results of our modeling, valid for the general family of boxy/peanut shaped bulges, are the following: (i) because of the spatial redistribution in the disk initiated at the epoch of bar formation, stars from the innermost regions to the outer Lindblad resonance of the stellar bar are mapped into a boxy bulge; (ii) the contribution of stars to the local bulge density depends on their birth radius: stars born in the innermost disk tend to dominate the innermost regions of the boxy bulge, while stars originating closer to the OLR are preferentially found in the outer regions of the boxy/peanut structure; (iii) stellar birth radii are imprinted in the bulge kinematics, the larger the birth radii of stars ending up in the bulge, the greater their rotational support and the higher their line-of- sight velocity dispersions (but note that this last trend depends on the bar viewing angle); (iv) the higher the classical bulge-over-disk ratio, the larger its fractional contribution of stars at large vertical distance from the galaxy mid-plane. (ABRIDGED) On the basis of their chemical and kinematic characteristics, the results of our modeling suggests that the populations A, B and C, as defined by the ARGOS survey, can be associated, respectively, with the inner thin disk, to the young thick and to the old thick disk, following the nomenclature recently suggested for stars in the solar neighborhood by Haywood et al. (2013). △ Less

Submitted 16 April, 2014; v1 submitted 1 April, 2014; originally announced April 2014.

Comments: 21 pages, 14 figures, A&A accepted

Journal ref: A&A 567, A122 (2014)

Light cone anisotropy in 21 cm signal from the epoch of reionization

Authors: Karolina Zawada, Benoit Semelin, Patrick Vonlanthen, Sunghye Baek, Yves Revaz

Abstract: Using a suite of detailed numerical simulations we estimate the level of anisotropy generated by the time evolution along the light cone of the 21cm signal from the epoch of reionization. Our simulations include the physics necessary to model the signal during both the late emission regime and the early absorption regime, namely X-ray and Lyman-band 3D radiative transfer in addition to the usual d… ▽ More Using a suite of detailed numerical simulations we estimate the level of anisotropy generated by the time evolution along the light cone of the 21cm signal from the epoch of reionization. Our simulations include the physics necessary to model the signal during both the late emission regime and the early absorption regime, namely X-ray and Lyman-band 3D radiative transfer in addition to the usual dynamics and ionizing UV transfer. The signal is analysed using correlation functions perpendicular and parallel to the line of sight (LOS). We reproduce general findings from previous theoretical studies: the overall amplitude of the correlations and the fact that the light cone anisotropy is visible only on large scales (100 cMpc). However, the detailed behaviour is different. At 3 different epochs, the amplitude of the correlations along and perpendicular to the LOS differ from each other, indicating anisotropy. These 3 epochs are associated with 3 events of the global reionization history: the overlap of ionized bubbles, the onset of mild heating by X-rays in regions around the sources, and the onset of efficient Lyman-alpha coupling in regions around the sources. A 20x20 deg^2 survey area may be necessary to mitigate sample variance when we use the directional correlation functions. On a 100 cMpc scale the light cone anisotropy dominates over the anisotropy generated by peculiar velocity gradients computed in the linear regime. By modelling instrumental noise and limited resolution, we find that the anisotropy should be easily detectable by the SKA, assuming perfect foreground removal, the limiting factor being a large enough survey size. In the case of the LOFAR, it is likely that only first anisotropy episode will fall in the observing frequency range and will be detectable only if sample variance is much reduced (i.e. a larger than 20x20 deg^2 survey, which is not presently planned). △ Less

Submitted 8 January, 2014; originally announced January 2014.

Comments: 13 pages, 12 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

Signatures of radial migration in barred galaxies: Azimuthal variations in the metallicity distribution of old stars

Authors: P. Di Matteo, M. Haywood, F. Combes, B. Semelin, O. N. Snaith

Abstract: By means of N-body simulations, we show that radial migration in galaxy disks, induced by bar and spiral arms, leads to significant azimuthal variations in the metallicity distribution of old stars at a given distance from the galaxy center. Metals do not show an axisymmetric distribution during phases of strong migration. Azimuthal variations are visible during the whole phase of strong bar phase… ▽ More By means of N-body simulations, we show that radial migration in galaxy disks, induced by bar and spiral arms, leads to significant azimuthal variations in the metallicity distribution of old stars at a given distance from the galaxy center. Metals do not show an axisymmetric distribution during phases of strong migration. Azimuthal variations are visible during the whole phase of strong bar phase, and tend to disappear as the effect of radial migration diminishes, together with a reduction in the bar strength. These results suggest that the presence of inhomogeneities in the metallicity distribution of old stars in a galaxy disk can be a probe of ongoing strong migration. Such signatures may be detected in the Milky Way by Gaia (and complementary spectroscopic data), as well as in external galaxies, by IFU surveys like CALIFA and ATLAS3D. Mixing - defined as the tendency toward a homogeneous, azimuthally symmetric, stellar distribution in the disk - and migration turns out to be two distinct processes, the effects of mixing starting to be visible when strong migration is over. △ Less

Submitted 11 January, 2013; originally announced January 2013.

Comments: 8 pages, 10 figures, accepted for publication on Astronomy and Astrophysics

Formation and Radiative Feedback of First Objects and First Galaxies

Authors: Masayuki Umemura, Hajime Susa, Kenji Hasegawa, Tamon Suwa, Benoit Semelin

Abstract: First, the formation of first objects driven by dark matter is revisited by high-resolution hydrodynamic simulations. It is revealed that dark matter haloes of ~10^4M_sun can produce first luminous objects with the aid of dark matter cusps. Therefore, the mass of first objects is smaller by roughly two orders of magnitude than in the previous prediction. This implies that the number of Pop III sta… ▽ More First, the formation of first objects driven by dark matter is revisited by high-resolution hydrodynamic simulations. It is revealed that dark matter haloes of ~10^4M_sun can produce first luminous objects with the aid of dark matter cusps. Therefore, the mass of first objects is smaller by roughly two orders of magnitude than in the previous prediction. This implies that the number of Pop III stars formed in the early universe could be significantly larger than hitherto thought. Secondly, the feedback by photo-ionization and photo-dissociation photons in the first objects is explored with radiation hydrodynamic simulations, and it is demonstrated that multiple stars can form in a 10^5M_sun halo. Thirdly, the fragmentation of an accretion disk around a primordial protostar is explored with photo-dissociation feedback. As a result, it is found that the photo-dissociation can reduce the mass accretion rate onto protostars. Also, protostars as small as 0.8M_sun may be ejected and evolve with keeping their mass, which might be detected as "real first stars" in the Galactic halo. Finally, state-of-the-art radiation hydrodynamic simulations are performed to investigate the internal ionization of first galaxies and the escape of ionizing photons. We find that UV feedback by forming massive stars enhances the escape fraction even in a halo as massive as > 6* 10^9M_sun, while it reduces the star formation rate significantly. This may have a momentous impact on the cosmic reionization. △ Less

Submitted 29 October, 2012; originally announced October 2012.

Comments: 26 pages, 14 figures

Journal ref: Progress of Theoretical and Experimental Physics, 2012, 01A306

Reionization and the Cosmic Dawn with the Square Kilometre Array

Authors: Garrelt Mellema, León Koopmans, Filipe Abdalla, Gianni Bernardi, Benedetta Ciardi, Soobash Daiboo, Ger de Bruyn, Kanan K. Datta, Heino Falcke, Andrea Ferrara, Ilian T. Iliev, Fabio Iocco, Vibor Jelić, Hannes Jensen, Ronniy Joseph, Hans-Rainer Kloeckner, Panos Labroupoulos, Avery Meiksin, Andrei Mesinger, Andre Offringa, V. N. Pandey, Jonathan R. Pritchard, Mario G. Santos, Dominik J. Schwarz, Benoit Semelin , et al. (3 additional authors not shown)

Abstract: The Square Kilometre Array (SKA) will have a low frequency component (SKA-low) which has as one of its main science goals the study of the redshifted 21cm line from the earliest phases of star and galaxy formation in the Universe. This 21cm signal provides a new and unique window on both the formation of the first stars and accreting black holes and the later period of substantial ionization of th… ▽ More The Square Kilometre Array (SKA) will have a low frequency component (SKA-low) which has as one of its main science goals the study of the redshifted 21cm line from the earliest phases of star and galaxy formation in the Universe. This 21cm signal provides a new and unique window on both the formation of the first stars and accreting black holes and the later period of substantial ionization of the intergalactic medium. The signal will teach us fundamental new things about the earliest phases of structure formation, cosmology and even has the potential to lead to the discovery of new physical phenomena. Here we present a white paper with an overview of the science questions that SKA-low can address, how we plan to tackle these questions and what this implies for the basic design of the telescope. △ Less

Submitted 24 March, 2013; v1 submitted 30 September, 2012; originally announced October 2012.

Comments: Accepted for publication in Experimental Astronomy, reformatted to 57 pages, some updated and improved figures and minor changes and updates to the text

The impacts of ultraviolet radiation feedback on galaxies during the epoch of reionization

Authors: Kenji Hasegawa, Benoit Semelin

Abstract: We explore the impacts of ultraviolet (UV) radiation feedback on galaxies during the epoch of reionisation by cosmological simulations in which hydrodynamics and the transfer of the H and He ionising photons are consistently coupled. Moreover we take into account H_2 non-equilibrium chemistry, including photo-dissociation. The most striking feature of the simulations is a high spatial resolution f… ▽ More We explore the impacts of ultraviolet (UV) radiation feedback on galaxies during the epoch of reionisation by cosmological simulations in which hydrodynamics and the transfer of the H and He ionising photons are consistently coupled. Moreover we take into account H_2 non-equilibrium chemistry, including photo-dissociation. The most striking feature of the simulations is a high spatial resolution for the radiative transfer (RT) calculation which enables us to start considering not only external UV feedback processes but also internal UV feedback processes in each galaxy. We find that the star formation is significantly suppressed due to the internal UV and supernova (SN) feedback. In low mass galaxies with M<10^9Msun, a large amount of gas is evacuated by photo-evaporation as previous studies have shown, which results in the suppression of star formation. Surprisingly, star formation in massive halos is also strongly suppressed despite the fact that these halos hardly lose any gas by photo-evaporation. The suppression of star formation in massive halos is mainly caused by following two factors; (i) small scale clumpy structures in the galaxies are smoothened by the internal feedback, (ii) although the dense gas in the galaxies is mostly neutral, the H_2 formation and cooling processes are disturbed by mild photo-heating. Photo-dissociating radiations actually suppress star formation, but the magnitude of the effect is not so large in massive galaxies. Even though our simulation volume is too small to be a representative patch of the Universe during reionisation, we find that our simulated star formation rate densities and HI fractions at z=6-7 are consistent with those found in observations. △ Less

Submitted 6 March, 2013; v1 submitted 18 September, 2012; originally announced September 2012.

Comments: 14 pages, 13 figures, 1 table, accepted for publication in MNRAS

Joint Lyα emitters - quasars reionization constraints

Authors: S. Baek, A. Ferrara, B. Semelin

Abstract: We present a novel method to investigate cosmic reionization, using joint spectral information on high redshift Lyman Alpha Emitters (LAE) and quasars (QSOs). Although LAEs have been proposed as reionization probes, their use is hampered by the fact their Lyα line is damped not only by intergalactic HI but also internally by dust. Our method allows to overcome such degeneracy. First, we carefully… ▽ More We present a novel method to investigate cosmic reionization, using joint spectral information on high redshift Lyman Alpha Emitters (LAE) and quasars (QSOs). Although LAEs have been proposed as reionization probes, their use is hampered by the fact their Lyα line is damped not only by intergalactic HI but also internally by dust. Our method allows to overcome such degeneracy. First, we carefully calibrate a reionization simulation with QSO absorption line experiments. Then we identify LAEs in two simulation boxes at z=5.7 and z=6.6 and we build synthetic images/spectra of a prototypical LAE. At redshift 5.7, we find that the Lyα transmissivity (T_LAE) ~ 0.25, almost independent of the halo mass. This constancy arises from the conspiracy of two effects: (i) the intrinsic Lyα line width and (ii) the infall peculiar velocity. At higher redshift, z=6.6, where the transmissivity is instead largely set by the local HI abundance and LAE transmissivity consequently increases with halo mass from 0.15 to 0.3. Although outflows are present, they are efficiently pressure-confined by infall in a small region around the LAE; hence they only marginally affect transmissivity. Finally, we cast LOS originating from background QSOs passing through foreground LAEs at different impact parameters, and compute the quasar transmissivity (T_QSO). At smaller impact parameters, d < 1 cMpc, a positive correlation between T_QSO and halo mass is found at z = 5.7, which tends to become less pronounced (i.e. flatter) at larger distances. By cross-correlating T_LAE and T_QSO, we can obtain a HI density estimate unaffected by dust. At z= 5.7, the cross-correlation is relatively weak,whereas at z = 6.6 we find a clear positive correlation. We conclude by briefly discussing the perspectives for the application of the method to existing and forthcoming data. △ Less

Submitted 13 February, 2012; originally announced February 2012.

Comments: 14 pages, 15 figures

Mass assembly of galaxies: Smooth accretion versus mergers

Authors: Benjamin L'Huillier, Francoise Combes, Benoit Semelin

Abstract: Galaxies accrete their mass by means of both smooth accretion from the cosmic web, and the mergers of smaller entities. We wish to quantify the respective role of these two modes of accretion, which could determine the morphological types of galaxies observed today. Multi-zoom cosmological simulations are used to estimate as a function of time the evolution of mass in bound systems, for dark matte… ▽ More Galaxies accrete their mass by means of both smooth accretion from the cosmic web, and the mergers of smaller entities. We wish to quantify the respective role of these two modes of accretion, which could determine the morphological types of galaxies observed today. Multi-zoom cosmological simulations are used to estimate as a function of time the evolution of mass in bound systems, for dark matter as well as baryons. The baryonic contents of dark matter haloes are studied. Merger histories are followed as a function of external density, and the different ways in which mass is assembled in galaxies and the stellar component accumulated are quantified. We find that most galaxies assemble their mass through smooth accretion, and only the most massive galaxies also grow significantly through mergers. The mean fraction of mass assembled by accretion is 77 %, and by mergers 23 %. We present typical accretion histories of hundreds of galaxies: masses of the most massive galaxies increase monotonically in time, mainly through accretion, many intermediate-mass objects also experience mass-loss events such as tidal stripping and evaporation. However, our simulations suffer from the overcooling of massive galaxies caused by the neglect of active galaxy nuclei (AGN) feedback. The time by which half of the galay mass has assembled, both in dark matter and baryons, is a decreasing function of mass, which is compatible with the observations of a so-called downsizing. At every epoch in the universe, there are low-mass galaxies actively forming stars, while more massive galaxies form their stars over a shorter period of time within half the age of the universe. △ Less

Submitted 21 June, 2012; v1 submitted 22 August, 2011; originally announced August 2011.

Comments: A&A Accepted, 19 pages, 17 figures

Journal ref: A&A (2012) 544, A68

Distinctive rings in the 21 cm signal of the epoch of reionization

Authors: P. Vonlanthen, B. Semelin, S. Baek, Y. Revaz

Abstract: It is predicted that sources emitting UV radiation in the Lyman band during the epoch of reionization (EoR) showed a series of discontinuities in their Ly-alpha flux radial profile as a consequence of the thickness of the Lyman line series in the primeval intergalactic medium. Through unsaturated Wouthuysen-Field coupling, these spherical discontinuities are also present in the 21 cm emission of t… ▽ More It is predicted that sources emitting UV radiation in the Lyman band during the epoch of reionization (EoR) showed a series of discontinuities in their Ly-alpha flux radial profile as a consequence of the thickness of the Lyman line series in the primeval intergalactic medium. Through unsaturated Wouthuysen-Field coupling, these spherical discontinuities are also present in the 21 cm emission of the neutral IGM. In this article, we study the effects these discontinuities have on the differential brightness temperature of the 21 cm signal of neutral hydrogen in a realistic setting including all other sources of fluctuations. We focus on the early phases of the EoR, and we address the question of the detectability by the planned Square Kilometre Array. Such a detection would be of great interest, because these structures could provide an unambiguous diagnostic for the cosmological origin of the signal remaining after the foreground cleaning procedure. Also, they could be used as a new type of standard rulers. We determine the differential brightness temperature of the 21 cm signal in the presence of inhomogeneous Wouthuysen-Field effect using simulations which include (hydro)dynamics and both ionizing and Lyman lines 3D radiative transfer with the code LICORICE. We find that the Lyman horizons are clearly visible on the maps and radial profiles around the first sources of our simulations, but for a limited time interval, typically Δz \approx 2 at z \sim 13. Stacking the profiles of the different sources of the simulation at a given redshift results in extending this interval to Δz \approx 4. When we take into account the implementation and design planned for the SKA (collecting area, sensitivity, resolution), we find that detection will be challenging. It may be possible with a 10 km diameter for the core, but will be difficult with the currently favored design of a 5 km core. △ Less

Submitted 27 June, 2011; v1 submitted 2 March, 2011; originally announced March 2011.

Comments: 10 pages, 10 figures; v2: Section 5.5 rewritten; some new references added; accepted for publication in Astronomy and Astrophysics

Gas inflows, star formation and metallicity evolution in galaxy pairs

Authors: P. Di Matteo, M. Montuori, M. D. Lehnert, F. Combes, B. Semelin

Abstract: It has been known since many decades that galaxy interactions can induce star formation (hereafter SF) enhancements and that one of the driving mechanisms of this enhancement is related to gas inflows into the central galaxy regions, induced by asymmetries in the stellar component, like bars. In the last years many evidences have been accumulating, showing that interacting pairs have central gas-p… ▽ More It has been known since many decades that galaxy interactions can induce star formation (hereafter SF) enhancements and that one of the driving mechanisms of this enhancement is related to gas inflows into the central galaxy regions, induced by asymmetries in the stellar component, like bars. In the last years many evidences have been accumulating, showing that interacting pairs have central gas-phase metallicities lower than those of field galaxies, by {\sim} 0.2-0.3 dex on average. These diluted ISM metallicities have been explained as the result of inflows of metal-poor gas from the outer disk to the galaxy central regions. A number of questions arises: What's the timing and the duration of this dilution? How and when does the SF induced by the gas inflow enrich the circumnuclear gas with re-processed material? Is there any correlation between the timing and strength of the dilution and the timing and intensity of the SF? By means of Tree-SPH simulations of galaxy major interactions, we have studied the effect that gas inflows have on the ISM dilution, and the effect that the induced SF has, subsequently, in re-enriching the nuclear gas. In this contribution, we present the main results of this study. △ Less

Submitted 10 February, 2011; originally announced February 2011.

Comments: Proceedings of the IAU Symposium 277 "Tracing the Ancestry of Galaxies", 4 pages, 2 figures

Galaxy Formation: Merger vs Gas Accretion

Authors: Benjamin L'Huillier, Francoise Combes, Benoit Semelin

Abstract: According to the hierarchical model, small galaxies form first and merge together to form bigger objects. In parallel, galaxies assemble their mass through accretion from cosmic filaments. Recently, the increased spatial resolution of the cosmological simulations have emphasised that a large fraction of cold gas can be accreted by galaxies. In order to compare the role of both phenomena and the co… ▽ More According to the hierarchical model, small galaxies form first and merge together to form bigger objects. In parallel, galaxies assemble their mass through accretion from cosmic filaments. Recently, the increased spatial resolution of the cosmological simulations have emphasised that a large fraction of cold gas can be accreted by galaxies. In order to compare the role of both phenomena and the corresponding star formation history, one has to detect the structures in the numerical simulations and to follow them in time, by building a merger tree. △ Less

Submitted 28 December, 2010; originally announced December 2010.

Comments: 2 pages, 2 figures, Proceedings SF2A 2010

The GalMer database: Galaxy Mergers in the Virtual Observatory

Authors: Igor Chilingarian, Paola Di Matteo, Francoise Combes, Anne-Laure Melchior, Benoit Semelin

Abstract: We present the GalMer database, a library of galaxy merger simulations, made available to users through tools compatible with the Virtual Observatory (VO) standards adapted specially for this theoretical database. To investigate the physics of galaxy formation through hierarchical merging, it is necessary to simulate galaxy interactions varying a large number of parameters: morphological types, ma… ▽ More We present the GalMer database, a library of galaxy merger simulations, made available to users through tools compatible with the Virtual Observatory (VO) standards adapted specially for this theoretical database. To investigate the physics of galaxy formation through hierarchical merging, it is necessary to simulate galaxy interactions varying a large number of parameters: morphological types, mass ratios, orbital configurations, etc. On one side, these simulations have to be run in a cosmological context, able to provide a large number of galaxy pairs, with boundary conditions given by the large-scale simulations, on the other side the resolution has to be high enough at galaxy scales, to provide realistic physics. The GalMer database is a library of thousands simulations of galaxy mergers at moderate spatial resolution and it is a compromise between the diversity of initial conditions and the details of underlying physics. We provide all coordinates and data of simulated particles in FITS binary tables. The main advantages of the database are VO access interfaces and value-added services which allow users to compare the results of the simulations directly to observations: stellar population modelling, dust extinction, spectra, images, visualisation using dedicated VO tools. The GalMer value-added services can be used as virtual telescope producing broadband images, 1D spectra, 3D spectral datacubes, thus making our database oriented towards the usage by observers. We present several examples of the GalMer database scientific usage obtained from the analysis of simulations and modelling their stellar population properties, including: (1) studies of the star formation efficiency in interactions; (2) creation of old counter-rotating components; (3) reshaping metallicity profiles in elliptical galaxies; (4) orbital to internal angular momentum transfer; (5) reproducing observed colour bimodality of galaxies. △ Less

Submitted 16 March, 2010; originally announced March 2010.

Comments: 15 pages, 11 figures, 10 tables accepted to A&A. Visualisation of GalMer simulations, access to snapshot files and value-added tools described in the paper are available at http://galmer.obspm.fr/

The dilution peak, metallicity evolution, and dating of galaxy interactions and mergers

Authors: M. Montuori, P. Di Matteo, M. D. Lehnert, F. Combes, B. Semelin

Abstract: Strong inflows of gas from the outer disk to the inner kiloparsecs are induced during the interaction of disk galaxies. This inflow of relatively low-metallicity gas dilutes the metallicity of the circumnuclear gas. We have investigated several aspects of the process as the timing and duration of the dilution and its correlation with the induced star formation. We analysed major (1:1) gas-rich int… ▽ More Strong inflows of gas from the outer disk to the inner kiloparsecs are induced during the interaction of disk galaxies. This inflow of relatively low-metallicity gas dilutes the metallicity of the circumnuclear gas. We have investigated several aspects of the process as the timing and duration of the dilution and its correlation with the induced star formation. We analysed major (1:1) gas-rich interactions and mergers, spanning a range of initial orbital characteristics. Star formation and metal enrichment from SNe are included in our model. Our results show that the strongest trend is between the star formation rate and the dilution of the metals in the nuclear region; i.e., the more intense the central burst of star formation, the more the gas is diluted. This trend comes from strong inflows of relatively metal-poor gas from the outer regions of both disks, which fuels the intense star formation and lowers the overall metallicity for a time. The strong inflows happen on timescales of about 10^8 years or less, and the most intense star formation and lowest gas phase metallicities are seen generally after the first pericentre passage. As the star formation proceeds and the merger advances, the dilution reduces and enrichment becomes dominant - ultimately increasing the metallicity of the circumnuclear gas to a level higher than the initial metallicities of the merging galaxies. The "fly-bys" - pairs that interact but do not merge - also cause some dilution. We even see some dilution early in the merger or in the "fly-bys" and thus do not observe a strong trend between the nuclear metallicities and separation in our simulations until the merger is well advanced. We also analyse the O and Fe enrichment of the ISM, and show that the evolution of the alpha/Fe ratios, as well as the dilution of the central gas metallicity, can be used as a clock for "dating" the interaction. △ Less

Submitted 1 July, 2010; v1 submitted 6 March, 2010; originally announced March 2010.

Comments: 13 pages, 11 figures, A&A in press. Final version after language editing

Reionization by UV or X-ray sources

Authors: S. Baek, B. Semelin, P. Di Matteo, Y. Revaz, F. Combes

Abstract: We present simulations of the 21-cm signal during the epoch of reionization. We focus on properly modeling the absorption regime in the presence of inhomogeneous Wouthuysen-Field effect and X-ray heating. We ran radiative transfer simulations for three bands in the source spectrum (Lyman, UV, and X-ray) to fully account for these processes. We find that the brightness temperature fluctuation of th… ▽ More We present simulations of the 21-cm signal during the epoch of reionization. We focus on properly modeling the absorption regime in the presence of inhomogeneous Wouthuysen-Field effect and X-ray heating. We ran radiative transfer simulations for three bands in the source spectrum (Lyman, UV, and X-ray) to fully account for these processes. We find that the brightness temperature fluctuation of the 21 cm signal has an amplitude greater than 100 mK during the early reionization, up to 10 times greater than the typical amplitude of a few 10 mK obtained during the later emission phase. More importantly, we find that even a rather high contribution from QSO-like sources only damps the absorption regime without erasing it. Heating the IGM with X-ray takes time. Our results show that observations of the early reionization will probably benefit from a higher signal-to-noise value than during later stages. After analyzing the statistical properties of the signal (power spectrum and PDF) we find three diagnostics to constrain the level of X-ray, hence the nature of the first sources. △ Less

Submitted 12 August, 2010; v1 submitted 3 March, 2010; originally announced March 2010.

Comments: 15 pages, 7 figures, 2 tables, accepted for publication in A&A. High resolution version available at http://homepage.sns.it/baek/quasar.pdf

How ubiquitous are massive starbursts in interacting galaxies?

Authors: P. Di Matteo, F. Bournaud, M. Martig, F. Combes, A. -L. Melchior, B. Semelin

Abstract: Many evidences exist for a connection between galaxy interactions and induced star formation. However, a large range of responses of galaxies to tidal interactions is found, both in observations and in numerical simulations. We will discuss some recent results obtained analysing a large sample (~ 1000) of simulations of interacting pairs and their agreement with the most recent observational wor… ▽ More Many evidences exist for a connection between galaxy interactions and induced star formation. However, a large range of responses of galaxies to tidal interactions is found, both in observations and in numerical simulations. We will discuss some recent results obtained analysing a large sample (~ 1000) of simulations of interacting pairs and their agreement with the most recent observational works. △ Less

Submitted 28 September, 2009; originally announced September 2009.

Comments: 4 pages, Proceedings of the conference "Galaxies in isolation: exploring Nature vs Nurture" (May 2009, Granada, Spain)

Generation of rotationally dominated galaxies by mergers of pressure-supported progenitors

Authors: P. Di Matteo, C. J. Jog, M. D. Lehnert, F. Combes, B. Semelin

Abstract: Through the analysis of a set of numerical simulations of major mergers between initially non-rotating, pressure supported progenitor galaxies with a range of central mass concentrations, we have shown that: (1) it is possible to generate elliptical-like galaxies, with v/sigma > 1 outside one effective radius, as a result of the conversion of orbital- into internal-angular momentum; (2) the oute… ▽ More Through the analysis of a set of numerical simulations of major mergers between initially non-rotating, pressure supported progenitor galaxies with a range of central mass concentrations, we have shown that: (1) it is possible to generate elliptical-like galaxies, with v/sigma > 1 outside one effective radius, as a result of the conversion of orbital- into internal-angular momentum; (2) the outer regions acquire part of the angular momentum first; (3) both the baryonic and the dark matter components of the remnant galaxy acquire part of the angular momentum, the relative fractions depend on the initial concentration of the merging galaxies. For this conversion to occur the initial baryonic component must be sufficiently dense and/or the encounter should take place on a orbit with high angular momentum. Systems with these hybrid properties have been recently observed through a combination of stellar absorption lines and planetary nebulae for kinematic studies of early-type galaxies. Our results are in qualitative agreement with such observations and demonstrate that even mergers composed of non-rotating, pressure-supported progenitor galaxies can produce early-type galaxies with significant rotation at large radii. △ Less

Submitted 29 May, 2009; originally announced June 2009.

Comments: 5 pages, 6 figures, 2 tables. Accepted for publication in A&A Letters

Cosmological Radiative Transfer Comparison Project II: The Radiation-Hydrodynamic Tests

Authors: Ilian T. Iliev, Daniel Whalen, Garrelt Mellema, Kyungjin Ahn, Sunghye Baek, Nickolay Y. Gnedin, Andrey V. Kravtsov, Michael Norman, Milan Raicevic, Daniel R. Reynolds, Daisuke Sato, Paul R. Shapiro, Benoit Semelin, Joseph Smidt, Hajime Susa, Tom Theuns, Masayuki Umemura

Abstract: The development of radiation hydrodynamical methods that are able to follow gas dynamics and radiative transfer self-consistently is key to the solution of many problems in numerical astrophysics. Such fluid flows are highly complex, rarely allowing even for approximate analytical solutions against which numerical codes can be tested. An alternative validation procedure is to compare different m… ▽ More The development of radiation hydrodynamical methods that are able to follow gas dynamics and radiative transfer self-consistently is key to the solution of many problems in numerical astrophysics. Such fluid flows are highly complex, rarely allowing even for approximate analytical solutions against which numerical codes can be tested. An alternative validation procedure is to compare different methods against each other on common problems, in order to assess the robustness of the results and establish a range of validity for the methods. Previously, we presented such a comparison for a set of pure radiative transfer tests (i.e. for fixed, non-evolving density fields). This is the second paper of the Cosmological Radiative Transfer (RT) Comparison Project, in which we compare 9 independent RT codes directly coupled to gasdynamics on 3 relatively simple astrophysical hydrodynamics problems: (5) the expansion of an H II region in a uniform medium; (6) an ionization front (I-front) in a 1/r^2 density profile with a flat core, and (7), the photoevaporation of a uniform dense clump. Results show a broad agreement between the different methods and no big failures, indicating that the participating codes have reached a certain level of maturity and reliability. However, many details still do differ, and virtually every code has showed some shortcomings and has disagreed, in one respect or another, with the majority of the results. This underscores the fact that no method is universal and all require careful testing of the particular features which are most relevant to the specific problem at hand. △ Less

Submitted 18 May, 2009; originally announced May 2009.

Comments: 36 pages, 44 figures (most in color), submitted to MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.400:1283-1316,2009

On the survival of metallicity gradients to major dry-mergers

Authors: P. Di Matteo, A. Pipino, M. D. Lehnert, F. Combes, B. Semelin

Abstract: Using a large suite of galaxies with a variety of concentrations and metallicity gradients, we study the evolution of non-dissipative ("dry") equal mass mergers. Our purpose in generating these simulations is to explore how the metallicity gradient in dry mergers depends on the structure and metallicity gradients of the galaxies involved in the merger. Specifically, we would like to answer: Coul… ▽ More Using a large suite of galaxies with a variety of concentrations and metallicity gradients, we study the evolution of non-dissipative ("dry") equal mass mergers. Our purpose in generating these simulations is to explore how the metallicity gradient in dry mergers depends on the structure and metallicity gradients of the galaxies involved in the merger. Specifically, we would like to answer: Could dry mergers lead to metallicity gradients as observed in elliptical galaxies in the local Universe? Do dry mergers always lead to a flattening of the initial (i.e., pre-merger) gradient? From this modeling, we conclude that: The ratio of the remnant and the initial galaxy slopes span a wide range of values, up to values greater than 1 (with values greater than one resulting only when companions have gradients twice the progenitor). For a merger between two ellipticals having identical initial metallicity slopes (i.e., equal companion and galaxy slopes), the metallicity profile of the remnant flattens, with a final gradient about 0.6 times the initial one. Ellipticals can maintain their original pre-merger metallicity gradient if the companion slope is sufficiently steep. The amount of flattening neither depends on the characteristics of the orbit of the progenitors or on their initial concentration. Given the diversity in outcomes of the mergers, we conclude that dry mergers do not violate any observational constraints on the systematic characteristics of metallicity gradients in local ellipticals. In fact, dry mergers may be important within the context of the results of our simulations and may explain the large scatter in gradients for massive ellipticals and the relative paucity of massive ellipticals with no or shallow metallicity gradients. △ Less

Submitted 16 March, 2009; originally announced March 2009.

Comments: 11 pages, 11 figures, 4 tables, A&A accepted

On the frequency, intensity and duration of starburst episodes triggered by galaxy interactions and mergers

Authors: P. Di Matteo, F. Bournaud, M. Martig, F. Combes, A. -L. Melchior, B. Semelin

Abstract: We investigate the intensity enhancement and the duration of starburst episodes, triggered by major galaxy interactions and mergers. To this aim, we analyze two large statistical datasets of numerical simulations. These have been obtained using two independent and different numerical techniques to model baryonic and dark matter evolution, that are extensively compared for the first time. One is… ▽ More We investigate the intensity enhancement and the duration of starburst episodes, triggered by major galaxy interactions and mergers. To this aim, we analyze two large statistical datasets of numerical simulations. These have been obtained using two independent and different numerical techniques to model baryonic and dark matter evolution, that are extensively compared for the first time. One is a Tree-SPH code, the other one is a grid-based N-body sticky-particles code. We show that, at low redshift, galaxy interactions and mergers in general trigger only moderate star formation enhancements. Strong starbursts where the star formation rate is increased by a factor larger than 5 are rare and found only in about 15% of major galaxy interactions and mergers. Merger-driven starbursts are also rather short-lived, with a typical duration of the activity of a few 10^8 yr. These conclusions are found to be robust, independent from the numerical techniques and star formation models. At higher redshifts where galaxies contain more gas, gas inflow-induced starbursts are neither stronger neither longer than their local counterparts. In turn, the formation of massive gas clumps, results of local Jeans instability that can occur spontaneously in gas-rich disks or be indirectly favored by galaxy interactions, could play a more important role in determining the duration and intensity of star formation episodes. △ Less

Submitted 15 September, 2008; originally announced September 2008.

Comments: 22 pages, 28 figures, A&A accepted. High resolution version available at http://aramis.obspm.fr/~paola/SFR_frequency/

The simulated 21 cm signal during the epoch of reionization : full modeling of the Ly-alpha pumping

Authors: S. Baek, P. Di Matteo, B. Semelin, F. Combes, Y. Revaz

Abstract: The 21 cm emission of neutral hydrogen is the most promising probe of the epoch of reionization(EoR). In the next few years, the SKA pathfinders will provide statistical measurements of this signal. Numerical simulations predicting these observations are necessary to optimize the design of the instruments. The main difficulty is the computation of the spin temperature of neutral hydrogen which d… ▽ More The 21 cm emission of neutral hydrogen is the most promising probe of the epoch of reionization(EoR). In the next few years, the SKA pathfinders will provide statistical measurements of this signal. Numerical simulations predicting these observations are necessary to optimize the design of the instruments. The main difficulty is the computation of the spin temperature of neutral hydrogen which depends on the gas kinetic temperature and on the level of the local Lyman-alpha flux. A T_s >> T_cmb assumption is usual. However, this assumption does not apply early in the reionization history, or even later in the history as long as the sources of X-rays are too weak to heat the intergalactic medium significantly. This work presents the first EoR numerical simulations including, beside dynamics and ionizing continuum radiative transfer, a self-consistent treatment of the Ly-alpha radiative transfer. This allows us to compute the spin temperature more accurately. We use two different box sizes, 20 Mpc/h and 100 Mpc/h, and a star source model. Using the redshift dependence of average quantities, maps, and power spectra, we quantify the effect of using different assumptions to compute the spin temperature and the influence of the box size. The first effect comes from allowing for a signal in absorption. The magnitude of this effect depends on the amount of heating by hydrodynamic shocks and X-rays in the intergalactic medium(IGM). The second effects comes from using the real, local, Lyman-alpha flux. This effect is important for an average ionization fraction of less than 10%: it changes the overall amplitude of the 21 cm signal, and adds its own fluctuations to the power spectrum. △ Less

Submitted 1 December, 2008; v1 submitted 6 August, 2008; originally announced August 2008.

Comments: 20 pages, 16 figures, 2 tables, To be publish A&A. High resolution version available at http://aramis.obspm.fr/~baek/21cm_Lya.pdf

Old stellar counter-rotating components in elliptical-spiral mergers: exploring the GalMer database

Authors: P. Di Matteo, F. Combes, A. -L. Melchior, B. Semelin

Abstract: We investigate, by means of numerical simulations, the kinematics of elliptical-spiral merger remnants. Counterrotation can appear both in coplanar and in non-coplanar retrograde mergers, and it is mostly associated to the presence of a disk component, which preserves part of its initial spin. In turn, the external regions of the two interacting galaxies acquire part of the orbital angular momen… ▽ More We investigate, by means of numerical simulations, the kinematics of elliptical-spiral merger remnants. Counterrotation can appear both in coplanar and in non-coplanar retrograde mergers, and it is mostly associated to the presence of a disk component, which preserves part of its initial spin. In turn, the external regions of the two interacting galaxies acquire part of the orbital angular momentum, due to the action of tidal forces. △ Less

Submitted 23 November, 2007; originally announced November 2007.

Comments: 2 pages, 2 figures. To appear in the proceedings of "Formation and Evolution of Galaxy Disks" (J. G. Funes, and E. M. Corsini eds)

Old stellar counter-rotating components in early-type galaxies from elliptical-spiral mergers

Authors: P. Di Matteo, F. Combes, A. -L. Melchior, B. Semelin

Abstract: We investigate, by means of numerical simulations, the possibility of forming counter-rotating old stellar components by major mergers between an elliptical and a spiral galaxy. We show that counter-rotation can appear both in dissipative and dissipationless retrograde mergers, and it is mostly associated to the presence of a disk component, which preserves part of its initial spin. In turn, the… ▽ More We investigate, by means of numerical simulations, the possibility of forming counter-rotating old stellar components by major mergers between an elliptical and a spiral galaxy. We show that counter-rotation can appear both in dissipative and dissipationless retrograde mergers, and it is mostly associated to the presence of a disk component, which preserves part of its initial spin. In turn, the external regions of the two interacting galaxies acquire part of the orbital angular momentum, due to the action of tidal forces exerted on each galaxy by the companion. △ Less

Submitted 5 October, 2007; originally announced October 2007.

Comments: 6 pages, 15 figures. Accepted on Astronomy & Astrophysics

Chemodynamical evolution of interacting galaxies: the GalMer view

Authors: P. Di Matteo, F. Combes, I. Chilingarian, A. -L. Melchior, B. Semelin

Abstract: We have undertaken a large set of simulations of galaxy interactions and mergers (GalMer Project) in order to study the physical processes related to galaxy encounters. All morphological types along the Hubble sequence are considered in the initial conditions of the two colliding galaxies, with varying bulge-to-disk ratios and gas mass fractions. Different types of orbits are simulated, direct a… ▽ More We have undertaken a large set of simulations of galaxy interactions and mergers (GalMer Project) in order to study the physical processes related to galaxy encounters. All morphological types along the Hubble sequence are considered in the initial conditions of the two colliding galaxies, with varying bulge-to-disk ratios and gas mass fractions. Different types of orbits are simulated, direct and retrograde, according to the initial relative energy and impact parameter. The self gravity of stars, gas and dark matter is taken into account through a tree-code algorithm, the gas hydrodynamics through SPH. Star formation is included adopting a density-dependent Schmidt law. This wide library of galaxy interactions and mergers, containing, at present, about 900 simulations of major encounters, represents an unique tool to investigate statistically the chemodynamical evolution of interacting systems. In the following, we present and discuss some results obtained exploring the dataset, together with some future perspectives. △ Less

Submitted 17 September, 2007; originally announced September 2007.

Comments: 4 pages, 5 figures. To appear in the proceedings of "SF2A-2007: Semaine de l'Astrophysique Francaise", (J. Bouvier, A. Chalabaev, and C. Charbonnel eds)

Lyman-alpha radiative transfer during the Epoch of Reionization: contribution to 21-cm signal fluctuations

Authors: B. Semelin, F. Combes, S. Baek

Abstract: During the epoch of reionization, Ly-alpha photons emitted by the first stars can couple the neutral hydrogen spin temperature to the kinetic gas temperature, providing the opportunity to observe the gas in emission or absorption in the 21-cm line. Given the bright foregrounds, it is of prime importance to determine precisely the fluctuations signature of the signal, to be able to extract it by… ▽ More During the epoch of reionization, Ly-alpha photons emitted by the first stars can couple the neutral hydrogen spin temperature to the kinetic gas temperature, providing the opportunity to observe the gas in emission or absorption in the 21-cm line. Given the bright foregrounds, it is of prime importance to determine precisely the fluctuations signature of the signal, to be able to extract it by its correlation power. LICORICE is a Monte-Carlo radiative transfer code, coupled to the dynamics via an adaptative Tree-SPH code. We present here the Ly-alpha part of the implementation, and validate it through three classical tests. Contrary to previous works, we do not assume that P_alpha, the number of scatterings of Ly-alpha photons per atom per second, is proportional to the Ly-alpha background flux, but take into account the scatterings in the Ly-alpha line wings. The latter have the effect to steepen the radial profile of P_alpha around each source, and re-inforce the contrast of the fluctuations. In the particular geometry of cosmic filaments of baryonic matter, Ly-alpha photons are scattered out of the filament, and the large scale structure of P_alpha is significantly anisotropic. This could have strong implications for the possible detection of the 21-cm signal. △ Less

Submitted 17 July, 2007; originally announced July 2007.

Comments: 13 pages, 9 figures. To be published in A&A

Star formation efficiency in galaxy interactions and mergers: a statistical study

Authors: P. Di Matteo, F. Combes, A. -L. Melchior, B. Semelin

Abstract: We investigate the enhancement of star formation efficiency in galaxy interactions and mergers, by numerical simulations of several hundred galaxy collisions. All morphological types along the Hubble sequence are considered in the initial conditions of the two colliding galaxies, with varying bulge-to-disk ratios and gas mass fractions. Different types of orbits are simulated, direct and retrogr… ▽ More We investigate the enhancement of star formation efficiency in galaxy interactions and mergers, by numerical simulations of several hundred galaxy collisions. All morphological types along the Hubble sequence are considered in the initial conditions of the two colliding galaxies, with varying bulge-to-disk ratios and gas mass fractions. Different types of orbits are simulated, direct and retrograde, according to the initial relative energy and impact parameter, and the resulting star formation history is compared to that occuring in the two galaxies when they are isolated. Our principal results are: (1) retrograde encounters have a larger star formation efficiency (SFE) than direct encounters; (2) the amount of gas available in the galaxy is not the main parameter governing the SFE in the burst phase; (3) there is an anticorrelation between the amplitude of the star forming burst and the tidal forces exerted per unit of time, which is due to the large amount of gas dragged outside the galaxy by tidal tails in strong interactions; (4) globally, the Kennicutt-Schmidt law is retrieved statistically for isolated galaxies, interacting pairs and mergers; (5) the enhanced star formation is essentially occurring in nuclear starbursts, triggered by inward gas flows driven by non-axisymmetries in the galaxy disks. Direct encounters develop more pronounced asymmetries than retrograde ones. Based on these statistical results, we derive general laws for the enhancement of star formation in galaxy interactions and mergers, as a function of the main parameters of the encounter. △ Less

Submitted 9 March, 2007; originally announced March 2007.

Comments: 22 pages, 37 figures, 4 tables. Accepted on Astronomy & Astrophysics

The life-time of galactic bars: central mass concentrations and gravity torques

Authors: F. Bournaud, F. Combes, B. Semelin

Abstract: Bars in gas-rich spiral galaxies are short-lived. They drive gas inflows through their gravity torques, and at the same time self-regulate their strength. Their robustness has been subject of debate, since it was thought that only the resulting central mass concentrations (CMCs) were weakening bars, and only relatively rare massive CMCs were able to completely destroy them. Through numerical sim… ▽ More Bars in gas-rich spiral galaxies are short-lived. They drive gas inflows through their gravity torques, and at the same time self-regulate their strength. Their robustness has been subject of debate, since it was thought that only the resulting central mass concentrations (CMCs) were weakening bars, and only relatively rare massive CMCs were able to completely destroy them. Through numerical simulations including gas dynamics, we find that with the gas parameters of normal spiral galaxies, the CMC is not sufficient to fully dissolve the bar. But another overlooked mechanism, the transfer of angular momentum from the infalling gas to the stellar bar, can also strongly weaken the bar. In addition, we show that gravity torques are correctly reproduced in simulations, and conclude that bars are transient features, with life-time of 1-2 Gyr in typical Sb-Sc galaxies, because of the combined effects of CMCs and gravity torques, while most existing works had focussed on the CMC effects alone. △ Less

Submitted 6 September, 2005; originally announced September 2005.

Comments: accepted for publication in MNRAS. The definitive version is available at www.blackwell-synergy.com

Journal ref: Mon.Not.Roy.Astron.Soc.Lett.364:L18-L22,2005