Variability Of Young Stellar Objects Research Articles

Submillimeter and Mid-Infrared Variability of Young Stellar Objects in the M17SWex Intermediate-mass Star-forming Region

We present a comprehensive analysis of young stellar object (YSO) variability within the M17 Southwest Extension (M17 SWex), using 3.5 yr of monitoring data from the James Clerk Maxwell Telescope (JCMT) Transient Survey at submillimeter and 9 yr from the NEOWISE mission at mid-infrared (mid-IR). Our study encompasses observations of 147 bright submillimeter peaks identified within our deep JCMT coadded map as well as 156 YSOs in NEOWISE W1 and 179 in W2 that were previously identified in Spitzer surveys. We find three robust submillimeter variables: two are candidate YSOs and one is a likely extragalactic source. At mid-IR wavelengths, our analysis reveals secular and stochastic variability in 47 YSOs, with the highest fraction of secular variability occurring at the earliest evolutionary stage. This is similar to what has previously been observed for low-mass YSO variability within the Gould Belt. However, we observe less overall variability in M17 SWex at both the submillimeter and mid-IR. We suspect that this lower fraction is due to the greater distance to M17 SWex. Our findings showcase the utility of multiwavelength observations to better capture the complex variability phenomena inherent to star formation processes and demonstrate the importance of years-long monitoring of a diverse selection of star-forming environments.

Mid-infrared Variability of Young Stellar Objects on Timescales of Days to Years

Variability in the brightness of young stellar objects (YSOs) is a common phenomenon that can be caused by changes in various factors, including accretion, extinction, disk morphology, interactions between the disk and the stellar photosphere, and the rotation of hot or cold magnetic spots on the stellar photosphere. Analyzing the variability on different timescales provides insight into the mechanisms driving the changes in the brightness of YSOs. We investigate the variability of YSOs on both long and short timescales using two mid-IR data sets: the Near-Earth Object Wide-field Infrared Survey Explorer (NEOWISE) 7.5 yr W2 (4.6 μm) data and the YSOVAR 40 day Spitzer/IRAC2 (4.5 μm) data, respectively. We classify the variability types in each timescale following Park et al. We find a higher detection rate of variable sources in the short term (77.6%) compared to the long term (43.0%) due to the higher sensitivity of the Spitzer observations. In addition, the higher cadence of the YSOVAR data results in the weeks-long short-term variability being mostly secular, while the years-long long-term variability explored with the coarsely sampled NEOWISE data is mostly stochastic. By crossmatching the two catalogs, we statistically analyze the variability types exhibited by YSOs across both timescales. The long-term variability amplitude is mostly three times (up to 10 times) greater than the short-term variability. Furthermore, we evaluate variability on very short (1–2 days) timescales and recover a trend of the increasing amplitude of variability as the timescales increase. By comprehensively analyzing the variability of YSOs over various timescales, we contribute to a deeper understanding of the underlying mechanisms driving their variability.

On the incidence of episodic accretion in Class I YSOs from VVV

ABSTRACT Episodic accretion is one of the competing models to explain the observed luminosity spread in young stellar clusters. These short-lived high accretion events could also have a strong impact on planet formation. Observations of high-amplitude variability in young stellar objects (YSOs) due to large changes in the accretion rate provide direct observational evidence for episodic accretion. However, there are still uncertainties in the frequency of these events and if episodic accretion is universal among YSOs. To determine the frequency of outbursts in Class I YSOs, we built a large and robust sample of objects at this evolutionary stage, and searched for high-amplitude near-infrared (ΔKS > 2 mag) variability in the VIRAC2 database of the Vista Variables in the Via Lactea survey. By complementing with near-IR (2MASS and DENIS) and mid-IR (WISE/Neo-WISE) data, we find that from ∼7000 Class I YSOs, 97 objects can be classified as eruptive variable YSOs. The duration of the outbursts vary from a few months to longer than 9 yr, and cover a similar range of amplitudes. Values of ΔKS > 5 mag, however, are only observed in outbursts with duration longer than 9 yr. When considering different effects of completeness and contamination, we estimate that the incidence of episodic accretion in Class I YSOs is between 2 and 3 per cent. Finally, we determine a recurrence time-scale of long-term outbursts (a.k.a FUors) of $\tau =1.75^{+1.12}_{-0.87}$ kyr. The latter value agrees with previous estimates and is in line with the expectations of higher frequency of FUor outbursts during younger stages of evolution.

