Properties of the Photometric Components of Lenticular Galaxies

After almost one century since their definition, the origin of lenticular galaxies is still a matter of debate. Several formation processes have been proposed in order to explain the wide variety of their observational properties. These properties could indicate that the class of lenticulars is a family formed by galaxies with different formations and evolutions. Here I review the main observational properties of the photometric components of lenticular galaxies reported in recent decades. I revise the main processes proposed in order to explain their origin and evolution. I also explain the different properties of the lenticular galaxies that have evolved through each of these formation processes. A unique opportunity for understanding the origin of S0 galaxies will shortly be forthcoming. This is due to the morphological classifications of large galaxy samples that have recently been published. These classifications have given us our first ever opportunity to study large samples of lenticulars within a wide range of masses and located in a great variety of environments. These large samples will provide us with a real census of nearby lenticular galaxies and could be crucial in finally helping us to understand the origin and evolution of these galaxies.


Introduction
Lenticular galaxies were introduced by Hubble in 1936 in order to fill the gap between elliptical and spiral galaxies.Their observational properties have been analyzed for decades.Several studies have shown that lenticular galaxies exhibit a great diversity of properties, being similar to both elliptical and spiral galaxies.This makes their origin and evolution still a matter of debate, and several questions still need to be answered.Are lenticular galaxies a well-defined and homogeneous class of galaxies, or are they on the contrary a morphological class of galaxies containing galaxies with different origin and evolution?
Concerning their formation, are lenticulars formed by major mergers of galaxies, similarly to bright elliptical galaxies, or did they form through slow galaxy processes including minor satellite accretion or other secular evolution processes?Are they the final steps in the evolution of late-type galaxies due to environmental mechanisms?Some light can be shed on these questions by analyzing the different observational properties of the structural components present in lenticulars.This is the main aim of this paper.I have compiled information from the literature about the main observational properties of the photometric components of lenticulars and I review the different formation mechanisms proposed.There is a wide diversity of formation mechanisms of S0 galaxies that act on different time scales and galaxy environments.For each formation mechanism, I review the expected properties of the lenticular galaxies formed.
This paper is organized as follows.In Section 2, I have revised visual and automatic classifications of lenticular galaxies.Multicomponent surface brightness decompositions of lenticular galaxies are reviewed in Section 3. The observational properties of the structural components of S0 galaxies are given in Section 4. The special case of dwarf lenticular galaxies (dS0s) is reviewed in Section 5.The different formation mechanisms of S0 galaxies are described in Section 6.Finally, some conclusions and future perspectives are given in Section 7.

Classification of Lenticular Galaxies:
Visual versus Automatic 2.1.Visual Galaxy Classification: Early Definitions of Lenticular Galaxies.Galaxy classification is a visual task.Several visual classifications have been proposed since the discovery of galaxies (extragalactic nebulae) at the beginning of the 20th century (see [1,2] for good reviews on galaxy classifications).One of the pioneering classification systems was proposed by Reynolds [3].He classified the spiral galaxies within a sequence of shapes that he defined as classes I to VII.Six years later, Hubble published the first version of his famous classification scheme.Galaxies were divided in two main classes: ellipticals (Es) and spirals.The elliptical galaxies were located in seven groups according to their flatness from E0 (round Es) to E7 (the most flattened Es).In addition, the spiral galaxies were put into adjacent classification groups depending on the bulge and the spiral arm relationships (Sa, Sb, and Sc groups).Hubble established a sequence of shapes from E0 to Sc with continuity between E7 and Sa galaxies [4].A decade later, Hubble published a new version of his classification system and introduced the type S0 (also called lenticulars by de Vaucouleurs [5]) to fill the gap between E7 and Sa galaxies [6].Lenticular galaxies were thus those showing a prominent central spheroidal component, similar to Es, but also a disc-like structure, similar to spirals, but without spiral arms.The properties of S0 galaxies were later updated in several works (see [7][8][9]).The Hubble classification is usually represented in a schematic way by the so-called Hubble tuning fork diagram (see Figure 1) Several galaxy catalogues show visual classifications of galaxies with a high level of detail.These collections contain examples of the different morphological types of galaxies, including lenticulars, and give the classification criteria for each galaxy type.Among the earlier generations of catalogues, we may mention the Hubble Atlas [10], the Revised Shapley-Ames Catalog [8], and the de Vaucouleurs Atlas [11,12].Although the definition of S0 galaxies is well established in these catalogues, visual differentiation between Es and S0s is not an easy task.Unless typical discs structures such as bars or dust lanes are present in the photometrical images.Figure 2 shows examples of four nearby S0 and E galaxies.In this case, the two lenticular galaxies show strong dust lanes.In several papers, van den Bergh has argued that not all S0 galaxies are properly classified as intermediate between E and Sa galaxies (see, e.g., [13][14][15][16]).He found that the frequency distribution of galaxy luminosities for S0 galaxies is not intermediate between E and Sa objects.Lenticulars are about one magnitude fainter than E and Sa galaxies.He proposed that S0 galaxies are a class of galaxies that includes different types of objects that only have a morphological similarity.The misclassification of E and S0 galaxies has been also pointed out in several recent papers from the SAURON and ATLAS-3D groups (e.g., [17][18][19]), who have studied the 2D kinematics of a large sample of nearby galaxies visually classified as Es and S0s.van den Bergh [20] even proposed a new classification system of galaxies in which spirals and lenticulars formed parallel sequences.This has recently been shown by different groups of researchers (see, e.g., [18,21,22]).
These studies indicate how difficult the classification of E and S0 galaxies is and how difficult it is to select galaxy samples containing only E or S0 galaxies using only photometric information (see [23] and references therein).Fortunately, the number of large samples of S0 galaxies including morphological and/or kinematical properties have increased during the last decade.I mention here the recent near-IR S0galaxy atlas of S0 galaxies (NIRSOS) by Laurikainen et al. (see [24]).It includes detailed morphological information of a sample of 185 S0 galaxies.The SAURON and the ATLAS-3D projects (e.g., [17][18][19]) have also provided detailed kinematic studies of large samples of Es and S0 galaxies, which have been used for a new classification of early-type galaxies according to their kinematical properties.Figure 3 shows the new classification scheme proposed by the ATLAS 3D project.In this classification lenticular galaxies are fast rotators that form a parallel sequence to spiral galaxies.

Automatic Galaxy Classifications.
The main problem with visual galaxy classifications is that they are not reproducible.Although the classifications are done by specialists, visual galaxy classifications are not universal due to the lack of objectivity of the human eye.Automatic algorithms can provide reproducible galaxy classifications and give important clues to the origin and evolution of galaxies.Automatic classifications also have the advantage that biases can be fully understood and characterized based on extensive simulations of galaxies treated as real ones.Nevertheless, automated visual galaxy classifications can only provide broad galaxy classes.They are not as accurate as visual classifications made by trained observers.Another problem with visual classifications is the time needed for detailed classifications of large galaxy samples.Present galaxy surveys such as Sloan Digital Sky Survey (SDSS; [25]) or Two Degree Field (2dF; [26]) have greatly increased the size of previous galaxy samples.Early galaxy catalogues contain several dozens, or even hundreds, of galaxies.In contrast, galaxy samples today can contain millions of galaxies.It is hard to give detailed morphological classifications of such large galaxy samples.Nevertheless, some contributions have been made in the last few years, and detailed visual morphological classifications of local galaxies are now available for small subsamples of the SDSS galaxy catalogue.Thus, Fukugita et al. [27] have visually classified a sample of 2253 galaxies, and Nair and Abraham [28] have also given detailed morphological types of 14 034 nearby galaxies.Special mention should be made of the Galaxy Zoo project [29], which has provided visual classifications (not so detailed as the previously mentioned ones) for more than 4 × 10 7 galaxies by involving about 10 5 classifiers [29,30].
The amount of time needed for detailed visual classifications of galaxies makes the use of automated algorithms a must.Automatic classifications of galaxies can be divided into two broad groups: parametric and nonparametric.Parametric classifications reproduce galaxy measurements using parametric laws.One of the most extended parametric classifications classifies galaxies according to some properties of their structural parameters obtained through the fitting of their surface brightness distribution.These algorithms assume that the luminosity distribution of a galaxy is the sum of the distributions of its individual components.The surface brightness of each component is modelled by a parametric law, which has to be strictly considered as an empirical fitting function (see the review of these algorithms in Section 3).Nonparametric galaxy classifications are based on the measurements of a set of galaxy parameters that are correlated with the Hubble sequence.The advantage of this method is that it assumes no analytic models and can classify regular and irregular galaxies.Several galaxy parameters have been used for decades for galaxy classification.The colour of galaxies was one of the first parameters used in such classifications since colours are easy to obtain and correlate with galaxy morphological types.Thus, E and S0 galaxies show redder colours than spirals.This correlation is due to colours of galaxies reflecting the dominant constituent stellar population and this correlates with morphology.de Vaucouleurs [31] established the dependence of galaxy colour and morphological type.Colours have also recently been used for the morphological classification of large galaxy databases (e.g., [32] and references therein).Spectra of galaxies also reflect their stellar populations and have been used for galaxy classification purposes.This property has been used in several spectroscopic galaxy classifications performed using certain features or full galaxy spectra in different wavelength and redshift ranges (e.g., [33][34][35][36][37][38][39][40][41][42][43]).All these classifications show correlations between morphological and spectral classes.However, these correlations show wide scatter, particularly for E and S0 galaxies (see Figure 4 and [44] for the quantification of the scatter between morphological and spectral galaxy classes).

