Source author record

Michele Trenti

Michele Trenti 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

53works
5topics
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

53 published item(s)

preprint2026arXiv

Accelerated evolution of galaxy host halo masses during Cosmic Dawn from deep JWST clustering

We present the deepest clustering analysis of early galaxies to date, analyzing $N_{\rm{g}} \simeq 6500$ photometrically-selected Lyman Break Galaxies from JWST's Advanced Deep Extragalactic Survey (JADES) to reveal how galaxies and dark matter evolved during cosmic dawn ($5 \leq z < 11$). Using halo occupation distribution (HOD) modeling of the two-point angular correlation function, we trace the galaxy-halo relationships across the first billion years of cosmic history. Our analysis reveals that galaxies at $z = 10.6$ reside in dark matter halos over an order of magnitude less massive ($M_{\rm{h}} \sim 10^{10.12} M_{\odot}$) than their counterparts at $z = 5.5$ ($M_{\rm{h}} \sim 10^{11.45} M_{\odot}$), while exhibiting correspondingly higher effective bias values ($b_{\rm{g}}^{\rm{eff}} = 8.13^{+0.04}_{-0.02}$ compared to $5.64^{+0.10}_{-0.13}$). Correspondingly, the satellite galaxy fraction hints at a declining trend with decreasing redshift, reaching $<1\%$ by $z \sim 5-6$. However, the significant systematic and random uncertainties in the data-model comparison prevent us from drawing robust conclusions on the evolution - if any - of the satellite fraction during the epoch of reionization. These results provide the first view of the coevolution between galaxies and dark matter evolved at redshift $\gtrsim 10$, offering additional and independent constraints on early galaxy formation models tuned to reproducing luminosity function evolution.

preprint2026arXiv

BEACON: JWST NIRCam Pure-parallel Imaging Survey. II. Physical Properties of $z=7-14$ Galaxies

We present photometric properties of 161 galaxy candidates at $z=7-14$ selected from the second data release (DR2) of BEACON, a JWST Cycle 2 pure-parallel NIRCam imaging program. Carefully selected from 36 independent pointings (corresponding to $\sim350$\,arcmin$^2$ sky coverage), and hence with reduced cosmic variance, our galaxy candidates provide an unbiased sample for investigating galaxy properties over a wide range of environments. We measure the physical properties, including UV continuum slope ($β_{\rm UV}$), stellar mass ($M_*$), star formation rate (SFR), and sizes. Our highest redshift galaxy candidate at $z=13.71\pm0.15$ has a remarkably bright UV luminosity of $M_{\rm UV}=-21.19\pm0.08$, making it the brightest galaxy at $z>12$ if spectroscopically confirmed. With an extremely blue UV slope, compact morphology, and high star formation rate surface density ($Σ_{\rm SFR}$), this candidate may have extremely low metallicity, high ionizing photon escape fraction, or contributions from an AGN. Among our multiple independent sightlines, we identify three fields of galaxy number overdensity with $>3σ$ significance. The properties of galaxies in various environments do not exhibit significant differences, implying either that accelerated galaxy evolution in overdense regions is not yet widespread at $z>7$, or that the current constraints are limited by sample size. Our simulations indicate that increasing the sample by an order of magnitude would allow such environmental trends to be robustly confirmed or ruled out, underscoring the importance of future pure-parallel observations.

preprint2026arXiv

Beyond UV: Rest-frame B-band and Apparent Luminosity Functions of z=5-9 Galaxies

We present new measurements of galaxy luminosity functions (LFs) from JWST/NIRCam imaging over the redshift range z=4.5-9.7, using photometric catalogs from JADES and public extragalactic fields. Our analysis includes rest-frame UV and B-band LFs, as well as apparent LFs in F090W, F115W, F200W, F356W, and F444W. We present the first constraints on the rest-frame B-band LF at z~7-8 and extend existing measurements at z~5 to M(B) = -18 mag. The B-band LFs evolve more strongly with redshift than UV LFs, though both decline more gradually than predicted by simulations at z>5. No single existing simulation reproduces all observed trends, with discrepancies likely driven by assumptions about binary evolution and stellar population synthesis models. The apparent LFs in F356W and F444W show hints of a bright-end excess at all redshifts, extending to fainter magnitudes at higher redshift. While extreme emission line galaxies may partially account for it, the excess may also indicate a population of moderately red, optically bright sources - potentially dusty star-forming galaxies or obscured AGNs. Finally, we find that rest-frame B-band luminosity correlates more tightly with stellar mass than UV, making it a powerful tracer of mass assembly and reinforcing the diagnostic value of rest-frame optical LFs in uncovering the physical processes that drive early galaxy formation.

preprint2026arXiv

The BoRG-$JWST$ Survey: Abundance and Mass-to-light Ratio of Luminous $z=7-9$ Galaxies from Independent Sight Lines with NIRSpec

We present new results on the rest-frame UV luminosity function (UVLF) and stellar mass-to-light (M/L) ratio of bright (M$_{\rm UV}\lesssim-20$ mag) spectroscopically-confirmed galaxies at $z=7-9$ derived from the BoRG-$JWST$ survey, a unique data set of NIRSpec prism follow up of $HST$-selected sources from random-pointing imaging. By selecting galaxies from over 200 independent sight lines, the survey minimizes cosmic variance ensuring a statistically robust sample of the bright-galaxy population during the epoch of reionization. The data is used to constrain, for the first time, the bright end of the UVLF at $z=7-9$ from spectroscopically-confirmed galaxies over eight independent fields. We find that the bright end of the UVLF is higher than found using imaging over $JWST$ legacy fields, suggesting the latter may be significantly affected by cosmic variance, and thus reducing the tension with recent findings from $JWST$ at $z>10$ and comparable to models invoking little dust attenuation and bursty star formation. Additionally, we use the galaxies' $JWST$ spectra to infer their stellar masses and M/L ratios relative to other $HST$ and $JWST$ studies. We show that the stellar mass scales almost linearly with UV luminosity (M$_* \propto L_{\rm UV}^{0.85\pm0.12}$), albeit with large ($\sim0.5$ dex) intrinsic scatter, consistent with stochastic bursts of star formation in early galaxy formation.

preprint2024arXiv

Early Results from GLASS-JWST XXIII: The transmission of Lyman-alpha from UV-faint z ~ 3-6 galaxies

Lyman-alpha (Ly$α$) emission from galaxies can be used to trace neutral hydrogen in the epoch of reionization, however, there is a degeneracy between the attenuation of Ly$α$ in the intergalactic medium (IGM) and the line profile emitted from the galaxy. Large shifts of Ly$α$ redward of systemic due to scattering in the interstellar medium can boost Ly$α$ transmission in the IGM during reionization. The relationship between Ly$α$ velocity offset from systemic and other galaxy properties is not well-established at high-redshift or low luminosities, due to the difficulty of observing emission lines which trace systemic redshift. Rest-frame optical spectroscopy with JWST/NIRSpec has opened a new window into understanding of Ly$α$ at z>3. We present a sample of 12 UV-faint galaxies ($-20 \lesssim$ MUV $\lesssim -16$) at $3 \lesssim z \lesssim 6$, with Ly$α$ velocity offsets, $Δv_{\mathrm{Ly}α}$, measured from VLT/MUSE and JWST/NIRSpec from the GLASS-JWST Early Release Program. We find median $Δv_{\mathrm{Ly}α}$ of 205 km s$^{-1}$ and standard deviation 75 km s$^{-1}$, compared to 320 and 170km s$^{-1}$ for MUV < -20 galaxies in the literature. Our new sample demonstrates the previously observed trend of decreasing Ly$α$ velocity offset with decreasing UV luminosity and optical line velocity dispersion, extends to MUV $\gtrsim$ -20, consistent with a picture where the Ly$α$ profile is shaped by gas close to the systemic redshift. Our results imply that during reionization Ly$α$ from UV-faint galaxies will be preferentially attenuated, but that detecting Ly$α$ with low $Δv_{\mathrm{Ly}α}$ can be an indicator of large ionized bubbles.

preprint2023arXiv

The brightest galaxies at Cosmic Dawn

Recent JWST observations suggest an excess of $z\gtrsim10$ galaxy candidates above most theoretical models. Here, we explore how the interplay between halo formation timescales, star formation efficiency and dust attenuation affects the properties and number densities of galaxies we can detect in the early universe. We calculate the theoretical upper limit on the UV luminosity function, assuming star formation is 100% efficient and all gas in halos is converted into stars, and that galaxies are at the peak age for UV emission (~10 Myr). This upper limit is ~4 orders of magnitude greater than current observations, implying these are fully consistent with star formation in $Λ$CDM cosmology. In a more realistic model, we use the distribution of halo formation timescales derived from extended Press-Schechter theory as a proxy for star formation rate (SFR). We predict that the galaxies observed so far at $z\gtrsim10$ are dominated by those with the fastest formation timescales, and thus most extreme SFRs and young ages. These galaxies can be upscattered by ~1.5 mag compared to the median UV magnitude vs halo mass relation. This likely introduces a selection effect at high redshift whereby only the youngest ($\lesssim$10 Myr), most highly star forming galaxies (specific SFR$\gtrsim$30 Gyr$^{-1}$) have been detected so far. Furthermore, our modelling suggests that redshift evolution at the bright end of the UV luminosity function is substantially affected by the build-up of dust attenuation. We predict that deeper JWST observations (reaching m~30) will reveal more typical galaxies with relatively older ages (~100 Myr) and less extreme specific SFRs (~10 Gyr$^{-1}$ for a $M_\mathrm{UV}$ ~ -20 galaxy at z~10).

preprint2022arXiv

A geostatistical analysis of multiscale metallicity variations in galaxies [II]: Predicting the metallicities of Hii and diffuse ionised gas regions via universal kriging

