Source author record

S. Michael Fall

S. Michael Fall 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

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

24 published item(s)

preprint2022arXiv

Galactic angular momentum in the IllustrisTNG simulation -- I. Connection to morphology, halo spin, and black hole mass

We use the TNG100 simulation of the IllustrisTNG project to investigate the stellar specific angular momenta ($j_{\ast}$) of $\sim$12,000 central galaxies at $z=0$ in a full cosmological context, with stellar masses ($M_{\ast}$) ranging from $10^{9}$ to $10^{12} \, {\rm M}_{\odot}$. We find that the $j_{\ast}$-$M_{\ast}$ relations for early-type and late-type galaxies in IllustrisTNG are in good overall agreement with observations, and that these galaxy types typically `retain' $\sim$10-20 and $\sim$50-60 per cent of their host haloes' specific angular momenta, respectively, with some dependence on the methodology used to measure galaxy morphology. We present results for kinematic as well as visual-like morphological measurements of the simulated galaxies. Next, we explore the scatter in the $j_{\ast}$-$M_{\ast}$ relation with respect to the spin of the dark matter halo and the mass of the supermassive black hole (BH) at the galactic centre. We find that galaxies residing in faster spinning haloes, as well as those hosting less massive BHs, tend to have a higher specific angular momentum. We also find that, at fixed galaxy or halo mass, halo spin and BH mass are anticorrelated with each other, probably as a consequence of more efficient gas flow toward the galactic centre in slowly rotating systems. Finally, we show that halo spin plays an important role in determining galaxy sizes - larger discs form at the centres of faster-rotating haloes - although the trend breaks down for massive galaxies with $M_{\ast} \gtrsim 10^{11} \, {\rm M}_{\odot}$, roughly the mass scale at which a galaxy's stellar mass becomes dominated by accreted stars.

preprint2020arXiv

Mass Functions of Giant Molecular Clouds and Young Star Clusters in Six Nearby Galaxies

We compare the mass functions of young star clusters (ages $\leq 10$ Myr) and giant molecular clouds (GMCs) in six galaxies that cover a large range in mass, metallicity, and star formation rate (LMC, M83, M51, NGC 3627, the Antennae, and NGC 3256). We perform maximum-likelihood fits of the Schechter function, $ψ(M) = dN/dM \propto M^β \exp(-M/M_*)$, to both populations. We find that most of the GMC and cluster mass functions in our sample are consistent with a pure power-law distribution ($M_* \rightarrow \infty$). M51 is the only galaxy that shows some evidence for an upper cutoff ($M_*$) in both populations. Therefore, physical upper mass cutoffs in populations of both GMCs and clusters may be the exception rather than the rule. When we perform power-law fits, we find a range of indices $β_{\rm PL}=-2.3\pm0.3$ for our GMC sample and $β_{\rm PL}=-2.0\pm0.3$ for the cluster sample. This result, that $β_{\rm Clusters} \approx β_{\rm GMC} \approx -2$, is consistent with theoretical predictions for cluster formation and suggests that the star-formation efficiency is largely independent of mass in the GMCs.

preprint2016arXiv

Energy budget of forming clumps in numerical simulations of collapsing clouds

We analyze the physical properties and energy balance of density enhancements in two SPH simulations of the formation, evolution, and collapse of giant molecular clouds. In the simulations, no feedback is included, so all motions are due either to the initial, decaying turbulence, or to gravitational contraction. We define clumps as connected regions above a series of density thresholds. The resultingfull set of clumps follows the generalized energy-equipartition relation $σ_{v}/R^{1/2} \propto Σ^{1/2}$, where $σ_{v}$ is the velocity dispersion, $R$ is the "radius", and $Σ$ is the column density. We interpret this as a natural consequence of gravitational contraction at all scales, rather than virial equilibrium. Nevertheless, clumps with low $Σ$ tend to show a large scatter around equipartition. In more than half of the cases, this scatter is dominated by external turbulent compressions that {\it assemble} the clumps, rather than by small-scale random motions that would disperse them. The other half does actually disperse. Moreover, clump sub-samples selected by means of different criteria exhibit different scalings. Sub-samples with narrow $Σ$ ranges follow Larson-like relations, although characterized by their respective value of $Σ$. Finally, we find that: i) clumps lying in filaments tend to appear sub-virial; ii) high-density cores ($n \ge 10^5$ cm$^3$) that exhibit moderate kinetic energy excesses often contain sink ("stellar") particles, and the excess disappears when the stellar mass is taken into account in the energy balance; iii) cores with kinetic energy excess but no stellar particles are truly in a state of dispersal.

preprint2016arXiv

Evolution of the Mass and Luminosity Functions of Globular Star Clusters

We reexamine the dynamical evolution of the mass and luminosity functions of globular star clusters (GCMF and GCLF). Fall & Zhang (2001, hereafter FZ01) showed that a power-law MF, as commonly seen among young cluster systems, would evolve by dynamical processes over a Hubble time into a peaked MF with a shape very similar to the observed GCMF in the Milky Way and other galaxies. To simplify the calculations, the semi-analytical FZ01 model adopted the "classical" theory of stellar escape from clusters, and neglected variations in the $M/L$ ratios of clusters. Kruijssen & Portegies Zwart (2009, hereafter KPZ09) modified the FZ01 model to include "retarded" and mass-dependent stellar escape, the latter causing significant $M/L$ variations. KPZ09 asserted that their model was compatible with observations whereas the FZ01 model was not. We show here that this claim is not correct; the FZ01 and KPZ09 models fit the observed Galactic GCLF equally well. We also show that there is no detectable correlation between $M/L$ and $L$ for GCs in the Milky Way and Andromeda galaxies, in contradiction with the KPZ09 model. Our comparisons of the FZ01 and KPZ09 models with observations can be explained most simply if stars escape at rates approaching the classical limit for high-mass clusters, as expected on theoretical grounds.

preprint2016arXiv

Version 1 of the Hubble Source Catalog

The Hubble Source Catalog is designed to help optimize science from the Hubble Space Telescope by combining the tens of thousands of visit-based source lists in the Hubble Legacy Archive into a single master catalog. Version 1 of the Hubble Source Catalog includes WFPC2, ACS/WFC, WFC3/UVIS, and WFC3/IR photometric data generated using SExtractor software to produce the individual source lists. The catalog includes roughly 80 million detections of 30 million objects involving 112 different detector/filter combinations, and about 160 thousand HST exposures. Source lists from Data Release 8 of the Hubble Legacy Archive are matched using an algorithm developed by Budavari & Lubow (2012). The mean photometric accuracy for the catalog as a whole is better than 0.10 mag, with relative accuracy as good as 0.02 mag in certain circumstances (e.g., bright isolated stars). The relative astrometric residuals are typically within 10 mas, with a value for the mode (i.e., most common value) of 2.3 mas. The absolute astrometric accuracy is better than $\sim$0.1 arcsec for most sources, but can be much larger for a fraction of fields that could not be matched to the PanSTARRS, SDSS, or 2MASS reference systems. In this paper we describe the database design with emphasis on those aspects that enable the users to fully exploit the catalog while avoiding common misunderstandings and potential pitfalls. We provide usage examples to illustrate some of the science capabilities and data quality characteristics, and briefly discuss plans for future improvements to the Hubble Source Catalog.

preprint2015arXiv

Galactic Angular Momentum in the Illustris Simulation: Feedback and the Hubble Sequence

We study the stellar angular momentum of thousands of galaxies in the Illustris cosmological simulation, which captures gravitational and gas dynamics within galaxies, as well as feedback from stars and black holes. We find that the angular momentum of the simulated galaxies matches observations well, and in particular two distinct relations are found for late-type versus early-type galaxies. The relation for late-type galaxies corresponds to the value expected from full conservation of the specific angular momentum generated by cosmological tidal torques. The relation for early-type galaxies corresponds to retention of only ~30% of that, but we find that those early-type galaxies with low angular momentum at z=0 nevertheless reside at high redshift on the late-type relation. Some of them abruptly lose angular momentum during major mergers. To gain further insight, we explore the scaling relations in simulations where the galaxy formation physics is modified with respect to the fiducial model. We find that galactic winds with high mass-loading factors are essential for obtaining the high angular momentum relation typical for late-type galaxies, while AGN feedback largely operates in the opposite direction. Hence, feedback controls the stellar angular momentum of galaxies, and appears to be instrumental for establishing the Hubble sequence.

preprint2015arXiv

Gravity or turbulence? -III. Evidence of pure thermal Jeans fragmentation at ~0.1 pc scale

We combine previously published interferometric and single-dish data of relatively nearby massive dense cores that are actively forming stars to test whether their `fragmentation level' is controlled by turbulent or thermal support. We find no clear correlation between the fragmentation level and velocity dispersion, nor between the observed number of fragments and the number of fragments expected when the gravitationally unstable mass is calculated including various prescriptions for `turbulent support'. On the other hand, the best correlation is found for the case of pure thermal Jeans fragmentation, for which we infer a core formation efficiency around 13 per cent, consistent with previous works. We conclude that the dominant factor determining the fragmentation level of star-forming massive dense cores at 0.1 pc scale seems to be thermal Jeans fragmentation.

