Search

Astrometry is one of the oldest branches of astronomy which measures the position, the proper motion and parallax of celestial objects. Following the Hipparcos and Gaia missions that have measured several billions of them using global astrometry, we propose to increase astrometry precision on pointed objects using differential astrometry in a large field in order to unravel rocky planets in habitable zones of stars in the Sun vicinity and investigate the nature of dark matter in galactic environments as recommended by the ESA Senior Committee in the Voyager 2050 prospective. Substantial technology developments in a number of critical areas is needed in order to reach the highest required precision of sub-micro-arcsecond. One of them is CMOS image sensors using the stitching technique to merge the multiple design structures on the wafer and produce array with very large number of pixels. Another one is to calibrate the pixel positions using projecting modulating interferometric laser fringes on the array. Finally, the distortion of the optical system can be monitored and compensated using reference stars as metrology sources. The final precision depends on the diameter and the field of view of the telescope that is used as well as the time spent on each target. We present here the science goals that can be achieved with such missions either within the framework of an ESA Medium-class mission or even in the NASA most challenging Habitable Worlds Observatory, a large space telescope recommended by the American Astronomy and Astrophysics prospective for the 2020s and designed specifically to search for signs of life on planets orbiting other stars.

Dwarf spheroidal galaxies are known to be dominated by old stellar populations. This has led to the assumption that their gas-rich progenitors lost their gas during their infall in the Milky Way (MW) halo at distant look-back times. Here, we report a discovery of a tiny but robustly detected population of possibly young ($\sim$ 1 Gyr old) and intermediate-mass ($\rm 1.8 M_{\odot} \le M < 3 M_{\odot}$) stars in MW dwarf spheroidal galaxies. This was established on the basis of their positions in color-magnitude diagrams, after filtering out the bulk of the foreground MW using Gaia DR3 proper motions. We have considered the possibility that this population is made of evolved blue stragglers. For Sculptor, it seems unlikely, because 95.5% of its stars are older than 8 Gyr, leading to masses smaller than 0.9 M$_{\odot}$. This would only allow blue straggler masses of less than 1.8 M$_{\odot}$, which is much lower than what we observed. Alternatively, it would require the merger of three turnoff stars, which appears even more unlikely. On the other hand, the recent Gaia proper motion measurements of MW dwarf galaxies infer their low binding energies and large angular momenta, pointing to a more recent, $\le$ 3 Gyr, infall. Although the nature of the newly discovered stars still needs further confirmation, we find that they are consistent with the recent infall of the dwarf galaxies into the MW halo, when star formation occurred from the ram pressurization of their gas content before its removal by the hot Galactic corona. The abundance of this plausibly young population of stars is similar to the expectations drawn from hydrodynamical simulations. These results point to a novel origin for MW dwarf spheroidal galaxies., Comment: A&A accepted. 20 pages, 11 figures

It is often believed that isolated, compact groups of galaxies (CGs) are special systems. But CGs have rarely been compared to regular groups, which are selected very differently. We study the global properties of a subsample of 80 groups of four members (CG4s), as well as their correlations, within the HMCG Hickson-like sample of compact groups, and compare them with those of the regular groups of Lim et al.. The latter are split into three control samples, complete in redshift and magnitude: one with the four brightest members, one with the four closest members to the brightest group galaxy (BGG) with less than three magnitudes in range (Control4Cs), and one with exactly four members (RG4s). The vast majority of CG4s are located within regular groups, and a large preponderance of the BGGs of these CG4s are the same as those of their host groups. CG4s are smaller than the groups of all other samples and more luminous than RG4s, both results as expected from their selection as high surface brightness systems. However, CG4s have similar luminosities as Control4Cs. CG4s also have higher velocity dispersions, probably because of a too permissive redshift accordance criterion. The BGGs of CG4s are not more dominant in luminosity than those of RG4s, but they are significantly more offset relative to the group size, because the Lim groups are built around their BGGs. In summary, compact groups have similar properties to the regular groups of four galaxies and to the cores of regular groups, once selection criteria of CGs are considered. A large fraction of CGs are the cores of regular groups, isolated on the sky by construction, but rarely isolated in real space (from simulations), indicating that they are often plagued by chance alignments of host groups galaxies along the line of sight., Comment: Submitted to A&A. 14 pages. Comments welcome! Latest version fixes 3rd author's name in metadata

We analyze four epochs of HST imaging over 18 years for the Draco dwarf spheroidal galaxy. We measure precise proper motions (PMs) for hundreds of stars and combine these with existing line-of-sight (LOS) velocities. This provides the first radially-resolved 3D velocity dispersion profiles for any dwarf galaxy. These constrain the intrinsic velocity anisotropy and resolve the mass-anisotropy degeneracy. We solve the Jeans equations in oblate axisymmetric geometry to infer the mass profile. We find the velocity dispersion to be radially anisotropic along the symmetry axis and tangentially anisotropic in the equatorial plane, with a globally-averaged value $\overline{\beta_{\mathrm B}}=-0.20^{+ 0.28}_{- 0.53}$, (where $1 - \beta_{\mathrm B} \equiv \langle v_{\mathrm{ tan}}^2 \rangle / \langle v_{\mathrm{ rad}}^2 \rangle$ in 3D). The logarithmic dark matter (DM) density slope over the observed radial range, $\Gamma_{\mathrm{ dark}}$, is $-0.83^{+ 0.32}_{- 0.37}$, consistent with the inner cusp predicted in $\Lambda$CDM cosmology. As expected given Draco's low mass and ancient star formation history, it does not appear to have been dissolved by baryonic processes. We rule out cores larger than 487, 717, 942 pc at respective 1-, 2-, 3-$\sigma$ confidence, thus imposing important constraints on the self-interacting DM cross-section. Spherical models yield biased estimates for both the velocity anisotropy and the inferred slope. The circular velocity at our outermost data point (900 pc) is $24.19^{+ 6.31}_{- 2.97} \ \mathrm{km~s^{-1}}s$. We infer a dynamical distance of $75.37^{+ 4.73}_{- 4.00}$ kpc, and show that Draco has a modest LOS rotation, with $\left = 0.22 \pm 0.09$. Our results provide a new stringent test of the so-called `cusp-core' problem that can be readily extended to other dwarfs., Comment: 35 pages, 18 figures, 3 Tables. Accepted for publication in ApJ. Journal version has better readability. Data is available at Zenodo: https://zenodo.org/records/11111113

