Source author record

Gerard Lemson

Gerard Lemson 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

15works
7topics
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

15 published item(s)

preprint2020arXiv

SciServer: a Science Platform for Astronomy and Beyond

We present SciServer, a science platform built and supported by the Institute for Data Intensive Engineering and Science at the Johns Hopkins University. SciServer builds upon and extends the SkyServer system of server-side tools that introduced the astronomical community to SQL (Structured Query Language) and has been serving the Sloan Digital Sky Survey catalog data to the public. SciServer uses a Docker/VM based architecture to provide interactive and batch mode server-side analysis with scripting languages like Python and R in various environments including Jupyter (notebooks), RStudio and command-line in addition to traditional SQL-based data analysis. Users have access to private file storage as well as personal SQL database space. A flexible resource access control system allows users to share their resources with collaborators, a feature that has also been very useful in classroom environments. All these services, wrapped in a layer of REST APIs, constitute a scalable collaborative data-driven science platform that is attractive to science disciplines beyond astronomy.

preprint2016arXiv

Galaxy formation in the Planck cosmology II. Star formation histories and post-processing magnitude reconstruction

We adapt the L-Galaxies semi-analytic model to follow the star-formation histories (SFH) of galaxies -- by which we mean a record of the formation time and metallicities of the stars that are present in each galaxy at a given time. We use these to construct stellar spectra in post-processing, which offers large efficiency savings and allows user-defined spectral bands and dust models to be applied to data stored in the Millennium data repository. We contrast model SFHs from the Millennium Simulation with observed ones from the VESPA algorithm as applied to the SDSS-7 catalogue. The overall agreement is good, with both simulated and SDSS galaxies showing a steeper SFH with increased stellar mass. The SFHs of blue and red galaxies, however, show poor agreement between data and simulations, which may indicate that the termination of star formation is too abrupt in the models. The mean star-formation rate (SFR) of model galaxies is well-defined and is accurately modelled by a double power law at all redshifts: SFR proportional to $1/(x^{-1.39}+x^{1.33})$, where $x=(t_a-t)/3.0\,$Gyr, $t$ is the age of the stars and $t_a$ is the loopback time to the onset of galaxy formation; above a redshift of unity, this is well approximated by a gamma function: SFR proportional to $x^{1.5}e^{-x}$, where $x=(t_a-t)/2.0\,$Gyr. Individual galaxies, however, show a wide dispersion about this mean. When split by mass, the SFR peaks earlier for high-mass galaxies than for lower-mass ones, and we interpret this downsizing as a mass-dependence in the evolution of the quenched fraction: the SFHs of star-forming galaxies show only a weak mass dependence.

preprint2016arXiv

The EAGLE simulations of galaxy formation: public release of halo and galaxy catalogues

We present the public data release of halo and galaxy catalogues extracted from the EAGLE suite of cosmological hydrodynamical simulations of galaxy formation. These simulations were performed with an enhanced version of the GADGET code that includes a modified hydrodynamics solver, time-step limiter and subgrid treatments of baryonic physics, such as stellar mass loss, element-by-element radiative cooling, star formation and feedback from star formation and black hole accretion. The simulation suite includes runs performed in volumes ranging from 25 to 100 comoving megaparsecs per side, with numerical resolution chosen to marginally resolve the Jeans mass of the gas at the star formation threshold. The free parameters of the subgrid models for feedback are calibrated to the redshift z=0 galaxy stellar mass function, galaxy sizes and black hole mass - stellar mass relation. The simulations have been shown to match a wide range of observations for present-day and higher-redshift galaxies. The raw particle data have been used to link galaxies across redshifts by creating merger trees. The indexing of the tree produces a simple way to connect a galaxy at one redshift to its progenitors at higher redshift and to identify its descendants at lower redshift. In this paper we present a relational database which we are making available for general use. A large number of properties of haloes and galaxies and their merger trees are stored in the database, including stellar masses, star formation rates, metallicities, photometric measurements and mock gri images. Complex queries can be created to explore the evolution of more than 10^5 galaxies, examples of which are provided in appendix. (abridged)

preprint2015arXiv

Galaxy formation in the Planck Cosmology - I. Matching the observed evolution of star formation rates, colours and stellar masses