Advances in Astronomy
Galaxies can also be classified on the basis of their light distribution.Morgan [45,46] used the central light concentration as a measure of the light distribution of galaxies.Some years later, Fraser [47] defined light concentration indices as the ratio of the radii containing certain light percentages.These indices became popular and were used in several galaxy classifications (see [48][49][50]).Abraham et al. [51,52] introduced the asymmetry of the galaxy light as another useful parameter in galaxy classification.They showed that different types of galaxies (early-type, latetype, and irregular) were located at different positions in the concentration-asymmetry (C-A) plane.Thus, early-type galaxies show a highly concentrated and symmetric light distribution.In contrast, the light distribution of irregular galaxies is less concentrated and homogeneous.A third parameter, smoothness (S), was proposed by Conselice et al. [53], given the CAS morphological classification system.Several other parameters such as the Gini coefficient (the Gini coefficient is a measures of the inequality of a distribution.This coefficient is a number between 0 for perfect equality and 1 for maximal inequality), galaxy ellipticity, or the M20 moment (the M20 moment is the second-order moment of the 20% brightest pixels of a galaxy.)have been included in nonparametric classification over the last decade [54][55][56].All these parameters give some information about the galaxy shape.
It is very difficult to classify galaxies with classical approaches using more than three of the previous parameters.Several classification methods take into account the information contained in more than three parameters, but the final calibrations are done in three dimensions (e.g., [57]).This problem was recently solved by Huertas-Company et al. [58], who developed a generalization of the nonparametric classification that uses an unlimited number of dimensions.The method is based on a particular class of learning machines called support vector machines.The algorithm finds the optimal decision regions in a multidimensional space using a set of visually classified galaxies as a training The sample of galaxies correspond to those galaxies classified visually by Nair and Abraham [28].See more details in [44].
sample (see [59][60][61]).This classification algorithm has recently been applied to ∼700 000 galaxies from the SDSS-DR7 spectroscopic galaxy catalogue [62].All the galaxies were classified into four morphological classes: E, S0, Sab, and Scd.The novelty of this classification is that they assign to each galaxy a probability of being in the four morphological classes instead of assigning a single class.This Bayesian approach is better adapted to nature, where a continuous transition between different morphological types is observed.Good agreement was found between this automatic classification and visually based ones.Thus, good correlations were observed between the probability values given by this automatic classification and visual galaxy classes.For example, galaxies visually classified as ellipticals have on average a probability of ∼0. 8 of being E and ∼0.2 of being S0.The other two probabilities of these galaxies turned to be almost zero (see [62]).This algorithm is one of the most appropriate for galaxy classifications of large galaxy samples due to the reasonable time required and the high level of agreement with other detailed visual classifications.
In the near future, all these new large galaxy samples will provide new insights in understanding the nature and evolution of lenticular and other morphological types of galaxies.

