This is the first paper in a series devoted to review the main properties of galaxies designated S0 in the Hubble classification system. Our aim is to gather abundant and, above all, robust information on the most relevant physical parameters of this poorly-understood morphological type and their possible dependence on the environment that could later be used to assess their possible formation channel(s). The adopted approach combines the characterisation of the fundamental features of the optical spectra of $68{,}043$ S0 with heliocentric $z\lesssim 0.1$ with the exploration of a comprehensive set of their global attributes. A principal component analysis is used to reduce the huge number of dimensions of the spectral data to a low-dimensional space facilitating a bias-free machine-learning-based classification of the galaxies. This procedure has revealed that objects bearing the S0 designation consist, despite their similar morphology, of two separate sub-populations with statistically inconsistent physical properties. Compared to the absorption-dominated S0, those with significant nebular emission are, on average, somewhat less massive, more luminous with less concentrated light profiles, have a younger, bluer and metal-poorer stellar component, and avoid high-galaxy-density regions. Noteworthy is the fact that the majority of members of this latter class, which accounts for at least a quarter of the local S0 population, show star formation rates and spectral characteristics entirely similar to those seen in late spirals. Our findings suggest that star-forming S0 might be less rare than hitherto believed and raise the interesting possibility of identifying them with plausible progenitors of their quiescent counterparts.

We present an automated statistical method that uses medium-resolution spectroscopic observations of a set of stars to select those that show evidence of possessing significant amounts of neutron-capture elements. Our tool was tested against a sample of $\sim 70,000$ F- and G-type stars distributed among $215$ plates from the Galactic Understanding and Exploration (SEGUE) survey, including $13$ that were directed at stellar Galaxy clusters. Focusing on five spectral lines of europium in the visible window, our procedure ranked the stars by their likelihood of having enhanced content of this atomic species and identifies the objects that exhibit signs of being rich in neutron-capture elements as those scoring in the upper $2.5\%$. We find that several of the cluster plates contain relatively large numbers of stars with significant absorption around at least three of the five selected lines. The most prominent is the globular cluster M3, where we measured a fraction of stars that are potentially rich in heavy nuclides, representing at least $15\%$.

A suite of 432 collisionless simulations of bound pairs of spiral galaxies with mass ratios 1:1 and 3:1, and global properties consistent with the $\Lambda$CDM paradigm, is used to test the conjecture that major mergers fuel the dual AGN (DAGN) of the local volume. Our analysis is based on the premise that the essential aspects of this scenario can be captured by replacing the physics of the central BH with restrictions on their relative separation in phase space. We introduce several estimates of the DAGN fraction and infer predictions for the activity levels and resolution limits usually involved in surveys of these systems, assessing their dependence on the parameters controlling the length of both mergers and nuclear activity. Given a set of constraints, we find that the values adopted for some of the latter factors often condition the outcomes from individual experiments. Still, the results do not reveal, in general, very tight correlations, being the tendency of the frequencies normalized to the merger time to anticorrelate with the orbital circularity the clearest effect. In agreement with other theoretical studies, our simulations predict intrinsic abundances of these systems that range from $\sim$few to $15\%$ depending on the maximum level of nuclear activity achieved. At the same time, we show that these probabilities are reduced by about an order of magnitude when they are filtered with the typical constraints applied by observational studies of the DAGN fraction at low redshift. As a whole, the results of the present work prove that the consideration of the most common limitations involved in the detection of close active pairs at optical wavelengths is sufficient by itself to reconcile the intrinsic frequencies envisaged in a hierarchical universe with the small fractions of double-peaked narrow-line systems which are often reported at kpc-scales.