No Arabic abstract
We present SimSpin, a new, public, software framework for generating integral field spectroscopy (IFS) data cubes from N-body/hydrodynamical simulations of galaxies, which can be compared directly with observational datasets. SimSpin provides a consistent method for studying a galaxys stellar component. It can be used to explore how observationally inferred measurements of kinematics, such as the spin parameter $lambda_R$, are impacted by the effects of, for example, inclination, seeing conditions, distance, etc. SimSpin is written in R and has been designed to be highly modular, flexible, and extensible. It is already being used by the astrophysics community to generate IFS-like cubes and FITS files for direct comparison of simulations to observations. In this paper, we explain the conceptual framework of SimSpin; how it is implemented in R; and we demonstrate SimSpins current capabilities, providing as an example a brief investigation of how numerical resolution affects how reliably we can recover the intrinsic stellar kinematics of a simulated galaxy.
Observers experience a series of limitations when measuring galaxy kinematics, such as variable seeing conditions and aperture size. These effects can be reduced using empirical corrections, but these equations are usually applicable within a restrictive set of boundary conditions (e.g. Sersic indices within a given range) which can lead to biases when trying to compare measurements made across a full kinematic survey. In this work, we present new corrections for two widely used kinematic parameters, $lambda_R$ and $V/sigma$, that are applicable across a broad range of galaxy shapes, measurement radii and ellipticities. We take a series of mock observations of N-body galaxy models and use these to quantify the relationship between the observed kinematic parameters, structural properties and different seeing conditions. Derived corrections are then tested using the full catalogue of galaxies, including hydro-dynamic models from the EAGLE simulation. Our correction is most effective for regularly-rotating systems, yet the kinematic parameters of all galaxies -- fast, slow and irregularly rotating systems -- are recovered successfully. We find that $lambda_R$ is more easily corrected than $V/sigma$, with relative deviations of 0.02 and 0.06 dex respectively. The relationship between $lambda_R$ and $V/sigma$, as described by the parameter $kappa$, also has a minor dependence on seeing conditions. These corrections will be particularly useful for stellar kinematic measurements in current and future integral field spectroscopic (IFS) surveys of galaxies.
We describe the construction of a suite of galaxy cluster mock catalogues from N-body simulations, based on the properties of the new ROSAT-ESO Flux-Limited X-Ray (REFLEX II) galaxy cluster catalogue. Our procedure is based on the measurements of the cluster abundance, and involves the calibration of the underlying scaling relation linking the mass of dark matter haloes to the cluster X-ray luminosity determined in the emph{ROSAT} energy band $0.1-2.4$ keV. In order to reproduce the observed abundance in the luminosity range probed by the REFLEX II X-ray luminosity function ($0.01<L_{X}/(10^{44}{rm erg},{rm s}^{-1}h^{-2})<10$), a mass-X ray luminosity relation deviating from a simple power law is required. We discuss the dependence of the calibration of this scaling relation on the X-ray luminosity and the definition of halo masses and analyse the one- and two-point statistical properties of the mock catalogues. Our set of mock catalogues provides samples with self-calibrated scaling relations of galaxy clusters together with inherent properties of flux-limited surveys. This makes them a useful tool to explore different systematic effects and statistical methods involved in constraining both astrophysical and cosmological information from present and future galaxy cluster surveys.
We compare the stellar motion around a spiral arm created in two different scenarios, transient/co-rotating spiral arms and density-wave-like spiral arms. We generate Gaia mock data from snapshots of the simulations following these two scenarios using our stellar population code, SNAPDRAGONS, which takes into account dust extinction and the expected Gaia errors. We compare the observed rotation velocity around a spiral arm similar in position to the Perseus arm, and find that there is a clear difference in the velocity features around the spiral arm between the co-rotating spiral arm and the density-wave-like spiral arm. Our result demonstrates that the volume and accuracy of the Gaia data are sufficient to clearly distinguish these two scenarios of the spiral arms.
Integral field spectroscopy can map astronomical objects spatially and spectroscopically. Due to instrumental and atmospheric effects, it is common for integral field instruments to yield a sampling of the sky image that is both irregular and wavelength-dependent. Most subsequent analysis procedures require a regular, wavelength independent sampling (for example a fixed rectangular grid), and thus an initial step of fundamental importance is to resample the data onto a new grid. The best possible resampling would produce a well-sampled image, with a resolution equal to that imposed by the intrinsic spatial resolution of the instrument, telescope, and atmosphere, and with no statistical correlations between neighboring pixels. A standard method in the field to produce a regular set of samples from an irregular set of samples is Shepards method, but Shepards method typically yields images with a degraded resolution and large statistical correlations between pixels. Here we introduce a new method, which improves on Shepards method in both these respects. We apply this method to data from the Mapping Nearby Galaxies at Apache Point Observatory survey, part of Sloan Digital Sky Survey IV, demonstrating a full-width half maximum close to that of the intrinsic resolution (and ~ 16% better than Shepards method) and low statistical correlations between pixels. These results nearly achieve the ideal resampling. This method can have broader applications to other integral field data sets and to other astronomical data sets (such as dithered images) with irregular sampling.
The MaNGA Survey (Mapping Nearby Galaxies at Apache Point Observatory) is one of three core programs in the Sloan Digital Sky Survey IV. It is obtaining integral field spectroscopy (IFS) for 10K nearby galaxies at a spectral resolution of R~2000 from 3,622-10,354A. The design of the survey is driven by a set of science requirements on the precision of estimates of the following properties: star formation rate surface density, gas metallicity, stellar population age, metallicity, and abundance ratio, and their gradients; stellar and gas kinematics; and enclosed gravitational mass as a function of radius. We describe how these science requirements set the depth of the observations and dictate sample selection. The majority of targeted galaxies are selected to ensure uniform spatial coverage in units of effective radius (Re) while maximizing spatial resolution. About 2/3 of the sample is covered out to 1.5Re (Primary sample), and 1/3 of the sample is covered to 2.5Re (Secondary sample). We describe the survey execution with details that would be useful in the design of similar future surveys. We also present statistics on the achieved data quality, specifically, the point spread function, sampling uniformity, spectral resolution, sky subtraction, and flux calibration. For our Primary sample, the median r-band signal-to-noise ratio is ~73 per 1.4A pixel for spectra stacked between 1-1.5 Re. Measurements of various galaxy properties from the first year data show that we are meeting or exceeding the defined requirements for the majority of our science goals.