Search

Relic galaxies, the oldest ultra-compact massive galaxies (UCMGs), contain almost exclusively "pristine" stars formed during an intense star formation (SF) burst at high redshift. As such, they allow us to study in detail the early mechanism of galaxy assembly in the Universe. Using the largest catalogue of spectroscopically confirmed UCMGs for which a degree of relicness (DoR) had been estimated, the INSPIRE catalogue, we investigate whether or not relics prefer dense environments. The objective of this study is to determine if the DoR, which measures how extreme the SF history was, and the surrounding environment are correlated. In order to achieve this goal, we employ the AMICO galaxy cluster catalogue to compute the probability for a galaxy to be a member of a cluster, and measure the local density around each UCMG using machine learning-based photometric redshifts. We find that UCMGs can reside both in clusters and in the field, but objects with very low DoR (< 0.3, i.e., a relatively extended SF history) prefer under-dense environments. We additionally report a correlation between the DoR and the distance from the cluster centre: more extreme relics, when located in clusters, tend to occupy the more central regions of them. We finally outline potential evolution scenarios for UCMGs at different DoR to reconcile their presence in both clusters and field environments, Comment: Accepted for publication on MNRAS, 11 pages, 8 figure, 1 table

Within the $\Lambda$CDM cosmology, dark matter haloes are expected to deviate from spherical symmetry. Constraining the halo shapes at large galactocentric distances is challenging due to the low density of luminous tracers. The well-studied early-type galaxy NGC 5128 (Centaurus A - CenA), has a large number of radial velocities for globular clusters (GCs) and planetary nebulae (PNe) of its extended stellar halo. In this work, we aim to determine the deviation from spherical symmetry of the dark matter halo of CenA at 5 $R_{\rm e}$ using its GCs as kinematic tracers. We used the largest photometric catalogue of GC candidates to accurately characterise the spatial distribution of the relaxed population and investigated the presence of non-relaxed structures in the kinematic catalogue of GCs using the relaxed point-symmetric velocity field as determined by the host's PNe population. We used anisotropic Jeans modelling under axisymmetric assumptions together with the Gaussian likelihood and GCs as discrete tracers. The gravitational potential is generated by flattened stellar and dark matter distributions. We leveraged different orbital properties of the blue and red GCs to model them separately. We find that discrete kinematics of the GCs are consistent with being drawn from an underlying relaxed velocity field determined from PNe. The best-fit parameters of the gravitational potential recovered from the blue and red GCs separately agree well and the joint results are: $M_{200} = 1.86^{1.61}_{-0.69}\times 10^{12}$ M$_\odot$, $M_\star/L_{\rm B} = 2.98^{+0.96}_{-0.78}$ and the flattening $q_{\rm DM} = 1.45^{+0.78}_{-0.53}$. Both GC populations show mild rotation, with red having a slightly stronger rotational signature and radially biased orbits, and blue GCs preferring negative velocity anisotropy. An oblate or a spherical dark matter halo of CenA is strongly disfavoured by our modelling., Comment: 22 pages, 17 figures, 8 tables, accepted for publication in A&A

Planetary nebulae (PNe) are excellent tracers of the metal-poor haloes of nearby early-type galaxies. They are commonly used to trace spatial distribution and kinematics of the halo and intracluster light at distances of up to 100 Mpcs. The results on the early-type galaxy M105 in the Leo I group represent a benchmark for the quantitative analysis of halo and intragroup light. Since the Leo I group lies at just a 10 Mpc distance, it is at the ideal location to compare results from resolved stellar populations with the homogeneous constraints over a much larger field of view from the PN populations. In M105, we have -- for the first time -- established a direct link between the presence of a metal-poor halo as traced by resolved red-giant branch stars and a PN population with a high specific frequency ($\alpha$-parameter). This confirms our inferences that the high $\alpha$-parameter PN population in the outer halo of M49 in the Virgo Cluster traces the metal-poor halo and intra-group light., Comment: 12 pages, 5 figures, accepted for publication in the proceedings of IAU Symposium 384: Planetary Nebulae: a Universal Toolbox in the Era of Precision Astrophysics

