Source author record

Benjamin P. Moster

Benjamin P. Moster 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

21works
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

21 published item(s)

preprint2022arXiv

Observational Measures of Halo Properties Beyond Mass

Different properties of dark matter haloes, including growth rate, concentration, interaction history, and spin, correlate with environment in unique, scale-dependent ways. While these halo properties are not directly observable, galaxies will inherit their host haloes' correlations with environment. In this paper, we show how these characteristic environmental signatures allow using measurements of galaxy environment to constrain which dark matter halo properties are most tightly connected to observable galaxy properties. We show that different halo properties beyond mass imprint distinct scale-dependent signatures in both the galaxy two-point correlation function and the distribution of distances to galaxies' $k$th nearest neighbours, with features strong enough to be accessible even with low-resolution (e.g., grism) spectroscopy at higher redshifts. As an application, we compute observed two-point correlation functions for galaxies binned by half-mass radius at $z=0$ from the Sloan Digital Sky Survey, showing that classic galaxy size models (i.e., galaxy size being proportional to halo spin) as well as other recent proposals show significant tensions with observational data. We show that the agreement with observed clustering can be improved with a simple empirical model in which galaxy size correlates with halo growth.

preprint2021arXiv

EMERGE: Empirical predictions of galaxy merger rates since $z\sim6$

We explore the galaxy-galaxy merger rate with the empirical model for galaxy formation, Emerge. On average, we find that between $2$ per cent and $20$ per cent of massive galaxies ($\log_{10}(m_{*}/M_{\odot}) \geq 10.3$) will experience a major merger per Gyr. Our model predicts galaxy merger rates that do not scale as a power-law with redshift when selected by descendant stellar mass, and exhibit a clear stellar mass and mass-ratio dependence. Specifically, major mergers are more frequent at high masses and at low redshift. We show mergers are significant for the stellar mass growth of galaxies $\log_{10}(m_{*}/M_{\odot}) \gtrsim 11.0$. For the most massive galaxies major mergers dominate the accreted mass fraction, contributing as much as $90$ per cent of the total accreted stellar mass. We reinforce that these phenomena are a direct result of the stellar-to-halo mass relation, which results in massive galaxies having a higher likelihood of experiencing major mergers than low mass galaxies. Our model produces a galaxy pair fraction consistent with recent observations, exhibiting a form best described by a power-law exponential function. Translating these pair fractions into merger rates results in an inaccurate prediction compared to the model intrinsic values when using published observation timescales. We find the pair fraction can be well mapped to the intrinsic merger rate by adopting an observation timescale that decreases linearly with redshift as $T_{\mathrm{obs}} = -0.36(1+z)+2.39$ [Gyr], assuming all observed pairs merge by $z=0$.

preprint2021arXiv

Exoplanet Detection using Machine Learning

We introduce a new machine learning based technique to detect exoplanets using the transit method. Machine learning and deep learning techniques have proven to be broadly applicable in various scientific research areas. We aim to exploit some of these methods to improve the conventional algorithm based approaches presently used in astrophysics to detect exoplanets. Using the time-series analysis library TSFresh to analyse light curves, we extracted 789 features from each curve, which capture the information about the characteristics of a light curve. We then used these features to train a gradient boosting classifier using the machine learning tool lightgbm. This approach was tested on simulated data, which showed that is more effective than the conventional box least squares fitting (BLS) method. We further found that our method produced comparable results to existing state-of-the-art deep learning models, while being much more computationally efficient and without needing folded and secondary views of the light curves. For Kepler data, the method is able to predict a planet with an AUC of 0.948, so that 94.8 per cent of the true planet signals are ranked higher than non-planet signals. The resulting recall is 0.96, so that 96 per cent of real planets are classified as planets. For the Transiting Exoplanet Survey Satellite (TESS) data, we found our method can classify light curves with an accuracy of 0.98, and is able to identify planets with a recall of 0.82 at a precision of 0.63.

preprint2020arXiv

Hot phase generation by supernovae: resolution, chemistry and thermal conduction