We have updated the Munich galaxy formation model to the Planck first-year cosmology, while modifying the treatment of baryonic processes to reproduce recent data on the abundance and passive fractions of galaxies from z= 3 down to z=0. Matching these more extensive and more precise observational results requires us to delay the reincorporation of wind ejecta, to lower the surface density threshold for turning cold gas into stars, to eliminate ram-pressure stripping in haloes less massive than ~10^14 Msun, and to modify our model for radio mode feedback. These changes cure the most obvious failings of our previous models, namely the overly early formation of low-mass galaxies and the overly large fraction of them that are passive at late times. The new model is calibrated to reproduce the observed evolution both of the stellar mass function and of the distribution of star formation rate at each stellar mass. Massive galaxies (M>10^11 [Msun]) assemble most of their mass before z=1 and are predominantly old and passive at z=0, while lower mass galaxies assemble later and, for M<10^9.5 (Msun), are still predominantly blue and star forming at z=0. This phenomenological but physically based model allows the observations to be interpreted in terms of the efficiency of the various processes that control the formation and evolution of galaxies as a function of their stellar mass, gas content, environment and time.

preprint2013arXiv

Galaxy formation in WMAP1 and WMAP7 cosmologies

Using the technique of Angulo & White (2010) we scale the Millennium and Millennium-II simulations of structure growth in a LCDM universe from the cosmological parameters with which they were carried out (based on first-year results from the Wilkinson Microwave Anisotropy Probe, WMAP1) to parameters consistent with the seven-year WMAP data (WMAP7). We implement semi-analytic galaxy formation modelling on both simulations in both cosmologies to investigate how the formation, evolution and clustering of galaxies are predicted to vary with cosmological parameters. The increased matter density Omega_m and decreased linear fluctuation amplitude sigma8 in WMAP7 have compensating effects, so that the abundance and clustering of dark halos are predicted to be very similar to those in WMAP1 for z <= 3. As a result, local galaxy properties can be reproduced equally well in the two cosmologies by slightly altering galaxy formation parameters. The evolution of the galaxy populations is then also similar. In WMAP7, structure forms slightly later. This shifts the peak in cosmic star formation rate to lower redshift, resulting in slightly bluer galaxies at z=0. Nevertheless, the model still predicts more passive low-mass galaxies than are observed. For rp< 1Mpc, the z=0 clustering of low-mass galaxies is weaker for WMAP7 than for WMAP1 and closer to that observed, but the two cosmologies give very similar results for more massive galaxies and on large scales. At z>1 galaxies are predicted to be more strongly clustered for WMAP7. Differences in galaxy properties, including, clustering, in these two cosmologies are rather small up to redshift 3. Given that there are still considerable residual uncertainties in galaxy formation models, it is very difficult to distinguish WMAP1 from WMAP7 through observations of galaxy properties or their evolution.

preprint2013arXiv

On the Spin Bias of Satellite Galaxies in the Local Group-like Environment

We utilize the Millennium-II simulation databases to study the spin bias of dark subhalos in the Local Group-like systems which have two prominent satellites with comparable masses. Selecting the group-size halos with total mass similar to that of the Local Group (LG) from the friends-of-friends halo catalog and locating their subhalos from the substructure catalog, we determine the most massive (main) and second to the most massive (submain) ones among the subhalos hosted by each selected halo. When the dimensionless spin parameter (lambda) of each subhalo is derived from its specific angular momentum and circular velocity at virial radius, a signal of correlation is detected between the spin parameters of the subhalos and the main-to-submain mass ratios of their host halos at z=0: The higher main-to-submain mass ratio a host halo has, the higher mean spin parameter its subhalos have. It is also found that the correlations exist even for the subhalo progenitors at z=0.5 and z=1. Our interpretation of this result is that the subhalo spin bias is not a transient effect but an intrinsic property of a LG-like system with higher main-to- submain mass ratio, caused by stronger anisotropic stress in the region. A cosmological implication of our result is also discussed.

preprint2012arXiv

Confronting theoretical models with the observed evolution of the galaxy population out to z=4