Nearby dwarf galaxies display a variety of effective radii (sizes) at a given stellar mass, suggesting different evolution scenarios according to their final "stellar" size. The TNG hydrodynamical simulations present a bimodality in the z = 0 size - mass relation (SMRz0) of dwarf galaxies, at $r_{1/2,\star}$ ~ 450 pc. Using the TNG50 simulation, we explored the evolution of the most massive progenitors of dwarf galaxies (z=0 $\log( M_\star / \mathrm{M}_\odot)$ between 8.4 and 9.2) that end up as central galaxies of their groups. We split these dwarfs into three classes of the SMRz0: "Normals" from the central spine of the main branch, and "Compacts" from the secondary branch as well as the lower envelope of the main branch. Both classes of Compacts see their stellar sizes decrease from z ~ 1 onwards in contrast to Normals, while the sizes of the gas and dark matter (DM) components continue to increase (as for Normals). A detailed analysis reveals that Compacts live in poorer environments, and thus suffer fewer major mergers from z = 0.8 onwards, which otherwise would pump angular momentum into the gas, allowing strong gas inflows, producing inner star formation, and thus leading to the buildup of a stellar core. Compacts are predicted to be rounder and to have bluer cores. Compact dwarfs of similar sizes are observed in the GAMA survey, but the bimodality in size is less evident and the most compact dwarfs tend to be passive rather than star forming, as in TNG50. Our conclusions should therefore be confirmed with future cosmological hydrodynamical simulations., Comment: 15 pages, 15 figures, 3 tables, Accepted for Publication in A&A. Final version

Compact groups of dwarf galaxies (CGDs) have been observed at low redshifts ($z<0.1$) and are direct evidence of hierarchical assembly at low masses. To understand the formation of CGDs and the galaxy assembly in the low-mass regime, we search for analogues of compact (radius $\leq 100$ kpc) groups of dwarfs ($7 \leq \log[M_{\ast}/{\rm M}_\odot] \leq 9.5$) in the IllustrisTNG highest-resolution simulation. Our analysis shows that TNG50-1 can successfully produce CGDs at $z=0$ with realistic total and stellar masses. We also find that the CGD number density decreases towards the present, especially at $z \lesssim 0.26$, reaching $n \approx 10^{-3.5}$ $\rm cMpc^{-3}$ at $z = 0$. This prediction can be tested observationally with upcoming surveys targeting the faint end of the galaxy population and is essential to constrain galaxy evolution models in the dwarf regime. The majority of simulated groups at $z \sim 0$ formed recently ($\lesssim 1.5 \ \rm Gyr$), and CGDs identified at $z \leq 0.5$ commonly take more than 1 Gyr to merge completely, giving origin to low- to intermediate-mass ($8 \leq \log[M_{\ast}/{\rm M}_\odot] \leq 10$) normally star-forming galaxies at $z=0$. We find that halos hosting CGDs at $z = 0$ formed later when compared to halos of similar mass, having lower stellar masses and higher total gas fractions. The simulations suggest that CGDs observed at $z \sim 0$ arise from a late hierarchical assembly in the last $\sim 3$ Gyr, producing rapid growth in total mass relative to stellar mass and creating dwarf groups with median halo masses of $\sim 10^{11.3}$ $\rm M_\odot$ and B-band mass-to-light ratios mostly in the range $10 \lesssim M/L \lesssim 100$, in agreement with previous theoretical and observational studies., Comment: 20 pages, 18 figures, 4 tables, accepted by MNRAS. Updated to reflect minor changes made during the referring process. Nearest neighbour catalogues available at http://www.tng-project.org/floresfreitas24

Most Milky Way dwarf galaxies are much less bound to their host than are relics of Gaia-Sausage-Enceladus and Sgr. These dwarfs are expected to have fallen into the Galactic halo less than 3 Gyr ago, and will therefore have undergone no more than one full orbit. Here, we have performed hydrodynamical simulations of this process, assuming that their progenitors are gas-rich, rotation-supported dwarfs. We follow their transformation through interactions with the hot corona and gravitational field of the Galaxy. Our dedicated simulations reproduce the structural properties of three dwarf galaxies: Sculptor, Antlia II and, with somewhat a lower accuracy, Crater II. This includes reproducing their large velocity dispersions, which are caused by ram-pressure stripping and Galactic tidal shocks. Differences between dwarfs can be interpreted as due to different orbital paths, as well as to different initial conditions for their progenitor gas and stellar contents. However, we failed to suppress in a single orbit the rotational support of our Sculptor analog if it is fully dark-matter dominated. In addition, we have found that classical dwarf galaxies like Sculptor may have stellar cores sufficiently dense to survive the pericenter passage through adiabatic contraction. On the contrary, our Antlia II and Crater II analogs are tidally stripped, explaining their large sizes, extremely low surface brightnesses, and velocity dispersion. This modeling explains differences between dwarf galaxies by reproducing them as being at different stages of out-of-equilibrium stellar systems., Comment: Paper III, 15 pages, 9 Figures, 2 Tables, accepted by MNRAS 23rd November, see also Paper II (arXiv:2311.05677), please see the video of a dwarf (Sculptor) transformation after its first infall in the MW halo at: https://www.youtube.com/watch?v=SwxSdmfQis4