Supernovae (SN) generate hot gas in the interstellar medium (ISM), help setting the ISM structure and support the driving of outflows. It is important to resolve the hot gas generation for galaxy formation simulations at solar mass and sub-parsec resolution which realise individual supernova (SN) explosions with ambient densities varying by several orders of magnitude in a realistic multi-phase ISM. We test resolution requirements by simulating SN blast waves at three metallicities ($Z = 0.01, 0.1$ and $1 Z_{\odot}$), six densities and their respective equilibrium chemical compositions ($n=0.001$ cm$^{-3}$ - $100$ cm$^{-3}$), and four mass resolutions ($0.1$ - $100$ M$_{\odot}$), in three dimensions. We include non-equilibrium cooling and chemistry, a homogenous interstellar radiation field, and shielding with a modern pressure-energy smoothed particle hydrodynamics (SPH) method including isotropic thermal conduction and a meshless-finite-mass (MFM) solver. We find stronger resolution requirements for chemistry and hot phase generation than for momentum generation. While at $10$ M$_{\odot}$ the radial momenta at the end of the Sedov phase start converging, the hot phase generation and chemistry require higher resolutions to represent the neutral to ionised hydrogen fraction at the end of the Sedov phase correctly. Thermal conduction typically reduces the hot phase by $0.2$ dex and has little impact on the chemical composition. In general, our $1$, and $0.1$ M$_{\odot}$ results agree well with previous numerical and analytic estimates. We conclude that for the thermal energy injection SN model presented here resolutions higher than $10$ M$_{\odot}$ are required to model the chemistry, momentum and hot phase generation in a multi-phase ISM.

preprint2020arXiv

The GRIFFIN project -- Formation of star clusters with individual massive stars in a simulated dwarf galaxy starburst

We describe a population of young star clusters (SCs) formed in a hydrodynamical simulation of a gas-rich dwarf galaxy merger resolved with individual massive stars at sub-parsec spatial resolution. The simulation is part of the GRIFFIN (Galaxy Realizations Including Feedback From INdividual massive stars) project. The star formation environment during the simulation spans seven orders of magnitude in gas surface density and thermal pressure, and the global star formation rate surface density ($Σ_\mathrm{SFR}$) varies by more than three orders of magnitude during the simulation. Young SCs more massive than $M_{\mathrm{*,cl}}\sim 10^{2.5}\,M_{\odot}$ form along a mass function with a power-law index $α\sim-1.7$ ($α\sim-2$ for $M_{\mathrm{*,cl}}\gtrsim10^{3}\,M_{\odot}$) at all merger phases, while the normalization and the highest SC masses (up to $\sim 10^6 M_{\odot}$) correlate with $Σ_\mathrm{SFR}$. The cluster formation efficiency varies from $Γ\sim20\%$ in early merger phases to $Γ\sim80\%$ at the peak of the starburst and is compared to observations and model predictions. The massive SCs ($\gtrsim10^4\,M_{\odot}$) have sizes and mean surface densities similar to observed young massive SCs. Simulated lower mass clusters appear slightly more concentrated than observed. All SCs form on timescales of a few Myr and lose their gas rapidly resulting in typical stellar age spreads between $σ\sim0.1-2$ Myr ($1σ$), consistent with observations. The age spreads increase with cluster mass, with the most massive cluster ($\sim10^6\, M_{\odot}$) reaching a spread of $5\, \mathrm{Myr}$ once its hierarchical formation finishes. Our study shows that it is now feasible to investigate the SC population of entire galaxies with novel high-resolution numerical simulations.

preprint2019arXiv

On the origin of magnetic driven winds and the structure of the galactic dynamo in isolated galaxies

We investigate the build-up of the galactic dynamo and subsequently the origin of a magnetic driven outflow. We use a setup of an isolated disc galaxy with a realistic circum-galactic medium (CGM). We find good agreement of the galactic dynamo with theoretical and observational predictions from the radial and toroidal components of the magnetic field as function of radius and disc scale height. We find several field reversals indicating dipole structure at early times and quadrupole structure at late times. Together with the magnetic pitch angle and the dynamo control parameters $R_α$, $R_ω$ and $D$ we present strong evidence for an $α^2$-$Ω$ dynamo. The formation of a bar in the centre leads to further amplification of the magnetic field via adiabatic compression which subsequently drives an outflow. Due to the Parker-Instability the magnetic field lines rise to the edge of the disc, break out and expand freely in the CGM driven by the magnetic pressure. Finally, we investigate the correlation between magnetic field and star formation rate. Globally, we find that the magnetic field is increasing as function of the star formation rate surface density with a slope between $0.3$ and $0.45$ in good agreement with predictions from theory and observations. Locally, we find that the magnetic field can decrease while star formation increases. We find that this effect is correlated with the diffusion of magnetic field from the spiral arms to the inter-arm regions which we explicitly include by solving the induction equation and accounting for non-linear terms.

