---
abstract: 'Recent studies adopting $\lambda_{\rm Re}$, a proxy for specific angular momentum, have highlighted how early-type galaxies (ETGs) are composed of two kinematical classes for which distinct formation mechanisms can be inferred. With upcoming surveys expected to obtain $\lambda_{\rm Re}$ from a broad range of environments (e.g. SAMI, MaNGA), we investigate in this numerical study how the $\lambda_{\rm Re}$-$\epsilon_{\rm e}$ distribution of fast-rotating dwarf satellite galaxies reflects their evolutionary state. By combining N-body/SPH simulations of progenitor disc galaxies (stellar mass $\simeq$10$^{\rm 9}$ M$_{\odot}$), their cosmologically-motivated sub-halo infall history and a characteristic group orbit/potential, we demonstrate the evolution of a satellite ETG population driven by tidal interactions (e.g. harassment). As a general result, these satellites remain intrinsically fast-rotating oblate stellar systems since their infall as early as $z=2$; mis-identifications as slow rotators often arise due to a bar/spiral lifecycle which plays an integral role in their evolution. Despite the idealistic nature of its construction, our mock $\lambda_{\rm Re}$-$\epsilon_{\rm e}$ distribution at $z<0.1$ reproduces its observational counterpart from the [ATLAS]{}$^{\rm 3D}$/SAURON projects. We predict therefore how the observed $\lambda_{\rm Re}$-$\epsilon_{\rm e}$ distribution of a group evolves according to these ensemble tidal interactions.'
author:
- |
    C. Yozin[^1] and K. Bekki\
    \
    ICRAR, M468, The University of Western Australia, 35 Stirling Highway, Crawley Western Australia, 6009, Australia
bibliography:
- 'bib.bib'
date: 'Accepted 2015 January. Received 2015 January; in original form 2015 January'
title: The global warming of group satellite galaxies
---

\[firstpage\]

