BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Benedetta Ciardi

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

ResearcherUnclaimed source record
astro-ph.COastro-ph.GAastro-phastro-ph.HEastro-ph.IMastro-ph.EPphysics.geo-phphysics.space-ph

Catalog footprint

What is connected

39works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

39 published item(s)

preprint2022arXiv

A CRASH simulation of the contribution of binary stars to the epoch of reionization

We use a set of 3D radiative transfer simulations to study the effect that a large fraction of binary stars in galaxies during the epoch of reionization has on the physical properties of the intergalactic medium (i.e. the gas temperature and the ionization state of hydrogen and helium), on the topology of the ionized bubbles and on the 21 cm power spectra. Consistently to previous literature, we find that the inclusion of binary stars can speed up the reionization process of HI and HeI, while HeII reionization is still dominated by more energetic sources, especially accreting black holes. The earlier ionization attained with binary stars allows for more time for cooling and recombination, so that gas fully ionized by binary stars is typically colder than that ionized by single stars at any given redshift. With the same volume averaged ionization fraction, the inclusion of binary stars results in fewer small ionized bubbles and more large ones, with visible effects also on the large scales of the 21 cm power spectrum.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Atomic and molecular gas from the epoch of reionization down to redshift 2

Cosmic gas makes up about 90% of baryonic matter in the Universe and H$_2$ is the closest molecule to star formation. In this work we study cold neutral gas and its H$_2$ component at different epochs, exploiting state-of-the-art hydrodynamic simulations that include time-dependent atomic and molecular non-equilibrium chemistry coupled to star formation, feedback effects, different UV backgrounds presented in the recent literature and a number of additional processes - such as gas self-shielding, H$_2$ dust grain catalysis, photoelectric and cosmic-ray heating - occurring during structure formation (ColdSIM). We find neutral-gas mass density parameters $ Ω_{\rm neutral} \simeq $10$^{-3}$ and increasing from lower to higher redshift, in agreement with available HI data. Resulting H$_2$ fractions can be as high as $\sim $50% at $z\sim $4-8, in line with the latest high-$z$ measurements. Albeit dependent on the adopted UV background, derived $ Ω_{\rm H_2} $ values agree with observations up to $z\sim$7 and both HI and H$_2$ trends are better reproduced by our non-equilibrium H$_2$-based star formation modelling. The predicted gas depletion timescales decrease towards lower $z$, with H$_2$ depletion times remaining below the Hubble time and comparable to the dynamical time at all considered redshifts. This implies that non-equilibrium molecular cooling is efficient at driving cold-gas collapse in a broad variety of environments and since the very early cosmic epochs. In appendix, we show detailed analyses of individual processes, as well as simple numerical parameterizations and fits to account for them. Our findings suggest that, in addition to HI, non-equilibrium H$_2$ observations are pivotal probes for assessing cold-gas abundances and the role of UV background radiation - Abridged

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A LOFAR Observation of Ionospheric Scintillation from Two Simultaneous Travelling Ionospheric Disturbances

This paper presents the results from one of the first observations of ionospheric scintillation taken using the Low-Frequency Array (LOFAR). The observation was of the strong natural radio source Cas A, taken overnight on 18-19 August 2013, and exhibited moderately strong scattering effects in dynamic spectra of intensity received across an observing bandwidth of 10-80MHz. Delay-Doppler spectra (the 2-D FFT of the dynamic spectrum) from the first hour of observation showed two discrete parabolic arcs, one with a steep curvature and the other shallow, which can be used to provide estimates of the distance to, and velocity of, the scattering plasma. A cross-correlation analysis of data received by the dense array of stations in the LOFAR "core" reveals two different velocities in the scintillation pattern: a primary velocity of ~30m/s with a north-west to south-east direction, associated with the steep parabolic arc and a scattering altitude in the F-region or higher, and a secondary velocity of ~110m/s with a north-east to south-west direction, associated with the shallow arc and a scattering altitude in the D-region. Geomagnetic activity was low in the mid-latitudes at the time, but a weak sub-storm at high latitudes reached its peak at the start of the observation. An analysis of Global Navigation Satellite Systems (GNSS) and ionosonde data from the time reveals a larger-scale travelling ionospheric disturbance (TID), possibly the result of the high-latitude activity, travelling in the north-west to south-east direction, and, simultaneously, a smaller--scale TID travelling in a north-east to south-west direction, which could be associated with atmospheric gravity wave activity. The LOFAR observation shows scattering from both TIDs, at different altitudes and propagating in different directions. To the best of our knowledge this is the first time that such a phenomenon has been reported.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Bubble mapping with the Square Kilometer Array -- I. Detecting galaxies with Euclid, JWST, WFIRST and ELT within ionized bubbles in the intergalactic medium at z>6

The Square Kilometer Array is expected to provide the first tomographic observations of the neutral intergalactic medium at redshifts z>6 and pinpoint the locations of individual ionized bubbles during early stages of cosmic reionization. In scenarios where star-forming galaxies provide most of the ionizing photons required for cosmic reionization, one expects the first ionized bubbles to be centered on overdensities of such galaxies. Here, we model the properties of galaxy populations within isolated, ionized bubbles that SKA-1 should be able to resolve at z=7-10, and explore the prospects for galaxy counts within such structures with various upcoming near-infrared telescopes. We find that, for the bubbles that are within reach of SKA-1 tomography, the bubble volume is closely tied to the number of ionizing photons that have escaped from the galaxies within. In the case of galaxy-dominated reionization, galaxies are expected to turn up above the spectroscopic detection threshold of JWST and ELT in even the smallest resolvable bubbles at redshifts z=10 or below. The prospects of detecting galaxies within these structures in purely photometric surveys with Euclid, WFIRST, JWST or ELT are also discussed. While spectroscopy is preferable towards the end of reionization to provide a robust sample of bubble members, multiband imaging may be a competitive option for bubbles at z~10, due to the very small number of line-of-sight interlopers expected at similar redshifts.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Extended X-ray emission from the z=4.26 radio galaxy 4C 63.20

