Source author record

Michael T. Busha

Michael T. Busha 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

18works
4topics
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

18 published item(s)

preprint2021arXiv

ADDGALS: Simulated Sky Catalogs for Wide Field Galaxy Surveys

We present a method for creating simulated galaxy catalogs with realistic galaxy luminosities, broad-band colors, and projected clustering over large cosmic volumes. The technique, denoted ADDGALS (Adding Density Dependent GAlaxies to Lightcone Simulations), uses an empirical approach to place galaxies within lightcone outputs of cosmological simulations. It can be applied to significantly lower-resolution simulations than those required for commonly used methods such as halo occupation distributions, subhalo abundance matching, and semi-analytic models, while still accurately reproducing projected galaxy clustering statistics down to scales of r ~ 100 kpc/h. We show that \addgals\ catalogs reproduce several statistical properties of the galaxy distribution as measured by the Sloan Digital Sky Survey (SDSS) main galaxy sample, including galaxy number densities, observed magnitude and color distributions, as well as luminosity- and color-dependent clustering. We also compare to cluster-galaxy cross correlations, where we find significant discrepancies with measurements from SDSS that are likely linked to artificial subhalo disruption in the simulations. Applications of this model to simulations of deep wide-area photometric surveys, including modeling weak-lensing statistics, photometric redshifts, and galaxy cluster finding are presented in DeRose et al (2019), and an application to a full cosmology analysis of Dark Energy Survey (DES) Year 3 like data is presented in DeRose etl al (2021). We plan to publicly release a 10,313 square degree catalog constructed using ADDGALS with magnitudes appropriate for several existing and planned surveys, including SDSS, DES, VISTA, WISE, and LSST.

preprint2014arXiv

Spectroscopic failures in photometric redshift calibration: cosmological biases and survey requirements

We use N-body-spectro-photometric simulations to investigate the impact of incompleteness and incorrect redshifts in spectroscopic surveys to photometric redshift training and calibration and the resulting effects on cosmological parameter estimation from weak lensing shear-shear correlations. The photometry of the simulations is modeled after the upcoming Dark Energy Survey and the spectroscopy is based on a low/intermediate resolution spectrograph with wavelength coverage of 5500Å < λ < 9500Å. The principal systematic errors that such a spectroscopic follow-up encounters are incompleteness (inability to obtain spectroscopic redshifts for certain galaxies) and wrong redshifts. Encouragingly, we find that a neural network-based approach can effectively describe the spectroscopic incompleteness in terms of the galaxies' colors, so that the spectroscopic selection can be applied to the photometric sample. Hence, we find that spectroscopic incompleteness yields no appreciable biases to cosmology, although the statistical constraints degrade somewhat because the photometric survey has to be culled to match the spectroscopic selection. Unfortunately, wrong redshifts have a more severe impact: the cosmological biases are intolerable if more than a percent of the spectroscopic redshifts are incorrect. Moreover, we find that incorrect redshifts can also substantially degrade the accuracy of training set based photo-z estimators. The main problem is the difficulty of obtaining redshifts, either spectroscopically or photometrically, for objects at z > 1.3. We discuss several approaches for reducing the cosmological biases, in particular finding that photo-z error estimators can reduce biases appreciably.

preprint2013arXiv

Gravitationally Consistent Halo Catalogs and Merger Trees for Precision Cosmology

We present a new algorithm for generating merger trees and halo catalogs which explicitly ensures consistency of halo properties (mass, position, and velocity) across timesteps. Our algorithm has demonstrated the ability to improve both the completeness (through detecting and inserting otherwise missing halos) and purity (through detecting and removing spurious objects) of both merger trees and halo catalogs. In addition, our method is able to robustly measure the self-consistency of halo finders; it is the first to directly measure the uncertainties in halo positions, halo velocities, and the halo mass function for a given halo finder based on consistency between snapshots in cosmological simulations. We use this algorithm to generate merger trees for two large simulations (Bolshoi and Consuelo) and evaluate two halo finders (ROCKSTAR and BDM). We find that both the ROCKSTAR and BDM halo finders track halos extremely well; in both, the number of halos which do not have physically consistent progenitors is at the 1-2% level across all halo masses. Our code is publicly available at http://code.google.com/p/consistent-trees . Our trees and catalogs are publicly available at http://hipacc.ucsc.edu/Bolshoi/ .

