Estimating and Analyzing Demographic Models Using the popbio Package in R

A complete assessment of population growth and viability from field census data often requires complex data manipulations, statistical routines, mathematical tools, programming environments, and graphical capabilities. We therefore designed an R package called popbio to facilitate both the construction and analysis of projection matrix models. The package consists primarily of the R translation of MATLAB code found in Caswell (2001) and Morris and Doak (2002) for the analysis of projection matrix models. The package also includes methods to estimate vital rates and construct projection matrix models from census data typically collected in plant demography studies. In these studies, vital rates can often be estimated directly from annual censuses of tagged individuals using transition frequency tables. Because the construction of projection matrix models requires careful management of census data, we describe the steps to construct a projection matrix in detail.

A matter of measurement: rotation velocities and the velocity function of dwarf galaxies

The velocity function derived from large scale surveys can be compared with the predictions of LCDM cosmology, by matching the measured rotation velocities Vrot of galaxies to the maximum circular velocity of dark matter (DM) halos Vmax. For Vrot<50km/s, a major discrepancy arises between the observed and LCDM velocity functions. However, the manner in which different observational measures of Vrot are associated with Vmax is not straight forward in dwarf galaxies. We instead relate galaxies to DM halos using the empirical baryon- mass to halo-mass relation, and show that different observational measures of Vrot result in very different velocity functions. We show how the W50 velocity function, i.e. using the HI profile line width at 50% of peak HI flux to measure Vrot, can be reconciled with a LCDM cosmology. Our semi-empirical methodology allows us to determine the region of rotation curves that are probed by HI measurements (RHI), and shows that the Vrot of dwarfs are generally measured at a fraction of Rmax, explaining their tendency to have rising rotation curves. We provide fitting formulae for relating RHI and Reff (the effective radius) to the virial radius of DM halos. To continue to use velocity functions as a probe of LCDM cosmology, it is necessary to be precise about how the different measures of rotation velocity are probing the mass of the DM halos, dropping the assumption that any measure of rotational velocity can be equally used as a proxy for Vmax.Comment: submitted to MNRAS. Comments very welcom

Galactic Halo Stars in Phase Space :A Hint of Satellite Accretion?

The present day chemical and dynamical properties of the Milky Way bear the imprint of the Galaxy's formation and evolutionary history. One of the most enduring and critical debates surrounding Galactic evolution is that regarding the competition between ``satellite accretion'' and ``monolithic collapse''; the apparent strong correlation between orbital eccentricity and metallicity of halo stars was originally used as supporting evidence for the latter. While modern-day unbiased samples no longer support the claims for a significant correlation, recent evidence has been presented by Chiba & Beers (2000,AJ,119,2843) for the existence of a minor population of high-eccentricity metal-deficient halo stars. It has been suggested that these stars represent the signature of a rapid (if minor) collapse phase in the Galaxy's history. Employing velocity- and integrals of motion-phase space projections of these stars, coupled with a series of N-body/Smoothed Particle Hydrodynamic (SPH) chemodynamical simulations, we suggest an alternative mechanism for creating such stars may be the recent accretion of a polar orbit dwarf galaxy.Comment: 12 pages(incl. figures). Accepted for publication in ApJ letters sectio

The Stellar Halo Metallicity - Luminosity Relationship for Spiral Galaxies

The stellar halos of spiral galaxies bear important chemo-dynamical signatures of galaxy formation. We present here the analysis of 89 semi-cosmological spiral galaxy simulations, spanning ~ 4 magnitudes in total galactic luminosity. These simulations sample a wide variety of merging histories and show significant dispersion in halo metallicity at a given total luminosity - more than a factor of ten in metallicity. Our preliminary analysis suggests that galaxies with a more extended merging history possess halos which have younger and more metal rich stellar populations than the stellar halos associated with galaxies with a more abbreviated assembly. A correlation between halo metallicity and its surface brightness has also been found, reflecting the correlation between halo metallicity and its stellar mass. Our simulations are compared with recent Hubble Space Telescope observations of resolved stellar halos in nearby spirals.Comment: 5 pages, 4 figures. MNRAS Letters, in pres

Signatures of dark matter halo expansion in galaxy populations

This article has been accepted for publication in Monthly Notices of the Royal Astronomical Society ©: 2015 The Authors. Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reservedDark matter cores within galaxy haloes can be formed by energy feedback from star-forming regions: an energy balance suggests that the maximum core formation efficiency arises in galaxies with Mstar ~ 108.5M⊙. We show that a model population of galaxies, in which the density profile has been modified by such baryonic feedback, is able to explain the observed galaxy velocity function and Tully-Fisher relations significantly better than a model in which a universal cuspy density profile is assumed. Alternative models, namely warm or self-interacting dark matter, also provide a better match to these observed relations than a universal profile model does, but make different predictions for how halo density profiles vary with mass compared to the baryonic feedback case. We propose that the expected signatures of the mass dependence of core formation generated by baryonic feedbackCB thanks the MICINN (Spain) for the financial support through the MINECO grant AYA2012-31101 and the Ramon y Cajal program. ADC is supported by the DARK independent fellowship program

Testing Feedback-Modified Dark Matter Haloes with Galaxy Rotation Curves: Estimation of Halo Parameters and Consistency with Λ\LambdaCDM

Cosmological $N$-body simulations predict dark matter (DM) haloes with steep central cusps (e.g. NFW, Navarro et al. 1996). This contradicts observations of gas kinematics in low-mass galaxies that imply the existence of shallow DM cores. Baryonic processes such as adiabatic contraction and gas outflows can, in principle, alter the initial DM density profile, yet their relative contributions to the halo transformation remain uncertain. Recent high resolution, cosmological hydrodynamic simulations (Di Cintio et al. 2014, DC14) predict that inner density profiles depend systematically on the ratio of stellar to DM mass (M$_*$/M$_{\text{halo}}$). Using a Markov Chain Monte Carlo approach, we test the NFW and the M$_*$/M$_{\text{halo}}$-dependent DC14 halo models against a sample of 147 galaxy rotation curves from the new {\it Spitzer} Photometry and Accurate Rotation Curves (SPARC) data set. These galaxies all have extended H{\small I} rotation curves from radio interferometry as well as accurate stellar mass density profiles from near-infrared photometry. The DC14 halo profile provides markedly better fits to the data compared to the NFW profile. Unlike NFW, the DC14 halo parameters found in our rotation curve fits naturally fall within two standard deviations of the mass-concentration relation predicted by $\Lambda$CDM and the stellar mass-halo mass relation inferred from abundance matching with few outliers. Halo profiles modified by baryonic processes are therefore more consistent with expectations from $\Lambda$ cold dark matter ($\Lambda$CDM) cosmology and provide better fits to galaxy rotation curves across a wide range of galaxy properties than do halo models that neglect baryonic physics. Our results offer a solution to the decade long cusp-core discrepancy.Comment: 23 Pages, 18 Figures, MNRAS Accepte