We develop and implement a model to analyze the internal kinematics of galaxy clusters that may contain subpopulations of galaxies that do not independently trace the cluster potential. The model allows for substructures within the cluster environment, disentangles cluster members from contaminating foreground and background galaxies, and includes an overall cluster rotation term as part of the cluster kinematics. We estimate the cluster velocity dispersion and/or mass while marginalizing over uncertainties in all of the above complexities. In a first application to our published data for Abell 267 (A267), we find no evidence for cluster rotation but we identify up to five distinct galaxy subpopulations. We use these results to explore the sensitivity of inferred cluster properties to the treatment of substructure. Compared to a model that assumes no substructure, our substructure model reduces the dynamical mass of A267 by $\sim 20\%$ and shifts the cluster mean velocity by $\sim 100$ km s$^{-1}$, approximately doubling the offset with respect to the velocity of A267's brightest cluster galaxy. Embedding the spherical Jeans equation within this framework, we infer for A267 a dark matter halo of mass $M_{200}=6.77\pm1.06\times10^{14}M_\odot/h$, concentration $\log_{10}c_{200}=0.61\pm0.39$, consistent with the mass-concentration relation found in cosmological simulations.

We report the results of a pilot program to use the Magellan/M2FS spectrograph to survey the galactic populations and internal kinematics of galaxy clusters. For this initial study, we present spectroscopic measurements for $223$ quiescent galaxies observed along the line of sight to the galaxy cluster Abell 267 ($z\sim0.23$). We develop a Bayesian method for modeling the integrated light from each galaxy as a simple stellar population, with free parameters that specify redshift ($v_\mathrm{los}/c$) and characteristic age, metallicity ($\mathrm{[Fe/H]}$), alpha-abundance ($[\alpha/\mathrm{Fe}]$), and internal velocity dispersion ($\sigma_\mathrm{int}$) for individual galaxies. Parameter estimates derived from our 1.5-hour observation of A267 have median random errors of $\sigma_{v_\mathrm{los}}=20\ \mathrm{km\ s^{-1}}$, $\sigma_{\mathrm{Age}}=1.2\ \mathrm{Gyr}$, $\sigma_{\mathrm{[Fe/H]}}=0.11\ \mathrm{dex}$, $\sigma_{[\alpha/\mathrm{Fe}]}=0.07\ \mathrm{dex}$, and $\sigma_{\sigma_\mathrm{int}}=20\ \mathrm{km\ s^{-1}}$. In a companion paper, we use these results to model the structure and internal kinematics of A267.

Studying the trajectories of objects like stars, globular clusters or satellite galaxies in the Milky Way allows to trace the dark matter halo but requires reliable models of its gravitational potential. Realistic, yet simple and fully analytical models have already been presented in the past. However, improved as well as new observational constraints have become available in the meantime calling for a recalibration of the respective model parameters. Three widely used model potentials are revisited. By a simultaneous least-squared fit to the observed rotation curve, in-plane proper motion of Sgr A*, local mass/surface density and the velocity dispersion in Baade's window, parameters of the potentials are brought up-to-date. The mass at large radii - and thus in particular that of the dark matter halo - is hereby constrained by imposing that the most extreme halo blue horizontal-branch star known has to be bound to the Milky Way. The Galactic mass models are tuned to yield a very good match to recent observations. The mass of the dark matter halo is - within the limitations of the applied models - estimated in a fully consistent way. As a first application, the trajectory of the hypervelocity star HE 0437-5439 is re-investigated to check its suggested origin in the Large Magellanic Cloud (LMC). Despite their simplicity, the presented Milky Way mass models are well able to reproduce all observational constraints. Their analytical and simple form makes them ideally suited for fast and accurate orbit calculations. The LMC cannot be ruled out as HE 0437-5439's birthplace.

(Abridged for arXiv)We investigate the properties of halos, galaxies and blackholes to z=0 in the high resolution hydrodynamical simulation MassiveBlack-II (MBII) which evolves a LCDM cosmology in a comoving volume Vbox=100(Mpc/h)^3. MBII is the highest resolution simulation of this size which includes a self-consistent model for star formation, black hole accretion and associated feedback. We provide a simulation browser web application which enables interactive search and tagging of halos, subhalos and their properties and publicly release our galaxy catalogs. Our analysis of the halo mass function (MF) in MBII reveals that baryons have strong effects, with changes in the halo abundance of 20-35% below the knee of the MF (Mhalo < 10^13.2 Msun/h at z=0) when compared to fits based on dark matter only simulations. We provide a fitting function for the halo MF out to redshift z=11 and discuss how the onset of non-universality in the MF limits the accuracy of our fit. We study the halo occupation distribution and clustering of galaxies, in particular the evolution and scale dependence of stochasticity and bias finding reasonable agreement with observational data. The shape of the cosmic spectral energy distribution predicted by MBII is consistent with observations, but lower in amplitude. The Galaxy Stellar Mass Function (GSMF) function is broadly consistent with observations at z>=2. At z<2, the population of passive low mass (for M*<10^9 Msun) galaxies in MBII makes the GSMF too steep compared to observations whereas at the high mass end (M*>10^11 Msun) galaxies hosting bright AGNs make significant contributions to the GSMF. The quasar bolometric luminosity function is also largely consistent with observations. We note however that more efficient AGN feedback (beyond simple thermal coupling used here) is likely necessary for the largest, rarest objects/clusters at low redshifts.

We present the results of a search for gravitational waves associated with 154 gamma-ray bursts (GRBs) that were detected by satellite-based gamma-ray experiments in 2009-2010, during the sixth LIGO science run and the second and third Virgo science runs. We perform two distinct searches: a modeled search for coalescences of either two neutron stars or a neutron star and black hole; and a search for generic, unmodeled gravitational-wave bursts. We find no evidence for gravitational-wave counterparts, either with any individual GRB in this sample or with the population as a whole. For all GRBs we place lower bounds on the distance to the progenitor, under the optimistic assumption of a gravitational-wave emission energy of 10^-2 M c^2 at 150 Hz, with a median limit of 17 Mpc. For short hard GRBs we place exclusion distances on binary neutron star and neutron star-black hole progenitors, using astrophysically motivated priors on the source parameters, with median values of 16 Mpc and 28 Mpc respectively. These distance limits, while significantly larger than for a search that is not aided by GRB satellite observations, are not large enough to expect a coincidence with a GRB. However, projecting these exclusions to the sensitivities of Advanced LIGO and Virgo, which should begin operation in 2015, we find that the detection of gravitational waves associated with GRBs will become quite possible.