[abridged] We construct lightcones for the semi-analytic galaxy formation simulation of Guo et al. (2011) and make mock catalogues for comparison with deep high-redshift surveys. Photometric properties are calculated with two different stellar population synthesis codes (Bruzual & Charlot 2003; Maraston 2005) in order to study sensitivity to this aspect of the modelling. The catalogues are publicly available and include photometry for a large number of observed bands from 4000°A to 6μm, as well as rest-frame photometry and intrinsic properties of the galaxies. Guo et al. (2011) tuned their model to fit the low-redshift galaxy population but noted that at z > 1 it overpredicts the abundance of galaxies below the "knee" of the stellar mass function. Here we extend the comparison to deep galaxy counts in the B, i, J, K and IRAC 3.6μm, 4.5μm and 5.8μm bands, to the redshift distributions of K and 5.8μm selected galaxies, and to the evolution of rest-frame luminosity functions in the B and K bands. The B, i and J counts are well reproduced, but at longer wavelengths the overabundant high-redshift galaxies produce excess faint counts. The predicted redshift distributions for K and 5.8μm selected samples highlight the effect of emission from thermally pulsing AGB stars. The full treatment of Maraston (2005) predicts three times as many z~2 galaxies in faint 5.8μm selected samples as the model of Bruzual & Charlot (2003), whereas the two models give similar predictions for K-band selected samples. Although luminosity functions are adequately reproduced out to z~3 in rest-frame B, the same is true at rest-frame K only if TP-AGB emission is included, and then only at high luminosity. Fainter than L* the two synthesis models agree but overpredict the number of galaxies, another reflection of the overabundance of ~10^10M\odot model galaxies at z > 1.

preprint2011arXiv

From dwarf spheroidals to cDs: Simulating the galaxy population in a LCDM cosmology

We apply updated semi-analytic galaxy formation models simultaneously to the stored halo/subhalo merger trees of the Millennium and Millennium-II simulations. These differ by a factor of 125 in mass resolution, allowing explicit testing of resolution effects on predicted galaxy properties. We have revised the treatments of the transition between the rapid infall and cooling flow regimes of gas accretion, of the sizes of bulges and of gaseous and stellar disks, of supernova feedback, of the transition between central and satellite status as galaxies fall into larger systems, and of gas and star stripping once they become satellites. Plausible values of efficiency and scaling parameters yield an excellent fit not only to the observed abundance of low-redshift galaxies over 5 orders of magnitude in stellar mass and 9 magnitudes in luminosity, but also to the observed abundance of Milky Way satellites. This suggests that reionisation effects may not be needed to solve the "missing satellite" problem except, perhaps, for the faintest objects. The same model matches the observed large-scale clustering of galaxies as a function of stellar mass and colour. The fit remains excellent down to ~30kpc for massive galaxies. For M* < 6 x 10^10Msun, however, the model overpredicts clustering at scales below 1 Mpc, suggesting that the sigma_8 adopted in the simulations (0.9) is too high. Galaxy distributions within rich clusters agree between the simulations and match those observed, but only if galaxies without dark matter subhalos (so-called orphans) are included. Our model predicts a larger passive fraction among low-mass galaxies than is observed, as well as an overabundance of ~10^10Msun galaxies beyond z~0.6, reflecting deficiencies in the way star-formation rates are modelled.

preprint2011arXiv

Local Gravity versus Local Velocity: Solutions for $β$ and nonlinear bias

(abridged) We perform a reconstruction of the cosmological large scale flows in the nearby Universe using two complementary observational sets. The first, the SFI++ sample of Tully-Fisher (TF) measurements of galaxies, provides a direct probe of the flows. The second, the whole sky distribution of galaxies in the 2MASS redshift survey (2MRS), yields a prediction of the flows given the cosmological density parameter, $Ω$, and a biasing relation between mass and galaxies. We aim at an unbiased comparison between the peculiar velocity fields extracted from the two data sets and its implication on the cosmological parameters and the biasing relation. We expand the fields in a set of orthonormal basis functions, each representing a plausible realization of a cosmological velocity field. Our analysis completely avoids the strong error covariance in the smoothed TF velocities by the use of orthonormal basis functions and employs elaborate realistic mock data sets to extensively calibrate the errors in 2MRS predicted velocities. We relate the 2MRS galaxy distribution to the mass density field by a linear bias factor, $b$, and include a luminosity dependent, $\propto L^α$, galaxy weighting. We assess the agreement between the fields as a function of $α$ and $β=f(Ω)/b$, where $f$ is the growth factor of linear perturbations. The agreement is excellent with a reasonable $χ^2$ per degree of freedom. For $α=0$, we derive $0.28<β<0.37$ and $0.24<β<0.43$, respectively, at the 68.3% and 95.4% confidence levels (CLs). For $β=0.33$, we get $α<0.25$ and $α<0.5$, respectively, at the 68.3% and 95.4% CLs. We set a constraint on the fluctuation normalization, finding $σ_8 = 0.73 \pm 0.1$, in very good agreement with the latest WMAP results.

