The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) will soon start to scan thousands of square degrees of the northern extragalactic sky with a unique set of $56$ optical filters from a dedicated $2.55$m telescope, JST, at the Javalambre Astrophysical Observatory. Before the arrival of the final instrument (a 1.2 Gpixels, 4.2deg$^2$ field-of-view camera), the JST was equipped with an interim camera (JPAS-Pathfinder), composed of one CCD with a 0.3deg$^2$ field-of-view and resolution of 0.23 arcsec pixel$^{-1}$. To demonstrate the scientific potential of J-PAS, with the JPAS-Pathfinder camera we carried out a survey on the AEGIS field (along the Extended Groth Strip), dubbed miniJPAS. We observed a total of $\sim 1$ deg$^2$, with the $56$ J-PAS filters, which include $54$ narrow band (NB, $\rm{FWHM} \sim 145$Angstrom) and two broader filters extending to the UV and the near-infrared, complemented by the $u,g,r,i$ SDSS broad band (BB) filters. In this paper we present the miniJPAS data set, the details of the catalogues and data access, and illustrate the scientific potential of our multi-band data. The data surpass the target depths originally planned for J-PAS, reaching $\rm{mag}_{\rm {AB}}$ between $\sim 22$ and $23.5$ for the NB filters and up to $24$ for the BB filters ($5\sigma$ in a $3$~arcsec aperture). The miniJPAS primary catalogue contains more than $64,000$ sources extracted in the $r$ detection band with forced photometry in all other bands. We estimate the catalogue to be complete up to $r=23.6$ for point-like sources and up to $r=22.7$ for extended sources. Photometric redshifts reach subpercent precision for all sources up to $r=22.5$, and a precision of $\sim 0.3$% for about half of the sample. (Abridged)

We present the first systematic study of the stellar populations of ultra-diffuse galaxies (UDGs) in the field, integrating the large area search and characterization of UDGs by the SMUDGes survey with the twelve-band optical photometry of the S-PLUS survey. Based on Bayesian modeling of the optical colors of UDGs, we determine the ages, metallicities and stellar masses of 100 UDGs distributed in an area of $\sim 330$ deg$^2$ in the Stripe 82 region. We find that the stellar masses and metallicities of field UDGs are similar to those observed in clusters and follow the trends previously defined in studies of dwarf and giant galaxies. However, field UDGs have younger luminosity-weighted ages than do UDGs in clusters. We interpret this result to mean that field UDGs have more extended star formation histories, including some that continue to form stars at low levels to the present time. Finally, we examine stellar population scaling relations that show that UDGs are, as a population, similar to other low-surface brightness galaxies.

From the approximately $\sim$3,500 planetary nebulae (PNe) discovered in our Galaxy, only 14 are known to be members of the Galactic halo. Nevertheless, a systematic search for halo PNe has never been performed. In this study, we present new photometric diagnostic tools to identify compact PNe in the Galactic halo by making use of the novel 12-filter system projects, J-PLUS (Javalambre Photometric Local Universe Survey) and S-PLUS (Southern-Photometric Local Universe Survey). We reconstructed the IPHAS (Isaac Newton Telescope (INT) Photometric H${\alpha}$ Survey of the Northern Galactic Plane) diagnostic diagram and propose four new ones using i) the J-PLUS and S-PLUS synthetic photometry for a grid of photo-ionisation models of halo PNe, ii) several observed halo PNe, as well as iii) a number of other emission-line objects that resemble PNe. All colour-colour diagnostic diagrams are validated using two known halo PNe observed by J-PLUS during the scientific verification phase and the first data release (DR1) of S-PLUS and the DR1 of J-PLUS. By applying our criteria to the DR1s ($\sim$1,190 deg$^2$), we identified one PN candidate. However, optical follow-up spectroscopy proved it to be a H II region belonging to the UGC 5272 galaxy. Here, we also discuss the PN and two H II galaxies recovered by these selection criteria. Finally, the cross-matching with the most updated PNe catalogue (HASH) helped us to highlight the potential of these surveys, since we recover all the known PNe in the observed area. The tools here proposed to identify PNe and separate them from their emission-line contaminants proved to be very efficient thanks to the combination of many colours, even when applied -like in the present work- to an automatic photometric search that is limited to compact PNe.

