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.

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

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

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

Expanded haloes, abundance matching and too-big-to-fail in the Local Group

Observed kinematical data of 40 Local Group (LG) members are used to derive the dark matter halo mass of such galaxies. Haloes are selected from the theoretically expected LG mass function and two different density profiles are assumed, a standard universal cuspy model and a mass dependent profile which accounts for the effects of baryons in modifying the dark matter distribution within galaxies. The resulting relations between stellar and halo mass are compared with expectations from abundance matching. Using a universal cuspy profile, the ensemble of LG galaxies is fit in relatively low mass haloes, leaving "dark" many massive haloes of \mhalo$\gtrsim$10$^{10}$\msun: this reflects the "too big to fail" problem and results in a \mstar-\mhalo\ relation that differs from abundance matching predictions. Moreover, the star formation efficiency of isolated LG galaxies increases with decreasing halo mass when adopting a cuspy model. By contrast, using the mass dependent density profile, dwarf galaxies with \mstar$\gtrsim$10$^{6}$\msun are assigned to more massive haloes, which have a central cored distribution of dark matter: the "too big to fail" problem is alleviated, the resultant \mstar-\mhalo\ relation follows abundance matching predictions down to the completeness limit of current surveys, and the star formation efficiency of isolated members decreases with decreasing halo mass, in agreement with theoretical expectations. Finally, the cusp/core space of LG galaxies is presented, providing a framework to understand the non-universality of their density profiles.Comment: Accepted in MNRAS, 15 pages 7 figures. section 3.3 has been added after report

The Emergence of the Thick Disk in a CDM Universe II: Colors and Abundance Patterns

The recently emerging conviction that thick disks are prevalent in disk galaxies, and their seemingly ubiquitous old ages, means that the formation of the thick disk, perhaps more than any other component, holds the key to unravelling the evolution of the Milky Way, and indeed all disk galaxies. In Paper I, we proposed that the thick disk was formed in an epoch of gas rich mergers, at high redshift. This hypothesis was based on comparing N-body/SPH simulations to a variety of Galactic and extragalactic observations, including stellar kinematics, ages and chemical properties.Here examine our thick disk formation scenario in light of the most recent observations of extragalactic thick disks. In agreement, our simulted thick disks are old and relatively metal rich, with V-I colors that do not vary significantly with distance from the plane. Further, we show that our proposal results in an enhancement of alpha-elements in thick disk stars as compared with thin disk stars, consistent with observations of the relevant populations of the Milky Way. We also find that our scenario naturally leads to the formation of an old metal weak stellar halo population with high alpha-element abundances.Comment: submitted to Ap

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