Source author record

Gary A. Mamon

Gary A. Mamon 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

38works
10topics
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

38 published item(s)

preprint2022arXiv

Stellar graveyards: Clustering of compact objects in globular clusters NGC 3201 and NGC 6397

We analyse Gaia EDR3 and re-calibrated HST proper motion data from the core-collapsed and non core-collapsed globular clusters NGC 6397 and NGC 3201, respectively, with the Bayesian mass-orbit modelling code MAMPOSSt-PM. We use Bayesian evidence and realistic mock data sets constructed with AGAMA to select between different mass models. In both clusters, the velocities are consistent with isotropy within the extent of our data. We robustly detect a dark central mass (DCM) of roughly 1000 solar masses in both clusters. Our MAMPOSSt-PM fits strongly prefer an extended DCM in NGC 6397, while only presenting a mild preference for it in NGC 3201, with respective sizes of a roughly one and a few per cent of the cluster effective radius. We explore the astrophysics behind our results with the CMC Monte Carlo N-body code, whose snapshots best matching the phase space observations lead to similar values for the mass and size of the DCM. The internal kinematics are thus consistent with a population of hundreds of massive white dwarfs in NGC 6397, and roughly 100 segregated stellar-mass black holes in NGC 3201, as previously found with CMC. Such analyses confirm the accuracy of both mass-orbit modelling and Monte Carlo N-body techniques, which together provide more robust predictions on the DCM of globular clusters (core-collapsed or not). This opens possibilities to understand a vast range of interesting astrophysical phenomena in clusters, such as fast radio bursts, compact object mergers, and gravitational waves.

preprint2022arXiv

Theia : science cases and mission profiles for high precision astrometry in the future

High-precision astrometry well beyond the capacities of Gaia will provide a unique way to achieve astrophysical breakthroughs, in particular on the nature of dark matter, and a complete survey of nearby habitable exoplanets. In this contribution, we present the scientific cases that require a flexibly-pointing instrument capable of high astrometric accuracy and we review the best mission profiles that can achieve such observations with the current space technology as well as within the boundary conditions defined by space agencies. We also describe the way the differential astrometric measurement is made using reference stars within the field. We show that the ultimate accuracy can be met without drastic constrains on the telescope stability.

preprint2021arXiv

Does NGC 6397 contain an intermediate-mass black hole or a more diffuse inner subcluster?

We analyze proper motions from the Hubble Space Telescope (HST) and the second Gaia data release along with line-of-sight velocities from the MUSE spectrograph to detect imprints of an intermediate-mass black hole (IMBH) in the center of the nearby, core-collapsed, globular cluster NGC 6397. For this, we use the new MAMPOSSt-PM Bayesian mass-modeling code, along with updated estimates of the surface density profile of NGC 6397. We consider different priors on velocity anisotropy and on the size of the central mass, and we also separate the stars into components of different mean mass to allow for mass segregation. The velocity ellipsoid is very isotropic throughout the cluster, as expected in post-core collapsed clusters subject to as strong a Galactic tidal field as NGC 6397. There is strong evidence for a central dark component of 0.8 to 2% of the total mass of the cluster. However, we find robust evidence disfavoring a central IMBH in NGC 6397, preferring instead a diffuse dark inner subcluster of unresolved objects with a total mass of 1000 to 2000 solar masses, half of which is concentrated within 6 arcsec (2% of the stellar effective radius). These results require the combination of HST and Gaia data: HST for the inner diagnostics and Gaia for the outer surface density and velocity anisotropy profiles. The small effective radius of the diffuse dark component suggests that it is composed of compact stars (white dwarfs and neutron stars) and stellar-mass black holes, whose inner locations are caused by dynamical friction given their high progenitor masses. We show that stellar-mass black holes should dominate the mass of this diffuse dark component, unless more than 25 per cent escape from the cluster. Their mergers in the cores of core-collapsed globular clusters could be an important source of the gravitational wave events detected by LIGO.

preprint2020arXiv

A Precise Analytical Approximation for the Deprojection of the Sérsic Profile

The Sérsic model is known to fit well the surface brightness (or surface density) profiles of elliptical galaxies and galaxy bulges, and possibly for dwarf spheroidal galaxies and globular clusters. The deprojected density and mass profiles are important for many astrophysical applications, in particular for mass-orbit modeling of these systems. However, the exact deprojection formula for the Sérsic model employs special functions not available in most computer languages. We show that all previous analytical approximations to the 3D density profile are imprecise at low Sérsic index ($n \lesssim 1.5$). We have derived a more precise analytical approximation to the deprojected Sérsic density profile by fitting two-dimensional 10th-order polynomials to the differences of the logarithms of the numerical deprojection and of the analytical approximation by Lima Neto et al. (1999, LGM) of the density profile on one hand and of the mass profile on the other. Our LGM-based polynomial fits have typical relative precision better than $0.2\%$ for both density and mass profiles, for Sérsic indices $0.5 \leq n \leq 10$ and radii $0.001 < r/R_{\rm e} < 1000$. Our approximation is much more precise than those of LGM, Simonneau & Prada (1999, 2004), Trujillo et al. (2002) for non-half-integer values of the index, and of Emsellem & van de Ven (2008) for non-one-tenth-integer values with $n \lesssim 3$, and are nevertheless more than $0.2\%$ precise for larger Sérsic indices, for both density and mass profiles. An appendix compares the deprojected Sérsic profiles with those of the popular simple models from Plummer (1911), Jaffe (1983), Hernquist (1990), Navarro et al. (1996), and Einasto (1965).

preprint2020arXiv

Regional analysis of COVID-19 in France from fit of hospital data with different evolutionary models