We report on deep Chandra X-ray Telescope imaging observations of 4C 63.20, one of the few known radio galaxies at z>3.5. The X-ray counterpart is resolved into a core plus two off-nuclear sources that (combined) account for close to 30% of the total X-ray flux. Their morphology and orientation are consistent with a diffuse, lobe-like nature, albeit compact hotspots cannot be ruled out. The broadband spectral energy distribution of 4C 63.20 can be reproduced with a jet model where the majority of the radio flux can be ascribed to synchrotron emission from the hotspots, whereas the (non-nuclear) X-ray emission is produced via Inverse Compton (IC) off of Cosmic Microwave Background (CMB) photons within the extended lobes. This scenario is broadly consistent with the expectation from highly magnetized lobes in a hotter CMB, and supports the view that IC/CMB may quench less extreme radio lobes at high redshifts.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Large scale simulations of H and He reionization and heating driven by stars and more energetic sources

We present simulations of cosmic reionization and reheating from $z=18$ to $z=5$, investigating the role of stars (emitting soft UV-photons), nuclear black holes (BHs, with power-law spectra), X-ray binaries (XRBs, with hard X-ray dominated spectra), and the supernova-associated thermal bremsstrahlung of the diffuse interstellar medium (ISM, with soft X-ray spectra). We post-process the hydrodynamical simulation Massive-Black II (MBII) with multifrequency ionizing radiative transfer. The source properties are directly derived from the physical environment of MBII, and our only real free parameter is the ionizing escape fraction $f_{\rm esc}$. We find that, among the models explored here, the one with an escape fraction that decreases with decreasing redshift yields results most in line with observations, such as of the neutral hydrogen fraction and the Thomson scattering optical depth. Stars are the main driver of hydrogen reionization and consequently of the thermal history of the intergalactic medium (IGM). We obtain $\langle x_{\rm HII} \rangle = 0.99998$ at $z=6$ for all source types, with volume averaged temperatures $\langle T \rangle \sim 20,000~{\rm K}$. BHs are rare and negligible to hydrogen reionization, but conversely they are the only sources which can fully ionize helium, increasing local temperatures by $\sim 10^4~{\rm K}$. The thermal and ionization state of the neutral and lowly ionized hydrogen differs significantly with different source combinations, with ISM and (to a lesser extent) XRBs, playing a significant role and, as a consequence, determining the transition from absorption to emission of the 21 cm signal from neutral hydrogen.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Predictions for the 21cm-galaxy cross-power spectrum observable with SKA and future galaxy surveys

In this paper we use radiative transfer + N-body simulations to explore the feasibility of measurements of cross-correlations between the 21cm field observed by the Square Kilometer Array (SKA) and high-z Lyman Alpha Emitters (LAEs) detected in galaxy surveys with the Subaru Hyper Supreme Cam (HSC), Subaru Prime Focus Spectrograph (PFS) and Wide Field Infrared Survey Telescope (WFIRST). 21cm-LAE cross-correlations are in fact a powerful probe of the epoch of reionization as they are expected to provide precious information on the progress of reionization and the typical scale of ionized regions at different redshifts. The next generation observations with SKA will have a noise level much lower than those with its precursor radio facilities, introducing a significant improvement in the measurement of the cross-correlations. We find that an SKA-HSC/PFS observation will allow to investigate scales below ~10 Mpc/h and ~60 Mpc/h at z=7.3 and 6.6, respectively. WFIRST will allow to access also higher redshifts, as it is expected to observe spectroscopically ~900 LAEs per square degree and unit redshift in the range 7.5<z<8.5. Because of the reduction of the shot noise compared to HSC and PFS, observations with WFIRST will result in more precise cross-correlations and increased observable scales.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Probing the high-z IGM with the hyperfine transition of $^3$He$^+$

The hyperfine transition of $^3$He$^+$ at 3.5cm has been thought as a probe of the high-z IGM since it offers a unique insight into the evolution of the helium component of the gas, as well as potentially give an independent constraint on the 21cm signal from neutral hydrogen. In this paper, we use radiative transfer simulations of reionization driven by sources such as stars, X-ray binaries, accreting black holes and shock heated interstellar medium, and simulations of a high-z quasar to characterize the signal and analyze its prospects of detection. We find that the peak of the signal lies in the range 1-50 $μ$K for both environments, but while around the quasar it is always in emission, in the case of cosmic reionization a brief period of absorption is expected. As the evolution of HeII is determined by stars, we find that it is not possible to distinguish reionization histories driven by more energetic sources. On the other hand, while a bright QSO produces a signal in 21cm that is very similar to the one from a large collection of galaxies, its signature in 3.5cm is very peculiar and could be a powerful probe to identify the presence of the QSO. We analyze the prospects of the signal's detectability using SKA1-mid as our reference telescope. We find that the noise power spectrum dominates over the power spectrum of the signal, although a modest S/N ratio can be obtained when the wavenumber bin width and the survey volume are sufficiently large.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Observing the redshifted 21 cm signal around a bright QSO at $z\sim 10$

