Source author record

Antonaldo Diaferio

Antonaldo Diaferio 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

57works
8topics
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

57 published item(s)

preprint2022arXiv

A Spectroscopic View of the JWST/GTO Strong Lensing Cluster A1489

We discuss a spectroscopic survey of the strong lensing cluster A1489 that includes redshifts for 195 cluster members along with central velocity dispersions for 188 cluster members. The caustic technique applied to the redshift survey gives the dynamical parameters $M_{200} = (1.25~\pm~0.09) \times 10^{15}~M_\odot$, $r_{200} = 1.97~\pm~{0.05}$ Mpc, and a cluster line-of sight velocity dispersion $1150~\pm~{72}~$km$~$s$^{-1}$ within $r_{200}$. These parameters are very similar to those of other strong lensing systems with comparably large Einstein radii. We use the spectroscopy and deep photometry to demonstrate that A1489 is probably dynamically active; its four BCGs have remarkably different rest frame radial velocities. Like other massive strong lensing clusters, the velocity dispersion function for members of A1489 shows an excess for dispersions $\geq~250~$km$~$s$^{-1}$. The central dispersions also provide enhanced constraints on future lensing models.

preprint2022arXiv

Dynamics of dwarf galaxies in $f(R)$ gravity

We use the kinematic data of the stars in eight dwarf spheroidal galaxies to assess whether $f(R)$ gravity can fit the observed profiles of the line-of-sight velocity dispersion of these systems without resorting to dark matter. Our model assumes that each galaxy is spherically symmetric and has a constant velocity anisotropy parameter $β$ and constant mass-to-light ratio consistent with stellar population synthesis models. We solve the spherical Jeans equation that includes the Yukawa-like gravitational potential appearing in the weak field limit of $f(R)$ gravity, and a Plummer density profile for the stellar distribution. The $f(R)$ velocity dispersion profiles depend on two parameters: the scale length $ξ^{-1}$, below which the Yukawa term is negligible, and the boost of the gravitational field $δ>-1$. $δ$ and $ξ$ are not universal parameters, but their variation within the same class of objects is expected to be limited. The $f(R)$ velocity dispersion profiles fit the data with a value $ξ^{-1}= 1.2^{+18.6}_{-0.9}$ Mpc for the entire galaxy sample. On the contrary, the values of $δ$ show a bimodal distribution that picks at $\barδ=-0.986\pm0.002$ and $\barδ=-0.92\pm0.01$. These two values disagree at $6σ$ and suggest a severe tension for $f(R)$ gravity. It remains to be seen whether an improved model of the dwarf galaxies or additional constraints provided by the proper motions of stars measured by future astrometric space missions can return consistent $δ$'s for the entire sample and remove this tension.

preprint2022arXiv

Probing modified Newtonian dynamics with hypervelocity stars

We show that measuring the velocity components of hypervelocity stars (HVSs) can discriminate between modified Newtonian dynamics (MOND) and Newtonian gravity. HVSs are ejected from the Galactic center on radial trajectories with a null tangential velocity component in the reference frame of the Galaxy. They acquire tangential components due to the nonspherical components of the Galactic gravitational potential. Axisymmetric potentials only affect the latitudinal components, $v_θ$, and non-null azimuthal components, $v_ϕ$, originate from non-axisymmetric matter distributions. For HVSs with sufficiently high ejection speed, $v_ϕ$ is proportionate to the deviation of the gravitational potential from axial symmetry. The ejection velocity threshold is $\sim$ 750 km/s for 4 $M_{\odot}$ stars and increases with decreasing HVS mass. We determine the upper limit of $v_ϕ$ as a function of the galactocentric distance for these high-speed HVSs if QUMOND, the quasi-linear formulation of MOND, is the correct theory of gravity and either the triaxial Galactic bulge or a nonspherical hot gaseous halo is the primary source of $v_ϕ$. In Newtonian gravity, the HVSs within 60 kpc of the Galactic center may easily have $v_ϕ$ values higher than the QUMOND upper limit if the dark matter (DM) halo is triaxial or if the DM halo and the baryonic components are axisymmetric but their two axes of symmetry are misaligned. Therefore, even a limited sample of high-speed HVSs could distinguish between QUMOND and the DM model. This test is currently limited by (i) the lack of a proper procedure to assess the HVS nature of a star in the model to be constrained; and (ii) the present uncertainties on $v_ϕ$, which are a factor of $\sim 10$ too large. A proper procedure to assess the HVS nature of the stars and astrometric measurements with microarcsecond precision would make this test feasible.

preprint2022arXiv

Probing the shape of the Milky Way dark matter halo with hypervelocity stars: a new method

We propose a new method to determine the shape of the gravitational potential of the dark matter (DM) halo of the Milky Way (MW) with the galactocentric tangential velocities of a sample of hypervelocity stars (HVSs). We compute the trajectories of different samples of HVSs in a MW where the baryon distribution is axisymmetric and the DM potential either is spherical or is spheroidal or triaxial with radial-dependent axis ratios. We determine the shape of the DM potential with the distribution of the latitudinal velocity $|v_{\vartheta}|$ in axisymmetric Galactic potentials, or with the distribution of $|v_{\vartheta}|$ and of a function $\bar v_φ$ of the azimuthal velocity in non-axisymmetric Galactic potentials. We recover the correct shape of the DM potential by comparing the distribution of $|v_{\vartheta}|$ and $\bar v_φ$ against the corresponding distributions of mock samples of HVSs that traveled in DM halos of different shapes. We use the largest possible sample of $\sim 800$ HVSs of $4~M_\odot$ ejected with the Hills mechanism at a rate $\sim 10^{-4}$ yr$^{-1}$, currently outgoing, and located at more than 10 kpc from the Galactic center. In our ideal case of galactocentric velocities with null uncertainties and no observational limitations, our method recovers the correct shape of the DM potential with a success rate $S\gtrsim 89\%$ in axisymmetric Galactic potentials, and $S > 96\%$ in the explored non-axisymmetric cases. The unsuccessful cases yield axis ratios of the DM potential that are off by $\pm 0.1$. The success rate decreases with decreasing sample size: for example, for a spherical DM halo, $S$ drops from $\sim 98\%$ to $\sim 38\%$ when the sample size decreases from $\sim 800$ to $\sim 40$ HVSs. A robust determination of the shape of the DM potential thus requires the measure of the galactocentric velocity of a few hundred genuine HVSs.

preprint2022arXiv

The dynamics of three nearby E0 galaxies in refracted gravity

We tested whether refracted gravity, a theory of modified gravity that describes the dynamics of galaxies without the aid of dark matter, can model the dynamics of the three massive elliptical galaxies, NGC 1407, NGC 4486, and NGC 5846, out to $\sim$$10R_{\rm e}$, where the baryonic mass component fades out and dark matter is required in Newtonian gravity. We probed these outer regions with the kinematics of the globular clusters provided by the SLUGGS survey. Refracted gravity mimics dark matter with the gravitational permittivity, a monotonic function of the local mass density depending on three parameters, $ε_0$, $ρ_{\rm c}$, and $Q$, which are expected to be universal. Refracted gravity satisfactorily reproduces the velocity dispersion profiles of the stars and red and blue globular clusters, with stellar mass-to-light ratios in agreement with stellar population synthesis models, and orbital anisotropy parameters consistent with previous results obtained in Newtonian gravity with dark matter. The sets of the three parameters of the gravitational permittivity found for each galaxy are consistent with each other within $\sim$2$σ$. We compare the mean $\{ε_0,Q,\log_{10}\left[ρ_{\rm c}\left(\mathrm{g}\, \mathrm{cm}^{-3}\right)\right]\} = \{0.089^{+0.038}_{-0.035}, 0.47^{+0.29}_{-0.21}, -24.25^{+0.28}_{-0.20}\}$ found here with the means of the parameters required to model the rotation curves and vertical velocity dispersion profiles of 30 disk galaxies from the DiskMass Survey (DMS): $ρ_{\rm c}$ and $Q$ agree within 1$σ$ with the DMS values, whereas $ε_0$ agrees within 3$σ$. This agreement suggests that ellipticals and disk galaxies allow for common values of the parameters of the permittivity and supports the universality of the permittivity function.

preprint2022arXiv

The proper motion of stars in dwarf galaxies: distinguishing central density cusps from cores

We show that measuring the proper motion of ${\sim 2000}$ stars within a dwarf galaxy, with an uncertainty of 1 km/s at most, can establish whether the Dark Matter (DM) density profile of the dwarf has a central core or cusp. We derive these limits by building mock star catalogues similar to those expected from future astrometric {\it Theia}-like missions and including celestial coordinates, radial velocity and proper motion of the stars. The density field of the DM halo of the dwarf is sampled from an extended Navarro-Frank-White (eNWF) spherical model, whereas the number density distribution of the stars is a Plummer sphere. The velocity field of the stars is set according to the Jeans equations. A Monte Carlo Markov Chain algorithm applied to a sample of $N\gtrsim 2000$ stars returns unbiased estimates of the eNFW DM parameters within $10\%$ of the true values and with $1σ$ relative uncertainties $\lesssim 20$\%. The proper motions of the stars lift the degeneracy among the eNFW parameters which appears when the line-of-sight velocities alone are available. {Our analysis demonstrates that, by estimating the log-slope of the mass density profile estimated at the half-light radius, a sample of $N=2000$ stars can distinguish between a core and a cusp at more than $8σ$.} Proper motions also return unbiased estimates of the dwarf mass profile with $1σ$ uncertainties that decrease, on average, from 2.65 dex to 0.15 dex when the size of the star sample increases from $N=100$ to $N=6000$ stars. The measure of the proper motions can thus strongly constrain the distribution of DM in nearby dwarfs and provides a fundamental contribution to understanding the nature and the properties of DM.

