A Simple Calculation in Service of Constraining the Rate of FU Orionis Outburst Events from Photometric Monitoring Surveys

An enigmatic and rare type of young stellar object is the FU Orionis class. The members are interpreted as "outbursting," that is, currently in a state of enhanced accretion by several orders of magnitude relative to the more modest disk-to-star accretion rates measured in typical T Tauri stars. They are key to our understanding of the history of stellar mass assembly and pre-main sequence evolution, as well as critical to consider in the chemical and physical evolution of the circumstellar environment -- where planets form. A common supposition is that *all* T Tauri stars undergo repeated such outbursts, more frequently in their earlier evolutionary stages when the disks are more massive, so as to build up the requisite amount of stellar mass on the required time scale. However, the actual data supporting this traditional picture of episodically enhanced disk accretion are limited, and the observational properties of the known sample of FU Ori objects quite diverse. To improve our understanding of these rare objects, we outline the logic for meaningfully constraining the rate of FU Ori outbursts and present numbers to guide parameter choices in the analysis of time domain surveys.Comment: accepted for publication in Ap

Simulated Performance of Timescale Metrics for Aperiodic Light Curves

Aperiodic variability is a characteristic feature of young stars, massive stars, and active galactic nuclei. With the recent proliferation of time domain surveys, it is increasingly essential to develop methods to quantify and analyze aperiodic variability. We develop three timescale metrics that have been little used in astronomy -- {\Delta}m-{\Delta}t plots, peak-finding, and Gaussian process regression -- and present simulations comparing their effectiveness across a range of aperiodic light curve shapes, characteristic timescales, observing cadences, and signal to noise ratios. We find that Gaussian process regression is easily confused by noise and by irregular sampling, even when the model being fit reflects the process underlying the light curve, but that {\Delta}m-{\Delta}t plots and peak-finding can coarsely characterize timescales across a broad region of parameter space. We make public the software we used for our simulations, both in the spirit of open research and to allow others to carry out analogous simulations for their own observing programs.Comment: 39 pages, 16 figures, 3 tables. To be published in The Astrophysical Journal. Simulation software is available online at http://ascl.net/1408.01

Disk-related Bursts and Fades in Young Stars

We present first results from a new, multiyear, time domain survey of young stars in the North America Nebula complex using the Palomar Transient Factory. Our survey is providing an unprecedented view of aperiodic variability in young stars on timescales of days to years. The analyzed sample covers R_(PTF) ≈ 13.5-18 and spans a range of mid-infrared color, with larger-amplitude optical variables (exceeding 0.4 mag root mean squared) more likely to have mid-infrared evidence for circumstellar material. This paper characterizes infrared excess stars with distinct bursts above or fades below a baseline of lower-level variability, identifying 41 examples. The light curves exhibit a remarkable diversity of amplitudes, timescales, and morphologies, with a continuum of behaviors that cannot be classified into distinct groups. Among the bursters, we identify three particularly promising sources that may represent theoretically predicted short-timescale accretion instabilities. Finally, we find that fading behavior is approximately twice as common as bursting behavior on timescales of days to years, although the bursting and fading duty cycle for individual objects often varies from year to year

CSI 2264: Simultaneous Optical and Infrared Light Curves of Young Disk-bearing Stars in NGC 2264 with CoRoT and Spitzer—Evidence for Multiple Origins of Variability

We present the Coordinated Synoptic Investigation of NGC 2264, a continuous 30 day multi-wavelength photometric monitoring campaign on more than 1000 young cluster members using 16 telescopes. The unprecedented combination of multi-wavelength, high-precision, high-cadence, and long-duration data opens a new window into the time domain behavior of young stellar objects. Here we provide an overview of the observations, focusing on results from Spitzer and CoRoT. The highlight of this work is detailed analysis of 162 classical T Tauri stars for which we can probe optical and mid-infrared flux variations to 1% amplitudes and sub-hour timescales. We present a morphological variability census and then use metrics of periodicity, stochasticity, and symmetry to statistically separate the light curves into seven distinct classes, which we suggest represent different physical processes and geometric effects. We provide distributions of the characteristic timescales and amplitudes and assess the fractional representation within each class. The largest category (>20%) are optical "dippers" with discrete fading events lasting ~1-5 days. The degree of correlation between the optical and infrared light curves is positive but weak; notably, the independently assigned optical and infrared morphology classes tend to be different for the same object. Assessment of flux variation behavior with respect to (circum)stellar properties reveals correlations of variability parameters with Hα emission and with effective temperature. Overall, our results point to multiple origins of young star variability, including circumstellar obscuration events, hot spots on the star and/or disk, accretion bursts, and rapid structural changes in the inner disk