preprint2013arXiv

Mergers and Mass Accretion for Infalling Halos Both End Well Outside Cluster Virial Radii

We find that infalling dark matter halos (i.e., the progenitors of satellite halos) begin losing mass well outside the virial radius of their eventual host halos. The peak mass occurs at a range of clustercentric distances, with median and 68th percentile range of 1.8 +2.3/-1.0 R_(vir,host) for progenitors of z=0 satellites. The peak circular velocity for infalling halos occurs at significantly larger distances (3.7 +3.3/-2.2 R_(vir,host) at z=0). This difference arises because different physical processes set peak circular velocity (typically, ~1:5 and larger mergers which cause transient circular velocity spikes) and peak mass (typically, smooth accretion) for infalling halos. We find that infalling halos also stop having significant mergers well before they enter the virial radius of their eventual hosts. Mergers larger than a 1:40 ratio in halo mass end for infalling halos at similar clustercentric distances (~ 1.9 R_(vir,host)) as the end of overall mass accretion. However, mergers larger than 1:3 typically end for infalling halos at more than 4 virial radial away from their eventual hosts. This limits the ability of mergers to affect quenching and morphology changes in clusters. We also note that the transient spikes which set peak circular velocity may lead to issues with abundance matching on that parameter, including unphysical galaxy stellar mass growth profiles near clusters; we propose a simple observational test to check if a better halo proxy for galaxy stellar mass exists.

preprint2012arXiv

A High Throughput Workflow Environment for Cosmological Simulations

The next generation of wide-area sky surveys offer the power to place extremely precise constraints on cosmological parameters and to test the source of cosmic acceleration. These observational programs will employ multiple techniques based on a variety of statistical signatures of galaxies and large-scale structure. These techniques have sources of systematic error that need to be understood at the percent-level in order to fully leverage the power of next-generation catalogs. Simulations of large-scale structure provide the means to characterize these uncertainties. We are using XSEDE resources to produce multiple synthetic sky surveys of galaxies and large-scale structure in support of science analysis for the Dark Energy Survey. In order to scale up our production to the level of fifty 10^10-particle simulations, we are working to embed production control within the Apache Airavata workflow environment. We explain our methods and report how the workflow has reduced production time by 40% compared to manual management.

preprint2012arXiv

Luminous Satellites II: Spatial Distribution, Luminosity Function and Cosmic Evolution

We infer the normalization and the radial and angular distributions of the number density of satellites of massive galaxies ($\log_{10}[M_{h}^*/M\odot]>10.5$) between redshifts 0.1 and 0.8 as a function of host stellar mass, redshift, morphology and satellite luminosity. Exploiting the depth and resolution of the COSMOS HST images, we detect satellites up to eight magnitudes fainter than the host galaxies and as close as 0.3 (1.4) arcseconds (kpc). Describing the number density profile of satellite galaxies to be a projected power law such that $P(R)\propto R^{\rpower}$, we find $\rpower=-1.1\pm 0.3$. We find no dependency of $\rpower$ on host stellar mass, redshift, morphology or satellite luminosity. Satellites of early-type hosts have angular distributions that are more flattened than the host light profile and are aligned with its major axis. No significant average alignment is detected for satellites of late-type hosts. The number of satellites within a fixed magnitude contrast from a host galaxy is dependent on its stellar mass, with more massive galaxies hosting significantly more satellites. Furthermore, high-mass late-type hosts have significantly fewer satellites than early-type galaxies of the same stellar mass, likely a result of environmental differences. No significant evolution in the number of satellites per host is detected. The cumulative luminosity function of satellites is qualitatively in good agreement with that predicted using subhalo abundance matching techniques. However, there are significant residual discrepancies in the absolute normalization, suggesting that properties other than the host galaxy luminosity or stellar mass determine the number of satellites.