\[3\][\#2]{}

galaxies: interactions – galaxies: dwarf

Introduction
============

In a major step beyond classifying early-type galaxies (ETGs) according to their morphology alone, the spatially-resolved velocity fields arising from recent volume-limited IFS (integral-field spectroscopy) surveys (SAURON, [[ATLAS]{}$^{\rm 3D}$]{}, CALIFA) have provided deeper insight into their formation. By consolidating this kinematic data into a proxy ($\lambda_{\rm Re}$) for the specific angular momentum, a metric tightly correlated with morphology [@obre14], it has become possible to distinguish two classes of ETGs: fast and slow rotators [FR and SR respectively; @emse07]. The pressure-supported SRs are predominantly round, old and massive (stellar mass [M$_{\ast}$]{}$>$[10$^{\rm 10.5}$]{} [M$_{\odot}$]{}), whilst FRs appear continuous (in both kinematics and shape) with late-type galaxies [LTGs; @emse11; @kraj11; @weij14].

This dichotomy has parallels with the red-blue classification of galaxies, and likewise, SR ETGs are conceived to be the natural products of their hierarchical assembly [@khoc06]. Myriad numerical studies have addressed this scenario; those recently motivated by [[ATLAS]{}$^{\rm 3D}$]{} have established that both binary major mergers and successive minor mergers can reproduce the kinematic peculiarities of SRs [@bour07; @jess09; @bois11; @mood14]. In practice, the cosmological assembly of a SR galaxy follows a non-trivial evolution in $\lambda_{\rm R}$ depending on merger mass ratios/gas fractions [@naab14].

As IFS surveys broaden their sample size, more focus will shift to the influence of different environments on forming SR/FRs. Early results in this context suggest SRs are preferentially located in dense locations [@capp13], with a clear radial trend in ellipticity among clusters [@deug15]. Mounting evidence also suggests the formation of FR ETGs [such as S0s, the underlying driver of the morphology-density relation; @dres80] most commonly follows the transformation of LTGs in groups [@wilm09; @just10].

The distinct kinematics between LTGs and S0s is a strong argument against fading alone in this transition, although it remains to be established which of the associated group mechanisms dominate. In a precursor to this paper, @bekk11 demonstrated that group tides/harassment can faciliate spiral-to-S0 transformation, and @foga15 showed how the corresponding evolution in $\lambda_{\rm Re}$ is consistent with those of the pilot SAMI survey. On the other hand, @quer15 proposed major mergers as a viable formation pathway, in the context of a group host, for those S0s observed in the CALIFA survey.

To address this issue (and build upon the success of adopting $\lambda_{\rm Re}$ as a discriminant of formation mechanisms), we introduce in this paper a new concept. Based on the premise that environmental mechanisms acting within a group or cluster can be imprinted on the [$\lambda_{\rm Re}$-$\epsilon_{\rm e}$]{} distribution of its satellite galaxies, we describe here the appearance of this [*global warming*]{} from LTGs to ETGs as a function of redshift.

Although tested against the SAURON and [[ATLAS]{}$^{\rm 3D}$]{} samples (the latter comprising 260 galaxies), the work is orchestrated for the purpose of comparison with the upcoming Sydney-AA0 Multi-object IFS survey [SAMI; @croo12]. The complete survey will catalogue 3400 galaxies in the local universe ($z$$<$0.1), with excellent coverage anticipated for galaxies of stellar mass [M$_{\ast}$]{}$=$[10$^{\rm 9-10}$]{} [M$_{\odot}$]{} in group hosts of dynamical mass [10$^{\rm 13}$]{} [M$_{\odot}$]{}. This particular combination incorporates some of the most common galaxies in the universe [@wein09], the most dominant environment to shape cosmic SFH since [$z=1$]{} [@pope15], and is therefore the focus in this present study. Additionally, this builds upon our previous work [@yozi15b hereafter YB15] concerning the long ($\sim$7 Gyr) quenching/transformation timescales of Magellanic Cloud-type dwarf group-satellites as inferred from near-field observations. The paper is organised as follows: in Section 2, we describe the numerical simulations and the post-processing that permits comparison with the [$\lambda_{\rm Re}$-$\epsilon_{\rm e}$]{} distributions of [[ATLAS]{}$^{\rm 3D}$]{}/SAMI galaxies; in Section 3, we describe the results of a parameter study with these models, and Section 4 concludes this study with a discussion on the merits of our proposed concept for understanding the dynamical state of a group galaxy population.

Method
======

![(Top panel) Probabilility distribution of lookback times at which a [M$_{\ast}$]{}$=$[10$^{\rm 9}$]{} [M$_{\odot}$]{} first becomes a satellite [reproduced from @delu12]; (Bottom panel) A schematic of the group model with the virial radius (r$_{\rm vir}$, shown as grey dashed line) bounding the group halo potential (with grey colour scale representing mass density), characteristic galaxy orbit (black solid line) and other satellites (red circles, with size$\propto$mass); (inset panel) The characterstic galaxy orbit, expressed with the group-centric radius as a function of time.](p1.eps){width="1.\columnwidth"}

Our methodology is comprised of two key steps, the first of which is the simulation of a self-consistent N-Body galaxy model with an interstellar medium (ISM) treated with smoothed particle hydrodynamics (SPH). This model is initially placed on a characteristic orbit within a fixed potential representing the group dark matter halo, as defined in YB15 for satellites with stellar mass ([M$_{\ast}$]{}) of [10$^{\rm 9-9.5}$]{} [M$_{\odot}$]{}.

The second step involves the statistical sampling of this simulated galaxy model evolution, such that we can construct a mock satellite population at a given redshift comprised of galaxies with nominally different infall times. We assume here that infall time ($t_{\rm inf}$) refers to the lookback time at which a satellite first crossed the virial radius of its host halo. Furthermore, our choice of $t_{\rm inf}$ is cosmologically motivated insofar as it is drawn from a probability distribution $p$($t_{\rm inf}$) of satellite infall times (with [M$_{\ast}$]{}$\simeq$[10$^{\rm 9}$]{} [M$_{\odot}$]{}; top panel of Fig. 1) corresponding to merger trees constructed from the Millennium Simulation [@delu12].

The main results of this paper, the mock [$\lambda_{\rm Re}$-$\epsilon_{\rm e}$]{} distribution at a given redshift (lookback time $t$), are thus constructed using the Monte-Carlo sampling of $p$($t_{\rm inf}$). Repeated for a minimum of 1000 samplings, we incorporate the simulated galaxy model properties into the mock distribution at a simulation timestep $t_{\rm i} = t_{\rm inf}-t$ if that timestep $t_{\rm i}\geq$0.

These steps are repeated for several galaxy models, differing in properties such as their gas fraction or orbital orientation. This scenario is also compared with those in which other group satellites, represented as point masses, are added to the group model to establish the influence of satellite-satellite tidal interactions. We emphasize therefore that our methodology is based on an idealised group model comprising a fixed dark halo and its satellite population (bottom panel, Fig. 1); the adopted orbit is sufficiently energetic (e.g. the pericentre radius or distance to the group halo centre is sufficiently large) that the explicit modelling of a central galaxy (which by definition lies at the halo centre) is not considered to be necessary.

Numerical models
----------------

  Parameter                                               Fiducial value
  ------------------------------------------------ -----------------------------
  N$^o$. DM particles                                    [10$^{\rm 6}$]{}
  N$^o$. Stellar particles                           3$\times$[10$^{\rm 5}$]{}
  N$^o$. Gas particles                              2.5$\times$[10$^{\rm 5}$]{}
  Stellar mass/[M$_{\odot}$]{} (M$_{\rm d}$)              10$^{\rm 9.25}$
  Gas/stellar mass (M$_{\rm g}$/M$_{\rm d}$)                     2
  DM halo/stellar mass (M$_{\rm h}$/M$_{\rm d}$)                85
  Stellar disc scalelength (r$_{\rm d}$/kpc)                    1.6
  Stellar disc scaleheight (z$_{\rm d}$/kpc)                   0.32
  Gas disc scalelength (r$_{\rm g}$/kpc)                        4.2
  Gas disc scaleheight (z$_{\rm g}$/kpc)                       0.84
  DM Halo virial radius (r$_{\rm vir}$/kpc)                     142
  DM Halo concentration ($c$)                                   6.5
  Metallicity (\[Fe/H\])                                       -0.72
  Metallicity gradient ($\Delta$\[Fe/H\]/kpc)                  -0.04
  SN expansion timescale ($\tau_{\rm SN}$/yr)           [10$^{\rm 5.5}$]{}

  : A summary of fiducial galaxy model properties

We present results from 22 galaxy model simulations, each run with an original parallelised chemodynamical code GRAPE-SPH [@bekk09], combining gravitational dynamics and an ISM with radiative cooling/stellar processing. Our numerical satellite/group models follow that of YB15 closely, to which we refer the interested reader.

To summarise, the galaxy commences at apocentre on a characteristic orbit (constrained to the $x-y$ plane) which is derived from simulations of collisionless satellite infall with a live model of the dark matter (DM) halo [@vill12]. The orientation of the galaxy model is specified by $\theta$ and $\phi$, corresponding to the angle between the $z-$axis and the vector of the angular momentum of the disc, and the azimuthal angle measured from the $x-$axis to the projection of the angular momentum vector of the disc on to the $x-y$ plane, respectively. Our fiducial model adopts $\theta$=45$^{\circ}$ and $\phi$=30$^{\circ}$.

Fig. 1 conveys this orbit and a schematic of the group model which consists of a spherically-symmetric potential alone in the fiducial case (Sections 3.1-2), and incorporates $\sim$50 group satellites (represented in the simulation as point masses) when considering the additional effect of tidal harassment (Section 3.3).

The group potential has total mass of [10$^{\rm 13}$]{} [M$_{\odot}$]{} and structural properties characteristic of [$z=1$]{} [virial radius 380 kpc and a NFW density profile with concentration 4.9; @vill12]. We use Monte-Carlo sampling to fit the point mass satellites to a Schecter function (slope -1.07) function, limited to a luminosity range 0.01 to 2.5 solar and with the assumption of a mean mass-to-light ratio $\sim$40. Their initial orbits are assigned according to a NFW profile with concentration 3.0.

### Fiducial galaxy model

Table 1 summarises the principle properties of our fiducial galaxy model, which consists of a gas-rich bulgeless (bulge-to-total ratio, B/T=0) galaxy as in YB15. The initial gas and halo mass with respect to the stellar mass is selected according to best-fit relations from the ALFALFA survey [@huan12] and @mill14 respectively, while DM halo properties are obtained from mass-redshift-concentration relations of @muno11. The stellar disc scalelength r$_{\rm d}$ is observationally motivated from the size-mass relations of @ichi12, while the gas disc scalelength r$_{\rm g}$ is assumed to be a factor $\sim$3 larger [@krav13]; stellar and gas scaleheights z$_{\rm d}$ and z$_{\rm g}$ are assigned as 0.2r$_{\rm d}$ and 2z$_{\rm d}$ respectively.

While we do not consider the star formation (SF) history here, gas inflow and bulge growth are crucial components to our model evolution, and therefore subgrid SF parametrisations and the efficiency of supernovae (SNe, controlled primarily with the adiabatic expansion timescale $\tau_{\rm SN}$) are constrained according to the Kennicutt-Schmidt law. To offset the idealistic nature of our scenario, we emphasise that our minimum baryonic resolution of $\sim$[10$^{\rm 4}$]{} [M$_{\odot}$]{} is two orders of magnitude less than recent cosmological simulations [e.g. Illustris; @gene15].

### Alternate galaxy models

The simulations constituting our parameter study are summarised in Table 2.

-   To account for pre-processing prior to group infall, we construct [$\lambda_{\rm Re}$-$\epsilon_{\rm e}$]{} distributions for a galaxy with an initial B/T of 0.2, where the bulge has a Hernquist density profile, isotropic velocity dispersion (with radial velocity dispersion assigned according the Jeans equation of a spherical system).

-   To further account for the possibility of efficient gas stripping prior to or upon infall to the group, we consider the evolution of a gas-poor galaxy (gas mass$=$$0.2$[M$_{\ast}$]{}).

-   The dependence on orbital orientation is addressed in three scenarios, including: Orbit:Face-on ($\theta$=90$^{\circ}$, $\phi$=0$^{\circ}$), Orbit:Edge-on ($\theta$=0$^{\circ}$, $\phi$=0$^{\circ}$) and Orbit:Retrograde ($\theta$=-45$^{\circ}$, $\phi$=-30$^{\circ}$).

-   The influence of satellite-satellite interactions (in a satellite population whose mean velocity dispersion is $\sim$150 kms$^{\rm -1}$) is established for the fiducial and B/T=0.2 models with two suites of simulations, each comprising 8 runs where the live galaxy is placed randomly on the locus defined by the group virial radius. This method accounts for the quasi-stochastic influence of the point-mass group satellites on the live galaxy model.

  Simulation            Deviation from fiducial
  --------------------- ------------------------------------
  B/T=0.2               Initial B/T of 0.2
  Orbit:Edge-on         Edge-on disc w.r.t. orbit plane
  Orbit:Face-on         Face-on disc w.r.t. orbit plane
  Orbit:Retrograde      Retrograde w.r.t. orbit plane
  Gas fraction:Low      Initial gas fraction of 0.2
  Fiducial+Harassment   Group satellites added to fiducial
  B/T=0.2+Harassment    Group satellites added to B/T=0.2

  : Summary of 21 non-fiducial models and their deviations from the fiducial; the two entries incorporating harassment are each comprised of eight simulations.

Analysis parameters
-------------------

We deduce both intrinsic and apparent specific angular momenta and ellipticities of our galaxy models. Intrinsic parameters are calculated from the 3D spatial/velocity data of all stellar particles within the half-mass radius, where the intrinsic ellipticity is established from decomposition of the inertial tensor of all stars (from which the ellipticity is 1-$c$/$a$ where $c$ and $a$ are the minor and major axes respectively).

The apparent specific angular momentum is computed according to the relation defined by @emse07: $$\lambda_{\rm R} = \frac{ \sum_{i=1}^{N_{\rm b}} F_{\rm i} R_{\rm i} |V_{\rm i}| } { \sum_{i=1}^{N_{\rm b}} F_{\rm i} R_{\rm i} \sqrt{ V_{\rm i}^{\rm 2} + \sigma_{\rm i}^{\rm 2} } }$$ where F$_{\rm i}$, R$_{\rm i}$, V$_{\rm i}$ and $\sigma_{\rm i}$ are the flux (computed here from the integrated surface brightness), radius, projected stellar velocity and stellar velocity dispersion respectively, in the $i$th spatial bin.

As in YB15, we adopt synthesis models for stellar evolution and an assumed gas-to-dust ratio of 5$\times$[10$^{\rm 21}$]{} cm$^{\rm -2}$A$_{\rm V}$ to establish stellar luminosities from the stellar age/mass of simulation particles. Surface brightness distributions of the galaxy models can be therefore constructed on 2D spatial grids. Ellipticities are calculated according to the [Kinemetry]{} method [@kraj06], in which concentric ellipses are fitted to binned line-of-sight kinematic maps.

@emse07 showed that the apparent properties taken at the effective radius, r$_{\rm e}$ (enclosing half the projected light), are sufficient to characterise inner disc kinematics. From the radial profile of $\lambda_{\rm R}$ and the ellipse ellipticities expressed as a function of their mean radii, we compute $\lambda_{\rm Re}$ and $\epsilon_{\rm e}$ respectively by linear interpolation at r$_{\rm e}$.

In general, our method in obtaining these values is consistent with those of previous studies [e.g. @emse07; @bois11; @naab14; @mood14]; however, we have not adopted their Voronoi binning in construction of the kinematic maps, as our simulation resolution is sufficient to provide the requisite signal-to-noise ratio up to r$_{\rm e}$ [i.e. at least $\sqrt{N}$$\ge$20 per bin; @capp03].

![(Top panel) Second Fourier mode amplitude (A$_{\rm 2}$, black solid line) and bar radius (r$_{\rm bar}$, normalised by stellar scalelength r$_{\rm d}$, black dashed line) as a function of time for the fiducial simulation. The red dashed line at A$_{\rm 2}=0.2$ shows our empirical threshold that separates ETGs and LTGs; (bottom panel) Bulge-to-Total ratio as a function of time for the fiducial and B/T=0.2 simulations, computed as the excess light with respect to an fitted exponential disc component. We qualify this excess as a pseudo-bulge if kinematic maps show it to be rotationally-supported.](p2.eps){width="1.\columnwidth"}

Our mock [$\lambda_{\rm Re}$-$\epsilon_{\rm e}$]{} distributions (at a given lookback time) combine the $\lambda_{\rm Re}$ and $\epsilon_{\rm e}$ values deduced from the corresponding galaxy model data when inclined and viewed at 100 equispaced increments of $i$, where $i$ is the stellar disc inclination $0^{\circ}<i<90^{\circ}$ with respect to face-on. This approach is based on the premise that observed galaxy populations have no preferential viewing inclination.

Furthermore, to compare this data with [[ATLAS]{}$^{\rm 3D}$]{} galaxies (and those of the SAMI pilot survey), whose criterion for ETGs is predicated on the absence of spiral arms, we exclude those $\lambda_{\rm Re}$ and $\epsilon_{\rm e}$ computed at timesteps where the face-on surface brightness distribution has a mean second Fourier amplitude (A$_{\rm 2}$, within 2r$_{\rm e}$) exceeding 0.2. We have confirmed that our criterion does indeed correspond to an apparent absence of spiral structure in mock B-band surface brightness imaging of our galaxy models down to $\sim$28 mag arcsec$^{\rm 2}$.

Fig. 2 compares this A$_{\rm 2}$ threshold with the computed A$_{\rm 2}$ of the fiducial model as a function of time. Also shown is the stellar bar radius, computed as in YB15. We find the model to be unstable to $m=2$ perturbations following tidal interactions at successive orbital pericentres. Subsequent gas inflow and angular momentum exchange to the bar can destroy it [@frie93; @bour05a], and therefore, the satellite undergoes multiple cycles between ETG and LTG classification.

Concurrently, the bulge-to-total ratio (which we calculate from the excess light with respect to a fitted exponential disc; see YB15 for more details) grows monotonically with time (Fig. 2, bottom panel). We assert that this constitutes a pseudo-bulge for the initially bulgeless Fiducial model because the corresponding kinematic fields strongly imply its rotational-support [@korm04].

Results
=======

![A mock distribution of apparent specific angular momentum, $\lambda_{\rm Re}$ vs. apparent ellipticity, $\epsilon_{\rm e}$ (both measured at at r$_{\rm e}$) for the Fiducial galaxy model at [$z=1.0$]{} (inset panel) and $<0.1$ (main panel). The 2D probability density plot in each panel shows the normalised distribution of ETGs (as classified with A$_{\rm 2}$). In the bottom panel, this is compared with the full LTG+ETG population (green contours), together with those FRs identified in [[ATLAS]{}$^{\rm 3D}$]{} [@emse11] to appear as SRs (white symbols). The white dashed line shows an empirical boundary between SR/FRs ($\lambda_{\rm Re}=0.31\epsilon_{\rm e}^{\rm 0.5}$). The red solid line shows show an isotropic stellar system viewed edge-on would appear in this plot; the red dashed lines show how this varies with viewing inclination (where face-on lies at the origin).](p3.eps){width="1.\columnwidth"}

Fig. 3 shows the mock probability distribution of [$\lambda_{\rm Re}$-$\epsilon_{\rm e}$]{} for our fiducial simulation at redshifts [$z=1$]{} and $<$0.1 (the priority coverage of the full SAMI survey). This distribution is comprised of galaxies with stellar masses in the range [10$^{\rm 9-9.5}$]{} [M$_{\odot}$]{} due to the combined effect of secular/triggered SF and gas stripping acting on a disc with initial mass [M$_{\ast}$]{}$=$[10$^{\rm 9.25}$]{} [M$_{\odot}$]{}. Overlaid upon each density plot is the theoretical properties of an edge-on stellar system (solid red line) assuming a two-integral dynamical Jeans model [@capp07; @emse11], and its corresponding projection towards face-on (dashed lines, towards the origin).

We find at an early epoch that our (predominantly Sc/d) satellite population is spread along along the upper-most dashed arcs (corresponding to an intrinsic $\epsilon_{\rm e}=0.7-0.8$). At lower redshifts, this distribution becomes more concentrated with a clear peak forming by $z<0.1$ at $\epsilon_{\rm e}=0.5$, $\lambda_{\rm Re}=0.5$. Fig. 2 (top panel) shows that in the intervening time the B/T ratio grows from 0 to 0.25; the corresponding (pseudo-)bulge is rotationally-supported and dominated by young stars.

An outlying proportion of the ETG population falls below an empirical criterion from [[ATLAS]{}$^{\rm 3D}$]{} that denotes slow rotation ($\lambda_{\rm Re}=0.31\sqrt{\epsilon_{\rm e}}$). Since we enforce a consistent definition of [$\lambda_{\rm Re}$]{} based on apparent properties at r$_{\rm e}$, we expect some apparent SR classifications to arise due to a stellar bar. This bar, the size of which diminishes with time after a peak at first pericentre (Fig. 2), distorts the ellipicity metric in viewing inclinations close to face-on [@capp07].

Our autonomous ETG criterion based on A$_{\rm 2}$ filters out many more such SR classifications arising from the full LTG+ETG population (Fig. 3, green contours); this distribution shows a [*tendril*]{} at $\epsilon_{\rm e}=0.6, \lambda_{\rm Re}$=0.3 which lies outside the theoretical bounds for an oblate axisymmetric system, and corresponds to strong tidal features formed during the first orbital pericentre.

The capacity for group tidal interactions to induce non-equilibrium kinematics in the satellite is demonstrated here also in terms of intrinsic properties (Fig. 4). Here we find in general a close concordance between our simulation of a single fiducial galaxy model (with infall at 10 Gyr ago) and the best-fit relation for a two-integral Jeans dynamical system [@capp07; @emse11]. This applies across a wide range of $\epsilon_{\rm e}$, albeit with some significant transient deviations at orbital pericentres (e.g. $\epsilon_{\rm e}=0.5, \lambda_{\rm Re}$=0.5).

The existence of apparent SRs are therefore traced to transient perturbations to an otherwise persistently disc-like kinematic field [@kraj11; @barr15], and also projection effects including that due to bars which we note are undetected by photometry in as many as 10 percent of galaxies in the [CALIFA]{} survey [@holm15]. The location of apparent SRs (at $\epsilon_{\rm e}=0.2$, $\lambda_{\rm Re}=0.1$; Fig. 3) coincides with the simulated remnants of late mergers [@naab14], and observed [[ATLAS]{}$^{\rm 3D}$]{} galaxies whose line-of-sight velocity distributions show kinematically-decoupled components (KDCs), non-regular rotation (NRR) or double-sigma components (2$\sigma$) as opposed to pressure-supported kinematics [white symbols; @emse11].

![The intrinsic evolution of the Fiducial simulation (black filled circles) and B/T=0.2 (red open squares)in the [$\lambda_{\rm Re}$-$\epsilon_{\rm e}$]{} plane; black labels give an approximate time since the start of the simulation. The approximate best-fit relation ($\lambda_{\rm Re}=0.93\epsilon_{\rm e}^{\rm 0.5}$) to the expected evolution of an isotropic stellar system is conveyed with the grey dashed line.](p4.eps){width="1.\columnwidth"}

Parameter dependencies
----------------------

Fig. 5 (top panel) shows how the $z<0.1$ [$\lambda_{\rm Re}$-$\epsilon_{\rm e}$]{} distribution depends on our initial galaxy model parameters. The key results are summarised as follows:

-   Our B/T=0.2 model displays an concentration lying closer to the origin than the fiducial case. This is to be expected from an inner disc hosting a classical (pressure-supported) bulge in addition to the rotation-supported pseudo-bulge formed also in the fiducial case (with a final B/T=0.4; Fig. 2). Though not shown in this figure, we note that fewer apparent SRs are detected due in part to an inner disc less favourable for the radial ($x_{\rm 1}$) stellar orbits that constitute a stellar bar [@norm96].

-   Gas fraction:Low does not form a pseudo-bulge as efficiently as the fiducial case, where in the latter case the loss in gas angular momentum by successive tidal torques fuels bulge growth [@barn91]. In spite of this, the [$\lambda_{\rm Re}$-$\epsilon_{\rm e}$]{} distribution shows a small deviation towards higher $e_{\rm e}$, a result similar to that found among the simulated remnants of dissipative/collisonless mergers [@jess09]. Of further note is that the relatively small dissipative component allows the triggered bar (and its associated distortion of $\epsilon_{\rm e}$) to persist on order of a Hubble time causing a higher incidence of SR mis-identifications.

-   The [$\lambda_{\rm Re}$-$\epsilon_{\rm e}$]{} distribution shows a clear dependence on the satellite orbit, a proxy to the dependence on tidal torque efficiency. Orbit:Face-on is not sufficiently perturbed to facilitate much morphological/kinematic evolution, and there exists negligible variation between prograde or retrograde orbits. In general, we find [$\lambda_{\rm Re}$]{} and $\epsilon_{\rm e}$ drop sharply at orbital pericentres, as a function of time, by values that scale with ${\sim}cos(i)$ (where $i$ is galaxy inclination to its orbit).

Dependence on tidal harassment
------------------------------

Fig. 5 (Bottom panel) shows the [$\lambda_{\rm Re}$-$\epsilon_{\rm e}$]{} distribution for Fiducial + Harassment at $z<0.1$. The addition of point-mass group satellites causes the dispersion of the highly concentrated distribution associated with our fiducial galaxy model at an intermediate orbital inclination (top panel). This dispersion, while largely remaining within theoretical limits (red lines), is due to the quasi-stochasticity of satellite-satellite interaction impact parameters and inclinations. As in the case of simulated minor mergers wherein kinematically-distinct stellar components often arise [@bois11], the proportion of SR misidentifications is more significant than in our fiducial case (Fig. 3).

![(Top panel) As Fig. 3, but where the [$\lambda_{\rm Re}$-$\epsilon_{\rm e}$]{} distribution of the fiducial simulation at $z<0.1$ (2D probability density plot) is compared with simulations commencing with differing galaxy model conditions (contours shown at $p=0.75$ and 0.875); (bottom panel) As Fig. 3, but showing the ensemble [$\lambda_{\rm Re}$-$\epsilon_{\rm e}$]{} distribution of Harassment at $z<0.1$. Green contours represent the mock distribution of 50,000 FR ETGs [@emse11] based on the Monte-Carlo method of @capp07.](p5.eps){width="1.\columnwidth"}

Discussion and conclusions
==========================

Although $\lambda_{\rm Re}$ has already proven a useful metric for distinguishing the formation paths of ETGs, galaxy evolution on the [$\lambda_{\rm Re}$-$\epsilon_{\rm e}$]{} plane remains a complex, non-linear process within a $\Lambda$CDM cosmology [@naab14]. In this study, we demonstrate how it is instructive to consider the global [$\lambda_{\rm Re}$-$\epsilon_{\rm e}$]{} distribution of many galaxies, analysed in a consistent manner and including projection effects. In environments like groups, a target for the upcoming SAMI survey among others, this method can overcome the non-monotonic evolution of individual satellites. Instead, we can describe the effective redshift-evolution of the group itself, which we euphemistically describe as undergoing a process of [*global warming*]{}.

This demonstration was performed with a suite of idealised simulations derived from those of @yozi15b and thus shares the same limitations, including a fixed spherical group potential. We argue that since the intrinsic evolution of our satellites generally follow that of an oblate FR (Fig. 4), they will tend towards a similar [$\lambda_{\rm Re}$-$\epsilon_{\rm e}$]{} distribution albeit at a rate determined by orbital parameters. This particular evolution is predicated on an absence of major mergers, which we deem rare for the [M$_{\ast}$]{}$=$[10$^{\rm 9}$]{} [M$_{\odot}$]{} satellites considered here, in accordance with the ILLUSTRIS simulation which predicts $\sim$0.1 such events per satellite since [$z=2$]{} [@rodr15].

Significant minor mergers (e.g. mass ratios of 1:10 to 1:6) will be more frequent, however; previous numerical studies [e.g. @bour05b] have suggested these events yield remnants analogous to S0s, and comparable therefore to the FR ETGs produced by our methodology (where bulge growth is likewise promoted by in-situ SF). However, this parallelism should be explicitly challenged in a future study.

Hydrodynamical interactions between the satellite and the intra-group medium are also ignored; our results and those of @bois11 have suggested that the extant gas fraction has only a small impact on [$\lambda_{\rm Re}$-$\epsilon_{\rm e}$]{}, although smooth accretion has been shown to raise [$\lambda_{\rm Re}$]{} if suitably aligned with the disc [@chri15]. We do not expect, however, the environmental mechanism of ram pressure stripping will significantly influence a kinematics-driven evolution of [$\lambda_{\rm Re}$-$\epsilon_{\rm e}$]{}; in fact, this can be exploited in efforts to discern galaxy formation mechanisms [e.g. S0s, see also @quer15].

Our models also lack an explicit model of galactic outflow, which recent cosmological simulations demonstrate as a key actor in regulating angular momentum [@gene15]. Our preliminary simulations suggest the effect of stronger feedback in our model (applied with a longer adiabatic timescale, e.g. $\tau_{\rm SN}={\rm 10}^{\rm 6}$ yr) merely controls the rate of [$\lambda_{\rm Re}$-$\epsilon_{\rm e}$]{} evolution rather than its shape, the latter being the focus of this present study.

Lastly, several studies have noted that SR formation is highly resolution dependent [@bois10], insofar as convergence requires a baryonic mass resolutions of at least $\sim$[10$^{\rm 4}$]{} [M$_{\odot}$]{}; indeed, the merger simulations of @mood14 adopt this threshold and find noise fluctuations cause $\lambda_{\rm Re}$ to deviate by only $\sim$0.01. Since our simulations also lie at this threshold, we maintain our claim that the group mechanisms modelled here produce only FR ETGs.

![As Fig. 3, but for mock [$\lambda_{\rm Re}$-$\epsilon_{\rm e}$]{} distributions derived from Fiducial + Harassment and B/T=0.2 + Harassment at [$z=1.0$]{} and $<0.1$. White contours are chosen at $p$=0.7 and 0.85 to illustrate the distinct evolutions associated with disc-dominated FR ETGs (B/T=0 and B/T=0.2) over $\sim$8 Gyr.](p6.eps){width="1.\columnwidth"}

In an attempt to verify our idealised numerical methodology, we compare our results in Fig. 5 (bottom panel) to the mock distribution of 50,000 FR ETGs galaxies obtained by @emse11. Their Monte-Carlo method follows from a linear relation found between intrinsic $\epsilon$ and anisotropy among [SAURON]{} galaxies [@capp07], for which the associated scatter is modelled as a Gaussian [see also @weij14]. The resulting mock distribution (green contours) shows qualitative agreement with both the [[ATLAS]{}$^{\rm 3D}$]{} sample of FRs (which we note is drawn from volume-limited surveys and therefore comprises satellite and central galaxies), and our $z<0.1$ snapshot of Fiducial + Harassment.

Specifically, Fiducial + Harasment reproduces the spread of [[ATLAS]{}$^{\rm 3D}$]{} galaxies and their concentration around $e_{\rm e}=0.5, \lambda_{\rm Re}=0.5$. Since the majority of the [[ATLAS]{}$^{\rm 3D}$]{} galaxies were identified by morphological type $T$ as S0s, we can claim therefore that our simulations capture the fundamental mechanisms that yield a S0 morphology. Fig. 6 (bottom panels) illustrates how those FR ETGs galaxies hosting an isotropic bulge (classically the remnant of merger events) will appear to be more concentrated towards the origin; this [$\lambda_{\rm Re}$-$\epsilon_{\rm e}$]{} distribution is therefore qualitatively consistent with that of the pilot-SAMI ETGs identified as cluster members by @foga15.

Of course, a comparison to the [[ATLAS]{}$^{\rm 3D}$]{}/SAURON and pilot-SAMI catalogs is not strictly correct (given their generally more massive stellar masses with respect to that represented by our models) and thus we look forward to the full SAMI survey for future comparisons.

However, this cursory verification of our low-redshift mock [$\lambda_{\rm Re}$-$\epsilon_{\rm e}$]{} distributions of ETG FRs allows us to propose, firstly, that the observed data (for a given stellar mass) represents the superposition of [$\lambda_{\rm Re}$-$\epsilon_{\rm e}$]{} distributions specific to particular morphological groups. @capp13b have already shown the bulge fraction can be more reliably inferred from kinematics than bulge-disc decompositions; our broader assertion is that the decomposition of [$\lambda_{\rm Re}$-$\epsilon_{\rm e}$]{} will in principle provide insight as to the evolutionary processes involved in those satellite populations [e.g. pre-processing, and the relative role of classical and pseudo-bulge building mechanisms; @delu12; @wetz13].

This concept complements that of the fundamental/mass plane (FP), a tight scaling relation for ETGs, which the SAMI survey has already proven capable in its construction [@scot15]. Unlike our concept, however, it is not conclusive if properties of the FP (e.g. slope, offset) vary with environment [e.g. @hou15], and indeed its utility as a distance indicator would suffer if it did [@mago12].

A second proposition is that the tidal-driven evolution of a satellite population drives an increasing spheroidal fraction with time, which manifests qualitatively as an increasing concentration of the [$\lambda_{\rm Re}$-$\epsilon_{\rm e}$]{} distribution. This can therefore be exploited to parametrise the [*global warming*]{} of the environment. While admittedly a non-trivial exercise given a cosmological assembly bias, the well-known environment-driven evolution of dwarf satellites since [$z=1$]{} makes them the ideal targets for this concept [@geha12]. Indeed, we advocate the comparison of [$\lambda_{\rm Re}$-$\epsilon_{\rm e}$]{} from high redshift IFU surveys [e.g. the KMOS Redshift One Spectroscopic Survey, presently conducted to compare star forming conditions within 800 star forming galaxies at $z=1$; @stot16] with local counterparts to verify the redshift-dependent behaviour illustrated in Fig. 6.

Acknowledgements {#acknowledgements .unnumbered}
================

We thank the anonymous referee for their comments which improved this paper. Part of this work was supported by the Australian Postgraduate Award Scholarship. The authors also wish to thank members of the Computational Theory Group at ICRAR for their useful comments.

\[lastpage\]

[^1]: [**E-mail**]{} 21101348@student.uwa.edu.au; kenji.bekki@uwa.edu.au