preprint2015arXiv

Simulating Deep Hubble Images With Semi-empirical Models of Galaxy Formation

We simulate deep images from the Hubble Space Telescope (HST) using semi-empirical models of galaxy formation with only a few basic assumptions and parameters. We project our simulations all the way to the observational domain, adding cosmological and instrumental effects to the images, and analyze them in the same way as real HST images ("forward modeling"). This is a powerful tool for testing and comparing galaxy evolution models, since it allows us to make unbiased comparisons between the predicted and observed distributions of galaxy properties, while automatically taking into account all relevant selection effects. Our semi-empirical models populate each dark matter halo with a galaxy of determined stellar mass and scale radius. We compute the luminosity and spectrum of each simulated galaxy from its evolving stellar mass using stellar population synthesis models. We calculate the intrinsic scatter in the stellar mass-halo mass relation that naturally results from enforcing a monotonically increasing stellar mass along the merger history of each halo. The simulated galaxy images are drawn from cutouts of real galaxies from the Sloan Digital Sky Survey, with sizes and fluxes rescaled to match those of the model galaxies. The distributions of galaxy luminosities, sizes, and surface brightnesses depend on the adjustable parameters in the models, and they agree well with observations for reasonable values of those parameters. Measured galaxy magnitudes and sizes have significant magnitude-dependent biases, with both being underestimated near the magnitude detection limit. The fraction of galaxies detected and fraction of light detected also depend sensitively on the details of the model.