preprint2015arXiv

A New Channel for the Formation of Kinematically Decoupled Cores in Early-type galaxies

We present the formation of a Kinematically Decoupled Core (KDC) in an elliptical galaxy, resulting from a major merger simulation of two disk galaxies. We show that although the two progenitor galaxies are initially following a prograde orbit, strong reactive forces during the merger can cause a short-lived change of their orbital spin; the two progenitors follow a retrograde orbit right before their final coalescence. This results in a central kinematic decoupling and the formation of a large-scale (~2 kpc radius) counter-rotating core (CRC) at the center of the final elliptical-like merger remnant (M*=1.3x10^11 Msun), while its outer parts keep the rotation direction of the initial orbital spin. The stellar velocity dispersion distribution of the merger remnant galaxy exhibits two symmetrical off-centered peaks, comparable to the observed "2-sigma galaxies". The KDC/CRC consists mainly of old, pre-merger population stars (older than 5 Gyr), remaining prominent in the center of the galaxy for more than 2 Gyr after the coalescence of its progenitors. Its properties are consistent with KDCs observed in massive elliptical galaxies. This new channel for the formation of KDCs from prograde mergers is in addition to previously known formation scenarios from retrograde mergers and can help towards explaining the substantial fraction of KDCs observed in early-type galaxies.

preprint2015arXiv

From Discs to Bulges: effect of mergers on the morphology of galaxies

We study the effect of mergers on the morphology of galaxies by means of the simulated merger tree approach first proposed by Moster et al. This method combines N-body cosmological simulations and semi-analytic techniques to extract realistic initial conditions for galaxy mergers. These are then evolved using high resolution hydrodynamical simulations, which include dark matter, stars, cold gas in the disc and hot gas in the halo. We show that the satellite mass accretion is not as effective as previously thought, as there is substantial stellar stripping before the final merger. The fraction of stellar disc mass transferred to the bulge is quite low, even in the case of a major merger, mainly due to the dispersion of part of the stellar disc mass into the halo. We confirm the findings of Hopkins et al., that a gas rich disc is able to survive major mergers more efficiently. The enhanced star formation associated with the merger is not localised to the bulge of galaxy, but a substantial fraction takes place in the disc too. The inclusion of the hot gas reservoir in the galaxy model contributes to reducing the efficiency of bulge formation. Overall, our findings suggest that mergers are not as efficient as previously thought in transforming discs into bulges. This possibly alleviates some of the tensions between observations of bulgeless galaxies and the hierarchical scenario for structure formation.

preprint2015arXiv

Star formation in mergers with cosmologically motivated initial conditions

We use semi-analytic models and cosmological merger trees to provide the initial conditions for multi-merger numerical hydrodynamic simulations, and exploit these simulations to explore the effect of galaxy interaction and merging on star formation (SF). We compute numerical realisations of twelve merger trees from z=1.5 to z=0. We include the effects of the large hot gaseous halo around all galaxies, following recent obervations and predictions of galaxy formation models. We find that including the hot gaseous halo has a number of important effects. Firstly, as expected, the star formation rate on long timescales is increased due to cooling of the hot halo and refuelling of the cold gas reservoir. Secondly, we find that interactions do not always increase the SF in the long term. This is partially due to the orbiting galaxies transferring gravitational energy to the hot gaseous haloes and raising their temperature. Finally we find that the relative size of the starburst, when including the hot halo, is much smaller than previous studies showed. Our simulations also show that the order and timing of interactions are important for the evolution of a galaxy. When multiple galaxies interact at the same time, the SF enhancement is less than when galaxies interact in series. All these effects show the importance of including hot gas and cosmologically motivated merger trees in galaxy evolution models.

preprint2015arXiv

The impact of mechanical AGN feedback on the formation of massive early-type galaxies