The SIR evolutionary model predicts too sharp a decrease of the fractions of people infected with COVID-19 in France after the start of the national lockdown, compared to what is observed. I fit the daily hospital data: arrivals in regular and critical care units, releases and deaths, using extended SEIR models. These involve ratios of evolutionary timescales to branching fractions, assumed uniform throughout a country, and the basic reproduction number, $R_0$, before and during the national lockdown, for each region of France. The joint-region Bayesian analysis allows precise evaluations of the time/fraction ratios and pre-hospitalized fractions. The hospital data are well fit by the models, except the arrivals in critical care, which decrease faster than predicted, indicating better treatment over time. Averaged over France, the analysis yields $R_0$= 3.4$\pm$0.1 before the lockdown and 0.65$\pm$0.04 (90% c.l.) during the lockdown, with small regional variations. On 11 May 2020, the Infection Fatality Rate in France was 4 $\pm$1% (90% c.l.), while the Feverish vastly outnumber the Asymptomatic, contrary to the early phases. Without the lockdown nor social distancing, over 2 million deaths from COVID-19 would have occurred throughout France, while a lockdown that would have been enforced 10 days earlier would have led to less than 1000 deaths. The fraction of immunized people reached a plateau below 1% throughout France (3% in Paris) by late April 2020 (95% c.l.), suggesting a lack of herd immunity. The widespread availability of face masks on 11 May, when the lockdown was partially lifted, should keep $R_0$ below unity if at least 46% of the population wear them outside their home. Otherwise, without enhanced other social distancing, a second wave is inevitable and cause the number of deaths to triple between early May and October (if $R_0$=1.2) or even late June (if $R_0$=2).

preprint2020arXiv

The frequency of very young galaxies in the local Universe: II. The view from SDSS spectra

Only a handful of galaxies in the local Universe appear to be very young. We estimate the fraction of very young galaxies (VYGs), defined as those with more than half their stellar masses formed within the last Gyr. We fit non-parametric star formation histories (SFHs) to ~280 000 galaxy spectra from a flux- and volume-limited subsample of the Main Galaxy Sample (MGS) of the SDSS, which is also complete in mass-to-light ratio, thus properly accounting for passive galaxies of a given mass. The VYG fractions decrease with increasing galaxy stellar mass, from ~50% at $m = 10^8\,\rm M_\odot$ to ~0.1% at $m = 10^{11.5}\,\rm M_\odot$, with differences of up to 1 dex between the different spectral models used to estimate the SFH and on how we treat aperture effects. But old stellar populations may hide in our VYGs despite our conservative VYG sample built with galaxies that are globally bluer than within the region viewed by the SDSS fibre. The VYG fractions versus mass decrease more gradually compared to the Tweed et al. predictions using analytical and semi-analytical models of galaxy formation, but agree better with the SIMBA hydrodynamical simulation. These discrepancies highlight the usefulness of VYGs in constraining the strong uncertainties in both galaxy formation models and spectral modelling of galaxy SFHs. Given the lognormal cosmic SFH, these mean VYG fractions suggest that galaxies with $m > 10^8\,\rm M_\odot$ undergo at most 4 major starbursts on average.

preprint2016arXiv

Anisotropic q-Gaussian velocity distributions in LambdaCDM halos

The velocity distribution function (VDF) of dark matter (DM) halos in $Λ$CDM dissipationless cosmological simulations, which must be non-separable in its radial and tangential components, is still poorly known. We present the first single-parameter, non-separable, anisotropic model for the VDF in $Λ$CDM halos, built from an isotropic $q$-Gaussian (Tsallis) VDF of the isotropic set of dimensionless spherical velocity components (after subtraction of streaming motions), normalized by the respective velocity dispersions. We test our VDF on 90 cluster-mass halos of a dissipationless cosmological simulation. Beyond the virial radius, $r_{\rm vir}$, our model VDF adequately reproduces that measured in the simulated halos, but no $q$-Gaussian model can adequately represent the VDF within $r_{\rm vir}$, as the speed distribution function is then flatter-topped than any $q$-Gaussian can allow. Nevertheless, our VDF fits significantly better the simulations than the commonly used Maxwellian (Gaussian) distribution, at virtually all radii within $5\,r_{\rm vir}$. Within 0.4 (1) $r_{\rm vir}$, the non-Gaussianity index $q$ is (roughly) linearly related to the slope of the density profile and also to the velocity anisotropy profile. We provide a parametrization of the modulation of $q$ with radius for both the median fits and the fit of the stacked halo. At radii of a few percent of $r_{\rm vir}$, corresponding to the Solar position in the Milky Way, our best-fit VDF, although fitting better the simulations than the Gaussian one, overproduces significantly the fraction of high velocity objects, indicating that one should not blindly use these $q$-Gaussian fits to make predictions on the direct detection rate of DM particles.

preprint2016arXiv

Dynamical signatures of infall around galaxy clusters: a generalized Jeans equation

We study the internal kinematics of galaxy clusters in the region beyond the sphere of virialization. Galaxies around a virialized cluster are infalling towards the cluster centre with a non-zero mean radial velocity. We develop a new formalism for describing the dynamical state of clusters, by generalizing the standard Jeans formalism with the inclusion of the peculiar infall motions of galaxies and the Hubble expansion as well as the contributions due to background cosmology. Using empirical fits to the radial profiles of density, mean radial velocity and velocity anisotropy of both a stacked cluster-mass halo and two isolated halos of a cosmological dark matter only simulation, we verify that our generalized Jeans equation correctly predicts the radial velocity dispersion out to 4 virial radii. We find that the radial velocity dispersion inferred from the standard Jeans equation is accurate up to 2 virial radii, but overestimated by \approx 20% for the stacked halo and by \approx 40% for the isolated halos, in the range \approx 2-3 virial radii. Our model depends on the logarithmic growth rate of the virial radius (function of halo mass or concentration), which we estimate in 7 different ways, and on the departure from self-similarity of the evolution of the peculiar velocity profile in virial units.

preprint2016arXiv

Outflows Driven by Quasars in High-Redshift Galaxies with Radiation Hydrodynamics

The quasar mode of Active Galactic Nuclei (AGN) in the high-redshift Universe is routinely observed in gas-rich galaxies together with large-scale AGN-driven winds. It is crucial to understand how photons emitted by the central AGN source couple to the ambient interstellar-medium to trigger large-scale outflows. By means of radiation-hydrodynamical simulations of idealised galactic discs, we study the coupling of photons with the multiphase galactic gas, and how it varies with gas cloud sizes, and the radiation bands included in the simulations, which are ultraviolet (UV), optical, and infrared (IR). We show how a quasar with a luminosity of $10^{46}$ erg/s can drive large-scale winds with velocities of $10^2-10^3$ km/s and mass outflow rates around $10^3$ M$_\odot$/yr for times of order a few million years. Infrared radiation is necessary to efficiently transfer momentum to the gas via multi-scattering on dust in dense clouds. However, IR multi-scattering, despite being extremely important at early times, quickly declines as the central gas cloud expands and breaks up, allowing the radiation to escape through low gas density channels. The typical number of multi-scattering events for an IR photon is only about a quarter of the mean optical depth from the center of the cloud. Our models account for the observed outflow rates of $\sim$500-1000 M$_\odot$/yr and high velocities of $\sim 10^3$ km/s, favouring winds that are energy-driven via extremely fast nuclear outflows, interpreted here as being IR-radiatively-driven winds.