We use hydrodynamics and radiative transfer simulations to study the 21~cm signal around a bright QSO at $z \sim 10$. Due to its powerful UV and X-ray radiation, the QSO quickly increases the extent of the fully ionized bubble produced by the pre-existing stellar type sources, in addition to partially ionize and heat the surrounding gas. As expected, a longer QSO lifetime, $t_{\rm QSO}$, results in a 21~cm signal in emission located at increasingly larger angular radii, $θ$, and covering a wider range of $θ$. Similar features can be obtained with a higher galactic emissivity efficiency, $f_{\rm UV}$, so that determining the origin of a large ionized bubble (i.e. QSO vs stars) is not straightforward. Such degeneracy could be reduced by taking advantage of the finite light traveltime effect, which is expected to affect an HII region produced by a QSO differently from one created by stellar type sources. From an observational point of view, we find that the 21 cm signal around a QSO at various $t_{\rm QSO}$ could be detected by SKA1-low with a high signal-noise ratio (S/N). As a reference, for $t_{\rm QSO} = 10\,\rm Myr$, a S/N $\sim 8$ is expected assuming that no pre-heating of the IGM has taken place due to high-$z$ energetic sources, while it can reach value above 10 in case of pre-heating. Observations of the 21~cm signal from the environment of a high-$z$ bright QSO could then be used to set constraints on its lifetime, as well as to reduce the degeneracy between $f_{\rm UV}$ and $t_{\rm QSO}$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Interpreting the evolution of galaxy colours from $z = 8$ to $z = 5$

We attempt to interpret existing data on the evolution of the UV luminosity function and UV colours, $β$, of galaxies at $5 \leq z \leq 8$, to improve our understanding of their dust content and ISM properties. To this aim, we post-process the results of a cosmological hydrodynamical simulation with a chemical evolution model, which includes dust formation by supernovae and intermediate mass stars, dust destruction in supernova shocks, and grain growth by accretion of gas-phase elements in dense gas. We find that observations require a steep, Small Magellanic Cloud-like extinction curve and a clumpy dust distribution, where stellar populations younger than 15 Myr are still embedded in their dusty natal clouds. Investigating the scatter in the colour distribution and stellar mass, we find that the observed trends can be explained by the presence of two populations: younger, less massive galaxies where dust enrichment is mainly due to stellar sources, and massive, more chemically evolved ones, where efficient grain growth provides the dominant contribution to the total dust mass. Computing the IR-excess - UV color relation we find that all but the dustiest model galaxies follow a relation shallower than the Meurer et al. (1999) one, usually adopted to correct the observed UV luminosities of high-$z$ galaxies for the effects of dust extinction. As a result, their total star formation rates might have been over-estimated. Our study illustrates the importance to incorporate a proper treatment of dust in simulations of high-$z$ galaxies, and that massive, dusty, UV-faint galaxies might have already appeared at $z \lesssim 7$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Systematic biases in low frequency radio interferometric data due to calibration: the LOFAR EoR case

The redshifted 21 cm line of neutral hydrogen is a promising probe of the Epoch of Reionization (EoR). However, its detection requires a thorough understanding and control of the systematic errors. We study two systematic biases observed in the LOFAR EoR residual data after calibration and subtraction of bright discrete foreground sources. The first effect is a suppression in the diffuse foregrounds, which could potentially mean a suppression of the 21 cm signal. The second effect is an excess of noise beyond the thermal noise. The excess noise shows fluctuations on small frequency scales, and hence it can not be easily removed by foreground removal or avoidance methods. Our analysis suggests that sidelobes of residual sources due to the chromatic point spread function and ionospheric scintillation can not be the dominant causes of the excess noise. Rather, both the suppression of diffuse foregrounds and the excess noise can occur due to calibration with an incomplete sky model containing predominantly bright discrete sources. We show that calibrating only on bright sources can cause suppression of other signals and introduce an excess noise in the data. The levels of the suppression and excess noise depend on the relative flux of sources which are not included in the model with respect to the flux of modeled sources. We discuss possible solutions such as using only long baselines to calibrate the interferometric gain solutions as well as simultaneous multi-frequency calibration along with their benefits and shortcomings.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

21cm Forest with the SKA

An alternative to both the tomography technique and the power spectrum approach is to search for the 21cm forest, that is the 21cm absorption features against high-z radio loud sources caused by the intervening cold neutral intergalactic medium (IGM) and collapsed structures. Although the existence of high-z radio loud sources has not been confirmed yet, SKA-low would be the instrument of choice to find such sources as they are expected to have spectra steeper than their lower-z counterparts. Since the strongest absorption features arise from small scale structures (few tens of physical kpc, or even lower), the 21cm forest can probe the HI density power spectrum on small scales not amenable to measurements by any other means. Also, it can be a unique probe of the heating process and the thermal history of the early universe, as the signal is strongly dependent on the IGM temperature. Here we show what SKA1-low could do in terms of detecting the 21cm forest in the redshift range z = 7.5-15.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Blue Wings and Bumps via Fermi-like Acceleration of Lyman-alpha photons across Shocks

We explore the impact of Fermi-like acceleration of Lyman-alpha (Lyα) photons across shock fronts on the observed Lyα spectral line shape. We first confirm the result of Neufeld & McKee (1988) that this mechanism gives rise to extended blue wings which may have been observed in some radio galaxies. Our Monte-Carlo radiative transfer calculations further show that in a minor modification of the shell-model, in which we add an additional static shell of hydrogen, this process can naturally explain the small blue bumps observed in a subset of Lyα emitting galaxies, which have been difficult to explain with conventional shell-models. Blue bumps can be produced with an additional column density of static hydrogen as small as $N_{HI}^{static} \ll N_{HI}^{shell}$, and typically occur at roughly the outflow velocity of the shell. In our model the spectra of so-called 'blue-bump objects' might reflect an evolutionary stage in which the outflows regulating the escape of Lyα photons are still engulfed within a static interstellar medium.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Bulk Flows and End of the Dark Ages with the SKA