We employ cosmological hydrodynamical simulations to investigate the effects of AGN feedback on the formation of massive galaxies with present-day stellar masses of $M_{stel} = 8.8 \times 10^{10} - 6.0 \times 10^{11} M_{sun}$. Using smoothed particle hydrodynamics simulations with a pressure-entropy formulation that allows an improved treatment of contact discontinuities and fluid mixing, we run three sets of simulations of 20 halos with different AGN feedback models: (1) no feedback, (2) thermal feedback, and (3) mechanical and radiation feedback. We assume that seed black holes are present at early cosmic epochs at the centre of emerging dark matter halos and trace their mass growth via gas accretion and mergers with other black holes. Both feedback models successfully recover the observed M_BH - sigma relation and black hole-to-stellar mass ratio for simulated central early-type galaxies. The baryonic conversion efficiencies are reduced by a factor of two compared to models without any AGN feedback at all halo masses. However, massive galaxies simulated with thermal AGN feedback show a factor of ~10-100 higher X-ray luminosities than observed. The mechanical/radiation feedback model reproduces the observed correlation between X-ray luminosities and velocity dispersion, e.g. for galaxies with sigma = 200 km/s, the X-ray luminosity is reduced from $10^{42}$ erg/s to $10^{40}$ erg/s. It also efficiently suppresses late time star formation, reducing the specific star formation rate from $10^{-10.5}$ $yr^{-1}$ to $10^{-14}$ $yr^{-1}$ on average and resulting in quiescent galaxies since z=2, whereas the thermal feedback model shows higher late time in-situ star formation rates than observed.

preprint2014arXiv

Consequences of Mechanical and Radiative Feedback from Black Holes in Disc Galaxy Mergers

We study the effect of AGN mechanical and radiation feedback on the formation of bulge dominated galaxies via mergers of disc galaxies. The merging galaxies have mass-ratios of 1:1 to 6:1 and include pre-existing hot gaseous halos to properly account for the global impact of AGN feedback. Using smoothed particle hydrodynamics simulation code (GADGET-3) we compare three models with different AGN feedback models: (1) no black hole and no AGN feedback; (2) thermal AGN feedback; and (3) mechanical and radiative AGN feedback. The last model is motivated by observations of broad line quasars which show winds with initial velocities of $v_w \ge$ 10,000 km/s and also heating associated with the central AGN X-ray radiation. The primary changes in gas properties due to mechanical AGN feedback are lower thermal X-ray luminosity from the final galaxy - in better agreement with observations - and galactic outflows with higher velocity $\sim 1000$ km/s similar to recent direct observations of nearby merger remnants. The kinetic energy of the outflowing gas is a factor of $\sim$ 20 higher than in the thermal feedback case. All merger remnants with momentum-based AGN feedback with $v_w \sim 10,000$ km/s and $ε_w=2 \times 10^{-3}$, independent of their progenitor mass-ratios, reproduce the observed relations between stellar velocity dispersion and black hole mass ($M_{\rm bh} - σ$) as well as X-ray luminosity ($L_X - σ$) with $10^{37.5}$ erg/s $\lesssim L_X (0.3-8~{\rm keV}) \lesssim 10^{39.2}$ erg/s for velocity dispersions in the range of 120 km/s $\lesssim σ\lesssim$ 190 km/s. In addition, the mechanical feedback produces a much greater AGN variability. We also show that gas is more rapidly and impulsively stripped from the galactic centres driving a moderate increase in galaxy size and decrease in central density with the mechanical AGN feedback model.

preprint2014arXiv

SPHGal: Smoothed Particle Hydrodynamics with improved accuracy for Galaxy simulations

We present the smoothed-particle hydrodynamics implementation SPHGal, which combines some recently proposed improvements in GADGET. This includes a pressure-entropy formulation with a Wendland kernel, a higher order estimate of velocity gradients, a modified artificial viscosity switch with a modified strong limiter, and artificial conduction of thermal energy. With a series of idealized hydrodynamic tests we show that the pressure-entropy formulation is ideal for resolving fluid mixing at contact discontinuities but performs conspicuously worse at strong shocks due to the large entropy discontinuities. Including artificial conduction at shocks greatly improves the results. In simulations of Milky Way like disk galaxies a feedback-induced instability develops if too much artificial viscosity is introduced. Our modified artificial viscosity scheme prevents this instability and shows efficient shock capturing capability. We also investigate the star formation rate and the galactic outflow. The star formation rates vary slightly for different SPH schemes while the mass loading is sensitive to the SPH scheme and significantly reduced in our favored implementation. We compare the accretion behavior of the hot halo gas. The formation of cold blobs, an artifact of simple SPH implementations, can be eliminated efficiently with proper fluid mixing, either by conduction and/or by using a pressure-entropy formulation.