preprint2020arXiv

Dark matters on the scale of galaxies

The cold dark matter model successfully explains both the emergence and evolution of cosmic structures on large scales and, when we include a cosmological constant, the properties of the homogeneous and isotropic Universe. However, the cold dark matter model faces persistent challenges on the scales of galaxies. {Indeed,} N-body simulations predict some galaxy properties that are at odds with the observations. These discrepancies are primarily related to the dark matter distribution in the innermost regions of the halos of galaxies and to the dynamical properties of dwarf galaxies. They may have three different origins: (1) the baryonic physics affecting galaxy formation is still poorly understood and it is thus not properly included in the model; (2) the actual properties of dark matter differs from those of the conventional cold dark matter; (3) the theory of gravity departs from General Relativity. Solving these discrepancies is a rapidly evolving research field. We illustrate some of the solutions proposed} within the cold dark matter model, and solutions when including warm dark matter, self-interacting dark matter, axion-like particles, or fuzzy dark matter. { We also illustrate some modifications of the theory of gravity: Modified Newtonian Dynamics (MOND), MOdified Gravity (MOG), and $f(R)$ gravity.

preprint2020arXiv

Unveiling the Hierarchical Structure of Open Star Clusters: the Perseus Double Cluster

We introduce a new kinematic method to investigate the structure of open star clusters. We adopt a hierarchical clustering algorithm that uses the celestial coordinates and the proper motions of the stars in the field of view of the cluster to estimate a proxy of the pairwise binding energy of the stars and arrange them in a binary tree. The cluster substructures and their members are identified by trimming the tree at two thresholds, according to the $σ$-plateau method. Testing the algorithm on 100 mock catalogs shows that, on average, the membership of the identified clusters is $(91.5\pm 3.5)$\% complete and the fraction of unrelated stars is $(10.4\pm 2.0)$\%. We apply the algorithm to the stars in the field of view of the Perseus double cluster from the Data Release 2 of Gaia. This approach identifies a single structure, Sub1, that separates into two substructures, Sub1-1 and Sub1-2. These substructures coincide with $h$ Per and $χ$ Per: the distributions of the proper motions and the color-magnitude diagrams of the members of Sub1-1 and Sub1-2 are fully consistent with those of $h$ Per and $χ$ Per reported in the literature. These results suggest that our hierarchical clustering algorithm can be a powerful tool to unveil the complex kinematic information of star clusters.

preprint2019arXiv

Velocity Dispersions of Brightest Cluster Galaxies and Their Host Clusters

We explore connections between brightest cluster galaxies (BCGs) and their host clusters. We first construct a HeCS-omnibus cluster sample including 227 galaxy clusters within $0.02 < z < 0.30$; the total number of spectroscopic members from MMT/Hectospec and SDSS observations is 52325. Taking advantage of the large spectroscopic sample, we compute physical properties of the clusters including the dynamical mass and cluster velocity dispersion ($σ_{cl}$). We also measure the central stellar velocity dispersion of the BCGs ($σ_{*,BCGs}$) to examine the relation between BCG velocity dispersion and cluster velocity dispersion for the first time. The observed relation between BCG velocity dispersion and the cluster velocity dispersion is remarkably tight. Interestingly, the $σ_{*, BCG} / σ_{cl}$ ratio decreases as a function of $σ_{cl}$ unlike the prediction from the numerical simulation of Dolag et al. (2010). The trend in $σ_{*, BCG} / σ_{cl}$ suggests that the BCG formation is more efficient in lower mass halos.

preprint2016arXiv

Compact Groups of Galaxies with Complete Spectroscopic Redshifts in the Local Universe

Dynamical analysis of compact groups provides important tests of models of compact group formation and evolution. By compiling 2066 redshifts from FLWO/FAST, from the literature, and from SDSS DR12 in the fields of compact groups in \citet{McC09}, we construct the largest sample of compact groups with complete spectroscopic redshifts in the redshift range $0.01 < z < 0.22$. This large redshift sample shows that the interloper fraction in the \citet{McC09} compact group candidates is $\sim 42\%$. A secure sample of 332 compact groups includes 192 groups with four or more member galaxies and 140 groups with three members. The fraction of early-type galaxies in these compact groups is 62\%, slightly higher than for the original Hickson compact groups. The velocity dispersions of early- and late-type galaxies in compact groups change little with groupcentric radius; the radii sampled are less than $100 ~h^{-1}$ kpc, smaller than the radii typically sampled by members of massive clusters of galaxies. The physical properties of our sample compact groups include size, number density, velocity dispersion, and local environment; these properties slightly differ from those derived for the original Hickson compact groups and for the DPOSS II compact groups. Differences result from subtle differences in the way the group candidates were originally selected. The space density of the compact groups changes little with redshift over the range covered by this sample. The approximate constancy of the space density for this sample is a potential constraint on the evolution of compact groups on a few Gigayear timescale.

preprint2016arXiv

Dynamics of galaxies and clusters in \textit{refracted gravity}

We investigate the proof of concept and the implications of \textit{refracted gravity}, a novel modified gravity aimed to solve the discrepancy between the luminous and the dynamical mass of cosmic structures without resorting to dark matter. Inspired by the behavior of electric fields in matter, refracted gravity introduces a gravitational permittivity that depends on the local mass density and modifies the standard Poisson equation. The resulting gravitational field can become more intense than the Newtonian field and can mimic the presence of dark matter. We show that the refracted gravitational field correctly describes (1) the rotation curves and the Tully-Fisher relation of disk galaxies; and (2) the observed temperature profile of the X-ray gas of galaxy clusters. According to these promising results, we conclude that refracted gravity deserves further investigation.

preprint2016arXiv

HectoMAP and Horizon Run 4: Dense Structures and Voids in the Real and Simulated Universe

HectoMAP is a dense redshift survey of red galaxies covering a 53 $deg^{2}$ strip of the northern sky. HectoMAP is 97\% complete for galaxies with $r<20.5$, $(g-r)>1.0$, and $(r-i)>0.5$. The survey enables tests of the physical properties of large-scale structure at intermediate redshift against cosmological models. We use the Horizon Run 4, one of the densest and largest cosmological simulations based on the standard $Λ$ Cold Dark Matter ($Λ$CDM) model, to compare the physical properties of observed large-scale structures with simulated ones in a volume-limited sample covering 8$\times10^6$ $h^{-3}$ Mpc$^3$ in the redshift range $0.22<z<0.44$. We apply the same criteria to the observations and simulations to identify over- and under-dense large-scale features of the galaxy distribution. The richness and size distributions of observed over-dense structures agree well with the simulated ones. Observations and simulations also agree for the volume and size distributions of under-dense structures, voids. The properties of the largest over-dense structure and the largest void in HectoMAP are well within the distributions for the largest structures drawn from 300 Horizon Run 4 mock surveys. Overall the size, richness and volume distributions of observed large-scale structures in the redshift range $0.22<z<0.44$ are remarkably consistent with predictions of the standard $Λ$CDM model.

preprint2016arXiv

Hydrostatic and Caustic Mass Profiles of Galaxy Clusters

We compare X-ray and caustic mass profiles for a sample of 16 massive galaxy clusters. We assume hydrostatic equilibrium in interpreting the X-ray data, and use large samples of cluster members with redshifts as a basis for applying the caustic technique. The hydrostatic and caustic masses agree to better than $\approx20\%$ on average across the radial range covered by both techniques $(\sim[0.2-1.25]R_{500})$. The mass profiles were measured independently and do not assume a common functional form. Previous studies suggest that, at $R_{500}$, the hydrostatic and caustic masses are biased low and high respectively. We find that the ratio of hydrostatic to caustic mass at $R_{500}$ is $1.20^{+0.13}_{-0.11}$; thus it is larger than 0.9 at $\approx3σ$ and the combination of under- and over-estimation of the mass by these two techniques is $\approx10\%$ at most. There is no indication of any dependence of the mass ratio on the X-ray morphology of the clusters, indicating that the hydrostatic masses are not strongly systematically affected by the dynamical state of the clusters. Overall, our results favour a small value of the so-called hydrostatic bias due to non-thermal pressure sources.

preprint2016arXiv

The mass accretion rate of galaxy clusters: a measurable quantity

We explore the possibility of measuring the mass accretion rate (MAR) of galaxy clusters from their mass profiles beyond the virial radius $R_{200}$. We derive the accretion rate from the mass of a spherical shell whose inner radius is $2R_{200}$, whose thickness changes with redshift, and whose infall velocity is assumed to be equal to the mean infall velocity of the spherical shells of dark matter halos extracted from $N$-body simulations. This approximation is rather crude in hierarchical clustering scenarios where both smooth accretion and aggregation of smaller dark matter halos contribute to the mass accretion of clusters.Nevertheless, in the redshift range $z=[0,2]$, our prescription returns an average MAR within $20-40 \%$ of the average rate derived from the merger trees of dark matter halos extracted from $N$-body simulations. The MAR of galaxy clusters has been the topic of numerous detailed numerical and theoretical investigations, but so far it has remained inaccessible to measurements in the real universe. Since the measurement of the mass profile of clusters beyond their virial radius can be performed with the caustic technique applied to dense redshift surveys of the cluster outer regions, our result suggests that measuring the mean MAR of a sample of galaxy clusters is actually feasible. We thus provide a new potential observational test of the cosmological and structure formation models.

preprint2016arXiv

The potential role of NGC 205 in generating Andromeda's vast thin co-rotating plane of satellite galaxies

The Andromeda galaxy is observed to have a system of two large dwarf ellipticals and ~13 smaller satellite galaxies that are currently co-rotating in a thin plane, in addition to 2 counter-rotating satellite galaxies. We explored the consistency of those observations with a scenario where the majority of the co-rotating satellite galaxies originated from a subhalo group, where NGC 205 was the host and the satellite galaxies occupied dark matter sub-subhalos. We ran N-body simulations of a close encounter between NGC 205 and M31. In the simulations, NGC 205 was surrounded by massless particles to statistically sample the distribution of the sub-subhalos expected in a subhalo that has a mass similar to NGC 205. We made Monte Carlo samplings and found that, using a set of reference parameters, the probability of producing a thinner distribution of sub-subhalos than the observed NGC 205 + 15 smaller satellites (thus including the 2 counter-rotators, but excluding M32) increased from <1e-8 for the initial distribution to ~0.01 at pericentre. The probability of the simulated sub-subhalos occupying the locations of the observed co-rotating satellites in the line of sight velocity versus projected on-sky distance plane is at most 0.002 for 11 out of 13 satellites. Increasing the mass of M31 and the extent of the initial distribution of sub-subhalos gives a maximum probability of 0.004 for all 13 co-rotating satellites, but the probability of producing the thinness would drop to ~ 0.001.

preprint2016arXiv

The unrelaxed dynamical structure of the galaxy cluster Abell 85

For the first time, we explore the dynamics of the central region of a galaxy cluster within $r_{500}\sim 600h^{-1}$~kpc from its center by combining optical and X-ray spectroscopy. We use (1) the caustic technique that identifies the cluster substructures and their galaxy members with optical spectroscopic data, and (2) the X-ray redshift fitting procedure that estimates the redshift distribution of the intracluster medium (ICM). We use the spatial and redshift distributions of the galaxies and of the X-ray emitting gas to associate the optical substructures to the X-ray regions. When we apply this approach to Abell 85 (A85), a complex dynamical structure of A85 emerges from our analysis: a galaxy group, with redshift $z=0.0509 \pm 0.0021$ is passing through the cluster center along the line of sight dragging part of the ICM present in the cluster core; two additional groups, at redshift $z=0.0547 \pm 0.0022$ and $z=0.0570 \pm 0.0020$, are going through the cluster in opposite directions, almost perpendicularly to the line of sight, and have substantially perturbed the dynamics of the ICM. An additional group in the outskirts of A85, at redshift $z=0.0561 \pm 0.0023$, is associated to a secondary peak of the X-ray emission, at redshift $z=0.0583^{+0.0039}_{-0.0047}$. Although our analysis and results on A85 need to be confirmed by high-resolution spectroscopy, they demonstrate how our new approach can be a powerful tool to constrain the formation history of galaxy clusters by unveiling their central and surrounding structures.

preprint2015arXiv

GLACE survey: OSIRIS/GTC Tuneable Filter H$α$ imaging of the rich galaxy cluster ZwCl 0024.0+1652 at z = 0.395. Part I -- Survey presentation, TF data reduction techniques and catalogue

The cores of clusters at 0 $\lesssim$ z $\lesssim$ 1 are dominated by quiescent early-type galaxies, whereas the field is dominated by star-forming late-type ones. Galaxy properties, notably the star formation (SF) ability, are altered as they fall into overdense regions. The critical issues to understand this evolution are how the truncation of SF is connected to the morphological transformation and the responsible physical mechanism. The GaLAxy Cluster Evolution Survey (GLACE) is conducting a study on the variation of galaxy properties (SF, AGN, morphology) as a function of environment in a representative sample of clusters. A deep survey of emission line galaxies (ELG) is being performed, mapping a set of optical lines ([OII], [OIII], H$β$ and H$α$/[NII]) in several clusters at z $\sim$ 0.40, 0.63 and 0.86. Using the Tunable Filters (TF) of OSIRIS/GTC, GLACE applies the technique of TF tomography: for each line, a set of images at different wavelengths are taken through the TF, to cover a rest frame velocity range of several thousands km/s. The first GLACE results target the H$α$/[NII] lines in the cluster ZwCl 0024.0+1652 at z = 0.395 covering $\sim$ 2 $\times$ r$_{vir}$. We discuss the techniques devised to process the TF tomography observations to generate the catalogue of H$α$ emitters of 174 unique cluster sources down to a SFR below 1 M$_{\odot}$/yr. The AGN population is discriminated using different diagnostics and found to be $\sim$ 37% of the ELG population. The median SFR is 1.4 M$_{\odot}$/yr. We have studied the spatial distribution of ELG, confirming the existence of two components in the redshift space. Finally, we have exploited the outstanding spectral resolution of the TF to estimate the cluster mass from ELG dynamics, finding M$_{200}$ = 4.1 $\times$ 10$^{14}$ M$_{\odot} h^{-1}$, in agreement with previous weak-lensing estimates.

preprint2015arXiv

HeCS-SZ: The Hectospec Survey of Sunyaev-Zeldovich Selected Clusters

We estimate cluster masses and velocity dispersions for 123 clusters from optical spectroscopy to compare the Sunyaev-Zeldovich (SZ) mass proxy and dynamical masses. Our new survey, HeCS-SZ (Hectospec Cluster Survey of SZ-selected clusters), includes 7,721 new or remeasured redshifts from MMT/Hectospec observations of 24 SZ-selected clusters at redshifts $z$=0.05-0.20 and not in previous surveys. We supplement the Hectospec data with spectra from the Sloan Digital Sky Survey (SDSS) and cluster data from the Cluster Infall Regions in SDSS (CIRS) project and the Hectospec Cluster Survey (HeCS), our Hectospec survey of clusters selected by X-ray flux. We measure the scaling relation between velocity dispersion and SZ mass estimates from the integrated Compton parameter for an SZ complete sample of 83 clusters. The observed relation agrees very well with a simple virial scaling from mass (based on SZ) to velocity dispersion. The SZ mass estimates (calibrated with hydrostatic X-ray mass estimates) are not significantly biased. Further, the velocity dispersion of cluster galaxies is consistent with the expected velocity dispersion of dark matter particles, indicating that galaxies are good dynamical tracers (i.e., velocity bias is small). Significant mass bias in SZ mass estimates could relieve tension between cosmological results from Planck SZ cluster counts and Planck CMB data. However, the excellent agreement between our measured velocity dispersions and those predicted from a virial scaling relation suggests that any SZ mass bias is too small to reconcile SZ and CMB results. In principle, SZ mass bias and velocity bias of galaxies could conspire to yield good agreement, but the required velocity bias is $σ_{galaxy}\approx 0.77σ_{DM}$, outside the range of plausible models of velocity bias in the literature.

preprint2015arXiv

Identification of galaxy cluster substructures with the Caustic method

We investigate the power of the caustic technique for identifying substructures of galaxy clusters from optical redshift data alone. The caustic technique is designed to estimate the mass profile of galaxy clusters to radii well beyond the virial radius, where dynamical equilibrium does not hold. Two by-products of this technique are the identification of the cluster members and the identification of the cluster substructures. We test the caustic technique as a substructure detector on two samples of 150 mock redshift surveys of clusters; the clusters are extracted from a large cosmological $N$-body simulation of a $Λ$CDM model and have masses of $M_{200} \sim 10^{14} h^{-1} M_{\odot}$ and $M_{200} \sim 10^{15} h^{-1} M_{\odot}$ in the two samples. We limit our analysis to substructures identified in the simulation with masses larger than $10^{13} h^{-1} M_{\odot}$. With mock redshift surveys with 200 galaxies within $3R_{200}$, (1) the caustic technique recovers $\sim 30-50$\% of the real substructures, and (2) $\sim 15-20$\% of the substructures identified by the caustic technique correspond to real substructures of the central cluster, the remaining fraction being low-mass substructures, groups or substructures of clusters in the surrounding region, or chance alignments of unrelated galaxies. These encouraging results show that the caustic technique is a promising approach for investigating the complex dynamics of galaxy clusters.

preprint2015arXiv

Mass models of disk galaxies from the DiskMass Survey in MOND

This article explores the agreement between the predictions of Modified Newtonian Dynamics (MOND) and the rotation curves and stellar velocity dispersion profiles measured by the DiskMass Survey. A bulge-disk decomposition was made for each of the thirty published galaxies, and a MOND Poisson solver was used to simultaneously compute, from the baryonic mass distributions, model rotation curves and vertical velocity dispersion profiles, which were compared to the measured values. The two main free parameters, the stellar disk's mass-to-light ratio ($M/L$) and its exponential scale-height ($h_z$), were estimated by Markov Chain Monte Carlo modelling. The average best-fit K-band stellar mass-to-light ratio was $M/L \simeq 0.55 \pm 0.15$. However, to match the DiskMass Survey data, the vertical scale-heights would have to be in the range $h_z=200$ to $400$ pc which is a factor of two lower than those derived from observations of edge-on galaxies with a similar scale-length. The reason is that modified gravity versions of MOND characteristically require a larger $M/L$ to fit the rotation curve in the absence of dark matter and therefore predict a stronger vertical gravitational field than Newtonian models. It was found that changing the MOND acceleration parameter, the shape of the velocity dispersion ellipsoid, the adopted vertical distribution of stars, as well as the galaxy inclination, within any realistic range, all had little impact on these results.

preprint2015arXiv

Measuring the growth of galaxy clusters

We suggest how we can use the mass profile of galaxy clusters beyond their virial radius to measure their mass accretion rate, a key prediction of structure formation models. The mass profile can be estimated by applying the caustic technique to dense redshift surveys of clusters and their outskirts, where dynamical equilibrium does not necessarily hold. An additional probe of the mass growth of clusters is their mass fraction in substructures. We show that the caustic technique, that identifies cluster substructures as a by-product, returns catalogs of substructures with mass larger than a few $10^{13}h^{-1}M_\odot$ that are between 60% and 80% complete, depending on the density of the redshift survey.

preprint2014arXiv

A Redshift Survey of the Strong Lensing Cluster Abell 383

Abell 383 is a famous rich cluster (z = 0.1887) imaged extensively as a basis for intensive strong and weak lensing studies. Nonetheless there are few spectroscopic observations. We enable dynamical analyses by measuring 2360 new redshifts for galaxies with r$_{petro} \leq 20.5$ and within 50$^\prime$ of the BCG (Brightest Cluster Galaxy: R.A.$_{2000} = 42.014125^\circ$, Decl$_{2000} = -03.529228^\circ$). We apply the caustic technique to identify 275 cluster members within 7$h^{-1}$ Mpc of the hierarchical cluster center. The BCG lies within $-11 \pm 110$ km s$^{-1}$ and 21 $\pm 56 h^{-1}$ kpc of the hierarchical cluster center; the velocity dispersion profile of the BCG appears to be an extension of the velocity dispersion profile based on cluster members. The distribution of cluster members on the sky corresponds impressively with the weak lensing contours of Okabe et al. (2010) especially when the impact of foreground and background structure is included. The values of R$_{200}$ = $1.22\pm 0.01 h^{-1}$ Mpc and M$_{200}$ = $(5.07 \pm 0.09)\times 10^{14} h^{-1}$ M$_\odot$ obtained by application of the caustic technique agree well with recent completely independent lensing measures. The caustic estimate extends direct measurement of the cluster mass profile to a radius of $\sim 5 h^{-1}$ Mpc.

preprint2014arXiv

Comparing Dense Galaxy Cluster Redshift Surveys with Weak Lensing Maps

We use dense redshift surveys of nine galaxy clusters at $z\sim0.2$ to compare the galaxy distribution in each system with the projected matter distribution from weak lensing. By combining 2087 new MMT/Hectospec redshifts and the data in the literature, we construct spectroscopic samples within the region of weak-lensing maps of high (70--89%) and uniform completeness. With these dense redshift surveys, we construct galaxy number density maps using several galaxy subsamples. The shape of the main cluster concentration in the weak-lensing maps is similar to the global morphology of the number density maps based on cluster members alone, mainly dominated by red members. We cross correlate the galaxy number density maps with the weak-lensing maps. The cross correlation signal when we include foreground and background galaxies at 0.5$z_{\rm cl}<z<2z_{\rm cl}$ is $10-23$% larger than for cluster members alone at the cluster virial radius. The excess can be as high as 30% depending on the cluster. Cross correlating the galaxy number density and weak-lensing maps suggests that superimposed structures close to the cluster in redshift space contribute more significantly to the excess cross correlation signal than unrelated large-scale structure along the line of sight. Interestingly, the weak-lensing mass profiles are not well constrained for the clusters with the largest cross correlation signal excesses ($>$20% for A383, A689 and A750). The fractional excess in the cross correlation signal including foreground and background structures could be a useful proxy for assessing the reliability of weak-lensing cluster mass estimates.

preprint2014arXiv

Galaxy evolution in the mid-infrared green valley: a case of the A2199 supercluster

We study the mid-infrared (MIR) properties of the galaxies in the A2199 supercluster at z = 0.03 to understand the star formation activity of galaxy groups and clusters in the supercluster environment. Using the Wide-field Infrared Survey Explorer data, we find no dependence of mass-normalized integrated SFRs of galaxy groups/clusters on their virial masses. We classify the supercluster galaxies into three classes in the MIR color-luminosity diagram: MIR blue cloud (massive, quiescent and mostly early-type), MIR star-forming sequence (mostly late-type), and MIR green valley galaxies. These MIR green valley galaxies are distinguishable from the optical green valley galaxies, in the sense that they belong to the optical red sequence. We find that the fraction of each MIR class does not depend on virial mass of each group/cluster. We compare the cumulative distributions of surface galaxy number density and cluster/group-centric distance for the three MIR classes. MIR green valley galaxies show the distribution between MIR blue cloud and MIR SF sequence galaxies. However, if we fix galaxy morphology, early- and late-type MIR green valley galaxies show different distributions. Our results suggest a possible evolutionary scenario of these galaxies: 1) Late-type MIR SF sequence galaxies -> 2) Late-type MIR green valley galaxies -> 3) Early-type MIR green valley galaxies -> 4) Early-type MIR blue cloud galaxies. In this sequence, star formation of galaxies is quenched before the galaxies enter the MIR green valley, and then morphological transformation occurs in the MIR green valley.

