Source author record

Jaime E. Forero-Romero

Jaime E. Forero-Romero appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

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

23 published item(s)

preprint2025arXiv

Galaxy-Multiplet Clustering from DESI DR2

We present an efficient estimator for higher-order galaxy clustering using small groups of nearby galaxies, or multiplets. Using the Luminous Red Galaxy (LRG) sample from the Dark Energy Spectroscopic Instrument (DESI) Data Release 2, we identify galaxy multiplets as discrete objects and measure their cross-correlations with the general galaxy field. Our results show that the multiplets exhibit stronger clustering bias as they trace more massive dark matter halos than individual galaxies. When comparing the observed clustering statistics with the mock catalogs generated from the N-body simulation AbacusSummit, we find that the mocks underpredict multiplet clustering despite reproducing the galaxy two-point auto-correlation reasonably well. This discrepancy indicates that the standard Halo Occupation Distribution (HOD) model is insufficient to describe the properties of galaxy multiplets, revealing the greater constraining power of this higher-order statistic on galaxy-halo connection and the possibility that multiplets are specific to additional assembly bias. We demonstrate that incorporating secondary biases into the HOD model improves agreement with the observed multiplet statistics, specifically by allowing galaxies to preferentially occupy halos in denser environments. Our results highlight the potential of utilizing multiplet clustering, beyond traditional two-point correlation measurements, to break degeneracies in models describing the galaxy-dark matter connection.

preprint2023arXiv

The Target-selection Pipeline for the Dark Energy Spectroscopic Instrument

In 2021 May, the Dark Energy Spectroscopic Instrument (DESI) began a 5 yr survey of approximately 50 million total extragalactic and Galactic targets. The primary DESI dark-time targets are emission line galaxies (ELGs), luminous red galaxies (LRGs) and quasars (QSOs). In bright time, DESI will focus on two surveys known as the Bright Galaxy Survey (BGS) and the Milky Way Survey (MWS). DESI also observes a selection of "secondary" targets for bespoke science goals. This paper gives an overview of the publicly available pipeline (desitarget) used to process targets for DESI observations. Highlights include details of the different DESI survey targeting phases, the targeting ID (TARGETID) used to define unique targets, the bitmasks used to indicate a particular type of target, the data model and structure of DESI targeting files, and examples of how to access and use the desitarget code base. This paper will also describe "supporting" DESI target classes, such as standard stars, sky locations, and random catalogs that mimic the angular selection function of DESI targets. The DESI target selection pipeline is complex and sizable; this paper attempts to summarize the most salient information required to understand and work with DESI targeting data.

preprint2022arXiv

Snowmass2021 Cosmic Frontier White Paper: Prospects for obtaining Dark Matter Constraints with DESI

Despite efforts over several decades, direct-detection experiments have not yet led to the discovery of the dark matter (DM) particle. This has led to increasing interest in alternatives to the Lambda CDM (LCDM) paradigm and alternative DM scenarios (including fuzzy DM, warm DM, self-interacting DM, etc.). In many of these scenarios, DM particles cannot be detected directly and constraints on their properties can ONLY be arrived at using astrophysical observations. The Dark Energy Spectroscopic Instrument (DESI) is currently one of the most powerful instruments for wide-field surveys. The synergy of DESI with ESA's Gaia satellite and future observing facilities will yield datasets of unprecedented size and coverage that will enable constraints on DM over a wide range of physical and mass scales and across redshifts. DESI will obtain spectra of the Lyman-alpha forest out to z~5 by detecting about 1 million QSO spectra that will put constraints on clustering of the low-density intergalactic gas and DM halos at high redshift. DESI will obtain radial velocities of 10 million stars in the Milky Way (MW) and Local Group satellites enabling us to constrain their global DM distributions, as well as the DM distribution on smaller scales. The paradigm of cosmological structure formation has been extensively tested with simulations. However, the majority of simulations to date have focused on collisionless CDM. Simulations with alternatives to CDM have recently been gaining ground but are still in their infancy. While there are numerous publicly available large-box and zoom-in simulations in the LCDM framework, there are no comparable publicly available WDM, SIDM, FDM simulations. DOE support for a public simulation suite will enable a more cohesive community effort to compare observations from DESI (and other surveys) with numerical predictions and will greatly impact DM science.