preprint2013arXiv

The Dynamics of Galaxy Pairs in a Cosmological Setting

We use the Millennium Simulation, and an abundance-matching framework, to investigate the dynamical behaviour of galaxy pairs embedded in a cosmological context. Our main galaxy-pair sample, selected to have separations under 250 kpc/h, consists of over 1.3 million pairs at redshift z = 0, with stellar masses greater than 10^9 Msun, probing mass ratios down to 1:1000. We use dark matter halo membership and energy to classify our galaxy pairs. In terms of halo membership, central-satellite pairs tend to be in isolation (in relation to external more massive galaxies), are energetically- bound to each other, and are also weakly-bound to a neighbouring massive galaxy. Satellite-satellite pairs, instead, inhabit regions in close proximity to a more massive galaxy, are energetically-unbound, and are often bound to that neighbour. We find that 60% of our paired galaxies are bound to both their companion and to a third external object. Moreover, only 9% of our pairs resemble the kind of systems described by idealised binary merger simulations in complete isolation. In sum, we demonstrate the importance of properly connecting galaxy pairs to the rest of the Universe.

preprint2012arXiv

Galactic star formation and accretion histories from matching galaxies to dark matter haloes

We present a new statistical method to determine the relationship between the stellar masses of galaxies and the masses of their host dark matter haloes over the entire cosmic history from z~4 to the present. This multi-epoch abundance matching (MEAM) model self-consistently takes into account that satellite galaxies first become satellites at times earlier than they are observed. We employ a redshift-dependent parameterization of the stellar-to-halo mass relation to populate haloes and subhaloes in the Millennium simulations with galaxies, requiring that the observed stellar mass functions at several redshifts be reproduced simultaneously. Using merger trees extracted from the dark matter simulations in combination with MEAM, we predict the average assembly histories of galaxies, separating into star formation within the galaxies (in-situ) and accretion of stars (ex-situ). The peak star formation efficiency decreases with redshift from 23% at z=0 to 9% at z=4 while the corresponding halo mass increases from 10^11.8M\odot to 10^12.5M\odot. The star formation rate of central galaxies peaks at a redshift which depends on halo mass; for massive haloes this peak is at early cosmic times while for low-mass galaxies the peak has not been reached yet. In haloes similar to that of the Milky-Way about half of the central stellar mass is assembled after z=0.7. In low-mass haloes, the accretion of satellites contributes little to the assembly of their central galaxies, while in massive haloes more than half of the central stellar mass is formed ex-situ with significant accretion of satellites at z<2. We find that our method implies a cosmic star formation history and an evolution of specific star formation rates which are consistent with those inferred directly. We present convenient fitting functions for stellar masses, star formation rates, and accretion rates as functions of halo mass and redshift.

preprint2012arXiv

Numerical hydrodynamic simulations based on semi-analytic galaxy merger trees: method and Milky-Way like galaxies

We present a new approach to study galaxy evolution in a cosmological context. We combine cosmological merger trees and semi-analytic models of galaxy formation to provide the initial conditions for multi-merger hydrodynamic simulations. In this way we exploit the advantages of merger simulations (high resolution and inclusion of the gas physics) and semi-analytic models (cosmological background and low computational cost), and integrate them to create a novel tool. This approach allows us to study the evolution of various galaxy properties, including the treatment of the hot gaseous halo from which gas cools and accretes onto the central disc, which has been neglected in many previous studies. This method shows several advantages over other methods. As only the particles in the regions of interest are included, the run time is much shorter than in traditional cosmological simulations, leading to greater computational efficiency. Using cosmological simulations, we show that multiple mergers are expected to be more common than sequences of isolated mergers, and therefore studies of galaxy mergers should take this into account. In this pilot study, we present our method and illustrate the results of simulating ten Milky Way-like galaxies since z=1. We find good agreement with observations for the total stellar masses, star formation rates, cold gas fractions and disc scale length parameters. We expect that this novel numerical approach will be very useful for pursuing a number of questions pertaining to the transformation of galaxy internal structure through cosmic time.

preprint2011arXiv

Interaction between dark matter sub-halos & galactic gaseous disk

