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]

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.

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.

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.

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.

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]

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.