preprint2020arXiv

MANTRA: A Machine Learning reference lightcurve dataset for astronomical transient event recognition

We introduce MANTRA, an annotated dataset of 4869 transient and 71207 non-transient object lightcurves built from the Catalina Real Time Transient Survey. We provide public access to this dataset as a plain text file to facilitate standardized quantitative comparison of astronomical transient event recognition algorithms. Some of the classes included in the dataset are: supernovae, cataclysmic variables, active galactic nuclei, high proper motion stars, blazars and flares. As an example of the tasks that can be performed on the dataset we experiment with multiple data pre-processing methods, feature selection techniques and popular machine learning algorithms (Support Vector Machines, Random Forests and Neural Networks). We assess quantitative performance in two classification tasks: binary (transient/non-transient) and eight-class classification. The best performing algorithm in both tasks is the Random Forest Classifier. It achieves an F1-score of 96.25% in the binary classification and 52.79% in the eight-class classification. For the eight-class classification, non-transients ( 96.83% ) is the class with the highest F1-score, while the lowest corresponds to high-proper-motion stars ( 16.79% ); for supernovae it achieves a value of 54.57% , close to the average across classes. The next release of MANTRA includes images and benchmarks with deep learning models.

preprint2020arXiv

The cosmic web through the lens of graph entropy

We explore the information theory entropy of a graph as a scalar to quantify the cosmic web. We find entropy values in the range between 1.5 and 3.2 bits. We argue that this entropy can be used as a discrete analogue of scalars used to quantify the connectivity in continuous density fields. After showing that the entropy clearly distinguishes between clustered and random points, we use simulations to gauge the influence of survey geometry, cosmic variance, redshift space distortions, redshift evolution, cosmological parameters and spatial number density. Cosmic variance shows the least important influence while changes from the survey geometry, redshift space distortions, cosmological parameters and redshift evolution produce larger changes on the order of $10^{-2}$ bits. The largest influence on the graph entropy comes from changes in the number density of clustered points. As the number density decreases, and the cosmic web is less pronounced, the entropy can diminish up to 0.2 bits. The graph entropy is simple to compute and can be applied both to simulations and observational data from large galaxy redshift surveys; it is a new statistic that can be used in a complementary way to other kinds of topological or clustering measurements.

preprint2016arXiv

Boosting Lya and HeII 1640A Line Fluxes from Pop III Galaxies: Stochastic IMF Sampling and Departures from Case-B

We revisit calculations of nebular hydrogen Lya and HeII 1640A line strengths for population III galaxies, undergoing continuous and bursts of star formation. We focus on initial mass functions (IMFs) motivated by recent theoretical studies, which generally span a lower range of stellar masses than earlier works. We also account for case-B departures and the stochastic sampling of the IMF. In agreement with previous works, we find that departures from case-B can enhance the Lya flux by a factor of a few, but we argue that this enhancement is driven mainly by collisional excitation and ionization, and not due to photoionization from the n = 2 state of atomic hydrogen. The increased sensitivity of the Lya flux to the high-energy end of the galaxy spectrum makes it more subject to stochastic sampling of the IMF. The latter introduces a dispersion in the predicted nebular line fluxes around the deterministic value by as much as a factor of ~4. In contrast, the stochastic sampling of the IMF has less impact on the emerging Lyman Werner (LW) photon flux. When case-B departures and stochasticity effects are combined, nebular line emission from population III galaxies can be up to one order of magnitude brighter than predicted by 'standard' calculations that do not include these effects. This enhances the prospects for detection with future facilities such as JWST and large, groundbased telescopes.