Despite numerous efforts, it is still unclear whether lenticular galaxies (S0s) evolve from spirals whose star formation was suppressed, or formed trough mergers or disk instabilities. In this paper we present a pilot study of 21 S0 galaxies in extreme environments (field and cluster), and compare their spatially-resolved kinematics and global stellar populations. Our aim is to identify whether there are different mechanisms that form S0s in different environments. Our results show that the kinematics of S0 galaxies in field and cluster are, indeed, different. Lenticulars in the cluster are more rotationally supported, suggesting that they are formed through processes that involve the rapid consumption or removal of gas (e.g. starvation, ram pressure stripping). In contrast, S0s in the field are more pressure supported, suggesting that minor mergers served mostly to shape their kinematic properties. These results are independent of total mass, luminosity, or disk-to-bulge ratio. On the other hand, the mass-weighted age, metallicity, and star formation time-scale of the galaxies correlate more with mass than with environment, in agreement with known relations from previous work such as the one between mass and metallicity. Overall, our results re-enforce the idea that there are multiple mechanisms that produce S0s, and that both mass $and$ environment play key roles. A larger sample is highly desirable to confirm or refute the results and the interpretation of this pilot study.

We present a star/galaxy classification for the Southern Photometric Local Universe Survey (S-PLUS), based on a Machine Learning approach: the Random Forest algorithm. We train the algorithm using the S-PLUS optical photometry up to $r$=21, matched to SDSS/DR13, and morphological parameters. The metric of importance is defined as the relative decrease of the initial accuracy when all correlations related to a certain feature is vanished. In general, the broad photometric bands presented higher importance when compared to narrow ones. The influence of the morphological parameters has been evaluated training the RF with and without the inclusion of morphological parameters, presenting accuracy values of 95.0\% and 88.1\%, respectively. Particularly, the morphological parameter {\rm FWHM/PSF} performed the highest importance over all features to distinguish between stars and galaxies, indicating that it is crucial to classify objects into stars and galaxies. We investigate the misclassification of stars and galaxies in the broad-band colour-colour diagram $(g-r)$ versus $(r-i)$. The morphology can notably improve the classification of objects at regions in the diagram where the misclassification was relatively high. Consequently, it provides cleaner samples for statistical studies. The expected contamination rate of red galaxies as a function of the redshift is estimated, providing corrections for red galaxy samples. The classification of QSOs as extragalactic objects is slightly better using photometric-only case. An extragalactic point-source catalogue is provided using the classification without any morphology feature (only the SED information) with additional constraints on photometric redshifts and {\rm FWHM/PSF} values.

The Southern Photometric Local Universe Survey (S-PLUS) is imaging ~9300 deg^2 of the celestial sphere in twelve optical bands using a dedicated 0.8 m robotic telescope, the T80-South, at the Cerro Tololo Inter-American Observatory, Chile. The telescope is equipped with a 9.2k by 9.2k e2v detector with 10 um pixels, resulting in a field-of-view of 2 deg^2 with a plate scale of 0.55"/pixel. The survey consists of four main subfields, which include two non-contiguous fields at high Galactic latitudes (8000 deg^2 at |b| > 30 deg) and two areas of the Galactic plane and bulge (for an additional 1300 deg^2). S-PLUS uses the Javalambre 12-band magnitude system, which includes the 5 u, g, r, i, z broad-band filters and 7 narrow-band filters centered on prominent stellar spectral features: the Balmer jump/[OII], Ca H+K, H-delta, G-band, Mg b triplet, H-alpha, and the Ca triplet. S-PLUS delivers accurate photometric redshifts (delta_z/(1+z) = 0.02 or better) for galaxies with r < 20 AB mag and redshift < 0.5, thus producing a 3D map of the local Universe over a volume of more than 1 (Gpc/h)^3. The final S-PLUS catalogue will also enable the study of star formation and stellar populations in and around the Milky Way and nearby galaxies, as well as searches for quasars, variable sources, and low-metallicity stars. In this paper we introduce the main characteristics of the survey, illustrated with science verification data highlighting the unique capabilities of S-PLUS. We also present the first public data release of ~336 deg^2 of the Stripe-82 area, which is available at http://datalab.noao.edu/splus.