Photometric and spectroscopic monitoring of YSOs in nearby star-forming regions – I. Eruptive YSOs

ABSTRACT Mid-infrared (mid-IR) variability in young stellar objects (YSOs) is driven by several physical mechanisms, which produce a variety of amplitudes and light-curve shapes. One of these mechanisms, variable disc accretion, is predicted by models of episodic accretion to drive secular variability, including in the mid-IR. Because the largest accretion bursts are rare, adding new objects to the YSO eruptive variable class aids our understanding of the episodic accretion phenomenon and its possible impact on stellar and planetary formation. A previous analysis of 6.5 yr of NeoWISE light curves (3–5 μm) of ∼7000 nearby YSOs found an increase in the fraction of variability and variability amplitude for objects at younger stages of evolution. To help interpret these light curves, we have obtained low- and high-resolution near-IR spectra of 78 objects from this sample of YSOs. In this work, we present the analysis of nine nearby YSOs (d<1 kpc) that show the characteristics of known classes of eruptive variable YSOs. We find one FU Orionis (FUor)-like source, one EX Lupi type object, and six YSOs with mixed characteristics or V1647 Ori like objects. The varied characteristics observed in our sample are consistent with recent discoveries of eruptive YSOs. We discuss how a wide range in YSO outburst parameters (central mass, maximum accretion rate during outburst, evolutionary stage, and/or instability leading to the outburst) may play a significant role in the observed spectrophotometric properties of YSO outbursts.

A systematic survey of millimetre-wavelength flaring variability of young stellar objects in the Orion Nebula Cluster

ABSTRACT High-energy processes are ubiquitous even in the earliest stages of protostellar evolution. Motivated by the results of our systematic search for intense centimetre radio flares in Young Stellar Objects (YSOs) and by rare findings of strong millimetre-wavelength variability, we have conducted a systematic search for such variability in the Orion Nebula Cluster (ONC) using Atacama Large Millimetre/submillimetre Array (ALMA). Rapid variability on time-scales of minutes to hours in the (centimetre) millimetre-wavelength range indicates (gyro)synchrotron radiation. Additionally, mass accretion will also affect the millimetre-wavelength luminosity but typically on longer time-scales. Beyond studies of individual YSOs, our characterization of strong millimetre-wavelength variability with ALMA in the ONC sets first systematic constraints on the occurrence of such variability in a large number of YSOs (∼130). We report the discovery of an order of magnitude millimetre-flare within just a few minutes from a known YSO previously reported as a radio flaring source at cm-wavelengths (the ‘ORBS’ source). We also present an assessment of the systematic variability effects caused by the use of time-sliced imaging of a complex region. These are mostly due to the impact of a changing synthesized beam throughout the observations. We use simulated ALMA observations to reproduce and quantify these effects, and set a lower limit for the variability that can be studied using our method in a complex region such as the ONC. Our results demonstrate that the utility of time domain analysis of YSOs extends into the millimetre-wavelength range, potentially interfering with the conversion of observed fluxes into dust masses.

Young Stellar Objects, Accretion Disks, and Their Variability with Rubin Observatory LSST