The early Universe is a precious probe of the birth of primordial objects, first star formation events and consequent production of photons and heavy elements. Higher-order corrections to the cosmological linear perturbation theory predicts the formation of coherent supersonic gaseous streaming motions at decoupling time. These bulk flows impact the gas cooling process and determine a cascade effect on the whole baryon evolution. By analytical estimates and N-body hydrodynamical chemistry numerical simulations including atomic and molecular evolution, gas cooling, star formation, feedback effects and metal spreading for individual species from different stellar populations according to the proper yields and lifetimes, we discuss the role of these primordial bulk flows at the end of the dark ages and their detectable impacts during the first Gyr in view of the upcoming SKA mission. Early bulk flows can inhibit molecular gas cooling capabilities, suppressing star formation, metal spreading and the abundance of small primordial galaxies in the infant Universe. This can determine a delay in the re-ionization process and in the heating of neutral hydrogen making the observable HI signal during cosmic evolution patchier and noisier. The planned SKA mission will represent a major advance over existing instruments, since it will be able to probe the effects on HI 21cm at z ~ 6-20 and on molecular line emissions from first collapsing sites at z ~ 20-40. Therefore, it will be optimal to address the effects of primordial streaming motions on early baryon evolution and to give constraints on structure formation in the first Gyr.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Effects of the sources of reionization on 21-cm redshift-space distortions

The observed 21-cm signal from the epoch of reionization will be distorted along the line-of-sight by the peculiar velocities of matter particles. These redshift-space distortions will affect the contrast in the signal and will also make it anisotropic. This anisotropy contains information about the cross-correlation between the matter density field and the neutral hydrogen field, and could thus potentially be used to extract information about the sources of reionization. In this paper, we study a collection of simulated reionization scenarios assuming different models for the sources of reionization. We show that the 21-cm anisotropy is best measured by the quadrupole moment of the power spectrum. We find that, unless the properties of the reionization sources are extreme in some way, the quadrupole moment evolves very predictably as a function of global neutral fraction. This predictability implies that redshift-space distortions are not a very sensitive tool for distinguishing between reionization sources. However, the quadrupole moment can be used as a model-independent probe for constraining the reionization history. We show that such measurements can be done to some extent by first-generation instruments such as LOFAR, while the SKA should be able to measure the reionization history using the quadrupole moment of the power spectrum to great accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Lyman-Alpha Emitting Galaxies as a Probe of Reionization: Large-Scale Bubble Morphology and Small-Scale Absorbers

The visibility of LyA emitting galaxies during the Epoch of Reionization is controlled by both diffuse HI patches in large-scale bubble morphology and small-scale absorbers. To investigate the impact on LyA photons, we apply a novel combination of analytic and numerical calculations to three scenarios: (i) the `bubble' model, where only diffuse HI outside ionized bubbles is present; (ii) the `web' model, where HI exists only in overdense self-shielded gas; and (iii) the more realistic 'web-bubble' model, which contains both. Our analysis confirms that there is a degeneracy between the ionization structure of the intergalactic medium (IGM) and the HI fraction inferred from LyA surveys, as the three models suppress LyA flux equally with very different HI fractions. We argue that a joint analysis of the LyA luminosity function and the rest-frame equivalent width distribution/LyA fraction can break this degeneracy and provide constraints on the reionization history and its topology. We further show that constraints can improve if we consider the full shape of the M_UV-dependent redshift evolution of the LyA fraction of Lyman break galaxies. Contrary to conventional wisdom, we find that (i) a drop of LyA fraction larger for UV-faint than for UV-bright galaxies can be reproduced with web and web-bubble models and therefore does not provide exclusive evidence of patchy reionization, and (ii) the IGM-transmission PDF is unimodal for bubble models and bimodal in web models. We further highlight the importance of galaxy-absorber cross-correlation. Comparing our models to observations, the neutral fraction at z~7 is likely to be of order of tens of per cent when interpreted with bubble or web-bubble models. Alternatively, we obtain a conservative lower limit ~1% in the web models, if we allow for a drop in the photoionization rate by a factor of ~100 from the post-reionized universe. [abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Predictions for the 21cm-galaxy cross-power spectrum observable with LOFAR and Subaru

The 21cm-galaxy cross-power spectrum is expected to be one of the promising probes of the Epoch of Reionization (EoR), as it could offer information about the progress of reionization and the typical scale of ionized regions at different redshifts. With upcoming observations of 21cm emission from the EoR with the Low Frequency Array (LOFAR), and of high redshift Lyalpha emitters (LAEs) with Subaru's Hyper Suprime Cam (HSC), we investigate the observability of such cross-power spectrum with these two instruments, which are both planning to observe the ELAIS-N1 field at z=6.6. In this paper we use N-body + radiative transfer (both for continuum and Lyalpha photons) simulations at redshift 6.68, 7.06 and 7.3 to compute the 3D theoretical 21cm-galaxy cross-power spectrum, as well as to predict the 2D 21cm-galaxy cross-power spectrum expected to be observed by LOFAR and HSC. Once noise and projection effects are accounted for, our predictions of the 21cm-galaxy cross-power spectrum show clear anti-correlation on scales larger than ~ 60 h$^{-1}$ Mpc (corresponding to k ~ 0.1 h Mpc$^{-1}$), with levels of significance p=0.04 at z=6.6 and p=0.048 at z=7.3. On smaller scales, instead, the signal is completely contaminated.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

SKA - EoR correlations and cross-correlations: kSZ, radio galaxies, and NIR background

The Universe's Cosmic Dawn (CD) and Epoch of Reionization (EoR) can be studied using a number of observational probes that provide complementary or corroborating information. Each of these probes suffers from its own systematic and statistical uncertainties. It is therefore useful to consider the mutual information that these data sets contain. In this paper, we discuss a potential of cross-correlations between the SKA cosmological 21 cm data with: (i) the kinetic Sunyaev- Zel'dovich (kSZ) effect in the CMB data; (ii) the galaxy surveys; and (iii) near infrared (NIR) backgrounds.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Clumping factors of HII, HeII and HeIII