preprint2016arXiv

Prospects for annihilating Dark Matter towards Milky Way's dwarf galaxies by the Cherenkov Telescope Array

We derive the Cherenkov Telescope Array (CTA) sensitivity to dark matter (DM) annihilation in several primary channels, over a broad range of DM masses. These sensitivities are estimated when CTA is pointed towards a large sample of Milky Way's dwarf spheroidal galaxies (dSphs) with promising $J$-factors and small statistical uncertainties. This analysis neglects systematic uncertainties, which we estimate at the level of at least 1 dex. We also present sensitivities on the annihilation cross section from a combined analysis of 4 dSphs. We assess the CTA sensitivity by: $i)$ using, for each dSph, a recent determination of the $J$-factor and its statistical error; $ii)$ considering the most up-to-date cosmic ray background; and $iii)$ including both spatial and spectral terms in the likelihood analysis. We find that a joint spectral and spatial analysis improves the CTA sensitivity, in particular for primary channels with sharp features in the $γ$-ray energy spectrum and for dSphs with steep $J$-factor profiles, as deduced from the internal kinematics. The greatest sensitivities are obtained for observations of Ursa Minor among the classical dSphs and of Ursa Major II for ultra-faint dSphs.

preprint2015arXiv

Black hole formation and growth with non-Gaussian primordial density perturbations

Quasars powered by massive black holes (BHs) with mass estimates above a billion solar masses have been identified at redshift 6 and beyond. The existence of such BHs requires almost continuous growth at the Eddington limit for their whole lifetime, of order of one billion years. In this paper, we explore the possibility that positively skewed scale-dependent non-Gaussian primordial fluctuations may ease the assembly of massive BHs. In particular, they produce more low-mass halos at high redshift, thus altering the production of metals and ultra-violet flux, believed to be important factors in BH formation. Additionally, a higher number of progenitors and of nearly equal-mass halo mergers would boost the mass increase provided by BH-BH mergers and merger-driven accretion. We use a set of two cosmological simulations, with either Gaussian or scale-dependent non-Gaussian primordial fluctuations to perform a proof-of-concept experiment to estimate how BH formation and growth are altered. We estimate the BH number density and the fraction of halos where BHs form, for both simulations and for two popular scenarios of BH formation (remnants of the first generation of stars and direct collapse in the absence of metals and molecular hydrogen). We find that the fractions of halos where BHs form are almost identical, but that non-Gaussian primordial perturbations increase the total number density of BHs for the both BH formation scenarios. We also evolve BHs using merger trees extracted from the simulations and find that non-Gaussianities increase both the BH masses and the number of the most massive BHs.

preprint2015arXiv

nIFTy Cosmology: Comparison of Galaxy Formation Models