preprint2011arXiv

Observing simulated galaxy clusters with PHOX: a novel X-ray photon simulator

We present a novel, virtual X-ray observatory designed to obtain synthetic observations from hydro-numerical simulations, named PHOX. In particular, we provide a description of the code constituting the photon simulator and of the new approach implemented. We apply PHOX to simulated galaxy clusters in order to demonstrate its capabilities. In fact, X-ray observations of clusters of galaxies continue to provide us with an increasingly detailed picture of their structure and of the underlying physical phenomena governing the gaseous component, which dominates their baryonic content. Therefore, it is fundamental to find the most direct and faithful way to compare such observational data with hydrodynamical simulations of cluster-like objects, which can currently include various complex physical processes. Here, we present and analyse synthetic Suzaku observations of two cluster-size haloes obtained by processing with PHOX the hydrodynamical simulation of the large-scale, filament-like region in which they reside. Taking advantage of the simulated data, we test the results inferred from the X-ray analysis of the mock observations against the underlying, known solution. Remarkably, we are able to recover the theoretical temperature distribution of the two haloes by means of the multi-temperature fitting of the synthetic spectra. Moreover, the shapes of the reconstructed distributions allow us to trace the different thermal structure that distinguishes the dynamical state of the two haloes.

preprint2011arXiv

The MultiDark Database: Release of the Bolshoi and MultiDark Cosmological Simulations

We present the online MultiDark Database -- a Virtual Observatory-oriented, relational database for hosting various cosmological simulations. The data is accessible via an SQL (Structured Query Language) query interface, which also allows users to directly pose scientific questions, as shown in a number of examples in this paper. Further examples for the usage of the database are given in its extensive online documentation (www.multidark.org). The database is based on the same technology as the Millennium Database, a fact that will greatly facilitate the usage of both suites of cosmological simulations. The first release of the MultiDark Database hosts two 8.6 billion particle cosmological N-body simulations: the Bolshoi (250/h Mpc simulation box, 1/h kpc resolution) and MultiDark Run1 simulation (MDR1, or BigBolshoi, 1000/h Mpc simulation box, 7/h kpc resolution). The extraction methods for halos/subhalos from the raw simulation data, and how this data is structured in the database are explained in this paper. With the first data release, users get full access to halo/subhalo catalogs, various profiles of the halos at redshifts z=0-15, and raw dark matter data for one time-step of the Bolshoi and four time-steps of the MultiDark simulation. Later releases will also include galaxy mock catalogs and additional merging trees for both simulations as well as new large volume simulations with high resolution. This project is further proof of the viability to store and present complex data using relational database technology. We encourage other simulators to publish their results in a similar manner.

preprint2009arXiv

Cosmic Cartography of the Large-Scale Structure with Sloan Digital Sky Survey Data Release 6

We present the largest Wiener reconstruction of the cosmic density field made to date. The reconstruction is based on the Sloan Digital Sky Survey data release 6 covering the northern Galactic cap. We use a novel supersampling algorithm to suppress aliasing effects and a Krylov-space inversion method to enable high performance with high resolution. These techniques are implemented in the ARGO computer code. We reconstruct the field over a 500 Mpc cube with Mpc grid-resolution while accounting both for the angular and radial selection functions of the SDSS, and the shot noise giving an effective resolution of the order of ~10 Mpc. In addition, we correct for the redshift distortions in the linear and nonlinear regimes in an approximate way. We show that the commonly used method of inverse weighting the galaxies by the corresponding selection function heads to excess noise in regions where the density of the observed galaxies is small. It is more accurate and conservative to adopt a Bayesian framework in which we model the galaxy selection/detection process to be Poisson-binomial. This results in heavier smoothing in regions of reduced sampling density. Our results show a complex cosmic web structure with huge void regions indicating that the recovered matter distribution is highly non-Gaussian. Filamentary structures are clearly visible on scales up to ~20 Mpc. We also calculate the statistical distribution of density after smoothing the reconstruction with Gaussian kernels of different radii r_S and find good agreement with a log-normal distribution for ~10 Mpc < r_S < ~30 Mpc.