We study how structural properties of globular clusters and dwarf galaxies are linked to their orbits in the Milky Way halo. From the inner to the outer halo, orbital energy increases and stellar-systems gradually move out of internal equilibrium: in the inner halo, high-surface brightness globular clusters are at pseudo-equilibrium, while further away, low-surface brightness clusters and dwarfs appear more tidally disturbed. Dwarf galaxies are the latest to arrive into the halo as indicated by their large orbital energies and pericenters, and have no time for more than one orbit. Their (gas-rich) progenitors likely lost their gas during their recent arrival in the Galactic halo. If dwarfs are at equilibrium with their dark matter (DM) content, the DM density should anti-correlate with pericenter. However, the transformation of DM dominated dwarfs from gas-rich rotation-supported into gas-poor dispersion-supported systems is unlikely accomplished during a single orbit. We suggest instead that the above anti-correlation is brought by the combination of ram-pressure stripping and of Galactic tidal shocks. Recent gas removal leads to an expansion of their stellar content caused by the associated gravity loss, making them sufficiently fragile to be transformed near pericenter passage. Out of equilibrium dwarfs would explain the observed anti-correlation of kinematics-based DM density with pericenter without invoking DM density itself, questioning its previous estimates. Ram-pressure stripping and tidal shocks may contribute to the dwarf velocity dispersion excess. It predicts the presence of numerous stars in their outskirts and a few young stars in their cores., Comment: Paper II, 18 pages, 15 Figures, MNRAS 527, 2178-2733 (2024), see also paper III with simulations (arXiv:2311.05687), as well as a video showing the transformation of a dwarf recently falling into the MW halo at: https://www.youtube.com/watch?v=SwxSdmfQis4

Using over 3000 compact groups (CGs) of galaxies extracted from mock catalogues built from semi-analytical models of galaxy formation (SAMs), we study whether the CG assembly channel affects the z=0 properties of galaxies and their evolution. The evolution of CG galaxy properties with time is a clear function of their stellar masses. For instance, high-stellar-mass CG galaxies have lived their last 8 Gyr with little cold gas content while maintaining their reservoir of hot gas, while low-mass CG galaxies still preserve some of their cold gas content at the present but they have completely drained their hot gas reservoir. Beyond that, we find that the evolution of CG galaxies is also a function of the assembly history of the CGs: with more extreme losses of gas content, faster mass gain rates for black holes and more marked suppression of star formation as a function of cosmic time as we go from recent to early CG assembly. Thus, CGs constitute another laboratory for galaxy assembly bias, as the later assembling groups have later star formation. Our results show that classifying CGs according to their assembly channel is a way of distinguishing different paths by which galaxies transform their properties throughout their history., Comment: 19 pages, 10 figures, 4 table, accepted for publication in MNRAS

Recent studies of nearby globular clusters have discovered excess dark mass in their cores, apparently in an extended distribution, and simulations indicate that this mass is composed mostly of white dwarfs (respectively stellar-mass black holes) in clusters that are core-collapsed (respectively with a flatter core). We perform mass-anisotropy modelling of the closest globular cluster, M4, with intermediate slope for the inner stellar density. We use proper-motion data from Gaia EDR3 and from observations by the Hubble Space Telescope. We extract the mass profile employing Bayesian Jeans modelling, and check our fits with realistic mock data. Our analyses return isotropic motions in the cluster core and tangential motions ($\beta\approx -0.4$$\pm$$0.1$) in the outskirts. We also robustly measure a dark central mass of roughly $800\pm300 \,$M$_{\odot}$, but it is not possible to distinguish between a point-like source, such as an intermediate-mass black hole (IMBH), or a dark population of stellar remnants of extent $\approx 0.016\,\rm pc \simeq 3300\,AU$. However, when removing a high-velocity star from the cluster centre, the same mass excess is found, but more extended ($\sim 0.034\, \rm{pc} \approx 7000\,\rm AU$). We use Monte Carlo $N$-body models of M4 to interpret the second outcome, and find that our excess mass is not sufficiently extended to be confidently associated with a dark population of remnants. Finally, we discuss the feasibility of these two scenarios (i.e., IMBH vs. remnants), and propose new observations that could help to better grasp the complex dynamics in M4's core., Comment: 19 page, 15 figures, 3 tables. Accepted for publication in MNRAS

Halo inhabitants are individual stars, stellar streams, star and globular clusters, and dwarf galaxies. Here we compare the two last categories that include objects of similar stellar mass, which are often studied as self-dynamical equilibrium systems. We discover that the half-light radius of globular clusters depends on their orbital pericenter and total energy, and that Milky Way (MW) tides may explain the observed correlation. We also suggest that the accretion epoch of stellar systems in the MW halo can be calibrated by the total orbital energy, and that such a relation is due to both the mass growth of the MW and dynamical friction affecting mostly satellites with numerous orbits. This calibration starts from the bulge, to Kraken, Gaia Sausage Enceladus, Sagittarius stellar systems, and finally to the new coming dwarfs, either or not linked to the vast-polar structure. The most eccentric globular clusters and dwarfs have their half-light radius scaling as the inverse of their binding energy, and this over more than two decades. This means that earlier arriving satellites are smaller due to the tidal effects of the MW. Therefore, most halo inhabitants appear to have their structural parameters shaped by MW tides and also by ram-pressure for the most recent arrivals, the dwarf galaxies. The correlations found in this study can be used as tools to further investigate the origin of globular clusters and dwarfs, as well as the assembly history of our Galaxy., Comment: 18 pages, 8 Figures, 2 Tables with orbital properties of globular clusters, accepted in MNRAS, December the 12, 2022, version similar to the printed version

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., Comment: SPIE Conference : "Space Telescopes and Instrumentation 2022: Optical, Infrared, and Millimeter Wave", 17-22 July 2022, Montr\'eal

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., Comment: 28 pages, 16 figures, 4 tables, accepted for publication in MNRAS

Many catalogues of isolated compact groups of galaxies (CGs) have been extracted using Hickson's criteria to identify isolated, dense systems of galaxies, with at least three or four galaxies concordant in magnitude and redshift. But is not clear to what extent the catalogues of CGs are complete and reliable, relative to 3D truly isolated, dense groups. Using five different semi-analytical models of galaxy formation (SAMs), we identify isolated dense groups in 3D real space, containing at least three galaxies. We then build mock redshift space galaxy catalogues and run a Hickson-like CG finder. We find that the Hickson-like algorithm in redshift space is poor at recovering 3D CGs of at least 3 galaxies, with a purity of $\sim 10\%$ and a completeness of $\sim 22\%$. Among the $\sim 90\%$ of spurious systems, typically $60\%$ are dense structures that failed the 3D isolation criteria, while the remaining $40\%$ are chance alignments of galaxies along the line of sight, nearly all of which are within regular groups, with some variation with the SAM used for the analysis. In other words, while only $10\%$ of CGs are isolated dense groups, as intended, half are dense structures embedded within larger groups, and one-third are chance alignments within larger groups. The low completeness of the extracted CG sample is mainly due to the flux limits of the selection criteria. Our results suggest that a new observational algorithm to identify compact groups in redshift space is required to obtain dense isolated galaxy systems., Comment: 13 pages, 5 figures, 5 tables, accepted for publication in MNRAS