preprint2015arXiv

The Link Between the Formation Rates of Clusters and Stars in Galaxies

The goal of this paper is to test whether the formation rate of star clusters is proportional to the star formation rate (SFR) in galaxies. As a first step, we present the mass functions of compact clusters younger than 10 Myr in seven star-forming galaxies of diverse masses, sizes, and morphologies: the Large and Small Magellanic Clouds, NGC 4214, NGC 4449, M83, M51, and the Antennae. These cluster mass functions (CMFs) are well represented by power laws, dN/dM~M^b, with similar exponents b=-1.92+/-0.27, but with amplitudes that differ by factors up to ~10^3, corresponding to vast differences in the sizes of the cluster populations in these galaxies. We then normalize these CMFs by the SFRs in the galaxies, derived from dust-corrected H-alpha luminosities, and find that the spread in the amplitudes collapses, with a remaining rms deviation of only sigma_(logA)= 0.2. This is close to the expected dispersion from random uncertainties in the CMFs and SFRs. Thus, the data presented here are consistent with exact proportionality between the formation rates of stars and clusters. However, the data also permit weak deviations from proportionality, at the factor of two level, within the statistical uncertainties. We find the same spread in amplitudes when we normalize the mass functions of much older clusters, with ages in the range 100 to 400 Myr, by the current SFR. This is another indication of the general similarity among the cluster populations of different galaxies.

preprint2013arXiv

Angular Momentum and Galaxy Formation Revisited: Effects of Variable Mass-to-light Ratios