preprint2014arXiv

Modified Baryonic Dynamics: two-component cosmological simulations with light sterile neutrinos

In this article we continue to test cosmological models centred on Modified Newtonian Dynamics (MOND) with light sterile neutrinos, which could in principle be a way to solve the fine-tuning problems of the standard model on galaxy scales while preserving successful predictions on larger scales. Due to previous failures of the simple MOND cosmological model, here we test a speculative model where the modified gravitational field is produced only by the baryons and the sterile neutrinos produce a purely Newtonian field (hence Modified Baryonic Dynamics). We use two component cosmological simulations to separate the baryonic N-body particles from the sterile neutrino ones. The premise is to attenuate the over-production of massive galaxy cluster halos which were prevalent in the original MOND plus light sterile neutrinos scenario. Theoretical issues with such a formulation notwithstanding, the Modified Baryonic Dynamics model fails to produce the correct amplitude for the galaxy cluster mass function for any reasonable value of the primordial power spectrum normalisation.

preprint2014arXiv

N-body simulations of the Carina dSph in MOND

The classical dwarf spheroidals (dSphs) provide a critical test for Modified Newtonian Dynamics (MOND) because they are observable satellite galactic systems with low internal accelerations and low, but periodically varying, external acceleration. This varying external gravitational field is not commonly found acting on systems with low internal acceleration. Using Jeans modelling, Carina in particular has been demonstrated to require a V-band mass-to-light ratio greater than 5, which is the nominal upper limit for an ancient stellar population. We run MOND N-body simulations of a Carina-like dSph orbiting the Milky Way to test if dSphs in MOND are stable to tidal forces over the Hubble time and if those same tidal forces artificially inflate their velocity dispersions and therefore their apparent mass-to-light ratio. We run many simulations with various initial total masses for Carina, and Galactocentric orbits (consistent with proper motions), and compare the simulation line of sight velocity dispersions (losVDs) with the observed losVDs of Walker et al. (2007). We find that the dSphs are stable, but that the tidal forces are not conducive to artificially inflating the losVDs. Furthermore, the range of mass-to-light ratios that best reproduces the observed line of sight velocity dispersions of Carina is 5.3 to 5.7 and circular orbits are preferred to plunging orbits. Therefore, some tension still exists between the required mass-to-light ratio for the Carina dSph in MOND and those expected from stellar population synthesis models. It remains to be seen whether a careful treatment of the binary population or triaxiality might reduce this tension.