preprint2012arXiv

Sample variance in photometric redshift calibration: cosmological biases and survey requirements

We use N-body/photometric galaxy simulations to examine the impact of sample variance of spectroscopic redshift samples on the accuracy of photometric redshift (photo-z) determination and calibration of photo-z errors. We estimate the biases in the cosmological parameter constraints from weak lensing and derive requirements on the spectroscopic follow-up for three different photo-z algorithms chosen to broadly span the range of algorithms available. We find that sample variance is much more relevant for the photo-z error calibration than for photo-z training, implying that follow-up requirements are similar for different algorithms. We demonstrate that the spectroscopic sample can be used for training of photo-zs and error calibration without incurring additional bias in the cosmological parameters. We provide a guide for observing proposals for the spectroscopic follow-up to ensure that redshift calibration biases do not dominate the cosmological parameter error budget. For example, assuming optimistically (pessimistically) that the weak lensing shear measurements from the Dark Energy Survey could obtain 1-sigma constraints on the dark energy equation of state w of 0.035 (0.055), implies a follow-up requirement of 150 (40) patches of sky with a telescope such as Magellan, assuming a 1/8^2 deg effective field of view and 400 galaxies per patch. Assuming (optimistically) a VVDS-like spectroscopic completeness with purely random failures, this could be accomplished with about 75 (20) nights of observation. For more realistic assumptions regarding spectroscopic completeness, or in the presence of other sources of systematics not considered here, further degradations to dark energy constraints are possible. We test several approaches for reducing the requirements. Abridged

preprint2011arXiv

Galaxies in X-ray Groups I: Robust Membership Assignment and the Impact of Group Environments on Quenching

Understanding the mechanisms that lead dense environments to host galaxies with redder colors, more spheroidal morphologies, and lower star formation rates than field populations remains an important problem. As most candidate processes ultimately depend on host halo mass, accurate characterizations of the local environment, ideally tied to halo mass estimates and spanning a range in halo mass and redshift are needed. In this work, we present and test a rigorous, probabalistic method for assigning galaxies to groups based on precise photometric redshifts and X-ray selected groups drawn from the COSMOS field. The groups have masses in the range 10^13 < M_200c/M_sun < 10^14 and span redshifts 0<z<1. We characterize our selection algorithm via tests on spectroscopic subsamples, including new data obtained at the VLT, and by applying our method to detailed mock catalogs. We find that our group member galaxy sample has a purity of 84% and completeness of 92% within 0.5 R200c. We measure the impact of uncertainties in redshifts and group centering on the quality of the member selection with simulations based on current data as well as future imaging and spectroscopic surveys. As a first application of our new group member catalog which will be made publicly available, we show that member galaxies exhibit a higher quenched fraction compared to the field at fixed stellar mass out to z~1, indicating a significant relationship between star formation and environment at group scales. We also address the suggestion that dusty star forming galaxies in such groups may impact the high-l power spectrum of the cosmic microwave background and find that such a population cannot explain the low power seen in recent SZ measurements.

preprint2011arXiv

How Common are the Magellanic Clouds?

We introduce a probabilistic approach to the problem of counting dwarf satellites around host galaxies in databases with limited redshift information. This technique is used to investigate the occurrence of satellites with luminosities similar to the Magellanic Clouds around hosts with properties similar to the Milky Way in the object catalog of the Sloan Digital Sky Survey. Our analysis uses data from SDSS Data Release 7, selecting candidate Milky-Way-like hosts from the spectroscopic catalog and candidate analogs of the Magellanic Clouds from the photometric catalog. Our principal result is the probability for a Milky-Way-like galaxy to host N_{sat} close satellites with luminosities similar to the Magellanic Clouds. We find that 81 percent of galaxies like the Milky Way are have no such satellites within a radius of 150 kpc, 11 percent have one, and only 3.5 percent of hosts have two. The probabilities are robust to changes in host and satellite selection criteria, background-estimation technique, and survey depth. These results demonstrate that the Milky Way has significantly more satellites than a typical galaxy of its luminosity; this fact is useful for understanding the larger cosmological context of our home galaxy.