Vera C. Rubin Observatory, through the Legacy Survey of Space and Time (LSST), will allow us to derive a panchromatic view of variability in young stellar objects (YSOs) across all relevant timescales. Indeed, both short-term variability (on timescales of hours to days) and long-term variability (months to years), predominantly driven by the dynamics of accretion processes in disk-hosting YSOs, can be explored by taking advantage of the multiband filters option available in Rubin LSST, in particular the u, g, r, i filters that enable us to discriminate between photospheric stellar properties and accretion signatures. The homogeneity and depth of sky coverage that will be achieved with LSST will provide us with a unique opportunity to characterize the time evolution of disk accretion as a function of age and varying environmental conditions (e.g., field crowdedness, massive neighbors, metallicity) by targeting different star-forming regions. In this contribution to the Rubin LSST Survey Strategy Optimization Focus Issue, we discuss how implementing a dense observing cadence to explore short-term variability in YSOs represents a key complementary effort to the Wide–Fast–Deep observing mode that will be used to survey the sky over the full duration of the main survey (≈10 yr). The combination of these two modes will be vital to investigate the connection between the inner-disk dynamics and longer-term eruptive variability behaviors, such as those observed on EX Lupi–type objects.

A Zwicky Transient Facility Look at Optical Variability of Young Stellar Objects in the North America and Pelican Nebulae Complex

We present a study of 323 photometrically variable young stellar objects that are likely members of the North America and Pelican nebulae star-forming region. To do so, we utilize over two years of data in the g and r photometric bands from the Zwicky Transient Facility. We first investigate periodic variability, finding 46 objects (∼15% of the sample) with significant periods that phase well and can be attributed to stellar rotation. We then use the quasiperiodicity (Q) and flux asymmetry (M) variability metrics to assign morphological classifications to the remaining aperiodic light curves. Another ∼39% of the variable star sample beyond the periodic (low Q) sources are also flux-symmetric, but with a quasiperiodic (moderate Q) or stochastic (high Q) nature. Concerning flux-asymmetric sources, our analysis reveals ∼14% bursters (high negative M) and ∼29% dippers (high positive M). We also investigate the relationship between variability slopes in the g versus g − r color–magnitude diagram, and the light-curve morphological classes. Burster-type objects have shallow slopes, while dipper-type variables tend to have higher slopes that are consistent with extinction-driven variability. Our work is one of the earliest applications of the Q and M metrics to ground-based data. We therefore contrast the Q values of high-cadence and high-precision space-based data, for which these metrics were designed, with Q determinations resulting from degraded space-based light curves that have the cadence and photometric precision characteristic of ground-based data.

Quantifying Variability of Young Stellar Objects in the Mid-infrared Over 6 Years with the Near-Earth Object Wide-field Infrared Survey Explorer

Abstract Variability in young stellar objects (YSOs) can be caused by various time-dependent phenomena associated with star formation, including accretion rates, geometric changes in the circumstellar disks, stochastic hydromagnetic interactions between stellar surfaces and inner-disk edges, reconnections within the stellar magnetosphere, and hot/cold spots on stellar surfaces. We uncover and characterize ∼1700 variables from a sample of ∼5400 YSOs in nearby low-mass star-forming regions using mid-IR light curves obtained from the 6.5 yr Near-Earth Object Wide-field Infrared Survey Explorer All Sky Survey. The mid-IR variability traces a wide range of dynamical, physical, and geometrical phenomenon. We classify six types of YSO mid-IR variability based on their light curves: secular variability (linear, curved, and periodic) and stochastic variability (burst, drop, and irregular). YSOs in earlier evolutionary stages have higher fractions of variables and higher amplitudes for the variability, with the recurrence timescale of FUor-type outbursts (defined here as ΔW1 or ΔW2 > 1 mag followed by inspection of candidates) of ∼1000 yr in the early embedded protostellar phase. Known eruptive young stars and subluminous objects show fractions of variables similar to the fraction (∼55%) found in typical protostars, suggesting that these two distinct types are not distinct in variability over the 6.5 yr timescale. Along with brightness variability, we also find a diverse range of secular color variations, which can be attributed to a competitive interplay between the variable accretion luminosity of the central source and the variable extinction by material associated with the accretion process.

High‐energy emission and its variability in young stellar objects

Young stars show a variety of highly energetic phenomena, from accretion and outflow processes to hot coronal plasmas confined in their outer atmosphere, all regulated by the intense stellar magnetic fields. Many aspects on each of these phenomena are debated, but, most notably, their complex mutual interaction remains obscure. In this work, I report how these phenomena are simultaneously responsible for the high‐energy emission from young stars, with a special focus on the expected and observed variability in the X‐ray band. Investigating variations in the X‐ray emission from young stars allows us to pose constraints on flare and coronal plasma properties, coronal heating, accretion stream properties, and accretion geometries. All these results are important building blocks for constructing a comprehensive picture of the complex magnetosphere of young stars.