preprint2016arXiv

Impact of Cosmic Variance on the Galaxy-Halo Connection for Lyman-$α$ Emitters

In this paper we study the impact of cosmic variance and observational uncertainties in constraining the mass and occupation fraction, $f_{\rm occ}$, of dark matter halos hosting Ly-$α$ Emitting Galaxies (LAEs) at high redshift. To this end, we construct mock catalogs from an N-body simulation to match the typical size of observed fields at $z=3.1$ ($\sim 1 {\rm deg^2}$). In our model a dark matter halo with mass in the range $M_{\rm min}<M_{\mathrm h}<M_{\rm max}$ can only host one detectable LAE at most. We proceed to explore the parameter space determined by $M_{\rm min}$,$M_{\rm max}$ and $f_{\rm occ}$ with a Markov Chain Monte-Carlo algorithm using the angular correlation function (ACF) and the LAEs number density as observational constraints. We find that the preferred minimum and maximum masses in our model span a wide range $10^{10.0}h^{-1}{\rm{M_{\odot}}}\leq M_{\rm min} \leq 10^{11.1}h^{-1}{\rm{M_{\odot}}}$ , $10^{11.0}h^{-1}{\rm{M_{\odot}}}\leq M_{\rm max} \leq 10^{13.0}h^{-1}{\rm{M_{\odot}}}$; followed by a wide range in the occupation fraction $0.02\leq f_{\rm occ} \leq 0.30$. As a consequence the median mass, $M_{50}$, of all the consistent models has a large uncertainty $M_{50} = 3.16^{+9.34}_{-2.37}\times 10^{10}$$h^{-1}{\rm{M_{\odot}}}$. However, we find that the same individual models have a relatively tight $1σ$ scatter around the median mass $ΔM_{1σ} = 0.55^{+0.11}_{-0.31}$ dex. We are also able to show that \focc\ is uniquely determined by $M_{\rm min}$, regardless of $M_{\rm max}$. We argue that upcoming large surveys covering at least $25$ deg$^{2}$ should be able to put tighter constraints on $M_{\rm min}$ and $f_{\rm occ}$ through the LAE number density distribution width constructed over several fields of $\sim 1$ deg$^{2}$.

preprint2016arXiv

The DESI Experiment Part I: Science,Targeting, and Survey Design

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$α$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

preprint2016arXiv

The DESI Experiment Part II: Instrument Design

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= λ/Δλ$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

preprint2015arXiv

Cosmological constraints from the redshift dependence of the Alcock-Paczynski test: galaxy density gradient field

We propose a method based on the redshift dependence of the Alcock-Paczynski (AP) test to measure the expansion history of the Universe. It uses the isotropy of the galaxy density gradient field to constrain cosmological parameters. If the density parameter $Ω_m$ or the dark energy equation of state $w$ are incorrectly chosen, the gradient field appears to be anisotropic with the degree of anisotropy varying with redshift. We use this effect to constrain the cosmological parameters governing the expansion history of the Universe. Although redshift-space distortions (RSD) induced by galaxy peculiar velocities also produce anisotropies in the gradient field, these effects are close to uniform in magnitude over a large range of redshift. This makes the redshift variation of the gradient field anisotropy relatively insensitive to the RSD. By testing the method on mock surveys drawn from the Horizon Run 3 cosmological N-body simulations, we demonstrate that the cosmological parameters can be estimated without bias. Our method is complementary to the baryon acoustic oscillation or topology methods as it depends on $D_AH$, the product of the angular diameter distance and the Hubble parameter.

preprint2015arXiv

Tensor anisotropy as a tracer of cosmic voids