CSI 2264: Characterizing Accretion-Burst Dominated Light Curves for Young Stars in NGC 2264

Based on more than four weeks of continuous high cadence photometric monitoring of several hundred members of the young cluster NGC 2264 with two space telescopes, NASA’s Spitzer and the CNES CoRoT (Convection, Rotation, and planetary Transits), we provide high quality, multi-wavelength light curves for young stellar objects (YSOs) whose optical variability is dominated by short duration flux bursts, which we infer are due to enhanced mass accretion rates. These light curves show many brief – several hour to one day – brightenings at optical and near-infrared (IR) wavelengths with amplitudes generally in the range 5-50% of the quiescent value. Typically, a dozen or more of these bursts occur in a thirty day period. We demonstrate that stars exhibiting this type of variability have large ultraviolet (UV) excesses and dominate the portion of the u − g vs. g − r color-color diagram with the largest UV excesses. These stars also have large Hɑ equivalent widths, and either centrally peaked, lumpy Hɑ emission profiles or profiles with blue-shifted absorption dips associated with disk or stellar winds. Light curves of this type have been predicted for stars whose accretion is dominated by Rayleigh-Taylor instabilities at the boundary between their magnetosphere and inner circumstellar disk, or where magneto-rotational instabilities modulate the accretion rate from the inner disk. Amongst the stars with the largest UV excesses or largest Hɑ equivalent widths, light curves with this type of variability greatly outnumber light curves with relatively smooth sinusoidal variations associated with long-lived hot spots. We provide quantitative statistics for the average duration and strength of the accretion bursts and for the fraction of the accretion luminosity associated with these bursts

The First Extensive Spectroscopic Study of Young Stars in the North America and Pelican Nebulae

We present a spectroscopic survey of over 3400 potential members in the North America and Pelican Nebulae (NAP) using several low-resolution (R ≈ 1300–2000) spectrographs: Palomar/Norris, WIYN/Hydra, Keck/DEep Imaging Multi-Object Spectrograph (DEIMOS), and the Multiple Mirror Telescope (MMT)/Hectospec. We identify 580 young stars as likely members of the NAP region based on criteria involving infrared excess, Li I 6708 Å absorption, X-ray emission, parallax, and proper motions. The spectral types of individual spectra are derived by fitting them with templates that are either empirical spectra of pre-main-sequence stars or model atmospheres. The templates are artificially veiled, and a best-fit combination of spectral type and veiling parameter is derived for each star. We use the spectral types with archival photometry to derive V-band extinction and stellar luminosity. From the Hertzsprung–Russell diagram, the median age of the young stars is about 1 Myr, with a luminosity dispersion of ~0.3–0.4 dex. We investigate the photometric variability of the spectroscopic member sample using Zwicky Transient Facility data and conclude that photometric variability, while present, does not significantly contribute to the luminosity dispersion. While larger than the formal errors, the luminosity dispersion is smaller than if veiling were not taken into account in our spectral typing process. The measured ages of the stellar kinematic groups, combined with the inferred ages for embedded stellar populations revealed by Spitzer, suggest a sequential history of star formation in the NAP region

CSI 2264: Characterizing Accretion-Burst Dominated Light Curves for Young Stars in NGC 2264

Based on more than four weeks of continuous high cadence photometric monitoring of several hundred members of the young cluster NGC 2264 with two space telescopes, NASA's Spitzer and the CNES CoRoT (Convection, Rotation, and planetary Transits), we provide high quality, multi-wavelength light curves for young stellar objects (YSOs) whose optical variability is dominated by short duration flux bursts, which we infer are due to enhanced mass accretion rates. These light curves show many brief -- several hour to one day -- brightenings at optical and near-infrared (IR) wavelengths with amplitudes generally in the range 5-50% of the quiescent value. Typically, a dozen or more of these bursts occur in a thirty day period. We demonstrate that stars exhibiting this type of variability have large ultraviolet (UV) excesses and dominate the portion of the u-g vs. g-r color-color diagram with the largest UV excesses. These stars also have large Halpha equivalent widths, and either centrally peaked, lumpy Halpha emission profiles or profiles with blue-shifted absorption dips associated with disk or stellar winds. Light curves of this type have been predicted for stars whose accretion is dominated by Rayleigh-Taylor instabilities at the boundary between their magnetosphere and inner circumstellar disk, or where magneto-rotational instabilities modulate the accretion rate from the inner disk. Amongst the stars with the largest UV excesses or largest Halpha equivalent widths, light curves with this type of variability greatly outnumber light curves with relatively smooth sinusoidal variations associated with long-lived hot spots. We provide quantitative statistics for the average duration and strength of the accretion bursts and for the fraction of the accretion luminosity associated with these bursts.Comment: Accepted for publication in AJ. 39 pages; 6 tables; 25 figures, many of which are highly degraded to meet size limits. Please download the regular resolution version at http://web.ipac.caltech.edu/staff/amc/staufferetal2014.pd