Constraining the Rate of Protostellar Accretion Outbursts in the Orion Molecular Clouds

Abstract Outbursts due to dramatic increases in the mass accretion rate are the most extreme type of variability in young stellar objects. We searched for outbursts among 319 protostars in the Orion molecular clouds by comparing 3.6, 4.5, and 24 μm photometry from the Spitzer Space Telescope to 3.4, 4.6, and 22 μm photometry from the Wide-field Infrared Survey Explorer (WISE) obtained ∼6.5 yr apart. Sources that brightened by more than two standard deviations above the mean variability at all three wavelengths were marked as burst candidates, and they were inspected visually to check for false positives due primarily to the reduced angular resolution of WISE compared to Spitzer. We recovered the known burst V2775 Ori (HOPS 223) as well as a previously unknown burst, HOPS 383, which we announced in an earlier paper. No other outbursts were found. With observations over 6.5 yr, we estimate an interval of about 1000 yr between bursts with a 90% confidence interval of 690–40,300 yr. The most likely burst interval is shorter than those found in studies of optically revealed young stellar objects, suggesting that outbursts are more frequent in protostars than in pre-main-sequence stars that lack substantial envelopes.

A survey for variable young stars with small telescopes: First results from HOYS-CAPS

Variability in Young Stellar Objects (YSOs) is one of their primary characteristics. Long-term, multi-filter, high-cadence monitoring of large YSO samples is the key to understand the partly unusual light-curves that many of these objects show. Here we introduce and present the first results of the HOYS-CAPS citizen science project which aims to perform such monitoring for nearby (d<1kpc) and young (age<10Myr) clusters and star forming regions, visible from the northern hemisphere, with small telescopes. We have identified and characterised 466 variable (413 confirmed young) stars in 8 young, nearby clusters. All sources vary by at least 0.2mag in V, have been observed at least 15 times in V, R and I in the same night over a period of about 2yrs and have a Stetson index of larger than 1. This is one of the largest samples of variable YSOs observed over such a time-span and cadence in multiple filters. About two thirds of our sample are classical T-Tauri stars, while the rest are objects with depleted or transition disks. Objects characterised as bursters show by far the highest variability. Dippers and objects whose variability is dominated by occultations from normal interstellar dust or dust with larger grains (or opaque material) have smaller amplitudes. We have established a hierarchical clustering algorithm based on the light-curve properties which allows the identification of the YSOs with the most unusual behaviour, and to group sources with similar properties. We discuss in detail the light-curves of the unusual objects V2492Cyg, V350Cep and 2MASSJ21383981+5708470.

YSOVAR: Mid-infrared Variability among YSOs in the Star Formation Region Serpens South

Abstract We present a time-variability study of young stellar objects (YSOs) in the Serpens South cluster performed at 3.6 and 4.5 μm with the Spitzer Space Telescope; this study is part of the Young Stellar Object VARiability project. We have collected light curves for more than 1500 sources, including 85 cluster members, over 38 days. This includes 44 class I sources, 19 sources with flat spectral energy distributions (SEDs), 17 class II sources, and five diskless YSO candidates. We find a high variability fraction among embedded cluster members of ∼70%, whereas young stars without a detectable disk display no variability. We detect periodic variability for 32 sources with periods primarily in the range of 0.2–14 days and a subset of fast rotators thought to be field binaries. The timescale for brightness changes are shortest for stars with the most photospheric SEDs and longest for those with flat or rising SEDs. While most variable YSOs become redder when fainter, as would be expected from variable extinction, about 10% get bluer as they get fainter. One source, SSTYSV J183006.13−020108.0, exhibits “cyclical” color changes.

Star–Disk Interactions in Multiband Photometric Monitoring of the Classical T Tauri Star GI Tau