We present a new method to find voids in cosmological simulations based on the tidal and the velocity shear tensors definitions of the cosmic web. We use the fractional anisotropy (FA) computed from the eigenvalues of each web scheme as a void tracer. We identify voids using a watershed transform based on the local minima of the FA field without making any assumption on the shape or structure of the voids. We test the method on the Bolshoi simulation and report on the abundance and radial averaged profiles for the density, velocity and fractional anisotropy. We find that voids in the velocity shear web are smaller than voids in the tidal web, with a particular overabundance of very small voids in the inner region of filaments/sheets. We classify voids as subcompensated/overcompansated depending on the absence/presence of an overdense matter ridge in their density profile, finding that close to $65\%$ and $35\%$ of the total population are classified into each category, respectively. Finally, we find evidence for the existence of universal profiles from the radially averaged profiles for density, velocity and fractional anisotropy. This requires that the radial coordinate is normalized to the effective radius of each void. Put together, all these results show that the FA is a reliable tracer for voids, which can be used in complementarity to other existing methods and tracers.

preprint2014arXiv

Cosmic web alignments with the shape, angular momentum and peculiar velocities of dark matter haloes

We study the alignment of dark matter haloes with the cosmic web characterized by the tidal and velocity shear fields. We focus on the alignment of their shape, angular momentum and peculiar velocities. We use a cosmological N-body simulation that allows to study dark matter halos spanning almost five orders of magnitude in mass ($10^{9}$-$10^{14}$) $h^{-1}$$M_{\odot}$ and spatial scales of $(0.5$-$1.0)$ $h^{-1}$ Mpc to define the cosmic web. We find that the halo shape presents the strongest alignment along the smallest tidal eigenvector, e.g. along filaments and walls, with a signal that gets stronger as the halo mass increases. In the case of the velocity shear field only massive halos $>10^{12}$ $h^{-1}$$M_{\odot}$ tend to have their shapes aligned along the largest tidal eigenvector; that is, perpendicular to filaments and walls. For the angular momentum we find alignment signals only for halos more massive than $10^{12}$ $h^{-1}$$M_{\odot}$ both in the tidal and velocity shear webs where the preferences are for it to be parallel to the middle eigenvector; perpendicular to filaments and parallel to walls. Finally, the peculiar velocities show a strong alignment along the smallest tidal eigenvector for all halo masses; halos move along filaments and walls. In the velocity shear the same alignment is present but weaker and only for haloes less massive than $10^{12}$ $h^{-1}$$M_{\odot}$. Our results clearly show that the two different algorithms we used to define the cosmic web describe different physical aspects of non-linear collapse and should be used in a complementary way to understand the effect of the cosmic web on galaxy evolution.

preprint2014arXiv

Systematic uncertainties from halo asphericity in dark matter searches

Although commonly assumed to be spherical, dark matter halos are predicted to be non-spherical by N-body simulations and their asphericity has a potential impact on the systematic uncertainties in dark matter searches. The evaluation of these uncertainties is the main aim of this work, where we study the impact of aspherical dark matter density distributions in Milky-Way-like halos on direct and indirect searches. Using data from the large N-body cosmological simulation Bolshoi, we perform a statistical analysis and quantify the systematic uncertainties on the determination of local dark matter density and the so-called $J$ factors for dark matter annihilations and decays from the galactic center. We find that, due to our ignorance about the extent of the non-sphericity of the Milky Way dark matter halo, systematic uncertainties can be as large as 35%, within the 95% most probable region, for a spherically averaged value for the local density of 0.3-0.4 GeV/cm$^3$. Similarly, systematic uncertainties on the $J$ factors evaluated around the galactic center can be as large as 10% and 15%, within the 95% most probable region, for dark matter annihilations and decays, respectively.

preprint2014arXiv

Target allocation yields for massively multiplexed spectroscopic surveys with fibers

