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.

The origin of the arc-shaped Sh2-296 nebula is still unclear. Mainly due to its morphology, the nebula has been suggested to be a 0.5 Myr-old supernova remnant (SNR) that could be inducing star formation in the CMa OB1 association. We aim to show, for the first time, that the nebula is part of a large, shell-like structure, which we have designated the ``CMa shell'', enclosing a bubble created by successive supernova (SN) explosions. We identified three runaway stars, associated with bow-shock structures, in the direction of the CMa shell and we investigate the possibility that they have originated in the center of the shell. By analyzing images of the CMa OB1 association at several wavelengths, we clearly see that the Sh2-296 nebula is in fact part of a large structure, which can be approximated by a large (with a diameter of ~60 pc) elliptical shell. Using the recent Gaia-DR2 astrometric data, we trace back the path of the three runaway stars, in order to find their original position in the past, with relation to the CMa shell. We also revise the heating and ionization of the Sh2-296 nebula, by comparing the photon budget provided by the O stars in the region with results from radio observations. We find that the runaway stars have likely been ejected from a Trapezium-like progenitor cluster on three successive SN explosions having taken place ~6, ~2 and ~1 Myr ago. We also show that the few late-type O stars in the region cannot explain the ionization of the Sh~2-296 nebula and other mechanisms need to be at work. We argue that, though we now have evidence for several SNe events in the CMa OB1 association, the SNe probably played a minor role in triggering star formation in these clouds. In contrast, the CMa OB1 association, as it is now, likely testifies to the last stages of a star-forming region.

The Canis Major OB1 Association has an intriguing scenario of star formation, especially in the Canis Major R1 (CMa R1) region traditionally assigned to a reflection nebula, but in reality an ionized region. This work is focused on the young stellar population associated to CMa R1, for which our previous results from ROSAT, optical and near-infrared data had revealed two stellar groups with different ages, suggesting a possible mixing of populations originated from distinct star-formation episodes. The X-ray data allow the detected sources to be characterized according to hardness ratios, light curves and spectra. Estimates of mass and age were obtained from the 2MASS catalogue, and used to define a complete subsample of stellar counterparts, for statistical purposes. A catalogue of 387 XMM-Newton sources is provided, 78% being confirmed as members or probable members of the CMa R1 association. Flares were observed for 13 sources, and the spectra of 21 bright sources could be fitted by a thermal plasma model. Mean values of fits parameters were used to estimate X-ray luminosities. We found a minimum value of log(L$_X$[erg/s]) = 29.43, indicating that our sample of low-mass stars (M$_\star$ $\leq$ 0.5 M$_\odot$), being faint X-ray emitters, is incomplete. Among the 250 objects selected as our complete subsample (defining our best sample), 171 are found to the East of the cloud, near Z CMa and dense molecular gas, 50% of them being young (< 5 Myr) and 30% being older (> 10 Myr). The opposite happens to the West, near GU CMa, in areas lacking molecular gas: among 79 objects, 30% are young and 50% are older. These findings confirm that a first episode of distributed star formation occurred in the whole studied region ~10 Myr ago and dispersed the molecular gas, while a second, localized episode (< 5 Myr) took place in the regions where molecular gas is still present.

A large sample of young stellar groups is analysed aiming to investigate their clustering properties and dynamical evolution. A comparison of the Q statistical parameter, measured for the clusters, with the fractal dimension estimated for the projected clouds shows that 52% of the sample has substructures and tends to follow the theoretically expected relation between clusters and clouds, according to calculations for artificial distribution of points. The fractal statistics was also compared to structural parameters revealing that clusters having radial density profile show a trend of parameter s increasing with mean surface stellar density. The core radius of the sample, as a function of age, follows a distribution similar to that observed in stellar groups of Milky Way and other galaxies. They also have dynamical age, indicated by their crossing time that is similar to unbound associations. The statistical analysis allowed us to separate the sample into two groups showing different clustering characteristics. However, they have the same dynamical evolution, since the whole sample has been revealed as expanding objects, for which the substructures seem to have not been erased. These results are in agreement with simulations that adopt low surface densities and models under supervirial conditions.

We studied a sample of stars associated with the Sh 2-296 nebula, part of the reflection nebulae complex in the region of Canis Major (CMa R1). Our sample corresponds to optical counterparts of X-ray sources detected from observations with the XMM-Newton satellite, which revealed dozens of possible low-mass young stars not yet known in this region. A sample of 58 young star candidates were selected based on optical spectral features, mainly H{\alpha} and lithium lines, observed with multi-objects spectroscopy performed by the Gemini South telescope. Among the candidates, we find 41 confirmed T Tauri and 15 very likely young stars. Based on the H{\alpha} emission, the T Tauri stars were distinguished between classical (17%) and weak-lined (83%), but no significant difference was found in the age and mass distribution of these two classes. The characterization of the sample was complemented by near- and mid-infrared data, providing an estimate of ages and masses from the comparison with pre-main-sequence evolutionary models. While half of the young stars have an age of 1-2 Myrs or less, only a small fraction (~25%) shows evidence of IR excess revealing the presence of circumstellar discs. This low fraction is quite rare compared to most young star-forming regions, suggesting that some external factor has accelerated the disc dissipation.

Four young star clusters were studied in order to characterize their anomalous extinction or variable reddening that could be due to a possible contamination by dense clouds or circumstellar effects. The extinction law (Rv) was evaluated by adopting two methods: (i) the use of theoretical expressions based on the colour-excess of stars with known spectral type, and (ii) the analysis of two-colour diagrams, where the slope of observed colours distribution is compared to the normal distribution. An algorithm to reproduce the zero age main sequence (ZAMS) reddened colours was developed in order to derive the average visual extinction (Av) that provides the best fitting of the observational data. The structure of the clouds was evaluated by means of statistical fractal analysis, aiming to compare their geometric structure with the spatial distribution of the cluster members. The cluster NGC 6530 is the only object of our sample showing anomalous extinction. In average, the other clusters are suffering normal extinction, but several of their members, mainly in NGC 2264, seem to have high Rv, probably due to circumstellar effects. The ZAMS fitting provides Av values that are in good agreement with those found in the literature. The fractal analysis shows that NGC 6530 has a centrally concentrated distribution of stars that is different of the sub-structures found in the density distribution of the cloud projected in the Av map, suggesting that the original cloud has been changed with the cluster formation. On the other hand, the fractal dimension and the statistical parameters of Berkeley 86, NGC 2244, and NGC 2264 indicate a good cloud-cluster correlation, when compared to other works based on artificial distribution of points.