Sky survey telescopes and powerful targeted telescopes play complementary roles in astronomy. In order to investigate the nature and characteristics of the motions of very faint objects, a flexibly-pointed instrument capable of high astrometric accuracy is an ideal complement to current astrometric surveys and a unique tool for precision astrophysics. Such a space-based mission will push the frontier of precision astrometry from evidence of Earth-mass habitable worlds around the nearest stars, to distant Milky Way objects, and out to the Local Group of galaxies. As we enter the era of the James Webb Space Telescope and the new ground-based, adaptive-optics-enabled giant telescopes, by obtaining these high precision measurements on key objects that Gaia could not reach, a mission that focuses on high precision astrometry science can consolidate our theoretical understanding of the local Universe, enable extrapolation of physical processes to remote redshifts, and derive a much more consistent picture of cosmological evolution and the likely fate of our cosmos. Already several missions have been proposed to address the science case of faint objects in motion using high precision astrometry missions: NEAT proposed for the ESA M3 opportunity, micro-NEAT for the S1 opportunity, and Theia for the M4 and M5 opportunities. Additional new mission configurations adapted with technological innovations could be envisioned to pursue accurate measurements of these extremely small motions. The goal of this White Paper is to address the fundamental science questions that are at stake when we focus on the motions of faint sky objects and to briefly review instrumentation and mission profiles., Comment: arXiv admin note: substantial text overlap with arXiv:1910.08028, arXiv:1707.01348

We study the formation of over 6000 compact groups (CGs) of galaxies identified in mock redshift-space galaxy catalogues built from semi-analytical models of galaxy formation (SAMs) run on the Millennium Simulations. We select CGs of 4 members in our mock SDSS galaxy catalogues and, for each CG, we trace back in time the real-space positions of the most massive progenitors of their 4 galaxies. By analysing the evolution of the distance of the galaxy members to the centre of mass of the group, we identify 4 channels of CG formation. The classification of these assembly channels is performed with an automatic recipe inferred from a preliminary visual inspection and based on the orbit of the galaxy with the fewest number of orbits. Most CGs show late assembly, with the last galaxy arriving on its first or second passage, while only 10-20 per cent form by the gradual contraction of their orbits by dynamical friction, and only a few per cent forming early with little subsequent contraction. However, a SAM from a higher resolution simulation leads to earlier assembly. Assembly histories of CGs also depend on cosmological parameters. At similar resolution, CGs assemble later in SAMs built on parent cosmological simulations of high density parameter. Several observed properties of mock CGs correlate with their assembly history: early-assembling CGs are smaller, with shorter crossing times, and greater magnitude gaps between their brightest two members, and their brightest galaxies have smaller spatial offsets and are more passive., Comment: 14 pages, 11 figures, 3 tables, accepted for publication in MNRAS

We use simulations to study the growth of a pseudobulge in an isolated thin exponential stellar disc embedded in a static spherical halo. We observe a transition from later to earlier morphological types and an increase in bar prominence for higher disc-to-halo mass ratios, for lower disc-to-halo size ratios, and for lower halo concentrations. We compute bulge-to-total stellar mass ratios $B/T$ by fitting a two-component S\'ersic-exponential surface-density distribution. The final $B/T$ is strongly related to the disc's fractional contribution $f_{\rm d}$ to the total gravitational acceleration at the optical radius. The formula $B/T=0.5\,f_{\rm }^{1.8}$ fits the simulations to an accuracy of $30\%$, is consistent with observational measurements of B/T and f_d as a function of luminosity, and reproduces the observed relation between $B/T$ and stellar mass when incorporated into the GalICS~2.0 semi-analytic model of galaxy formation., Comment: 19 pages, 11 figures, accepted for publication in Astronomy & Astrophysics

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., Comment: 27 pages, 20 figures, accepted by Astronomy & Astrophysics

Mergers can be detected as double-peak narrow emission line galaxies but they are difficult to disentangle from disc rotations and gas outflows. We aim to properly detect such galaxies and distinguish the underlying mechanisms. Relying on RCSED, we developed an automated selection procedure and found 5663 double-peak emission line galaxies at z<0.34 corresponding to 0.8% of the parent database. To characterise these galaxies, we built a single-peak no-bias control sample (NBCS) with the same redshift and stellar mass distributions as the double-peak sample (DPS). These two samples are indeed very similar in terms of absolute magnitude, [OIII] luminosity, colour-colour diagrams, age and specific star formation rate, metallicity, and environment. We find an important excess of S0 galaxies in the DPS, not observed in the NBCS, and which cannot be accounted for by the environment, as most of these galaxies are isolated or in poor groups. Similarly, we find a relative deficit of pure discs in the DPS late-type galaxies, that are preferentially of Sa type. In parallel, we observe a systematic central excess of star formation and extinction for DP galaxies. Finally, there are noticeable differences in the kinematics: the gas velocity dispersion is correlated with the galaxy inclination in the NBCS, whereas this relation does not hold for the DPS. Furthermore, the DP galaxies show larger stellar velocity dispersions and they deviate from the Tully-Fisher relation for both late-type and S0 galaxies. These discrepancies can be reconciled if one considers the two peaks as two different components. Considering the morphological biases in favour, bulge-dominated galaxies and star-formation central enhancement, we suggest a scenario of multiple sequential minor mergers driving the increase of the bulge size, leading to larger fractions of S0 galaxies and a deficit of pure disc galaxies., Comment: 30 pages, 27 figures, article accepted to A&A

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)., Comment: 21 pages. Comments welcome. This version 4 has a different title in the PDF to match that of arXiv, and a retouch of the last sentence of the last section and of the Acknowledgements

The S\'ersic 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\'ersic 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\'ersic index ($n \lesssim 1.5$). We have derived a more precise analytical approximation to the deprojected S\'ersic 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\'ersic 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\'ersic indices, for both density and mass profiles. An appendix compares the deprojected S\'ersic profiles with those of the popular simple models from Plummer (1911), Jaffe (1983), Hernquist (1990), Navarro et al. (1996), and Einasto (1965)., Comment: 8 pages, 6 figures, 1 table, accepted by Astronomy & Astrophysics

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., Comment: Published version, consolidated with Erratum on caption of Table A1 (the table is for all 404 931 galaxies in the Clean sample, not just those that are candidate VYGs according to at least one spectral model). New is in green and omitted is in strikeout-red