Planetary nebulae (PNe) are essential tracers of the kinematics of the diffuse halo and intracluster light where stellar spectroscopy is unfeasible, due to their strong emission lines. However, that is not all they can reveal about the underlying stellar population. In recent years, it has also been found that PNe in the metal-poor halos of galaxies have different properties (specific frequency, luminosity function), than PNe in the more metal-rich galaxy centers. A more quantitative understanding of the role of age and metallicity in these relations would turn PNe into valuable stellar-population tracers. In order to do that, a full characterization of PNe in regions where the stellar light can also be analysed in detail is necessary. In this work, we make use of integral-field spectroscopic data covering the central regions of galaxies, which allow us to measure both stellar ages and metallicities as well as to detect PNe. This analysis is fundamental to calibrate PNe as stellar population tracers and to push our understanding of galaxy properties at unprecedented galactocentric distances., Comment: Submitted to the Proceedings of IAU Symposium 384: Planetary Nebulae: a Universal Toolbox in the Era of Precision Astrophysics. Eds: O. De Marco, A. Zijlstra, R. Szczerba

This paper presents the third data release of the INvestigating Stellar Population In RElics (INSPIRE) project, comprising 52 ultra-compact massive galaxies (UCMGs) observed with the X-Shooter spectrograph. We measure integrated stellar velocity dispersion, [Mg/Fe] abundances, ages, and metallicities for all the INSPIRE objects. We thus infer star formation histories and confirm the existence of a degree of relicness (DoR), defined in terms of the fraction of stellar mass formed by $z=2$, the time at which a galaxy has assembled 75\% of its mass, and the final assembly time. Objects with a high DoR assembled their stellar mass at early epochs, while low-DoR objects show a non-negligible fraction of later-formed populations and hence a spread in ages and metallicities. A higher DoR correlates with larger [Mg/Fe], super-solar metallicity, and larger velocity dispersion values. The 52 UMCGs span a large range of DoR from 0.83 to 0.06, with 38 of them having formed more than 75\% of their mass by $z=2$. Of these, nine are extreme relics (DoR$>0.7$), since they formed the totality ($>99\%$) of their stellar mass by redshift $z=2$. The remaining 14 UCMGs cannot be considered relics, as they are characterised by more extended star formation histories. With INSPIRE, we built the first sizeable sample of relics outside the local Universe, up to $z\sim0.4$, increasing the number of confirmed relics by a factor of $>10$, and opening up an important window to explain the mass assembly of massive galaxies in the high-z Universe., Comment: accepted for publication on MNRAS, 20 pages, 15 figures, 3 tables

Cosmological simulations predict dark matter shapes that deviate from spherical symmetry. The exact shape depends on the prescription of the simulation and the interplay between dark matter and baryons. This signature is most pronounced in the diffuse galactic haloes that can be observationally probed with planetary nebulae and globular clusters (GCs). The kinematic observations of these halo tracers support intrinsic triaxial shape for the mass generating the gravitational potential. With discrete axisymmetric modelling of GCs as the halo tracers of NGC 5128 we investigate the overall mass distribution of this nearby giant elliptical galaxy. Our modelling approach constrains $c_{200}$, $(M/L)_{\star, B}$ and inclination. We derive a preliminary $M_{200}\sim 1 \times 10^{12}$ M$_\odot$ and flattening $q_{\mathrm{DM}}\sim 1.3$ indicative of prolate/triaxial halo for NGC 5128., Comment: 6 pages, 2 figures, Accepted for publication in Proceedings of the IAUS 379: Dynamical Masses of Local Group Galaxies, Potsdam, March 20-24, 2023