We rederive the relation between the specific angular momentum j_* and the mass M_* of the stellar matter in galaxies of different morphological types. This is a revision of the j_*--M_* diagram presented in our recent comprehensive study of galactic angular momentum. In that work, we estimated j_* from kinematic and photometric data that extended to large radii and M_* from near-infrared luminosities L_K with an assumed universal mass-to-light ratio M_*/L_K. However, recent stellar population models show large variations in M_*/L_K correlated with B-V color. In the present work, we use this correlation to estimate M_*/L_K and hence M_* from the measured B-V and L_K. Our revised j_*--M_* diagram is similar to our previous one; both disk-dominated and elliptical galaxies follow nearly parallel sequences with j_* \propto M_*^α and α= 0.6 +/- 0.1. However, the offset between the sequences is now a factor of about 5, some 30% larger than before (and close to the offset found by Fall in 1983). Thus, our new results place even tighter constraints on the loss of specific angular momentum by galactic disks over their lifetimes.

preprint2012arXiv

Angular momentum and galaxy formation revisited

Motivated by new kinematic data in the outer parts of early-type galaxies (ETGs), we re-examine angular momentum (AM) in all galaxy types. We present methods for estimating the specific AM j, focusing on ETGs, to derive relations between stellar j_* and mass M_* (after Fall 1983). We perform analyses of 8 galaxies out to ~10 R_e, finding that data at 2 R_e are sufficient to estimate total j_*. Our results contravene suggestions that ellipticals (Es) harbor large reservoirs of hidden j_* from AM transport in major mergers. We carry out a j_*-M_* analysis of literature data for ~100 nearby bright galaxies of all types. The Es and spirals form parallel j_*-M_* tracks, which for spirals is like the Tully-Fisher relation, but for Es derives from a mass-size-rotation conspiracy. The Es contain ~3-4 times less AM than equal-mass spirals. We decompose the spirals into disks+bulges and find similar j_*-M_* trends to spirals and Es overall. The S0s are intermediate, and we propose that morphological types reflect disk/bulge subcomponents following separate j_*-M_* scaling relations -- providing a physical motivation for characterizing galaxies by mass and bulge/disk ratio. Next, we construct idealized cosmological models of AM content, using a priori estimates of dark matter halo spin and mass. We find that the scatter in halo spin cannot explain the spiral/E j_* differences, but the data are matched if the galaxies retained different fractions of initial j (~60% and ~10%). We consider physical mechanisms for j_* and M_* evolution (outflows, stripping, collapse bias, merging), emphasizing that the vector sum of such processes must produce the observed j_*-M_* relations. A combination of early collapse and multiple mergers (major/minor) may account for the trend for Es. More generally, the observed AM variations represent fundamental constraints for any galaxy formation model.

preprint2012arXiv

Similarities in Populations of Star Clusters

We compare the observed mass functions and age distributions of star clusters in six well-studied galaxies: the Milky Way, Magellanic Clouds, M83, M51, and Antennae. In combination, these distributions span wide ranges of mass and age: $10^2\lea M/M_{\odot}\lea10^6$ and $10^6\leaτ/yr \lea10^9$. We confirm that the distributions are well represented by power laws: $dN/dM\propto M^β$ with $β\approx-1.9$ and $dN/dτ\proptoτ^γ$ with $γ\approx -0.8$. The mass and age distributions are approximately independent of each other, ruling out simple models of mass-dependent disruption. As expected, there are minor differences among the exponents, at a level close to the true uncertainties, $ε_β\simε_γ\sim$~0.1--0.2. However, the overwhelming impression is the similarity of the mass functions and age distributions of clusters in these different galaxies, including giant and dwarf, quiescent and interacting galaxies. This is an important empirical result, justifying terms such as "universal" or "quasi-universal." We provide a partial theoretical explanation for these observations in terms of physical processes operating during the formation and disruption of the clusters, including star formation and feedback, subsequent stellar mass loss, and tidal interactions with passing molecular clouds. A full explanation will require additional information about the molecular clumps and star clusters in galaxies beyond the Milky Way.

preprint2012arXiv

The Radius of Baryonic Collapse in Disc Galaxy Formation

