KMOS study of the mass accretion rate from Class I to Class II in NGC 1333
Authors: Fiorellino Eleonora , Manara Carlo F. , Nisini Brunella , Ramsay Suzanne , Antoniucci Simone , Giannini Teresa , Biazzo Katia , Alcalà Juan , Fedele Davide 
Published: March 2021

Context. The mass accretion rate (Ṁacc) is the fundamental parameter to understand the process of mass assembly that results in the formation of a low-mass star. This parameter has been largely studied in Classical T Tauri stars in star-forming regions with ages of ∼1 − 10 Myr. However, little is known about the accretion properties of young stellar objects (YSOs) in younger regions and early stages of star formation, such as in the Class 0/I phases.
Aims. We present new near-infrared spectra of 17 Class I/Flat and 35 Class II sources located in the young (< 1 Myr) NGC 1333 cluster, acquired with the KMOS instrument at the Very Large Telescope. Our goal is to study whether the mass accretion rate evolves with age, as suggested by the widely adopted viscous evolution model, by comparing the properties of the NGC 1333 members with samples of older regions.
Methods. For the Class II sources in our sample, we measured the stellar parameters (SpT, AV, and L⋆) through a comparison of the IR spectra with a grid of non-accreting Class III stellar templates. We then computed the accretion luminosity by using the known correlation between Lacc and the luminosity of HI lines (Paβ and Brγ). For the Class I sample, where the presence of a large IR excess makes it impossible to use the same spectral typing method, we applied a procedure that allowed us to measure the stellar and accretion luminosity in a self-consistent way. Mass accretion rates Ṁacc were then measured once masses and radii were estimated adopting suitable evolutionary tracks.
Results. The NGC 1333 Class II sources of our sample have Lacc ∼ 10−4 − 1 L⊙ and Ṁacc ∼ 10−11 − 10−7 M⊙ yr−1. We find a correlation between accretion and stellar luminosity in the form of log Lacc = (1.5 ± 0.2)log L⋆ + ( − 1.0 ± 0.1), and a correlation between the mass accretion rate and stellar mass in the form of log Ṁacc = (2.6 ± 0.9) log M⋆ + (−7.3 ± 0.7). Both correlations are compatible within the errors with the older Lupus star-forming region, while only the latter is consistent with results from Chamaeleon I. The Class I sample shows larger accretion luminosities (∼10−2 − 102 L⊙) and mass accretion rates (∼10−9 − 10−6 M⊙ yr−1) with respect to the Class II stars of the same cloud. However, the derived mass accretion rates are not sufficiently high to build up the inferred stellar masses, assuming steady accretion during the Class I lifetime. This suggests that the sources are not in their main accretion phase and that most of their mass has already been accumulated during a previous stage and/or that the accretion is an episodic phenomenon. We show that some of the targets originally classified as Class I through Spitzer photometry are in fact evolved or low accreting objects. This evidence can have implications for the estimated protostellar phase lifetimes.
Conclusions. The accretion rates of our sample are larger in more embedded and early stage YSOs. Further observations of larger samples in young star-forming regions are needed to determine if this is a general result. In addition, we highlight the importance of spectroscopic surveys of YSOs to confirm their classification and perform a more correct estimate of their lifetime.