We present a comparison of 14 galaxy formation models: 12 different semi-analytical models and 2 halo-occupation distribution models for galaxy formation based upon the same cosmological simulation and merger tree information derived from it. The participating codes have proven to be very successful in their own right but they have all been calibrated independently using various observational data sets, stellar models, and merger trees. In this paper we apply them without recalibration and this leads to a wide variety of predictions for the stellar mass function, specific star formation rates, stellar-to- halo mass ratios, and the abundance of orphan galaxies. The scatter is much larger than seen in previous comparison studies primarily because the codes have been used outside of their native environment within which they are well tested and calibrated. The purpose of the `nIFTy comparison of galaxy formation models' is to bring together as many different galaxy formation modellers as possible and to investigate a common approach to model calibration. This paper provides a unified description for all participating models and presents the initial, uncalibrated comparison as a baseline for our future studies where we will develop a common calibration framework and address the extent to which that reduces the scatter in the model predictions seen here.

preprint2015arXiv

Playing with positive feedback: external pressure-triggering of a star-forming disc galaxy

Feedback in massive galaxies generally involves quenching of star formation, a favored candidate being outflows from a central supermassive black hole. At high redshifts however, explanation of the huge rates of star formation often found in galaxies containing AGN may require a more vigorous mode of star formation than attainable by simply enriching the gas content of galaxies in the usual gravitationally-driven mode that is associated with the nearby Universe. Using hydrodynamical simulations, we demonstrate that AGN-pressure-driven star formation potentially provides the positive feedback that may be required to generate the accelerated star formation rates observed in the distant Universe.

preprint2014arXiv

Dark matter in elliptical galaxies: I. Is the total mass density profile of the NFW form or even steeper?

Elliptical galaxies are modelled as Sersic luminosity distributions with density profiles (DPs) for the total mass adopted from the DPs of haloes within dissipationless LambdaCDM N-body simulations. Ellipticals turn out to be inconsistent with cuspy low-concentration NFW models representing the total mass, nor are they consistent with a steeper -1.5 inner slope, nor with the shallower models proposed by Navarro et al. 04, nor with NFW models 10 times more concentrated than predicted, as deduced from several X-ray observations: the mass models, extrapolated inwards, lead to local mass-to-light ratios that are smaller than the stellar value inside an effective radius (R_e), and to central aperture velocity dispersions that are much smaller than observed. This conclusion remains true as long as there is no sharp steepening (slope < -2) of the dark matter (DM) DPs just inside 0.01 virial radii. The too low total mass and velocity dispersion produced within R_e by an NFW-like total mass profile suggests that the stellar component should dominate the DM one out to at least R_e. It should then be difficult to kinematically constrain the inner slope of the dark matter DP of ellipticals. The high concentration parameters deduced from X-ray observations appear to be a consequence of fitting an NFW model to the total mass DP made up of a stellar component that dominates inside and a DM component that dominates outwards. An appendix gives the virial mass dependence of the concentration parameter, central density, and total mass of the Navarro et al. model. In a 2nd appendix are given single integral expressions for the velocity dispersions averaged along the line-of-sight, in circular apertures and in thin slits, for general luminosity density and mass distributions, with isotropic orbits.

preprint2014arXiv

Dark matter in elliptical galaxies: II. Estimating the mass within the virial radius

Elliptical galaxies are modelled with a a 4-component model: Sersic stars, LCDM dark matter (DM), hot gas and central black hole. DM is negligible in the inner regions, which are dominated by stars and the central black hole. This prevents any kinematical estimate (using a Jeans analysis) of the inner slope of the DM density profile. The gas fraction rises, but the baryon fraction decreases with radius, at least out to 10 effective radii (R_e). Even with line-of-sight velocity dispersion (VD) measurements at 4 to 6 R_e with 20 km/s accuracy and perfectly known velocity anisotropy, the total mass within the virial radius (r_v) is uncertain by a factor over 3. The DM distributions found in LCDM simulations are consistent with the stellar VD profiles, but appear inconsistent with the low VDs measured by Romanowsky et al. (2003) of planetary nebulae between 2 and 5 R_e, which imply such low M/Ls that the baryon fraction within r_v must be greater than the universal value. Replacing the NFW DM model by the new model of Navarro et al. (2004) decreases slightly the VD at a given radius. So, given the observed VD measured at 5 R_e, the inferred M/L within r_v is 40% larger than predicted with the NFW model. Folding in the slight (strong) radial anisotropy found in LCDM (merger) simulations, which is well modelled (much better than with the Osipkov-Merritt formula) with beta(r) = 1/2 r/(r+a), the inferred M/L within r_v is another 1.6 (2.4) times higher than for the isotropic NFW model. Thus, the DM model and radial anisotropy can partly explain the low PN VDs, but not in full. In an appendix, single integral expressions are derived for the VDs in terms of the tracer density and total mass profiles, for 3 anisotropic models: radial, Osipkov-Merritt, and the model above, for general radial profiles of luminosity density and mass.

preprint2014arXiv

Effect of primordial non-Gaussianities on the far-UV luminosity function of high-redshift galaxies: implications for cosmic reionization

[Abridged] Understanding how the intergalactic medium (IGM) was reionized at z > 6 is one of the big challenges of current high redshift astronomy. It requires modelling the collapse of the first astrophysical objects (Pop III stars, first galaxies) and their interaction with the IGM, while at the same time pushing current observational facilities to their limits. The observational and theoretical progress of the last few years have led to the emergence of a coherent picture in which the budget of hydrogen-ionizing photons is dominated by low-mass star-forming galaxies, with little contribution from Pop III stars and quasars. The reionization history of the Universe therefore critically depends on the number density of low-mass galaxies at high redshift. In this work, we explore how changes in the statistical properties of initial density fluctuations affect the formation of early galaxies. Following Habouzit et al. (2014), we run 5 N-body simulations with Gaussian and (scale-dependent) non-Gaussian initial conditions, all consistent with Planck constraints. By appealing to a galaxy formation model and to a population synthesis code, we compute the far-UV galaxy luminosity function down to M_UV = -14 at redshift 7 < z < 15. We find that models with strong primordial non-Gaussianities on < Mpc scales show a far-UV luminosity function significantly enhanced in low-mass galaxies. We adopt a reionization model calibrated from state-of-the-art hydrodynamical simulations and show that such non-Gaussianities leave a clear imprint on the Universe reionization history and electron Thomson scattering optical depth tau_E. Although current uncertainties in the physics of reionization and on the determination of tau_E still dominate the signatures of non-Gaussianities, our results suggest that tau_E could ultimately be used to constrain the statistical properties of initial density fluctuations.

preprint2014arXiv

Testing primordial non-Gaussianities on galactic scales at high redshift

Primordial non-Gaussianities provide an important test of inflationary models. Although the Planck CMB experiment has produced strong limits on non-Gaussianity on scales of clusters, there is still room for considerable non-Gaussianity on galactic scales. We have tested the effect of local non-Gaussianity on the high redshift galaxy population by running five cosmological N-body simulations down to z=6.5. For these simulations, we adopt the same initial phases, and either Gaussian or scale-dependent non-Gaussian primordial fluctuations, all consistent with the constraints set by Planck on clusters scales. We then assign stellar masses to each halo using the halo - stellar mass empirical relation of Behroozi et al. (2013). Our simulations with non-Gaussian initial conditions produce halo mass functions that show clear departures from those obtained from the analogous simulations with Gaussian initial conditions at z>~10. We observe a >0.3 dex enhancement of the low-end of the halo mass function, which leads to a similar effect on the galaxy stellar mass function, which should be testable with future galaxy surveys at z>10. As cosmic reionization is thought to be driven by dwarf galaxies at high redshift, our findings may have implications for the reionization history of the Universe.

preprint2013arXiv

MAMPOSSt: Modelling Anisotropy and Mass Profiles of Observed Spherical Systems. I. Gaussian 3D velocities

Mass modelling of spherical systems through internal motions is hampered by the mass/velocity anisotropy (VA) degeneracy inherent in the Jeans equation, as well as the lack of techniques that are both fast and adaptable to realistic systems. A new fast method, called MAMPOSSt, which performs a maximum likelihood fit of the distribution of observed tracers in projected phase space, is developed and thoroughly tested. MAMPOSSt assumes a shape for the gravitational potential, but instead of postulating a shape for the distribution function in terms of energy and angular momentum, or supposing Gaussian line-of-sight velocity distributions, MAMPOSSt assumes a VA profile and a shape for the 3D velocity distribution, here Gaussian. MAMPOSSt requires no binning, differentiation, nor extrapolation of the observables. Tests on cluster-mass haloes from LambdaCDM cosmological simulations show that, with 500 tracers, MAMPOSSt is able to jointly recover the virial radius, tracer scale radius, dark matter scale radius and outer or constant VA with small bias (<10% on scale radii and <2% on the two other quantities) and inefficiencies of 10%, 27%, 48% and 20%, respectively. MAMPOSSt does not perform better when some parameters are frozen, and even worse when the virial radius is set to its true value, which appears to be the consequence of halo triaxiality. The accuracy of MAMPOSSt depends weakly on the adopted interloper removal scheme, including an efficient iterative Bayesian scheme that we introduce here, which can directly obtain the virial radius with as good precision as MAMPOSSt. Our tests show that MAMPOSSt with Gaussian 3D velocities is very competitive with, and up to 1000x faster than other methods. Hence, MAMPOSSt is a very powerful and rapid tool for the mass and anisotropy modeling of systems such as clusters and groups of galaxies, elliptical and dwarf spheroidal galaxies.

preprint2013arXiv

Physical properties underlying observed kinematics of satellite galaxies

We study the kinematics of satellites around isolated galaxies selected from the Sloan Digital Sky Survey (SDSS) spectroscopic catalog. Using a model of the phase-space density previously measured for the halos of LCDM dark matter cosmological simulations, we determine the properties of the halo mass distribution and the orbital anisotropy of the satellites as a function of the colour-based morphological type and the stellar mass of the central host galaxy. We place constraints on the halo mass and the concentration parameter of dark matter and the satellite number density profiles. We obtain a concentration-mass relation for galactic dark matter haloes that is consistent with predictions of a standard LCDM cosmological model. At given halo or stellar mass, red galaxies have more concentrated halos than their blue counterparts. The fraction of dark matter within a few effective radii is minimal for 11.25<log M_star<11.5. The number density profile of the satellites appears to be shallower than of dark matter, with the scale radius typically 60 per cent larger than of dark matter. The orbital anisotropy around red hosts exhibits a mild excess of radial motions, in agreement with the typical anisotropy profiles found in cosmological simulations, whereas blue galaxies are found to be consistent with an isotropic velocity distribution. Our new constraints on the halo masses of galaxies are used to provide analytic approximations of the halo-to-stellar mass relation for red and blue galaxies.

preprint2012arXiv

The Planetary Nebulae Population in the Nuclear Regions of M31: the SAURON view

Following a first study of the central regions of M32 that illustrated the power of integral-field spectroscopy (IFS) in detecting and measuring the [O III]λ5007 emission of PNe against a strong stellar background, we turn to the very nuclear PN population of M31, within 80 pc of its centre. We show that PNe can also be found in the presence of emission from diffuse gas and further illustrate the excellent sensitivity of IFS in detecting extragalactic PNe through a comparison with narrowband images obtained with the Hubble Space Telescope. We find the nuclear PNe population of M31 is only marginally consistent with the generally adopted form of the PNe luminosity function (PNLF). In particular, this is due to a lack of PNe with absolute magnitude M5007 brighter than -3, which would only result from a rather unfortunate draw from such a model PNLF. We suggest that the observed lack of bright PNe in the nuclear regions of M31 is due to a horizontal-branch population that is more tilted toward less massive and hotter He-burning stars, so that its progeny consists mostly of UV-bright stars that fail to climb back up the asymptotic giant branch (AGB) and only of few, if any, bright PNe powered by central post-AGB stars. These results are also consistent with recent reports on a dearth of bright post-AGB stars towards the nucleus of M31, and lend further support to the idea that the metallicity of a stellar population has an impact on the way the horizontal branch is populated and to the loose anticorrelation between the strength of the UV-upturn and the specific number of PNe that is observed in early-type galaxies. Finally, our investigation also serves to stress the importance of considering the same spatial scales when comparing the PNe population of galaxies with the properties of their stellar populations.

preprint2011arXiv

Improved Modeling of the Mass Distribution of Disk Galaxies by the Einasto Halo Model

(Abridged) The analysis of the rotation curves (RCs) of spiral galaxies provides an efficient diagnostic for studying the properties of dark matter halos and their relations with the baryonic material. We have modeled the RCs of galaxies from The HI Nearby Galaxy Survey (THINGS) with the Einasto halo model, which has emerged as the best-fitting model of the halos arising in dissipationless cosmological N-body simulations. We find that the RCs are significantly better fit with the Einasto halo than with either a pseudo-isothermal sphere (Iso) or Navarro-Frenk-White (NFW) halo models. In our best-fit models, the radius of density slope -2 and the density at this radius are highly correlated. The Einasto index, which controls the overall shape of the density profile, is near unity on average for intermediate and low mass halos. This is not in agreement with the predictions from LCDM simulations. The indices of the most massive halos are in rough agreement with those of cosmological simulations and appear correlated with the halo virial mass. We find that a typical Einasto density profile declines more strongly in its outermost parts than any of the Iso or NFW models whereas it is relatively shallow in its innermost regions. The core nature of those regions of halos thus extends the cusp-core controversy found for the NFW model with low surface density galaxies to the Einasto halo with more massive galaxies like those of THINGS. We thus find that the Einasto halo model provides, so far, the best match to the observed RCs, and can therefore be considered as a new standard model for dark matter halos.

preprint2011arXiv

Is the 2MASS clustering dipole convergent?

There is a long-standing controversy about the convergence of the dipole moment of the galaxy angular distribution (the so-called clustering dipole). We study the growth of the clustering dipole of galaxies as a function of the limiting flux of the sample from the Two Micron All Sky Survey (2MASS). Contrary to some earlier claims, we find that the dipole does not converge before the completeness limit of the 2MASS Extended Source Catalog, i.e. up to 13.5 mag in the near-infrared K_s band (equivalent to an effective distance of 300 Mpc/h). We compare the observed growth of the dipole with the theoretically expected, conditional one (i.e., given the velocity of the Local Group relative to the CMB), for the LambdaCDM power spectrum and cosmological parameters constrained by WMAP. The observed growth turns out to be within 1-sigma confidence level of its theoretical counterpart once the proper observational window of the 2MASS flux-limited catalog is included. For a contrast, if the adopted window is a top-hat, then the predicted dipole grows significantly faster and converges to its final value for a distance of about 300 Mpc/h. By comparing the observational windows, we show that for a given flux limit and a corresponding distance limit, the 2MASS flux-weighted window passes less large-scale signal than the top-hat one. We conclude that the growth of the 2MASS dipole for effective distances greater than 200 Mpc/h is only apparent. On the other hand, for a distance of 80 Mpc/h (mean depth of the 2MASS Redshift Survey) and the LambdaCDM power spectrum, the true dipole is expected to reach only ~80% of its final value. Eventually, since for the window function of 2MASS the predicted growth is consistent with the observed one, we can compare the two to evaluate beta = (Omega_m)^0.55 / b. The result is beta = 0.38+-0.04, which leads to an estimate of the density parameter Omega_m = 0.20+-0.08.

preprint2011arXiv

Predicting the frequencies of young and of tiny galaxies

A simple, 1-equation, galaxy formation model is applied to both the halo merger tree derived from a high-resolution dissipationless cosmological simulation and to 1/4 million Monte-Carlo halo merger trees. The galaxy formation model involves a sharp entropy barrier against the accretion of gas onto low-mass halos, the shock heating of infalling gas far from the central regions of massive halos, and supernova feedback that drives the gas out of shallow halo potential wells. With the first approach, we show that the large majority of galaxies within group- and cluster-mass halos, known to be mainly dwarf ellipticals, have acquired the bulk of their stellar mass through gas accretion and not via galaxy mergers. With the second approach, we qualitatively reproduce the downsizing trend of greater ages at greater masses in stars and predict an upsizing trend of greater ages as one proceeds to masses lower than 10^10 M_Sun. We find that the fraction of galaxies with very young stellar populations (more than half the stellar mass formed within the last 1.5 Gyr) is a function of present-day stellar mass, which peaks at 0.5% at m_crit=10^7.5-9.5 M_Sun, roughly corresponding to the masses of blue compact dwarfs. We predict that the stellar mass function of galaxies should not show a maximum at m_stars > 10^{5.5}, M_Sun, with a power-law stellar mass function with slope approx -1.6 if the IGM temperature in the outskirts of halos before reionization is set by H2 cooling. We speculate on the nature of the lowest mass galaxies.

preprint2011arXiv

The velocity modulation of galaxy properties in and near clusters: quantifying the decrease in star formation in backsplash galaxies

The efficiency of recent star formation (SF) in galaxies increases with increasing projected distance from the centre of a cluster out to several times its virial radius (R_v). Using a complete sample of galaxies in 268 clusters from the SDSS DR4, we investigate how, at a given projected radius from the cluster centre, M* and SF properties of a galaxy depend on its absolute line-of-sight velocity in the cluster rest frame, |v_LOS|. We find that for R<0.5 R_v, the fraction of high mass non-BCG galaxies increases towards the centre for low |v_LOS|. At a given projected radius, the fraction of Galaxies with Ongoing or Recent (<1-3 Gyr) Efficient Star Formation (GORES, with EW(H_delta)>2 ang & D_4000>1.5) is slightly but significantly lower for low |v_LOS| galaxies than for their high velocity counterparts. We study these observational trends with the help of a dark matter (DM) cosmological simulation. We find that the backsplash particles account for at least one-third (half) of all particles at projected radii slightly greater than the virial radius and |v_LOS|<sigma_v. The deprojection of the GORES fraction leads to a saturated linear increase with radius. We fit simple models of the fraction of GORES as a function of class only or class and distance to the cluster centre (as in our deprojected fraction). In our best-fitting model GORES account for 13% of galaxies within the virial sphere, 11% of the virial population, 34% of the distant (for projected radii R<2 R_v) infall population and 19% of the backsplash galaxies. Given the 1-3 Gyr lookback time of our GORES indicators, these results suggest that SF in a galaxy is almost completely quenched in a single passage through the cluster.

preprint2010arXiv

Is the 2MASS dipole convergent?

We study the growth of the clustering dipole of galaxies from the Two Micron All Sky Survey (2MASS). We find that the dipole does not converge before the completeness limit of the 2MASS Extended Source Catalog, i.e. up to about 300 Mpc/h. We compare the observed growth of the dipole with the theoretically expected, conditional growth for the LambdaCDM power spectrum and cosmological parameters constrained by WMAP. The observed growth turns out to be within 1-sigma confidence level of the theoretical one, once the proper observational window of the 2MASS flux dipole is included. For a contrast, if the adopted window is a top hat, then the predicted dipole grows significantly faster and converges to its final value at a distance of about 200 Mpc/h. We study the difference between the top-hat window and the window for the flux-limited 2MASS survey and we conclude that the growth of the 2MASS dipole at effective distances greater than 200 Mpc/h is only apparent. Eventually, since for the window function of 2MASS the predicted growth is consistent with the observed one, we can compare the two to evaluate beta = (Omega_m)^0.55 / b. The result is beta = 0.38+-0.05, which gives a rough estimate of Omega_m = 0.2+-0.1.

preprint2010arXiv

The universal distribution of halo interlopers in projected phase space. Bias in galaxy cluster concentration and velocity anisotropy?

When clusters of galaxies are viewed in projection, one cannot avoid picking up foreground/background interlopers (FBIs), that lie within the virial cone (VC), but outside the virial sphere. Structural & kinematic deprojection equations are not known for an expanding Universe, where the Hubble flow (HF) stretches the line-of-sight (LOS) distribution of velocities. We analyze 93 mock relaxed clusters, built from a cosmological simulation. The stacked mock cluster is well fit by an m=5 Einasto DM density profile (but only out to 1.5 virial radii [r_v]), with velocity anisotropy (VA) close to the Mamon-Lokas model with VA radius equal to that of density slope -2. The surface density of FBIs is nearly flat out to r_v, while their LOS velocity distribution shows a dominant gaussian cluster-outskirts component and a flat field component. This distribution of FBIs in projected phase space is nearly universal in mass. A local k=2.7 sigma velocity cut returns the LOS velocity dispersion profile (LOSVDP) expected from the NFW density and VA profiles measured in 3D. The HF causes a shallower outer LOSVDP that cannot be well matched by the Einasto model for any k. After this velocity cut, FBIs still account for 23% of DM particles within the VC (close to the observed fraction of cluster galaxies lying off the Red Sequence). The best-fit projected NFW/Einasto models underestimate the 3D concentration by 6+/-6% (16+/-7%) after (before) the velocity cut, unless a constant background is included in the fit. Assuming the correct mass profile, the VA profile is well recovered from the measured LOSVDP, with a slight bias towards more radial orbits in the outer regions. These small biases are overshadowed by large cluster-cluster variations caused by cosmic variance. An appendix provides an analytical approximation to the surface density, projected mass and tangential shear profiles of the Einasto model.

preprint2009arXiv

Kinematic deprojection and mass inversion of spherical systems of known velocity anisotropy

Traditionally, the mass / velocity anisotropy degeneracy (MAD) inherent in the spherical, stationary, non-streaming Jeans equation has been handled by assuming a mass profile and fitting models to the observed kinematical data. Here, the opposite approach is considered: the equation of anisotropic kinematic projection is inverted for known arbitrary anisotropy to yield the space radial velocity dispersion profile in terms of an integral involving the radial profiles of anisotropy and isotropic dynamical pressure. Then, through the Jeans equation, the mass profile is derived in terms of double integrals of observable quantities. Single integral formulas for both deprojection and mass inversion are provided for several simple anisotropy models (isotropic, radial, circular, general constant, Osipkov-Merritt, Mamon-Lokas and Diemand-Moore-Stadel). Tests of the mass inversion on NFW models with these anisotropy models yield accurate results in the case of perfect observational data, and typically better than 70% (in 4 cases out of 5) accurate mass profiles for the sampling errors expected from current observational data on clusters of galaxies. For the NFW model with mildly increasing radial anisotropy, the mass is found to be insensitive to the adopted anisotropy profile at 7 scale radii and to the adopted anisotropy radius at 3 scale radii. This anisotropic mass inversion method is a useful complementary tool to analyze the mass and anisotropy profiles of spherical systems. It provides the practical means to lift the MAD in quasi-spherical systems such as globular clusters, round dwarf spheroidal and elliptical galaxies, as well as groups and clusters of galaxies, when the anisotropy of the tracer is expected to be linearly related to the slope of its density.

preprint2009arXiv

Suppressed radio emission in supercluster galaxies: enhanced ram pressure in merging clusters?

The environmental influence on the 1.4 GHz continuum radio emission of galaxies is analyzed in a 600 deg2 region of the local Universe containing the Shapley Supercluster (SSC). Galaxies in the FLASH and 6dFGS redshift surveys are cross-identified with NVSS radio sources, selected in a subsample doubly complete in volume and luminosity. Environmental effects are studied through a smoothed density field (normalized with random catalogs with the same survey edges and redshift selection function) and the distance to the nearest cluster (R/r200, where r200 is the virial radius, whose relation to the aperture velocity dispersion is quantified). The fraction of high radio loudness (R_K=L_radio/L_K) galaxies in the 10 Mpc Abell 3558 cluster complex at the core of the SSC (SSC-CR) is half as large than elsewhere. In the SSC-CR, R_K is anti-correlated with the density of the large-scale environment and correlated with R/r200: central brightest cluster galaxies (BCGs) in the SSC-CR are 10x less radio-loud than BCGs elsewhere, with signs of suppressed radio loudness in the SSC-CR also present beyond the BCGs, out to at least 0.3 r200. This correlation is nearly as strong as the tight correlation of L_K with R/r200 (K-luminosity segregation), inside the SSC-CR. The suppression of radio loudness in SSC-CR BCGs can be attributed to cluster-cluster mergers that destroy the cool core and thus the supply of gas to the central AGN. We analytically demonstrate that the low radio loudness of non-BCG galaxies within SSC-CR clusters cannot be explained by direct major galaxy mergers or rapid galaxy flyby collisions, but by the loss of gas supply through the enhanced ram pressure felt when these galaxies cross the shock front between the 2 merging clusters and are later subjected to the stronger wind from the 2nd cluster.

preprint2009arXiv

The mass and anisotropy profiles of galaxy clusters from the projected phase space density: testing the method on simulated data

We present a new method of constraining the mass and velocity anisotropy profiles of galaxy clusters from kinematic data. The method is based on a model of the phase space density which allows the anisotropy to vary with radius between two asymptotic values. The characteristic scale of transition between these asymptotes is fixed and tuned to a typical anisotropy profile resulting from cosmological simulations. The model is parametrized by two values of anisotropy, at the centre of the cluster and at infinity, and two parameters of the NFW density profile, the scale radius and the scale mass. In order to test the performance of the method in reconstructing the true cluster parameters we analyze mock kinematic data for 20 relaxed galaxy clusters generated from a cosmological simulation of the standard LCDM model. We use Bayesian methods of inference and the analysis is carried out following the Markov Chain Monte Carlo approach. The parameters of the mass profile are reproduced quite well, but we note that the mass is typically underestimated by 15 percent, probably due to the presence of small velocity substructures. The constraints on the anisotropy profile for a single cluster are in general barely conclusive. Although the central asymptotic value is determined accurately, the outer one is subject to significant systematic errors caused by substructures at large clustercentric distance. The anisotropy profile is much better constrained if one performs joint analysis of at least a few clusters. In this case it is possible to reproduce the radial variation of the anisotropy over two decades in radius inside the virial sphere.

preprint2009arXiv

The orientation and kinematics of inner tidal tails around dwarf galaxies orbiting the Milky Way

Using high-resolution collisionless N-body simulations we study the properties of tidal tails formed in the immediate vicinity of a two-component dwarf galaxy evolving in a static potential of the Milky Way (MW). The stellar component of the dwarf is initially in the form of a disk and the galaxy is placed on an eccentric orbit motivated by CDM-based cosmological simulations. We measure the orientation, density and velocity distribution of the stars in the tails. Due to the geometry of the orbit, in the vicinity of the dwarf, where the tails are densest and therefore most likely to be detectable, they are typically oriented towards the MW and not along the orbit. We report on an interesting phenomenon of `tidal tail flipping': on the way from the pericentre to the apocentre the old tails following the orbit are dissolved and new ones pointing towards the MW are formed over a short timescale. We also find a tight linear relation between the velocity of stars in the tidal tails and their distance from the dwarf. Using mock data sets we demonstrate that if dwarf spheroidal (dSph) galaxies in the vicinity of the MW are tidally affected their kinematic samples are very likely contaminated by tidally stripped stars which tend to artificially inflate the measured velocity dispersion. The effect is stronger for dwarfs on their way from the peri- to the apocentre due to the formation of new tidal tails after pericentre. Realistic mass estimates of dSph galaxies thus require removal of these stars from kinematic samples.