In the standard picture of disc galaxy formation, baryons and dark matter receive the same tidal torques, and therefore approximately the same initial specific angular momentum. However, observations indicate that disc galaxies typically have only about half as much specific angular momentum as their dark matter haloes. We argue this does not necessarily imply that baryons lose this much specific angular momentum as they form galaxies. It may instead indicate that galaxies are most directly related to the inner regions of their host haloes, as may be expected in a scenario where baryons in the inner parts of haloes collapse first. A limiting case is examined under the idealised assumption of perfect angular momentum conservation. Namely, we determine the density contrast Delta, with respect to the critical density of the Universe, by which dark matter haloes need to be defined in order to have the same average specific angular momentum as the galaxies they host. Under the assumption that galaxies are related to haloes via their characteristic rotation velocities, the necessary Delta is ~600. This Delta corresponds to an average halo radius and mass which are ~60% and ~75%, respectively, of the virial values (i.e., for Delta = 200). We refer to this radius as the radius of baryonic collapse R_BC, since if specific angular momentum is conserved perfectly, baryons would come from within it. It is not likely a simple step function due to the complex gastrophysics involved, therefore we regard it as an effective radius. In summary, the difference between the predicted initial and the observed final specific angular momentum of galaxies, which is conventionally attributed solely to angular momentum loss, can more naturally be explained by a preference for collapse of baryons within R_BC, with possibly some later angular momentum transfer.

preprint2011arXiv

Disruption of Star Clusters in the Interacting Antennae Galaxies

We reexamine the age distribution of star clusters in the Antennae in the context of N-body+hydrodynamical simulations of these interacting galaxies. All of the simulations that account for the observed morphology and other properties of the Antennae have star formation rates that vary relatively slowly with time, by factors of only 1.3 - 2.5 in the past 10^8 yr. In contrast, the observed age distribution of the clusters declines approximately as a power law, dN/dt \propto t^{gamma} with gamma = -1.0, for ages 10^6 yr \la t \la 10^9 yr. These two facts can only be reconciled if the clusters are disrupted progressively for at least 10^8 yr and possibly 10^9 yr. When we combine the simulated formation rates with a power-law model, f_surv \propto t^{delta}, for the fraction of clusters that survive to each age t, we match the observed age distribution with exponents in the range -0.9 \la delta \la -0.6 (with a slightly different delta for each simulation). The similarity between delta and gamma indicates that dN/dt is shaped mainly by the disruption of clusters rather than variations in their formation rate. Thus, the situation in the interacting Antennae resembles that in relatively quiescent galaxies such as the Milky Way and the Magellanic Clouds.

preprint2010arXiv

A Comparison of Methods for Determining the Age Distribution of Star Clusters: Application to the Large Magellanic Cloud

The age distribution of star clusters in nearby galaxies plays a crucial role in evaluating the lifetimes and disruption mechanisms of the clusters. Two very different results have been found recently for the age distribution chi(t) of clusters in the Large Magellanic Cloud (LMC). We found that chi(t) can be described approximately by a power law chi(t) propto t^{gamma}, with gamma -0.8, by counting clusters in the mass-age plane, i.e., by constructing chi(t) directly from mass-limited samples. Gieles & Bastian inferred a value of gamma~, based on the slope of the relation between the maximum mass of clusters in equal intervals of log t, hereafter the M_max method, an indirect technique that requires additional assumptions about the upper end of the mass function. However, our own analysis shows that the M_max method gives a result consistent with our direct counting method for clusters in the LMC, namely chi(t) propto t^-0.8 for t<10^9 yr. The reason for the apparent discrepancy is that our analysis includes many massive (M>1.5x10^3 M_sol), recently formed (t<10^7 yr) clusters, which are known to exist in the LMC, whereas Gieles & Bastian are missing such clusters. We compile recent results from the literature showing that the age distribution of young star clusters in more than a dozen galaxies, including dwarf and giant galaxies, isolated and interacting galaxies, irregular and spiral galaxies, has a similar declining shape. We interpret this approximately "universal" shape as due primarily to the progressive disruption of star clusters over their first ~few x 10^8 yr, starting soon after formation, and discuss some observational and physical implications of this early disruption for stellar populations in galaxies.