In this paper we present a thorough discussion about the photometric redshift (photo-z) performance of the Southern Photometric Local Universe Survey (S-PLUS). This survey combines a 7 narrow + 5 broad passband filter system, with a typical photometric-depth of r$\sim$21 AB. For this exercise, we utilize the Data Release 1 (DR1), corresponding to 336 deg$^{2}$ from the Stripe-82 region. We rely on the \texttt{BPZ2} code to compute our estimates, using a new library of SED models, which includes additional templates for quiescent galaxies. When compared to a spectroscopic redshift control sample of $\sim$100k galaxies, we find a precision of $\sigma_{z}<$0.8\%, $<$2.0\% or $<$3.0\% for galaxies with magnitudes r$<$17, $<$19 and $<$21, respectively. A precision of 0.6\% is attained for galaxies with the highest \texttt{Odds} values. These estimates have a negligible bias and a fraction of catastrophic outliers inferior to 1\%. We identify a redshift window (i.e., 0.26$<z<$0.32) where our estimates double their precision, due to the simultaneous detection of two emission-lines in two distinct narrow-bands; representing a window opportunity to conduct statistical studies such as luminosity functions. We forecast a total of $\sim$2M, $\sim$16M and $\sim$32M galaxies in the S-PLUS survey with a photo-z precision of $\sigma_{z}<$1.0\%, $<$2.0\% and $<$2.5\% after observing 8000 $deg^{2}$. We also derive redshift Probability Density Functions, proving their reliability encoding redshift uncertainties and their potential recovering the $n(z)$ of galaxies at $z<0.4$, with an unprecedented precision for a photometric survey in the southern hemisphere.

Galaxies within galaxy clusters are known to be subject to a wide variety of environmental effects, both gravitational and hydrodynamical. In this study, we examine the purely gravitational interaction of idealised galaxy models falling into a galaxy cluster in the context of Modified Newtonian Dynamics (MOND). This modification of gravity gives rise to an external field effect (EFE), where the internal dynamics of a system are affected by the presence of external gravitational fields. We examine the consequences of the EFE on low and high mass disk galaxies in time-evolving analytic background cluster potentials, considering orbits with weak and strong tidal fields. By varying the orbital plane of the galaxies we also test the effect of having the tidal interaction orthogonal or parallel to the disk. Furthermore, we consider as a control sample models where the EFE has been removed and they are only affected by tides. Our results suggest that MOND cluster galaxies should exhibit clear asymmetries in their isophotes, suffer increased mass loss and a reduction in their rotation curves due to the combined effect of cluster tides and the external field. In particular, low mass galaxies are hit hard by the EFE, becoming dominated by dispersion rather than rotation even in the absence of tides.

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

Lenticular galaxies are generally thought to have descended from spirals via morphological transformation, although recent numerical simulations have shown that minor or even major merger can also lead to an S0-like remnant. These mechanisms, however, are active in a dense environment such as a group or a cluster of galaxies - making it harder to explain the remarkable fraction of S0s found in the field. Here, we propose a new mechanism to form such lenticular galaxies. We show that an isolated cold disk settled into rotational equilibrium becomes violently unstable - leading to fragmentation and formation of stellar clumps that, in turn, not only grow the bulge, but also increase the stellar disk velocity dispersion optimally in less than a billion year. Subsequently, the galaxy evolves passively without any conspicuous spiral structure. The final galaxy models resemble remarkably well the morphology and stellar kinematics of the present-day S0s observed by the Planetary Nebulae spectrograph. Our findings suggest a natural link between the high-redshift clumpy progenitors to the present-day S0 galaxies.

Recovering the origins of lenticular galaxies can shed light on the understanding of galaxy formation and evolution, since they present properties that can be found in both elliptical and spiral galaxies. In this work we study the kinematics of the globular cluster (GC) systems of three lenticular galaxies located in low density environments (NGC 2768, NGC 3115 and NGC 7457), and compare them with the kinematics of planetary nebulae (PNe). The PNe and GC data come from the Planetary Nebulae Spectrograph and the SLUGGS Surveys. Through photometric spheroid-disc decomposition and PNe kinematics we find the probability for a given GC to belong to either the spheroid or the disc of its host galaxy or be rejected from the model. We find that there is no correlation between the components that the GCs are likely to belong to and their colours. Particularly, for NGC 2768 we find that its red GCs display rotation preferentially at inner radii (Re<1). In the case of the GC system of NGC 3115 we find a group of GCs with similar kinematics that are not likely to belong to either its spheroid nor disc. For NGC 7457 we find that 70% of its GCs are likely to belong to the disc. Overall, our results suggest that these galaxies assembled into S0s through different evolutionary paths. Mergers seem to have been very important for NGC 2768 and NGC 3115 while NGC 7457 is more likely to have experienced secular evolution.