The metallicity of diffuse ionised gas (DIG) cannot be determined using strong emission line diagnostics, which are calibrated to calculate the metallicity of Hii regions. Because of this, resolved metallicity maps from integral field spectroscopy (IFS) data remain largely incomplete. In this paper (the second of a series), we introduce the geostatistical technique of universal kriging, which allows the complete 2D metallicity distribution of a galaxy to be reconstructed from metallicities measured at Hii regions, accounting for spatial correlations between nearby data points. We apply this method to construct high-fidelity metallicity maps of the local spiral galaxy NGC 5236 using data from the TYPHOON/PrISM survey. We find significant correlation in the metallicity of Hii regions separated by up to 0.4-1.2 kpc. Predictions constructed using this method were tested using cross-validation in Hii regions, and we show that they outperform significantly interpolation based on metallicity gradients. Furthermore, we apply kriging to predict the metallicities in regions dominated by DIG emission, considering seven additional spiral galaxies with high resolution (<100pc) metallicity maps. We compare kriging maps to DIG metallicities computed with novel ionisation corrections, and find that such corrections introduce a systematic offset of up to $\pm0.1$ dex for any individual galaxy, with a scatter of 0.02-0.07 dex for the sample. Overall we recommend universal kriging, together with a calibrated geostatistical model, as the superior method for inferring the metallicities of DIG-dominated regions in local spiral galaxies, demonstrating further the potential of applying geostatistical methods to spatially resolved galaxy observations.

preprint2022arXiv

SkyHopper mission science case I: Identification of high redshift Gamma-Ray Bursts through space-based near-infrared afterglow observations

Long-duration gamma-ray burst (GRB) afterglow observations offer cutting-edge opportunities to characterise the star formation history of the Universe back to the epoch of reionisation, and to measure the chemical composition of interstellar and intergalactic gas through absorption spectroscopy. The main barrier to progress is the low efficiency in rapidly and confidently identifying which bursts are high redshift ($z > 5$) candidates before they fade, as this requires low-latency follow-up observations at near-infrared wavelengths (or longer) to determine a reliable photometric redshift estimate. So far this task has been performed by instruments on the ground, but sky visibility and weather constraints limit the number of GRB targets that can be observed and the speed at which follow-up is possible. In this work we develop a Monte Carlo simulation framework to investigate an alternative approach based on the use of a rapid-response near-infrared nano-satellite, capable of simultaneous imaging in four bands from $0.8$ to $1.7μ$m (a mission concept called SkyHopper). We find that such a nano-satellite is capable of detecting in the H band (1.6 $μ$m) $72.5\% \pm 3.1\%$ of GRBs concurrently observable with the Swift satellite via its UVOT instrument (and $44.1\% \pm 12.3\%$ of high redshift ($z>5$) GRBs) within 60 minutes of the GRB prompt emission. This corresponds to detecting $\sim 55$ GRB afterglows per year, of which 1-3 have $z > 5$. These rates represent a substantial contribution to the field of high-$z$ GRB science, as only 23 $z > 5$ GRBs have been collectively discovered by the entire astronomical community over the last $\sim 24$ years. Additionally, we find that launching a mini-constellation of 3 near-infrared nano-satellites would increase the detection fraction of afterglows to $\sim 83\%$ and substantially reduce the latency in the photometric redshift determination.

preprint2022arXiv

The size-luminosity relation of lensed galaxies at $z=6-9$ in the Hubble Frontier Fields

We measure the size-luminosity relation of photometrically-selected galaxies within the redshift range $z\sim6-9$, using galaxies lensed by six foreground Hubble Frontier Fields (HFF) clusters. The power afforded by strong gravitational lensing allows us to observe fainter and smaller galaxies than in blank fields. We select our sample of galaxies and obtain their properties, e.g., redshift, magnitude, from the photometrically-derived ASTRODEEP catalogues. The intrinsic size is measured with the Lenstruction software, and completeness maps are created as a function of size and luminosity via the GLACiAR2 software. We perform a Bayesian analysis to estimate the intrinsic and incompleteness-corrected size-luminosity distribution, with parameterization $r_e \propto L^β$. We find slopes of $β\sim0.48\pm0.08$ at $z\sim6-7$ and $β\sim0.68\pm0.14$ at $z\sim8.5$, adopting the Bradac lens model. The slope derived by lensed galaxies is steeper than that obtained in blank fields and is consistent with other independent determinations of the size-luminosity relation from the HFF dataset. We also investigate the systematic uncertainties correlated with the choice of lens models, finding that the slopes of size-luminosity relations derived from different models are consistent with each other, i.e. the modeling errors are not a significant source of uncertainty in the size-luminosity relation.

preprint2021arXiv

RELICS-DP7: Spectroscopic Confirmation of a Dichromatic Primeval Galaxy at z ~ 7

We report the discovery of a spectroscopically-confirmed strong Lyman-$α$ emitter at $z=7.0281\pm0.0003$, observed as part of the Reionization Cluster Lensing Survey (RELICS). This galaxy, dubbed "Dichromatic Primeval Galaxy" at $z\sim7$ (DP7), shows two distinct components. While fairly unremarkable in terms of its ultraviolet (UV) luminosity ($\sim0.3L^{\ast}_{UV}$, where $L^{\ast}_{UV}$ is the characteristic luminosity), DP7 has one of the highest observed Lyman-$α$ equivalent widths (EWs) among Lyman-$α$ emitters at $z>6$ ($>200$ Angstrom in the rest frame). The strong Lyman-$α$ emission generally suggests a young metal-poor, low-dust galaxy; however, we find that the UV slope $β$ of the galaxy as a whole is redder than typical star-forming galaxies at these redshifts, $-1.13\pm 0.84$, likely indicating, on average, a considerable amount of dust obscuration, or an older stellar population. When we measure $β$ for the two components separately, however, we find evidence of differing UV colors, suggesting two separate stellar populations. Also, we find that Lyman-$α$ is spatially extended and likely larger than the galaxy size, hinting to the possible existence of a Lyman-$α$ halo. Rejuvenation or merging events could explain these results. Either scenario requires an extreme stellar population, possibly including a component of Population III stars, or an obscured Active Galactic Nucleus. DP7, with its low UV luminosity and high Lyman-$α$ EW, represents the typical galaxies that are thought to be the major contribution to the reionization of the Universe, and for this reason DP7 is an excellent target for follow-up with the James Webb Space Telescope.

preprint2021arXiv

The Physical Properties of Luminous $z\gtrsim8$ Galaxies and Implications for the Cosmic Star Formation Rate Density From $\sim$0.35 deg$^{2}$ of (Pure-)Parallel HST Observations

We present the largest systematic, HST-based search to date for luminous $z\gtrsim8$ galaxy candidates using $\sim$1267 arcmin$^{2}$ of (pure-)parallel observations from a compilation of 288 random sightlines with ACS and WFC3 observations, derived from the SuperBoRG data set and together representing a factor $\sim1.12\times$ larger than existing space-based data sets. Using NIR color cuts and careful photo-$z$ analyses, we find 31 $z\gtrsim8$ galaxy candidates over 29 unique sightlines, and derive global galaxy properties such as UV magnitudes and continuum slopes, sizes, and rest-frame optical properties (e.g., SFRs, stellar masses, $A_{\rm v}$). Taking advantage of the (pure-)parallel nature of our data set - making it one of the most representative thus far - and derived SFRs, we evaluate the cosmic star formation rate density for the bright end of the UV luminosity function at $z\sim8-10$ and test the validity of luminosity function-derived results using a conversion factor. We find our method yields comparable results to those derived with luminosity functions. Furthermore, we present follow up observations of 5 (Super)BoRG targets with Keck/MOSFIRE, finding no evidence of Ly$α$ in $>3$ hrs of $Y-$band observations in either, consistent with a largely neutral medium at $z\sim8$. Our results offer a definitive HST legacy on the bright end of the luminosity function and provide a valuable benchmark as well as targets for follow up with JWST.

preprint2020arXiv

One star, two stars, or both? Investigating metallicity-dependant models for Gamma-Ray Burst progenitors with the IllustrisTNG simulation

The rate of long-duration gamma ray bursts (GRBs) has been identified as a potential proxy for the star formation rate (SFR) across redshift, but the exact relationship depends on GRB progenitor models (single versus binary). The single-progenitor collapsar model accounts for the preference towards low-metallicity GRB progenitors, but is in apparent tension with some high-metallicity GRB host galaxy measurements. As a possible solution, we consider the scenario where high-metallicity GRB hosts harbour low metallicity regions in which GRB progenitors form. For this, we use the IllustrisTNG cosmological hydrodynamical simulation to investigate the internal metallicity distribution of GRB hosts, implementing in post-processing different GRB formation models. Predictions (GRB rate, host metallicities and stellar masses) are compared to the high-completeness GRB legacy surveys BAT6 and SHOALS and a sample of high-redshift GRB-DLA metallicities, allowing us to compute their relative likelihoods. When the internal metallicity distribution of galaxies is ignored, the best-fitting model requires a metallicity-independent channel, as previously proposed by Trenti, Perna & Jimenez. However, when the internal metallicity distribution is considered, a basic metallicity bias model with a cutoff at $Z_{max}=0.35Z_\odot$ is the best fitting one. Current data are insufficient to discriminate among more detailed metallicity bias models, such as weak metallicity dependence of massive binaries vs stronger metallicity bias of collapsars. An increased sample of objects, and direct measurements of host stellar masses at redshift $z>2$ would allow to further constrain the origin of long GRBs.

preprint2020arXiv