preprint2014arXiv

The galaxy population of the complex cluster system Abell 3921

We present a spectrophotometric analysis of the galaxy pop. in the area of the merging cluster Abell 3921 at redshift 0.093. We investigate the impact of the complex cluster environment on galaxy properties such as morphology or star formation rate. We combine multi-object spectroscopy from the 2dF spectrograph with optical imaging taken with the ESO WFI. We carry out a redshift analysis and determine cluster velocity dispersions using biweight statistics. Applying a Dressler-Shectman (DS-)test we seek evidence for cluster substructure. Cluster and field galaxies are investigated with respect to [OII] and Hα equivalent width, SFR and morphological descriptors such as concentration index and Gini coefficient. We study these cluster galaxy properties as a function of clustercentric distance and investigate the spatial distribution of various galaxy types. Applying the DS-test we find a 3rd component (A3921-C) in addition to the two main subclusters (A3921-A and A3921-B) already known. The re-determined mass ratio between the main components A and B is approx. 2:1. Similar to prev. studies of galaxy clusters, we find that a large fraction of the disk galaxies close to the cluster core show no detectable star formation. These are likely systems that are quenched due to ram pressure stripping. We also find quenched spirals at rather large distances of 3 to 4 Mpc from the cluster core. A3921-C might be a group of galaxies falling onto the main cluster components. We speculate that the unexpected population of quenched spirals at large clustercentric radii in A3921-A and A3921-B might be an effect of the ongoing cluster merger: shocks in the ICM might give raise to enhanced ram pressure stripping and at least in part be the cause for the quenching of star formation. These quenched spirals might be an interm. stage in the morphological transformation of field spirals into cluster S0s.

preprint2013arXiv

Cosmological simulations in MOND: the cluster scale halo mass function with light sterile neutrinos

We use our Modified Newtonian Dynamics (MOND) cosmological particle-mesh N-body code to investigate the feasibility of structure formation in a framework involving MOND and light sterile neutrinos in the mass range 11 - 300 eV: always assuming that Ω_{ν_s}=0.225 for H_o=72 \kms Mpc^{-1}. We run a suite of simulations with variants on the expansion history, cosmological variation of the MOND acceleration constant, different normalisations of the power spectrum of the initial perturbations and interpolating functions. Using various box sizes, but typically with ones of length 256 Mpc/h, we compare our simulated halo mass functions with observed cluster mass functions and show that (i) the sterile neutrino mass must be larger than 30 eV to account for the low mass (M_{200}<10^{14.6} solar masses) clusters of galaxies in MOND and (ii) regardless of sterile neutrino mass or any of the variations we mentioned above, it is not possible to form the correct number of high mass (M_{200}>10^{15.1} solar masses) clusters of galaxies: there is always a considerable over production. This means that the ansatz of considering the weak-field limit of MOND together with a component of light sterile neutrinos to form structure from z ~ 200 fails. If MOND is the correct description of weak-field gravitational dynamics, it could mean that subtle effects of the additional fields in covariant theories of MOND render the ansatz inaccurate, or that the gravity generated by light sterile neutrinos (or by similar hot dark matter particles) is different from that generated by the baryons.