Estimating the intergalactic medium ionization level of a region needs proper treatment of the reionization process for a large representative volume of the universe. The clumping factor, a parameter which accounts for the effect of recombinations in unresolved, small-scale structures, aids in achieving the required accuracy for the reionization history even in simulations with low spatial resolution. In this paper, we study for the first time the redshift evolution of clumping factors of different ionized species of H and He in a small but very high resolution simulation of the reionization process. We investigate the dependence of the value and redshift evolution of clumping factors on their definition, the ionization level of the gas, the grid resolution, box size and mean dimensionless density of the simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Constraining the epoch of reionization with the variance statistic: simulations of the LOFAR case

Several experiments are underway to detect the cosmic redshifted 21-cm signal from neutral hydrogen from the Epoch of Reionization (EoR). Due to their very low signal-to-noise ratio, these observations aim for a statistical detection of the signal by measuring its power spectrum. We investigate the extraction of the variance of the signal as a first step towards detecting and constraining the global history of the EoR. Signal variance is the integral of the signal's power spectrum, and it is expected to be measured with a high significance. We demonstrate this through results from a simulation and parameter estimation pipeline developed for the Low Frequency Array (LOFAR)-EoR experiment. We show that LOFAR should be able to detect the EoR in 600 hours of integration using the variance statistic. Additionally, the redshift ($z_r$) and duration ($Δz$) of reionization can be constrained assuming a parametrization. We use an EoR simulation of $z_r = 7.68$ and $Δz = 0.43$ to test the pipeline. We are able to detect the simulated signal with a significance of 4 standard deviations and extract the EoR parameters as $z_r = 7.72^{+0.37}_{-0.18}$ and $Δz = 0.53^{+0.12}_{-0.23}$ in 600 hours, assuming that systematic errors can be adequately controlled. We further show that the significance of detection and constraints on EoR parameters can be improved by measuring the cross-variance of the signal by cross-correlating consecutive redshift bins.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Reionization and the Cosmic Dawn with the Square Kilometre Array

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Simulating extremely metal-poor gas and DLA metal content at redshift z=7

We present the first theoretical study of metals in damped Ly-alpha (DLA) systems at redshift z=7. The features of cold, primordial gas are studied by means of N-body, hydro, chemistry simulations, including atomic and molecular non-equilibrium chemistry, cooling, star formation for population III and population II-I regimes, stellar evolution, cosmic metal spreading according to proper yields (for He, C, O, Si, Fe, Mg, S, etc.) and lifetimes, and feedback effects. Theoretical expectations are then compared to recently available constraints from DLA observations. We find that DLA galaxies at z=7 account for 10 per cent of the whole galaxy population and for most of the metal-poor galaxies at these epochs. About 7 per cent of these DLA galaxies contain purely pristine material and 34 per cent of them consist of very weakly polluted gas, being, therefore, suitable candidate s as population III sites. The remaining 59 per cent are enriched above ~10^{-4} Zsun. Additionally, DLA candidates appear to have: gas masses <~2x10^8 Msun; very low star formation rate, ~10^{-3} - 10^{-2} Msun/yr (significantly weaker than late-time counterparts); mean molecular fractions covering a fairly wide range, x_mol~10^{-3}- 10^{-6}; typical metallicities Z<~3x10^{-3}Zsun and HI column densities NHI>~3x10^{20} cm^{-2}. They present no or weak correlations between their gas mass and Z, NHI, or x_mol; a moderate correlation between x_mol and Z, linked to the ongoing molecular-driven star formation and metal pollution processes; a mild anti-correlation between NHI and x_mol, due to H depletion into molecules; and a chemical content that is subject to environmental dependencies.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Simulating the assembly of galaxies at redshifts z = 6 - 12

We use state-of-the-art simulations to explore the physical evolution of galaxies in the first billion years of cosmic time. First, we demonstrate that our model reproduces the basic statistical properties of the observed Lyman-break galaxy (LBG) population at z = 6 - 8, including the evolving ultra-violet (UV) luminosity function (LF), the stellar-mass density (SMD), and the average specific star-formation rates (sSFR) of LBGs with M_{UV} < -18 (AB mag). Encouraged by this success we present predictions for the behaviour of fainter LBGs extending down to M_{UV} <= -15 (as will be probed with the James Webb Space Telescope) and have interrogated our simulations to try to gain insight into the physical drivers of the observed population evolution. We find that mass growth due to star formation in the mass-dominant progenitor builds up about 90% of the total z ~ 6 LBG stellar mass, dominating over the mass contributed by merging throughout this era. Our simulation suggests that the apparent "luminosity evolution" depends on the luminosity range probed: the steady brightening of the bright end of the LF is driven primarily by genuine physical luminosity evolution and arises due to a fairly steady increase in the UV luminosity (and hence star-formation rates) in the most massive LBGs. However, at fainter luminosities the situation is more complex, due in part to the more stochastic star-formation histories of lower-mass objects; at this end, the evolution of the UV LF involves a mix of positive and negative luminosity evolution (as low-mass galaxies temporarily brighten then fade) coupled with both positive and negative density evolution (as new low-mass galaxies form, and other low-mass galaxies are consumed by merging). We also predict the average sSFR of LBGs should rise from sSFR = 4.5 Gyr^-1 at z = 6 to about 11 Gyr^-1 by z = 9.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Counts of high-redshift GRBs as probe of primordial non-Gaussianities

