Search

This study focuses on the analysis of mid-infrared variability in a sample of high-mass young stellar objects (YSOs) associated with the cataloged sources from the APEX Telescope Large Area Survey of the Galaxy (ATLASGAL). The Near-Earth Object Wide-field Infrared Survey Explorer Reactivation Mission (NEOWISE) database was used to explore the long-term mid-infrared variability of these high-mass YSOs at a half-year scale. After matching with NEOWISE photometric measurements, a total of 2230 ATLASGAL sources were selected for the variability analysis, out of which 717 were identified as variables. The derived proportions of variables at different evolutionary stages show that the variability rate of high-mass YSOs is highest during the YSO stage and decreases with evolution toward the H ii region stage, resembling the behavior of low-mass YSOs. The variables can be classified into six types based on their light curves, divided into two categories: secular (linear, sin, sin+linear) and stochastic variables (burst, drop, and irregular). The magnitude–color variations observed in ∼160 secular variables can be mainly divided into “bluer when brighter/redder when dimming” and “redder when brighter/bluer when dimming,” likely originating from changes in accretion rate or the effect of extinction due to obscuration. Moreover, several episodic accretion candidates were selected for further observational studies.

We present a statistical study of infrared variability using the Near-Earth Object Wide-field Infrared Survey Explorer (NEOWISE) database for a sample consisting of 1085 high-mass young stellar objects (YSOs) related to 6.7 GHz methanol masers. A total of 383 maser sources were identified as NEOWISE variables and classified in two variability behavior classifications: 204 secular (linear, curved, and periodic) and 179 stochastic (burst, drop, and irregular) variables. Statistical analysis of the properties of these variables (e.g., the dust temperature, bolometric luminosity, hydrogen column density, W4 luminosity, and W1−W2 color) has revealed a potential evolutionary sequence among different light-curve types of variables. There is a possible general evolutionary (from less to more evolved) trend between the three variable types from secular to stochastic to nonvariable. For the specific classifications, the evolutionary trend for secular variables is linear to sin to sin+linear, and for stochastic variables it is burst to irregular to drop. These sequences may reflect the evolution of the envelope or accretion disk of high-mass YSOs, from large to small radii due to gravitational collapse. Although no significant variability correlation was established between the 6.7 GHz methanol maser and the W2-band emission based on the data collected so far, a number of candidates were found for further investigating the accretion burst events via future variability monitoring programs of both mid-infrared emission and masers.

Light variability is an essential attribute of young stellar objects (YSOs) at the pre-main-sequence (PMS) stage, which can be treated as an indicator to trace star formation. We present a statistical study on the infrared NEOWISE variability of new and large-sized PMS stars as well as Herbig Ae/Be (HAeBe) and classical Be (CBe) stars identified from Gaia Data Release 2. A total of 3792 (including 714 HAeBe), 417, and 473 NEOWISE variable sources were detected from the catalog of 8470 PMS (including 1361 HAeBe), 693 CBe, and 1309 uncertain (possibly belonging to either HAeBe or CBe) stars, respectively. We found that 80% of the PMS variables are dominated by irregular infrared variability, whereas variables in the CBe and uncertain categories exhibit comparable contributions from secular (linear, curved, and periodic) and stochastic (burst, drop, and irregular) variables. This, in combination with infrared W1–W2 colors, W4 luminosities, and the spectral energy distribution slope at near- and mid-infrared wavelengths, suggests different evolutionary stages of the PMS (including HAeBe), uncertain, and CBe stars, from less to more evolved. A statistical analysis of variable source properties (e.g., period and fractional amplitude) further reveals a potential sequence possibly associated with disk evolution for different types of light curves (from linear, to curved, periodic, stochastic, and nonvariable stars). Possible mechanisms for the variable sources are suggested to be related to either extinction via circumstellar matter, or contribution of multiple systems surrounding the YSOs, as opposed to variable accretion (e.g., episodic accretion or accretion burst).