The Javalambre Photometric Local Universe Survey (J-PLUS) provides wide field-of-view images in 12 narrow, intermediate and broad-band filters optimized for stellar photometry. Here we have applied J-PLUS data for the first time for the study of Galactic GCs using science verification data obtained for the very metal-poor GC M\,15. Our J-PLUS data provide low-resolution spectral energy distributions covering the near-UV to the near-IR, allowing us to search for MPs based on pseudo-spectral fitting diagnostics. J-PLUS CMDs are found to be particularly useful to search for splits in the sequences formed by the upper red giant branch (RGB) and asymptotic giant branch (AGB) stars. We interpret these split sequences as evidence for the presence of MPs. This demonstrates that the J-PLUS survey will have sufficient spatial coverage and spectral resolution to perform a large statistical study of GCs through multi-band photometry in the coming years.

We aim to use multi-band imaging from the Phase-3 Verification Data of the J-PLUS survey to derive accurate photometric redshifts (photo-z) and look for potential new members in the surroundings of the nearby galaxy clusters A2589 (z=0.0414) & A2593 (z=0.0440), using redshift probability distribution functions. The ultimate goal is to demonstrate the usefulness of a 12-band filter system in the study of large-scale structure in the local universe. We present an optimized pipeline for the estimation of photo-z in clusters of galaxies. We tested our photo-z with a sample of 296 spectroscopically confirmed cluster members with a magnitude of <r>= 16.6 and redshift <z>=0.041. The combination of seven narrow and five broadband filters with a typical photometric-depth of r<21.5 provides dz/(1+z)=0.01 photo-z estimates. A precision of dz/(1+z)=0.005 is obtained for the 177 galaxies brighter than magnitude r<17. To foresee the precision beyond the spectroscopic sample, we designed a set of simulations in which real cluster galaxies are modeled and reinjected inside the images at different signal-to-noise. A precision of dz/(1+z)=0.02 and dz/(1+z)=0.03 is expected at <r>= 18-22, respectively. Complementarily, we used SDSS/DR12 data to derive photo-z estimates for the same galaxy sample, demonstrating that the wavelength-resolution of the J-PLUS can double the precision achieved by SDSS for galaxies with a high S/N. We find as much as 170 new candidates across the entire field. The spatial distribution of these galaxies may suggest an overlap between the systems with no evidence of a clear filamentary structure connecting the clusters. These preliminary results show the potential of J-PLUS data to revisit membership of groups and clusters from nearby galaxies, important for the determination of luminosity and mass functions and environmental studies at the intermediate and low-mass regime.

The spatial variations of stellar population properties within a galaxy are intimately related to their formation process. Therefore, spatially resolved studies of galaxies are essential to uncover their formation and assembly. The Javalambre Photometric Local Universe Survey (J-PLUS) is a dedicated multi-filter designed to observed ~8500 deg2 using twelve narrow-, intermediate- and broad-band filters in the optical range. In this study, we test the potential of the multi-filter observation carried out with J-PLUS to investigate the properties of spatially-resolved nearby galaxies. We present detailed 2D maps of stellar population properties (age, metallicity, extinction, and stellar mass surface density) for two early-type galaxies observed in both, J-PLUS and CALIFA surveys: NGC 5473 and NGC 5485. Radial structures are also compared and luminosity- and mass-weighted profiles are derived. We use MUFFIT to process the J-PLUS observations, and two different techniques (STARLIGHT and STECKMAP) to analyze IFU CALIFA data. We demonstrate that this novel technique delivers radial stellar population gradients in good agreement with the IFU technique CALIFA/STECKMAP although comparison of the absolute values reveals the existence of intrinsic systematic differences. Radial stellar population gradients differ when CALIFA/STARLIGHT methodology is used. Age and metallicity radial profiles derived from J-PLUS/MUFFIT are very similar when luminosity- or mass-weighted properties are used, suggesting that the contribution of a younger component is small. Comparison between the three methodologies reveals some discrepancies suggesting that the specific characteristics of each method causes important differences. We conclude that the ages, metallicities and extinction derived for individual galaxies not only depend on the chosen models but also depend on the method used.