Quantitative Morphology of Lenticular Galaxies
From a structural point of view, galaxies can be simple systems formed by one main structural component, such as elliptical galaxies, or they can be more complex systems, such as spiral galaxies, showing structures such as spheroids, discs, bars, spiral arms, rings, and lenses.The luminosity distribution of a galaxy is the sum of the distributions of its individual components.The surface brightness distribution of the galaxies then provides information about their different structural components.It is assumed that each galaxy component can be modelled by a parametric law.Thus, the fit of the surface brightness distribution of galaxies by one or more components provides their structural parameters (size, scale, shape, and luminosity).Comparison of these structural parameters with different types of galaxies or galaxies of the same type but at different redshifts provides useful information for classifying galaxies and understanding their formation and evolution.One of the basic problems of this quantitative morphology is prior knowledge, gained from solely photometrical information, of the number of structural components of the galaxies.This is crucial because the values of the resulting structural parameters of the galaxies depend on the number of structures used in the fit of their surface brightness distributions (e.g., [63]).This problem has special relevance in the case of S0 galaxies as they are objects between ellipticals and spirals.Usually, studies of large E and S0 datasets fit their surface brightness distributions by using only one component, obtaining global parameters such as scales, effective surface brightness, and the total luminosity of the galaxies.Nevertheless, S0 galaxies have at least two structural components (bulge and disc).Thus, their surface brightness should be fitted by using multicomponent models.The only way to solve this problem is through prior knowledge of the different structural components of the galaxies.This can be done by detailed studies of the photometric information of each galaxy before and/or after performing the decomposition of their surface brightness distribution.After the Advances in Astronomy decompositions it should also be considered whether the parameters obtained are physically meaningful (see [64][65][66]).
Several multicomponent fitting algorithms have been developed for fitting the surface brightness distribution of galaxies (e.g., [64,[67][68][69][70][71][72][73]).Prieto et al. [64] developed one of these multicomponent techniques for fitting the surface brightness profiles of disc galaxies, including lenticulars.The process begins by analysing colour and colour-index images, which provides a qualitative idea of the numbers of components present in the galaxies.The scale lengths of these components were estimated taking into account the colourindex profiles of the galaxies.Confirmation of the assumed components was then obtained in the ellipticity and position angle radial profiles of the isophotes of the galaxies where the geometry of the different components, projected onto the sky plane, are well reflected.In this study and others using multicomponents fits, the surface brightness distribution of lenticular galaxies was fitted with a bulge, disc, bar, and lens.Below, I review the most popular mathematical laws used for fitting the surface brightness distribution of S0 galaxies.Figure 5 shows an example of 2D decompositions of the barred galaxy NGC 4643 using the multicomponent fitting method developed by [73].

Fitting Function of the Structural Components.
Several mathematical functions have been used in order to fit the distribution of the stellar surface brightness of the different structural components of galaxies.These fitting functions include those for bulges, discs, lenses, and bars.

Bulge Fitting Functions.
Historically, the first surface brightness models were proposed for fitting spheroidal galaxy components.The reason for this was that the first galaxies analysed were ellipticals.We can mention the surface-density models proposed by Reynolds [74], Hubble [75], King [76], Rood et al. [77], and Oemler [78].For a good recent review of the mathematical laws proposed for fitting the surface brightness of galaxies see Graham [79].
One of the most popular fitting functions for describing the surface brightness of Es and bulges of S0s and spirals has been the de Vaucouleurs r 1/4 profile [80,81].This model was proposed by de Vaucouleurs, who argued that it fitted the surface brightness of galaxies better than previous models for E galaxies.Indeed, since the work by de Vaucouleurs [82] for the bulge of M31, the r 1/4 law has been used for fitting the surface brightness of bulges of S0 and spiral galaxies.The radial variation of the surface intensity (I dV (r)) of this model is given by I dV (r) = I dV ,e 10 −7.669((r/re) 1/4 −1) , where I dV ,e is the effective surface intensity of the galaxy and r e the effective radius, defined as the radius enclosing half of the total galaxy light.This model has been extensively used in surface brightness decompositions of ellipticals and bulges of S0 and spiral galaxies (e.g., [83][84][85][86]).More recent studies have found deviations from the r 1/4 law of the surface brightness of the spheroidal components of individual galaxies (see [67,[87][88][89]), including the bulge of the Milky Way [90].The universality of the r 1/4 profile for spheroidals began to be questioned after the work by Andredakis and Sanders [91], who found that bulges of latetype galaxies were better modelled by the Sérsic [92] law.This work was extended for a larger sample of galaxies of all morphological types, revealing that their bulges were accurately fitted by Sérsic profiles [93].This was also confirmed by other researchers.Today, the Sérsic law is extensively used for modelling the surface brightness distribution of bulges in galaxies in the optical and near-infrared (e.g., [21,63,64,[94][95][96][97][98][99][100][101][102][103][104][105]).The radial surface intensity (I S (r)) of the Sérsic law is given by where I S,e is the intensity at the effective radius r e that encloses half of the total light of the model.The parameter n describes the shape of the profile, and the coefficient b n is given by ≈ 1.9992n − 0.327 [106].Notice that this profile is a generalization of the r 1/4 law.Thus, for n = 1 the Sérsic profile is exponential, while n = 4 is the de Vaucouleurs profile.Figure 6 shows the Sérsic radial profiles with different values of the n parameter.For useful Sérsic reviews see Ciotti [107], Ciotti and Bertin [108], Simonneau and Prada [109], and Graham and Driver [110].

The Disc Component.
Freeman [111] found that the surface brightness distribution of discs (I D (r)) in S0 and spiral galaxies could be modelled with an exponential function given by where I D,0 is the central disc intensity and h is the scale length of the profile.This exponential profile has been extensively used in surface brightness decompositions of S0 and spiral galaxies.Although some studies have used different fitting functions for the discs (e.g., [84]), they were not as popular as the exponential one.
Recently, detailed analyses of the external regions of disc galaxies have shown that not all galaxy discs show pure exponential surface brightness radial profiles.Most galaxies show truncations (downbending) or antitruncations (upbending) in their outermost regions.Thus, about 60% of galaxies have a break followed by downbending.Moreover, another 30% of galaxies have breaks followed by upbending.The shape of the surface brightness profiles of discs correlates with the Hubble type.Thus, downbending breaks are more frequent in late-type galaxies.In contrast, upbending breaks are more frequent in early-type galaxies (see [112]).

Bar Surface Brightness
Laws.Bars are elliptical structures located in discs of galaxies.It is well known since pioneering morphological classifications of galaxies that strong bars are common in discs of S0 and spiral galaxies (see e.g., [12]).Early studies showed that the surface brightness profiles of bars in early-type galaxies are constant within the bar region while late-type galaxies show exponential-like bar profiles (see [113]).Nevertheless, it is still unclear which is  the best mathematical fitting function describing the surface brightness radial profile of bars.
There are not many studies of the surface-brightness decomposition of barred galaxies.The large number of free parameters needed in these kind of fits makes them difficult.Fortunately, the number of studies including bars in the photometrical decomposition is increasing in the last years.Thus, some studies have used the Sérsic functions for fitting the observed surface brightness profile of bars (see [114]).Nevertheless, the most popular functions for fitting the observed surface brightness profiles of bars are Ferrers (see [65,70]), Freeman (see [63,115]), and flat (see [63,64,115,116]) profiles.
The radial surface intensity profile of a Ferrers ellipsoid (I Ferrers (r)) is given by where I Fe,0 , r bar , and n bar are the central surface intensity, length and a shape parameter of the bar, respectively.The total luminosity for a perfect ellipse is given by In the case of n bar = 2, the previous integral has a primitive, and the total luminosity (L Ferrers ) is given by The radial surface intensity profile (I Freeman (r)) of a Freeman bar is given by I Fr,0 and r bar are the central surface intensity and length of the bar, respectively (see [117]).The Freeman bar is a particular case of the Ferrers bar with n bar = 0.In this case, the total luminosity for a Freeman bar is given by The radial surface intensity profile of a flat bar (I Flat (r)) is where I Fl,0 and r bar are the central surface intensity and the length of the bar.For radii greater than r bar the profile falls with a scale length r s [116].In this case, the total luminosity of this bar is given by where Li 2 is the dilogarithm function (also called the Jonquiere function).
The surface radial profiles of Ferrers, Freeman and flat bars with the same bar radius can be seen in Figure 6.It is not clear which of the previous functions better fit the observed surface brightness distribution of bars.Probably, the type of bar is correlated with the morphological type of the galaxy (see [118]).Nevertheless, more investigation is needed in this subject.
3.1.4.Lens Component.Lenses are photometric components commonly observed in lenticular and some early-type galaxies.In contrast, few late-type galaxies show lenses.Kormendy [119] found that 54% of SB0-SBa galaxies have lenses.A similar fraction of lenses has been observed in nonbarred S0 galaxies [65,66,70].Their surface brightness radial intensity (I l (r)) is characterized by smooth luminosity gradients with a very sharp cut-off.Duval and Athanassoula [120] found a lens in the luminosity profile of the galaxy NGC 5383.This lens was fitted by the expression where I l,0 is the surface intensity of the lens at the centre and r l is the size of the lens.This expression was also used by Prieto et al. [64] and Aguerri et al. [63] for fitting the surface brightness radial profiles for some lenticular galaxies.In multicomponent fits, lenses represent a small percent (less than 10%) of the total light of the galaxies (see [63,64]).The origin and evolution of lenses are unclear.Nevertheless, Kormendy [119] has argued that lenses in lenticulars originate as bars.He suggested that bars evolve into lenses through some internal process.

Properties of the Structural Components of Lenticular Galaxies
The relations between the structural parameters have been extensively studied in the literature in order to understand the formation and evolution processes of different galaxy components (for a recent review see [79]).I review below some of the structural relations given in the literature for samples of galaxies including S0s.In this paper, I took into account galaxy samples with multicomponent surface brightness decompositions.Thus, the structural parameters of all galaxies have been obtained in a homogeneous way.

Bulge Structural Parameters
4.1.1.Correlations with Hubble Type.Andredakis et al. [93], using an extension of the model-independent decomposition method introduced by Kent [121], found correlations of certain structural parameters of bulges with galaxy Hubble type.In particular, using near-infrared images, they found that the bulges of early-type spiral galaxies had a Sérsic shape parameter of n ≈ 4, while bulges from late-type galaxies showed n ≈ 1.This continuous relation indicated a wide diversity concerning the shape parameter in galaxy bulges, similar to ellipticals (e.g., [122][123][124]).This diversity also indicated that the different types of bulges were not well fitted using a single law (such as the de Vaucouleurs law).This trend was later confirmed by other works (e.g., [94,97,[125][126][127][128][129]).
Recent surface brightness decomposition of large galaxy samples, including S0s, show that the correlation between n and Hubble type is not so evident.For example, Mendez-Abreu et al. [103] found no correlations between the bulge parameters (r e , n, and μ e ) and Hubble type for a sample of 148 unbarred S0-Sb bulges (see Figure 7).Laurikainen et al. [21] found that n was almost constant from S0 − to Sa bulges of galaxies and was only slightly larger for S0-Sa galaxies than for late-type bulges.Similar results have been found in other samples with few S0 galaxies (see [63,130]).All these studies indicate that there is a wide scatter of n values for S0 galaxies.Thus, some bulges of S0 galaxies show n values similar to those of Es (n ≈ 4).But there is also a significant fraction of bulges of S0 galaxies showing similar values of the shape parameter (n < 2) to those of late-type (see, e.g., [63,65,73,131] and Figure 7).

The Fundamental Planes and Their
Projections.Dynamically hot systems are not scattered randomly in the parameter space defined by their effective surface brightness (μ e ), half-light radius (r e ), and central velocity dispersion (σ 0 ).On the contrary, they follow a tight correlation called the fundamental plane (FP; [132,133]).The bulges of S0 galaxies are located in the same FPs as E galaxies [63,103] or with a slight offset [134].No statistically significant difference has been observed between the FP defined by bulges of S0s or early-type spirals [103].These studies have also shown that the FP of bulges of S0 galaxies are independent of their n Sérsic values.Thus, S0 bulges with n ≈ 1 and n ≈ 4 are located in the same FP (see [63]).
Kormendy [84] found a strong correlation between the surface central brightness of E galaxies and their effective radius.This relation is also kept when the mean effective surface brightness ( μ e ) is used.This is the so-called Kormendy relation (KR) and is one of the projections of the FP.Several authors have shown that bulges of S0 galaxies follow a similar KR as E galaxies.In contrast, bulges of latetype galaxies show a large offset with respect to the KR of Es [85,88,93,101,135].For a large sample of S0 galaxies, Laurikainen et al. [21] found that the slopes of the KR of bulges of S0 and E galaxies were similar.But, the zero points of the relations were not the same.They reported a small offset of about 0.3 mag.This offset was not observed when bulges brighter than M K = −22. 7were considered (see Figure 8).Similarly, Barway et al. [136] also found that bright S0 galaxies (M K < −24.5) followed the same KR as Es in the Coma cluster, while the KR of faint S0s showed an offset with respect to E galaxies.
Another projection of the FP is the so-called Faber-Jackson (FJ) relation ([see 144]), which relates the luminosity (L) of spheroidal systems to their central velocity dispersion (σ 0 ).Several studies have found that this relation depends strongly on the luminosity of the galaxies (see [137] and references therein).Thus, Matković and Guzmán [138] found that bright elliptical galaxies show a tight FJ relation given by L ∝ σ 4 0 , while faint early-type galaxies follow a significantly shallower relation (L ∝ σ 2.01 0 ).Bulges of S0 and Sa galaxies follow a similar FJ relation to that of bright Es (L ∝ σ 4.2 0 ; [103]).The deviation from the elliptical FJ relation is due to bulges fainter than M B = −20 (see [103,130]).
4. 1.3.Other Relations between Bulge Effective Radius, Luminosity, and Sérsic Shape Parameter.The Sérsic shape index n has a correlation with the total luminosity of bulges (M bulge ) for spiral galaxies (see [21,96,97,123,139]). Thus, brighter bulges show larger values of n.This correlation is also shown for E galaxies (e.g., [124]).In contrast, bulges of S0 galaxies do not show a statistically significant relation between n and M bulge (see [21]).
The effective radius (r e ) of bulges of early-type spirals and S0 correlates with the bulge luminosity (see [65,130]).Nevertheless, bulges of spirals and lenticulars show a different r e -M bulge relation from that for elliptical galaxies (see [140]).In particular, for a given bulge luminosity, bulges of S0 galaxies are more compact than those Es with the same luminosity.This difference was argued by Laurikainen et al. [21], who concluded that bulges of S0s galaxies were photometrically different from those of Es (see Figure 9).Classical bulges can be distinguished from pseudobulges (Many bulges show disk-like features that do not resemble hot stellar systems.These kind of bulges are usually called pseudobulges (see [141]).) by their scaling relations (e.g., [141,142]).In a series of papers, Fisher and Drory [131,143] have shown that classical bulges (defined as those with n > 2) follow the same relations as E galaxies between the effective radius, luminosity and Sérsic shape parameter.In contrast, these relations are not followed for pseudobulges (those with n < 2).They conclude that the Sérsic shape index can be used in order to distinguish between classical bulges and pseudobulges (see [141] for a definition of pseudobulge).This classification is independent of the Hubble type.Thus, their galaxy sample contains lenticular bulges, most of them being of the classical type.But a few of their lenticular bulges are catalogued as pseudobulges.This indicates the wide variety of bulges of lenticular galaxies.

The Intrinsic Shape of Bulges.
In comparison with the large number of works investigating the correlations between the bulge parameters of S0 and spiral galaxies, little attention has been given to the three-dimensional shape of bulges (e.g., [144][145][146][147][148]).Independently of their internal structure, the bulges of lenticulars and spiral galaxies have been traditionally considered as ellipsoidal stellar systems located at the centre of galaxies.Bulges and discs are concentric with the galactic centre and share the same polar axis.In this framework, the twisting of bulge isophotes [144,149] and misalignment between the major axes of the bulge and disk [145,147] are not possible if the bulge and disk are both oblate.They were interpreted as a signature of bulge triaxiality.This idea has been supported by the presence of noncircular gas motions (e.g., [150][151][152][153]) and a velocity gradient along the galaxy minor axis [154][155][156].
In a series of papers, Méndez-Abreu et al. [103,104] have studied the 3D shape of the bulges of a large sample of S0 and   early-type spiral galaxies.They found that 80% of bulges in unbarred lenticulars and early-type galaxies were not oblate but triaxial ellipsoids.Their mean axial ratio in the equatorial plane was B/A = 0.85.This was in agreement with early findings by Bertola et al. [145] and Fathi and Peletier [157].
In addition, Méndez-Abreu et al. [103] found no significant dependence on the probability distribution function of B/A and other galaxy properties, such as morphological type, light concentration, or luminosity of bulges.Méndez-Abreu et al. [104] have found that most of the bulges of S0 and early-type spirals are flattened along their polar axis.Only 18% of the bulges had a probability larger than 50% of being elongated along their polar axis.In this study, the triaxiality parameter (T = (1 − (B/A) 2 )/(1 − (C/A) 2 ), where A, B, and C were the lengths of the longest, intermediate and shortest semiaxes of the bulges; see [158]), was also analysed with the result that the triaxiality distribution of the bulges was strongly bimodal.This bimodality was driven by the Sérsic shape parameter and the B/T ratio of the bulges.Thus, bulges with n < 2 and B/T < 0.3 exhibit a larger fraction of oblate axisymmetric bulges, a smaller fraction of triaxial bulges, and fewer prolate axisymmetric bulges with respect to bulges with n > 2 and B/T > 0.3 (see [104] and Figure 10).No correlation was found between the intrinsic shape of bulges and other galaxy properties, such as the total luminosity or central velocity dispersion.

Correlations with Hubble
Type.Several works found no correlation between the disc scale length and the Hubble type over a wide range of morphological types, including S0 discs (see [94,97,98,103,111,129,130]).Thus, the scale length of the discs was independent of the Hubble type of the galaxy.Freeman [111] discovered that the central surface brightness of the discs in spiral and S0 galaxies has a constant value of (μ 0 = 21.65 mag/arcsec 2 ) with only 0.3 mag/arcsec 2 of standard deviation.Nevertheless, this was not confirmed in later studies.Thus, late-type spirals show fainter central surface brightness than early-type ones (see shows the T distributions for bulges of galaxies with B/T < 0.3 (dotted line) and B/T > 0.3 (dashed line).See [104] for more details.[94,96,129,130]).The correlations of the disc structural parameters for a large sample of bulges from S0 and earlyype galaxies are shown in Figure 11.

4.2.2.
Correlations between h, μ 0 , and Disc Luminosity.The scale length of discs of spiral and lenticular galaxies correlates with the disc luminosity.Thus, larger discs are located in brighter ones (see [21,103,130]).The discs of S0 and spiral galaxies also show a correlation between μ 0 and h in the sense that larger discs have a fainter central surface brightness [21,63,103,130].In both relations the discs of spiral galaxies and lenticulars are located in the same regions, indicating that they are similar (see [21] and Figure 11).

The Tully-Fisher
Relation.The Tully-Fisher relation (TFR; [159]) is a correlation between the asymptotic circular velocity (V c ) of disc galaxies and their integrated stellar luminosity (L).It is well known that the scatter and the zero point of this relation depend on the Hubble type.Thus, late-type spirals show a tight TFR with low dispersion (σ ≈ 0.3).In contrast, Sa and Sb galaxies show TFRs with wider scatter and different zero points from those of Sc ones (e.g., [160,161]).The TFR of S0 galaxies has been less investigated in the literature owing to the difficulty in obtaining the asymptotic circular velocity of these gas-poor galaxies.Nevertheless, several studies indicate that the TFR of S0 galaxies shows a similar slope to that of late-type ones but with a different zero point (see Figure 12).This offset indicates that, for a given circular velocity, S0 galaxies are fainter than late-type ones (see [162][163][164][165][166][167]).The offset is the expected one if the progenitors of present S0 galaxies were late-type galaxies that have stopped their star formation a few Gyr ago.

The Rotational
Velocity-Size Plane.The distribution of disc galaxies in the rotational velocity-size (V c -R) plane provides information about galaxy formation.These two parameters are related with the total angular momentum (J) and mass (M) of galaxies, which have been repeatedly invoked as the fundamental parameters in galaxy formation.The dotted line indicates the TF for late-type spirals and the full line shows the TF of S0 galaxies from the Coma cluster.Credit: Hinz et al. [162].Reproduced by permission from the AAS.
The location of galaxies in the V c -R plane was analysed by Prieto et al. [168] and Campos-Aguilar et al. [169] for a sample of spiral and S0 galaxies.They found that for a given mass Sc galaxies showed larger angular momentum than Sa galaxies.In contrast, lenticular and irregular galaxies were located at the same position as spiral galaxies in the V-R plane.This result indicates that, dynamically, S0 galaxies are no different from spirals of any type.This could point toward the old idea from van den Bergh [20] that lenticulars and spirals form parallel sequences as shown lately by different groups of researchers [18,21,22].

The Prominence of the Bulge in the Hubble Sequence.
Traditionally, the bulge-to-disc-luminosity ratio (B/D) and the bulge-to-total-luminosity ratio (B/T) of galaxies have been used for morphological galaxy classifications.The early work by Burstein [88] shows that for S0 galaxies B/D > 0.3 (or B/T > 0.75).Some years later, Simien and de Vaucouleurs [86] found that the Hubble classification sequence was mainly driven by the decreasing B/D ratio.In this picture, S0 galaxies were located between the Es and spirals, as Hubble claimed.Several researchers have used the previous results and have classified their galaxies according to the values of the B/D or B/T parameters.Thus, S0 galaxies have been considered as those disc galaxies dominated by the bulge component (B/T > 0.6-0.7).
Recent surface brightness decomposition of samples of galaxies including S0s challenge the previous results.In particular, Balcells et al. [130], using a small galaxy sample (17) including only 4 S0s, found no correlation between B/D and the Hubble type but that S0 galaxies has a mean value of B/D ≈ 0.25.The large B/T dispersion of S0 galaxies have also been observed in other larger samples.Thus, Laurikainen et al. [65] and Weinzirl et al. [73] show that B/T increases from late to early-type galaxies, but that, nevertheless, S0/a galaxies have similar B/T values (≈ 0.2-0.3) to Sa ones.A similar result was obtained by Laurikainen et al. [21], who shows that S0 galaxies have a wide spread of B/T values.Thus, 13% of their S0 galaxies have B/T values similar to those for late-type galaxies (B/T < 0.15).All these new results show that some galaxies visually classified as S0 have small bulges, as do late-type ones.These galaxies were classified as S0s due to their smooth extended disc, in spite of their bulge prominence.Notice also that galaxy classifications using B/D or B/T could have been misclassified a significant fraction of S0 galaxies.
The large spread found in the B/T values of S0 galaxies could also be in agreement with other results suggesting that spiral and S0 galaxies formed parallel sequences within which early-and late-type galaxies are distinguished by their B/T values (see e.g., [18,[20][21][22].

Correlation between r e and h.
The correlation between the bulge effective radius and disc scale length has been extensively studied in the literature (see [66,94,95,97,103,126]).These studies show that larger bulges reside in larger discs.This correlation is shown in galaxy samples including (or not including) S0 galaxies.It has been argued that this relation indicates a link between bulge and disc formation.Barway et al. [170] found that for S0 galaxies the correlation between r e and h depends on the galaxy luminosity.Thus, faint lenticulars show positive correlation between r e and h while this correlation was not present for bright S0s.Nevertheless, this dichotomy was not confirmed by other researches in larger galaxy samples (see [66]).
The relation between the ratio of the bulge effective radius to disk scale length (r e /h) and the Hubble type has also been studied.Several investigations have found a lack of correlation between those two quantities [94,95,130], indicating that the Hubble sequence is scale-free.Thus, Graham and Worley [140] and McArthur et al. [126] found a mean value r e /h ≈ 0.2 for spirals.In the case of S0 galaxies, the mean values of r e /h reported in the literature goes from 0.2 (see [21]) to 0.36 (see [103]).

Bar Parameters.
The presence of bars in disc galaxies has been extensively studied in the literature.Bars have been identified through visual inspection (e.g., [12,28,171,172]), by analysing the shape and orientation of the galaxy isophotes [118,[173][174][175][176], by studying the Fourier modes of the light distribution [113,[177][178][179]) and by decomposing the surface brightness of the galaxies into their structural components including bars [21,63,64,70,73,114,116,179,180].These studies indicate that bars are observed in optical images of roughly half of all the nearby disc galaxies.This fraction rises to about 70% when near-infrared images are analysed.
The fraction of barred galaxies is not constant with the Hubble type.Neither is there agreement on whether the bar fraction is larger in S0 galaxies or in late-type ones.Thus, Aguerri et al. [118] found a bar in 29% of lenticular galaxies and ∼55% in early-and late-type spirals.In addition, the bar fraction in galaxy discs is higher in those galaxies with bulges with n < 2 and B/T < 0.2 (see [73]), suggesting that bars avoid galaxies with large and prominent bulges (see [118]; but see also [171]).This is consistent with the findings by Laurikainen et al. [65] who found that B/T in barred galaxies was smaller than in nonbarred ones.A smaller fraction of bars in lenticulars than in spirals was also found by Laurikainen et al. [66] for a low redshift sample and by Nair and Abraham [28] for a sample of galaxies located at z < 0.1.Similarly, Barway et al. [181] also found a small bar fraction in lenticulars (15%), this fraction being a function of galaxy luminosity.Thus, less luminous galaxies have a higher fraction of bars in their sample.Giordano et al. [182] also observed a smaller fraction of bars in S0s than in spirals but, in contrast, less massive galaxies show smaller bar fractions.This was also found by Méndez-Abreu et al. [104] for a sample of bars in the Coma cluster spanning a large mass range.Larger fraction of bars in S0 galaxies, similar to those reported for late-type ones, has been found in other studies (see e.g., [183]).
There are three main observational parameters of bars: length, strength, and pattern speed.I review below the most important results concerning these three parameters for S0 galaxies.

The Bar Length.
The bar length has been obtained by eye estimates from galaxy images [119,184], by locating the maximum ellipticity of the galaxy isophotes [173,175,176], by looking for variations in the isophotal position angle [118,185,186] or in the phase angle of the Fourier modes of galaxy light distributions [115,187], or by photometric decomposition of surface brightness distributions [21,63,64,70,73,114,116].These studies have reported a typical bar length of 3-4 kpc correlated with the disc scale length or galaxy size [63,173,180].The bar length also depends on other galaxy properties such as colour or prominence of the bulge [188].
The bar length normalized by galaxy size depends on the Hubble type.Thus, lenticular galaxies show larger bars than late-type ones [113,118,186,189].Laurikainen et al. [65] also found that the ratio between the bar radius and the disc scale depends on the Hubble type, being larger for S0 galaxies than for late-types.

The Bar Strength.
The bar strength represents the contribution of the bar to the total potential of the galaxy.Several methods have been developed to measure it.The most commonly used is the Q g parameter.This parameter estimates the bar torque and was suggested by Combes and Sanders [191].It has been calculated using different methods (see, [187,[192][193][194]).This parameter can be accurately estimated by using near-infrared galaxies (see [65,195,196]).Other bar strength measurements are related to the bar ellipticity.This method assumes that bar strength is a function of their ellipticity.Thus, more elliptical bars are stronger than rounder ones (see [197]).
The two methods give similar strength measurements of the bars due to the fact that bar strength computed with bar ellipticities correlates with Q g (e.g., [198]).Note that bar ellipticity and Q g do not give similar strengths as a function of the Hubble type (e.g., [65]).Nevertheless, Laurikainen et al. [65] showed that S0 bars are weaker than late-type ones using the Q g parameter or the bar ellipticity for measuring the bar strength.In this sense, some studies have shown that bars in lenticular galaxies are weaker than in late-type spirals measuring the bar strength by the Q g parameter (see [65,190] and Figure 13) or by the bar ellipticity [118].

The Bar Pattern Speed.
Bar pattern speed is the main kinematic observable in barred galaxies and measures the tumbling pattern of the bar.The pattern speed is related to the so-called corotation radius of the galaxy.This is the radius at which the angular velocity of the disc is equal to the pattern speed of the bar.Barred galaxies are classified according to the distance-independent ratio R = D L /a B , where D L is the corotation radius of the bar and a B is the bar semi-axis length.Self-consistent bars cannot exist with R < 1 (e.g., [199]).Bars showing 1.0 < R < 1.4 are called fast bars and those with R > 1.4 are slow bars.
Several methods have been developed for measuring the pattern speed of bars, including hydrodynamic simulations [201][202][203][204][205][206][207][208] and the shape of dust lanes in the bar regions of some barred galaxies [209,210].However, the best way to measure the bar pattern speed of galaxies is through the model-independent method developed by Tremaine and Weinberg ( [211]; hereafter the TW method).
The TW method has already been applied to a large sample of S0 barred galaxies [115,154,[212][213][214][215][216].All these measurements show that bars of S0 galaxies are fast.According to numerical simulations of Debattista and Sellwood [217,218], strong fast bars are located in galaxies with maximal discs.This means that inside the bar region the baryonic matter dominates.Otherwise, the bar dynamical friction with the dark matter halo of the galaxy slows down the bar efficiently and would produce a slow bar.Other researches claim that the bar pattern speed does not only depend on the mass concentration (see, e.g., [219]).Figure 14 shows the corotation radius versus the bar radius for all galaxies with pattern speed measured using the TW method.This figure shows that, within the error bars, all bars are compatible with being fast (i.e., R CR /a bar < 1.4).
The dependence of the bar pattern speed on the Hubble type is not well understood yet.Few measurements of bar pattern speeds for late-type galaxies are available in the literature.Nevertheless, there are some hints that S0 and latetype galaxies show different pattern speed.Thus, Aguerri et al. [178] found a weak correlation between the parameter R and the Hubble type for a sample of ten barred galaxies.They showed that lenticular galaxies have smaller values of R than late-type ones.Recently, Rautiainen et al. [208] have found that all the early-type galaxies of their sample show fast bars,  d) of bar strengths in S0 (dotted line) and spiral galaxies (full line).This figure shows that bars are weaker in S0 than in spiral galaxies.Credit: Buta et al. [190].Reproduced by permission from the AAS.whereas bars of the late-type galaxies included both fast and slow rotators (see Figure 15).
It is important to notice that all these results concerning bars indicate that discs of late-type and lenticular galaxies are different.Thus, lenticulars show a smaller fraction of bars than late-type ones.In addition, bars in lenticular galaxies are larger, weaker, and faster than late-type bars.

The Case of dS0 Galaxies
Dwarf galaxies are the most numerous type of galaxies in the universe.They are defined as those with low central surface brightness and luminosities below M B = −18.According to their colours we can distinguish two main classes of dwarf galaxies: dwarf ellipticals (dEs) and dwarf irregulars (dIs).In general, their surface-brightness profiles can be well fitted by the Sérsic profiles with small shape parameter values n ≈ 1-2 [102,124,220].Nevertheless, there are some dE galaxies with surface brightness profiles that are not well fitted by a single component.Sandage and Binggeli [221] noticed this in some dE galaxies in the Virgo Cluster.Later, Binggeli and Cameron [122] fitted an exponential profile to the outermost regions of the surface brightness of these galaxies, but this profile clearly cannot fit their surface brightness at all radii.These results indicate that these objects have two structural components.These types of dwarf galaxies were called dwarf lenticulars (dS0s).These galaxies are also gas-poor systems like dEs and have evolved with old stellar components but Figure 14: The corotation radius (R CR ) versus the bar radius (a bar ) for barred galaxies from the literature with bar pattern speed measured by applying the TW method.Full and open stars refer to barred and unbarred lenticular galaxies, respectively.Open circles show spiral barred galaxies.The full, dashed, and dotted lines represent R CR /a bar = 1.0, 1.4, and 2.0, respectively.Fast bars are usually defined as those with R CR /a bar < 1.4 (see text for more details).The data for this figure was taken from [200].they have higher flattening [222].Some of these galaxies in the Virgo cluster have disc-like structures due to the presence of bars and/or spiral arms patterns [223][224][225][226][227][228].
In the Coma cluster, 56% of the dwarf galaxies in the magnitude interval −16 < M B < −18 show surface brightness radial profiles well fitted by only one structural component.In contrast, the surface brightness radial profiles of 29% of the dwarfs in the same magnitude interval cannot be fitted by a single Sérsic component (see [102]).
The structural relations of dwarf galaxies in the Coma cluster were analysed by Aguerri et al. [102], who found that the dS0 population in the Coma cluster show similar structural parameters to those of the dS0s located in the Virgo cluster.In addition, the structural parameters of the disc components in the dS0 galaxies showed similar relations to the discs of bright spirals.The scale length of the bulges of dS0 galaxies was similar to the scale length of bulges of latetype galaxies.Nevertheless, the mean scale of dEs was found to be higher than that of the bulges of all morphological types of bright galaxies.These results indicate that dS0 galaxies are similar to late-type spirals.Thus, dS0s seem to be low-luminosity late-type spiral galaxies (see Figure 16).The dichotomy in the formation between elliptical and spheroidal galaxies has been recently pointed out by other authors (see e.g., [229]).

Formation of Bright and Dwarf Lenticular Galaxies: Different Formation Scenarios
In the previous sections I have reviewed the properties of the structural parameters of lenticular and spiral galaxies.
It should be noticed that the scaling relations shown before give a wide variety of results.Thus, some results indicate that S0 bulges are similar to E galaxies (e.g., same FP, FJ, and size-mass relations, large n Sérsic values), while other scaling relations show that S0 galaxies are similar to late-type ones (e.g., small values of the n parameter, similar slope of the TF relation but different zero point).This could reinforce the idea of S. van den Bergh that the class of S0 galaxies is formed by galaxies with different origins and evolutions that only share a similar morphological shape.In this section I review the different formation scenarios proposed for lenticular galaxies: I indicate which S0s could be formed by each of these scenarios according to their properties.

Major Merger Formation.
According to the current paradigm of galaxy formation, spheroidal systems (ellipticals and bulges of galaxies) have been formed by major or minor galaxy mergers.In contrast, discs do not survive after major mergers.They were formed by cooling of gas in rotating dark matter haloes.In this paradigm the central bulge of disc galaxies forms prior to the disc as a result of early merging (see the reviews by [141,230]).In contrast to the bulge, the thin disc forms later from surrounding gas (e.g., [231]).Galaxy major mergers produce remnants with similar shape and scaling relations to those observed in bulges of S0 galaxies.In particular, the gas content of the progenitors determines the shape and structural parameters of the remnants.Thus, dissipationless numerical simulations of major mergers produce triaxial remnants with a tendency to become prolate systems.In contrast, dissipational remnants are weakly triaxial, but close to oblate [232][233][234][235][236]. Dissipationless numerical simulations of bulge + disc galaxy mergers with different mass ratios produce disc galaxies with B/T > 0.2 and bulges with a Sérsic shape parameter mainly in the interval 3 < n < 4 [237].In addition, Hopkins et al. [238] has found that high B/T disc galaxies tend to form in major mergers, whereas low B/T objects tend to form from minor mergers.The remnants produced by galaxy major mergers also follow typical scaling relations to those of elliptical galaxies (see [239]).The major merger scenario could be applicable to some of the S0 galaxies, especially the most luminous and with prominent bulges (n > 2 and B/T > 0.3).The scaling relations previously shown indicate that the bulges of these S0 galaxies follow the typical scaling relations of elliptical galaxies (fundamental plane, Kormendy relation, size-luminosity relation).In addition, the observed bimodal distribution of the triaxiality parameter shown in a large sample of bulges of S0-Sa galaxies has been compared with the predictions of numerical simulations of spheroid formation (see [104]).In this context, the variety of intrinsic shapes of bulges of S0-Sa galaxies with n > 2 and B/T > 0.3 required both major dissipational and dissipationless mergers [104].
Recently, it was discovered that a small fraction of nearby E/S0 galaxies show blue colours.These galaxies, called blue ellipticals, constitute only about 10% of the total population of nearby E/S0 galaxies (e.g., [61,[240][241][242]).This fraction changes with redshift, constituting up to 30% at z ≈ 1 (see [61]).Blue ellipticals with small masses could be the results of merger remnants that triggered the central star formation and built a central bulge component.The gas surrounding this bulge could build a disc on a longer time scale (see [61]).These small blue ellipticals observed at high redshift could be the progenitors of some nearby S0 galaxies.
Major mergers are rare in rich galaxy clusters.They are more frequent in field or galaxy groups.Thus, S0 galaxies formed by this mechanism should be located preferentially in less dense galaxy environments.

Formation of S0 by Satellite Accretion.
Satellite accretion must have occurred several times in a disc galaxy over a Hubble time.Mergers of galaxies with similar mass fully destroy the disc and leave remnants similar to E galaxies.In contrast, minor mergers damage the disc less and can produce bulge growth.The growth of the central bulge of the main galaxy depends on the density and relative mass of the accreted satellite.The remnants of these kind of mergers retain significant amounts of rotation (see [244,245]).Nevertheless, the resulting disc is heated by the Figure 17: The time evolution of the luminous matter of a minor merger (see details in [243]).
tidal interactions between the satellite and the main galaxy ( [246][247][248]).The final disc is similar to S0 galaxies ( [249]).Figure 17 shows the time evolution of the luminous matter of simulated minor mergers.It can be seen that after the merger the main galaxy retains a thick disc (see [243]).The transformation of late-type galaxies by minor mergers was analysed in a series of papers by Aguerri et al. [243] and Eliche-Moral et al. [250].
The most important galaxy transformations were produced by the accretion of dense satellites (see [243]).Dissipationless simulations of minor mergers with different mass ratios and dense satellites show that the accreted satellite can reach the galaxy centre and drive significant growth of the central bulge.The initial conditions of the main galaxy in these simulations were bulge + disc galaxies with an exponential bulge, similar to observed bulges of  [243].Each arrow starts at the location of the original model and ends at the n and B/D derived from the twocomponent fit of the surface mass density profile of the models.The points represent the distribution of n versus log(B/D) for bulges of galaxies from [93].See more details is [243].late-type galaxies.After the accretion of these types of satellites the Sérsic shape parameter can rise from 1 to 4 dependent on the mass merger rate.Moreover, the B/D ratio of the final galaxy was in all cases higher than 0.43, similar to the B/D ratios observed in some S0 and Sa galaxies.These simulations show that the efficiency of the process suggests that present day bulges of late-type spirals showing exponential profiles cannot have grown significantly by collisionless mergers.Figure 18 shows the evolution of n and B/D due to the accretion of dense satellites of different masses (see [243]).These simulations also show that S0 galaxies with bulges with large n Sérsic values n > 2-3 and large B/D ratios can be formed by the transformation of latetype galaxies through minor mergers (see [243]).
The accretion of less dense satellites was studied by Eliche-Moral et al. [250].The primary and secondary galaxies were spiral galaxies with their particles distributed in bulge, disc, and halo components.Both galaxies follow the Tully-Fisher relation similar to real ones.This scaling means that the material of the satellites did not reach the galaxy centre but was instead distributed through the disc of the main galaxy.Thus, particles from the satellite disc deposit their mass onto the primary disc.In contrast, particles from the satellite bulge settle into a rather cold inner disc or ring (see also Eliche-Moral et al. [251]).In addition, an inward flow of material from the disc to the bulge region is produced during the satellite accretion due to the presence of transitory nonasymmetric distortions in the main disc.The satellite accretion drives the initial main galaxy towards an earliertype secular evolution.In particular, the primary galaxy evolves toward higher B/D, n, central velocity dispersion, scale length of the disc, and lower central surface brightness of the disc.In particular, n increases by a factor of ≈ 2.3 and B/D doubles its initial value (see Figure 19).
In both cases the accretion of high-and low-density satellites hits the thin disc of the primary galaxy.This vertical heating of the disc depends on the orbit and mass of the satellite.The thick disc created by these types of mergers is reminiscent of the thick discs seen in S0 and early-type galaxies.Bekki [252] and Bournaud et al. [253] also showed that unequal-mass gas-rich spirals (mass ratio ≈ 0.3) can also produce remnants similar to S0 galaxies.The mergers produce an enhanced star formation, thicken the disc of the main galaxy, and produce a subsequent transformation of the main galaxy into an S0 type galaxy.This formation mechanism has been proposed in order to explain the formation of some S0 galaxies [254,255].
The minor accretion models proposed by Aguerri et al. [243] and Eliche-Moral et al. [250] have certain limitations and caveats that should be discussed.Neither model includes gas hydrodynamics or star formation.We may expect some changes in the dynamics of the remnants due to gas concentration in the galaxy nucleus.This could produce changes in the M/L of the central regions of the galaxies and changes in some structural relations.It is also important to notice that dense satellites easily produce counterrotating bulges in the primary galaxy when the satellite has retrograde orbits.This is an important limitation because counterrotating bulges are not very common in galaxies (see [256]).This is not the case in the simulations with less dense satellites [250].Nevertheless, in this case the Sérsic shape parameter only grows from n = 1 to n ≈ 2-2.5.Thus, lenticular galaxies with n ≈ 4 and large B/D cannot be explained.

Evolution of Galaxies in High-Density Environments.
Since the first observations of galaxies in clusters it has been well known that the galaxies located in high-density environments have different properties from those located in the field [257][258][259].Early studies claimed that galaxies in clusters and the field were different because they were born differently.Today, it is thought that these differences could be due to the evolution of galaxies being dependent on the environment where galaxies live (e.g., [260] But, see also [18]).
One of the differences between galaxies in clusters and the field is morphology.In dense environments the fraction of ellipticals and S0 galaxies is larger than in the field.In contrast, the fraction of spiral and irregular galaxies is higher in less dense environments.This is the so-called morphologydensity relation observed in nearby [261,262] and high redshift [263] clusters.Figure 20 shows the morphologydensity relation for a sample of nearby galaxy clusters (see [262]).This morphological mixing is not constant with redshift.Thus, it was observed that the S0 population of galaxies is smaller in more distant clusters.This decrement is compensated by an increase in the spiral population (see [263,264]).The morphology-density relation has been one of the strongest arguments in favour of the existence of an evolutionary trend between spiral and lenticular galaxies.In  this picture, S0 galaxies are descendent from spiral galaxies that have stopped their star formation due to environmental processes.Galaxies in high-density environments evolve through many different physical processes.We may mention harassment [260], ram pressure stripping [265][266][267], tidal effects and mergers [268,269], or starvation [267,270].I review below the morphological galaxy transformations produced by these environmental processes.6.3.1.Gas Stripping.X-ray observations of galaxy clusters reveal the presence of a large amount of hot gas (T ≈ 10 7 -10 8 K) in the intracluster medium.This hot gas contains most of the baryons of galaxy clusters.The motion of the galaxies through this hot intracluster medium produces a ram pressure that can strip a significant fraction of the gas content of galaxies in clusters (see [265,266]).This ram pressure is proportional to ρ ICM v 2 gal , where ρ ICM and v 2 gal are the gas density of the intracluster medium (ICM) and the  relative velocity between the galaxy and the ICM.This effect is thus stronger in the central regions of massive clusters where the gas density and the relative velocity of the galaxies are higher.The efficiency of the stripping depends on several properties of clusters and galaxies, such as the location of the stripped gas in the galaxy, the masses of the host cluster and galaxies, and the density of the gas.The gas in galaxies can be found as cold or hot gas.The cold gas is located inside galaxy discs, while the hot gas is located in the galaxy halo.Both types of gas can be stripped from galaxies.The gas from galaxy haloes can be stripped more easily than gas from the discs.The required relative velocity between the galaxy and cluster and the gas density of the intracluster gas are thus significantly much lower than those from disc stripping (see [267]).The gas stripping also depends on the mass of the host cluster and the mass of the galaxy.Stripping is more efficient in more massive clusters and in less luminous galaxies [267,271,272].The stripping also depends on the density of the gas inside galaxies.Thus, low-density gas is more quickly stripped at all galaxy radii than high-density gas.
Gas stripping has been observed in several bright dwarf galaxies located in nearby galaxy clusters (e.g., [273][274][275][276][277][278][279]).It is also supported by the HI deficiency observed in some spiral galaxies in clusters [280][281][282][283][284]. Both disc and halo gas stripping produce similar effects.This is the truncation of the star formation and reddening of galaxy discs producing strong morphological transformations.Thus, late-type spirals can be transformed into lenticular ones (see [270] and Figure 21).The gas stripping and the stopping of the gas accretion onto the disc of the galaxies also stop the formation of bars in galaxies.In addition, it is also expected that gas stripping does not produce significant disk heating.This mechanism could explain several observational properties of lenticular galaxies such as the observed low values of n and B/T for some lenticular galaxies, the lower fraction of bars in S0 galaxies, the different Tully-Fisher relation shown by S0 galaxies, and the specific number of globular clusters observed in S0 galaxies (e.g., [21,22,118,167,285]).

Galaxy Interactions.
Galaxies in high-density environments like galaxy clusters can suffer several gravitational interactions with a wide variety of impact parameters, relative velocity between galaxies and mass ratios.Interactions in rich galaxy clusters are expected to be more abundant and faster than in small galaxy groups.The harassment scenario proposes that the cumulative effects of fast tidal encounters between galaxies and with the gravitational potential can produce dramatic morphological transformations in galaxies.Thus, bright galaxies can lose a significant fraction of their stellar and gas content.The fraction of mass lost can be high enough in order to transform bright galaxies into dwarfs ( [260,[286][287][288][289][290]).
Observational indications of this mass loss have been found for the bright spiral galaxies located in the Coma cluster.Thus, Aguerri et al. [101] and Gutierrez et al. [100] found that the scale lengths of bright spiral galaxies in the Coma cluster were significantly smaller than those for spirals with similar luminosities located in less dense environments.The disrupted stars could form part of the so-called diffuse light detected in some nearby galaxy clusters such as Virgo [291][292][293][294][295][296], Coma [297,298], and Fornax [299].
Mastropietro et al. [288] and Aguerri and González-Garcia [289] have studied the properties of the remnants after several fast tidal encounters.The initial conditions of the simulations done by Aguerri and González-Garca [289] consist of disc galaxies with different B/T ratios.These initial systems undergo several tidal interactions, transforming the initial spiral galaxy into a dwarf-like system.The simulations show that fast tidal interactions are very efficient mechanisms for removing stars and dark matter particles from galaxies, especially from the disc and halo components.In contrast, the bulge particles were not removed.After four or five tidal interactions, the initial galaxy has lost up to 80% of its total mass (see Figure 22).The kinematic analysis of the remnants shows that some of them still conserve an external rotating disc.In contrast, other remnants are systems supported by pressure.Nevertheless, the scale length of the resulting disks is about 40-50% shorter than the initial ones.In addition, the tidal interaction easily creates tidal features in the discs such as stable bars.
The final position of the resulting galaxies in the fundamental plane depends on the initial conditions of the simulations.Thus, models with initially large B/D ratios are closer to the face-on view of the fundamental plane defined Figure 22: Change of total (a) and luminous (b) mass of simulated galaxies after different fast tidal events (see details in [289]).
by bright Es and bulges of early-type galaxies.In contrast, galaxies with initially small B/D ratios are located at the position of the dwarf galaxies in the face-on view of the FP.These simulations show that dE and dS0 galaxies can be formed by the tidal disruption of late-type spirals.In contrast, the tidal disruption of early-type galaxies would produce a final galaxy similar to elliptical or S0 galaxies.One of the most prominent features that tidal interactions induce in the discs of galaxies is the presence of stable and prominent bars (see [288,289]).Although the properties of these bars are similar to those observed in some dwarf galaxies, bars are not common in dwarf galaxies (see [172]).This bar formation is one of the caveats of these simulations and could be related to their initial conditions.Other initial conditions of these models should be taken into account, in particular initial models with thicker discs, which could reduce the bar formation (see [300]).
Aguerri and González-Garcia [289] have analysed the correlations between certain structural parameters of the remnants.In particular, they analysed the correlation between Sérsic parameters r e /h and n observed in spiral galaxies (e.g., [94,103,125,126,128,129]). Thus, those bulges from spiral galaxies with larger values of r e /h also show higher n values.The simulated dwarf objects formed by fast tidal interactions do not follow the trend followed by real galaxies.Thus, for a fixed value of n, tidally formed dwarf galaxies have larger r e /h ratios.This is because discs were strongly truncated due to tidal encounters.In contrast, the effective radius of the bulges does not change.The search for outliers in the r e /h versus n relation could provide dwarf systems (dS0s) formed through the fast interactions of disc galaxies.
It is well known that early-type galaxies follow the luminosity-metallicity and velocity dispersion-metallicity relations (e.g., [301]).Thus, more metal-rich galaxies have higher central velocity dispersions and are more luminous.Tidal stripping should not change the metallicity of the galaxies.This would imply that low-mass systems formed by tidal stripping of bright disc galaxies could well deviate from the luminosity-metallicity relations displayed by bright galaxies.This point is still not well analysed in the literature.Nevertheless, some approaches have recently been undertaken.For example, it has recently been observed that some dwarf galaxies are located outside to the metallicityluminosity relation (e.g., [302]).

Cluster Tides.
Galaxies in high-density environments do not only interact with one another but also interact with the cluster gravitational potential.Thus, cluster and/or group tides can also produce strong morphological transformations in spiral galaxies.The tidal field of the clusters as a whole can trigger cluster spiral, nuclear, and disc activities.This effect is more important near the centre of clusters or groups.A higher fraction of starburst galaxies might be expected in the cores of clusters or groups of galaxies.The long-term effects of this activity would be to clear the gas from disc spirals and changing them into S0 galaxies [303,304].Semianalytic models suggest that tidal fields can produce a segregation of galaxies by mass in clusters [305,306].The evolution of the galaxies inside clusters is determined by their trajectories and due to the substructure presented in galaxy clusters.These tidal interactions are insufficient to transform a spiral galaxy into an elliptical one.Nevertheless, they can produce the transformation of spiral galaxies into S0s [307].
Evidence for the morphological transformation of galaxies by the cluster environment has been observed in the S0 population of galaxies in the Coma cluster.Poggianti et al. [308] studied the stellar populations of elliptical and S0 galaxies in the Coma cluster.They found that a large fraction of the S0 population has undergone star formation during the last 5 Gyr, not observed in elliptical galaxies.They concluded that the S0 galaxies with recent star formation are the descendants of star-forming spirals whose star formation has been halted due to the effect of the cluster environment.Nevertheless, the Coma cluster could be an exception to the fact that most of the morphological transformations in clusters since z = 1 is not confined to massive clusters.On the contrary, it was most pronounced in low-mass clusters (see [309]).

Environmental Tides.
The picture of galaxies falling into static galaxy clusters is very simplistic and far removed from reality.In the hierarchical structure formation theory dark matter haloes grow through the accretion of smaller units.This formation is scale-free.Thus, clusters and groups of galaxies grow through the accretion of smaller groups.Optical and X-ray studies of galaxy clusters show that a large fraction of them have recently accreted other groups or clusters [310][311][312][313][314].This can be seen in the substructure detected in some nearby galaxy clusters (see [315] and references therein).
Mergers of galaxy clusters could enhance star formation, as suggested by numerical simulations (e.g., [316]) and is claimed to have been seen in some observational studies [317,318].These mergers therefore produce strong gravitational tides acting across galaxies and transforming spiral galaxies into S0s [307,319].Nevertheless, it is not easy to observe a cluster merger.The transformation of galaxies due to cluster mergers could depend, among other considerations, on the phase in which the merger is observed, or on the viewing angle.There are very few well-identified examples of cluster mergers in the literature.They have been identified by a clear segregation between the dark matter and the intracluster gas.These mergers are the bullet cluster [320,321], Abell 520 [322], and MACS J0025.4-1225[323].This last cluster has recently been studied by Ma et al. [324], who found a large fraction of S0 galaxies and interpreted this result as a transformation of spiral galaxies into S0s due to ram pressure and/or tidal forces during the cluster merger.6.4.Passive Evolution of Spiral Galaxies.The distribution of galaxies in colour-magnitude and colour-mass diagrams shows that they are located in two main regions: the blue cloud and the red sequence.Study of the colour-magnitude diagram of galaxies at different redshifts indicates that galaxies stop their star formation and pass from the blue cloud to the red sequence.Moreover, this migration is strongly mass dependent.Thus, more massive galaxies reach the red sequence at higher redshift than less massive ones (e.g., [325,326]).The fact that galaxies are not homogeneously distributed throughout the colour-magnitude diagram indicates that galaxies pass from the blue cloud to the red sequence on short time scales.
The reddening of spiral galaxies could produce S0s by fading the spiral arms of the progenitors.This reddening could be due to processes related to the galaxy environment, or they could be processes related only to the internal galaxy evolution.One of the internal processes that could regulate the star formation in galaxies is feedback processes heating the cold gas in galaxies and stopping their star formation.Dekel and Silk [327] proposed that feedback from supernovae could stop the star formation in galaxies and drive their evolution.This feedback scenario could drive the star formation history and evolution of dwarf galaxies, but it cannot heat the gas of massive haloes and stop their star formation.Such additional heating could come from feedback of active galactic nuclei (AGN; see [328]).Thus, outflows and jets from the central supermassive black holes of galaxies are expected to heat the cold gas and stop further star formation in the galaxies (e.g., [329][330][331]).In a series of papers, Schawinski et al. [332,333] proposed that AGN feedback is mainly responsible for stopping the star formation in star-forming early-type galaxies.This process would transform early-type spiral galaxies located in the blue cloud into galaxies located close to or in the red sequence.These transformed early-type galaxies could be the progenitors of later S0 galaxies.This mechanism would produce S0 galaxies with similar properties to those of earlytype spiral galaxies.
The red sequence is not only formed by ellipticals and early-type galaxies.It has recently been discovered that a small fraction of late-type spirals (about 6%) show red colour and are located in the red sequence (see [334]).The red colours of these late-type galaxies is not due to large dust extinction.In addition, these galaxies are mainly located in intermediate galaxy density environments.This implies that the environment alone has not produced the reddening of these galaxies.Masters et al. [334] have also demonstrated that these red spiral galaxies are not recent poststarbursts.Rather, their star formation has ceased gradually.They found that most red spirals (about 70%) exhibit bars.They argued that the cessation of the star formation in these galaxies should be related to the presence of bars.These passive spirals could be precursors of some S0 galaxies (see also [335]).
These two internal processes (AGN feedback and passive evolution) could explain the similarities of S0s and spiral galaxies.The S0 galaxies formed by these mechanisms should be in intermediate or less dense galaxy environments.However, these internal processes could not explain the different bar properties observed between spiral and lenticular galaxies: different bar fractions (e.g, [66,118,182]) or pattern speed (e.g., [178,194,208]).

Conclusions
In this paper I have reviewed the observational properties of the photometric components of lenticular galaxies.I have selected samples of galaxies containing S0s and with structural parameters obtained in a homogeneous way by multicomponent surface brightness decompositions.The properties of the different components show that some lenticular galaxies are similar to ellipticals whereas others show similar relations to spiral galaxies.In particular, the most luminous S0s with B/T > 0.3 and n > 2 follow similar relations to elliptical galaxies.In contrast, less luminous lenticulars with B/T < 0.3 and n < 2 have structural relations closer to spirals.These properties could indicate that the class of lenticular galaxies could be formed by galaxies with different origins and evolutions with only the shape as the common link.
I have reviewed the different formation mechanisms proposed in the literature for understanding the formation of lenticular galaxies.There is a large variety of formation mechanisms that produce S0s over very different time scales.Each of these mechanisms produces S0 galaxies with different properties.Thus, lenticular galaxies with high n and B/T values should be produced by one major or several minor mergers or passive evolution of early-type galaxies.Those S0 galaxies with low n and B/T values should be produced by transformation of late-type galaxies within high-density environments or by passive evolution of late-type spirals.
The morphological classifications of large galaxy samples that have recently appeared in the literature will in the near future provide important clues to understanding the formation and evolution of galaxies, including lenticulars.These classifications provide a unique opportunity for selecting large samples of lenticular galaxies located in a large diversity of environments in a homogeneous way, from field to galaxy clusters.These large samples will permit us to solve the biases of present small galaxy samples.Detailed analysis of the structural parameters, together with the analysis of the environments where lenticulars live, will in the near future provide new and important constrains to the formation and evolution of lenticular galaxies.

Figure 1 :
Figure 1: The Hubble tuning fork diagram.A schematic view of the Hubble galaxy classification.Credit: NASA and ESA.

Figure 3 :
Figure3: Galaxy classification scheme proposed by the ATLAS 3D project.Galaxies are classified as fast and slow rotators.Fast rotators form a parallel sequence to spiral galaxies.Credit: Cappellari et al.[18].Reproduced by permission from WILEY.

Figure 4 :
Figure4: Scatter plot of Hubble versus spectral types.The sample of galaxies correspond to those galaxies classified visually by Nair and Abraham[28].See more details in[44].

Figure 5 :
Figure 5: (left panels) r-band images the NGC4643.Sérsic (middle upper panel), Sérsic plus disc (middle middle panel), and Sérsic plus disc plus bar (middle-lower panel) models of the galaxy.(right panels) the corresponding residuals.Credit: Weinzirl et al. [73].Reproduced by permission from the AAS.

Figure 7 :
Figure 7: Correlations between the bulge parameters for a sample of lenticulars and early-type spirals.Correlation between the Hubble type and effective radius (a), effective surface brightness (b), and Sérsic parameter (c).Correlations between the effective radius of the bulge and the Sérsic parameter (d), effective surface brightness (e), and absolute magnitude (f).The full line represents the best linear fit.The Pearson correlation coefficient (r) is also indicated.See [103] for more details.

Figure 8 :
Figure8: The Kormendy relation for bulges of S0, early-type and late-type spirals.In all panels the full line represents the KR of elliptical galaxies in the Coma cluster.The dotted black line corresponds to the extension of the KR of ellipticals towards galaxies with smaller effective radius.The red dashed line of the upper panel represents the KR for the S0 bulges.The red and green dashed lines in the lower panel show the KR for S0 bulges brighter and fainter than M k = −22.7,respectively.Credit: Laurikainen et al.[21].Reproduced by permision from WILEY.

Figure 9 :Figure 10 :
Figure 9: Effective radius (a) and effective surface brightness (b) of bulges of S0 (red points), early-type (crosses), and late-type (squares) spirals as a function of the bulge absolute magnitude.The dashed line in the left panel corresponds to the relation followed by ellipticals galaxies from the Coma cluster.The continuous lines in the left panel correspond to the mean values of the data points in one-magnitude bins: S0s (red line) and Sa-Sbc (blue line).Credit: Laurikainen et al. [21].Reproduced by permision from WILEY.

Figure 11 :
Figure 11: Correlations between the disc parameters for a sample of lenticulars and early-type spirals.Correlation between the Hubble type and disc scale length (a) and the central surface brightness (b).Correlations between the disc scale length and the absolute magnitude of the disc (c) and the central surface brightness (d).The full line represents the best linear fit (y = a + bx).The Pearson correlation coefficient (r) is also indicated.See [103] for more details.

Figure 12 :
Figure 12: Tully-Fisher diagram for S0 Coma galaxies in I-band.The dotted line indicates the TF for late-type spirals and the full line shows the TF of S0 galaxies from the Coma cluster.Credit: Hinz et al.[162].Reproduced by permission from the AAS.

Figure 13 :
Figure13: Histograms of the distribution (a) and (c) and cumulative distribution (b) and (d) of bar strengths in S0 (dotted line) and spiral galaxies (full line).This figure shows that bars are weaker in S0 than in spiral galaxies.Credit: Buta et al.[190].Reproduced by permission from the AAS.

5 RFigure 15 :
Figure 15: The R = R CR /a bar parameter as a function of the Hubble type.The horizontal line (R = 1.4) represents the separation between fast and slow bars.Credit: Rautiainen et al. [208].Reproduced by permision from WILEY.

Figure 16 :
Figure 16: Effective radius (a), Sérsic shape parameter (b), and central surface brightness (c) as functions of the R-band absolute magnitude of dE (asterisks) and bulges of dS0 (diamonds), late-type spirals (circles), early-type spirals (crosses), and S0 galaxies (squares) in the Coma cluster.See for more details [102].

Figure 18 :
Figure18: Growth vectors of the bulges in the n-log(B/D) plane for the accretion models by[243].Each arrow starts at the location of the original model and ends at the n and B/D derived from the twocomponent fit of the surface mass density profile of the models.The points represent the distribution of n versus log(B/D) for bulges of galaxies from[93].See more details is[243].

Figure 19 :
Figure19: The dependence of bulge and disc parameters, and B/D ratio, on the increment of the bulge magnitude (see details in[250]).

Figure 20 :
Figure 20: Galaxy population gradients in galaxy clusters as a function of radial distance from the cluster center (a) and local galaxy density (b) Credit: Dressler [262].Reproduced by permission from the AAS.

Figure 21 :
Figure 21: Morphological properties for different models from [270] with no gas accretion at T = 0.4 Gyr (a) and T = 3.4 Gyr (b) Models with high (c) and low (d) gas accretion rates.Credit: Bekki et al. [270].Reproduced by permission from the AAS.