We present Simulated Annealing fiber-to-target allocation simulations for the proposed DESI and 4MOST massively multiplexed spectroscopic surveys, and for both Poisson and realistically clustered mock target samples. We simulate both Echidna and theta-phi actuator designs, including the restrictions caused by the physical actuator characteristics during repositioning. For DESI, with theta-phi actuators, used in 5 passes over the sky for a mock ELG/LRG/QSO sample, with matched fiber and target densities, a total target allocation yield of 89.3% was achieved, but only 83.7% for the high-priority Ly-alpha QSOs. If Echidna actuators are used with the same pitch and number of passes, the yield increases by 5.7% and 16% respectively. Echidna also allows a factor-of-two increase in the number of close Ly-alpha QSO pairs that can be observed. Echidna spine tilt causes a variable loss of throughput, with average loss being the same as the loss at the rms tilt. With a natural tilt minimization scheme, we find an rms tilt always close to 0.58 x maximum. There is an additional but much smaller defocus loss, equivalent to an average defocus of 30microns. These tilt losses offset the gains in yield for Echidna, but because the survey strategy is driven by the higher priority targets, a clear survey speed advantage remains. For 4MOST, high and low latitude sample mock catalogs were supplied by the 4MOST team, and allocations were carried out with the proposed Echidna-based positioner geometry. At high latitudes, the resulting target completeness was 85.3% for LR targets and 78.9% for HR targets. At low latitude, the target completeness was 93.9% for LR targets and 71.2% for HR targets.

preprint2014arXiv

The impact of gas bulk rotation on the lyman-alpha line

We present results of radiative transfer calculations to measure the impact of gas bulk rotation on the morphology of the Lyman $α$ emission line in distant galaxies. We model a galaxy as a sphere with an homogeneous mixture of dust and hydrogen at a constant temperature. These spheres undergo solid-body rotation with maximum velocities in the range 0-300 \kms and neutral hydrogen optical depths in the range $τ_{\rm H}=10^{5}-10^{7}$. We consider two types of source distributions in the sphere: central and homogeneous. Our main result is that rotation introduces a dependence of the line morphology with viewing angle and rotational velocity. Observations with a line of sight parallel to the rotation axis yield line morphologies similar to the static case. For lines of sight perpendicular to the rotation axis both the intensity at the line center and the line width increase with rotational velocity. Along the same line of sight, the line becomes single peaked at rotational velocities close to half the line width in the static case. Notably, we find that rotation does not induce any spatial anisotropy in the integrated line flux, the escape fraction or the average number of scatterings. This is because Lyman α scattering through a rotating solid-body proceeds identical as in the static case. The only difference is the doppler shift from the different regions in the sphere that move with respect to the observer. This allows us to derive an analytic approximation for the viewing-angle dependence of the emerging spectrum, as a function of rotational velocity.

preprint2014arXiv

The Local Group in the cosmic web

We explore the characteristics of the cosmic web around Local Group(LG) like pairs using a cosmological simulation in the $Λ$CDM cosmology. We use the Hessian of the gravitational potential to classify regions on scales of $\sim 2$ Mpc as a peak, sheet, filament or void. The sample of LG counterparts is represented by two samples of halo pairs. The first is a general sample composed by pairs with similar masses and isolation criteria as observed for the LG. The second is a subset with additional observed kinematic constraints such as relative pair velocity and separation. We find that the pairs in the LG sample with all constraints are: (i) Preferentially located in filaments and sheets, (ii) Located in in a narrow range of local overdensity $0<δ<2$, web ellipticity $0.1<e<1.0$ and prolateness $-0.4<p<0.4$. (iii) Strongly aligned with the cosmic web. The alignments are such that the pair orbital angular momentum tends to be perpendicular to the smallest tidal eigenvector, $\hat{e}_3$, which lies along the filament direction or the sheet plane. A stronger alignment is present for the vector linking the two halos with the vector $\hat{e}_3$. Additionally, we fail to find a strong correlation of the spin of each halo in the pair with the cosmic web. All these trends are expected to a great extent from the selection on the LG total mass on the general sample. Applied to the observed LG, there is a potential conflict between the alignments of the different planes of satellites and the numerical evidence for satellite accretion along filaments; the direction defined by $\hat{e}_3$. This highlights the relevance of achieving a precise characterization of the place of the LG in the cosmic web in the cosmological context provided by $Λ$CDM.