We propose to use high-redshift long $γ$-ray bursts (GRBs) as cosmological tools to constrain the amount of primordial non-Gaussianity in the density field. By using numerical, N-body, hydrodynamic, chemistry simulations of different cosmological volumes with various Gaussian and non-Gaussian models, we self-consistently relate the cosmic star formation rate density to the corresponding GRB rate. Assuming that GRBs are fair tracers of cosmic star formation, we find that positive local non-Gaussianities, described in terms of the non-linear parameter, \fnl, might boost significantly the GRB rate at high redshift, $z \gg 6$. Deviations with respect to the Gaussian case account for a few orders of magnitude if \fnl$\sim 1000$, one order of magnitude for \fnl$\sim 100$, and a factor of $\sim 2$ for \fnl$\sim 50$. These differences are found only at large redshift, while at later times the rates tend to converge. Furthermore, a comparison between our predictions and the observed GRB data at $z > 6$ allows to exclude large negative \fnl, consistently with previous works. Future detections of any long GRB at extremely high redshift ($z\sim 15-20$) could favor non-Gaussian scenarios with positive \fnl. More stringent constraints require much larger high-$z$ GRB complete samples, currently not available in literature. By distinguishing the contributions to the GRB rate from the metal-poor population III regime, and the metal-enriched population II-I regime, we conclude that the latter is a more solid tracer of the underlying matter distribution, while the former is strongly dominated by feedback mechanisms from the first, massive, short-lived stars, rather than by possible non-Gaussian fluctuations. This holds quite independently of the assumed population III initial mass function.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dark Matter Halo Environment for Primordial Star Formation

We study the statistical properties (such as shape and spin) of high-z halos likely hosting the first (PopIII) stars with cosmological simulations including detailed gas physics. In the redshift range considered ($11 < z < 16$) the average sphericity is $<s> = 0.3 \pm 0.1$, and for more than 90% of halos the triaxiality parameter is $T \lesssim 0.4$, showing a clear preference for oblateness over prolateness. Larger halos in the simulation tend to be both more spherical and prolate: we find $s \propto M_h^{α_s}$ and $T \propto M_h^{α_T}$, with $α_s \approx 0.128$ and $α_T= 0.276$ at z = 11. The spin distributions of dark matter and gas are considerably different at $z=16$, with the baryons rotating slower than the dark matter. At lower redshift, instead, the spin distributions of dark matter and gas track each other almost perfectly, as a consequence of a longer time interval available for momentum redistribution between the two components. The spin of both the gas and dark matter follows a lognormal distribution, with a mean value at z=16 of $<λ> =0.0184$, virtually independent of halo mass. This is in good agreement with previous studies. Using the results of two feedback models (MT1 and MT2) by McKee & Tan (2008) and mapping our halo spin distribution into a PopIII IMF, we find that at high-$z$ the IMF closely tracks the spin lognormal distribution. Depending on the feedback model, though, the distribution can be centered at $\approx 65 M_\odot$ (MT1) or $\approx 140 M_\odot$ (MT2). At later times, model MT1 evolves into a bimodal distribution with a second prominent peak located at $35-40 M_\odot$ as a result of the non-linear relation between rotation and halo mass. We conclude that the dark matter halo properties might be a key factor shaping the IMF of the first stars.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Effect of Intergalactic Medium on the Observability of Lyman Alpha Emitters during Cosmic Reionization

We perform a systematic study of how the inhomogeneities in the Inter-Galactic Medium (IGM) affect the observability of Lyman Alpha Emitters (LAEs) around the Epoch of Reionization. We focus on the IGM close to the galaxies as the detailed ionization distribution and velocity fields of this region could significantly influence the scattering of Ly-alpha photons off neutral H atoms as they traverse the IGM after escaping from the galaxy. We simulate the surface brightness (SB) maps and spectra of more than 100 LAEs at z=7.7 as seen by an observer at z=0. To achieve this, we extract the source properties of galaxies and their surrounding IGM from cosmological simulations of box sizes 5-30 Mpc/h and follow the coupled radiative transfer of ionizing and Ly-alpha radiation through the IGM using CRASH-alpha. We find that the simulated SB profiles are extended and their detailed structure is affected by inhomogeneities in the IGM, especially at high neutral fractions. The detectability of LAEs and the fraction of the flux observed depend heavily on the shape of the SB profile and the SB threshold (SB_th) of the observational campaign. Only ultradeep observations (e.g. SB_th ~ 10^-23 ergs/s/cm^2/arcsec^2) would be able to obtain the true underlying mass-luminosity relation and luminosity functions of LAEs. The details of our results depend on whether Ly-alpha photons are significantly shifted in the galaxy to longer wavelengths, the mean ionization fraction in the IGM and the clustering of ionizing sources. These effects can lead to an easier escape of Ly-alpha photons with less scattering in the IGM and a concentrated SB profile similar to the one of a point source. Finally, we show that the SB profiles are steeper at high ionization fraction for the same LAE sample which can potentially be observed from the stacked profile of a large number of LAEs.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Searching for the first stars with the Gaia mission

We construct a theoretical model to predict the number of orphan afterglows (OA) from gamma-ray bursts (GRBs) triggered by primordial metal-free (Pop III) stars expected to be observed by the Gaia mission. In particular, we consider primordial metal-free stars that were affected by radiation from other stars (Pop III.2) as a possible target. We use a semi-analytical approach that includes all relevant feedback effects to construct cosmic star formation history and its connection with the cumulative number of GRBs. The OA events are generated using the Monte Carlo method, and realistic simulations of Gaia's scanning law are performed to derive the observation probability expectation. We show that Gaia can observe up to 2.28 $\pm$ 0.88 off-axis afterglows and 2.78 $\pm$ 1.41 on-axis during the five-year nominal mission. This implies that a nonnegligible percentage of afterglows that may be observed by Gaia ($\sim 10%$) could have Pop III stars as progenitors.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Enabling parallel computing in CRASH

We present the new parallel version (pCRASH2) of the cosmological radiative transfer code CRASH2 for distributed memory supercomputing facilities. The code is based on a static domain decomposition strategy inspired by geometric dilution of photons in the optical thin case that ensures a favourable performance speed-up with increasing number of computational cores. Linear speed-up is ensured as long as the number of radiation sources is equal to the number of computational cores or larger. The propagation of rays is segmented and rays are only propagated through one sub-domain per time step to guarantee an optimal balance between communication and computation. We have extensively checked pCRASH2 with a standardised set of test cases to validate the parallelisation scheme. The parallel version of CRASH2 can easily handle the propagation of radiation from a large number of sources and is ready for the extension of the ionisation network to species other than hydrogen and helium.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The interplay between chemical and mechanical feedback from the first generation of stars