Our goal is to morphologically classify the sources identified in the images of the J-PLUS early data release (EDR) into compact (stars) or extended (galaxies) using a suited Bayesian classifier. J-PLUS sources exhibit two distinct populations in the r-band magnitude vs. concentration plane, corresponding to compact and extended sources. We modelled the two-population distribution with a skewed Gaussian for compact objects and a log-normal function for the extended ones. The derived model and the number density prior based on J-PLUS EDR data were used to estimate the Bayesian probability of a source to be star or galaxy. This procedure was applied pointing-by-pointing to account for varying observing conditions and sky position. Finally, we combined the morphological information from g, r, and i broad bands in order to improve the classification of low signal-to-noise sources. The derived probabilities are used to compute the pointing-by-pointing number counts of stars and galaxies. The former increases as we approach to the Milky Way disk, and the latter are similar across the probed area. The comparison with SDSS in the common regions is satisfactory up to r ~ 21, with consistent numbers of stars and galaxies, and consistent distributions in concentration and (g - i) colour spaces. We implement a morphological star/galaxy classifier based on PDF analysis, providing meaningful probabilities for J-PLUS sources to one magnitude deeper (r ~ 21) than a classical boolean classification. These probabilities are suited for the statistical study of 150k stars and 101k galaxies with 15 < r < 21 present in the 31.7 deg2 of the J-PLUS EDR. In a future version of the classifier, we will include J-PLUS colour information from twelve photometric bands.

J-PLUS is an ongoing 12-band photometric optical survey, observing thousands of square degrees of the Northern hemisphere from the dedicated JAST/T80 telescope at the Observatorio Astrof\'isico de Javalambre. T80Cam is a 2 sq.deg field-of-view camera mounted on this 83cm-diameter telescope, and is equipped with a unique system of filters spanning the entire optical range. This filter system is a combination of broad, medium and narrow-band filters, optimally designed to extract the rest-frame spectral features (the 3700-4000\AA\ Balmer break region, H$\delta$, Ca H+K, the G-band, the Mgb and Ca triplets) that are key to both characterize stellar types and to deliver a low-resolution photo-spectrum for each pixel of the sky observed. With a typical depth of AB $\sim 21.25$ mag per band, this filter set thus allows for an indiscriminate and accurate characterization of the stellar population in our Galaxy, it provides an unprecedented 2D photo-spectral information for all resolved galaxies in the local universe, as well as accurate photo-z estimates ($\Delta\,z\sim 0.01-0.03$) for moderately bright (up to $r\sim 20$ mag) extragalactic sources. While some narrow band filters are designed for the study of particular emission features ([OII]/$\lambda$3727, H$\alpha$/$\lambda$6563) up to $z < 0.015$, they also provide well-defined windows for the analysis of other emission lines at higher redshifts. As a result, J-PLUS has the potential to contribute to a wide range of fields in Astrophysics, both in the nearby universe (Milky Way, 2D IFU-like studies, stellar populations of nearby and moderate redshift galaxies, clusters of galaxies) and at high redshifts (ELGs at $z\approx 0.77, 2.2$ and $4.4$, QSOs, etc). With this paper, we release $\sim 36$ sq.deg of J-PLUS data, containing about $1.5\times 10^5$ stars and $10^5$ galaxies at $r<21$ mag.

Globular clusters (GCs) can be considered discrete, long-lived, dynamical tracers that retain crucial information about the assembly history of their parent galaxy. In this paper, we present a new catalogue of GC velocities and colours for the lenticular galaxy NGC 1023, we study their kinematics and spatial distribution, in comparison with the underlying stellar kinematics and surface brightness profile, and we test a new method for studying GC properties. Specifically, we decompose the galaxy light into its spheroid (assumed to represent the bulge + halo components) and disk components and use it to assign to each GC a probability of belonging to one of the two components. Then we model the galaxy kinematics, assuming a disk and spheroidal component, using planetary nebulae (PNe) and integrated stellar light. We use this kinematic model and the probability previously obtained from the photometry to recalculate for each GC its likelihood of being associated with the disk, the spheroid, or neither. We find that the reddest GCs are likely to be associated with the disk, as found for faint fuzzies in this same galaxy, suggesting that the disk of this S0 galaxy originated at z ~ 2. The majority of blue GCs are found likely to be associated with the spheroidal (hot) component. The method also allows us to identify objects that are unlikely to be in equilibrium with the system. In NGC1023 some of the rejected GCs form a substructure in phase space that is connected with NGC 1023 companion galaxy.