preprint2009arXiv

Tidal evolution of disky dwarf galaxies in the Milky Way potential: the formation of dwarf spheroidals

We conduct high-resolution collisionless N-body simulations to investigate the tidal evolution of dwarf galaxies on an eccentric orbit in the Milky Way (MW) potential. The dwarfs originally consist of a low surface brightness stellar disk embedded in a cosmologically motivated dark matter halo. During 10 Gyr of dynamical evolution and after 5 pericentre passages the dwarfs suffer substantial mass loss and their stellar component undergoes a major morphological transformation from a disk to a bar and finally to a spheroid. The bar is preserved for most of the time as the angular momentum is transferred outside the galaxy. A dwarf spheroidal (dSph) galaxy is formed via gradual shortening of the bar. This work thus provides a comprehensive quantitative explanation of a potentially crucial morphological transformation mechanism for dwarf galaxies that operates in groups as well as in clusters. We compare three cases with different initial inclinations of the disk and find that the evolution is fastest when the disk is coplanar with the orbit. Despite the strong tidal perturbations and mass loss the dwarfs remain dark matter dominated. For most of the time the 1D stellar velocity dispersion, σ, follows the maximum circular velocity, V_{\rm max}, and they are both good tracers of the bound mass. Specifically, we find that M_{\rm bound} \propto V_{\rm max}^{3.5} and V_{\rm max} \sim \sqrt{3} σin agreement with earlier studies based on pure dark matter simulations. The latter relation is based on directly measuring the stellar kinematics of the simulated dwarf and may thus be reliably used to map the observed stellar velocity dispersions of dSphs to halo circular velocities when addressing the missing satellites problem.