preprint2014arXiv

The place of the Local Group in the cosmic web

We use the Bolshoi Simulation to find the most probable location of the Local Group (LG) in the cosmic web. Our LG simulacra are pairs of halos with isolation and kinematic properties consistent with observations. The cosmic web is defined using a tidal tensor approach. We find that the LG's preferred location is regions with a dark matter overdensity close to the cosmic average. This makes filaments and sheets the preferred environment. We also find a strong alignment between the LG and the cosmic web. The orbital angular momentum is preferentially perpendicular to the smallest tidal eigenvector, while the vector connecting the two halos is strongly aligned along the smallest tidal eigenvector and perpendicular to the largest tidal eigenvector; the pair lies and moves along filaments and sheets. We do not find any evidence for an alignment between the spin of each halo in the pair and the cosmic web.

preprint2013arXiv

The kinematics of the Local Group in a cosmological context

Recent observations constrained the tangential velocity of M31 with respect to the Milky Way (MW) to be v_tan<34.4 km/s and the radial velocity to be in the range v_rad=-109+/- 4.4 km/s (van der Marel et al. 2012). In this study we use a large volume high resolution N-body cosmological simulation (Bolshoi) together with three constrained simulations to statistically study this kinematics in the context of the LCDM. The comparison of the ensembles of simulated pairs with the observed LG at the 1-sigma level in the uncertainties has been done with respect to the radial and tangential velocities, the reduced orbital energy (e_tot), angular momentum (l_orb) and the dimensionless spin parameter, lambda. Our main results are: (i) the preferred radial and tangential velocities for pairs in LCDM are v_rad=-80+/-20 km/s, v_tan=50+/-10 km/s, (ii) pairs around that region are 3 to 13 times more common than pairs within the observational values, (iii) 15%to 24% of LG-like pairs in LCDM have energy and angular momentum consistent with observations while (iv) 9% to 13% of pairs in the same sample show similar values in the inferred dimensionless spin parameter. It follows that within current observational uncertainties the quasi-conserved quantities that characterize the orbit of the LG, i.e. e_tot, r_orb and lambda, do not challenge the standard LCDM model, but the model is in tension with regard to the actual values of the radial and tangential velocities. This might hint to a problem of the LCDM model to reproduce the observed LG.

preprint2012arXiv

A kinematic classification of the cosmic web

A new approach for the classification of the cosmic web is presented. In extension of the previous work of Hahn et al. (2007) and Forero-Romero et al. (2009) the new algorithm is based on the analysis of the velocity shear tensor rather than the gravitational tidal tensor. The procedure consists of the construction of the the shear tensor at each (grid) point in space and the evaluation of its three eigenvectors. A given point is classified to be either a void, sheet, filament or a knot according to the number of eigenvalues above a certain threshold, 0, 1, 2, or 3 respectively. The threshold is treated as a free parameter that defines the web. The algorithm has been applied to a dark matter only, high resolution simulation of a box of side-length 64$h^{-1}$Mpc and N = $1024^3$ particles with the framework of the WMAP5/LCDM model. The resulting velocity based cosmic web resolves structures down to <0.1$h^{-1}$Mpc scales, as opposed to the ~1$h^{-1}$Mpc scale of the tidal based web. The under-dense regions are made of extended voids bisected by planar sheets, whose density is also below the mean. The over-dense regions are vastly dominated by the linear filaments and knots. The resolution achieved by the velocity based cosmic web provides a platform for studying the formation of halos and galaxies within the framework of the cosmic web.

preprint2012arXiv

Effects of Star Formation Stochasticity on the Ly-alpha & Lyman Continuum Emission from Dwarf Galaxies during Reionization