Compact groups (CGs) of galaxies are defined as isolated and dense galaxy systems that appear to be a unique site of multiple galaxy interactions. Semi-analytical models of galaxy formation (SAMs) are a prime tool to understand CGs. We investigate how the frequency and the three-dimensional nature of CGs depends on the SAM and its underlying cosmological parameters. Extracting 9 lightcones of galaxies from 5 different SAMs and selecting CGs as in observed samples, we find that the frequency and nature of CGs depends strongly on the cosmological parameters. Moving from the WMAP1 to the WMAP7 and Planck cosmologies (increasing density of the Universe and decreasing normalisation of the power spectrum), the space density of CGs is decreased by a factor 2.5, while the fraction of CGs that are physically dense falls from 50 to 35 percent. The lower $\sigma_8$ leads to fewer dense groups, while the higher $\Omega_{\rm m}$ causes more chance alignments. However, with increased mass and spatial resolution, the fraction of CGs that are physically dense is pushed back up to 50 percent. The intrinsic differences in the SAM recipes also lead to differences in the frequency and nature of CGs, particularly those related to how SAMs treat orphan galaxies. We find no dependence of CG properties on the flux limit of the mock catalogues nor on the waveband in which galaxies are selected. One should thus be cautious when interpreting a particular SAM for the frequency and nature of CGs., Comment: 17 pages, 10 figures, 4 tables, accepted for publication in MNRAS

We describe a physical model of the outflows produced as a result of gas accretion onto a black hole, and the resultant changes to star formation rates and efficiencies in galaxies, using the Radio-SAGE semi-analytic galaxy formation model. We show that the ratio of outflow rate to SFR of galaxies is mainly driven by black hole mass and virial halo mass, and show that the SFR is higher than the outflow rate at low black hole masses. The model consistently reproduces the observed evolution of star formation rate density from z = 6 to z = 0, as well as the trend of the stellar mass - halo mass relations. We show the characteristic growth of massive galaxies influenced by AGN feedback at different redshifts. We find feedback to be prevalent in the most massive galaxy halos, inhibiting the cooling catastrophe., Comment: 15 pages, 8 figures. Accepted for publication in MNRAS

We present a comparison of the observed evolving galaxy stellar mass functions with the predictions of eight semi-analytic models and one halo occupation distribution model. While most models are able to fit the data at low redshift, some of them struggle to simultaneously fit observations at high redshift. We separate the galaxies into 'passive' and 'star-forming' classes and find that several of the models produce too many low-mass star-forming galaxies at high redshift compared to observations, in some cases by nearly a factor of 10 in the redshift range $2.5 < z < 3.0$. We also find important differences in the implied mass of the dark matter haloes the galaxies inhabit, by comparing with halo masses inferred from observations. Galaxies at high redshift in the models are in lower mass haloes than suggested by observations, and the star formation efficiency in low-mass haloes is higher than observed. We conclude that many of the models require a physical prescription that acts to dissociate the growth of low-mass galaxies from the growth of their dark matter haloes at high redshift., Comment: 16 pages, 8 figures, accepted for publication in MNRAS

We investigate the dynamical evolution of galaxies in groups with different formation epochs. Galaxy groups have been selected to be in different dynamical states, namely dynamically old and dynamically young, which reflect their early and late formation times, respectively, based on their halo mass assembly. Brightest galaxies in dynamically young groups have suffered their last major galaxy merger typically $\sim 2$ Gyr more recently than their counterparts in dynamically old groups. Furthermore, we study the evolution of velocity dispersion in these two classes and compare them with the analytic models of isolated halos. The velocity dispersion of dwarf galaxies in high mass, dynamically young groups increases slowly in time, while the analogous dispersion in dynamically old high-mass groups is constant. In contrast, the velocity dispersion of giant galaxies in low mass groups decreases rapidly at late times. This increasing velocity bias is caused by dynamical friction, and starts much earlier in the dynamically old groups. The recent {\sc Radio-SAGE} model of galaxy formation suggests that radio luminosities of central galaxies, considered to be tracers of AGN activity, are enhanced in halos that assembled more recently, independent of the time since the last major merger., Comment: 11 pages, 7 figures. Accepted for publication in ApJ

In the local Universe, the existence of very young galaxies (VYGs), having formed at least half their stellar mass in the last 1 Gyr, is debated. We predict the present-day fraction of VYGs among central galaxies as a function of galaxy stellar mass. For this, we apply to high mass resolution Monte-Carlo halo merger trees (MCHMTs) three (one) analytical models of galaxy formation, where the ratio of stellar to halo mass (mass growth rate) is a function of halo mass and redshift. Galaxy merging is delayed until orbital decay by dynamical friction. With starbursts associated with halo mergers, our models predict typically one percent of VYGs up to galaxy masses of $10^{10}$ M$_\odot$, falling rapidly at higher masses, and VYGs are usually associated with recent major mergers of their haloes. Without these starbursts, two of the models have VYG fractions reduced by 1 or 2 dex at low or intermediate stellar masses, and VYGs are rarely associated with major halo mergers. In comparison, the state-of-the-art semi-analytical model (SAM) of Henriques et al. produces only 0.01% of VYGs at intermediate masses. Finally, the Menci et al. SAM run on MCMHTs with Warm Dark Matter cosmology generates 10 times more VYGs at masses below $10^8$ M$_\odot$ than when run with Cold Dark Matter. The wide range in these VYG fractions illustrates the usefulness of VYGs to constrain both galaxy formation and cosmological models., Comment: Accepted in MNRAS. 26 pages, 20 figures