preprint2007arXiv

SDSSJ124155.33+114003.7 -- a Missing Link Between Compact Elliptical and Ultracompact Dwarf Galaxies

We report the discovery of a compact object (R_e = 32 pc, M_B = -12.34 mag) at a projected distance of 9 kpc from Messier 59, a giant elliptical in the Virgo cluster. Using HST imaging and SDSS spectroscopy, both available in the Virtual Observatory, we find that this object has a blue core containing one-quarter of the light, and a redder n=1 Sersic envelope, as well as luminosity-weighted age of 9.3 +/- 1.4 Gyr, a metallicity of -0.03 +/- 0.04 dex and a velocity dispersion of 48 +/- 5 km/s. While ultra-compact dwarfs (UCDs) in the face-on view of the Fundamental Plane are found to form a sequence connecting the highest-luminosity globular clusters with the lowest-luminosity dwarf ellipticals, the compact object near M59 lies in between this UCD sequence and the positions of compact ellipticals. Its stellar age, metallicity, and effective surface brightness are similar to low-luminosity ellipticals and lenticulars, suggesting that SDSSJ124155.33+114003.7 is a result of the tidal stripping of such an object.

preprint2006arXiv

The evolution of galaxy groups and of galaxies therein

Properties of groups of galaxies depend sensitively on the algorithm for group selection, and even the most recent catalogs of groups built from redshift-space selection should suffer from projections and infalling galaxies. The cosmo-dynamical evolution of groups from initial Hubble expansion to collapse and virialization leads to a fundamental track (FT) in virial-theorem-M/L vs crossing time. The increased rates of mergers, both direct and after dynamical friction, in groups relative to clusters, explain the higher fraction of elliptical galaxies at given local number density in X-ray selected groups, relative to clusters, even when the hierarchical evolution of groups is considered. Galaxies falling into groups and clusters should later travel outwards to typically 2 virial radii, which is somewhat less than the outermost radius where observed galaxy star formation efficiencies are enhanced relative to field galaxies of same morphological type. An ongoing analysis of the internal kinematics of X-ray selected groups suggests that the radial profiles of line of sight velocity dispersion are consistent with isotropic NFW distributions for the total mass density, with higher (lower) concentrations than LambdaCDM predictions in groups of high (low) mass. The critical mass, at M200 ~ 10^13 M_sun is consistent with possible breaks in the X-ray luminosity-temperature and Fundamental Plane relations. The internal kinematics of groups indicate that the M-T relation of groups should agree with that extrapolated from clusters with no break at the group scale. The analyses of observed velocity dispersion profiles and of the FT both suggest that low velocity dispersion groups (compact and loose, X-ray emitting or undetected) are quite contaminated by chance projections.