preprint2010arXiv

New Tests for Disruption Mechanisms of Star Clusters: The Large and Small Magellanic Clouds

We compare the observed bivariate distribution of masses(M) and ages(t) of star clusters in the LMC with the predicted distributions g(M,t) from 3 idealized models for the disruption of star clusters: (1)sudden mass-dependent disruption;(2)gradual mass-dependent disruption; and (3)gradual mass-independent disruption. The model with mass-{\em in}dependent disruption provides a good, first-order description of these cluster populations, with g(M,t) propto M^{beta} t^{gamma}, beta=-1.8+/-0.2 and gamma=-0.8+/-0.2, at least for clusters with ages t<10^9 yr and masses M<10^3 M_sol (more specifically, t<10^7(M/10^2 M_sol)^{1.3} yr). This model predicts that the clusters should have a power-law luminosity function, dN/dL propto L^-1.8, in agreement with observations. The first two models, on the other hand, fare poorly when describing the observations, refuting previous claims that mass-dependent disruption of star clusters is observed in the LMC over the studied M-t domain. Clusters in the SMC can be described by the same g(M,t) distribution as for the LMC, but with smaller samples and hence larger uncertainties. The successful g(M,t) model for clusters in the Magellanic Clouds is virtually the same as the one for clusters in the merging Antennae galaxies, but extends the domain of validity to lower masses and to older ages. This indicates that the dominant disruption processes are similar in these very different galaxies over at least t<10^8 yr and possibly t<10^9 yr. The mass functions for young clusters in the LMC are power-laws, while that for ancient globular clusters is peaked. We show that the observed shapes of these mass functions are consistent with expectations from the simple evaporation model presented by McLaughlin & Fall.

preprint2010arXiv

Stellar Feedback in Molecular Clouds and its Influence on the Mass Function of Young Star Clusters

We investigate how the removal of interstellar material by stellar feedback limits the efficiency of star formation in molecular clouds and how this determines the shape of the mass function of young star clusters. In particular, we derive relations between the power-law exponents of the mass functions of the clouds and clusters in the limiting regimes in which the feedback is energy-driven and momentum-driven, corresponding to minimum and maximum radiative losses and likely to bracket all realistic cases. We find good agreement between the predicted and observed exponents, especially for momentum-driven feedback, provided the protoclusters have roughly constant mean surface density, as indicated by observations of the star-forming clumps within molecular clouds. We also consider a variety of specific feedback mechanisms, concluding that HII regions inflated by radiation pressure predominate in massive protoclusters, a momentum-limited process when photons can escape after only a few interactions with dust grains. We then present a first estimate of the star formation efficiency in this case, finding that it depends on the masses and sizes of the protoclusters only through their mean surface density, thus ensuring consistency between the observed power-law exponents of the mass functions of the clouds and clusters. The numerical value of this efficiency is also consistent with observations.

preprint2009arXiv

Large Area Survey for z=7 Galaxies in SDF and GOODS-N: Implications for Galaxy Formation and Cosmic Reionization

We present results of our large-area survey for z'-band dropout galaxies at z=7 in a 1568 arcmin^2 sky area covering the SDF and GOODS-N fields. Combining our ultra-deep Subaru/Suprime-Cam z'- and y-band (lambda_eff=1um) images with legacy data of Subaru and HST, we have identified 22 bright z-dropout galaxies down to y=26, one of which has a spectroscopic redshift of z=6.96 determined from Lya emission. The z=7 luminosity function (LF) yields the best-fit Schechter parameters of phi*=0.69 +2.62/-0.55 x10^(-3) Mpc^(-3), Muv*=-20.10 +/-0.76 mag, and alpha=-1.72 +/-0.65, and indicates a decrease from z=6 at a >95% confidence level. This decrease is beyond the cosmic variance in our two fields, which is estimated to be a factor of <~2. We have found that the cosmic star formation rate density drops from the peak at z=2-3 to z=7 roughly by a factor of ~10 but not larger than ~100. A comparison with the reionization models suggests either that the Universe could not be totally ionized by only galaxies at z=7, or more likely that properties of galaxies at z=7 are different from those at low redshifts having, e.g., a larger escape fraction (>~0.2), and/or a flatter IMF. Our SDF z-dropout galaxies appear to form 60-Mpc long filamentary structures, and the z=6.96 galaxy with Lya emission is located at the center of an overdense region consisting of four UV bright dropout candidates, which might suggest an existence of a well-developed ionized bubble at z=7.