We present a comparison of nine galaxy formation models, eight semi-analytical and one halo occupation distribution model, run on the same underlying cold dark matter simulation (cosmological box of co-moving width 125$h^{-1}$ Mpc, with a dark-matter particle mass of $1.24\times 10^9 h^{-1}$ Msun) and the same merger trees. While their free parameters have been calibrated to the same observational data sets using two approaches, they nevertheless retain some 'memory' of any previous calibration that served as the starting point (especially for the manually-tuned models). For the first calibration, models reproduce the observed z = 0 galaxy stellar mass function (SMF) within 3-{\sigma}. The second calibration extended the observational data to include the z = 2 SMF alongside the z~0 star formation rate function, cold gas mass and the black hole-bulge mass relation. Encapsulating the observed evolution of the SMF from z = 2 to z = 0 is found to be very hard within the context of the physics currently included in the models. We finally use our calibrated models to study the evolution of the stellar-to-halo mass (SHM) ratio. For all models we find that the peak value of the SHM relation decreases with redshift. However, the trends seen for the evolution of the peak position as well as the mean scatter in the SHM relation are rather weak and strongly model dependent. Both the calibration data sets and model results are publicly available., Comment: 20 pages, 14 figures (+5 supplementary figures in the Appendix), accepted for publication in MNRAS

We carry out a systematic investigation of the total mass density profile of massive (Mstar~3e11 Msun) early-type galaxies and its dependence on redshift, specifically in the range 0

Dwarf galaxy anomalies, such as their abundance and cusp-core problems, remain a prime challenge in our understanding of galaxy formation. The inclusion of baryonic physics could potentially solve these issues, but the efficiency of stellar feedback is still controversial. We analytically explore the possibility of feedback from Active Galactic Nuclei (AGN) in dwarf galaxies and compare AGN and supernova (SN) feedback. We assume the presence of an intermediate mass black hole within low mass galaxies and standard scaling relations between the relevant physical quantities. We model the propagation and properties of the outflow and explore the critical condition for global gas ejection. Performing the same calculation for SNe, we compare the ability of AGN and SNe to drive gas out of galaxies. We find that a critical halo mass exists below which AGN feedback can remove gas from the host halo and that the critical halo mass for AGN is greater than the equivalent for SNe in a significant part of the parameter space, suggesting that AGN could provide an alternative and more successful source of negative feedback than SNe, even in the most massive dwarf galaxies., Comment: 6 pages, 5 figures

[Abridged] We use the WINGS database to select a sample of 67 nearby galaxy clusters with at least 30 spectroscopic members each. 53 of these clusters do not show evidence of substructures in phase-space, while 14 do. We estimate the virial radii and circular velocities of the 67 clusters by a variety of proxies (velocity dispersion, X-ray temperature, and richness) and use these estimates to build stack samples from these 53 and 14 clusters ('Reg' and 'Irr' stacks, respectively). We determine the number-density and velocity-dispersion profiles (VDPs) of E, S0, and Sp+Irr (S) galaxies in the Reg and Irr samples, separately, and fit models to these profiles. The number density profiles of E, S0, and S galaxies are adequately described by either a NFW or a cored King model, both for the Reg and Irr samples, with a slight preference for the NFW model. The spatial distribution concentration increases from the S to the S0 and to the E populations, both in the Reg and the Irr stacks, reflecting the well-known morphology-radius relation. Reg clusters have a more concentrated spatial distribution of E and S0 galaxies than Irr clusters, while the spatial distributions of S galaxies in Reg and Irr clusters are similar. We propose a new phenomenological model that provides acceptable fits to the VDP of all our galaxy samples. The VDPs become steeper and with a higher normalization from E to S0 to S galaxies. The S0 VDP is close to that of E galaxies in Reg clusters, and intermediate between those of E and S galaxies in Irr clusters. Our results suggest that S galaxies are a recently accreted cluster population, that take less than 3 Gyr to evolve into S0 galaxies after accretion, and in doing so modify their phase-space distribution, approaching that of cluster ellipticals. While in Reg clusters this evolutionary process is mostly completed, it is still ongoing in Irr clusters., Comment: A&A, in press - 11 pages, 9 figures, 4 tables

We estimate the stellar mass that satellite galaxies lose once they enter groups (and clusters) by identifying groups in a high-resolution cosmological N-body simulation, assigning entry masses to satellite galaxies with halo abundance matching at the entry time, and comparing the predicted conditional stellar mass function of satellite galaxies at $z\simeq0$ with observations. Our results depend on the mass of the stars that form in satellite galaxies after the entry time. A model in which star formation shuts down completely as soon a galaxy enters a group environment is ruled out because it underpredicts the stellar masses of satellite galaxies even in the absence of tidal stripping. The greater is the stellar mass that is allowed to form, the greater the fraction that needs to be tidally stripped. The stellar mass fraction lost by satellite galaxies after entering a group or cluster environment is consistent with any value in the range $0-40\%$. To place stronger constraints, we consider a more refined model of tidal stripping of galaxies on elongated orbits (where stripping occurs at orbit pericentres). Our model predicts less tidal stripping: satellite galaxies lose $\sim 20-25\%$ of their stellar mass since their entry into the group. This finding is consistent with a slow-starvation delayed-quenching picture, in which galaxies that enter a group or cluster environment keep forming stars until at least the first pericentric passage., Comment: Accepted for publication in MNRAS. 27 pages, 22 figures

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision., Comment: http://theia.phyip3.dur.ac.uk