We report high-quality H${\alpha}$/CO, imaging spectroscopy of nine massive (log median stellar mass = 10.65 $M_{\odot}$), disk galaxies on the star-forming, main sequence (henceforth `SFGs'), near the peak of cosmic galaxy evolution ($z\sim$1.1-2.5), taken with the ESO-Very Large Telescope, IRAM-NOEMA and Atacama Large Millimeter/submillimeter Array. We fit the major axis position-velocity cuts with beam-convolved, forward models with a bulge, a turbulent rotating disk, and a dark matter (DM) halo. We include priors for stellar and molecular gas masses, optical light effective radii and inclinations, and DM masses from our previous rotation curve analyses of these galaxies. We then subtract the inferred 2D model-galaxy velocity and velocity dispersion maps from those of the observed galaxies. We investigate whether the residual velocity and velocity dispersion maps show indications for radial flows. We also carry out kinemetry, a model-independent tool for detecting radial flows. We find that all nine galaxies exhibit significant non-tangential flows. In six SFGs, the inflow velocities ($v_r\sim$30-90 km s$^{-1}$, 10%-30% of the rotational component) are along the minor axis of these galaxies. In two cases the inflow appears to be off the minor axis. The magnitudes of the radial motions are in broad agreement with the expectations from analytic models of gravitationally unstable, gas-rich disks. Gravitational torques due to clump and bar formation, or spiral arms, drive gas rapidly inward and result in the formation of central disks and large bulges. If this interpretation is correct, our observations imply that gas is transported into the central regions on ~10 dynamical time scales., Comment: 31 pages, 27 figures, 3 tables. The Astrophysical Journal in press

Mass and angular momentum are key parameters of galaxies. Their co-evolution establishes an empirical relation between the specific stellar angular momentum j* and the stellar mass M* that depends on morphology. In this work, we measure j* in a sample of 32 early type galaxies (ETGs) from the ePN.S survey, using full 2D kinematic information out to a mean 6Re. We present lambda and j* profiles. We derive the distribution of these galaxies on the total j*-M* plane and determine the ratio between the stellar j* and the specific angular momentum of the host dark matter halo. The radially extended, 2D kinematic data show that the stellar halos of ETGs do not contain large stellar mass fractions of high j*. The j*-profiles of fast-rotator ETGs are largely converged within the range of the data. For slow rotators, j* is still rising and is estimated to increase beyond 6Re by up to 40%. More than 60% of their halo angular momentum is in misaligned rotation. We find that the ePN.S ETG sample displays the well-known correlation between j*, M*, and morphology: elliptical galaxies have systematically lower j* than similar mass S0 galaxies. However, fast and slow rotators lie on the same relation within errors with the slow rotators falling at the high M* end. A power-law fit to the j*-M* relation gives a slope of 0.55+-0.17 for the S0s and 0.76+-0.23 for the ellipticals, with normalisation about 4 and 9 times lower than spirals, respectively. The estimated retained fraction of angular momentum at 10^10-10^10.5 Msun is 25% for S0s and >10% for ellipticals, and decreases by 1.5 orders of magnitude at M*~10^12 Msun. Our results show that ETGs have substantially lower j* than spiral galaxies with similar M*. Their j* must be lost during their evolution, and/or retained in the hot gas component and the satellite galaxies that have not yet merged with the central galaxy. [abridged], Comment: 26 pages, 14 figures, accepted for publication in A&A

In the last decade, growing evidence has emerged supporting a non-universal stellar Initial Mass Function (IMF) in massive galaxies, with a larger number of dwarf stars with respect to the Milky-Way (bottom-heavy IMF). However, a consensus about the mechanisms that cause IMF variations is yet to be reached. Recently, it has been suggested that stars formed early-on in cosmic time, via a star formation burst, could be characterised by a bottom-heavy IMF. A promising way to confirm this is to use relics, ultra-compact massive galaxies, almost entirely composed by these "pristine" stars. The INSPIRE Project aims at assembling a large sample of confirmed relics, that can serve as laboratory to investigate on the conditions of star formation in the first 1-3 Gyr of the Universe. In this third INSPIRE paper, we build a high signal-to-noise spectrum from five relics and one from five galaxies with similar sizes, masses, and kinematical properties, but characterised by a more extended star formation history (non-relics). Our detailed stellar population analysis suggests a systematically bottom-heavier IMF slope for relics than for non-relics, adding new observational evidence for the non-universality of the IMF at various redshifts and further supporting the above proposed physical scenario., Comment: Accepted for publication in MNRAS. 8 pages, 5 figures

This is the second data release (DR2) of the INvestigating Stellar Population In RElics (INSPIRE) project, comprising 21 new systems with observations completed before March 2022. For each system, we release four one-dimensional (1D) spectra to the ESO Science Archive, one spectrum for each arm of the X-Shooter spectrograph. In this paper, we focus on the line-of-sight velocity distribution, measuring integrated stellar velocity dispersions from the spectra, and assessing their robustness and the associated uncertainties. For each of the 21 new systems, we systematically investigated the effect of the parameters and set-ups of the full spectral fitting on the stellar velocity dispersion ($\sigma$) measurements. In particular, we tested how $\sigma$ changes when several parameters of the fit as well as the resolution and spectral coverage of the input spectra are varied. We found that the effect that causes the largest systematic uncertainties on $\sigma$ is the wavelength range used for the fit, especially for spectra with a lower signal-to-noise ratio (S/N $\leq$ 30). When using blue wavelengths (UVB arm) one generally underestimates the velocity dispersion (by $\sim$15 km/s). The values obtained from the near-IR (NIR) arm present a larger scatter because the quality of the spectra is lower. We finally compared our results with those in literature, finding a very good agreement overall. Joining results obtained in DR1 with those presented here, INSPIRE contains 40 ultra-compact massive galaxies, corresponding to 75% of the whole survey. By plotting these systems in a stellar mass-velocity dispersion diagram, we identify at least four highly reliable relic candidates among the new systems. Their velocity dispersion is larger than that of normal-sized galaxies of similar stellar mass., Comment: Accepted for publication, and in press on A&A

We present a high-resolution kinematic study of the massive main-sequence star-forming galaxy (SFG) SDSS J090122.37+181432.3 (J0901) at z=2.259, using 0.36 arcsec ALMA CO(3-2) and 0.1-0.5 arcsec SINFONI/VLT H-alpha observations. J0901 is a rare, strongly-lensed but otherwise normal massive (log(M_star/M_sun)~11) main sequence SFG, offering a unique opportunity to study a typical massive SFG under the microscope of lensing. Through forward dynamical modeling incorporating lensing deflection, we fit the CO and H-alpha kinematics in the image plane out to about one disk effective radius (R_e ~ 4 kpc) at a ~600pc delensed physical resolution along the kinematic major axis. Our results show high intrinsic dispersions of the cold molecular and warm ionized gas (sig0_mol ~ 40 km/s and sig0_ion ~ 66 km/s) that remain constant out to R_e; a moderately low dark matter fraction (f_DM(R_e) ~ 0.3-0.4) within R_e; and a centrally-peaked Toomre Q-parameter -- agreeing well with the previously established sig0 vs. z, f_DM vs. Sig_baryon, and Q's radial trends using large-sample non-lensed main sequence SFGs. Our data further reveal a high stellar mass concentration within ~1-2 kpc with little molecular gas, and a clumpy molecular gas ring-like structure at R ~ 2-4 kpc, in line with the inside-out quenching scenario. Our further analysis indicates that J0901 had assembled half of its stellar mass only ~400 Myrs before its observed cosmic time, and cold gas ring and dense central stellar component are consistent with signposts of a recent wet compaction event of a highly turbulent disk found in recent simulations., Comment: 29 pages and 21 figures in total (14 pages and 10 figures in main text and the rest in appendix). Accepted for publication in ApJ

We discuss the statistical distribution of galaxy shapes and viewing angles under the assumption of triaxiality by deprojecting observed Surface Brightness (SB) profiles of 56 Brightest Cluster Galaxies coming from a recently published large deep-photometry sample. For the first time, we address this issue by directly measuring axis ratio profiles without limiting ourselves to a statistical analysis of average ellipticities. We show that these objects are strongly triaxial, with triaxiality parameters 0.39 $ \leq T \leq $ 0.72, have on average axis ratios $< p(r) > = $ 0.84 and $< q(r) > =$ 0.68, and are more spherical in the central regions but flatten out at large radii. Measured shapes in the outskirts agree well with the shapes found for simulated massive galaxies and their dark matter halos from both the IllustrisTNG and the Magneticum simulations, possibly probing the nature of dark matter. In contrast, both simulations fail to reproduce the observed inner regions of BCGs, producing too flattened objects., Comment: v2, 27 pages, 9 Figures. Published in ApJ, vol. 933, p. 215

M105 is an early-type galaxy in the nearby Leo I group, the closest galaxy group to contain all galaxy types and therefore an excellent environment to explore the low-mass end of intra-group light (IGL) assembly. We present a new extended kinematic survey of planetary nebulae (PNe) in M105 and the surrounding 30'x30' in the Leo I group with the Planetary Nebula Spectrograph. We use PNe as kinematic tracers of the diffuse stellar light in the halo and IGL and employ Gaussian mixture models to separate contributions from the companion galaxy NGC 3384 and associate PNe with halo and IGL components around M105. We present a catalogue of 314 PNe and firmly associate 93 with NGC 3384 and 169 with M105. The PNe in M105 are further associated with its halo and the surrounding exponential envelope. We construct smooth velocity and velocity dispersion fields and calculate projected rotation, velocity dispersion, and $\lambda_R$ profiles for each component. Halo PNe exhibit declining velocity dispersion and rotation profiles, while the velocity dispersion and rotation of the exponential envelope increase notably at large radii. We identify three regimes with distinct kinematics that are linked to distinct stellar population properties: (i) the rotating core (within $1~R_\mathrm{eff}$) formed in situ and dominated by metal-rich ([M/H]~0) stars likely formed in situ, (ii) the halo from 1 to $7.5~R_\mathrm{eff}$ consisting of intermediate-metallicity stars ([M/H]>-1), either formed in situ or brought in through major mergers, and (iii) the exponential envelope reaching beyond our farthest data point at 16 $R_\mathrm{eff}$, predominately composed of metal-poor ([M/H]<-1) stars. The high velocity dispersion and moderate rotation of the latter are consistent with that measured for dwarf satellite galaxies in the Leo I group, indicating that the exponential envelope traces the transition to the IGL., Comment: Submitted to A&A, comments welcome. 20 pages, 14 figures, abstract abridged

The INvestigating Stellar Population In RElics is an on-going project targeting 52 ultra-compact massive galaxies at 0.12) through a short and intense star formation burst, that evolved passively and undisturbed until the present-day. Relics provide a unique opportunity to study the mechanisms of star formation at high-z. In this paper, we present the first INSPIRE Data Release, comprising 19 systems with observations completed in 2020. We use the methods already presented in the INSPIRE Pilot, but revisiting the 1D spectral extraction. For these 19 systems, we obtain an estimate of the stellar velocity dispersion, fitting separately the two UVB and VIS XSH arms at their original resolution. We estimate [Mg/Fe] abundances via line-index strength and mass-weighted integrated stellar ages and metallicities with full spectral fitting on the combined spectrum. Ages are generally old, in agreement with the photometric ones, and metallicities are almost always super-solar, confirming the mass-metallicity relation. The [Mg/Fe] ratio is also larger than solar for the great majority of the galaxies, as expected. We find that 10 objects have formed more than 75% of their stellar mass (M*) within 3 Gyr from the Big Bang and classify them as relics. Among these, we identify 4 galaxies which had already fully assembled their M* by that time. They are therefore `extreme relics' of the ancient Universe. The INSPIRE DR1 catalogue of 10 known relics to-date augment by a factor of 3.3 the total number of confirmed relics, also enlarging the redshift window. It is therefore the largest publicly available collection. Thanks to the larger number of systems, we can also better quantify the existence of a 'degree of relicness', already hinted at the Pilot Paper., Comment: (Abstract abridged) 21 pages, 12 figures and 5 tables in the main body, plus 3 figure and 1 table in the appendix, accepted for publication on A&A. The associated data are available via the ESO Phase 3 Science Portal

Modern galaxy formation theory suggests that the misalignment between stellar and gaseous components usually results from an external gas accretion and/or interaction with other galaxies. The extreme case of the kinematic misalignment is demonstrated by so-called galaxies with counterrotation that possess two distinct components rotating in opposite directions with respect to each other. We provide an in-deep analysis of galaxies with counterrotation from IllustrisTNG100 cosmological simulations. We have found $25$ galaxies with substantial stellar counterrotation in the stellar mass range of $2\times10^{9}-3\times10^{10}$~\Msun. In our sample the stellar counterrotation is a result of an external gas infall happened $\approx 2-8$~Gyr ago. The infall leads to the initial removal of pre-existing gas, which is captured and mixed together with the infalling component. The gas mixture ends up in the counterrotating gaseous disc. We show that $\approx 90\%$ of the stellar counterrotation formed in-situ, in the counterrotating gas. During the early phases of the infall, gas can be found in inclined extended and rather thin disc-like structures, and in some galaxies they are similar to (nearly-)~polar disc or ring-like structures. We discuss a possible link between the gas infall, AGN activity and the formation of misaligned components. In particular, we suggest that the AGN activity does not cause the counterrotation, although it is efficiently triggered by the retrograde gas infall, and it correlates well with the misaligned component appearance. We also find evidence of the stellar disc heating visible as an increase of the vertical-to-radial velocity dispersion ratio above unity in both co- and counterrotating components, which implies the importance of the kinematical misalignment in shaping the velocity ellipsoids in disc galaxies., Comment: 20 pages, 20 figures, accepted for publication in MNRAS

Stellar halos in early-type galaxies (ETGs) are shaped by their accretion and merger histories. We use a sample of 1114 ETGs in the TNG100 simulation with stellar masses $10^{10.3}\leq M_{*}/M_\odot\leq 10^{12}$, selected at z=0 within the range of g-r colour and lambda-ellipticity diagram populated by observed ETGs. We study how the rotational support and intrinsic shapes of the stellar halos depend on the fraction of accreted stars, overall and separately by major, minor, and mini mergers. Accretion histories in TNG100 ETGs as well as the radial distributions of ex-situ stars $f_{ex}(R)$ strongly correlate with stellar mass. Low-mass ETGs have characteristic peaked rotation profiles and near-oblate shapes with rounder halos that are completely driven by the in-situ stars. At high $f_{ex}$ major mergers decrease the in-situ peak in rotation velocity, flatten the $V_{*}/\sigma_{*}(R)$ profiles, and increase the triaxiality of the stellar halos. Kinematic transition radii do not trace the transition between in-situ and ex-situ dominated regions, but for systems with $M_{*}>10^{10.6}M_\odot$ the local rotational support decreases with the local ex-situ fraction $f_{ex}(R)$ and triaxiality increases with $f_{ex}$. These correlations are followed by fast and slow rotators alike with a continuous and overlapping sequence of properties. Merger events dynamically couple stars and dark matter: in high mass ETGs and at large radii where $f_{ex}\gtrsim0.5$, both components tend to have similar intrinsic shapes and rotational support, and nearly aligned principal axes and spin directions. Based on these results we suggest that extended photometry and kinematics of massive ETGs ($M_{*}>10^{10.6}M_\odot$) can be used to estimate the local fraction of ex-situ stars and to approximate the intrinsic shapes and rotational support of the co-spatial dark matter component. [abridged], Comment: 22 pages, 17 figures, accepted for publication on A&A update version including a new section on high-redshift compact progenitor galaxies and a new Appendix with with tables reporting relevant physical relations derived throughout this paper

M105 (NGC 3379) is an early-type galaxy in the Leo I group. This group is the nearest group that contains all main galaxy types and can thus be used as a benchmark to study the properties of the intra-group light (IGL) in low-mass groups. We use PNe as discrete stellar tracers of the diffuse light around M105. PNe were identified on the basis of their bright [OIII]5007 AA emission and the absence of a broad-band continuum. We compare the PN number density profile with the galaxy surface-brightness profile decomposed into metallicity components using published HST photometry in two halo fields. We identify 226 PNe candidates within a limiting magnitude of mlim = 28.1 from our Subaru-SuprimeCam imaging, covering 67.6 kpc along the major axis of M105 and the halos of NGC 3384 and NGC 3398. We find an excess of PNe at large radii compared to the stellar surface brightness profile from broad-band surveys. This excess is related to a variation in the luminosity-specific PN number $\alpha$ with radius. The $\alpha$-parameter value of the extended halo is more than 7 times higher than that of the inner halo. We also measure an increase in the slope of the PN luminosity function at fainter magnitudes with radius. We infer that the radial variation of the PN population properties is due to a diffuse population of metal-poor stars ([M/H] < -1.0) following an exponential profile, in addition to the M105 halo. The spatial coincidence between the number density profile of these metal-poor stars and the increase in the $\alpha$-parameter value with radius establishes the missing link between metallicity and the post-AGB phases of stellar evolution. We estimate that the total bolometric luminosity associated with the exponential IGL population is 2.04x10^9 Lsun as a lower limit, corresponding to an IGL fraction of 3.8%. This work sets the stage for kinematic studies of the IGL in low-mass groups., Comment: 22 pages, 18 figures, accepted for publication in Astronomy & Astrophysics. Abridged abstract

We characterize the photometric and kinematic properties of simulated early-type galaxy (ETG) stellar halos, and compare them to observations. We select a sample of ~1200 ETGs in the TNG100 and TNG50 simulations, spanning a stellar mass range of $10^{10.3}-10^{12}M_{\odot}$ and within the range of (g-r) colour and lambda-ellipticity diagram populated by observed ETGs. We determine photometric parameters, intrinsic shapes, and kinematic observables in their extended stellar halos. We study the variation in kinematics from center to halo and connect it to a change in the intrinsic shape of the galaxies. We find that the simulated galaxy sample reproduces the diversity of kinematic properties observed in ETG halos. Simulated fast rotators (FRs) divide almost evenly in one third having flat lambda profiles and high halo rotational support, a third with gently decreasing profiles, and another third with low halo rotation. Slow rotators (SRs) tend to have increased rotation in the outskirts, with half of them exceeding lambda=0.2. For $M_{*}>10^{11.5}M_{\odot}$ halo rotation is unimportant. A similar variety of properties is found for the stellar halo intrinsic shapes. Rotational support and shape are deeply related: the kinematic transition to lower rotational support is accompanied by a change towards rounder intrinsic shape. Triaxiality in the halos of FRs increases outwards and with stellar mass. Simulated SRs have relatively constant triaxiality profiles. Simulated stellar halos show a large variety of structural properties, with quantitative but no clear qualitative differences between FRs and SRs. At the same stellar mass, stellar halo properties show a gradual transition and significant overlap between the two families, despite the clear bimodality in the central regions. This is in agreement with observations of extended photometry and kinematics. [abridged], Comment: accepted for publication in A&A, 25 pages, 22 figures

Cosmological simulations predict that early-type galaxies (ETGs) are the results of extended mass accretion histories. The latter are characterized by different numbers of mergers, mergers mass ratios and gas fractions, and timing. Depending on the sequence and nature of these mergers that follow the first phase of the in-situ star formation, these accretion histories may lead to ETGs that have low or high mass halos, and that rotate fast or slow. Since the stellar halos maintain the fossil records of the events that led to their formation, a discontinuity may be in place between the inner regions of ETGs and their outer halos, because the time required for the halos stars to exchange their energies and momenta is very long compared with the age of these systems. Exquisite deep photometry and extended spectroscopy for significant samples of ETGs are then used to quantify the occurrence and significance of such a transition in the galaxies structural and kinematical parameters. Once this transition radius is measured, its dependency with the effective radius of the galaxies light distribution and total stellar masses can be investigated. Such correlations can then be compared with the predictions of accreted, i.e. ex-situ vs. in-situ components from cosmological simulations to validate such models., Comment: 7 pages, 3 figures, Conference Proc. IAU Symp. 353, M. Valluri and J.A. Sellwood eds

The age-velocity dispersion relation is an important tool to understand the evolution of the disc of the Andromeda galaxy (M31) in comparison with the Milky Way. We use Planetary Nebulae (PNe) to obtain the age-velocity dispersion relation in different radial bins of the M31 disc. We separate the observed PNe sample based on their extinction values into two distinct age populations. The observed velocities of our high- and low-extinction PNe, which correspond to higher and lower mass progenitors respectively, are fitted in de-projected elliptical bins to obtain their rotational velocities, $V_{\phi}$, and corresponding dispersions, $\rm\sigma_{\phi}$. We assign ages to the two PNe populations by comparing central-star properties of an archival sub-sample of PNe, having models fitted to their observed spectral features, to stellar evolution tracks. For the high- and low-extinction PNe, we find ages of $\sim2.5$ Gyr and $\sim4.5$ Gyr respectively, with distinct kinematics beyond a deprojected radius R$\rm_{GC}= 14$ kpc. At R$\rm_{GC}$=17--20 kpc, which is the equivalent distance in disc scale lengths of the Sun in the Milky Way disc, we obtain $\rm\sigma_{\phi,~2.5~Gyr}= 61\pm 14$ km s$^{-1}$ and $\rm\sigma_{\phi,~4.5~Gyr}= 101\pm 13$ km s$^{-1}$. The age-velocity dispersion relation for the M31 disc is obtained in two radial bins, R$\rm_{GC}$=14--17 and 17--20 kpc. The high- and low-extinction PNe are associated with the young thin and old thicker disc of M31 respectively, whose velocity dispersion values increase with age. These values are almost twice and thrice that of the Milky Way disc stellar population of corresponding ages. From comparison with simulations of merging galaxies, we find that the age-velocity dispersion relation in the M31 disc measured using PNe is indicative of a single major merger that occurred 2.5 -- 4.5 Gyr ago with an estimated merger mass ratio $\approx$ 1:5., Comment: 8 pages, 9 figures, accepted for publication in A&A

In this contribution we report on a kinematic study for 33 early type galaxies (ETGs) into their outer halos (average 6 effective radii, Re). We use planetary nebulae (PNe) as tracers of the main stellar population at large radii, where absorption line spectroscopy is no longer feasible. The ePN.S survey is the largest survey to-date of ETG kinematics with PNe, based on data from the Planetary Nebula Spectrograph (PN.S), counter-dispersed imaging, and high-resolution PN spectroscopy. We find that ETGs typically show a kinematic transition between inner regions and halos. Slow rotators have increased rotational support at large radii. Most of the ePN.S fast rotators show a decrease in rotation, due to the fading of the stellar disk in the outer, more slowly rotating spheroid. 30% of these fast rotators are dominated by rotation also at large radii, 40% show kinematic twists or misalignments, indicating a transition from oblate to triaxial in the halo. Despite this variety of kinematic behaviors, the ePN.S ETG halos have similar angular momentum content, independently of fast/slow rotation of the central regions. Estimated kinematic transition radii in units of Re are ~1-3 Re and anti-correlate with stellar mass. These results are consistent with cosmological simulations and support a two-phase formation scenario for ETGs., Comment: 5 pages, 1 figure, proceeding of the 6th Focus Meeting (FM6) at the XXXth IAU GA 2018

We present a kinematic study of a sample of 298 planetary nebulas (PNs) in the outer halo of the central Virgo galaxy M87 (NGC~4486). The line-of-sight velocities of these PNs are used to identify sub-components, to measure the angular momentum content of the smooth M87 halo, and to constrain the orbital distribution of the stars at these large radii. We find that the surface density profile of the ICPNS at $100\, \mathrm{kpc}$ radii has a shallow logarithmic slope, $-\alpha_{\rm ICL}\simeq -0.8$, dominating the light at the largest radii. The angular momentum-related $\lambda_{\mathrm{R}}$ profile for the smooth halo remains in the slow rotator regime, out to 135 kpc average ellipse radius. The velocity dispersione profile is measured to decrease sharply to $\sigma_{\mathrm{halo}} \sim 100\, \mathrm{kms}^{-1}$ at $R_{\rm avg}=135$ kpc. Such a steeply decline of the velocity dispersion and the surface density profile of the smooth halo can be reconciled with the circular velocity curve inferred from assuming hydrostatic equilibrium for the hot X-ray gas if the orbit distribution of the smooth M87 halo changes strongly from approximately isotropic within $R_{\rm avg}\sim 60$ kpc to very radially anisotropic at the largest distances probed. The extended LOSVD of the PNs in the M87 halo allows the identification of several subcomponents: the ICPNs, the "crown" accretion event, and the smooth M87 halo. In galaxies like M87, the presence of these sub-components needs to be taken into account to avoid systematic biases in estimating the total enclosed mass. The dynamical structure inferred from the velocity dispersion profile indicates that the smooth halo of M87 steepens beyond $R_{\rm avg}=60$ kpc and becomes strongly radially anisotropic, and that the velocity dispersion profile is consistent with the X-ray circular velocity curve at these radii without non-thermal pressure effects., Comment: 24 pages. Accepted on A\&A

M49 (NGC 4472) is the dominant galaxy in subcluster B of the Virgo Cluster, and a benchmark for studying the build-up of the extended halos of brightest group galaxies in the outskirts of galaxy clusters. We investigate the kinematics in the outer halo of M49, look for substructures, and describe the transition to the surrounding intra-group light. As kinematic tracers we use planetary nebulae (PNe), combining kinematics from the extended Planetary Nebula Spectrograph (PN.S) early-type galaxy survey with our recent deep photometric sample. We study the position-velocity-plane for bright and faint PN populations out to 95 kpc radius, and employ a multi-Gaussian model for the velocity distribution to identify stellar populations with distinct kinematics and histories. We report the detection of stellar-kinematic substructure associated with the interaction of M49 with the dwarf irregular galaxy VCC 1249. We find two kinematically distinct PN populations associated with the main M49 halo and the extended intra-group light (IGL). These have velocity dispersions $\sigma_{halo}\simeq 170$ km/s and $\sigma_{IGL} \simeq400$ km/s at 10-80 kpc radii. The overall luminosity profile and velocity dispersion at $\sim80$ kpc are consistent with a flat circular velocity curve extrapolated from X-ray observations. The dispersion of the PNe associated with the IGL joins onto that of the satellite galaxies in subcluster B at $\sim100$ kpc radius. This is the first time that the transition from halo to IGL is observed based on the velocities of individual stars. Therefore the halo of M49, consisting of at least three distinct components, has undergone an extended accretion history within its parent group potential. The blue colours of the IGL component are consistent with a population of stars formed in low-mass galaxies at redshift $\sim0.5$ that has since evolved passively, as suggested by other data., Comment: 10 pages, 7 figures, accepted for publication in Astronomy & Astrophysics

In the hierarchical two-phase formation scenario, the extended halos of early type galaxies (ETGs) are expected to have different physical properties from those of the galaxies' central regions. This work aims at characterizing the kinematic properties of ETG halos using planetary nebulae (PNe) as tracers, which allow us to overcome the limitations of absorption line spectroscopy of continuum at low surface brightness. We present two-dimensional velocity and velocity dispersion fields for 33 ETGs, including both fast (FRs) and slow rotators (SRs), making this the largest kinematic survey to-date of extragalactic PNe. The velocity fields extend out to a median 5.6 effective radii (Re), with a range [3Re-13Re] for the PN.S ETGs. We complemented the PN kinematics with absorption line data from the literature. We find that ETGs typically show a kinematic transition between inner regions and halo. Estimated transition radii in units of Re anticorrelate with stellar mass. SRs have increased but still modest rotational support at large radii, while most of the FRs show a decrease in rotation, due to the fading of the stellar disk in the outer, more slowly rotating spheroid. 30% of the FRs are dominated by rotation also at large radii. Most ETGs have flat or slightly falling halo velocity dispersion profiles, but 15% of the sample have steeply falling profiles. 40% of the FRs show kinematic twists, misalignments, or rotation along two axes, indicating that they turn from oblate at the center to triaxial in the halo, consistently with the distribution of isophote twists from extended photometry. ETGs have more diverse kinematic properties in their halos than in the central regions. The observed kinematic transition to the halo and its dependence on stellar mass is consistent with LambdaCDM simulations and supports a two-phase formation scenario. [abridged], Comment: Accepted for publication in A&A, 52 pages, 30 figures