preprint2009arXiv

New Tests for Disruption Mechanisms of Star Clusters: Methods and Application to the Antennae Galaxies

We present new tests for disruption mechanisms of star clusters based on the bivariate mass-age distribution g(M,τ). In particular, we derive formulae for g(M,τ) for two idealized models in which the rate of disruption depends on the masses of the clusters and one in which it does not. We then compare these models with our Hubble Space Telescope observations of star clusters in the Antennae galaxies over the mass-age domain in which we can readily distinguish clusters from individual stars: τ\la10^7(M/10^4 M_{\odot})^{1.3} yr. We find that the models with mass-dependent disruption are poor fits to the data, even with complete freedom to adjust several parameters, while the model with mass-independent disruption is a good fit. The successful model has the simple form g(M,τ) \propto M^{-2} τ^{-1}, with power-law mass and age distributions, dN/dM propto M^{-2} and dN/dτ\proptoτ^{-1}. The predicted luminosity function is also a power law, dN/dL \propto L^{-2}, in good agreement with our observations of the Antennae clusters. The similarity of the mass functions of star clusters and molecular clouds indicates that the efficiency of star formation in the clouds is roughly independent of their masses. The age distribution of the massive young clusters is plausibly explained by the following combination of disruption mechanisms: (1) removal of interstellar material by stellar feedback, τ\la 10^7$ yr; (2) continued stellar mass loss, 10^7 yr \la τ\la 10^8 yr; (3), tidal disturbances by passing molecular clouds, τ\ga 10^8 yr. None of these processes is expected to have a strong dependence on mass, consistent with our observations of the Antennae clusters. We speculate that this simple picture also applies--at least approximately--to the clusters in many other galaxies.

preprint2008arXiv

Shaping the Globular Cluster Mass Function by Stellar-Dynamical Evaporation

We show that the globular cluster mass function (GCMF) in the Milky Way depends on cluster half-mass density (rho_h) in the sense that the turnover mass M_TO increases with rho_h while the width of the GCMF decreases. We argue that this is the expected signature of the slow erosion of a mass function that initially rose towards low masses, predominantly through cluster evaporation driven by internal two-body relaxation. We find excellent agreement between the observed GCMF -- including its dependence on internal density rho_h, central concentration c, and Galactocentric distance r_gc -- and a simple model in which the relaxation-driven mass-loss rates of clusters are approximated by -dM/dt = mu_ev ~ rho_h^{1/2}. In particular, we recover the well-known insensitivity of M_TO to r_gc. This feature does not derive from a literal ``universality'' of the GCMF turnover mass, but rather from a significant variation of M_TO with rho_h -- the expected outcome of relaxation-driven cluster disruption -- plus significant scatter in rho_h as a function of r_gc. Our conclusions are the same if the evaporation rates are assumed to depend instead on the mean volume or surface densities of clusters inside their tidal radii, as mu_ev ~ rho_t^{1/2} or mu_ev ~ Sigma_t^{3/4} -- alternative prescriptions that are physically motivated but involve cluster properties (rho_t and Sigma_t) that are not as well defined or as readily observable as rho_h. In all cases, the normalization of mu_ev required to fit the GCMF implies cluster lifetimes that are within the range of standard values (although falling towards the low end of this range). Our analysis does not depend on any assumptions or information about velocity anisotropy in the globular cluster system.

preprint2006arXiv

Connection between the Age Distributions of Star Clusters and Field Stars: A First Application to the Small Magellanic Cloud