preprint2013arXiv

Identification of members in the central and outer regions of galaxy clusters

The caustic technique measures the mass of galaxy clusters in both their virial and infall regions and, as a byproduct, yields the list of cluster galaxy members. Here we use 100 galaxy clusters with mass M200>=1E14 Msun/h extracted from a cosmological N-body simulation of a LambdaCDM universe to test the ability of the caustic technique to identify the cluster galaxy members. We identify the true three-dimensional members as the gravitationally bound galaxies. The caustic technique uses the caustic location in the redshift diagram to separate the cluster members from the interlopers. We apply the technique to mock catalogues containing 1000 galaxies in the field of view of 12 Mpc/h on a side at the cluster location. On average, this sample size roughly corresponds to 180 real galaxy members within 3r200, similar to recent redshift surveys of cluster regions. The caustic technique yields a completeness, the fraction of identified true members, fc=0.95 (+- 0.03) within 3r200. The contamination increases from fi=0.020 (+0.046;-0.015) at r200 to fi=0.08 (+0.11;-0.05) at 3r200. No other technique for the identification of the members of a galaxy cluster provides such large completeness and small contamination at these large radii. The caustic technique assumes spherical symmetry and the asphericity of the cluster is responsible for most of the spread of the completeness and the contamination. By applying the technique to an approximately spherical system obtained by stacking the individual clusters, the spreads decrease by at least a factor of two. We finally estimate the cluster mass within 3r200 after removing the interlopers: for individual clusters, the mass estimated with the virial theorem is unbiased and within 30 per cent of the actual mass; this spread decreases to less than 10 per cent for the spherically symmetric stacked cluster.

preprint2013arXiv

Measuring the Ultimate Mass of Galaxy Clusters: Redshifts and Mass Profiles from the Hectospec Cluster Survey (HeCS)

The infall regions of galaxy clusters represent the largest gravitationally bound structures in a $Λ$CDM universe. Measuring cluster mass profiles into the infall regions provides an estimate of the ultimate mass of these haloes. We use the caustic technique to measure cluster mass profiles from galaxy redshifts obtained with the Hectospec Cluster Survey (HeCS), an extensive spectroscopic survey of galaxy clusters with MMT/Hectospec. We survey 58 clusters selected by X-ray flux at 0.1$<$$z$$<$0.3. The survey includes 21,314 unique MMT/Hectospec redshifts for individual galaxies; 10,275 of these galaxies are cluster members. For each cluster we acquired high signal-to-noise spectra for $\sim 200$ cluster members and a comparable number of foreground/background galaxies. The cluster members trace out infall patterns around the clusters. The members define a very narrow red sequence. The velocity dispersions decline with radius; we demonstrate that the determination of the velocity dispersion is insensitive to the inclusion of bluer members (a small fraction of the cluster population). We apply the caustic technique to define membership and estimate the mass profiles to large radii. The ultimate halo mass of clusters (the mass that remains bound in the far future of a $Λ$CDM universe) is on average (1.99$\pm$0.11)$M_{200}$, a new observational cosmological test in essential agreement with simulations. Summed profiles binned in $M_{200}$ and in $L_X$ demonstrate that the predicted NFW form of the density profile is a remarkably good representation of the data in agreement with weak lensing results extending to large radius. The concentration of these summed profiles is also consistent with theoretical predictions.

preprint2013arXiv

The dynamics of the bulge dominated galaxy NGC 7814 in MOND

The bulge dominated galaxy NGC 7814 provides one of the strongest dynamical tests possible for Modified Newtonian Dynamics (MOND). Spitzer 3.6 micron photometry fixes the bulge parameterisation and strongly constrains the properties of the sub-dominant stellar disk. Furthermore, the distance is known to better than 5 percent, virtually eliminating it as a free parameter. The rotation curve is easily measured, since the H I (and stellar) disks are edge on, and both the receding and approaching sides agree very well. In this paper we explore the agreement between the model and observed rotation curves in MOND given that the only two free parameters available are the mass-to-light ratios of the bulge and disk. We use a grid based MOND Poisson solver that accurately solves for the MOND gravity and produces our model rotation curves from a given mass distribution. The input to the Poisson solver is a 3D distribution of N particles which is generated from modelling the observed distribution of stars and gas in the galaxy. By ensuring a superior fit to the radial surface brightness profile than previous works, by virtue of a double Sersic fit to the bulge, we were able to produce excellent fits to the rotation curve with typical values for both mass-to-light ratios. We conclude that the model rotation curve of a mass distribution in MOND is extremely sensitive to the bulge-disk decomposition and even slight deviation from the observed mass distribution can produce large differences in the model rotation curve.

preprint2012arXiv

A WISE View of a Nearby Supercluster A2199