preprint2011arXiv

New constraints on the evolution of the stellar-to-dark matter connection: a combined analysis of galaxy-galaxy lensing, clustering, and stellar mass functions from z=0.2 to z=1

Using data from the COSMOS survey, we perform the first joint analysis of galaxy-galaxy weak lensing, galaxy spatial clustering, and galaxy number densities. Carefully accounting for sample variance and for scatter between stellar and halo mass, we model all three observables simultaneously using a novel and self-consistent theoretical framework. Our results provide strong constraints on the shape and redshift evolution of the stellar-to-halo mass relation (SHMR) from z=0.2 to z=1. At low stellar mass, we find that halo mass scales as Mh M*^0.46 and that this scaling does not evolve significantly with redshift to z=1. We show that the dark-to-stellar ratio, Mh/M*, varies from low to high masses, reaching a minimum of Mh/M*~27 at M*=4.5x10^10 Msun and Mh=1.2x10^12 Msun. This minimum is important for models of galaxy formation because it marks the mass at which the accumulated stellar growth of the central galaxy has been the most efficient. We describe the SHMR at this minimum in terms of the "pivot stellar mass", M*piv, the "pivot halo mass", Mhpiv, and the "pivot ratio", (Mh/M*)piv. Thanks to a homogeneous analysis of a single data set, we report the first detection of mass downsizing trends for both Mhpiv and M*piv. The pivot stellar mass decreases from M*piv=5.75+-0.13x10^10 Msun at z=0.88 to M*piv=3.55+-0.17x10^10 Msun at z=0.37. Intriguingly, however, the corresponding evolution of Mhpiv leaves the pivot ratio constant with redshift at (Mh/M*)piv~27. We use simple arguments to show how this result raises the possibility that star formation quenching may ultimately depend on Mh/M* and not simply Mh, as is commonly assumed. We show that simple models with such a dependence naturally lead to downsizing in the sites of star formation. Finally, we discuss the implications of our results in the context of popular quenching models, including disk instabilities and AGN feedback.

preprint2011arXiv

Statistics of Satellite Galaxies Around Milky Way-Like Hosts

We calculate the probability that a Milky-Way-like halo in the standard cosmological model has the observed number of Magellanic Clouds (MCs). The statistics of the number of MCs in the LCDM model are in good agreement with observations of a large sample of SDSS galaxies. Under the sub-halo abundance matching assumption of a relationship with small scatter between galaxy r-band luminosities and halo internal velocities v_max, we make detailed comparisons to similar measurements using SDSS DR7 data by Liu et al. (2010). Models and observational data give very similar probabilities for having zero, one, and two MC-like satellites. In both cases, Milky Way-luminosity hosts have just a \sim 10% chance of hosting two satellites similar to the Magellanic Clouds. In addition, we present a prediction for the probability for a host galaxy to have Nsats satellite galaxies as a function of the magnitudes of both the host and satellite. This probability and its scaling with host properties is significantly different from that of mass-selected objects because of scatter in the mass- luminosity relation and because of variations in the star formation efficiency with halo mass.

preprint2011arXiv

The Mass Distribution and Assembly of the Milky Way from the Properties of the Magellanic Clouds