CSI 2264: Simultaneous optical and infrared light curves of young disk-bearing stars in NGC 2264 with CoRoT and Spitzer-- evidence for multiple origins of variability

Anne Marie Cody, et al, 'CSI 2264: SIMULTANEOUS OPTICAL AND INFRARED LIGHT CURVES OF YOUNG DISK-BEARING STARS IN NGC 2264 WITH CoRoT and SPITZER—EVIDENCE FOR MULTIPLE ORIGINS OF VARIABILITY', The Astronomical Journal, 147:82 (47pp), 2014 April doi:10.1088/0004-6256/147/4/82 © 2014. The American Astronomical Society.We present the Coordinated Synoptic Investigation of NGC 2264, a continuous 30-day multi-wavelength photometric monitoring campaign on more than 1000 young cluster members using 16 telescopes. The unprecedented combination of multi-wavelength, high-precision, high-cadence, and long-duration data opens a new window into the time domain behavior of young stellar objects. Here we provide an overview of the observations, focusing on results from Spitzer and CoRoT. The highlight of this work is detailed analysis of 162 classical T Tauri stars for which we can probe optical and mid-infrared flux variations to 1% amplitudes and sub-hour timescales. We present a morphological variability census and then use metrics of periodicity, stochasticity, and symmetry to statistically separate the light curves into seven distinct classes, which we suggest represent different physical processes and geometric effects. We provide distributions of the characteristic timescales and amplitudes, and assess the fractional representation within each class. The largest category (>20%) are optical "dippers" having discrete fading events lasting ~1-5 days. The degree of correlation between the optical and infrared light curves is positive but weak; notably, the independently assigned optical and infrared morphology classes tend to be different for the same object. Assessment of flux variation behavior with respect to (circum)stellar properties reveals correlations of variability parameters with H$\alpha$ emission and with effective temperature. Overall, our results point to multiple origins of young star variability, including circumstellar obscuration events, hot spots on the star and/or disk, accretion bursts, and rapid structural changes in the inner disk.Peer reviewe

New Insights from Aperiodic Variability of Young Stars

Nearly all young stars are variable, with the variability traditionally divided into two classes: periodic variables and aperiodic or "irregular" variables. Periodic variables have been studied extensively, typically using periodograms, while aperiodic variables have received much less attention due to a lack of standard statistical tools. However, aperiodic variability can serve as a powerful probe of young star accretion physics and inner circumstellar disk structure. For my dissertation, I analyzed data from a large-scale, long-term survey of the nearby North America Nebula complex, using Palomar Transient Factory photometric time series collected on a nightly or every few night cadence over several years. This survey is the most thorough exploration of variability in a sample of thousands of young stars over time baselines of days to years, revealing a rich array of lightcurve shapes, amplitudes, and timescales. I have constrained the timescale distribution of all young variables, periodic and aperiodic, on timescales from less than a day to ~100 days. I have shown that the distribution of timescales for aperiodic variables peaks at a few days, with relatively few (~15%) sources dominated by variability on tens of days or longer. My constraints on aperiodic timescale distributions are based on two new tools, magnitude- vs. time-difference (Δm-Δt) plots and peak-finding plots, for describing aperiodic lightcurves; this thesis provides simulations of their performance and presents recommendations on how to apply them to aperiodic signals in other time series data sets. In addition, I have measured the error introduced into colors or SEDs from combining photometry of variable sources taken at different epochs. These are the first quantitative results to be presented on the distributions in amplitude and time scale for young aperiodic variables, particularly those varying on timescales of weeks to months.</p