We investigate the idea that the interaction of Dark Matter (DM) sub-halos with the gaseous disk of galaxies can be the origin for the observed holes and shells found in their neutral hydrogen (HI) distributions. We use high resolution hydrodynamic simulations to show that pure DM sub-halos impacting a galactic disk are not able to produce holes; on the contrary, they result in high density regions in the disk. However, sub-halos containing a small amount of gas (a few percent of the total DM mass of the sub-halo) are able to displace the gas in the disk and form holes and shells. The size and lifetime of these holes depend on the sub-halo gas mass, density and impact velocity. A DM sub-halo of mass $10^8$ $M_{\odot}$ and a gas mass fraction of $\sim 3$%, is able to create a kpc scale hole, with a lifetime similar to those observed in nearby galaxies. We also register an increase in the star formation rate at the rim of the hole, again in agreement with observations. Even though the properties off these simulated structures resemble those found in observations we find that the number of predicted holes (based on mass and orbital distributions of DM halos derived from cosmological N-body simulations) falls short compared to the observations. Only a handful of holes are produced per Giga year. This leads us to conclude that DM halo impact is not the major channel through which these holes are formed.

preprint2011arXiv

The effects of a hot gaseous halo in galaxy major mergers

Cosmological hydrodynamical simulations as well as observations indicate that spiral galaxies are comprised of five different components: dark matter halo, stellar disc, stellar bulge, gaseous disc and gaseous halo. While the first four components have been extensively considered in numerical simulations of binary galaxy mergers, the effect of a hot gaseous halo has usually been neglected even though it can contain up to 80% of the total gas within the galaxy virial radius. We present a series of hydrodynamic simulations of major mergers of disc galaxies, that for the first time include a diffuse, rotating, hot gaseous halo. Through cooling and accretion, the hot halo can dissipate and refuel the cold gas disc before and after a merger. This cold gas can subsequently form stars, thus impacting the morphology and kinematics of the remnant. Simulations of isolated systems with total mass M~10^12Msun show a nearly constant star formation rate of ~5Msun/yr if the hot gaseous halo is included, while the star formation rate declines exponentially if it is neglected. We conduct a detailed study of the star formation efficiency during mergers and find that the presence of a hot gaseous halo reduces the starburst efficiency (e=0.5) compared to simulations without a hot halo (e=0.68). Moreover we find cases where the stellar mass of the merger remnant is lower than the sum of the stellar mass of the two progenitor galaxies when evolved in isolation. This suggests a revision to semi-analytic galaxy formation models which assume that a merger always leads to enhanced star formation. We show that adding the hot gas component has a significant effect on the kinematics and internal structure of the merger remnants, like an increased abundance of fast rotators and an r^(1/4) surface brightness profile at small scales.

preprint2011arXiv

The effects of a hot gaseous halo on disc thickening in galaxy minor mergers

We employ hydrodynamical simulations to study the effects of dissipational gas physics on the vertical heating and thickening of disc galaxies during minor mergers. For the first time we present a suite of simulations that includes a diffuse, rotating, cooling, hot gaseous halo, as predicted by cosmological hydrodynamical simulations as well as models of galaxy formation. We study the effect of this new gaseous component on the vertical structure of a Milky Way-like stellar disc during 1:10 and 1:5 mergers. For 1:10 mergers we find no increased final thin disc scale height compared to the isolated simulation, leading to the conclusion that thin discs can be present even after a 1:10 merger if a reasonable amount of hot gas is present. The reason for this is the accretion of new cold gas, leading to the formation of a massive new thin stellar disc that dominates the surface brightness profile. In a previous study, in which we included only cold gas in the disk, we showed that the presence of cold gas decreased the thickening by a minor merger relative to the no-gas case. Here, we show that the evolution of the scale height in the presences of a cooling hot halo is dominated by the formation of the new stellar disc. In this scenario, the thick disc is the old stellar disc that has been thickened in a minor merger at z>1, while the thin disc is the new stellar disc that reforms after this merger. In addition, we study the evolution of the scale height during a 1:5 merger and find that a thin disc can be present even after this merger, provided enough hot gas is available. The final scale height in our simulations depends on the mass of the hot gaseous halo, the efficiency of the winds and the merger mass ratio. We find post-merger values in the range 0.5<z0<1.0 kpc in good agreement with observational constraints by local galaxies.

preprint2010arXiv

A Cosmic Variance Cookbook