The Javalambre-Physics of the Accelerated Universe Astrophysical Survey (J-PAS) is a narrow band, very wide field Cosmological Survey to be carried out from the Javalambre Observatory in Spain with a purpose-built, dedicated 2.5m telescope and a 4.7 sq.deg. camera with 1.2Gpix. Starting in late 2015, J-PAS will observe 8500sq.deg. of Northern Sky and measure $0.003(1+z)$ photo-z for $9\times10^7$ LRG and ELG galaxies plus several million QSOs, sampling an effective volume of $\sim 14$ Gpc$^3$ up to $z=1.3$ and becoming the first radial BAO experiment to reach Stage IV. J-PAS will detect $7\times 10^5$ galaxy clusters and groups, setting constrains on Dark Energy which rival those obtained from its BAO measurements. Thanks to the superb characteristics of the site (seeing ~0.7 arcsec), J-PAS is expected to obtain a deep, sub-arcsec image of the Northern sky, which combined with its unique photo-z precision will produce one of the most powerful cosmological lensing surveys before the arrival of Euclid. J-PAS unprecedented spectral time domain information will enable a self-contained SN survey that, without the need for external spectroscopic follow-up, will detect, classify and measure $\sigma_z\sim 0.5\%$ redshifts for $\sim 4000$ SNeIa and $\sim 900$ core-collapse SNe. The key to the J-PAS potential is its innovative approach: a contiguous system of 54 filters with $145\AA$ width, placed $100\AA$ apart over a multi-degree FoV is a powerful "redshift machine", with the survey speed of a 4000 multiplexing low resolution spectrograph, but many times cheaper and much faster to build. The J-PAS camera is equivalent to a 4.7 sq.deg. "IFU" and it will produce a time-resolved, 3D image of the Northern Sky with a very wide range of Astrophysical applications in Galaxy Evolution, the nearby Universe and the study of resolved stellar populations.

Faint fuzzies are metal-rich apparently-old star clusters with unusually large radii (7-15 pc), found mostly in S0 galaxies, whose source remain obscure. To identify their origins, we compare planetary nebulae and neutral hydrogen with faint fuzzy positions and line-of-sight velocities in NGC1023. In this way, we rule out scenarios in which these objects are associated with an on-going merger or with a spheroid population in NGC1023. Their kinematics are indistinguishable from the stellar disk population in this galaxy, and we conclude that faint fuzzies are most likely just remnant open clusters. Their observed association with S0s then simply reflects the difficulty of identifying such objects in later-type disk galaxies.

The stellar kinematics of the spheroids and discs of S0 galaxies contain clues to their formation histories. Unfortunately, it is difficult to disentangle the two components and to recover their stellar kinematics in the faint outer parts of the galaxies using conventional absorption line spectroscopy. This paper therefore presents the stellar kinematics of six S0 galaxies derived from observations of planetary nebulae (PNe), obtained using the Planetary Nebula Spectrograph. To separate the kinematics of the two components, we use a maximum-likelihood method that combines the discrete kinematic data with a photometric component decomposition. The results of this analysis reveal that: the discs of S0 galaxies are rotationally supported; however, the amount of random motion in these discs is systematically higher than in comparable spiral galaxies; and the S0s lie around one magnitude below the Tully--Fisher relation for spiral galaxies, while their spheroids lie nearly one magnitude above the Faber--Jackson relation for ellipticals. All of these findings are consistent with a scenario in which spirals are converted into S0s through a process of mild harassment or "pestering," with their discs somewhat heated and their spheroid somewhat enhanced by the conversion process. In such a scenario, one might expect the properties of S0s to depend on environment. We do not see such an effect in this fairly small sample, although any differences would be diluted by the fact that the current location does not necessarily reflect the environment in which the transformation occurred. Similar observations of larger samples probing a broader range of environments, coupled with more detailed modelling of the transformation process to match the wide range of parameters that we have shown can now be measured, should take us from these first steps to the definitive answer as to how S0 galaxies form.