Abstract The variability of young stellar objects is mostly driven by star–disk interactions. In long-term photometric monitoring of the accreting T Tauri star GI Tau, we detect extinction events with typical depths of mag that last for days to months and often appear to occur stochastically. In 2014–2015, extinctions that repeated with a quasi-period of 21 days over several months are the first empirical evidence of slow warps predicted by magnetohydrodynamic simulations to form at a few stellar radii away from the central star. The reddening is consistent with and, along with an absence of diffuse interstellar bands, indicates that some dust processing has occurred in the disk. The 2015–2016 multiband light curve includes variations in spot coverage, extinction, and accretion, each of which results in different traces in color–magnitude diagrams. This light curve is initially dominated by a month-long extinction event and a return to the unocculted brightness. The subsequent light curve then features spot modulation with a 7.03 day period, punctuated by brief, randomly spaced extinction events. The accretion rate measured from U-band photometry ranges from to M ⊙ yr−1 (excluding the highest and lowest 5% of high- and low- accretion rate outliers), with an average of M ⊙ yr−1. A total of 50% of the mass is accreted during bursts of M ⊙ yr , which indicates limitations on analyses of disk evolution using single-epoch accretion rates.

HIGHLY VARIABLE YOUNG MASSIVE STARS IN ATLASGAL CLUMPS

High-amplitude variability in Young Stellar Objects (YSOs) is usually associated with episodic accretion events. It has not been observed so far in massive YSOs. Here, the high-amplitude variable star sample of ContrerasPe\~{n}a et al.(2016) has been used to search for highly-variable($\Delta$K$\ge$1\,mag) sources coinciding with dense clumps mapped using the 850\mum continuum emission by the ATLASGAL survey. 18 variable sources are centred on the sub-mm clump peaks, and coincide ($<$1") with a 24$\mu$m point or compact ($<$10") source. 13 of these 18 sources can be fit by YSO models. The 13 variable YSOs(VYSO) have luminosities of $\sim$10$^3$ L$_{\odot}$, an average mass of 8 M$_{\odot}$ and a range of ages up to 10$^6$ yr. 11 of these 13 VYSOs are located in the midst of infrared dark clouds. 9 of the 13 sources have $\Delta$K$>$2 mag, significantly higher compared to the mean variability of the entire VVV sample. The light curves of these objects sampled between 2010-2015 display rising, declining, or quasi-periodic behaviour but no clear periodicity. Light-curve analysis using Plavchan method show that the most prominent phased signals have periods of a few hundred days. The nature and time-scale of variations found in 6.7 Ghz methanol maser emission (MME) in massive stars are similar to that of the VYSO light curves. We argue that the origin of the observed variability is episodic accretion. We suggest that the timescale of a few hundred days may represent the frequency at which a spiralling disk feeds dense gas to the young massive star.

YSOVAR: MID-INFRARED VARIABILITY IN NGC 1333

As part of the Young Stellar Object VARiability (YSOVAR) program, we monitored NGC 1333 for ~35 days at 3.6 and 4.5 um using the Spitzer Space Telescope. We report here on the mid-infrared variability of the point sources in the ~ 10 X ~ 20 area centered on 03:29:06, +31:19:30 (J2000). Out of 701 light curves in either channel, we find 78 variables over the YSOVAR campaign. About half of the members are variable. The variable fraction for the most embedded SEDs (Class I, flat) is higher than that for less embedded SEDs (Class II), which is in turn higher than the star-like SEDs (Class III). A few objects have amplitudes (10-90th percentile brightness) in [3.6] or [4.5]>0.2 mag; a more typical amplitude is 0.1-0.15 mag. The largest color change is >0.2 mag. There are 24 periodic objects, with 40% of them being flat SED class. This may mean that the periodic signal is primarily from the disk, not the photosphere, in those cases. We find 9 variables likely to be 'dippers', where texture in the disk occults the central star, and 11 likely to be 'bursters', where accretion instabilities create brightness bursts. There are 39 objects that have significant trends in [3.6]-[4.5] color over the campaign, about evenly divided between redder-when-fainter (consistent with extinction variations) and bluer-when-fainter. About a third of the 17 Class 0 and/or jet-driving sources from the literature are variable over the YSOVAR campaign, and a larger fraction (~half) are variable between the YSOVAR campaign and the cryogenic-era Spitzer observations (6-7 years), perhaps because it takes time for the envelope to respond to changes in the central source. The NGC 1333 brown dwarfs do not stand out from the stellar light curves in any way except there is a much larger fraction of periodic objects (~60% of variable brown dwarfs are periodic, compared to ~30% of the variables overall).