We present a new measurement of the mass of the Milky Way (MW) based on observed properties of its largest satellite galaxies, the Magellanic Clouds (MCs), and an assumed prior of a ΛCDM universe. The large, high-resolution Bolshoi cosmological simulation of this universe provides a means to statistically sample the dynamical properties of bright satellite galaxies in a large population of dark matter halos. The observed properties of the MCs, including their circular velocity, distance from the center of the MW, and velocity within the MW halo, are used to evaluate the likelihood that a given halo would have each or all of these properties; the posterior PDF for any property of the MW system can thus be constructed. This method provides a constraint on the MW virial mass, 1.2 +0.7 -0.3(stat.) +0.3 -0.4 (sys.) x 10^12 M\odot (68% confidence), which is consistent with recent determinations that involve very different assumptions. In addition, we calculate the posterior PDF for the density profile of the MW and its satellite accretion history. Although typical satellites of 10^12 M\odot halos are accreted over a wide range of epochs over the last 10 Gyr, we find a \sim72% probability that the Magellanic Clouds were accreted within the last Gyr, and a 50% probability that they were accreted together.

preprint2010arXiv

The Voronoi Tessellation cluster finder in 2+1 dimensions

We present a detailed description of the Voronoi Tessellation (VT) cluster finder algorithm in 2+1 dimensions, which improves on past implementations of this technique. The need for cluster finder algorithms able to produce reliable cluster catalogs up to redshift 1 or beyond and down to $10^{13.5}$ solar masses is paramount especially in light of upcoming surveys aiming at cosmological constraints from galaxy cluster number counts. We build the VT in photometric redshift shells and use the two-point correlation function of the galaxies in the field to both determine the density threshold for detection of cluster candidates and to establish their significance. This allows us to detect clusters in a self consistent way without any assumptions about their astrophysical properties. We apply the VT to mock catalogs which extend to redshift 1.4 reproducing the $Λ$CDM cosmology and the clustering properties observed in the SDSS data. An objective estimate of the cluster selection function in terms of the completeness and purity as a function of mass and redshift is as important as having a reliable cluster finder. We measure these quantities by matching the VT cluster catalog with the mock truth table. We show that the VT can produce a cluster catalog with completeness and purity $>80%$ for the redshift range up to $\sim 1$ and mass range down to $\sim 10^{13.5}$ solar masses.

preprint2009arXiv

ArborZ: Photometric Redshifts Using Boosted Decision Trees

Precision photometric redshifts will be essential for extracting cosmological parameters from the next generation of wide-area imaging surveys. In this paper we introduce a photometric redshift algorithm, ArborZ, based on the machine-learning technique of Boosted Decision Trees. We study the algorithm using galaxies from the Sloan Digital Sky Survey and from mock catalogs intended to simulate both the SDSS and the upcoming Dark Energy Survey. We show that it improves upon the performance of existing algorithms. Moreover, the method naturally leads to the reconstruction of a full probability density function (PDF) for the photometric redshift of each galaxy, not merely a single "best estimate" and error, and also provides a photo-z quality figure-of-merit for each galaxy that can be used to reject outliers. We show that the stacked PDFs yield a more accurate reconstruction of the redshift distribution N(z). We discuss limitations of the current algorithm and ideas for future work.

preprint2009arXiv

The Impact of Inhomogeneous Reionization on the Satellite Galaxy Population of the Milky Way

We use the publicly available subhalo catalogs from the Via Lactea simulation along with a Gpc-scale N-body simulation to understand the impact of inhomogeneous reionization on the satellite galaxy population of the Milky Way. The large-volume simulation is combined with a model for reionization that allows us to predict the distribution of reionization times for Milky Way mass halos. Motivated by this distribution, we identify candidate satellite galaxies in the simulation by requiring that any subhalo must grow above a specified mass threshold before it is reionized; after this time the photoionizing background will suppress both the formation of stars and the accretion of gas. We show that varying the reionization time over the range expected for Milky Way mass halos can change the number of satellite galaxies by roughly two orders of magnitude. This conclusion is in contradiction with a number of studies in the literature, and we conclude that this is a result of inconsistent application of the results of Gnedin (2000). We compare our satellite galaxies to observations using both abundance matching and stellar population synthesis methods to assign luminosities to our subhalos and account for observational completeness effects. Additionally, if we assume that the mass threshold is set by the virial temperature Tvir = 8e3K we find that our model accurately matches the vmax distribution, radial distribution, and luminosity function of observed Milky Way satellites for a reionization time zreion = 9.6^{1.0}_{-2.1}, assuming that the Via Lacteasubhalo distribution is representative of the Milky Way. This results in the presence of 119^{+202}_{-50} satellite galaxies.

