Draft version June 18, 2021 Typeset using LATEX twocolumn style in AASTeX63 Stars with Photometrically Young Gaia Luminosities Around the Solar System (SPYGLASS) I: Mapping Young Stellar Structures and their Star Formation Histories Ronan Kerr,1 Aaron C. Rizzuto,1, ∗ Adam L. Kraus,1 and Stella S. R. Offner1 1Department of Astronomy, University of Texas at Austin 2515 Speedway, Stop C1400 Austin, Texas, USA 78712-1205 (Accepted May 16, 2021) ABSTRACT Young stellar associations hold a star formation record that can persist for millions of years, revealing the progression of star formation long after the dispersal of the natal cloud. To identify nearby young stellar populations that trace this progression, we have designed a comprehensive framework for the identification of young stars, and use it to identify ∼3×104 candidate young stars within a distance of 333 pc using Gaia DR2. Applying the HDBSCAN clustering algorithm to this sample, we identify 27 top-level groups, nearly half of which have little to no presence in previous literature. Ten of these groups have visible substructure, including notable young associations such as Orion, Perseus, Taurus, and Sco-Cen. We provide a complete subclustering analysis on all groups with substructure, using age estimates to reveal each region's star formation history. The patterns we reveal include an apparent star formation origin for Sco-Cen along a semicircular arc, as well as clear evidence for sequential star formation moving away from that arc with a propagation speed of ∼4 km s−1 (∼4 pc Myr−1). We also identify earlier bursts of star formation in Perseus and Taurus that predate current, kinematically identical active star-forming events, suggesting that the mechanisms that collect gas can spark multiple generations of star formation, punctuated by gas dispersal and cloud regrowth. The large spatial scales and long temporal scales on which we observe star formation offer a bridge between the processes within individual molecular clouds and the broad forces guiding star formation at galactic scales. Keywords: young stars | open clusters and associations | star formation | T Tauri stars | Sco-Cen | Taurus | Perseus | Orion 1. INTRODUCTION 2005; Tobin et al. 2009; Kirk et al. 2013, 2017; Kerr Most young stars are not found in isolation, instead et al. 2019), studies of clusters and associations provide a residing in co-moving star clusters or associations (e.g., much longer-term view of star formation that has unique Lada & Lada 2003; Krumholz et al. 2019). These stellar value. Rather than providing a snapshot in time, a de- overdensities are direct remnants of the molecular cloud tailed study of a young association can trace back tens of that collapsed to create them, and as such they preserve millions of years of star formation, enabling a complete study of a star-forming event from its onset to present arXiv:2105.09338v2 [astro-ph.GA] 16 Jun 2021 significant information on the structure and dynamics of those parent clouds (Elmegreen & Lada 1977; Krause day, and by extension the processes that drove the evolu- et al. 2020). While studies of sites of active star forma- tion of the population (e.g.,de Zeeuw et al. 1999; Pecaut tion are popular due to the presence of both young stars & Mamajek 2016; Wright & Mamajek 2018). Further- and dense gas for which exquisite dynamical studies can more, as these associations persist long after gas disper- be performed (e.g., Palla & Stahler 1999; Hatchell et al. sal, the study of stellar populations through associations may also enable the identification of more rapid forms of star formation that do not have known equivalents in Corresponding author: Ronan Kerr active star-forming sites. [email protected] Some processes in star formation require a complete ∗ 51 Pegasi b Fellow record of the star-forming event, not just a snapshot, to 2 Kerr et al. be properly studied in nature. One example is sequential the young stars together and distinguish a young associ- star formation, a process in which previous generations ation from the field. Improvements in measurements of of stars compress the cloud beside them, which can then parallax and proper motions can therefore significantly collapse to form stars, producing a self-sustaining cy- expand our knowledge of these associations, and such cle of star formation that can slowly propagate across developments frequently result in the addition of new an entire molecular cloud (Elmegreen & Lada 1977). members and the revision of associations' known spatial Most cases where a sequential process has been sug- extents (e.g., de Zeeuw et al. 1999; Preibisch & Mamajek gested include just two generations of star formation: 2008; Rizzuto et al. 2015; Zari et al. 2018; Kraus et al. one recently-formed generation powering an H II region, 2017; Luhman 2018). With its unprecedented breadth and one site of active star formation triggered in a shell of spatial, kinematic, and photometric data, Gaia Data that the previous generation compressed (e.g., Lee et al. Release 2 provides a data set capable of dramatically 2005; Maaskant et al. 2011; Nony et al. 2021). Given expanding our knowledge of young stellar populations that these processes are capable of continuing without (Gaia Collaboration et al. 2018a). Gaia photometry can limit as long as unused gas remains, the current view of be combined with parallaxes to identify stars based on sequential star formation has yet to explore large scales their high locations on the HR diagram, and those stars in both time and space. can then be clustered according to their spatial coordi- One of the greatest strengths in using clusters and as- nates and transverse motions. sociations as a record of past star formation lies in the There have been multiple recent searches for stellar time over which detailed information can be extracted. populations in the solar neighborhood enabled by Gaia Having such an extensive star formation record allows DR2, including Sim et al.(2019), which reports the dis- for the study of long-lived star-forming processes, while covery of 207 new open clusters, and Kounkel & Covey also revealing unexpected anomalies in time, such as (2019), which notes the identification of 1901 stellar periods of dormancy in the star formation record. As overdensities in Gaia DR2. However, these surveys in- simulations become increasingly sophisticated and be- cluded stars of all ages in their clustering, so young pop- gin to include more physical processes (e.g., Grudi´cet al. ulations may not always stand out above the field den- 2020a), having a robust record of star formation, com- sity of the older stars. The work of Zari et al.(2018) plete with currently unexplained features can provide does focus on these younger populations, separating out critical comparisons capable of testing new theoretical pre-main sequence stars quite effectively and revealing frameworks. multiple associations in the form of stellar overdensi- While more spatially compact open clusters have ties. While Zari et al.(2018) revealed substantial young long ago been discovered and catalogued (e.g., Trum- stellar populations, it did not include a clustering anal- pler 1921; Klein Wassink 1927; Mermilliod 1995), much ysis, and as such the structures present within those sparser stellar associations are considerably more chal- populations have yet to be rigorously defined. Some re- lenging to identify due to relatively sparse on-sky densi- cent investigations have focused on young stellar popula- ties and large spatial extents, often barring identification tions while also performing a complete clustering analy- from on-sky density. For the nearest populations, stars sis, however these have all been restricted to individual with strong indicators of youth such as protoplanetary associations, such as Zari et al.(2019), which studies disks have typically been used as signposts for the iden- Orion, and Cantat-Gaudin et al.(2019b), which focuses tification of associated populations nearby, such as TW on Vela. No paper to date has performed a spatially un- Hydrae and β Pictoris (e.g., Kastner et al. 1997; Zucker- biased, all-sky search for young stellar populations with man et al. 2001). Beyond about 50-100 pc, however, rec- both robust youth diagnostics and a broad clustering ognizing stars with youth indicators and finding poten- analysis. tial companions becomes increasingly difficult, leaving In this paper, we present the deepest comprehensive the population of low-mass, more distant associations study of young stellar populations in the solar neighbor- largely unexplored. Due to their unbound nature, small hood to date. In Section2, we present the Gaia DR2 velocities relative to a typical field star, and geometric data set used for this work. We outline our methods effects from their wide spatial distributions, the kine- for the identification of young stars, age estimation, and matics of associations are often difficult to disentangle clustering in Section3, while a detailed analysis into the from the field. The effective identification of young stel- populations we identify is performed in Section4. Pos- lar populations therefore requires both the suppression sible implications of this work with respect to the study of older background populations and the use of accurate of star formation progression in the solar neighborhood kinematics and 3-d spatial coordinates to properly group Young Stellar Populations 3 are discussed in Section5, while we conclude in Section In addition, we impose a restriction on the BP/RP 6. flux excess factor3 E, which is defined as the sum of the fluxes in the Gaia BP and RP photometric bands 2. DATA divided by the flux in the G band. The configuration of 2.1. Gaia Astrometry and Photometry the Gaia photometric bands implies that E should be In this paper, we search for nearby young stars among slightly larger than 1 for stars with good photometric a large sample drawn from Gaia Data Release 2 (Gaia measurements, however the BP and RP bands are vul- Collaboration et al.