We study cosmological simulations of early structure formation, including non-equilibrium molecular chemistry, metal pollution from stellar evolution, transition from population III (popIII) to population II (popII) star formation, regulated by a given critical metallicity, and feedback effects. We investigate the properties of early metal spreading from the different stellar populations and its interplay with primordial molecular gas. We find that, independently of the details about popIII modeling, after the onset of star formation, regions enriched below the critical level are mostly found in isolated environments, while popII star formation regions are much more clumped. Typical star forming haloes show average SN driven outflow rates of up to 10^{-4} Msun/yr in enriched gas, initially leaving the original star formation regions almost devoid of metals. The polluted material, which is gravitationally incorporated in over-dense environments on timescales of 10^7 yr, is mostly coming from external, nearby star forming sites ("gravitational enrichment"). In parallel, the pristine-gas inflow rates are between 10^{-3} - 10^{-1} Msun/yr. However, thermal feedback from SN generates turbulence and destroys molecules within the pristine gas, and only the polluted material, incorporated via gravitational enrichment, can continue to cool by atomic metal fine-structure transitions on time scales short enough to end the initial popIII regime within less than 10^8 yr.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Power spectrum extraction for redshifted 21-cm epoch of reionization experiments: the LOFAR case

One of the aims of the Low Frequency Array (LOFAR) Epoch of Reionization (EoR) project is to measure the power spectrum of variations in the intensity of redshifted 21-cm radiation from the EoR. The sensitivity with which this power spectrum can be estimated depends on the level of thermal noise and sample variance, and also on the systematic errors arising from the extraction process, in particular from the subtraction of foreground contamination. We model the extraction process using realistic simulations of the cosmological signal, the foregrounds and noise, and so estimate the sensitivity of the LOFAR EoR experiment to the redshifted 21-cm power spectrum. Detection of emission from the EoR should be possible within 360 hours of observation with a single station beam. Integrating for longer, and synthesizing multiple station beams within the primary (tile) beam, then enables us to extract progressively more accurate estimates of the power at a greater range of scales and redshifts. We discuss different observational strategies which compromise between depth of observation, sky coverage and frequency coverage. A plan in which lower frequencies receive a larger fraction of the time appears to be promising. We also study the nature of the bias which foreground fitting errors induce on the inferred power spectrum, and discuss how to reduce and correct for this bias. The angular and line-of-sight power spectra have different merits in this respect, and we suggest considering them separately in the analysis of LOFAR data.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Probing intergalactic radiation fields during cosmic reionization through gamma-ray absorption

We discuss expectations for the absorption of high-energy gamma-rays by gamma-gamma pair production with intergalactic radiation fields (IRFs) at very high redshifts (z~5-20), and the prospects thereof for probing the cosmic reionization era. For the evolving IRF, a semi-analytical model incorporating both Population II and Population III stars is employed, which is consistent with a wide variety of existing high-z observations including QSO spectral measurements, WMAP Thomson depth constraints, near-IR source count limits, etc. We find that the UV IRF below the Lyman edge energy with intensities in the range of a few times 10^{-19} erg cm^{-2} s^{-1} Hz^{-1} sr^{-1} can cause appreciable attenuation above ~12 GeV at z~5, down to ~6-8 GeV at z>~8-10. This may be observable in the spectra of blazars or gamma-ray bursts by the Fermi Gamma-ray Space Telescope or next generation facilities such as the Cherenkov Telescope Array, Advanced Gamma-ray Imaging System or 5@5, providing invaluable insight into early star formation and cosmic reionization.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The impact of primordial supersonic flows on early structure formation, reionization and the lowest-mass dwarf galaxies