RELICS: A Very Large ($θ_{E}\sim40"$) Cluster Lens -- RXC J0032.1+1808

Extensive surveys with the \textit{Hubble Space Telescope} (HST) over the past decade, targeting some of the most massive clusters in the sky, have uncovered dozens of galaxy-cluster strong lenses. The massive cluster strong-lens scale is typically $θ_{E}\sim10\arcsec$ to $\sim30-35\arcsec$, with only a handful of clusters known with Einstein radii $θ_{E}\sim40\arcsec$ or above (for $z_{source}=2$, nominally). Here we report another very large cluster lens, RXC J0032.1+1808 ($z=0.3956$), the second richest cluster in the redMapper cluster catalog and the 85th most massive cluster in the Planck Sunyaev-Zel'dovich catalog. With our Light-Traces-Mass and fully parametric (dPIEeNFW) approaches, we construct strong lensing models based on 18 multiple images of 5 background galaxies newly identified in the \textit{Hubble} data mainly from the \textit{Reionization Lensing Cluster Survey} (RELICS), in addition to a known sextuply imaged system in this cluster. Furthermore, we compare these models to Lenstool and GLAFIC models that were produced independently as part of the RELICS program. All models reveal a large effective Einstein radius of $θ_{E}\simeq40\arcsec$ ($z_{source}=2$), owing to the obvious concentration of substructures near the cluster center. Although RXC J0032.1+1808 has a very large critical area and high lensing strength, only three magnified high-redshift candidates are found within the field targeted by RELICS. Nevertheless, we expect many more high-redshift candidates will be seen in wider and deeper observations with \textit{Hubble} or \emph{JWST}. Finally, the comparison between several algorithms demonstrates that the total error budget is largely dominated by systematic uncertainties.

preprint2020arXiv

Spectroscopically Confirmed Lyman-Alpha Emitters from Redshift 5 to 7 Behind Ten Galaxy Cluster Lenses

We present 36 spectroscopically confirmed intrinsically UV-faint Ly$α$ emitting galaxies from follow-up observations with Keck/DEIMOS of gravitationally lensed high-redshift candidates. Candidates were selected to be between $5\lesssim z \lesssim 7$ from photometric data using \textit{HST} and \textit{Spitzer} imaging surveys. We used photometric redshift information to perform an integrated photometric redshift probability cut $>1\%$ between $5<z<7$ to construct a sample of 198 high-redshift objects. Our high-redshift sample spans intrinsic UV luminosities from a few $L^*$ down to $0.001L^*$. We identified 19 high-confidence detections of Ly$α$ and an additional 17 likely detections. We divided our sample into lower-redshift ($z\sim5.5$) and higher-redshift ($z\sim6.5$) bins and ran Monte Carlo trials, incorporating the strength of the Ly$α$ emission and the photometric redshift of the non-detections. Considering only objects where Ly$α$ could be detected at EW(Ly$α$)$>$25Å at $3σ$ at the fiducial depth of our survey, and only those galaxies with EW(Ly$α$)$>$25Å as true LAEs, and finally, only objects with $m_{AB}<26.8$, we found the LAE fraction to be flat, or modestly increase from 0.26$\pm0.04$ to 0.30$\pm0.04$. These values relative to those for lower-redshift samples are consistent with a rising LAE fraction with redshift out to $z\sim6$, but at $z\sim6.5$ there is some tension between our results and results from surveys at intrinsically brighter luminosities. We conclude intrinsically fainter galaxies have Ly$α$ emission, and there is a steep drop in the LAE fraction from our high-redshift sample at $z\sim6.5$ and from similar galaxies at $z\sim7.5$. This likely indicates we are witnessing the tail end of the epoch of reionization, as such a drop is not expected due to changes of intrinsic galaxy properties between these redshifts.

preprint2020arXiv

Stochastic Processes as the Origin of the Double-Power Law Shape of the Quasar Luminosity Function

The Quasar Luminosity Function (QLF) offers insight into the early co-evolution of black holes and galaxies. It has been characterized observationally up to redshift $z\sim6$ with clear evidence of a double power-law shape, in contrast to the Schechter-like form of the underlying dark-matter halo mass function. We investigate a physical origin for the difference in these distributions by considering the impact of stochasticity induced by the processes that determine the quasar luminosity for a given host halo and redshift. We employ a conditional luminosity function and construct the relation between median quasar magnitude versus halo mass $M_{UV,\rm{c}}(M_{\rm{h}})$ with log-normal in luminosity scatter $Σ$, and duty-cycle $ε_{\rm{DC}}$, and focus on high redshift $z\gtrsim4$. We show that, in order to reproduce the observed QLF, the $Σ=0$ abundance matching requires all of the brightest quasars to be hosted in the rarest most massive dark-matter halos (with an increasing $M_{UV,\rm{c}}/M_{\rm{h}}$ in halo mass). Conversely, for $Σ>0$ the brightest quasars can be over-luminous outliers hosted in relatively common dark-matter halos. In this case, the median quasar magnitude versus halo mass relation, $M_{UV,\rm{c}}$, flattens at the high-end, as expected in self-regulated growth due to feedback. We sample the parameter space of $Σ$ and $ε_{\rm{DC}}$ and show that $M_{UV,\rm{c}}$ flattens above $M_{\rm{h}}\sim 10^{12}M_{\odot}$ for $ε_{\rm{DC}}<10^{-2}$. Models with $ε_{\rm{DC}}\sim1$ instead require a high mass threshold close to $M_{\rm{h}}\gtrsim10^{13}M_{\odot}$. We investigate the impact of $ε_{\rm{DC}}$ and $Σ$ on measurements of clustering and find there is no luminosity dependence on clustering for $Σ>0.3$, consistent with recent observations from Subaru HSC.

preprint2020arXiv

The Impact of Scatter in the Galaxy UV Luminosity to Halo Mass Relation on Ly$α$ Visibility During the Epoch of Reionization

The reionization of hydrogen is closely linked to the first structures in the universe, so understanding the timeline of reionization promises to shed light on the nature of these early objects. In particular, transmission of Lyman alpha (Ly$α$) from galaxies through the intergalactic medium (IGM) is sensitive to neutral hydrogen in the IGM, so can be used to probe the reionization timeline. In this work, we implement an improved model of the galaxy UV luminosity to dark matter halo mass relation to infer the volume-averaged fraction of neutral hydrogen in the IGM from Ly$α$ observations. Many models assume that UV-bright galaxies are hosted by massive dark matter haloes in overdense regions of the IGM, so reside in relatively large ionized regions. However, observations and N-body simulations indicate that scatter in the UV luminosity-halo mass relation is expected. Here, we model the scatter (though we assume the IGM topology is unaffected) and assess the impact on Ly$α$ visibility during reionization. We show that UV luminosity-halo mass scatter reduces Ly$α$ visibility compared to models without scatter, and that this is most significant for UV-bright galaxies. We then use our model with scatter to infer the neutral fraction, $\overline{x}_{\mathrm{HI}}$, at $z \sim 7$ using a sample of Lyman-break galaxies in legacy fields. We infer $\overline{x}_{\mathrm{HI}} = 0.55_{-0.13}^{+0.11}$ with scatter, compared to $\overline{x}_{\mathrm{HI}} = 0.59_{-0.14}^{+0.12}$ without scatter, a very slight decrease and consistent within the uncertainties. Finally, we place our results in the context of other constraints on the reionization timeline and discuss implications for future high-redshift galaxy studies.

preprint2018arXiv

The Reionization Lensing Cluster Survey (RELICS) and the Brightest High-z Galaxies

Massive foreground galaxy clusters magnify and distort the light of objects behind them, permitting a view into both the extremely distant and intrinsically faint galaxy populations. We present here the z ~ 6 - 8 candidate high-redshift galaxies from the Reionization Lensing Cluster Survey (RELICS), a Hubble and Spitzer Space Telescope survey of 41 massive galaxy clusters spanning an area of ~200 arcmin^2. These clusters were selected to be excellent lenses and we find similar high-redshift sample sizes and magnitude distributions as CLASH. We discover 321 candidate galaxies with photometric redshifts between z ~ 6 to z ~ 8, including extremely bright objects with H-band magnitudes of m_AB ~ 23 mag. As a sample, the observed (lensed) magnitudes of these galaxies are among the brightest known at z> 6, comparable to much wider, blank-field surveys. RELICS demonstrates the efficiency of using strong gravitational lenses to produce high-redshift samples in the epoch of reionization. These brightly observed galaxies are excellent targets for follow-up study with current and future observatories, including the James Webb Space Telescope.

preprint2016arXiv

Detection of Lyman-Alpha Emission From a Triple Imaged z=6.85 Galaxy Behind MACS J2129.4-0741

We report the detection of Ly$α$ emission at $\sim9538$Å in the Keck/DEIMOS and \HST WFC3 G102 grism data from a triply-imaged galaxy at $z=6.846\pm0.001$ behind galaxy cluster MACS J2129.4$-$0741. Combining the emission line wavelength with broadband photometry, line ratio upper limits, and lens modeling, we rule out the scenario that this emission line is \oii at $z=1.57$. After accounting for magnification, we calculate the weighted average of the intrinsic Ly$α$ luminosity to be $\sim1.3\times10^{42}~\mathrm{erg}~\mathrm{s}^{-1}$ and Ly$α$ equivalent width to be $74\pm15$Å. Its intrinsic UV absolute magnitude at 1600Å is $-18.6\pm0.2$ mag and stellar mass $(1.5\pm0.3)\times10^{7}~M_{\odot}$, making it one of the faintest (intrinsic $L_{UV}\sim0.14~L_{UV}^*$) galaxies with Ly$α$ detection at $z\sim7$ to date. Its stellar mass is in the typical range for the galaxies thought to dominate the reionization photon budget at $z\gtrsim7$; the inferred Ly$α$ escape fraction is high ($\gtrsim 10$\%), which could be common for sub-$L^*$ $z\gtrsim7$ galaxies with Ly$α$ emission. This galaxy offers a glimpse of the galaxy population that is thought to drive reionization, and it shows that gravitational lensing is an important avenue to probe the sub-$L^*$ galaxy population.

preprint2016arXiv

The Close Stellar Companions to Intermediate Mass Black Holes

When embedded in dense cluster cores, intermediate mass black holes (IMBHs) acquire close stellar or stellar-remnant companions. These companions are not only gravitationally bound, they tend to hierarchically isolate from other cluster stars through series of multibody encounters. In this paper, we study the demographics of IMBH companions in compact star clusters through direct $N$-body simulation. We study clusters initially composed of $10^5$ or $2\times 10^5$ stars with IMBHs of 75 and 150 solar masses, and follow their evolution for 6-10 Gyr. A tight innermost binary pair of IMBH and stellar object rapidly forms. The IMBH has a companion with orbital semi-major axis at least three times tighter than the second-most bound object over 90% of the time. These companionships have typical periods of order years and are subject to cycles of exchange and destruction. The most frequently observed, long-lived pairings persist for $\sim 10^7$ yr. The demographics of IMBH companions in clusters are diverse; they include both main sequence, giant stars, and stellar remnants. Companion objects may reveal the presence of an IMBH in a cluster in one of several ways. Most-bound companion stars routinely suffer grazing tidal interactions with the IMBH, offering a dynamical mechanism to produce repeated flaring episodes like those seen in the IMBH candidate HLX-1. Stellar winds of companion stars provide a minimum quiescent accretion rate for IMBHs, with implications for radio searches for IMBH accretion in globular clusters. Finally, gravitational wave inspirals of compact objects are found to occur with promising frequency.

preprint2016arXiv

The Grism Lens-Amplified Survey from Space (GLASS). VI. Comparing the Mass and Light in MACSJ0416.1-2403 using Frontier Field imaging and GLASS spectroscopy

We present a strong and weak gravitational lens model of the galaxy cluster MACSJ0416.1-2403, constrained using spectroscopy from the Grism Lens-Amplified Survey from Space (GLASS) and Hubble Frontier Fields (HFF) imaging data. We search for emission lines in known multiply imaged sources in the GLASS spectra, obtaining secure spectroscopic redshifts of 31 multiple images belonging to 16 distinct source galaxies. The GLASS spectra provide the first spectroscopic measurements for 6 of the source galaxies. The weak lensing signal is acquired from 884 galaxies in the F606W HFF image. By combining the weak lensing constraints with 15 multiple image systems with spectroscopic redshifts and 9 multiple image systems with photometric redshifts, we reconstruct the gravitational potential of the cluster on an adaptive grid. The resulting total mass density map is compared with a stellar mass density map obtained from the deep Spitzer Frontier Fields imaging data to study the relative distribution of stellar and total mass in the cluster. We find that the projected stellar mass to total mass ratio, $f_{\star}$, varies considerably with the stellar surface mass density. The mean projected stellar mass to total mass ratio is $\langle f_{\star} \rangle= 0.009 \pm 0.003 $ (stat.), but with a systematic error as large as $0.004-0.005$, dominated by the choice of the IMF. We find agreement with several recent measurements of $f_{\star}$ in massive cluster environments. The lensing maps of convergence, shear, and magnification are made available to the broader community in the standard HFF format.

preprint2016arXiv

The Grism lens-amplified survey from space (GLASS). VII. The diversity of the distribution of star formation in cluster and field galaxies at 0.3<z<0.7

Exploiting the slitless spectroscopy taken as part of the Grism Lens-Amplified Survey from Space (GLASS), we present an extended analysis of the spatial distribution of star formation in 76 galaxies in 10 clusters at 0.3< z <0.7. We use 85 foreground and background galaxies in the same redshift range as a field sample. The samples are well matched in stellar mass (10^8-10^11 M_sun) and star formation rate (0.5-50 M_sun/yr). We visually classify galaxies in terms of broad-band morphology, Halpha morphology and likely physical process acting on the galaxy. Most Halpha emitters have a spiral morphology (41+/-8% in clusters, 51+/-8% in the field), followed by mergers/interactions (28+/-8%, 31+/-7%, respectively) and early-type galaxies (remarkably as high as 29+/-8% in clusters and 15+/-6% in the field). A diversity of Halpha morphologies is detected, suggesting a diversity of physical processes. In clusters, 30+/-8% of the galaxies present a regular morphology, mostly consistent with star formation diffused uniformly across the stellar population (mostly in the disk component, when present). The second most common morphology (28+/-8%) is asymmetric/jellyfish, consistent with ram pressure stripping or other non-gravitational processes in 18+/-8% of the cases. Ram pressure stripping appears significantly less prominent in the field (2+/-2%), where the most common morphology/mechanism appears to be consistent with minor gas rich mergers or clump accretion. This work demonstrates that while environment specific mechanisms affect galaxy evolution at this redshift, they are diverse and their effects subtle. A full understanding of this complexity requires larger samples and detailed and spatially resolved physical models.

preprint2015arXiv

Compact Stellar Binary Assembly in the First Nuclear Star Clusters and r-Process Synthesis in the Early Universe

Investigations of elemental abundances in the ancient and most metal deficient stars are extremely important because they serve as tests of variable nucleosynthesis pathways and can provide critical inferences of the type of stars that lived and died before them. The presence of r-process elements in a handful of carbon-enhanced metal-poor (CEMP-r) stars, which are assumed to be closely connected to the chemical yield from the first stars, is hard to reconcile with standard neutron star mergers. Here we show that the production rate of dynamically assembled compact binaries in high-z nuclear star clusters can attain a sufficient high value to be a potential viable source of heavy r-process material in CEMP-r stars. The predicted frequency of such events in the early Galaxy, much lower than the frequency of Type II supernovae but with significantly higher mass ejected per event, can naturally lead to a high level of scatter of Eu as observed in CEMP-r stars.

preprint2015arXiv

Connecting Faint End Slopes of the Lyman-$α$ emitter and Lyman-break Galaxy Luminosity Functions

We predict Lyman-$α$ (Ly$α$) luminosity functions (LFs) of Ly$α$-selected galaxies (Ly$α$ emitters, or LAEs) at $z=3-6$ using the phenomenological model of Dijkstra & Wyithe (2012). This model combines observed UV-LFs of Lyman-break galaxies (LBGs, or drop out galaxies), with constraints on their distribution of Ly$α$ line strengths as a function of UV-luminosity and redshift. Our analysis shows that while Ly$α$ LFs of LAEs are generally not Schechter functions, these provide a good description over the luminosity range of $\log_{10}( L_α/{\rm erg}\,{\rm s}^{-1})=41-44$. Motivated by this result, we predict Schechter function parameters at $z=3-6$. Our analysis further shows that (i) the faint end slope of the Ly$α$ LF is steeper than that of the UV-LF of Lyman-break galaxies, (with a median $α_{Lyα} < -2.0$ at $z\gtrsim 4$), and (ii) a turn-over in the Ly$α$ LF of LAEs at Ly$α$ luminosities $10^{40}$ erg s$^{-1}<L_α\lesssim 10^{41}$ erg s$^{-1}$ may signal a flattening of UV-LF of Lyman-break galaxies at $-12>M_{\rm UV}>-14$. We discuss the implications of these results - which can be tested directly with upcoming surveys - for the Epoch of Reionization.

preprint2015arXiv

Correcting the z~8 Galaxy Luminosity Function for Gravitational Lensing Magnification Bias

We present a Bayesian framework to account for the magnification bias from both strong and weak gravitational lensing in estimates of high-redshift galaxy luminosity functions. We illustrate our method by estimating the $z\sim8$ UV luminosity function using a sample of 97 Y-band dropouts (Lyman break galaxies) found in the Brightest of Reionizing Galaxies (BoRG) survey and from the literature. We find the luminosity function is well described by a Schechter function with characteristic magnitude of $M^\star = -19.85^{+0.30}_{-0.35}$, faint-end slope of $α= -1.72^{+0.30}_{-0.29}$, and number density of $\log_{10} Ψ^\star [\textrm{Mpc}^{-3}] = -3.00^{+0.23}_{-0.31}$. These parameters are consistent within the uncertainties with those inferred from the same sample without accounting for the magnification bias, demonstrating that the effect is small for current surveys at $z\sim8$, and cannot account for the apparent overdensity of bright galaxies compared to a Schechter function found recently by Bowler et al. (2014a,b) and Finkelstein et al. (2014). We estimate that the probability of finding a strongly lensed $z\sim8$ source in our sample is in the range $\sim 3-15 \%$ depending on limiting magnitude. We identify one strongly-lensed candidate and three cases of intermediate lensing in BoRG (estimated magnification $μ>1.4$) in addition to the previously known candidate group-scale strong lens. Using a range of theoretical luminosity functions we conclude that magnification bias will dominate wide field surveys -- such as those planned for the Euclid and WFIRST missions -- especially at $z>10$. Magnification bias will need to be accounted for in order to derive accurate estimates of high-redshift luminosity functions in these surveys and to distinguish between galaxy formation models.

preprint2015arXiv

Gamma-Ray Bursts Trace UV Metrics of Star Formation over 3 < z < 5

We present the first uniform treatment of long duration gamma-ray burst (GRB) host galaxy detections and upper limits over the redshift range 3<z<5, a key epoch for observational and theoretical efforts to understand the processes, environments, and consequences of early cosmic star formation. We contribute deep imaging observations of 13 GRB positions yielding the discovery of eight new host galaxies. We use this dataset in tandem with previously published observations of 31 further GRB positions to estimate or constrain the host galaxy rest-frame ultraviolet (UV; 1600 A) absolute magnitudes M_UV. We then use the combined set of 44 M_UV estimates and limits to construct the M_UV luminosity function (LF) for GRB host galaxies over 3<z<5 and compare it to expectations from Lyman break galaxy (LBG) photometric surveys with the Hubble Space Telescope. Adopting standard prescriptions for the luminosity dependence of galaxy dust obscuration (and hence, total star formation rate), we find that our LF is compatible with LBG observations over a factor of 600x in host luminosity, from M_UV = -22.5 mag to >-15.6 mag, and with extrapolations of the assumed Schechter-type LF well beyond this range. We review proposed astrophysical and observational biases for our sample, and find they are for the most part minimal. We therefore conclude, as the simplest interpretation of our results, that GRBs successfully trace UV metrics of cosmic star formation over the range 3<z<5. Our findings suggest GRBs are providing an accurate picture of star formation processes from z ~3 out to the highest redshifts.

preprint2015arXiv

Multiple Images of a Highly Magnified Supernova Formed by an Early-Type Cluster Galaxy Lens

In 1964, Refsdal hypothesized that a supernova whose light traversed multiple paths around a strong gravitational lens could be used to measure the rate of cosmic expansion. We report the discovery of such a system. In Hubble Space Telescope imaging, we have found four images of a single supernova forming an Einstein cross configuration around a redshift z=0.54 elliptical galaxy in the MACS J1149.6+2223 cluster. The cluster's gravitational potential also creates multiple images of the z=1.49 spiral supernova host galaxy, and a future appearance of the supernova elsewhere in the cluster field is expected. The magnifications and staggered arrivals of the supernova images probe the cosmic expansion rate, as well as the distribution of matter in the galaxy and cluster lenses.

preprint2015arXiv

Primordial star clusters at extreme magnification

Gravitationally lensed galaxies with magnification ~10-100 are routinely detected at high redshifts, but magnifications significantly higher than this are hampered by a combination of low probability and large source sizes. Magnifications of ~1000 may nonetheless be relevant in the case of intrinsically small, high-redshift objects with very high number densities. Here, we explore the prospects of detecting compact (< 10 pc), high-redshift (z > 7) Population III star clusters at such extreme magnifications in large-area surveys with planned telescopes like Euclid, WFIRST and WISH. We find that the planned WISH 100 sq. deg ultradeep survey may be able to detect a small number of such objects, provided that the total stellar mass of these star clusters is > 10000 solar masses. If candidates for such lensed Population III star clusters are found, follow-up spectroscopy of the surrounding nebula with the James Webb Space Telescope or groundbased Extremely Large Telescopes should be able to confirm the Population III nature of these objects. Multiband photometry of these objects with the James Webb Space Telescope also has the potential to confirm that the stellar initial mass function in these Population III star clusters is top-heavy, as supported by current simulations.

preprint2015arXiv

The environment of bright QSOs at z ~ 6: Star forming galaxies and X-ray emission

We employ cosmological hydrodynamical simulations to investigate models in which the supermassive black holes (BHs) powering luminous z ~ 6 QSOs grow from massive seeds. We simulate at high resolution 18 fields sampling regions with densities ranging from the mean cosmic density all the way to the highest sigma peaks in the Millennium simulation volume. Only in the most massive halos, BHs can grow to masses up to ~ 10^9 Msun by z ~ 6 without invoking super-Eddington accretion. Accretion onto the most massive BHs becomes limited by thermal AGN feedback by z ~ 9-8 with further BH growth proceeding in short Eddington limited bursts. Our modelling suggests that current flux-limited surveys of QSOs at high redshift preferentially detect objects at their peak luminosity and therefore miss a substantial population of QSOs powered by similarly massive BHs but with low accretion rates. To test whether the required host halo masses are consistent with the observed galaxy environment of z ~ 6 QSOs, we produce realistic rest-frame UV images of our simulated galaxies. Without strong stellar feedback, our simulations predict numbers of bright galaxies larger than observed by a factor ten or more. Supernova-driven galactic winds reduce the predicted numbers to a level consistent with observations indicating that stellar feedback was already very efficient at high redshifts. We have further investigated the effect of thermal AGN feedback on the surrounding gas. Our adopted AGN feedback prescription drives mostly energy-driven highly anisotropic outflows with gas speeds of >= 1000 km/s to distances of >= 10 kpc consistent with observations. The spatially extended thermal X-ray emission around bright QSOs powered by these outflows can exceed by large factors the emission expected without AGN feedback and is an important diagnostic of the mechanism whereby AGN feedback energy couples to surrounding gas.

preprint2015arXiv

The Galaxy UV Luminosity Function Before the Epoch of Reionization

We present a model for the evolution of the galaxy ultraviolet (UV) luminosity function (LF) across cosmic time where star formation is linked to the assembly of dark matter halos under the assumption of a mass dependent, but redshift independent, efficiency. We introduce a new self-consistent treatment of the halo star formation history, which allows us to make predictions at $z>10$ (lookback time $\lesssim500$ Myr), when growth is rapid. With a calibration at a single redshift to set the stellar-to-halo mass ratio, and no further degrees of freedom, our model captures the evolution of the UV LF over all available observations ($0\lesssim z\lesssim10$). The significant drop in luminosity density of currently detectable galaxies beyond $z\sim8$ is explained by a shift of star formation toward less massive, fainter galaxies. Assuming that star formation proceeds down to atomic cooling halos, we derive a reionization optical depth $τ= 0.056^{+0.007}_{-0.010}$, fully consistent with the latest Planck measurement, implying that the universe is fully reionized at $z=7.84^{+0.65}_{-0.98}$. In addition, our model naturally produces smoothly rising star formation histories for galaxies with $L\lesssim L_*$ in agreement with observations and hydrodynamical simulations. Before the epoch of reionization at $z>10$ we predict the LF to remain well-described by a Schechter function, but with an increasingly steep faint-end slope ($α\sim-3.5$ at $z\sim16$). Finally, we construct forecasts for surveys with \JWST~and \WFIRST and predict that galaxies out to $z\sim14$ will be observed. Galaxies at $z>15$ will likely be accessible to JWST and WFIRST only through the assistance of strong lensing magnification.

preprint2015arXiv

The Grism Lens-Amplified Survey from Space (GLASS). II. Gas-phase metallicity and radial gradients in an interacting system at z~2

We present spatially resolved gas-phase metallicity for a system of three galaxies at z=1.85 detected in the Grism Lens-Amplified Survey from Space (GLASS). The combination of HST's diffraction limit and strong gravitational lensing by the cluster MACS J0717+3745 results in a spatial resolution of ~200-300 pc, enabling good spatial sampling despite the intrinsically small galaxy sizes. The galaxies in this system are separated by 50-200 kpc in projection and are likely in an early stage of interaction, evidenced by relatively high specific star formation rates. Their gas-phase metallicities are consistent with larger samples at similar redshift, star formation rate, and stellar mass. We obtain a precise measurement of the metallicity gradient for one galaxy and find a shallow slope compared to isolated galaxies at high redshift, consistent with a flattening of the gradient due to gravitational interaction. An alternative explanation for the shallow metallicity gradient and elevated star formation rate is rapid recycling of metal-enriched gas, but we find no evidence for enhanced gas-phase metallicities which should result from this effect. Notably, the measured stellar masses log(M/Msun) = 7.2-9.1 probe to an order of magnitude below previous mass-metallicity studies at this redshift. The lowest mass galaxy has properties similar to those expected for Fornax at this redshift, indicating that GLASS is able to directly study the progenitors of local group dwarf galaxies on spatially resolved scales. Larger samples from the full GLASS survey will be ideal for studying the effects of feedback, and the time evolution of metallicity gradients. These initial results demonstrate the utility of HST spectroscopy combined with gravitational lensing for characterizing resolved physical properties of galaxies at high redshift.

preprint2015arXiv

The Grism Lens-Amplified Survey from Space (GLASS). V. Extent and spatial distribution of star formation in z~0.5 cluster galaxies

We present the first study of the spatial distribution of star formation in z~0.5 cluster galaxies. The analysis is based on data taken with the Wide Field Camera 3 as part of the Grism Lens-Amplified Survey from Space (GLASS). We illustrate the methodology by focusing on two clusters (MACS0717.5+3745 and MACS1423.8+2404) with different morphologies (one relaxed and one merging) and use foreground and background galaxies as field control sample. The cluster+field sample consists of 42 galaxies with stellar masses in the range 10^8-10^11 M_sun, and star formation rates in the range 1-20 M_sun/yr. Both in clusters and in the field, Hα is more extended than the rest-frame UV continuum in 60% of the cases, consistent with diffuse star formation and inside out growth. In ~20% of the cases, the Hα emission appears more extended in cluster galaxies than in the field, pointing perhaps to ionized gas being stripped and/or star formation being enhanced at large radii. The peak of the Hα emission and that of the continuum are offset by less than 1 kpc. We investigate trends with the hot gas density as traced by the X-ray emission, and with the surface mass density as inferred from gravitational lens models and find no conclusive results. The diversity of morphologies and sizes observed in H_alpha illustrates the complexity of the environmental process that regulate star formation. Upcoming analysis of the full GLASS dataset will increase our sample size by almost an order of magnitude, verifying and strengthening the inference from this initial dataset.

preprint2015arXiv

The luminosity and stellar mass functions of GRB host galaxies: Insight into the metallicity bias

Long-Duration Gamma-Ray Bursts (GRBs) are powerful probes of the Universe star formation history, but correlation between the two depends on the highly debated presence/strength of a metallicity bias. To investigate this correlation, we use a phenomenological model that successfully describes star formation rates, luminosities and stellar masses of star forming galaxies, applying it to GRB production. We predict luminosities, stellar masses, and metallicities of host galaxies depending on the metallicity bias. Our best-fitting model includes a moderate metallicity bias, broadly consistent with the large majority of long-duration GRBs in metal-poor environments originating from collapsars (probability ~83%), but with a secondary contribution (~17%) from metal-independent production channels, such as binary evolution. Because of the mass-metallicity relation of galaxies, the maximum likelihood model predicts that the metal-independent channel becomes dominant at z<2, where hosts have higher metallicities and collapsars are suppressed. This possibly explains why some studies find no clear evidence of a metal-bias based on low-z samples. However, while metallicity predictions match observations well at high redshift, there is tension with low redshift observations, since a significant fraction of GRB hosts are predicted to have (near-)solar metallicity. This is in contrast to observations, unless obscured, metal-rich hosts are preferentially missed in current datasets, and suggests that lower efficiencies of the metal-independent GRB channel might be preferred following a comprehensive fit from complete samples. Overall, we are able to establish the presence of a metallicity bias for GRB production, but continued characterization of GRB host galaxies is needed to quantify its strength.

preprint2015arXiv

The relative and absolute ages of old globular clusters in the LCDM framework

Old Globular Clusters (GCs) in the Milky Way have ages of about 13 Gyr, placing their formation time in the reionization epoch. We propose a novel scenario for the formation of these systems based on the merger of two or more atomic cooling halos at high-redshift (z>6). First generation stars are formed as an intense burst in the center of a minihalo that grows above the threshold for hydrogen cooling (halo mass M_h~10^8 Msun) by undergoing a major merger within its cooling timescale (~150 Myr). Subsequent minor mergers and sustained gas infall bring new supply of pristine gas at the halo center, creating conditions that can trigger new episodes of star formation. The dark-matter halo around the GC is then stripped during assembly of the host galaxy halo. Minihalo merging is efficient only in a short redshift window, set by the LCDM parameters, allowing us to make a strong prediction on the age distribution for old GCs. From cosmological simulations we derive an average merging redshift <z>=9 and narrow distribution Dz=2, implying average GC age <t_age>=13.0+/-0.2 Gyr including ~0.2 Gyr of star formation delay. Qualitatively, our scenario reproduces other general old GC properties (characteristic masses and number of objects, metallicity versus galactocentric radius anticorrelation, radial distribution), but unlike age, these generally depend on details of baryonic physics. In addition to improved age measurements, direct validation of the model at z~10 may be within reach of ultradeep gravitationally lensed observations with the James Webb Space Telescope.

preprint2014arXiv

The Luminosity Function at z~8 from 97 Y-band dropouts: Inferences About Reionization

[Abbreviated] We present the largest search to date for $z\sim8$ Lyman break galaxies (LBGs) based on 350 arcmin$^2$ of HST observations in the V-, Y-, J- and H-bands from the Brightest of Reionizing Galaxies (BoRG) survey. The BoRG dataset includes $\sim$50 arcmin$^2$ of new data and deeper observations of two previous BoRG pointings, from which we present 9 new $z\sim8$ LBG candidates, bringing the total number of BoRG LBGs to 38 with $25.5\leqslant m_{J} \leqslant 27.6$ (AB system). We introduce a new Bayesian formalism for estimating the galaxy luminosity function (LF), which does not require binning (and thus smearing) of the data and includes a likelihood based on the formally correct binomial distribution as opposed to the often used approximate Poisson distribution. We demonstrate the utility of the new method on a sample of $97$ LBGs that combines the bright BoRG galaxies with the fainter sources published in Bouwens et al. (2012) from the HUDF and ERS programs. We show that the $z\sim8$ LF is well described by a Schechter function with a characteristic magnitude $M^\star = -20.15^{+0.29}_{-0.38}$, a faint-end slope of $α= -1.87^{+0.26}_{-0.26}$, and a number density of $\log_{10} ϕ^\star [\textrm{Mpc}^{-3}] = -3.24^{+0.25}_{-0.24}$. Integrated down to $M=-17.7$ this LF yields a luminosity density, $\log_{10} ε[\textrm{erg}/\textrm{s/Hz/Mpc}^{3}] = 25.52^{+0.05}_{-0.05}$. Our LF analysis is consistent with previously published determinations within 1$σ$. We discuss the implication of our study for the physics of reionization. By assuming theoretically motivated priors on the clumping factor and the photon escape fraction we show that the UV LF from galaxy samples down to $M=-17.7$ can ionize only 10-50% of the neutral hydrogen at $z\sim8$. Full reionization would require extending the LF down to $M=-15$.

preprint2013arXiv

A physical model for the 0 < z < 8 redshift evolution of the galaxy UV luminosity and stellar mass functions

We present a model to understand the redshift evolution of the UV luminosity and stellar mass functions of Lyman Break Galaxies. Our approach is based on the assumption that the luminosity and stellar mass of a galaxy is related to its dark matter halo assembly and gas infall rate. Specifically, galaxies experience a burst of star formation at the halo assembly time, followed by a constant star formation rate, representing a secular star formation activity sustained by steady gas accretion. Star formation from steady gas accretion is the dominant contribution to the galaxy UV luminosity at all redshifts. The model is calibrated by constructing a galaxy luminosity versus halo mass relation at $z=4$ via abundance matching. After this luminosity calibration, the model naturally fits the $z=4$ stellar mass function, and correctly predicts the evolution of both luminosity and stellar mass functions from $z=0$ to $z=8$. While the details of star formation efficiency and feedback are hidden within our calibrated luminosity versus halo mass relation, our study highlights that the primary driver of galaxy evolution across cosmic time is the build-up of dark matter halos, without the need to invoke a redshift dependent efficiency in converting gas into stars.

preprint2013arXiv

Local Group and Star Cluster Dynamics from HSTPROMO (The Hubble Space Telescope Proper Motion Collaboration)

The Hubble Space Telescope (HST) has proven to be uniquely suited for the measurement of proper motions (PMs) of stars and galaxies in the nearby Universe. Here we summarize the main results and ongoing studies of the HSTPROMO collaboration, which over the past decade has executed some two dozen observational and theoretical HST projects on this topic. This is continuing to revolutionize our dynamical understanding of many objects, including: globular clusters; young star clusters; stars and stellar streams in the Milky Way halo; Local Group galaxies, including dwarf satellite galaxies, the Magellanic Clouds, and the Andromeda galaxy; and AGN Black Hole Jets.

preprint2013arXiv

Tracing the Mass Growth and Star Formation Rate Evolution of Massive Galaxies from z~6 to z~1 in the Hubble Ultra-Deep Field

We present an analysis of $\sim$1500 H160-selected photometric galaxies detected to a limiting magnitude of 27.8 in the HUDF, using imaging from the HST WFC3/IR camera in combination with archival UV, optical, and NIR imaging. We fit photometric redshifts and stellar population estimates for all galaxies with well-determined Spitzer IRAC fluxes, allowing for the determination of the cumulative mass function within the range $1<z<6$. By selecting samples of galaxies at a constant cumulative number density, we explore the co-evolution of stellar masses and star formation rates from z$\sim$6. We find a steady increase in the SFRs of galaxies at constant number density from z$\sim$6 to z$\sim$3. The peak epoch of star formation is found to shift to later times for galaxies with increasing number densities, in agreement with the expectations from cosmic downsizing. The observed SFRs can fully account for the mass growth to z$\sim$2 amongst galaxies with cumulative number densities greater than 10$^{-3.5}$ Mpc$^{-3}$. For galaxies with a lower constant number density we find the observed stellar masses are $\sim$3 times greater than that which may be accounted for by the observed star formation alone at late times, implying that growth from mergers plays an important role at $z<2$. We additionally observe a decreasing sSFR, equivalent to approximately one order of magnitude, from z$\sim$6 to z$\sim$2 amongst galaxies with number densities less than 10$^{-3.5}$ Mpc$^{-3}$ along with significant evidence that at any redshift the sSFR is higher for galaxies at higher number density. The combination of these findings can qualitatively explain the previous findings of a sSFR plateau at high redshift. Tracing the evolution of the fraction of quiescent galaxies for samples matched in cumulative number density over this redshift range, we find no unambiguous examples of quiescent galaxies at $z>4$.

preprint2012arXiv

Constraints on the Ionizing Efficiency of the First Galaxies

Observations of the Lyman-alpha forest and of high-redshift galaxies at z~5-10 imply that there were just enough photons to maintain the universe in an ionized state at z~5-6, indicating a "photon-starved" end to reionization. The ionizing emissivity must have been larger at earlier times in order to yield the extended reionization history implied by the electron scattering optical depth constraint from WMAP. Here we address the possibility that a faint population of galaxies with host halo masses of ~1e8-1e9 Msun dominated the ionizing photon budget at redshifts of about z>9, due to their much higher escape fractions. Such faint, early galaxies, would not have formed in ionized regions due to suppression by heating from the UV background (UVB), and would therefore not contribute to the ionizing background at z<6, after reionization is complete. Our model matches: (1) the low escape fractions observed for high-redshift galaxies, (2) the WMAP constraint of tau_es~0.09, (3) the low values for the UVB at z<6, and (4) the observed star formation rate density inferred from Lyman-break galaxies. A top heavy IMF from Pop III stars is not required in this scenario. We compare our model to recent ones in the literature that were forced to introduce an escape fraction that increases strongly towards high redshift, and show that a similar evolution occurs naturally if low mass galaxies possess high escape fractions.

preprint2012arXiv

Gamma-Ray-Burst Host Galaxy Surveys at Redshift z>4: Probes of Star Formation Rate and Cosmic Reionization

Measuring the star formation rate (SFR) at high redshift is crucial for understanding cosmic reionization and galaxy formation. Two common complementary approaches are Lyman-Break-Galaxy (LBG) surveys for large samples and Gamma-Ray-Burst (GRB) observations for sensitivity to SFR in small galaxies. The z>4 GRB-inferred SFR is higher than the LBG rate, but this difference is difficult to understand, as both methods rely on several modeling assumptions. Using a physically motivated galaxy luminosity function model, with star formation in dark-matter halos with virial temperature Tvir>2e4 K (M_DM>2e8 M_sun), we show that GRB and LBG-derived SFRs are consistent if GRBs extend to faint galaxies (M_AB<-11). To test star formation below the detection limit L_lim~0.05L^*_{z=3} of LBG surveys, we propose to measure the fraction f_det(L>L_lim,z) of GRB hosts with L>L_lim. This fraction quantifies the missing star formation fraction in LBG surveys, constraining the mass-suppression scale for galaxy formation, with weak dependence on modeling assumptions. Because f_det(L>L_lim,z) corresponds to the ratio of star formation rates derived from LBG and GRB surveys, if these estimators are unbiased, measuring f_det(L>L_lim,z) also constrains the redshift evolution of the GRB production rate per unit mass of star formation. Our analysis predicts significant success for GRB host detections at z~5 with f_det(L>L_lim,z)~0.4, but rarer detections at z>6. By analyzing the upper limits on host-galaxy luminosities of six z>5 GRBs from literature data, we infer that galaxies with M_AB>-15 were present at z>5 at 95% confidence, demonstrating the key role played by very faint galaxies during reionization.

preprint2012arXiv

The Evolution of Mass-size Relation for Lyman Break Galaxies From z=1 to z=7

For the first time, we study the evolution of the stellar mass-size relation for star-forming galaxies from z ~ 4 to z ~ 7 from Hubble-WFC3/IR camera observations of the HUDF and Early Release Science (ERS) field. The sizes are measured by determining the best fit model to galaxy images in the rest-frame 2100 Å\ with the stellar masses estimated from SED fitting to rest-frame optical (from Spitzer/IRAC) and UV fluxes. We show that the stellar mass-size relation of Lyman-break galaxies (LBGs) persists, at least to z ~ 5, and the median size of LBGs at a given stellar mass increases towards lower redshifts. For galaxies with stellar masses of 9.5<Log(M*/Msun)<10.4 sizes evolve as $(1+z)^{-1.20\pm0.11}$. This evolution is very similar for galaxies with lower stellar masses of 8.6<Log(M*/Msun)<9.5 which is $r_{e} \propto (1+z)^{-1.18\pm0.10}$, in agreement with simple theoretical galaxy formation models at high z. Our results are consistent with previous measurements of the LBGs mass-size relation at lower redshifts (z ~ 1-3).

preprint2012arXiv

The stellar mass structure of massive galaxies from z=0 to z=2.5; surface density profiles and half-mass radii

We present stellar mass surface density profiles of a mass-selected sample of 177 galaxies at 0.5 < z < 2.5, obtained using very deep HST optical and near-infrared data over the GOODS-South field, including recent CANDELS data. Accurate stellar mass surface density profiles have been measured for the first time for a complete sample of high-redshift galaxies more massive than 10^10.7 M_sun. The key advantage of this study compared to previous work is that the surface brightness profiles are deconvolved for PSF smoothing, allowing accurate measurements of the structure of the galaxies. The surface brightness profiles account for contributions from complex galaxy structures such as rings and faint outer disks. Mass profiles are derived using radial rest-frame u-g color profiles and a well-established empirical relation between these colors and the stellar mass-to-light ratio. We derive stellar half-mass radii from the mass profiles, and find that these are on average ~25% smaller than rest-frame g band half-light radii. This average size difference of 25% is the same at all redshifts, and does not correlate with stellar mass, specific star formation rate, effective surface density, Sersic index, or galaxy size. Although on average the difference between half-mass size and half-light size is modest, for approximately 10% of massive galaxies this difference is more than a factor two. These extreme galaxies are mostly extended, disk-like systems with large central bulges. These results are robust, but could be impacted if the central dust extinction becomes high. ALMA observations can be used to explore this possibility. These results provide added support for galaxy growth scenarios wherein massive galaxies at these epochs grow by accretion onto their outer regions.

preprint2011arXiv

Critical Star-Formation Rates for Reionization: Full Reionization occurs at z = 7

We assess the probable redshift (z_rei ~ 7) for full reionization of the intergalactic medium (IGM) using a prescription for the co-moving star-formation-rate (SFR) density (rho_SFR) required to maintain photoionization against recombination. Our newly developed on-line reionization simulator allows users to assess the required SFR and ionization histories, using a variety of assumptions for galactic and stellar populations, IGM clumping factor and temperature, and LyC escape fraction. The decline in high-redshift galaxy candidates and Lya emitters at z = 6-8 suggests a rising neutral fraction, with reionization at z > 7 increasingly difficult owing to increased recombination rates and constraints from the ionizing background and LyC mean free path. The required rate is rho_SFR = (0.018 M_sun/yr/Mpc^3) [(1+z)/8]^3 (C_H/3)(0.2/f_esc) T_4^(-0.845) scaled to fiducial values of clumping factor C_H = 3, escape fraction f_esc = 0.2, electron temperature T_e = 10^4 K, and low-metallicity initial mass functions (IMF) and stellar atmospheres. Our hydrodynamical + N-body simulations find a mean clumping factor C_H = (2.9)[(1+z)/6]^-1.1 in the photoionized, photoheated filaments at z = 5-9. The critical SFR could be reduced by increasing the minimum stellar mass, invoking a top-heavy IMF, or systematically increasing f_esc at high z. The CMB optical depth, tau_e = 0.088 +/- 0.015, can be explained by full reionization, producing tau_e = 0.050 back to z_rei = 7, augmented by Delta-tau_e = 0.01-0.04 in a UV/X-ray partially ionized IGM at z > 7. In this scenario, the strongest 21-cm signal should occur at redshifted frequencies 124-167 MHz owing to IGM heating over an interval Delta z ~ 3 from z = 7.5-10.5.

preprint2011arXiv

Galaxy Formation in Heavily Overdense Regions at z~10: the Prevalence of Disks in Massive Halos

Using a high-resolution cosmological numerical simulation, we have analyzed the evolution of galaxies at z~10 in a highly overdense region of the universe. These objects could represent the high redshift galaxies recently observed by the Hubble's WFC3, and be as well possible precursors of QSOs at z~6-7. To overcome the sampling and resolution problems in cosmological simulations, we have used the Constrained Realizations method. Our main result for z~10 shows the region of 3.5h^{-1}Mpc radius in comoving coordinates completely dominated by disk galaxies in the total mass range of >=10^9h^{-1}Mo. We have verified that the gaseous and stellar disks we identify are robust morphological features, capable of surviving the ongoing merger process at these redshifts. Below this mass range, we find a sharp decline in the disk fraction to negligible numbers. At this redshift, the disks appear to be gas-rich and the dark matter halos baryon-rich, by a factor of ~2-3 above the average fraction of baryons in the universe. The prevalence of disk galaxies in the high density peaks during the epoch of reionization is contrary to the morphology-density trend observed at low redshifts.

preprint2011arXiv

Morphological evolution of galaxies from ultradeep HST WFC3 imaging: the Hubble sequence at z~2

We use ultradeep HST WFC3/IR imaging of the HUDF to investigate the rest-frame optical morphologies of a mass-selected sample of galaxies at z~2. We find a large variety of galaxy morphologies, ranging from large, blue, disk-like galaxies to compact, red, early-type galaxies. We derive rest-frame u-g color profiles for these galaxies and show that most z~2 galaxies in our sample have negative color gradients such that their cores are red. Although these color gradients may partly be caused by radial variations in dust content, they point to the existence of older stellar populations in the centers of z~2 galaxies. This result is consistent with an "inside-out" scenario of galaxy growth. We find that the median color gradient is fairly constant with redshift: [Delta(u-g_{rest})/Delta(log r)]_{median} = -0.46, -0.44 and -0.49 for z~2, z~1 and z=0, respectively. Using structural parameters derived from surface brightness profiles we confirm that at z~2 galaxy morphology correlates well with specific star formation rate. At the same mass, star forming galaxies have larger effective radii, bluer rest-frame u-g colors and lower Sersic indices than quiescent galaxies. These correlations are very similar to those at lower redshift, suggesting that the relations that give rise to the Hubble sequence at z=0 are already in place for massive galaxies at this early epoch.

preprint2011arXiv

The Dark Side of QSO Formation at High Redshifts

Observed high-redshift QSOs, at z~6, may reside in massive dark matter (DM) halos of more than 10^{12} Msun and are thus expected to be surrounded by overdense regions. In a series of 10 constrained simulations, we have tested the environment of such QSOs. Comparing the computed overdensities with respect to the unconstrained simulations of regions empty of QSOs, assuming there is no bias between the DM and baryon distributions, and invoking an observationally-constrained duty-cycle for Lyman Break Galaxies, we have obtained the galaxy count number for the QSO environment. We find that a clear discrepancy exists between the computed and observed galaxy counts in the Kim et al. (2009) samples. Our simulations predict that on average eight z~6 galaxies per QSO field should have been observed, while Kim et al. detect on average four galaxies per QSO field compared to an average of three galaxies in a control sample (GOODS fields). While we cannot rule out a small number statistics for the observed fields to high confidence, the discrepancy suggests that galaxy formation in the QSO neighborhood proceeds differently than in the field. We also find that QSO halos are the most massive of the simulated volume at z~6 but this is no longer true at z~3. This implies that QSO halos, even in the case they are the most massive ones at high redshifts, do not evolve into most massive galaxy clusters at z=0.

preprint2011arXiv

The UDF05 Follow-up of the Hubble Ultra Deep Field. III. The Luminosity Function at z~6

In this paper, we present a derivation of the rest-frame 1400A luminosity function (LF) at redshift six from a new application of the maximum likelihood method by exploring the five deepest HST/ACS fields, i.e., the HUDF, two UDF05 fields, and two GOODS fields. We work on the latest improved data products, which makes our results more robust than those of previous studies. We use un-binned data and thereby make optimal use of the information contained in the dataset. We focus on the analysis to a magnitude limit where the completeness is larger than 50% to avoid possibly large errors in the faint end slope that are difficult to quantify. We also take into account scattering in and out of the dropout sample due to photometric errors by defining for each object a probability that it belongs to the dropout sample. We find the best fit Schechter parameters to the z~6 LF are: alpha = 1.87 +/- 0.14, M* = -20.25 +/- 0.23, and phi*=1.77^{+0.62}_{-0.49} * 10^{-3} Mpc^{-3}. Such a steep slope suggests that galaxies, especially the faint ones, are possibly the main sources of ionizing photons in the universe at redshift six. We also combine results from all studies at z~6 to reach an agreement in 95% confidence level that -20.45<M*<-20.05 and -1.90<alpha<-1.55. The luminosity density has been found not to evolve significantly between z~6 and z~5, but considerable evolution is detected from z~6 to z~3.

preprint2010arXiv

Confirmation of the compactness of a z=1.91 quiescent galaxy with Hubble Space Telescope's Wide Field Camera 3

We present very deep Wide Field Camera 3 (WFC3) photometry of a massive, compact galaxy located in the Hubble Ultra Deep Field. This quiescent galaxy has a spectroscopic redshift z=1.91 and has been identified as an extremely compact galaxy by Daddi et al. 2005. We use new H-F160W imaging data obtained with Hubble Space Telescope/WFC3 to measure the deconvolved surface brightness profile to H = 28 mag arcsec**-2. We find that the surface brightness profile is well approximated by an n=3.7 Sersic profile. Our deconvolved profile is constructed by a new technique which corrects the best-fit Sersic profile with the residual of the fit to the observed image. This allows for galaxy profiles which deviate from a Sersic profile. We determine the effective radius of this galaxy: r_e=0.42 +- 0.14 kpc in the observed H-F160W-band. We show that this result is robust to deviations from the Sersic model used in the fit. We test the sensitivity of our analysis to faint "wings" in the profile using simulated galaxy images consisting of a bright compact component and a faint extended component. We find that due to the combination of the WFC3 imaging depth and our method's sensitivity to extended faint emission we can accurately trace the intrinsic surface brightness profile, and that we can therefore confidently rule out the existence of a faint extended envelope around the observed galaxy down to our surface brightness limit. These results confirm that the galaxy lies a factor of 10 off from the local mass-size relation.

preprint2010arXiv

Extremely-Metal Poor Stars in the Milky Way: A Second Generation Formed after Reionization

Cosmological simulations of Population III star formation suggest an initial mass function (IMF) biased toward very massive stars (M>100Msun) formed in minihalos at redshift z>20, when the cooling is driven by molecular hydrogen. However, this result conflicts with observations of extremely-metal poor (EMP) stars in the Milky Way halo, whose r-process elemental abundances appear to be incompatible with those expected from very massive Population III progenitors. We propose a new solution to the problem in which the IMF of second-generation stars formed at z>10, before reionization, is deficient in sub-solar mass stars, owing to the high cosmic microwave background temperature floor. The observed EMP stars are formed preferentially at z<10 in pockets of gas enriched to metallicity Z>10^{-3.5} Zsun by winds from Population II stars. Our cosmological simulations of dark matter halos like the Milky Way show that current samples of EMP stars can only constrain the IMF of late-time Population III stars, formed at z<13 in halos with virial temperature Tvir~10^4 K. This suggests that pair instability supernovae were not produced primarily by this population. To begin probing the IMF of Population III stars formed at higher redshift will require a large survey, with at least 500 and probably several thousand EMP stars of metallicities Z~10^{-3.5} Zsun.

preprint2010arXiv

Population III Star Formation During and After the Reionization Epoch

Population III star formation during the dark ages shifted from minihalos (~10^6 Msun) cooled via molecular hydrogen to more massive halos (~10^8 Msun) cooled via Ly-alpha as Lyman-Werner backgrounds progressively quenched molecular hydrogen cooling. Eventually, both modes of primordial star formation were suppressed by the chemical enrichment of the IGM. We present a comprehensive model for following the modes of Population III star formation that is based on a combination of analytical calculations and cosmological simulations. We characterize the properties of the transition from metal-free star formation to the first Population II clusters for an average region of the Universe and for the progenitors of the Milky Way. Finally, we highlight the possibility of observing the explosion of Population III stars within Ly-alpha cooled halos at redshift z~6 in future deep all sky surveys such as LSST.

preprint2010arXiv

Widespread presence of shallow cusps in the surface-brightness profile of globular clusters

Surface brightness profiles of globular clusters with shallow central cusps (Sigma ~ R^v with -0.3<~ v <~ -0.05) have been associated by several recent studies with the presence of a central intermediate mass black hole (IMBH). Such shallow slopes are observed in several globular clusters thanks to the high angular resolution of Hubble Space Telescope imaging. In this Letter we evaluate whether shallow cusps are a unique signature of a central IMBH by analyzing a sample of direct N-body simulations of star clusters with and without a central IMBH. We ``observe'' the simulations as if they were HST images. Shallow cusps are common in our simulation sample: star clusters without an IMBH have v >~ -0.3 in the pre-core-collapse and core-collapse phases. Post-core-collapse clusters without an IMBH transition to steeper cusps, -0.7<~ v <~ -0.4, only if the primordial binary fraction is very small, f_{bin}< 3 per cent, and if there are few stellar-mass black holes remaining. Otherwise v values overlap the range usually ascribed to the presence of an IMBH throughout the entire duration of the simulations. In addition, measuring v is intrinsically prone to significant uncertainty, therefore typical measurement errors may lead to v > -0.3 even when <v> <~ -0.4. Overall our analysis shows that a shallow cusp is not an unequivocal signature of a central IMBH and casts serious doubts on the usefulness of measuring v in the context of the hunt for IMBHs in globular clusters.

preprint2009arXiv

Mass Segregation in NGC 2298: limits on the presence of an Intermediate Mass Black Hole

[abridged] Theoretical investigations have suggested the presence of Intermediate Mass Black Holes (IMBHs, with masses in the 100-10000 Msun range) in the cores of some Globular Clusters (GCs). In this paper we present the first application of a new technique to determine the presence or absence of a central IMBH in globular clusters that have reached energy equipartition via two-body relaxation. The method is based on the measurement of the radial profile for the average mass of stars in the system, using the fact that a quenching of mass segregation is expected when an IMBH is present. Here we measure the radial profile of mass segregation using main-sequence stars for the globular cluster NGC 2298 from resolved source photometry based on HST-ACS data. The observations are compared to expectations from direct N-body simulations of the dynamics of star clusters with and without an IMBH. The mass segregation profile for NGC 2298 is quantitatively matched to that inferred from simulations without a central massive object over all the radial range probed by the observations, that is from the center to about two half-mass radii. Profiles from simulations containing an IMBH more massive than ~ 300-500 Msun (depending on the assumed total mass of NGC 2298) are instead inconsistent with the data at about 3 sigma confidence, irrespective of the IMF and binary fraction chosen for these runs. While providing a null result in the quest of detecting a central black hole in globular clusters, the data-model comparison carried out here demonstrates the feasibility of the method which can also be applied to other globular clusters with resolved photometry in their cores.

preprint2009arXiv

Star Clusters in Pseudo-Bulges of Spiral Galaxies

We present a study of the properties of the star-cluster systems around pseudo-bulges of late-type spiral galaxies using a sample of 11 galaxies with distances from 17 to 37 Mpc. Star clusters are identified from multiband HST ACS and WFPC2 imaging data by combining detections in 3 bands (F435W and F814W with ACS and F606W with WFPC2). The photometric data are then compared to population synthesis models to infer the masses and ages of the star clusters. Photometric errors and completeness are estimated by means of artificial source Monte Carlo simulations. Dust extinction is estimated by considering F160W NICMOS observations of the central regions of the galaxies, augmenting our wavelength coverage. In all galaxies we identify star clusters with a wide range of ages, from young (age < 8 Myr) blue clusters, with typical mass of 10^3 Msun to older (age > 100-250 Myr), more massive, red clusters. Some of the latter might likely evolve into objects similar to the Milky Way's globular clusters. We compute the specific frequencies for the older clusters with respect to the galaxy and bulge luminosities. Specific frequencies relative to the galaxy light appear consistent with the globular cluster specific frequencies of early-type spirals. We compare the specific frequencies relative to the bulge light with the globular cluster specific frequencies of dwarf galaxies, which have a surface-brightness profile that is similar to that of the pseudo-bulges in our sample. The specific frequencies we derive for our sample galaxies are higher than those of the dwarf galaxies, supporting an evolutionary scenario in which some of the dwarf galaxies might be the remnants of harassed late-type spiral galaxies which hosted a pseudo-bulge.

preprint2009arXiv

Tidal disruption, global mass function and structural parameters evolution in star clusters

[abridged] We present a unified picture for the evolution of star clusters on the two-body relaxation timescale. We use direct N-body simulations of star clusters in a galactic tidal field starting from different multi-mass King models, up to 10% of primordial binaries and up to Ntot=65536 particles. An additional run also includes a central Intermediate Mass Black Hole. We find that for the broad range of initial conditions we have studied the stellar mass function of these systems presents a universal evolution which depends only on the fractional mass loss. The structure of the system, as measured by the core to half mass radius ratio, also evolves toward a universal state, which is set by the efficiency of heating on the visible population of stars induced by dynamical interactions in the core of the system. Interactions with dark remnants are dominant over the heating induced by a moderate population of primordial binaries (3-5%), especially under the assumption that most of the neutron stars and black holes are retained in the system. All our models without primordial binaries undergo a deep gravothermal collapse in the radial mass profile. However their projected light distribution can be well fitted by medium concentration King models (with parameter W0 ~ 8), even though there tends to be an excess over the best fit for the innermost points of the surface brightness. This excess is consistent with a shallow cusp in the surface brightness (mu(R) ~ R^{-v} with v ~ 0.4-0.7), like it has been observed for many globular clusters from high-resolution HST imaging. Classification of core-collapsed globular clusters based on their surface brightness profile is likely to fail in systems that have already bounced back to lower concentrations.