The origins of S0 galaxies remain obscure, with various mechanisms proposed for their formation, likely depending on environment. These mechanisms would imprint different signatures in the galaxies' stellar kinematics out to large radii, offering a method for distinguishing between them. We aim to study a sample of six S0 galaxies from a range of environments, and use planetary nebulae (PNe) as tracers of their stellar populations out to very large radii, to determine their kinematics in order to understand their origins. Using a special-purpose instrument, the Planetary Nebula Spectrograph, we observe and extract PNe catalogues for these six systems*. We show that the PNe have the same spatial distribution as the starlight, that the numbers of them are consistent with what would be expected in a comparable old stellar population in elliptical galaxies, and that their kinematics join smoothly onto those derived at smaller radii from conventional spectroscopy. The high-quality kinematic observations presented here form an excellent set for studying the detailed kinematics of S0 galaxies, in order to unravel their formation histories. We find that PNe are good tracers of stellar kinematics in these systems. We show that the recovered kinematics are largely dominated by rotational motion, although with significant random velocities in most cases.

There are only a few tracers available to probe the kinematics of individual early-type galaxies beyond one effective radius. Here we directly compare a sample of planetary nebulae (PNe), globular clusters (GCs) and galaxy starlight velocities out to ~4 effective radii, in the S0 galaxy NGC 2768. Using a bulge-to-disk decomposition of a K-band image we assign PNe and starlight to either the disk or the bulge. We show that the bulge PNe and bulge starlight follow the same radial density distribution as the red subpopulation of GCs, whereas the disk PNe and disk starlight are distinct components. We find good kinematic agreement between the three tracers to several effective radii (and with stellar data in the inner regions). Further support for the distinct nature of the two galaxy components come from our kinematic analysis. After separating the tracers into bulge and disk components we find the bulge to be a slowly rotating pressure-supported system, whereas the disk reveals a rapidly rising rotation curve with a declining velocity dispersion profile. The resulting V/sigma ratio for the disk resembles that of a spiral galaxy and hints at an origin for NGC 2768 as a transformed late-type galaxy. A two-component kinematic analysis for a sample of S0s will help to elucidate the nature of this class of galaxy.

To investigate the origins of S0 galaxies, we present a new method of analyzing their stellar kinematics from discrete tracers such as planetary nebulae. This method involves binning the data in the radial direction so as to extract the most general possible non-parametric kinematic profiles, and using a maximum likelihood fit within each bin in order to make full use of the information in the discrete kinematic tracers. Both disk and spheroid kinematic components are fitted, with a two-dimensional decomposition of imaging data used to attribute to each tracer a probability of membership in the separate components. Likelihood clipping also allows us to identify objects whose properties are not consistent with the adopted model, rendering the technique robust against contaminants and able to identify additional kinematic features. The method is first tested on an N-body simulated galaxy to assess possible sources of systematic error associated with the structural and kinematic decomposition, which are found to be small. It is then applied to the S0 system NGC~1023, for which a planetary nebula catalogue has already been released and analyzed by (Noordermeer et al., 2008). The correct inclusion of the spheroidal component allows us to show that, contrary to previous claims, the stellar kinematics of this galaxy are indistinguishable from those of a normal spiral galaxy, indicating that it may have evolved directly from such a system via gas stripping or secular evolution. The method also successfully identifies a population of outliers whose kinematics are different from those of the main galaxy; these objects can be identified with a stellar stream associated with the companion galaxy NGC~1023A.