YSOVAR: MID-INFRARED VARIABILITY AMONG YSOs IN THE STAR FORMATION REGION GGD12-15

We present an IR-monitoring survey with the $Spitzer$ Space Telescope of the star forming region GGD 12-15. Over 1000 objects were monitored including about 350 objects within the central 5 arcminutes which is found to be especially dense in cluster members. The monitoring took place over 38 days and is part of the Young Stellar Object VARiability (YSOVAR) project. The region was also the subject of a contemporaneous 67ks $Chandra$ observation. The field includes 119 previously identified pre-main sequence star candidates. X-rays are detected from 164 objects, 90 of which are identified with cluster members. Overall, we find that about half the objects in the central 5 arcminutes are young stellar objects based on a combination of their spectral energy distribution, IR variability and X-ray emission. Most of the stars with IR excess relative to a photosphere show large amplitude (>0.1 mag) mid-IR variability. There are 39 periodic sources, all but one of these is found to be a cluster member. Almost half of the periodic sources do not show IR excesses. Overall, more than 85% of the Class I, flat spectrum, and Class II sources are found to vary. The amplitude of the variability is larger in more embedded young stellar objects. Most of the Class~I/II objects exhibit redder colors in a fainter state, compatible with time-variable extinction. A few become bluer when fainter, which can be explained with significant changes in the structure of the inner disk. A search for changes in the IR due to X-ray events is carried out, but the low number of flares prevented an analysis of the direct impact of X-ray flares on the IR lightcurves. However, we find that X-ray detected Class II sources have longer timescales for change in the mid-IR than a similar set of non-X-ray detected Class IIs.

NEAR-INFRARED VARIABILITY IN THE ORION NEBULA CLUSTER

Using the United Kingdom Infrared Telescope on Mauna Kea, we have carried out a new near-infrared J, H, K monitoring survey of almost a square degree of the star-forming Orion Nebula Cluster with observations on 120 nights over three observing seasons, spanning a total of 894 days. We monitored ~15,000 stars down to J=20 using the WFCAM instrument, and have extracted 1203 significantly variable stars from our data. By studying variability in young stellar objects (YSOs) in the H-K, K color-magnitude diagram, we are able to distinguish between physical mechanisms of variability. Many variables show color behavior indicating either dust-extinction or disk/accretion activity, but we find that when monitored for longer periods of time, a number of stars shift between these two variability mechanisms. Further, we show that the intrinsic timescale of disk/accretion variability in young stars is longer than that of dust-extinction variability. We confirm that variability amplitude is statistically correlated with evolutionary class in all bands and colors. Our investigations of these 1203 variables have revealed 73 periodic AA Tau type variables, many large-amplitude and long-period (P > 15 day) YSOs, including three stars showing widely-spaced periodic brightening events consistent with circumbinary disk activity, and four new eclipsing binaries. These phenomena and others indicate the activity of long-term disk/accretion variability processes taking place in young stars. We have made the light curves and associated data for these 1203 variables available online.

YSOVAR: MID-INFRARED VARIABILITY OF YOUNG STELLAR OBJECTS AND THEIR DISKS IN THE CLUSTER IRAS 20050+2720