preprint2003arXiv

The Selection of Tenured Astronomers in France

The organization of the recruitment of tenured astronomers and astrophysicists in France is presented and compared with the setup in other countries. The ages for getting tenure have increased from 27-28 in 1980 to 31 today. Foreign scientists constitute at least 11% of the recruits and the delay in their hiring is quantified. The large reliance on national tenure committees is justified, while the increased targeting of positions is questioned and a compromise proposed.

preprint2000arXiv

Tangential Large Scale Structure as a Standard Ruler: Curvature Parameters from Quasars

Several observational analyses suggest that matter is spatially structured at $\approx 130h^{-1}Mpc$ at low redshifts. This peak in the power spectrum provides a standard ruler in comoving space which can be used to compare the local geometry at high and low redshifts, thereby constraining the curvature parameters. It is shown here that this power spectrum peak is present in the observed quasar distribution at $z\sim 2$: qualitatively, via wedge diagrams which clearly show a void-like structure, and quantitatively, via one-dimensional Fourier analysis of the quasars' tangential distribution. The sample studied here contains 812 quasars. The method produces strong constraints (68% confidence limits) on the density parameter $Ω_0$ and weaker constraints on the cosmological constant $λ_0$, which can be expressed by the relation $ Ω_0 = (0.24 \pm0.15) + (0.10\pm0.08) λ_0 $. Independently of $λ_0$ (in the range $λ_0 \in [0,1]$), the constraint is $0.1 < Ω_0 < 0.45$. Combination of the present results with SN Type Ia results yields $Ω_0 = (0.30\pm0.11) + (0.57\pm0.11) (λ_0-0.7),$ $0.55 < λ_0 < 0.95,$ (68% confidence limits). This strongly supports the possibility that the observable universe satisfies a nearly flat, perturbed Friedmann-Lema\^ıtre-Robertson-Walker model, independently of any cosmic microwave background observations.

preprint1995arXiv

The DENIS 2 Micron Survey and its Cosmological Applications

Near-infrared surveys have very promising cosmological applications. The DENIS 2 micron survey, currently in camera-testing, is scheduled to begin in early 1995, and to map the entire southern sky in the $I$, $J$, and $K'$ bands. The steps for the near-homogeneous extraction of $\sim 10^5$ galaxies are outlined, some pitfalls are given, and the estimated limiting capabilities of the survey are presented. A preliminary $J$-band image illustrates these capabilities and some of the hidden pitfalls. The necessary spectroscopic followup is shown to be feasible on a reasonably short timescale.