As part of our current programme to test LCDM predictions for dark matter (DM) haloes using extended kinematical observations of early-type galaxies, we present a dynamical analysis of the bright elliptical galaxy NGC 4374 (M84) based on ~450 Planetary Nebulae (PNe) velocities from the PN.Spectrograph, along with extended long-slit stellar kinematics. This is the first such analysis of a galaxy from our survey with a radially constant velocity dispersion profile. We find that the spatial and kinematical distributions of the PNe agree with the field stars in the region of overlap. The velocity kurtosis is consistent with zero at almost all radii. We construct a series of Jeans models, fitting both velocity dispersion and kurtosis to help break the mass-anisotropy degeneracy. Our mass models include DM halos either with shallow cores or with central cusps as predicted by cosmological simulations - along with the novel introduction in this context of adiabatic halo contraction from baryon infall. Both classes of models confirm a very massive dark halo around NGC 4374, demonstrating that PN kinematics data are well able to detect such haloes when present. Considering the default cosmological mass model, we confirm earlier suggestions that bright galaxies tend to have halo concentrations higher than LCDM predictions, but this is found to be solved if either a Salpeter IMF or adiabatic contraction with a Kroupa IMF is assumed. Thus for the first time a case is found where the PN dynamics may well be consistent with a standard dark matter halo. A cored halo can also fit the data, and prefers a stellar mass consistent with a Salpeter IMF. The less dramatic dark matter content found in lower-luminosity "ordinary" ellipticals suggests a bimodality in the halo properties which may be produced by divergent baryonic effects during their assembly histories.

We present new planetary nebulae (PNe) positions, radial velocities, and magnitudes for 6 early-type galaxies obtained with the Planetary Nebulae Spectrograph, their two-dimensional velocity and velocity dispersion fields. We extend this study to include an additional 10 early-type galaxies with PNe radial velocity measurements available from the literature, to obtain a broader description of the outer-halo kinematics in early-type galaxies. These data extend the information derived from stellar kinematics to typically up to ~8 Re. The combination of photometry, stellar and PNe kinematics shows: i) good agreement between the PNe number density and the stellar surface brightness in the region where the two data sets overlap; ii) good agreement between PNe and stellar kinematics; iii) that the mean rms velocity profiles fall into two groups: with of the galaxies characterized by slowly decreasing profiles and the remainder having steeply falling profiles; iv) a larger variety of velocity dispersion profiles; v) that twists and misalignments in the velocity fields are more frequent at large radii, including some fast rotators; vi) that outer haloes are characterised by more complex radial profiles of the specific angular momentum-related lambda_R parameter than observed within 1Re; vii) that many objects are more rotationally dominated at large radii than in their central parts; and viii) that the halo kinematics are correlated with other galaxy properties, such as total luminosity, isophotal shape, total stellar mass, V/sigma, and alpha parameter, with a clear separation between fast and slow rotators.

We present first results of a study of the halo kinematics for a sample of early type galaxies using planetary nebulae (PNe) as kinematical tracers. PNe allow to extend up to several effective radii (Re) the information from absorption line kinematics (confined to within 1 or 2 Re), providing valuable information and constraints for merger simulations and galaxy formation models. We find that the specific angular momentum per unit mass has a more complex radial dependence when the halo region is taken into account and that the halo velocity dispersion is related to the total galaxy luminosity, isophotal shape, and number of PNe per unit of luminosity

We present new Planetary Nebula Spectrograph observations of the ordinary elliptical galaxy NGC 4494, resulting in positions and velocities of 255 PNe out to 7 effective radii (25 kpc). We also present new wide-field surface photometry from MMT/Megacam, and long-slit stellar kinematics from VLT/FORS2. The spatial and kinematical distributions of the PNe agree with the field stars in the region of overlap. The mean rotation is relatively low, with a possible kinematic axis twist outside 1 Re. The velocity dispersion profile declines with radius, though not very steeply, down to ~70 km/s at the last data point. We have constructed spherical dynamical models of the system, including Jeans analyses with multi-component LCDM-motivated galaxies as well as logarithmic potentials. These models include special attention to orbital anisotropy, which we constrain using fourth-order velocity moments. Given several different sets of modelling methods and assumptions, we find consistent results for the mass profile within the radial range constrained by the data. Some dark matter (DM) is required by the data; our best-fit solution has a radially anisotropic stellar halo, a plausible stellar mass-to-light ratio, and a DM halo with an unexpectedly low central density. We find that this result does not substantially change with a flattened axisymmetric model. Taken together with other results for galaxy halo masses, we find suggestions for a puzzling pattern wherein most intermediate-luminosity galaxies have very low concentration halos, while some high-mass ellipticals have very high concentrations. We discuss some possible implications of these results for DM and galaxy formation.