We present a time-variability study of young stellar objects in the cluster IRAS 20050+2720, performed at 3.6 and 4.5 micron with the Spitzer Space Telescope; this study is part of the Young Stellar Object VARiability project (YSOVAR). We have collected light curves for 181 cluster members over 40 days. We find a high variability fraction among embedded cluster members of ca. 70%, whereas young stars without a detectable disk display variability less often (in ca. 50% of the cases) and with lower amplitudes. We detect periodic variability for 33 sources with periods primarily in the range of 2-6 days. Practically all embedded periodic sources display additional variability on top of their periodicity. Furthermore, we analyze the slopes of the tracks that our sources span in the color-magnitude diagram (CMD). We find that sources with long variability time scales tend to display CMD slopes that are at least partially influenced by accretion processes, while sources with short variability time scales tend to display extinction-dominated slopes. We find a tentative trend of X-ray detected cluster members to vary on longer time scales than the X-ray undetected members.

Circumstellar disks revealed byH/Kflux variation gradients

The variability of young stellar objects (YSO) changes their brightness and color preventing a proper classification in traditional color-color and color magnitude diagrams. We have explored the feasibility of the flux variation gradient (FVG) method for YSOs, using $H$ and $K$ band monitoring data of the star forming region RCW\,38 obtained at the University Observatory Bochum in Chile. Simultaneous multi-epoch flux measurements follow a linear relation $F_{H}=\alpha + \beta \cdot F_{K}$ for almost all YSOs with large variability amplitude. The slope $\beta$ gives the mean $HK$ color temperature $T_{var}$ of the varying component. Because $T_{var}$ is hotter than the dust sublimation temperature, we have tentatively assigned it to stellar variations. If the gradient does not meet the origin of the flux-flux diagram, an additional non- or less-varying component may be required. If the variability amplitude is larger at the shorter wavelength, e.g. $\alpha < 0$, this component is cooler than the star (e.g. a circumstellar disk); vice versa, if $\alpha > 0$, the component is hotter like a scattering halo or even a companion star. We here present examples of two YSOs, where the $HK$ FVG implies the presence of a circumstellar disk; this finding is consistent with additional data at $J$ and $L$. One YSO shows a clear $K$-band excess in the $JHK$ color-color diagram, while the significance of a $K$-excess in the other YSO depends on the measurement epoch. Disentangling the contributions of star and disk it turns out that the two YSOs have huge variability amplitudes ($\sim 3-5$\,mag). The $HK$ FVG analysis is a powerful complementary tool to analyze the varying components of YSOs and worth further exploration of monitoring data at other wavelengths.

YOUNG STELLAR OBJECT VARIABILITY (YSOVAR): LONG TIMESCALE VARIATIONS IN THE MID-INFRARED

The YSOVAR (Young Stellar Object VARiability) Spitzer Space Telescope observing program obtained the first extensive mid-infrared (3.6 & 4.5 um) time-series photometry of the Orion Nebula Cluster plus smaller footprints in eleven other star-forming cores (AFGL490, NGC1333, MonR2, GGD 12-15, NGC2264, L1688, Serpens Main, Serpens South, IRAS 20050+2720, IC1396A, and Ceph C). There are ~29,000 unique objects with light curves in either or both IRAC channels in the YSOVAR data set. We present the data collection and reduction for the Spitzer and ancillary data, and define the "standard sample" on which we calculate statistics, consisting of fast cadence data, with epochs about twice per day for ~40d. We also define a "standard sample of members", consisting of all the IR-selected members and X-ray selected members. We characterize the standard sample in terms of other properties, such as spectral energy distribution shape. We use three mechanisms to identify variables in the fast cadence data--the Stetson index, a chi^2 fit to a flat light curve, and significant periodicity. We also identified variables on the longest timescales possible of ~6 years, by comparing measurements taken early in the Spitzer mission with the mean from our YSOVAR campaign. The fraction of members in each cluster that are variable on these longest timescales is a function of the ratio of Class I/total members in each cluster, such that clusters with a higher fraction of Class I objects also have a higher fraction of long-term variables. For objects with a YSOVAR-determined period and a [3.6]-[8] color, we find that a star with a longer period is more likely than those with shorter periods to have an IR excess. We do not find any evidence for variability that causes [3.6]-[4.5] excesses to appear or vanish within our data; out of members and field objects combined, at most 0.02% may have transient IR excesses.