preprint2006arXiv

The Asymptotic Form of Cosmic Structure: Small Scale Power and Accretion History

We explore the effects of small scale structure on the formation and equilibrium of dark matter halos in a universe dominated by vacuum energy. We present the results of a suite of four N-body simulations, two with a LCDM initial power spectrum and two with WDM-like spectra that suppress the early formation of small structures. All simulations are run into to far future when the universe is 64Gyr/h old, long enough for halos to essentially reach dynamical equilibrium. We quantify the importance of hierarchical merging on the halo mass accretion history, the substructure population, and the equilibrium density profile. We modify the mass accretion history function of Wechsler et al. (2002) by introducing a parameter, γ, that controls the rate of mass accretion, dln(M) / dln(a) ~ a^(-γ), and find that this form characterizes both hierarchical and monolithic formation. Subhalo decay rates are exponential in time with a much shorter time scale for WDM halos. At the end of the simulations, we find truncated Hernquist density profiles for halos in both the CDM and WDM cosmologies. There is a systematic shift to lower concentration for WDM halos, but both cosmologies lie on the same locus relating concentration and formation epoch. Because the form of the density profile remains unchanged, our results indicate that the equilibrium halo density profile is set independently of the halo formation process.

preprint2004arXiv

The Ultimate Halo Mass in a LCDM Universe

In the far future of an accelerating LCDM cosmology, the cosmic web of large-scale structure consists of a set of increasingly isolated halos in dynamical equilibrium. We examine the approach of collisionless dark matter to hydrostatic equilibrium using a large N-body simulation evolved to scale factor a = 100, well beyond the vacuum--matter equality epoch, a_eq ~ 0.75, and 53/h Gyr into the future for a concordance model universe (Omega_m ~ 0.3, Omega_Lambda ~ 0.7). The radial phase-space structure of halos -- characterized at a < a_eq by a pair of zero-velocity surfaces that bracket a dynamically active accretion region -- simplifies at a > 10 a_eq when these surfaces merge to create a single zero-velocity surface, clearly defining the halo outer boundary, rhalo, and its enclosed mass, mhalo. This boundary approaches a fixed physical size encompassing a mean interior density ~ 5 times the critical density, similar to the turnaround value in a classical Einstein-deSitter model. We relate mhalo to other scales currently used to define halo mass (m200, mvir, m180b) and find that m200 is approximately half of the total asymptotic cluster mass, while m180b follows the evolution of the inner zero velocity surface for a < 2 but becomes much larger than the total bound mass for a > 3. The radial density profile of all bound halo material is well fit by a truncated Hernquist profile. An NFW profile provides a somewhat better fit interior to r200 but is much too shallow in the range r200 < r < rhalo.

preprint2003arXiv

Future Evolution of Structure in an Accelerating Universe

Current cosmological data indicate that our universe contains a substantial component of dark vacuum energy that is driving the cosmos to accelerate. We examine the immediate and longer term consequences of this dark energy (assumed here to have a constant density). Using analytic calculations and supporting numerical simulations, we present criteria for test bodies to remain bound to existing structures. We show that collapsed halos become spatially isolated and dynamically relax to a particular density profile with logarithmic slope steeper than -3 at radii beyond r_200. The asymptotic form of the space-time metric is then specified. We develop this scenario further by determining the effects of the accelerating expansion on the background radiation fields and individual particles. In an appendix, we generalize these results to include quintessence.