Deep pencil beam surveys (<1 deg^2) are of fundamental importance for studying the high-redshift universe. However, inferences about galaxy population properties are in practice limited by 'cosmic variance'. This is the uncertainty in observational estimates of the number density of galaxies arising from the underlying large-scale density fluctuations. This source of uncertainty can be significant, especially for surveys which cover only small areas and for massive high-redshift galaxies. Cosmic variance for a given galaxy population can be determined using predictions from cold dark matter theory and the galaxy bias. In this paper we provide tools for experiment design and interpretation. For a given survey geometry we present the cosmic variance of dark matter as a function of mean redshift z and redshift bin size Dz. Using a halo occupation model to predict galaxy clustering, we derive the galaxy bias as a function of mean redshift for galaxy samples of a given stellar mass range. In the linear regime, the cosmic variance of these galaxy samples is the product of the galaxy bias and the dark matter cosmic variance. We present a simple recipe using a fitting function to compute cosmic variance as a function of the angular dimensions of the field, z, Dz and stellar mass m*. We also provide tabulated values and a software tool. We find that for GOODS at z=2 and with Dz=0.5 the relative cosmic variance of galaxies with m*>10^11 Msun is ~38%, while it is ~27% for GEMS and ~12% for COSMOS. For galaxies of m*~10^10 Msun the relative cosmic variance is ~19% for GOODS, ~13% for GEMS and ~6% for COSMOS. This implies that cosmic variance is a significant source of uncertainty at z=2 for small fields and massive galaxies, while for larger fields and intermediate mass galaxies cosmic variance is less serious.

preprint2009arXiv

Can Gas prevent the Destruction of Thin Stellar Discs by Minor Mergers?

We study the effect of dissipational gas physics on the vertical heating and thickening of disc galaxies during minor mergers. We produce a suite of minor merger simulations for Milky Way-like galaxies. This suite consists of collisionless simulations as well as hydrodynamical runs including a gaseous component in the galactic disc. We find that in dissipationless simulations minor mergers cause the scale height of the disc to increase by up to a factor of ~2. When the presence of gas in the disc is taken into account this thickening is reduced by 25% (50%) for an initial disc gas fraction of 20% (40%), leading to a final scale height z0 between 0.6 and 0.7 kpc, regardless of the initial scale height. We argue that the presence of gas reduces disc heating via two mechanisms: absorption of kinetic impact energy by the gas and/or formation of a new thin stellar disc that can cause heated stars to recontract towards the disc plane. We show that in our simulations most of the gas is consumed during the merger and thus the regrowth of a new thin disc has a negligible impact on the z0 of the post merger galaxy. Final disc scale heights found in our simulations are in good agreement with studies of the vertical structure of spiral galaxies where the majority of the systems are found to have scale heights of 0.4 kpc < z0 < 0.8 kpc. We also found no tension between recent measurements of the scale height of the Milky Way thin disc and results coming from our hydrodynamical simulations. We conclude that the existence of a thin disc in the Milky Way and in external galaxies is not in obvious conflict with the predictions of the CDM model.

preprint2009arXiv

Constraints on the relationship between stellar mass and halo mass at low and high redshift

We use a statistical approach to determine the relationship between the stellar masses of galaxies and the masses of the dark matter halos in which they reside. We obtain a parameterized stellar-to-halo mass (SHM) relation by populating halos and subhalos in an N-body simulation with galaxies and requiring that the observed stellar mass function be reproduced. We find good agreement with constraints from galaxy-galaxy lensing and predictions of semi-analytic models. Using this mapping, and the positions of the halos and subhalos obtained from the simulation, we find that our model predictions for the galaxy two-point correlation function (CF) as a function of stellar mass are in excellent agreement with the observed clustering properties in the SDSS at z=0. We show that the clustering data do not provide additional strong constraints on the SHM function and conclude that our model can therefore predict clustering as a function of stellar mass. We compute the conditional mass function, which yields the average number of galaxies with stellar masses in the range [m, m+dm] that reside in a halo of mass M. We study the redshift dependence of the SHM relation and show that, for low mass halos, the SHM ratio is lower at higher redshift. The derived SHM relation is used to predict the stellar mass dependent galaxy CF and bias at high redshift. Our model predicts that not only are massive galaxies more biased than low mass ones at all redshifts, but the bias increases more rapidly with increasing redshift for massive galaxies than for low mass ones. We present convenient fitting functions for the SHM relation as a function of redshift, the conditional mass function, and the bias as a function of stellar mass and redshift.