While the density profiles (DPs) of $\Lambda$CDM haloes obey the NFW law out to roughly one virial radius, $r_{\rm vir}$, the structure of their outer parts is still poorly understood, since the 1-halo term describing the halo itself is dominated by the 2-halo term representing the other haloes picked up. Using a semi-analytical model, we measure the real-space `1-halo' number DP of groups out to $20\,r_{\rm vir}$ by assigning each galaxy to its nearest group with mass above $M_{\rm a}$, in units of the group $r_{\rm vir}$. If $M_{\rm a}$ is small (large), the outer DP of groups falls rapidly (slowly). We find that there is an optimal $M_{\rm a}$ for which the stacked DP resembles the NFW model to $0.1$ dex accuracy out to $\simeq 13\,r_{\rm vir}$. We find similar long-range NFW surface DPs (out to $\simeq 10\,r_{\rm vir}$) in the SDSS observations using a galaxy assignment scheme that combines the non-linear virialized regions of groups with their linear outer parts. The optimal $M_{\rm a}$ scales as the minimum mass of the groups that are stacked to the power $0.25-0.3$. Our results suggest that the NFW model does not solely originate from violent relaxation. Moreover, populating haloes with galaxies using HOD models must proceed out to larger radii than usually done., Comment: 5 pages, 4 figures. Accepted for publication in MNRAS Letters. Final version including discussion on the backsplash radius

The gap between first and second ranked galaxy magnitudes in groups is often considered a tracer of their merger histories, which in turn may affect galaxy properties, and also serves to test galaxy luminosity functions (LFs). We remeasure the conditional luminosity function (CLF) of the Main Galaxy Sample of the SDSS in an appropriately cleaned subsample of groups from the Yang catalog. We find that, at low group masses, our best-fit CLF have steeper satellite high ends, yet higher ratios of characteristic satellite to central luminosities in comparison with the CLF of Yang et al. (2008). The observed fractions of groups with large and small magnitude gaps as well as the Tremaine & Richstone (1977) statistics, are not compatible with either a single Schechter LF or with a Schechter-like satellite plus lognormal central LF. These gap statistics, which naturally depend on the size of the subsamples, and also on the maximum projected radius, $R_{\rm max}$, for defining the 2nd brightest galaxy, can only be reproduced with two-component CLFs if we allow small gap groups to preferentially have two central galaxies, as expected when groups merge. Finally, we find that the trend of higher gap for higher group velocity dispersion, $\sigma_{\rm v}$, at given richness, discovered by Hearin et al. (2013), is strongly reduced when we consider $\sigma_{\rm v}$ in bins of richness, and virtually disappears when we use group mass instead of $\sigma_{\rm v}$. This limits the applicability of gaps in refining cosmographic studies based on cluster counts., Comment: 18 pages, 13 figures, 5 tables, submitted to MNRAS, comments are welcome

We compare the mean mass assembly histories of compact and fossil galaxy groups in the Millennium dark matter simulation and an associated semi-analytic galaxy formation model. Tracing the halo mass of compact groups (CGs) from z=0 to z=1 shows that, on average, 55 per cent of the halo mass in compact groups is assembled since z~1, compared to 40 per cent of the halo mass in fossil groups (FGs) in the same time interval, indicating that compared to FGs, CGs are relatively younger galaxy systems. At z=0, for a given halo mass, fossil groups tend to have a larger concentration than compact groups. Investigating the evolution of CG's parameters show that they become more compact with time. CGs at z=0.5 see their magnitude gaps increase exponentially, but it takes ~10 Gyr for them to reach a magnitude gap of 2 magnitudes. The slow growth of the magnitude gap leads to only a minority (~41 per cent) of CGs selected at z=0.5 turning into a FG by z=0. Also, while three-quarters of FGs go through a compact phase, most fail to meet the CG isolation criterion, leaving only ~30 per cent of FGs fully satisfying the CG selection criteria. Therefore, there is no strong link of CGs turning into FGs or FGs originating from CGs. The relation between CGs and FGs is thus more complex, and in most cases, FGs and CGs follow different evolutionary tracks., Comment: 11 pages, 11 figures, Accepted for publication in ApJ

We carry out a systematic investigation of the total mass density profile of massive (Mstar>2e11 Msun) early-type galaxies and its dependence on galactic properties and host halo mass with the aid of a variety of lensing/dynamical data and large mock galaxy catalogs. The latter are produced via semi-empirical models that, by design, are based on just a few basic input assumptions. Galaxies, with measured stellar masses, effective radii and S\'{e}rsic indices, are assigned, via abundance matching relations, host dark matter halos characterized by a typical LCDM profile. Our main results are as follows: (i) In line with observational evidence, our semi-empirical models naturally predict that the total, mass-weighted density slope at the effective radius gamma' is not universal, steepening for more compact and/or massive galaxies, but flattening with increasing host halo mass. (ii) Models characterized by a Salpeter or variable initial mass function and uncontracted dark matter profiles are in good agreement with the data, while a Chabrier initial mass function and/or adiabatic contractions/expansions of the dark matter halos are highly disfavored. (iii) Currently available data on the mass density profiles of very massive galaxies (Mstar>1e12 Msun), with Mhalo>3e14 Msun, favor instead models with a stellar profile flatter than a S\'{e}rsic one in the very inner regions (r<3-5 kpc), and a cored NFW or Einasto dark matter profile with median halo concentration a factor of ~2 or <1.3, respectively, higher than those typically predicted by N-body numerical simulations., Comment: 23 pages, 10 figures, 3 Appendices (with an extra 7 figures). ApJ, accepted. Main results in Figures 3, 5, 6, 8

We present a clustering comparison of 12 galaxy formation models (including Semi-Analytic Models (SAMs) and Halo Occupation Distribution (HOD) models) all run on halo catalogues and merger trees extracted from a single {\Lambda}CDM N-body simulation. We compare the results of the measurements of the mean halo occupation numbers, the radial distribution of galaxies in haloes and the 2-Point Correlation Functions (2PCF). We also study the implications of the different treatments of orphan (galaxies not assigned to any dark matter subhalo) and non-orphan galaxies in these measurements. Our main result is that the galaxy formation models generally agree in their clustering predictions but they disagree significantly between HOD and SAMs for the orphan satellites. Although there is a very good agreement between the models on the 2PCF of central galaxies, the scatter between the models when orphan satellites are included can be larger than a factor of 2 for scales smaller than 1 Mpc/h. We also show that galaxy formation models that do not include orphan satellite galaxies have a significantly lower 2PCF on small scales, consistent with previous studies. Finally, we show that the 2PCF of orphan satellites is remarkably different between SAMs and HOD models. Orphan satellites in SAMs present a higher clustering than in HOD models because they tend to occupy more massive haloes. We conclude that orphan satellites have an important role on galaxy clustering and they are the main cause of the differences in the clustering between HOD models and SAMs., Comment: 17 pages, 7 figures

We report the discovery of a very diverse set of five low-surface brightness (LSB) dwarf galaxy candidates in Hickson Compact Group 90 (HCG 90) detected in deep U- and I-band images obtained with VLT/VIMOS. These are the first LSB dwarf galaxy candidates found in a compact group of galaxies. We measure spheroid half-light radii in the range $0.7\!\lesssim\! r_{\rm eff}/{\rm kpc}\! \lesssim\! 1.5$ with luminosities of $-11.65\!\lesssim\! M_U\! \lesssim\! -9.42$ and $-12.79\!\lesssim\! M_I\! \lesssim\! -10.58$ mag, corresponding to a color range of $(U\!-\!I)_0\!\simeq\!1.1\!-\!2.2$ mag and surface brightness levels of $\mu_U\!\simeq\!28.1\,{\rm mag/arcsec^2}$ and $\mu_I\!\simeq\!27.4\,{\rm mag/arcsec^2}$. Their colours and luminosities are consistent with a diverse set of stellar population properties. Assuming solar and 0.02 Z$_\odot$ metallicities we obtain stellar masses in the range $M_*|_{Z_\odot} \simeq 10^{5.7-6.3} M_{\odot}$ and $M_*|_{0.02\,Z_\odot}\!\simeq\!10^{6.3-8}\,M_{\odot}$. Three dwarfs are older than 1 Gyr, while the other two significantly bluer dwarfs are younger than $\sim 2$ Gyr at any mass/metallicity combination. Altogether, the new LSB dwarf galaxy candidates share properties with dwarf galaxies found throughout the Local Volume and in nearby galaxy clusters such as Fornax. We find a pair of candidates with $\sim\!2$ kpc projected separation, which may represent one of the closest dwarf galaxy pairs found. We also find a nucleated dwarf candidate, with a nucleus size of $r_{\rm eff}\!\simeq\!46\!-\!63$ pc and magnitude M$_{U,0}=-7.42$ mag and $(U\!-\!I)_0\!=\!1.51$ mag, which is consistent with a nuclear stellar disc with a stellar mass in the range $10^{4.9-6.5}\,M_\odot$., Comment: Accepted by MNRAS. 7 pages, 5 figures, 2 tables

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., Comment: submitted to MNRAS

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 $\gamma$-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., Comment: 22 pages + references, 4 figures. v2: a few comments and references added, version published in JCAP

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., Comment: 11 pages, 10 figures, MNRAS accepted

Feedback from active galactic nuclei (AGN) has often been invoked both in simulations and in interpreting observations for regulating star formation and quenching cooling flows in massive galaxies. AGN activity can, however, also over-pressurise the dense star-forming regions of galaxies and thus enhance star formation, leading to a positive feedback effect. To understand this pressurisation better, we investigate the effect of an ambient external pressure on gas fragmentation and triggering of starburst activity by means of hydrodynamical simulations. We find that moderate levels of over-pressurisation of the galaxy boost the global star formation rate of the galaxy by an order of magnitude, turn stable discs unstable, and lead to significant fragmentation of the gas content of the galaxy, similar to what is observed in high redshift galaxies., Comment: submitted to MNRAS

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., Comment: 35 pages, 20 figures, accepted for publication in MNRAS

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., Comment: 5 pages, 4 figures. Submitted to ApJ

Combining our knowledge of halo structure and internal kinematics from cosmological dark matter simulations and the distribution of halo interlopers in projected phase space measured in cosmological galaxy simulations, we develop MAGGIE, a prior- and halo-based, probabilistic, abundance matching (AM) grouping algorithm for doubly complete subsamples (in distance and luminosity) of flux-limited samples. We test MAGGIE-L and MAGGIE-M (in which group masses are derived from AM applied to the group luminosities and stellar masses, respectively) on groups of at least three galaxies extracted from a mock Sloan Digital Sky Survey Legacy redshift survey, incorporating realistic observational errors on galaxy luminosities and stellar masses. In comparison with the optimal Friends-of-Friends (FoF) group finder, groups extracted with MAGGIE are much less likely to be secondary fragments of true groups; in primary fragments, its galaxy memberships (relative to the virial sphere of the real-space group) are much more complete and usually more reliable, and its masses are much less biased and usually with less scatter, as are its group luminosities and stellar masses (computed in MAGGIE using the membership probabilities as weights). FoF outperforms MAGGIE only for high-mass clusters: for the reliability of the galaxy population and the dispersion of its total mass. In comparison with our implementation of the Yang et al. group finder, MAGGIE reaches much higher completeness and slightly lower group fragmentation and dispersion on group total masses, luminosities and stellar masses, but slightly greater bias in the latter two and lower reliabilities. MAGGIE should therefore lead to sharper trends of environmental effects on galaxies and more accurate mass/orbit modelling., Comment: Version published in MNRAS consolidated with Erratum and additional fixes

[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., Comment: 18 pages, 12 figures, accepted for publication in MNRAS

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., Comment: 6 pages, 3 figures, 1 table, MNRAS (Letters) in press

We test the Friends-of-Friends (FoF) grouping algorithm, which depends on two linking lengths (LLs), plane-of-sky and line-of-sight (LOS), normalized to the mean nearest neighbor separation of field galaxies. We argue, on theoretical grounds, that LLs should be $b_\perp \simeq 0.11$, and $b_\parallel \simeq 1.3$ to recover 95% of all galaxies with projected radii within the virial radius $r_{200}$ and 95% of the galaxies along the LOS. We then predict that 80 to 90% of the galaxies in FoF groups should lie within their parent real-space groups (RSGs), defined within their virial spheres. We test the FoF extraction for 16x16 pairs of LLs, using subsamples of galaxies, doubly complete in distance and luminosity, of a flux-limited mock SDSS galaxy catalog. We find that massive RSGs are more prone to fragmentation, while the fragments typically have low estimated mass, with typically 30% of groups of low and intermediate estimated mass being fragments. Group merging rises drastically with estimated mass. For groups of 3 or more galaxies, galaxy completeness and reliability are both typically better than 80% (after discarding the fragments). Estimated masses of extracted groups are biased low, by up to a factor 4 at low richness, while the inefficiency of mass estimation improves from 0.85 dex to 0.2 dex when moving from low to high multiplicity groups. The optimal LLs depend on the scientific goal for the group catalog. We propose $b_\perp \simeq 0.07$, with $b_\parallel \simeq 1.1$ for studies of environmental effects, $b_\parallel \simeq 2.5$ for cosmographic studies and $b_\parallel \simeq 5$ for followups of individual groups., Comment: Revised version resubmitted to MNRAS. 16 pages. Smaller optimal linking lengths. Two extra figures. Comments welcome!