Observations of distant galaxies play a key role in improving our understanding of the Epoch of Reionization (EoR). The observed Ly-alpha emission line strength - quantified by its restframe equivalent width (EW) - provides a valuable diagnostic of stellar populations and dust in galaxies during and after the EoR. In this paper we quantify the effects of star formation stochasticity on the predicted Ly-alpha EW in dwarf galaxies, using the publicly available code SLUG ('Stochastically Light Up Galaxies'). We compute the number of hydrogen ionizing photons, as well as flux in the Far UV for a set of models with star formation rates (SFR) in the range 10-3-1 Msol/yr. From these fluxes we compute the luminosity, L-alpha, and the EW of the Ly-alpha line. We find that stochasticity alone induces a broad distribution in L-alpha and EW at a fixed SFR, and that the widths of these distributions decrease with increasing SFR. We parameterize the EW probability density function (PDF) as an SFR-dependent double power law. We find that it is possible to have EW as low as ~EW0/4 and as high as ~3 times the EW0, where EW0 denotes the expected EW in the absence of stochasticity. We argue that stochasticity may therefore be important when linking drop-out and narrow-band selected galaxies, when identifying population III galaxies, and that it may help to explain the large EW (EW > 100 - 200 A) observed for a fraction of Ly-alpha- selected galaxies. Finally, we show that stochasticity can also affect the inferred escape fraction of ionizing photons from galaxies. In particular, we argue that stochasticity may simultaneously explain the observed anomalous ratios of the Lyman continuum flux density to the (non-ionizing) UV continuum density in so-called Lyman-Bump galaxies at z = 3.1, as well as the absence of such objects among a sample of z = 1.3 drop-out galaxies.

preprint2012arXiv

Towards a panchromatic picture of galaxy evolution during the reionization epoch

There are thousands of confirmed detections of star forming galaxies at high redshift (z > 4). These observations rely primarily on the detection of the spectral Lyman Break and the Lyman-alpha emission line. Theoretical modelling of these sources helps to interpret the observations in the framework of the standard cosmological paradigm. We present results from the High-z MareNostrum Project, aimed at constructing a panchromatic picture of the high redshift galaxy evolution that will improve our understanding of young star forming galaxies. Our simulation successfully reproduces the observational constraints from Lyman Break Galaxies and Lyman-alpha emitters at 5 < z < 7 . Based on this model we make predictions on the expected Far Infrared (FIR) emission that should be observed for LAEs. These predictions will help to settle down the question on the dust content of massive high-z galaxies, an issue that will be feasible to probe observationally with the Atacama Large Millimetre Array (ALMA).

preprint2011arXiv

Visualising Matter and Cosmologies: A Transhistorical Example

We propose a connection between the visualisation of cosmic matter and structure formation in the Cartesian tradition and that used by contemporary astrophysics. More precisely, we identify cosmological simulations of large scale structure in the Universe with the system of vortices in Descartes physics. This connection operates at different levels of the images: their representational purpose; the theoretical systems behind their use; and, finally, their function and materiality as visual productions. A skilled use of image analysis is necessary to stress the continuities and peculiarities between different epochs and disciplines.

preprint2008arXiv

The Coarse Geometry of Merger Trees in ΛCDM

We introduce the contour process to describe the geometrical properties of merger trees. The contour process produces a one-dimensional object, the contour walk, which is a translation of the merger tree. We portray the contour walk through its length and action. The length is proportional to to the number of progenitors in the tree, and the action can be interpreted as a proxy of the mean length of a branch in a merger tree. We obtain the contour walk for merger trees extracted from the public database of the Millennium Run and also for merger trees constructed with a public Monte-Carlo code which implements a Markovian algorithm. The trees correspond to halos of final masses between 10^{11} h^{-1} M_sol and 10^{14} h^{-1} M_sol. We study how the length and action of the walks evolve with the mass of the final halo. In all the cases, except for the action measured from Markovian trees, we find a transitional scale around 3 \times 10^{12} h^{-1} M_sol. As a general trend the length and action measured from the Markovian trees show a large scatter in comparison with the case of the Millennium Run trees.