We present the age distributions for star clusters and individual stars in the Small Magellanic Cloud (SMC) based on data from the Magellanic Clouds Photometric Survey by Zaritsky and collaborators. The age distribution of the SMC clusters shows a steep decline, dN_{cluster}/dt \propto t^{-0.85\pm0.15}, over the period 10^7 < t <10^9 yr. This decline is essentially identical to that observed previously for more massive clusters in the merging Antennae galaxies, and also for lower-mass embedded clusters in the solar neighborhood. The SMC cluster age distribution therefore provides additional evidence for the rapid disruption of star clusters (``infant mortality''). These disrupted clusters deliver their stars to the general field population, implying that the field star age distribution, dN_{fld star}/dt, should have an inverse relation to dN_{cluster}/dt if most stars form initially in clusters. We make specific predictions for dN_{fldstar}/dt based on our cluster disruption models, and compare them with current data available for stars in the SMC. While these data do not extend to sufficiently young ages for a definitive test, they are consistent with a scenario wherein most SMC stars formed in clusters. Future analyses of dN_{fldstar}/dt that extend down to ages of approximately few million years are needed to verify the age relationship between stars residing in clusters and in the field.

preprint2004arXiv

Hubble Space Telescope Observations of Element Abundances in Low-redshift Damped Lyman-alpha Galaxies and Implications for the Global Metallicity-Redshift Relation

Most models of cosmic chemical evolution predict that the mass-weighted mean interstellar metallicity of galaxies should rise with time from a low value $\sim 1/30$ solar at $z \sim 3$ to a nearly solar value at $z = 0$. In the absence of any selection effects, the damped Lyman-alpha absorbers (DLAs) in quasar spectra are expected to show such a rise in global metallicity. However, it has been difficult to determine whether or not DLAs show this effect, primarily because of the very small number of DLA metallicity measurements at low redshifts. In an attempt to put tighter constraints on the low-redshift end of the DLA metallicity-redshift relation, we have observed Zn II and Cr II lines in four DLAs at $0.09 < z < 0.52$, using the Space Telescope Imaging Spectrograph (STIS) onboard the Hubble Space Telescope (HST). These observations have provided the first constraints on Zn abundances in DLAs with $z < 0.4$. In all the three DLAs for which our observations offer meaningful constraints on the metallicity, the data suggest that the metallicities are much lower than the solar value. These results are consistent with recent imaging studies indicating that these DLAs may be associated with dwarf or low surface brightness galaxies. We combine our results with higher redshift data from the literature to estimate the global mean metallicity-redshift relation for DLAs. We find that the global mean metallicity shows at most a slow increase with decreasing redshift. ...(Please see the paper for the complete abstract).

preprint2003arXiv

The Size Evolution of High Redshift Galaxies

Hubble Space Telescope images of high-redshift galaxies selected via color and photometric redshifts are used to examine the size and axial-ratio distribution of galaxies as a function of redshift at lookback times t > 8 Gyr. These parameters are measured at rest-frame UV wavelengths (1200-2000 Angstroms) on images with a rest-frame resolution of less than 0.8 kpc. Galaxy radii are found to scale with redshift approximately as the inverse of the Hubble parameter H(z). This is in accord with the theoretical expectation that the typical sizes of the luminous parts of galaxies should track the expected evolution in the virial radius of dark-matter halos. The mean ratio of semi-major to semi-minor axis for a bright well-resolved sample of galaxies at z ~ 4 is b/a = 0.65, suggesting that these Lyman break galaxies are not drawn from a spheroidal population. However the median concentration index of this sample is C = 3.5, which is closer to the typical concentration indices, C ~ 4, of nearby elliptical galaxies than to the values, C<2 for local disk galaxies of type Sb and later.

preprint1999arXiv

The Mass Function of Young Star Clusters in the "Antennae" Galaxies

We determine the mass function of young star clusters in the merging galaxies the ``Antennae'' (NGC 4038/39) from deep images taken with the WFPC2 on the refurbished Hubble Space Telescope (HST). This is accomplished by means of reddening-free parameters and a comparison with stellar population synthesis tracks to estimate the intrinsic luminosity and age and hence the mass of each cluster. We find that the mass function of the young star clusters (with ages <~ 160 Myr) is well represented by a power law of the form psi(M) ~ M^{-2} over the range 10^4 < M < 10^6 solar masses. This result may have important implications for our understanding of the origin of globular clusters during the early phases of galactic evolution.