We use Wide-field Infrared Survey Explorer (WISE) data covering the entire region (~130 deg^2) of the A2199 supercluster at z=0.03 to study the mid-infrared (MIR) properties of supercluster galaxies. We identify a `MIR star-forming sequence' in the WISE [3.4]-[12] color-12 μm luminosity diagram, consisting of late-type, star-forming galaxies. At a fixed star formation rate (SFR), the MIR-detected galaxies at 22 μm or 12 μm tend to be more metal rich and to have higher surface brightness than those without MIR detection. Using these MIR-detected galaxies, we construct the IR luminosity function (LF) and investigate its environmental dependence. Both total IR (TIR) and 12 μm LFs are dominated by late-type, star-forming galaxies. The contribution of active galactic nuclei (AGN)-host galaxies increases with both TIR and 12 μm luminosities. The contribution of early-type galaxies to the 12 μm LFs increases with decreasing luminosity. The faint-end slope of the TIR LFs does not change with environment, but the change of faint-end slope in the 12 μm LFs with the environment is significant: there is a steeper faint-end slope in the cluster core than in the cluster outskirts. This steepening results primarily from the increasing contribution of early-type galaxies toward the cluster. These galaxies are passively evolving, and contain old stellar populations with weak MIR emission from the circumstellar dust around asymptotic giant branch stars.

preprint2012arXiv

An updated analysis of two classes of f(R) theories of gravity

The observed accelerated cosmic expansion can be a signature of fourth\,-\,order gravity theories, where the acceleration of the Universe is a consequence of departures from Einstein General Relativity, rather than the sign of the existence of a fluid with negative pressure. In the fourth\,-\,order gravity theories, the gravity Lagrangian is described by an analytic function $f(R)$ of the scalar curvature $R$ subject to the demanding conditions that no detectable deviations from standard GR is observed on the Solar System scale. Here we consider two classes of $f(R)$ theories able to pass Solar System tests and investigate their viability on cosmological scales. To this end, we fit the theories to a large dataset including the combined Hubble diagram of Type Ia Supernovae and Gamma Ray Bursts, the Hubble parameter $H(z)$ data from passively evolving red galaxies, Baryon Acoustic Oscillations extracted from the seventh data release of the Sloan Digital Sky Survey (SDSS) and the distance priors from the Wilkinson Microwave Anisotropy Probe seven years (WMAP7) data. We find that both classes of $f(R)$ fit very well this large dataset with the present\,-\,day values of the matter density, Hubble constant and deceleration parameter in agreement with previous estimates; however, the strong degeneracy among the $f(R)$ parameters prevents us from strongly constraining their values. We also derive the growth factor $g = d\lnδ/d\ln{a}$, with $δ= δρ_M/ρ_M$ the matter density perturbation, and show that it can still be well approximated by $g(z) \propto Ω_M(z)^γ$. We finally constrain $γ$ (on some representative scales) and investigate its redshift dependence to see whether future data can discriminate between these classes of $f(R)$ theories and standard dark energy models.

preprint2012arXiv

CLASH: Precise New Constraints on the Mass Profile of Abell 2261

We precisely constrain the inner mass profile of Abell 2261 (z=0.225) for the first time and determine this cluster is not "over-concentrated" as found previously, implying a formation time in agreement with ΛCDM expectations. These results are based on strong lensing analyses of new 16-band HST imaging obtained as part of the Cluster Lensing and Supernova survey with Hubble (CLASH). Combining this with revised weak lensing analyses of Subaru wide field imaging with 5-band Subaru + KPNO photometry, we place tight new constraints on the halo virial mass M_vir = 2.2\pm0.2\times10^15 M\odot/h70 (within r \approx 3 Mpc/h70) and concentration c = 6.2 \pm 0.3 when assuming a spherical halo. This agrees broadly with average c(M,z) predictions from recent ΛCDM simulations which span 5 <~ <c> <~ 8. Our most significant systematic uncertainty is halo elongation along the line of sight. To estimate this, we also derive a mass profile based on archival Chandra X-ray observations and find it to be ~35% lower than our lensing-derived profile at r2500 ~ 600 kpc. Agreement can be achieved by a halo elongated with a ~2:1 axis ratio along our line of sight. For this elongated halo model, we find M_vir = 1.7\pm0.2\times10^15 M\odot/h70 and c_vir = 4.6\pm0.2, placing rough lower limits on these values. The need for halo elongation can be partially obviated by non-thermal pressure support and, perhaps entirely, by systematic errors in the X-ray mass measurements. We estimate the effect of background structures based on MMT/Hectospec spectroscopic redshifts and find these tend to lower Mvir further by ~7% and increase cvir by ~5%.

preprint2012arXiv

Measuring the Mass Distribution in Galaxy Clusters

Cluster mass profiles are tests of models of structure formation. Only two current observational methods of determining the mass profile, gravitational lensing and the caustic technique, are independent of the assumption of dynamical equilibrium. Both techniques enable determination of the extended mass profile at radii beyond the virial radius. For 19 clusters, we compare the mass profile based on the caustic technique with weak lensing measurements taken from the literature. This comparison offers a test of systematic issues in both techniques. Around the virial radius, the two methods of mass estimation agree to within about 30%, consistent with the expected errors in the individual techniques. At small radii, the caustic technique overestimates the mass as expected from numerical simulations. The ratio between the lensing profile and the caustic mass profile at these radii suggests that the weak lensing profiles are a good representation of the true mass profile. At radii larger than the virial radius, the lensing mass profile exceeds the caustic mass profile possibly as a result of contamination of the lensing profile by large-scale structures within the lensing kernel. We highlight the case of the closely neighboring clusters MS0906+11 and A750 to illustrate the potential seriousness of contamination of the the weak lensing signal by unrelated structures.

preprint2012arXiv

The Acceleration Scale, Modified Newtonian Dynamics, and Sterile Neutrinos

General Relativity is able to describe the dynamics of galaxies and larger cosmic structures only if most of the matter in the Universe is dark, namely it does not emit any electromagnetic radiation. Intriguingly, on the scale of galaxies, there is strong observational evidence that the presence of dark matter appears to be necessary only when the gravitational field inferred from the distribution of the luminous matter falls below an acceleration of the order of 10^(-10) m/s^2. In the standard model, which combines Newtonian gravity with dark matter, the origin of this acceleration scale is challenging and remains unsolved. On the contrary, the full set of observations can be neatly described, and were partly predicted, by a modification of Newtonian dynamics, dubbed MOND, that does not resort to the existence of dark matter. On the scale of galaxy clusters and beyond, however, MOND is not as successful as on the scale of galaxies, and the existence of some dark matter appears unavoidable. A model combining MOND with hot dark matter made of sterile neutrinos seems to be able to describe most of the astrophysical phenomenology, from the power spectrum of the cosmic microwave background anisotropies to the dynamics of dwarf galaxies. Whether there exists a yet unknown covariant theory that contains General Relativity and Newtonian gravity in the weak field limit, and MOND as the ultra-weak field limit is still an open question.

preprint2012arXiv

The Faint End of the Luminosity Function and Low Surface Brightness Galaxies

SHELS (Smithsonian Hectospec Lensing Survey) is a dense redshift survey covering a 4 square degree region to a limiting R = 20.6. In the construction of the galaxy catalog and in the acquisition of spectroscopic targets, we paid careful attention to the survey completeness for lower surface brightness dwarf galaxies. Thus, although the survey covers a small area, it is a robust basis for computation of the slope of the faint end of the galaxy luminosity function to a limiting M_R = -13.3 + 5logh. We calculate the faint end slope in the R-band for the subset of SHELS galaxies with redshif ts in the range 0.02 <= z < 0.1, SHELS_{0.1}. This sample contains 532 galaxies with R< 20.6 and with a median surface brightness within the half light radius of SB_{50,R} = 21.82 mag arcsec^{-2}. We used this sample to make one of the few direct measurements of the dependence of the faint end of the galaxy luminosity function on surface brightness. For the sample as a whole the faint end slope, alpha = -1.31 +/- 0.04, is consistent with both the Blanton et al. (2005b) analysis of the SDSS and the Liu et al. (2008) analysis of the COSMOS field. This consistency is impressive given the very different approaches of th ese three surveys. A magnitude limited sample of 135 galaxies with optical spectroscopic reds hifts with mean half-light surface brightness, SB_{50,R} >= 22.5 mag arcsec^{-2} is unique to SHELS_{0.1}. The faint end slope is alpha_{22.5} = -1.52+/- 0.16. SHELS_{0.1} shows that lower surface brightness objects dominate the faint end slope of the l uminosity function in the field, underscoring the importance of surface brightness limits in evaluating measurements of the faint end slope and its evolution.

preprint2011arXiv

Gamma-ray bursts as cosmological probes: LambdaCDM vs. conformal gravity

LambdaCDM, for the currently preferred cosmological density Omega_0 and cosmological constant Omega_Lambda, predicts that the Universe expansion decelerates from early times to redshift z~0.9 and accelerates at later times. On the contrary, the cosmological model based on conformal gravity predicts that the cosmic expansion has always been accelerating. To distinguish between these two very different cosmologies, we resort to gamma-ray bursts (GRBs), which have been suggested to probe the Universe expansion history at z>1, where identified type Ia supernovae (SNe) are rare. We use the full Bayesian approach to infer the cosmological parameters and the additional parameters required to describe the GRB data available in the literature. For the first time, we use GRBs as cosmological probes without any prior information from other data. In addition, when we combine the GRB samples with SNe, our approach neatly avoids all the inconsistencies of most numerous previous methods that are plagued by the so-called circularity problem. In fact, when analyzed properly, current data are consistent with distance moduli of GRBs and SNe that can respectively be, in a variant of conformal gravity, ~15 and ~3 magnitudes fainter than in LambdaCDM. Our results indicate that the currently available SN and GRB samples are accommodated equally well by both LambdaCDM and conformal gravity and do not exclude a continuous accelerated expansion. We conclude that GRBs are currently far from being effective cosmological probes, as they are unable to distinguish between these two very different expansion histories.

preprint2011arXiv

How unusual are the Shapley Supercluster and the Sloan Great Wall?

We use extreme value statistics to assess the significance of two of the most dramatic structures in the local Universe: the Shapley supercluster and the Sloan Great Wall. If we assume that Shapley (volume ~ 1.2 x 10^5 (Mpc/h)^3) evolved from an overdense region in the initial Gaussian fluctuation field, with currently popular choices for the background cosmological model and the shape and amplitude sigma8 of the initial power spectrum, we estimate that the total mass of the system is within 20 percent of 1.8 x 10^16 Msun/h. Extreme value statistics show that the existence of this massive concentration is not unexpected if the initial fluctuation field was Gaussian, provided there are no other similar objects within a sphere of radius 200 Mpc/h centred on our Galaxy. However, a similar analysis of the Sloan Great Wall, a more distant (z ~ 0.08) and extended concentration of structures (volume ~ 7.2 x 10^5 (Mpc/h)^3) suggests that it is more unusual. We estimate its total mass to be within 20 percent of 1.2 x 10^17 Msun/h; even if it is the densest such object of its volume within z=0.2, its existence is difficult to reconcile with Gaussian initial conditions if sigma8 < 0.9. This tension can be alleviated if this structure is the densest within the Hubble volume. Finally, we show how extreme value statistics can be used to address the likelihood that an object like Shapley exists in the same volume which contains the Great Wall, finding, again, that Shapley is not particularly unusual. It is straightforward to incorporate other models of the initial fluctuation field into our formalism.

preprint2011arXiv

Mapping the Universe: The 2010 Russell Lecture

Redshift surveys are a powerful tool of modern cosmology. We discuss two aspects of their power to map the distribution of mass and light in the universe: (1) measuring the mass distribution extending into the infall regions of rich clusters and (2) applying deep redshift surveys to the selection of clusters of galaxies and to the identification of very large structures (Great Walls). We preview the HectoMAP project, a redshift survey with median redshift z = 0.34 covering 50 square degrees to r= 21. We emphasize the importance and power of spectroscopy for exploring and understanding the nature and evolution of structure in the universe.

preprint2011arXiv

Measuring Unified Dark Matter with 3D Cosmic Shear

We present parameter estimation forecasts for future 3D cosmic shear surveys for a class of Unified Dark Matter (UDM) models, where a single scalar field mimics both Dark Matter (DM) and Dark Energy (DE). These models have the advantage that they can describe the dynamics of the Universe with a single matter component providing an explanation for structure formation and cosmic acceleration. A crucial feature of the class of UDM models we use in this work is characterised by a parameter, c_inf (in units of the speed of light c=1), that is the value of the sound speed at late times, and on which structure formation depends. We demonstrate that the properties of the DM-like behaviour of the scalar field can be estimated with very high precision with large-scale, fully 3D weak lensing surveys. We found that 3D weak lensing significantly constrains c_inf, and we find minimal errors 0.00003, for the fiducial value c_inf=0.001, and 0.000026, for c_inf=0.012. Moreover, we compute the Bayesian evidence for UDM models over the LCDM model as a function of c_inf. For this purpose, we can consider the LCDM model as a UDM model with c_inf=0. We find that the expected evidence clearly shows that the survey data would unquestionably favour UDM models over the LCDM model, for the values c_inf>0.001.

preprint2011arXiv

Testing a Phenomenologically Extended DGP Model with Upcoming Weak Lensing Surveys

A phenomenological extension of the well-known brane-world cosmology of Dvali, Gabadadze and Porrati (eDGP) has recently been proposed. In this model, a cosmological-constant-like term is explicitly present as a non-vanishing tension sigma on the brane, and an extra parameter alpha tunes the cross-over scale r_c, the scale at which higher dimensional gravity effects become non negligible. Since the Hubble parameter in this cosmology reproduces the same LCDM expansion history, we study how upcoming weak lensing surveys, such as Euclid and DES (Dark Energy Survey), can confirm or rule out this class of models. We perform Markov Chain Monte Carlo simulations to determine the parameters of the model, using Type Ia Supernovæ, H(z) data, Gamma Ray Bursts and Baryon Acoustic Oscillations. We also fit the power spectrum of the temperature anisotropies of the Cosmic Microwave Background to obtain the correct normalisation for the density perturbation power spectrum. Then, we compute the matter and the cosmic shear power spectra, both in the linear and non-linear regimes. The latter is calculated with the two different approaches of Hu and Sawicki (2007) (HS) and Khoury and Wyman (2009) (KW). With the eDGP parameters coming from the Markov Chains, KW reproduces the LCDM matter power spectrum at both linear and non-linear scales and the LCDM and eDGP shear signals are degenerate. This result does not hold with the HS prescription: Euclid can distinguish the eDGP model from LCDM because their expected power spectra roughly differ by the 3sigma uncertainty in the angular scale range 700<l<3000; on the contrary, the two models differ at most by the 1sigma uncertainty over the range 500<l<3000 in the DES experiment and they are virtually indistinguishable.

preprint2011arXiv

The abundance of galaxy clusters in MOND: Cosmological simulations with massive neutrinos

We present a new Particle-Mesh cosmological N-body code for accurately solving the modified Poisson equation of the Quasi Linear formulation of MOND. We generate initial conditions for the Angus (2009) cosmological model, which is identical to LCDM except that the cold dark matter is switched for a single species of thermal sterile neutrinos. We set the initial conditions at z=250 for a (512 Mpc/h)^3 box with 256^3 particles and we evolve them down to z=0. We clearly demonstrate the necessity of MOND for developing the large scale structure in a hot dark matter cosmology and contradict the naive expectation that MOND cannot form galaxy clusters. We find that the correct order of magnitude of X-ray clusters (with T_X > 4.5 keV) can be formed, but that we overpredict the number of very rich clusters and seriously underpredict the number of lower mass clusters. The latter is a shortcoming of the resolution of our simulations, whereas we suggest that the over production of very rich clusters might be prevented by incorporating a MOND acceleration constant that varies with redshift and an expansion history that cannot be described by the usual Friedmann models. We present evidence that suggests the density profiles of our simulated clusters are compatible with those of observed X-ray clusters in MOND. It remains to be seen if the low mass end of the cluster mass function can be reproduced and if the high densities of dark matter in the central 20 kpc of groups and clusters of galaxies, measured in the MOND framework, can be achieved. As a last test, we computed the relative velocity between pairs of halos within 10 Mpc and find that pairs with velocities larger than 3000 km/s like the bullet cluster, can form without difficulty.

preprint2011arXiv

The Power of Cosmic Flexion in Testing Modified Matter and Gravity

Flexion is the weak lensing effect responsible for the weakly skewed and arc-like appearance of lensed galaxies. The flexion signal-to-noise ratio can be an order of magnitude larger than that of shear. For the first time, we show how this makes flexion an invaluable tool for discriminating among alternative cosmological models. We analyse a scalar field model of unified dark matter and dark energy, a brane-world cosmology and two f(R) modified-action theories. We show that these models can be distinguished from LCDM at several standard deviations by measuring the power spectrum of cosmic flexion.

preprint2011arXiv

Tomography from the Next Generation of Cosmic Shear Experiments for Viable f(R) Models

We present the cosmic shear signal predicted by two viable cosmological models in the framework of modified-action f(R) theories. We use f(R) models where the current accelerated expansion of the Universe is a direct consequence of the modified gravitational Lagrangian rather than Dark Energy (DE), either in the form of vacuum energy/cosmological constant or of a dynamical scalar field (e.g. quintessence). We choose Starobinsky's (St) and Hu & Sawicki's (HS) f(R) models, which are carefully designed to pass the Solar System gravity tests. In order to further support - or rule out - f(R) theories as alternative candidates to the DE hypothesis, we exploit the power of weak gravitational lensing, specifically of cosmic shear. We calculate the tomographic shear matrix as it would be measured by the upcoming ESA Cosmic Vision Euclid satellite. We find that in the St model the cosmic shear signal is almost completely degenerate with LCDM, but it is easily distinguishable in the HS model. Moreover, we compute the corresponding Fisher matrix for both the St and HS models, thus obtaining forecasts for their cosmological parameters. Finally, we show that the Bayes factor for cosmic shear will definitely favour the HS model over LCDM if Euclid measures a value larger than ~0.02 for the extra HS parameter n_HS.

preprint2011arXiv

Using dwarf satellite proper motions to determine their origin

The highly organised distribution of satellite galaxies surrounding the Milky Way is a serious challenge to the concordance cosmological model. Perhaps the only remaining solution, in this framework, is that the dwarf satellite galaxies fall into the Milky Way's potential along one or two filaments, which may or may not plausibly reproduce the observed distribution. Here we test this scenario by making use of the proper motions of the Fornax, Sculptor, Ursa Minor and Carina dwarf spheroidals, and trace their orbits back through several variations of the Milky Way's potential and account for dynamical friction. The key parameters are the proper motions and total masses of the dwarf galaxies. Using a simple model we find no tenable set of parameters that can allow Fornax to be consistent with filamentary infall, mainly because the 1 sigma error on its proper motion is relatively small. The other three must walk a tightrope between requiring a small pericentre (less than 20 kpc) to lose enough orbital energy to dynamical friction and avoiding being tidally disrupted. We then employed a more realistic model with host halo mass accretion, and found that the four dwarf galaxies must have fallen in at least 5 Gyrs ago. This time interval is longer than organised distribution is expected to last before being erased by the randomisation of the satellite orbits.

preprint2010arXiv

A sub-resolution multiphase interstellar medium model of star formation and SNe energy feedback

We present a new multi-phase sub-resolution model for star formation and feedback in SPH numerical simulations of galaxy formation. Our model, called MUPPI (MUlti-Phase Particle Integrator), describes each gas particle as a multi-phase system, with cold and hot gas phases, coexisting in pressure equilibrium, and a stellar component. Cooling of the hot tenuous gas phase feeds the cold gas phase. Stars are formed out of molecular gas with a given efficiency, which scales with the dynamical time of the cold phase. Our prescription for star formation is not based on imposing the Schmidt-Kennicutt relation, which is instead naturally produced by MUPPI. Energy from supernova explosions is deposited partly into the hot phase of the gas particles, and partly to that of neighboring particles. Mass and energy flows among the different phases of each particle are described by a set of ordinary differential equations which we explicitly integrate for each gas particle, instead of relying on equilibrium solutions. This system of equations also includes the response of the multi-phase structure to energy changes associated to the thermodynamics of the gas. We apply our model to two isolated disk galaxy simulations and two spherical cooling flows. MUPPI is able to reproduce the Schmidt-Kennicutt relation for disc galaxies. It also reproduces the basic properties of the inter-stellar medium in disc galaxies, the surface densities of cold and molecular gas, of stars and of star formation rate, the vertical velocity dispersion of cold clouds and the flows connected to the galactic fountains. Quite remarkably, MUPPI also provides efficient stellar feedback without the need to include a scheme of kinetic energy feedback. [abridged]

preprint2010arXiv

Comparison of Hectospec Virial Masses with SZE Measurements

We present the first comparison of virial masses of galaxy clusters with their Sunyaev-Zel'dovich Effect (SZE) signals. We study 15 clusters from the Hectospec Cluster Survey (HeCS) with MMT/Hectospec spectroscopy and published SZE signals. We measure virial masses of these clusters from an average of 90 member redshifts inside the radius $r_{100}$. The virial masses of the clusters are strongly correlated with their SZE signals (at the 99% confidence level using a Spearman rank-sum test). This correlation suggests that $Y_{SZ}$ can be used as a measure of virial mass. Simulations predict a powerlaw scaling of $Y_{SZ}\propto M_{200}^α$ with $α\approx$1.6. Observationally, we find $α$=1.11$\pm$0.16, significantly shallower (given the formal uncertainty) than the theoretical prediction. However, the selection function of our sample is unknown and a bias against less massive clusters cannot be excluded (such a selection bias could artificially flatten the slope). Moreover, our sample indicates that the relation between velocity dispersion (or virial mass estimate) and SZE signal has significant intrinsic scatter, comparable to the range of our current sample. More detailed studies of scaling relations are therefore needed to derive a robust determination of the relation between cluster mass and SZE.

preprint2010arXiv

Implications for dwarf spheroidal mass content from interloper removal

Using the caustic method, we identify the member stars of five dwarf spheroidal (dSph) galaxies of the Milky Way, the smallest dark matter (DM) dominated systems in the Universe. After our interloper rejection, we compute line-of-sight velocity dispersion profiles that are substantially smoother than previous results. Moreover, two dSphs have line-of-sight velocity dispersions 20% smaller than previous calculations suggested. Our Jeans modelling confirms that the DM content interior to 300pc is roughly constant with satellite luminosity. Finally, if we assume that MOND provides the true law of gravity, our identification of interlopers implies that four dSphs have mass-to-light ratios in agreement with stellar population synthesis models, whereas Carina still has a mass-to-light ratio a factor of two too large and remains a problem for MOND.

preprint2010arXiv

Infall Regions and Scaling Relations of X-ray Selected Groups

We use the Fifth Data Release of the Sloan Digital Sky Survey to study X-ray-selected galaxy groups and compare their properties to clusters. We search for infall patterns around the groups and use these to measure group mass profiles to large radii. In previous work, we analyzed infall patterns for an X-ray-selected sample of 72 clusters from the ROSAT All-Sky Survey. Here, we extend this approach to a sample of systems with smaller X-ray fluxes selected from the 400 Square Degree serendipitous survey of clusters and groups in ROSAT pointed observations. We identify 16 groups with SDSS DR5 spectroscopy, search for infall patterns, and compute mass profiles out to 2-6 Mpc from the group centers with the caustic technique. No other mass estimation methods are currently available at such large radii for these low-mass groups, because the virial estimate requires dynamical equilibrium and the gravitational lensing signal is too weak. Despite the small masses of these groups, most display recognizable infall patterns. We use caustic and virial mass estimates to measure the scaling relations between different observables, extending these relations to smaller fluxes and luminosities than many previous surveys. Close inspection reveals that three of the groups are subclusters in the outskirts of larger clusters. A fourth group is apparently undergoing a group-group merger. These four merging groups represent the most extreme outliers in the scaling relations. Excluding these groups, we find $L_X\proptoσ_p^{3.4\pm1.6}$, consistent with previous determinations for both clusters and groups. Understanding cluster and group scaling relations is crucial for measuring cosmological parameters from clusters.

preprint2010arXiv

Measuring the escape velocity and mass profiles of galaxy clusters beyond their virial radius

The caustic technique uses galaxy redshifts alone to measure the escape velocity and mass profiles of galaxy clusters to clustrocentric distances well beyond the virial radius, where dynamical equilibrium does not necessarily hold. We provide a detailed description of this technique and analyse its possible systematic errors. We apply the caustic technique to clusters with mass M_200>=10^{14}h^{-1} M_sun extracted from a cosmological hydrodynamic simulation of a LambdaCDM universe. With a few tens of redshifts per squared comoving megaparsec within the cluster, the caustic technique, on average, recovers the profile of the escape velocity from the cluster with better than 10 percent accuracy up to r~4 r_200. The caustic technique also recovers the mass profile with better than 10 percent accuracy in the range (0.6-4) r_200, but it overestimates the mass up to 70 percent at smaller radii. This overestimate is a consequence of neglecting the radial dependence of the filling function F_beta(r). The 1-sigma uncertainty on individual escape velocity profiles increases from ~20 to ~50 percent when the radius increases from r~0.1 r_200 to ~4 r_200. Individual mass profiles have 1-sigma uncertainty between 40 and 80 percent within the radial range (0.6-4) r_200. We show that the amplitude of these uncertainties is completely due to the assumption of spherical symmetry, which is difficult to drop. Alternatively, we can apply the technique to synthetic clusters obtained by stacking individual clusters: in this case, the 1-sigma uncertainty on the escape velocity profile is smaller than 20 percent out to 4 r_200. The caustic technique thus provides reliable average profiles which extend to regions difficult or impossible to probe with other techniques.

preprint2009arXiv

Equilibrium configurations of 11eV sterile neutrinos in MONDian galaxy clusters

(Abridged) Modified Newtonian Dynamics (MOND) can fit a broad range of galaxy kinematic data, but struggles with clusters of galaxies. MONDian clusters need dark matter, and here we test the 11 eV sterile neutrino - used to fit the first three acoustic peaks of the cosmic microwave background - by investigating their equilibrium distributions in 30 groups and clusters over a wide range of temperatures. Knowing both the sterile neutrino density and velocity dispersion tells us the Tremaine-Gunn phase-space limit at all radii. We find that all 30 systems serendipitously reach the Tremaine-Gunn limit by the centre, which means a portion of the dynamical mass must always be covered by the brightest cluster galaxy. Interestingly, the typical fitted K-band mass-to-light ratio is unity and at most 1.2, which is very consistent - although leaving no margin for error - with stellar population synthesis models. Amidst the sample there are several special cases including the Coma cluster (for which dark matter was first proposed), NGC 720 (where geometrical evidence for dark matter was found) and the bullet cluster (where dark matter - of some kind - in clusters was directly proven to exist). We demonstrate that 11 eV sterile neutrinos are unlikely to influence spiral galaxy rotation curves, as they don't influence even some very massive early-types (NGC 4125 and NGC 6482). Finally, we conclude that it is intriguing that the minimum mass of sterile neutrino particle that can match the cosmic microwave background is identical to the minimum mass found here to be consistent with equilibrium configurations of MONDian clusters of galaxies.

preprint2009arXiv

Resolving the timing problem of the globular clusters orbiting the Fornax dwarf galaxy

We re-investigate the old problem of the survival of the five globular clusters orbiting the Fornax dwarf galaxy in both standard and modified Newtonian dynamics. For the first time in the history of the topic, we use accurate mass models for the Fornax dwarf, obtained through Jeans modelling of the recently published line of sight velocity dispersion data, and we are also not resigned to circular orbits for the globular clusters. Previously conceived problems stem from fixing the starting distances of the globulars to be less than half the tidal radius. We relax this constraint since there is absolutely no evidence for it and show that the dark matter paradigm, with either cusped or cored dark matter profiles, has no trouble sustaining the orbits of the two least massive globular clusters for a Hubble time almost regardless of their initial distance from Fornax. The three most massive globulars can remain in orbit as long as their starting distances are marginally outside the tidal radius. The outlook for modified Newtonian dynamics is also not nearly as bleak as previously reported. Although dynamical friction inside the tidal radius is far stronger in MOND, outside dynamical friction is negligible due to the absence of stars. This allows highly radial orbits to survive, but more importantly circular orbits at distances more than 85% of Fornax's tidal radius to survive indefinitely. The probability of the globular clusters being on circular orbits at this distance compared with their current projected distances is discussed and shown to be plausible. Finally, if we ignore the presence of the most massive globular (giving it a large line of sight distance) we demonstrate that the remaining four globulars can survive within the tidal radius for the Hubble time with perfectly sensible orbits.

preprint2009arXiv

Velocity Dispersion Profile of the Milky Way Halo

We present a spectroscopic sample of 910 distant halo stars from the Hypervelocity Star survey from which we derive the velocity dispersion profile of the Milky Way halo. The sample is a mix of 74% evolved horizontal branch stars and 26% blue stragglers. We estimate distances to the stars using observed colors, metallicities, and stellar evolution tracks. Our sample contains twice as many objects with R>50 kpc as previous surveys. We compute the velocity dispersion profile in two ways: with a parametric method based on a Milky Way potential model, and with a non-parametric method based on the caustic technique originally developed to measure galaxy cluster mass profiles. The resulting velocity dispersion profiles are remarkably consistent with those found by two independent surveys based on other stellar populations: the Milky Way halo exhibits a mean decline in radial velocity dispersion of -0.38+-0.12 km/s/kpc over 15<R<75 kpc. This measurement is a useful basis for calculating the total mass and mass distribution of the Milky Way halo.

preprint2009arXiv

Weak lensing signal in Unified Dark Matter models

Weak gravitational lensing is a powerful tool for studying both the geometry and the dynamics of the Universe. Its power spectrum contains information on the sources emitting photons and on the large--scale structures that these hotons cross. We calculate the weak lensing cosmic convergence and shear power spectra, in linear theory and Limber's approximation, for two different classes of cosmological models: the standard \LambdaCDM and Unified Dark Matter (UDM) models. The latter models attempt to solve the problems of the dark matter in the dynamics of galaxies and galaxy clusters and of the late-time acceleration of the Universe expansion by introducing a scalar field that mimics both dark matter and dark energy. A crucial feature of the UDM models is the speed of sound c_\infty, that is the value of the sound speed at late times, on which structure formation depends. In this paper, we provide the predictions of the UDM models for the weak lensing signal,with various values of c_\infty. We consider both the Cosmic Microwave Background and background galaxies at different redshifts as sources for the weak lensing power spectra. We find that the power spectra in UDM models are more sensitive to the variations of c_\infty for sources located at low redshifts. Moreover, we find that for c_\infty>10^{-3} (in units of the speed of light), the UDM weak lensing convergence power spectrum C^{κκ}(l) for background galaxies is strongly suppressed with respect to the \LambdaCDM spectrum, particularly at small angular scales l\gtrsim100.

preprint2008arXiv

How the Scalar Field of Unified Dark Matter Models Can Cluster

We use scalar-field Lagrangians with a non-canonical kinetic term to obtain unified dark matter models where both the dark matter and the dark energy, the latter mimicking a cosmological constant, are described by the scalar field itself. In this framework, we propose a technique to reconstruct models where the effective speed of sound is small enough that the scalar field can cluster. These models avoid the strong time evolution of the gravitational potential and the large Integrated Sachs-Wolfe effect which have been a serious drawback of previously considered models. Moreover, these unified dark matter scalar field models can be easily generalized to behave as dark matter plus a dark energy component behaving like any type of quintessence fluid.

preprint2005arXiv

CAIRNS: The Cluster And Infall Region Nearby Survey III. Environmental Dependence of H-alpha Properties of Galaxies

We investigate the environmental dependence of star formation in cluster virial regions and infall regions as part of CAIRNS (Cluster And Infall Region Nearby Survey), a large spectroscopic survey of the infall regions surrounding nine nearby rich clusters of galaxies. Our long-slit spectroscopy yields estimates of star formation rates in environments from cluster cores to the general large-scale structure. The fraction of galaxies with current star formation in their inner disks as traced by H-alpha emission increases with distance from the cluster and converges to the ``field'' value only at 2-3 virial radii, in agreement with other investigations. However, among galaxies with significant current star formation (EW[Ha]geq2Å), there is no difference in the distribution of EW[Ha] inside and outside the virial radius. This surprising result, first seen by Carter et al., suggests that (1) star formation is truncated on either very short timescales or only at moderate and high redshifts or (2) that projection effects contaminate the measurement. The number density profiles of star-forming and non-star-forming galaxies indicate that, among galaxies projected inside the virial radius, at least half of the former and 20% of the latter are ``infall interlopers,'' galaxies in the infall region but outside the virial region. The kinematics of star-forming galaxies in the infall region closely match those of absorption-dominated galaxies. This result shows that the star forming galaxies in the infall regions are not interlopers from the field and excludes one model of the backsplash scenario of galaxy transformation. Finally, we quantify systematic uncertainties in estimating the global star formation in galaxies from their inner disks.