preprint2008arXiv

Quantifying the Cosmic Web I: The large-scale halo ellipticity-ellipticity and ellipticity-direction correlations

The formation of dark matter halos tends to occur anisotropically along the filaments of the Cosmic Web, which induces both ellipticity-ellipticity (EE) correlations between the shapes of halos, as well as ellipticity-direction (ED) cross-correlations between halo shapes and the directions to neighboring halos. We analyze the halo catalogue and the semi-analytic galaxy catalogue of the recent Millennium Run Simulation to measure the EE and ED correlations numerically at four different redshifts (z=0, 0.5, 1 and 2). For the EE correlations, we find that (i) the major-axis correlation is strongest while the intermediate-axis correlation is weakest; (ii) the signal is significant at distances out to 10 Mpc/h; (iii) the signal decreases as z decreases; (iv) and its behavior depends strongly on the halo mass scale, with larger masses showing stronger correlations at large distances. For the ED correlations, we find that (i) the correlations are much stronger than the EE correlations, and are significant even out to distances of 50 Mpc/h; (ii) the signal also decreases as z decreases; (iii) and it increases with halo mass at all distances. We also provide empirical fitting functions for the EE and ED correlations. The EE correlations are found to scale linearly with the linear density correlation function, xi(r). While the ED cross-correlation is found to scale as xi^{1/2}(r) at large distances beyond 10 Mpc/h. The best-fit values of the fitting parameters for the EE and the ED correlations are all determined through chi^{2}-statistics. Our results may be useful for quantifying the filamentary distribution of dark matter halos over a wide range of scales.

preprint2008arXiv

Quenching of Star Formation

In the last decade we have seen an enormous increase in the size and quality of spectroscopic galaxy surveys, both at low and high redshift. New statistical techniques to analyse large portions of galaxy spectra are now finding favour over traditional index based methods. Here we will review a new robust and iterative Principal Component Analysis (PCA) algorithm, which solves several common issues with classic PCA. Application to the 4000AA break region of galaxies in the VIMOS VLT Deep Survey (VVDS) and Sloan Digital Sky Survey (SDSS) gives new high signal-to-noise ratio spectral indices easily interpretable in terms of recent star formation history. In particular, we identify a sample of post-starburst galaxies at z~0.7 and z~0.07. We quantify for the first time the importance of post-starburst galaxies, consistent with being descendants of gas-rich major mergers, for building the red sequence. Finally, we present a comparison with new low and high redshift "mock spectroscopic surveys" derived from a Millennium Run semi-analytic model.

preprint2007arXiv

Bursty stellar populations and obscured AGN in galaxy bulges

[Abridged] We investigate trends between the recent star formation history and black hole growth in galaxy bulges in the Sloan Digital Sky Survey (SDSS). The galaxies lie at 0.01<z<0.07 where the fibre aperture covers only the central 0.6-4.0kpc diameter of the galaxy. We find strong trends between black hole growth, as measured by dust-attenuation-corrected OIII luminosity, and the recent star formation history of the bulges. We conclude that our results support the popular hypothesis for black hole growth occurring through gas inflow into the central regions of galaxies, followed by a starburst and triggering of the AGN. However, while this is a significant pathway for the growth of black holes, it is not the dominant one in the present-day Universe. More unspectacular processes are apparently responsible for the majority of this growth. In order to arrive at these conclusions we have developed a set of new high signal-to-noise ratio (SNR) optical spectral indicators, designed to allow a detailed study of stellar populations which have undergone recent enhanced star formation. Working in the rest-frame wavelength range 3750-4150AA, ideally suited to many recent and ongoing spectroscopic surveys at low and high redshift, the first two indices are equivalent to the previously well studied 4000AA break strength and Hdelta equivalent width. The primary advantage of this new method is a greatly improved SNR for the latter index, allowing the present study to use spectra with SNR-per-pixel as low as 8.