Tseliakhovich & Hirata recently discovered that higher-order corrections to the cosmological linear-perturbation theory lead to supersonic coherent baryonic flows just after recombination (i.e. z~1020), with rms velocities of ~30 km/s relative to the underlying dark-matter distribution, on comoving scales of ~<3 Mpc/h. To study the impact of these coherent flows we performed high-resolution N-body plus SPH simulations in boxes of 5.0 and 0.7 Mpc/h, for bulk-flow velocities of 0 (as reference), 30 and 60 km/s. The simulations follow the evolution of cosmic structures by taking into account detailed, primordial, non-equilibrium gas chemistry (i.e. H, He, H2, HD, HeH, etc.), cooling, star formation, and feedback effects from stellar evolution. We find that these bulk flows suppress star formation in low-mass haloes (i.e. Mvir<~10^8Msun until z~13, lower the abundance of the first objects by ~1%-20%, and, as consequence, delay cosmic star formation history by ~2e7yr. The gas fractions in individual objects can change up to a factor of two at very early times. Coherent bulk flow, therefore, has implications for (i) the star-formation in the lowest-mass haloes (e.g. dSphs), (ii) the start of reionization by suppressing it in some patches of the Universe, and (iii) the heating (i.e. spin temperature) of neutral hydrogen. We speculate that the patchy nature of reionization and heating on several Mpc scales could lead to enhanced differences in the HI spin-temperature, giving rise to stronger variations in the HI brightness temperatures during the late dark ages.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Lyalpha versus X-ray heating in the high-z IGM

In this paper we examine the effect of X-ray and Lyalpha photons on the intergalactic medium temperature. We calculate the photon production from a population of stars and micro-quasars in a set of cosmological hydrodynamic simulations which self-consistently follow the dark matter dynamics, radiative processes as well as star formation, black hole growth and associated feedback processes. We find that, (i) IGM heating is always dominated by X-rays unless the Lyalpha photon contribution from stars in objects with mass M<10^8 Msun becomes significantly enhanced with respect to the X-ray contribution from BHs in the same halo (which we do not directly model). (ii) Without overproducing the unresolved X-ray background, the gas temperature becomes larger than the CMB temperature, and thus an associated 21 cm signal should be expected in emission, at z<11.5. We discuss how in such a scenario the transition redshift between a 21 cm signal in absorption and in emission could be used to constraint BHs accretion and associated feedback processes.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The contribution of the IGM and minihalos to the 21 cm signal of reionization

We study the statistical properties of the cosmological 21 cm signal from both the intergalactic medium (IGM) and minihalos, using a reionization simulation that includes a self--consistent treatment of minihalo photoevaporation. We consider two models for minihalo formation and three typical thermal states of the IGM -- heating purely by ionization, heating from both ionizing and $\rm Lyα$ photons, and a maximal "strong heating" model. We find that the signal from the IGM is almost always dominant over that from minihalos. In our calculation, the differential brightness temperature, $δT_b,$ of minihalos is never larger than 2 mK. Although there are indeed some differences in the signals from the minihalos and from the IGM, even with the planned generation of radio telescopes it will be unfeasible to detect them. However, minihalos significantly affect the ionization state of the IGM and the corresponding 21 cm flux.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The onset of star formation in primordial haloes

Star formation remains an unsolved problem in astrophysics. Numerical studies of large-scale structure simulations cannot resolve the whole process and their approach usually assumes that only gas denser than a typical threshold can host and form stars. We investigate the onset of cosmological star formation and compare several very-high-resolution, three-dimensional, N-body/SPH simulations that include non-equilibrium, atomic and molecular chemistry, star formation prescriptions, and feedback effects. We study how primordial star formation depends on gas density thresholds, cosmological parameters and initial set-ups. For mean-density initial conditions, we find that standard low-density star-formation threshold (0.2 h^2/cm3) models predict the onset of star formation at z~25-31, depending on the adopted cosmology. In these models stars are formed regardless of the time between the moment when the threshold is reached and the effective runaway collapse. At high redshift, this time interval represents a significant fraction of the Hubble time and thus this assumption can induce large artificial off-sets to the onset of star formation. Choosing higher density thresholds (135 h^2/cm3) allows the entire cooling process to be followed, and the onset of star formation is then estimated to be at redshift z~12-16. When isolated, rare, high-density peaks are considered, the chemical evolution is much faster and the first star formation episodes occur at z > 40, almost regardless of the choice for the density threshold. These results could have implications for the formation redshift of the first cosmological objects, as inferred from direct numerical simulations of mean-density environments, and on the studies of the reionization history of the universe.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

How rapidly do supermassive black hole "seeds" grow at early times?

We investigate the physical conditions for the growth of intermediate mass seed black holes assumed to have formed from remnants of the first generation of massive stars. We follow the collapse of high-sigma halos with Tvir > 1e4 K using cosmological, smooth-particle hydrodynamic (SPH) simulations in the standard LCDM model. During collapse of the parent halo the seed holes are incorporated through mergers into larger systems and accrete mass from the surrounding gas. We include a self-consistent treatment of star formation, black hole accretion and associated feedback processes. Even under optimistic assumptions for the seed black hole mass and for efficient merger rates, we find that seed holes in halos M<1e10 Msun never reach the conditions for critical Eddington growth. Most of the black hole growth in this regime is determined by the initial mass and the merger rates. Critical accretion rates are reached, albeit only after a significant delay, at the time of collapse z~7) for 3-4 sigma halos of M~1e11 Msun. Our results imply M_BH = 5e6 Msun (M_halo/1e11 Msun)^0.78 at the time of collapse. The required conditions of Eddington growth to explain the build-up of supermassive black holes (~1e9 Msun), as implied by Sloan quasars at z>6, are therefore hard to meet in such a scenario. Without a 'jump-start' these conditions may be only achieved in extremely rare halos with M_halo > 1e13 Msun that collapsed before z~6. The sub-Eddington regime in which black holes holes accrete at early time implies a small contribution to the reionization by miniquasar but still sufficient to cause appreciable heating of the IGM at z<15-18.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Detecting the First Objects in the Mid-IR with NGST

We calculate the expected Mid-IR molecular hydrogen line emission from the first objects in the universe. Due to their low masses, the stellar feedback from massive stars is able to blowaway their gas content and collect it into a cooling shell where H_2 rapidly forms and IR roto-vibrational (as for example the restframe 2.12 micron) lines carry away a large fraction (up to 10%) of the explosion energy. The fluxes from these sources are in the range 10^{-21}-10^{-17} erg/s/cm^2. The highest number counts are expected in the 20 micron band, where about 10^5 sources/deg^2 are predicted at the limiting flux of 3 \times 10^{-18} erg/s/cm^2. Among the planned observational facilities, we find that the best detection perspectives are offered by NGST, which should be able to reveal about 200 first objects in one hour observation time at its limiting flux in the above band. Therefore, Mid-IR instruments appear to represent perfect tools to trace star formation and stellar feedback in the high (z>5) redshift universe.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Ly-alpha Clouds Associated with Pop III Objects in Cold Dark Matter Models

We present a semi-analytical model to test the hypothesis that Ly-alpha QSO absorption lines originate in gaseous halos produced by multiple supernova explosions occurring in Pop III objects in a CDM cosmological scenario hoping to assess the validity of CDM models and/or constrain their parameters. The preliminary results indicate that the range of $N_{HI}$, redshift distribution and metallicity of clouds are well reproduced by CDM if they are associated with galaxy halos or groups. Firm conclusions on clouds with $N_{HI}\le 10^{14} cm^{-2}$ need a more refined study.

0 citations0 reviewsNo code linkedOpen access