Showing 1–50 of 55 results for author: Flaherty, K

Evidence for Non-zero Turbulence in the Protoplanetary disc around IM Lup

Authors: Kevin Flaherty, A. Meredith Hughes, Jacob B. Simon, Alicia Smith Reina, Chunhua Qi, Xue-Ning Bai, Sean M. Andrews, David J. Wilner, Agnes Kospal

Abstract: The amount of turbulence in protoplanetary discs around young stars is critical for determining the efficiency, timeline, and outcomes of planet formation. It is also difficult to measure. Observations are still limited, but direct measurements of the non-thermal, turbulent gas motion are possible with the Atacama Large Millimeter/submillimeter Array (ALMA). Using CO(2-1)/$^{13}$CO(2-1)/C$^{18}$O(… ▽ More The amount of turbulence in protoplanetary discs around young stars is critical for determining the efficiency, timeline, and outcomes of planet formation. It is also difficult to measure. Observations are still limited, but direct measurements of the non-thermal, turbulent gas motion are possible with the Atacama Large Millimeter/submillimeter Array (ALMA). Using CO(2-1)/$^{13}$CO(2-1)/C$^{18}$O(2-1) ALMA observations of the disc around IM Lup at ~0.4" (~60 au) resolution we find evidence of significant turbulence, at the level of $δv_{\rm turb}=(0.18-0.30)$c$_s$. This result is robust against systematic uncertainties (e.g., amplitude flux calibration, midplane gas temperature, disc self-gravity). We find that gravito-turbulence as the source of the gas motion is unlikely based on the lack of an imprint on the rotation curve from a massive disc, while magneto-rotational instabilities and hydrodynamic instabilities are still possible, depending on the unknown magnetic field strength and the cooling timescale in the outer disc. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Comments: Accepted by MNRAS, 17 pages, 12 figures

Mapping Protoplanetary Disk Vertical Structure with CO Isotopologue Line Emission

Authors: Charles J. Law, Richard Teague, Karin I. Öberg, Evan A. Rich, Sean M. Andrews, Jaehan Bae, Myriam Benisty, Stefano Facchini, Kevin Flaherty, Andrea Isella, Sheng Jin, Jun Hashimoto, Jane Huang, Ryan A. Loomis, Feng Long, Carlos E. Muñoz-Romero, Teresa Paneque-Carreño, Laura M. Pérez, Chunhua Qi, Kamber R. Schwarz, Jochen Stadler, Takashi Tsukagoshi, David J. Wilner, Gerrit van der Plas

Abstract: High spatial resolution observations of CO isotopologue line emission in protoplanetary disks at mid-inclinations (${\approx}$30-75°) allow us to characterize the gas structure in detail, including radial and vertical substructures, emission surface heights and their dependencies on source characteristics, and disk temperature profiles. By combining observations of a suite of CO isotopologues, we… ▽ More High spatial resolution observations of CO isotopologue line emission in protoplanetary disks at mid-inclinations (${\approx}$30-75°) allow us to characterize the gas structure in detail, including radial and vertical substructures, emission surface heights and their dependencies on source characteristics, and disk temperature profiles. By combining observations of a suite of CO isotopologues, we can map the 2D (r, z) disk structure from the disk upper atmosphere, as traced by CO, to near the midplane, as probed by less abundant isotopologues. Here, we present high angular resolution (${\lesssim}$0."1 to ${\approx}$0."2; ${\approx}$15-30 au) observations of CO, $^{13}$CO, and C$^{18}$O in either or both J=2-1 and J=3-2 lines in the transition disks around DM Tau, Sz 91, LkCa 15, and HD 34282. We derived line emission surfaces in CO for all disks and in $^{13}$CO for the DM Tau and LkCa 15 disks. With these observations, we do not resolve the vertical structure of C$^{18}$O in any disk, which is instead consistent with C$^{18}$O emission originating from the midplane. Both the J=2-1 and J=3-2 lines show similar heights. Using the derived emission surfaces, we computed radial and vertical gas temperature distributions for each disk, including empirical temperature models for the DM Tau and LkCa 15 disks. After combining our sample with literature sources, we find that $^{13}$CO line emitting heights are also tentatively linked with source characteristics, e.g., stellar host mass, gas temperature, disk size, and show steeper trends than seen in CO emission surfaces. △ Less

Submitted 16 December, 2022; originally announced December 2022.

Comments: 31 pages, 15 figures, accepted for publication in ApJ. Image cubes available at https://zenodo.org/record/7430257

Millimeter Dust Emission and Planetary Dynamics in the HD 106906 System

Authors: Anna Fehr, A. Meredith Hughes, Rebekah I. Dawson, Rachel E. Marino, Matan Ackelsberg, Jamar Kittling, Kevin M. Flaherty, Erika Nesvold, John Carpenter, Sean M. Andrews, Brenda Matthews, Katie Crotts, Paul Kalas

Abstract: Debris disks are dusty, optically thin structures around main sequence stars. HD 106906AB is a short-period stellar binary, host to a wide separation planet, HD 106906b, and a debris disk. Only a few known systems include a debris disk and a directly imaged planet, and HD 106906 is the only one in which the planet is exterior to the disk. The debris disk is edge-on and highly asymmetric in scatter… ▽ More Debris disks are dusty, optically thin structures around main sequence stars. HD 106906AB is a short-period stellar binary, host to a wide separation planet, HD 106906b, and a debris disk. Only a few known systems include a debris disk and a directly imaged planet, and HD 106906 is the only one in which the planet is exterior to the disk. The debris disk is edge-on and highly asymmetric in scattered light. Here we resolve the disk structure at a resolution of 0.38" (39 au) with the Atacama Large Millimeter/submillimeter Array (ALMA) at a wavelength of 1.3 mm. We model the disk with both a narrow and broad ring of material, and find that a radially broad, axisymmetric disk between radii of $\sim$50$-$100 au is able to capture the structure of the observations without evidence of any asymmetry or eccentricity, other than a tentative stellocentric offset. We place stringent upper limits on both the gas and dust content of a putative circumplanetary disk. We interpret the ALMA data in concert with scattered light observations of the inner ring and astrometric constraints on the planet's orbit, and find that the observations are consistent with a large-separation, low-eccentricity orbit for the planet. A dynamical analysis indicates that the central binary can efficiently stabilize planetesimal orbits interior to $\sim$100 au, which relaxes the constraints on eccentricity and semimajor axis somewhat. The observational constraints are consistent with in situ formation via gravitational instability, but cannot rule out a scattering event as the origin for HD 106906b's current orbit. △ Less

Submitted 23 September, 2022; originally announced September 2022.

Multiwavelength Vertical Structure in the AU Mic Debris Disk: Characterizing the Collisional Cascade

Authors: David Vizgan, A. Meredith Hughes, Evan S. Carter, Kevin M. Flaherty, Margaret Pan, Eugene Chiang, Hilke Schlichting, David J. Wilner, Sean M. Andrews, John M. Carpenter, Atilla Moór, Meredith A. MacGregor

Abstract: Debris disks are scaled-up analogs of the Kuiper Belt in which dust is generated by collisions between planetesimals. In the "collisional cascade" model of debris disks, dust lost to radiation pressure and winds is constantly replenished by grinding collisions between planetesimals. The model assumes that collisions are destructive and involve large velocities; this assumption has not been tested… ▽ More Debris disks are scaled-up analogs of the Kuiper Belt in which dust is generated by collisions between planetesimals. In the "collisional cascade" model of debris disks, dust lost to radiation pressure and winds is constantly replenished by grinding collisions between planetesimals. The model assumes that collisions are destructive and involve large velocities; this assumption has not been tested beyond our Solar System. We present 0"25 ($\approx$2.4 au) resolution observations of the $λ$ = 450 $μ$m dust continuum emission from the debris disk around the nearby M dwarf AU Microscopii with the Atacama Large Millimeter/submillimeter Array. We use parametric models to describe the disk structure, and an MCMC algorithm to explore the posterior distributions of the model parameters; we fit the structure of the disk to both our data and archival $λ= 1.3$ mm data (Daley et al. 2019), from which we obtain two aspect ratio measurements at 1.3 mm ($h_{1300}$ = 0.025$^{+0.008}_{-0.002}$) and at 450 $μ$m ($h_{450}$ = 0.019$^{+0.006}_{-0.001}$), as well as the grain size distribution index $q =$ 3.03 $\pm$ 0.02. Contextualizing our aspect ratio measurements within the modeling framework laid out in Pan & Schlichting (2012), we derive a power law index of velocity dispersion as a function of grain size $p = 0.28 \pm 0.06$ for the AU Mic debris disk. This result implies that smaller bodies are more easily disrupted than larger bodies by collisions, which is inconsistent with the strength regime usually assumed for such small bodies. Possible explanations for this discrepancy are discussed. △ Less

Submitted 11 July, 2022; originally announced July 2022.

Comments: 17 pages, 6 figures, 3 tables. Accepted for publication in ApJ

CO Line Emission Surfaces and Vertical Structure in Mid-Inclination Protoplanetary Disks

Authors: Charles J. Law, Sage Crystian, Richard Teague, Karin I. Öberg, Evan A. Rich, Sean M. Andrews, Jaehan Bae, Kevin Flaherty, Viviana V. Guzmán, Jane Huang, John D. Ilee, Joel H. Kastner, Ryan A. Loomis, Feng Long, Laura M. Pérez, Sebastián Pérez, Chunhua Qi, Giovanni P. Rosotti, Dary Ruíz-Rodríguez, Takashi Tsukagoshi, David J. Wilner

Abstract: High spatial resolution CO observations of mid-inclination (30-75°) protoplanetary disks offer an opportunity to study the vertical distribution of CO emission and temperature. The asymmetry of line emission relative to the disk major axis allows for a direct mapping of the emission height above the midplane, and for optically-thick, spatially-resolved emission in LTE, the intensity is a measure o… ▽ More High spatial resolution CO observations of mid-inclination (30-75°) protoplanetary disks offer an opportunity to study the vertical distribution of CO emission and temperature. The asymmetry of line emission relative to the disk major axis allows for a direct mapping of the emission height above the midplane, and for optically-thick, spatially-resolved emission in LTE, the intensity is a measure of the local gas temperature. Our analysis of ALMA archival data yields CO emission surfaces, dynamically-constrained stellar host masses, and disk atmosphere gas temperatures for the disks around: HD 142666, MY Lup, V4046 Sgr, HD 100546, GW Lup, WaOph 6, DoAr 25, Sz 91, CI Tau, and DM Tau. These sources span a wide range in stellar masses (0.50-2.10 M$_{\odot}$), ages (${\sim}$0.3-23 Myr), and CO gas radial emission extents (${\approx}$200-1000 au). This sample nearly triples the number of disks with mapped emission surfaces and confirms the wide diversity in line emitting heights ($z/r\approx0.1$ to ${\gtrsim}0.5$) hinted at in previous studies. We compute radial and vertical CO gas temperature distributions for each disk. A few disks show local temperature dips or enhancements, some of which correspond to dust substructures or the proposed locations of embedded planets. Several emission surfaces also show vertical substructures, which all align with rings and gaps in the millimeter dust. Combining our sample with literature sources, we find that CO line emitting heights weakly decline with stellar mass and gas temperature, which, despite large scatter, is consistent with simple scaling relations. We also observe a correlation between CO emission height and disk size, which is due to the flared structure of disks. Overall, CO emission surfaces trace ${\approx}2$-$5\times$ gas pressure scale heights (H$_{\rm{g}}$) and could potentially be calibrated as empirical tracers of H$_{\rm{g}}$. △ Less

Submitted 3 May, 2022; originally announced May 2022.

Comments: 31 pages, 17 figures, accepted for publication in ApJ. Image cubes available at https://zenodo.org/record/6410045

Kinematic Structures in Planet-Forming Disks

Authors: Christophe Pinte, Richard Teague, Kevin Flaherty, Cassandra Hall, Stefano Facchini, Simon Casassus

Abstract: The past 5 years have dramatically changed our view of the disks of gas and dust around young stars. Observations with the Atacama Large Millimeter/submillimeter Array (ALMA) and extreme adaptive optics systems have revealed that disks are dynamical systems. Most disks contain resolved structures, both in gas and dust, including rings, gaps, spirals, azimuthal dust concentrations, shadows cast by… ▽ More The past 5 years have dramatically changed our view of the disks of gas and dust around young stars. Observations with the Atacama Large Millimeter/submillimeter Array (ALMA) and extreme adaptive optics systems have revealed that disks are dynamical systems. Most disks contain resolved structures, both in gas and dust, including rings, gaps, spirals, azimuthal dust concentrations, shadows cast by misaligned inner disks, as well as deviations from Keplerian rotation. The origin of these structures and how they relate to the planet formation process remain poorly understood. Spatially resolved kinematic studies offer a new and necessary window to understand and quantify the physical processes (turbulence, winds, radial and meridional flows, stellar multiplicity, instabilities) at play during planet formation and disk evolution. Recent progress, driven mainly by resolved ALMA observations, includes the detection and mass determination of embedded planets, the mapping of the gas flow around the accreting planets, the confirmation of tidal interactions and warped disk geometries, and stringent limits on the turbulent velocities. In this chapter, we will review our current understanding of these dynamical processes and highlight how kinematic mapping provides new ways to observe planet formation in action. △ Less

Submitted 17 March, 2022; originally announced March 2022.

Comments: Protostars and Planets VII chapter, 35 pages, 14 figures. Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura

A Deep Search for Five Molecules in the 49 Ceti Debris Disk

Authors: Jessica Klusmeyer, A. Meredith Hughes, Luca Matra, Kevin Flaherty, Agnes Kospal, Attila Moor, Aki Roberge, Karin Oberg, Aaron Boley, Jacob White, David Wilner, Peter Abraham

Abstract: Surprisingly strong CO emission has been observed from more than a dozen debris disks around nearby main-sequence stars. The origin of this CO is unclear, in particular whether it is left over from the protoplanetary disk phase or is second-generation material released from collisions between icy bodies like debris dust. The primary unexplored avenue for distinguishing the origin of the material i… ▽ More Surprisingly strong CO emission has been observed from more than a dozen debris disks around nearby main-sequence stars. The origin of this CO is unclear, in particular whether it is left over from the protoplanetary disk phase or is second-generation material released from collisions between icy bodies like debris dust. The primary unexplored avenue for distinguishing the origin of the material is understanding its molecular composition. Here we present a deep search for five molecules (CN, HCN, HCO+, SiO, and CH3OH) in the debris disk around 49 Ceti. We take advantage of the high sensitivity of the Atacama Large Millimeter/submillimeter Array (ALMA) at Band 7 to integrate for 3.2 hours at modest spatial (1") and spectral (0.8 km/s) resolution. Our search yields stringent upper limits on the flux of all surveyed molecular lines, which imply abundances relative to CO that are orders of magnitude lower than those observed in protoplanetary disks and Solar System comets, and also those predicted in outgassing models of second-generation material. However, if CI shielding is responsible for extending the lifetime of any CO produced in second-generation collisions, as proposed by Kral et al. (2018), then the line ratios do not reflect true ice phase chemical abundances, but rather imply that CO is shielded by its own photodissociation product, CI, but other molecules are rapidly photodissociated by the stellar and interstellar radiation field. △ Less

Submitted 15 July, 2021; originally announced July 2021.

Comments: 26 pages, 4 figures, 4 tables, accepted for publication in ApJ

Resolving Structure in the Debris Disk around HD 206893 with ALMA

Authors: Ava Nederlander, A. Meredith Hughes, Anna J. Fehr, Kevin M. Flaherty, Kate Y. L. Su, Attila Moor, Eugene Chiang, Sean M. Andrews, David J. Wilner, Sebastian Marino

Abstract: Debris disks are tenuous, dusty belts surrounding main sequence stars generated by collisions between planetesimals. HD 206893 is one of only two stars known to host a directly imaged brown dwarf orbiting interior to its debris ring, in this case at a projected separation of 10.4 au. Here we resolve structure in the debris disk around HD 206893 at an angular resolution of 0.6" (24 au) and waveleng… ▽ More Debris disks are tenuous, dusty belts surrounding main sequence stars generated by collisions between planetesimals. HD 206893 is one of only two stars known to host a directly imaged brown dwarf orbiting interior to its debris ring, in this case at a projected separation of 10.4 au. Here we resolve structure in the debris disk around HD 206893 at an angular resolution of 0.6" (24 au) and wavelength of 1.3 mm with the Atacama Large Millimeter/submillimeter Array (ALMA). We observe a broad disk extending from a radius of <51 au to 194^{+13}_{-2} au. We model the disk with a continuous, gapped, and double power-law model of the surface density profile, and find strong evidence for a local minimum in the surface density distribution near a radius of 70 au, consistent with a gap in the disk with an inner radius of 63^{+8}_{-16} au and width 31^{+11}_{-7} au. Gapped structure has been observed in four other debris disks -- essentially every other radially resolved debris disk observed with sufficient angular resolution and sensitivity with ALMA -- and could be suggestive of the presence of an additional planetary-mass companion. △ Less

Submitted 21 January, 2021; originally announced January 2021.

Comments: 24 pages, 12 figures, Accepted for publication in ApJ

Measuring turbulent motion in planet-forming disks with ALMA: A detection around DM Tau and non-detections around MWC 480 and V4046 Sgr

Authors: Kevin Flaherty, A. Meredith Hughes, Jacob B. Simon, Chunhua Qi, Xue-ning Bai, Alyssa Bulatek, Sean M. Andrews, David J. Wilner, Agnes Kospal

Abstract: Turbulence is a crucial factor in many models of planet formation, but it has only been directly constrained among a small number of planet forming disks. Building on the upper limits on turbulence placed in disks around HD 163296 and TW Hya, we present ALMA CO J=2-1 line observations at $\sim0.3"$ (20-50 au) resolution and 80 ms$^{-1}$ channel spacing of the disks around DM Tau, MWC 480, and V404… ▽ More Turbulence is a crucial factor in many models of planet formation, but it has only been directly constrained among a small number of planet forming disks. Building on the upper limits on turbulence placed in disks around HD 163296 and TW Hya, we present ALMA CO J=2-1 line observations at $\sim0.3"$ (20-50 au) resolution and 80 ms$^{-1}$ channel spacing of the disks around DM Tau, MWC 480, and V4046 Sgr. Using parametric models of disk structure, we robustly detect non-thermal gas motions around DM Tau of between 0.25 c$_s$ and 0.33 c$_s$, with the range dominated by systematic effects, making this one of the only systems with directly measured non-zero turbulence. Using the same methodology, we place stringent upper limits on the non-thermal gas motion around MWC 480 ($<$0.08c$_s$) and V4046 Sgr ($<$0.12 c$_s$). The preponderance of upper limits in this small sample, and the modest turbulence levels consistent with dust studies, suggest that weak turbulence ($α\lesssim10^{-3}$) may be a common, albeit not universal, feature of planet-forming disks. We explore the particular physical conditions around DM Tau that could lead this system to be more turbulent than the others. △ Less

Submitted 28 April, 2020; v1 submitted 25 April, 2020; originally announced April 2020.

Comments: 21 pages, 10 figures, accepted to ApJ

From scattered-light to millimeter emission: A comprehensive view of the Gyr-old system of HD 202628 and its eccentric debris ring

Authors: Virginie Faramaz, John Krist, Karl R. Stapelfeldt, Geoffrey Bryden, Eric E. Mamajek, Luca Matrà, Mark Booth, Kevin Flaherty, Antonio S. Hales, A. Meredith Hughes, Amelia Bayo, Simon Casassus, Jorge Cuadra, Johan Olofsson, Kate Y. L. Su, David J. Wilner

Abstract: We present here new observations of the eccentric debris ring surrounding the Gyr-old solar-type star HD 202628: at millimeter wavelengths with ALMA, at far-infrared wavelengths with \textit{Herschel}, and in scattered light with \textit{HST}. The ring inner edge is found to be consistent between ALMA and \textit{HST} data. As radiation pressure affects small grains seen in scattered-light, the ri… ▽ More We present here new observations of the eccentric debris ring surrounding the Gyr-old solar-type star HD 202628: at millimeter wavelengths with ALMA, at far-infrared wavelengths with \textit{Herschel}, and in scattered light with \textit{HST}. The ring inner edge is found to be consistent between ALMA and \textit{HST} data. As radiation pressure affects small grains seen in scattered-light, the ring appears broader at optical than at millimeter wavelengths. The best fit to the ring seen with ALMA has inner and outer edges at $143.1 \pm 1.7$ AU and $165.5 \pm 1.4$, respectively, and an inclination of $57.4^\circ \pm 0.4$ from face-on. The offset of the ring centre of symmetry from the star allows us to quantify its eccentricity to be $e=0.09_{-0.01}^{+0.02}$. This eccentric feature is also detected in low resolution \textit{Herschel}/PACS observations, under the form of a pericenter-glow. Combining the infrared and millimeter photometry, we retrieve a disk grain size distribution index of $\sim -3.4$, and therefore exclude in-situ formation of the inferred belt-shaping perturber, for which we provide new dynamical constraints. Finally, ALMA images show four point-like sources that exceed 100$\,μ$Jy, one of them being just interior to the ring. Although the presence of a background object cannot be excluded, we cannot exclude either that this source is circumplanetary material surrounding the belt-shaper, in which case degeneracies between its mass and orbital parameters could be lifted, allowing us to fully characterize such a distant planet in this mass and age regime for the very first time. △ Less

Submitted 9 September, 2019; originally announced September 2019.

Comments: 29 pages, 15 figures

Modelling the Spatial Distribution and Origin of CO Gas in Debris Disks

Authors: Antonio Hales, Uma Gorti, John Carpenter, Meredith Hughes, Kevin Flaherty

Abstract: The detection of gas in debris disks raises the question of whether this gas is a remnant from the primordial protoplanetary phase, or released by the collision of secondary bodies. In this paper we analyze ALMA observations at 1-1.5" resolution of three debris disks where the $^{12}$CO(2-1) rotational line was detected: HD131835, HD138813, and HD156623. We apply the iterative Lucy-Richardson deco… ▽ More The detection of gas in debris disks raises the question of whether this gas is a remnant from the primordial protoplanetary phase, or released by the collision of secondary bodies. In this paper we analyze ALMA observations at 1-1.5" resolution of three debris disks where the $^{12}$CO(2-1) rotational line was detected: HD131835, HD138813, and HD156623. We apply the iterative Lucy-Richardson deconvolution technique to the problem of circumstellar disks to derive disk geometries and surface brightness distributions of the gas. The derived disk parameters are used as input for thermochemical models to test both primordial and cometary scenarios for the origin of the gas. We favor a secondary origin for the gas in these disks and find that the CO gas masses ($\sim 3\times10^{-3}$ M$_{\oplus}$) require production rates ($\sim 5\times 10^{-7}$ M$_{\oplus}$~yr$^{-1}$) similar to those estimated for the bona-fide gas rich debris disk $β$ Pic. △ Less

Submitted 9 May, 2019; originally announced May 2019.

Comments: accepted for publication in ApJ

The Inner Disk and Accretion Flow of the Close Binary DQ Tau

Authors: James Muzerolle, Kevin Flaherty, Zoltan Balog, Tracy Beck, Robert Gutermuth

Abstract: We present multi-epoch optical and near-infrared (NIR) photometry and spectroscopy of the spectroscopic binary T Tauri star DQ Tau. The photometric monitoring, obtained using SMARTS ANDICAM, recovers the previously-seen correlation between optical flux and the 15.8-day binary orbital period, with blue flux peaks occurring close to most observed periastron passages. For the first time, we find an e… ▽ More We present multi-epoch optical and near-infrared (NIR) photometry and spectroscopy of the spectroscopic binary T Tauri star DQ Tau. The photometric monitoring, obtained using SMARTS ANDICAM, recovers the previously-seen correlation between optical flux and the 15.8-day binary orbital period, with blue flux peaks occurring close to most observed periastron passages. For the first time, we find an even more consistent correlation between orbital period and NIR brightness and color. The onset of pulse events in the NIR on average precedes those in the optical by a few days, with the rise usually starting near apastron orbital phase. We further obtained five epochs of spectroscopy using IRTF SpeX, with a wavelength range of 0.8 to 5 microns, and derived spectra of the infrared excess emission. The shape and strength of the excess varies with time, with cooler and weaker characteristic dust emission (T ~ 1100-1300 K) over most of the binary orbit, and stronger/warmer dust emission (T ~ 1600 K, indicative of dust sublimation) just before periastron passage. We suggest our results are broadly consistent with predictions of simulations of disk structure and accretion flows around close binaries, with the varying dust emission possibly tracing the evolution of accretion streams falling inwards through a circumbinary disk cavity and feeding the accretion pulses traced by the optical photometry and NIR emission lines. However, our results also show more complicated behavior that is not fully explained by this simple picture, and will require further observations and modeling to fully interpret. △ Less

Submitted 12 April, 2019; originally announced April 2019.

Comments: accepted to ApJ

The Mass of Stirring Bodies in the AU Mic Debris Disk Inferred from Resolved Vertical Structure

Authors: Cail Daley, A. Meredith Hughes, Evan S. Carter, Kevin Flaherty, Zachary Lambros, Margaret Pan, Hilke Schlichting, Eugene Chiang, Mark Wyatt, David Wilner, Sean Andrews, John Carpenter

Abstract: The vertical distribution of dust in debris disks is sensitive to the number and size of large planetesimals dynamically stirring the disk, and is therefore well-suited for constraining the prevalence of otherwise unobservable Uranus and Neptune analogs. Information regarding stirring bodies has previously been inferred from infrared and optical observations of debris disk vertical structure, but… ▽ More The vertical distribution of dust in debris disks is sensitive to the number and size of large planetesimals dynamically stirring the disk, and is therefore well-suited for constraining the prevalence of otherwise unobservable Uranus and Neptune analogs. Information regarding stirring bodies has previously been inferred from infrared and optical observations of debris disk vertical structure, but theoretical works predict that the small particles traced by short-wavelength observations will be `puffed up' by radiation pressure, yielding only upper limits. The large grains that dominate the disk emission at millimeter wavelengths are much less sensitive to the effects of stellar radiation or stellar winds, and therefore trace the underlying mass distribution more directly. Here we present ALMA 1.3 mm dust continuum observations of the debris disk around the nearby M star AU Mic. The 3 au spatial resolution of the observations, combined with the favorable edge-on geometry of the system, allows us to measure the vertical thickness of the disk. We report a scale height-to-radius aspect ratio of $h = 0.031_{-0.004}^{+0.005}$ between radii of $\sim 23$ au and $\sim 41$ au. Comparing this aspect ratio to a theoretical model of size-dependent velocity distributions in the collisional cascade, we find that the perturbing bodies embedded in the local disk must be larger than about 400 km, and the largest perturbing body must be smaller than roughly 1.8 M$_\odot$. These measurements rule out the presence of a gas giant or Neptune analog near the $\sim 40$ au outer edge of the debris ring, but are suggestive of large planetesimals or an Earth-sized planet stirring the dust distribution. △ Less

Submitted 8 April, 2019; v1 submitted 29 March, 2019; originally announced April 2019.

Comments: 21 pages, 8 figures, 5 tables. Accepted for publication in ApJ

Realizing the Unique Potential of ALMA to Probe the Gas Reservoir of Planet Formation

Authors: L. Ilsedore Cleeves, Ryan Loomis, Richard Teague, Ke Zhang, Edwin Bergin, Karin Oberg, Crystal Brogan, Todd Hunter, Yuri Aikawa, Sean Andrews, Jaehan Bae, Jennifer Bergner, Kevin Flaherty, Viviana Guzman, Jane Huang, Michiel Hogerheijde, Shih-Ping Lai, Laura Perez, Luca Ricci, Colette Salyk, Kamber Schwarz, Jonathan Williams, David Wilner, Al Wootten

Abstract: Understanding the origin of the astonishing diversity of exoplanets is a key question for the coming decades. ALMA has revolutionized our view of the dust emission from protoplanetary disks, demonstrating the prevalence of ring and spiral structures that are likely sculpted by young planets in formation. To detect kinematic signatures of these protoplanets and to probe the chemistry of their gas a… ▽ More Understanding the origin of the astonishing diversity of exoplanets is a key question for the coming decades. ALMA has revolutionized our view of the dust emission from protoplanetary disks, demonstrating the prevalence of ring and spiral structures that are likely sculpted by young planets in formation. To detect kinematic signatures of these protoplanets and to probe the chemistry of their gas accretion reservoir will require the imaging of molecular spectral line emission at high angular and spectral resolution. However, the current sensitivity of ALMA limits these important spectral studies to only the nearest protoplanetary disks. Although some promising results are emerging, including the identification of the snowlines of a few key molecules and the first attempt at detecting a protoplanet's spiral wake, it is not yet possible to search for these important signatures in a population of disks in diverse environments and ages. Harnessing the tremendous power of (sub)mm observations to pinpoint and characterize the chemistry of planets in formation will require a major increase of ALMA's spectral sensitivity (5-10x), increase in instantaneous bandwidth (2x) at high spectral resolution, and improved angular resolution (2x) in the 2030 era. △ Less

Submitted 27 March, 2019; originally announced March 2019.

Comments: Astro2020 White Paper

The Planet Formation Potential Around a 45 Myr old Accreting M Dwarf

Authors: Kevin M. Flaherty, A. Meredith Hughes, Eric E. Mamajek, Simon J. Murphy

Abstract: Debris disk detections around M dwarfs are rare, and so far no gas emission has been detected from an M dwarf debris disk. This makes the 45 Myr old M dwarf WISEJ080822.18-644357.3 a bit of a curiosity; it has a strong infrared excess at an age beyond the lifetime of a typical planet-forming disk, and also exhibits broad H$α$ emission consistent with active accretion from a gaseous disk. To better… ▽ More Debris disk detections around M dwarfs are rare, and so far no gas emission has been detected from an M dwarf debris disk. This makes the 45 Myr old M dwarf WISEJ080822.18-644357.3 a bit of a curiosity; it has a strong infrared excess at an age beyond the lifetime of a typical planet-forming disk, and also exhibits broad H$α$ emission consistent with active accretion from a gaseous disk. To better understand the cold gas and dust properties of this system, we obtained ALMA observations of the 1.3mm continuum and the CO/$^{13}$CO/C$^{18}$O J=2-1 emission lines. No cold CO gas is detected from this system, ruling out a gas-rich protoplanetary disk. Unresolved dust continuum emission is detected at a flux of 198$\pm$15 $μ$Jy, consistent with 0.057$\pm$0.006 M$_{\oplus}$ worth of optically thin dust, and consistent with being generated through a collisional cascade induced by large bodies at radii $<$16 au. With a sufficiently strong stellar wind, dust grains released in the outer disk can migrate inwards via PR drag, potentially serving as a source of grains for the strong infrared excess. △ Less

Submitted 10 December, 2018; originally announced December 2018.

Comments: 9 pages, 4 figures, accepted to ApJ

The Leaky Pipeline for Postdocs: A study of the time between receiving a PhD and securing a faculty job for male and female astronomers

Authors: Kevin Flaherty

Abstract: The transition between receiving a PhD and securing a tenure track faculty position is challenging for nearly every astronomer interested in working in academia. Here we use a publicly available database of recently hired faculty (the Astrophysics Job Rumor Mill) to examine the amount of time astronomers typically spend in this transitory state. Using these data as a starting point to examine the… ▽ More The transition between receiving a PhD and securing a tenure track faculty position is challenging for nearly every astronomer interested in working in academia. Here we use a publicly available database of recently hired faculty (the Astrophysics Job Rumor Mill) to examine the amount of time astronomers typically spend in this transitory state. Using these data as a starting point to examine the experiences of astronomy postdocs, we find that the average time spent between receiving a PhD and being hired into a faculty position is 4.9$\pm$0.3 years, with female astronomers hired on average 4.2$\pm$0.4 years after receiving a PhD while male astronomers are typically hired after 5.3$\pm$0.4 years. Using a simple model of the labor market, we attempt to recreate this gendered difference in time spent as a postdoc. We can rule out the role of the increasing representation of women among astronomy PhDs, as well as any bias in favor of hiring female astronomers in response to efforts to diversify the faculty ranks. Instead the most likely explanation is that female astronomers are leaving the academic labor market, at a rate that is 3-4 times higher than male astronomers. This scenario explains the distinct hiring time distributions between male and female astronomers, as well as the measured percentage of female assistant professors, and the fraction of female applicants within a typical faculty search. These results provide evidence that more work needs to be done to support and retain female astronomers during the postdoctoral phase of their careers. △ Less

Submitted 2 October, 2018; originally announced October 2018.

Comments: 16 pages, 7 figures, comment welcome

Turbulence in the TW Hya Disk

Authors: Kevin M. Flaherty, A. Meredith Hughes, Richard Teague, Jacob B. Simon, Sean M. Andrews, David J. Wilner

Abstract: Turbulence is a fundamental parameter in models of grain growth during the early stages of planet formation. As such, observational constraints on its magnitude are crucial. Here we self-consistently analyze ALMA CO(2-1), SMA CO(3-2), and SMA CO(6-5) observations of the disk around TW Hya and find an upper limit on the turbulent broadening of $<$0.08c$_s$ ($α<$0.007 for $α$ defined only within 2-3… ▽ More Turbulence is a fundamental parameter in models of grain growth during the early stages of planet formation. As such, observational constraints on its magnitude are crucial. Here we self-consistently analyze ALMA CO(2-1), SMA CO(3-2), and SMA CO(6-5) observations of the disk around TW Hya and find an upper limit on the turbulent broadening of $<$0.08c$_s$ ($α<$0.007 for $α$ defined only within 2-3 pressure scale heights above the midplane), lower than the tentative detection previously found from an analysis of the CO(2-1) data. We examine in detail the challenges of image plane fitting vs directly fitting the visibilities, while also considering the role of the vertical temperature gradient, systematic uncertainty in the amplitude calibration, and assumptions about the CO abundance, as potential sources of the discrepancy in the turbulence measurements. These tests result in variations of the turbulence limit between $<$0.04c$_s$ and $<$0.13c$_s$, consistently lower than the 0.2-0.4c$_s$ found previously. Having ruled out numerous factors, we restrict the source of the discrepancy to our assumed coupling between temperature and density through hydrostatic equilibrium in the presence of a vertical temperature gradient and/or the confinement of CO to a thin molecular layer above the midplane, although further work is needed to quantify the influence of these prescriptions. Assumptions about hydrostatic equilibrium and the CO distribution are physically motivated, and may have a small influence on measuring the kinematics of the gas, but they become important when constraining small effects such as the strength of the turbulence within a protoplanetary disk. △ Less

Submitted 10 March, 2018; originally announced March 2018.

Comments: 14 pages, 4 figures, accepted to ApJ

Origin of Weak Turbulence in the Outer Regions of Protoplanetary Disks

Authors: Jacob B. Simon, Xue-Ning Bai, Kevin M. Flaherty, A. Meredith Hughes

Abstract: The mechanism behind angular momentum transport in protoplanetary disks, and whether this transport is turbulent in nature, is a fundamental issue in planet formation studies. Recent ALMA observations have suggested that turbulent velocities in the outer regions of these disks are less than ~5-10% of the sound speed, contradicting theoretical predictions of turbulence driven by the magnetorotation… ▽ More The mechanism behind angular momentum transport in protoplanetary disks, and whether this transport is turbulent in nature, is a fundamental issue in planet formation studies. Recent ALMA observations have suggested that turbulent velocities in the outer regions of these disks are less than ~5-10% of the sound speed, contradicting theoretical predictions of turbulence driven by the magnetorotational instability (MRI). These observations have generally been interpreted to be consistent with a large-scale laminar magnetic wind driving accretion. Here, we carry out local, shearing box simulations with varying ionization levels and background magnetic field strengths in order to determine which parameters produce results consistent with observations. We find that even when the background magnetic field launches a strong largely laminar wind, significant turbulence persists and is driven by localized regions of vertical magnetic field (the result of zonal flows) that are unstable to the MRI. The only conditions for which we find turbulent velocities below the observational limits are weak background magnetic fields and ionization levels well below that usually assumed in theoretical studies. We interpret these findings within the context of a preliminary model in which a large scale magnetic field, confined to the inner disk, hinders ionizing sources from reaching large radial distances, e.g., through a sufficiently dense wind. Thus, in addition to such a wind, this model predicts that for disks with weakly turbulent outer regions, the outer disk will have significantly reduced ionization levels compared to standard models and will harbor only a weak vertical magnetic field. △ Less

Submitted 6 August, 2018; v1 submitted 13 November, 2017; originally announced November 2017.

Comments: Accepted to ApJ, 13 pages, 6 figures, modified from original version to include more extensive discussion and analysis and a title change

ALMA 1.3 Millimeter Map of the HD 95086 System

Authors: Kate Y. L. Su, Meredith A. Macgregor, Mark Booth, David J. Wilner, Kevin Flaherty, A. Meredith Hughes, Neil M. Phillips, Renu Malhotra, Antonio S. Hales, Sarah Morrison, Steve Ertel, Brenda C. Matthews, William R. F. Dent, Simon Casassus

Abstract: Planets and minor bodies such as asteroids, Kuiper-belt objects and comets are integral components of a planetary system. Interactions among them leave clues about the formation process of a planetary system. The signature of such interactions is most prominent through observations of its debris disk at millimeter wavelengths where emission is dominated by the population of large grains that stay… ▽ More Planets and minor bodies such as asteroids, Kuiper-belt objects and comets are integral components of a planetary system. Interactions among them leave clues about the formation process of a planetary system. The signature of such interactions is most prominent through observations of its debris disk at millimeter wavelengths where emission is dominated by the population of large grains that stay close to their parent bodies. Here we present ALMA 1.3 mm observations of HD 95086, a young early-type star that hosts a directly imaged giant planet b and a massive debris disk with both asteroid- and Kuiper-belt analogs. The location of the Kuiper-belt analog is resolved for the first time. The system can be depicted as a broad ($ΔR/R \sim$0.84), inclined (30\arcdeg$\pm$3\arcdeg) ring with millimeter emission peaked at 200$\pm$6 au from the star. The 1.3 mm disk emission is consistent with a broad disk with sharp boundaries from 106$\pm$6 to 320$\pm$20 au with a surface density distribution described by a power law with an index of --0.5$\pm$0.2. Our deep ALMA map also reveals a bright source located near the edge of the ring, whose brightness at 1.3 mm and potential spectral energy distribution are consistent with it being a luminous star-forming galaxy at high redshift. We set constraints on the orbital properties of planet b assuming co-planarity with the observed disk. △ Less

Submitted 28 September, 2017; originally announced September 2017.

Comments: accepted for publication in AJ

A Three-Dimensional View of Turbulence: Constraints on Turbulent Motions in the HD 163296 Protoplanetary Disk using DCO$^+$

Authors: Kevin M. Flaherty, A. Meredith Hughes, Sanaea C. Rose, Jacob B. Simon, Chunhua Qi, Sean M. Andrews, Agnes Kospal, David J. Wilner, Eugene Chiang, Philip J. Armitage, Xue-ning Bai

Abstract: Gas kinematics are an important part of the planet formation process. Turbulence influences planetesimal growth and migration from the scale of sub-micron dust grains through gas-giant planets. Radio observations of resolved molecular line emission can directly measure this non-thermal motion and, taking advantage of the layered chemical structure of disks, different molecular lines can be combine… ▽ More Gas kinematics are an important part of the planet formation process. Turbulence influences planetesimal growth and migration from the scale of sub-micron dust grains through gas-giant planets. Radio observations of resolved molecular line emission can directly measure this non-thermal motion and, taking advantage of the layered chemical structure of disks, different molecular lines can be combined to map the turbulence throughout the vertical extent of a protoplanetary disk. Here we present ALMA observations of three molecules (DCO$^+$(3-2), C$^{18}$O(2-1) and CO(2-1)) from the disk around HD 163296. We are able to place stringent upper limits ($v_{\rm turb}<$0.06c$_s$, $<$0.05c$_s$ and $<$0.04c$_s$ for CO(2-1), C$^{18}$O(2-1) and DCO$^+$(3-2) respectively), corresponding to $α\lesssim$3$\times$10$^{-3}$, similar to our prior limit derived from CO(3-2). This indicates that there is little turbulence throughout the vertical extent of the disk, contrary to theoretical predictions based on the magneto-rotational instability and gravito-turbulence. In modeling the DCO$^+$ emission we also find that it is confined to three concentric rings at 65.7$\pm$0.9 au, 149.9$^{+0.5}_{-0.7}$ au and 259$\pm$1 au, indicative of a complex chemical environment. △ Less

Submitted 14 June, 2017; originally announced June 2017.

Comments: Accepted to ApJ

Radial Surface Density Profiles of Gas and Dust in the Debris Disk around 49 Ceti

Authors: A. M. Hughes, J. Lieman-Sifry, K. M. Flaherty, C. M. Daley, A. Roberge, A. Kospal, Attila Moor, Inga Kamp, D. J. Wilner, S. M. Andrews, J. H. Kastner, P. Abraham

Abstract: We present ~0.4 resolution images of CO(3-2) and associated continuum emission from the gas-bearing debris disk around the nearby A star 49 Ceti, observed with the Atacama Large Millimeter/Submillimeter Array (ALMA). We analyze the ALMA visibilities in tandem with the broad-band spectral energy distribution to measure the radial surface density profiles of dust and gas emission from the system. Th… ▽ More We present ~0.4 resolution images of CO(3-2) and associated continuum emission from the gas-bearing debris disk around the nearby A star 49 Ceti, observed with the Atacama Large Millimeter/Submillimeter Array (ALMA). We analyze the ALMA visibilities in tandem with the broad-band spectral energy distribution to measure the radial surface density profiles of dust and gas emission from the system. The dust surface density decreases with radius between ~100 and 310 au, with a marginally significant enhancement of surface density at a radius of ~110 au. The SED requires an inner disk of small grains in addition to the outer disk of larger grains resolved by ALMA. The gas disk exhibits a surface density profile that increases with radius, contrary to most previous spatially resolved observations of circumstellar gas disks. While ~80% of the CO flux is well described by an axisymmetric power-law disk in Keplerian rotation about the central star, residuals at ~20% of the peak flux exhibit a departure from axisymmetry suggestive of spiral arms or a warp in the gas disk. The radial extent of the gas disk (~220 au) is smaller than that of the dust disk (~300 au), consistent with recent observations of other gas-bearing debris disks. While there are so far only three broad debris disks with well characterized radial dust profiles at millimeter wavelengths, 49 Ceti's disk shows a markedly different structure from two radially resolved gas-poor debris disks, implying that the physical processes generating and sculpting the gas and dust are fundamentally different. △ Less

Submitted 6 April, 2017; originally announced April 2017.

Comments: 20 pages, 8 figures, accepted for publication in ApJ March 31, 2017 (submitted Nov 2016)

ALMA Observations of Asymmetric Molecular Gas Emission from a Protoplanetary Disk in the Orion Nebula

Authors: Samuel M. Factor, A. M. Hughes, Kevin M. Flaherty, Rita K. Mann, James Di Francesco, Jonathan P. Williams, Luca Ricci, Brenda C. Matthews, John Bally, Doug Johnstone

Abstract: We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of molecular line emission from d216-0939, one of the largest and most massive protoplanetary disks in the Orion Nebula Cluster (ONC). We model the spectrally resolved HCO$^+$ (4--3), CO (3--2), and HCN (4--3) lines observed at 0\farcs5 resolution to fit the temperature and density structure of the disk. We also weakly det… ▽ More We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of molecular line emission from d216-0939, one of the largest and most massive protoplanetary disks in the Orion Nebula Cluster (ONC). We model the spectrally resolved HCO$^+$ (4--3), CO (3--2), and HCN (4--3) lines observed at 0\farcs5 resolution to fit the temperature and density structure of the disk. We also weakly detect and spectrally resolve the CS (7--6) line but do not model it. The abundances we derive for CO and HCO$^+$ are generally consistent with expected values from chemical modeling of protoplanetary disks, while the HCN abundance is higher than expected. We dynamically measure the mass of the central star to be $2.17\pm0.07\,M_\odot$ which is inconsistent with the previously determined spectral type of K5. We also report the detection of a spatially unresolved high-velocity blue-shifted excess emission feature with a measurable positional offset from the central star, consistent with a Keplerian orbit at $60\pm20\,\mathrm{au}$. Using the integrated flux of the feature in HCO$^+$ (4--3), we estimate the total H$_2$ gas mass of this feature to be at least $1.8-8\,M_\mathrm{Jupiter}$, depending on the assumed temperature. The feature is due to a local temperature and/or density enhancement consistent with either a hydrodynamic vortex or the expected signature of the envelope of a forming protoplanet within the disk. △ Less

Submitted 6 April, 2017; originally announced April 2017.

Comments: 19 pages, 12 figures, accepted for publication in AJ

IN-SYNC. V. Stellar kinematics and dynamics in the Orion A Molecular Cloud

Authors: Nicola Da Rio, Jonathan C. Tan, Kevin R. Covey, Michiel Cottaar, Jonathan B. Foster, Nicholas C. Cullen, John J. Tobin, Jinyoung S. Kim, Michael R. Meyer, David L. Nidever, Keivan G. Stassun, S. Drew Chojnowski, Kevin M. Flaherty, Steve Majewski, Michael F. Skrutskie, Gail Zasowski, Kaike Pan

Abstract: The kinematics and dynamics of young stellar populations enable us to test theories of star formation. With this aim, we continue our analysis of the SDSS-III/APOGEE IN-SYNC survey, a high resolution near infrared spectroscopic survey of young clusters. We focus on the Orion A star-forming region, for which IN-SYNC obtained spectra of $\sim2700$ stars. In Paper IV we used these data to study the y… ▽ More The kinematics and dynamics of young stellar populations enable us to test theories of star formation. With this aim, we continue our analysis of the SDSS-III/APOGEE IN-SYNC survey, a high resolution near infrared spectroscopic survey of young clusters. We focus on the Orion A star-forming region, for which IN-SYNC obtained spectra of $\sim2700$ stars. In Paper IV we used these data to study the young stellar population. Here we study the kinematic properties through radial velocities ($v_r$). The young stellar population remains kinematically associated with the molecular gas, following a $\sim10\:{\rm{km\:s}}^{-1}$ gradient along filament. However, near the center of the region, the $v_r$ distribution is slightly blueshifted and asymmetric; we suggest that this population, which is older, is slightly in foreground. We find evidence for kinematic subclustering, detecting statistically significant groupings of co-located stars with coherent motions. These are mostly in the lower-density regions of the cloud, while the ONC radial velocities are smoothly distributed, consistent with it being an older, more dynamically evolved cluster. The velocity dispersion $σ_v$ varies along the filament. The ONC appears virialized, or just slightly supervirial, consistent with an old dynamical age. Here there is also some evidence for on-going expansion, from a $v_r$--extinction correlation. In the southern filament, $σ_v$ is $\sim2$--$3$ times larger than virial in the L1641N region, where we infer a superposition along the line of sight of stellar sub-populations, detached from the gas. On the contrary, $σ_v$ decreases towards L1641S, where the population is again in agreement with a virial state. △ Less

Submitted 2 August, 2017; v1 submitted 14 February, 2017; originally announced February 2017.

Comments: 14 pages, 13 figures, ApJ accepted

A "Rosetta Stone" for protoplanetary disks: The synergy of multi-wavelength observations

Authors: A. Sicilia-Aguilar, A. Banzatti, A. Carmona, T. Stolker, M. Kama, I Mendigutía, A. Garufi, K. Flaherty, N. van der Marel, J. Greaves

Abstract: The recent progress in instrumentation and telescope development has brought us different ways to observe protoplanetary disks, including interferometers, space missions, adaptive optics, polarimetry, and time- and spectrally-resolved data. While the new facilities have changed the way we can tackle the existing open problems in disk structure and evolution, there is a substantial lack of intercon… ▽ More The recent progress in instrumentation and telescope development has brought us different ways to observe protoplanetary disks, including interferometers, space missions, adaptive optics, polarimetry, and time- and spectrally-resolved data. While the new facilities have changed the way we can tackle the existing open problems in disk structure and evolution, there is a substantial lack of interconnection between different observing techniques and their user communities. Here, we explore the complementarity of some of the state-of-the-art observing techniques, and how they can be brought together in a collective effort to understand how disks evolve and disperse at the time of planet formation. This paper was born at the "Protoplanetary Discussions" meeting in Edinburgh, 2016. Its goal is to clarify where multi-wavelength observations of disks converge in unveiling disk structure and evolution, and where they diverge and challenge our current understanding. We discuss caveats that should be considered when linking results from different observations, or when drawing conclusions based on limited datasets (in terms of wavelength or sample). We focus on disk properties that are currently being revolutionized by multi-wavelength observations. Specifically: the inner disk radius, holes and gaps and their link to large-scale disk structures, the disk mass, and the accretion rate. We discuss how the links between them, as well as the apparent contradictions, can help us to disentangle the disk physics and to learn about disk evolution. △ Less

Submitted 10 November, 2016; v1 submitted 6 November, 2016; originally announced November 2016.

Comments: Accepted for publication in PASA. 37 pages, 9 figures, 1 table. Revised version: corrected problem in Fig 2

Spitzer Observations of Long Term Infrared Variability Among Young Stellar Objects in Chamaeleon I

Authors: Kevin M. Flaherty, Lindsay DeMarchi, James Muzerolle, Zoltan Balog, William Herbst, S. Thomas Megeath, Elise Furlan, Robert Gutermuth

Abstract: Infrared variability is common among young stellar objects, with surveys finding daily to weekly fluctuations of a few tenths of a magnitude. Space-based observations can produce highly sampled infrared light curves, but are often limited to total baselines of about a month due to the orientation of the spacecraft. Here we present observations of the Chameleon I cluster whose low declination makes… ▽ More Infrared variability is common among young stellar objects, with surveys finding daily to weekly fluctuations of a few tenths of a magnitude. Space-based observations can produce highly sampled infrared light curves, but are often limited to total baselines of about a month due to the orientation of the spacecraft. Here we present observations of the Chameleon I cluster whose low declination makes it observable by the Spitzer space telescope over a 200 day period. We observe 30 young stellar objects with a daily cadence to better sample variability on timescales of months. We find such variability is common, occurring in ~80% of the detected cluster members. The change in [3.6]-[4.5] color over 200 days for many of the sources falls between that expected for extinction and fluctuations in disk emission. With our high cadence and long baseline we can derive power spectral density curves covering two orders of magnitude in frequency and find significant power at low frequencies, up to the boundaries of our 200 day survey. Such long timescales are difficult to explain with variations driven by the interaction between the disk and stellar magnetic field, which has a dynamical timescale of days to weeks. The most likely explanation is either structural or temperature fluctuations spread throughout the inner ~0.5 au of the disk, suggesting that the intrinsic dust structure is highly dynamic. △ Less

Submitted 28 September, 2016; originally announced September 2016.

Comments: Accepted to ApJ

Debris Disks in the Scorpius-Centaurus OB Association Resolved by ALMA

Authors: Jesse Lieman-Sifry, A. Meredith Hughes, John M. Carpenter, Uma Gorti, Antonio Hales, Kevin M. Flaherty

Abstract: We present a CO(2-1) and 1240 um continuum survey of 23 debris disks with spectral types B9-G1, observed at an angular resolution of 0.5-1 arcsec with the Atacama Large Millimeter/Submillimeter Array (ALMA). The sample was selected for large infrared excess and age ~10 Myr, to characterize the prevalence of molecular gas emission in young debris disks. We identify three CO-rich debris disks, plus… ▽ More We present a CO(2-1) and 1240 um continuum survey of 23 debris disks with spectral types B9-G1, observed at an angular resolution of 0.5-1 arcsec with the Atacama Large Millimeter/Submillimeter Array (ALMA). The sample was selected for large infrared excess and age ~10 Myr, to characterize the prevalence of molecular gas emission in young debris disks. We identify three CO-rich debris disks, plus two additional tentative (3-sigma) CO detections. Twenty disks were detected in the continuum at the >3-sigma level. For the 12 disks in the sample that are spatially resolved by our observations, we perform an independent analysis of the interferometric continuum visibilities to constrain the basic dust disk geometry, as well as a simultaneous analysis of the visibilities and broad-band spectral energy distribution to constrain the characteristic grain size and disk mass. The gas-rich debris disks exhibit preferentially larger outer radii in their dust disks, and a higher prevalence of characteristic grain sizes smaller than the blowout size. The gas-rich disks do not exhibit preferentially larger dust masses, contrary to expectations for a scenario in which a higher cometary destruction rate would be expected to result in a larger mass of both CO and dust. The three debris disks in our sample with strong CO detections are all around A stars: the conditions in disks around intermediate-mass stars appear to be the most conducive to the survival or formation of CO. △ Less

Submitted 22 June, 2016; originally announced June 2016.

Comments: 16 pages, 6 figures, accepted for publication in ApJ

ALMA Observations of HD141569's Circumstellar Disk

Authors: J. A. White, A. C. Boley, A. M. Hughes, K. M. Flaherty, E. Ford, D. Wilner, S. Corder, M. Payne

Abstract: We present ALMA band 7 (345 GHz) continuum and $^{12}$CO(J = 3-2) observations of the circumstellar disk surrounding HD141569. At an age of about 5 Myr, the disk has a complex morphology that may be best interpreted as a nascent debris system with gas. Our $870\rm~μm$ ALMA continuum observations resolve a dust disk out to approximately $ 56 ~\rm au$ from the star (assuming a distance of 116 pc) wi… ▽ More We present ALMA band 7 (345 GHz) continuum and $^{12}$CO(J = 3-2) observations of the circumstellar disk surrounding HD141569. At an age of about 5 Myr, the disk has a complex morphology that may be best interpreted as a nascent debris system with gas. Our $870\rm~μm$ ALMA continuum observations resolve a dust disk out to approximately $ 56 ~\rm au$ from the star (assuming a distance of 116 pc) with $0."38$ resolution and $0.07 ~ \rm mJy~beam^{-1}$ sensitivity. We measure a continuum flux density for this inner material of $3.8 \pm 0.4 ~ \rm mJy$ (including calibration uncertainties). The $^{12}$CO(3-2) gas is resolved kinematically and spatially from about 30 to 210 au. The integrated $^{12}$CO(3-2) line flux density is $15.7 \pm 1.6~\rm Jy~km~s^{-1}$. We estimate the mass of the millimeter debris and $^{12}$CO(3-2) gas to be $\gtrsim0.04~\rm M_{\oplus}$ and $\sim2\times 10^{-3}~\rm M_{\oplus}$, respectively. If the millimeter grains are part of a collisional cascade, then we infer that the inner disk ($<50$ au) has $\sim 160~\rm M_{\oplus}$ contained within objects less than 50 km in radius, depending on the planetesimal size distribution and density assumptions. MCMC modeling of the system reveals a disk morphology with an inclination of $53.4^{\circ}$ centered around a $\rm M=2.39~ M_{\odot}$ host star ($\rm Msin(i)=1.92~ M_{\odot}$). We discuss whether the gas in HD141569's disk may be second generation. If it is, the system can be used to study the clearing stages of planet formation. △ Less

Submitted 1 June, 2016; originally announced June 2016.

Comments: 35 pages, 7 figures, submitted to ApJ

Journal ref: The Astrophysical Journal, 829:6 (11pp), 2016 September 20

Resolved Gas Interior to the Dust Rings of the HD 141569 Disk

Authors: Kevin M. Flaherty, A. Meredith Hughes, Sean M. Andrews, Chunhua Qi, David J Wilner, Aaron C. Boley, Jacob A. White, Will Harney, Julia Zachary

Abstract: The disk around HD 141569 is one of a handful of systems whose weak infrared emission is consistent with a debris disk, but still has a significant reservoir of gas. Here we report spatially resolved mm observations of the CO(3-2) and CO(1-0) emission as seen with the SMA and CARMA. We find that the excitation temperature for CO is lower than expected from cospatial blackbody grains, similar to pr… ▽ More The disk around HD 141569 is one of a handful of systems whose weak infrared emission is consistent with a debris disk, but still has a significant reservoir of gas. Here we report spatially resolved mm observations of the CO(3-2) and CO(1-0) emission as seen with the SMA and CARMA. We find that the excitation temperature for CO is lower than expected from cospatial blackbody grains, similar to previous observations of analogous systems, and derive a gas mass that lies between that of gas-rich primordial disks and gas-poor debris disks. The data also indicate a large inner hole in the CO gas distribution and an outer radius that lies interior to the outer scattered light rings. This spatial distribution, with the dust rings just outside the gaseous disk, is consistent with the expected interactions between gas and dust in an optically thin disk. This indicates that gas can have a significant effect on the location of the dust within debris disks. △ Less

Submitted 11 January, 2016; originally announced January 2016.

Comments: Accepted to ApJ

IN-SYNC IV - The Young Stellar Population in the Orion A Molecular Cloud

Authors: Nicola Da Rio, Jonathan C. Tan, Kevin R. Covey, Michiel Cottaar, Jonathan B. Foster, Nicholas C. Cullen, John J. Tobin, Jinyoung S. Kim, Michael R. Meyer, David L. Nidever, Keivan G. Stassun, S. Drew Chojnowski, Kevin M. Flaherty, Steve Majewski, Michael F. Skrutskie, Gail Zasowski, Kaike Pan

Abstract: We present the results of the SDSS APOGEE INfrared Spectroscopy of Young Nebulous Clusters program (IN-SYNC) survey of the Orion A molecular cloud. This survey obtained high resolution near infrared (NIR) spectroscopy of about 2700 young pre-main sequence stars throughout the region, acquired across five distinct fields spanning 6deg field of view (FOV). With these spectra, we have measured accura… ▽ More We present the results of the SDSS APOGEE INfrared Spectroscopy of Young Nebulous Clusters program (IN-SYNC) survey of the Orion A molecular cloud. This survey obtained high resolution near infrared (NIR) spectroscopy of about 2700 young pre-main sequence stars throughout the region, acquired across five distinct fields spanning 6deg field of view (FOV). With these spectra, we have measured accurate stellar parameters (T_eff, log g, v sin i) and extinctions, and placed the sources in the Hertzsprung-Russel Diagram (HRD). We have also extracted radial velocities for the kinematic characterization of the population. We compare our measurements with literature results for a sub-sample of targets in order to assess the performances and accuracy of the survey. Source extinction shows evidence for dust grains that are larger than those in the diffuse interstellar medium (ISM): we estimate an average R_V=5.5 in the region. Importantly, we find a clear correlation between HRD inferred ages and spectroscopic surface-gravity inferred ages. This clearly indicates a real spread of stellar radii at fixed temperature, and together with additional correlations with extinction and with disk presence, strongly suggests a real spread of ages large than a few Myr. Focussing on the young population around NGC1980 iota Ori, which has previously been suggested to be a separate, foreground, older cluster, we confirm its older (5Myr) age and low A_V, but considering that its radial velocity distribution is indistinguishable from the Orion A's population, we suggest that NGC1980 is part of Orion A's star formation activity. Based on their stellar parameters and kinematic properties, we identify 383 new candidate members of Orion A, most of which are diskless sources in areas of the region poorly studied by previous works. △ Less

Submitted 10 December, 2015; v1 submitted 12 November, 2015; originally announced November 2015.

Comments: ApJ accepted

The Spitzer Space Telescope Survey of the Orion A and B Molecular Clouds II: the Spatial Distribution and Demographics of Dusty Young Stellar Objects

Authors: S. T. Megeath, R. Gutermuth, J. Muzerolle, E. Kryukova, J. L. Hora, L. E. Allen, K. Flaherty, L. Hartmann, P. C. Myers, J. L. Pipher, J. Stauffer, E. T. Young, G. G. Fazio

Abstract: We analyze the spatial distribution of dusty young stellar objects (YSOs) identified in the Spitzer Survey of the Orion Molecular clouds, augmenting these data with Chandra X-ray observations to correct for incompleteness in dense clustered regions. We also devise a scheme to correct for spatially varying incompleteness when X-ray data are not available. The local surface densities of the YSOs ran… ▽ More We analyze the spatial distribution of dusty young stellar objects (YSOs) identified in the Spitzer Survey of the Orion Molecular clouds, augmenting these data with Chandra X-ray observations to correct for incompleteness in dense clustered regions. We also devise a scheme to correct for spatially varying incompleteness when X-ray data are not available. The local surface densities of the YSOs range from 1 pc$^{-2}$ to over 10,000 pc$^{-2}$, with protostars tending to be in higher density regions. This range of densities is similar to other surveyed molecular clouds with clusters, but broader than clouds without clusters. By identifying clusters and groups as continuous regions with surface densities $\ge10$ pc$^{-2}$, we find that 59% of the YSOs are in the largest cluster, the Orion Nebular Cluster (ONC), while 13% of the YSOs are found in a distributed population. A lower fraction of protostars in the distributed population is evidence that it is somewhat older than the groups and clusters. An examination of the structural properties of the clusters and groups show that the peak surface densities of the clusters increase approximately linearly with the number of members. Furthermore, all clusters with more than 70 members exhibit asymmetric and/or highly elongated structures. The ONC becomes azimuthally symmetric in the inner 0.1 pc, suggesting that the cluster is only $\sim 2$ Myr in age. We find the star formation efficiency (SFE) of the Orion B cloud is unusually low, and that the SFEs of individual groups and clusters are an order of magnitude higher than those of the clouds. Finally, we discuss the relationship between the young low mass stars in the Orion clouds and the Orion OB 1 association, and we determine upper limits to the fraction of disks that may be affected by UV radiation from OB stars or by dynamical interactions in dense, clustered regions. △ Less

Submitted 3 November, 2015; originally announced November 2015.

Comments: Accepted to AJ. Higher resolution and two column version of paper as well as photometry data available at http://astro1.physics.utoledo.edu/~megeath/Orion/The_Spitzer_Orion_Survey.html

Weak Turbulence in the HD 163296 Protoplanetary Disk Revealed by ALMA CO Observations

Authors: Kevin M. Flaherty, A. Meredith Hughes, Katherine A. Rosenfeld, Sean M. Andrews, Eugene Chiang, Jacob B. Simon, Skylar Kerzner, David J. Wilner

Abstract: Turbulence can transport angular momentum in protoplanetary disks and influence the growth and evolution of planets. With spatially and spectrally resolved molecular emission line measurements provided by (sub)millimeter interferometric observations, it is possible to directly measure non-thermal motions in the disk gas that can be attributed to this turbulence. We report a new constraint on the t… ▽ More Turbulence can transport angular momentum in protoplanetary disks and influence the growth and evolution of planets. With spatially and spectrally resolved molecular emission line measurements provided by (sub)millimeter interferometric observations, it is possible to directly measure non-thermal motions in the disk gas that can be attributed to this turbulence. We report a new constraint on the turbulence in the disk around HD 163296, a nearby young A star, determined from ALMA Science Verification observations of four CO emission lines (the CO(3-2), CO(2-1), 13CO(2-1), and C18O(2-1) transitions). The different optical depths for these lines permit probes of non-thermal line-widths at a range of physical conditions (temperature and density) and depths into the disk interior. We derive stringent limits on the non-thermal motions in the upper layers of the outer disk such that any contribution to the line-widths from turbulence is <3% of the local sound speed. These limits are approximately an order of magnitude lower than theoretical predictions for full-blown MHD turbulence driven by the magneto-rotational instability, potentially suggesting that this mechanism is less efficient in the outer (R>30AU) disk than has been previously considered. △ Less

Submitted 5 October, 2015; originally announced October 2015.

Comments: Accepted to ApJ

IN-SYNC III: The dynamical state of IC 348 - A super-virial velocity dispersion and a puzzling sign of convergence

Authors: Michiel Cottaar, Kevin R. Covey, Jonathan B. Foster, Michael R. Meyer, Jonathan C. Tan, David L. Nidever, S. Drew Chojnowski, Nicola da Rio, Kevin M. Flaherty, Peter M. Frinchaboy, Steve Majewski, Michael F. Skrutskie, John C. Wilson, Gail Zasowski

Abstract: Most field stars will have encountered the highest stellar density and hence the largest number of interactions in their birth environment. Yet the stellar dynamics during this crucial phase are poorly understood. Here we analyze the radial velocities measured for 152 out of 380 observed stars in the 2-6 Myr old star cluster IC 348 as part of the SDSS-III APOGEE. The radial velocity distribution o… ▽ More Most field stars will have encountered the highest stellar density and hence the largest number of interactions in their birth environment. Yet the stellar dynamics during this crucial phase are poorly understood. Here we analyze the radial velocities measured for 152 out of 380 observed stars in the 2-6 Myr old star cluster IC 348 as part of the SDSS-III APOGEE. The radial velocity distribution of these stars is fitted with one or two Gaussians, convolved with the measurement uncertainties including binary orbital motions. Including a second Gaussian improves the fit; the high-velocity outliers that are best fit by this second component may either (1) be contaminants from the nearby Perseus OB2 association, (2) be a halo of ejected or dispersing stars from IC 348, or (3) reflect that IC 348 has not relaxed to a Gaussian velocity distribution. We measure a velocity dispersion for IC 348 of $0.72 \pm 0.07$ km s$^{-1}$ (or $0.64 \pm 0.08$ km s$^{-1}$ if two Gaussians are fitted), which implies a supervirial state, unless the gas contributes more to the gravitational potential than expected. No evidence is found for a dependence of this velocity dispersion on distance from the cluster center or stellar mass. We also find that stars with lower extinction (in the front of the cloud) tend to be redshifted compared with stars with somewhat higher extinction (towards the back of the cloud). This data suggests that the stars in IC 348 are converging along the line of sight. We show that this correlation between radial velocity and extinction is unlikely to be spuriously caused by the small cluster rotation of $0.024 \pm 0.013$ km s$^{-1}$ arcmin$^{-1}$ or by correlations between the radial velocities of neighboring stars. This signature, if confirmed, will be the first detection of line-of-sight convergence in a star cluster(...) △ Less

Submitted 27 May, 2015; originally announced May 2015.

Comments: online data available at http://www.astro.ufl.edu/insync/perseus.html

Signatures of MRI-Driven Turbulence in Protoplanetary Disks: Predictions for ALMA Observations

Authors: Jacob B. Simon, A. Meredith Hughes, Kevin M. Flaherty, Xue-Ning Bai, Philip J. Armitage

Abstract: Spatially resolved observations of molecular line emission have the potential to yield unique constraints on the nature of turbulence within protoplanetary disks. Using a combination of local non-ideal magnetohydrodynamic simulations and radiative transfer calculations, tailored to properties of the disk around HD 163296, we assess the ability of ALMA to detect turbulence driven by the magnetorota… ▽ More Spatially resolved observations of molecular line emission have the potential to yield unique constraints on the nature of turbulence within protoplanetary disks. Using a combination of local non-ideal magnetohydrodynamic simulations and radiative transfer calculations, tailored to properties of the disk around HD 163296, we assess the ability of ALMA to detect turbulence driven by the magnetorotational instability (MRI). Our local simulations show that the MRI produces small-scale turbulent velocity fluctuations that increase in strength with height above the mid-plane. For a set of simulations at different disk radii, we fit a Maxell-Boltzmann distribution to the turbulent velocity and construct a turbulent broadening parameter as a function of radius and height. We input this broadening into radiative transfer calculations to quantify observational signatures of MRI-driven disk turbulence. We find that the ratio of the peak line flux to the flux at line center is a robust diagnostic of turbulence that is only mildly degenerate with systematic uncertainties in disk temperature. For the CO(3-2) line, which we expect to probe the most magnetically active slice of the disk column, variations in the predicted peak-to-trough ratio between our most and least turbulent models span a range of approximately 15%. Additional independent constraints can be derived from the morphology of spatially resolved line profiles, and we estimate the resolution required to detect turbulence on different spatial scales. We discuss the role of lower optical depth molecular tracers, which trace regions closer to the disk mid-plane where velocities in MRI-driven models are systematically lower. △ Less

Submitted 22 June, 2015; v1 submitted 12 January, 2015; originally announced January 2015.

Comments: accepted to ApJ after minor revisions

The Eleventh and Twelfth Data Releases of the Sloan Digital Sky Survey: Final Data from SDSS-III

Authors: Shadab Alam, Franco D. Albareti, Carlos Allende Prieto, F. Anders, Scott F. Anderson, Brett H. Andrews, Eric Armengaud, Éric Aubourg, Stephen Bailey, Julian E. Bautista, Rachael L. Beaton, Timothy C. Beers, Chad F. Bender, Andreas A. Berlind, Florian Beutler, Vaishali Bhardwaj, Jonathan C. Bird, Dmitry Bizyaev, Cullen H. Blake, Michael R. Blanton, Michael Blomqvist, John J. Bochanski, Adam S. Bolton, Jo Bovy, A. Shelden Bradley , et al. (249 additional authors not shown)

Abstract: The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11… ▽ More The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra. △ Less

Submitted 21 May, 2015; v1 submitted 5 January, 2015; originally announced January 2015.

Comments: DR12 data are available at http://www.sdss3.org/dr12. 30 pages. 11 figures. Accepted to ApJS

IN-SYNC II: Virial Stars from Sub-Virial Cores -- The Velocity Dispersion of Embedded Pre-Main-Sequence Stars in NGC 1333

Authors: Jonathan B. Foster, Michiel Cottaar, Kevin R. Covey, Héctor G. Arce, Michael R. Meyer, David L. Nidever, Keivan G. Stassun, Jonathan C. Tan, S. Drew Chojnowski, Nicola da Rio, Kevin M. Flaherty, Luisa Rebull, Peter M. Frinchaboy, Steven R. Majewski, Michael Skrutskie, John C. Wilson, Gail Zasowski

Abstract: The initial velocity dispersion of newborn stars is a major unconstrained aspect of star formation theory. Using near-infrared spectra obtained with the APOGEE spectrograph, we show that the velocity dispersion of young (1-2 Myr) stars in NGC 1333 is 0.92+/-0.12 km/s after correcting for measurement uncertainties and the effect of binaries. This velocity dispersion is consistent with the virial ve… ▽ More The initial velocity dispersion of newborn stars is a major unconstrained aspect of star formation theory. Using near-infrared spectra obtained with the APOGEE spectrograph, we show that the velocity dispersion of young (1-2 Myr) stars in NGC 1333 is 0.92+/-0.12 km/s after correcting for measurement uncertainties and the effect of binaries. This velocity dispersion is consistent with the virial velocity of the region and the diffuse gas velocity dispersion, but significantly larger than the velocity dispersion of the dense, star-forming cores, which have a sub-virial velocity dispersion of 0.5 km/s. Since the NGC 1333 cluster is dynamically young and deeply embedded, this measurement provides a strong constraint on the initial velocity dispersion of newly-formed stars. We propose that the difference in velocity dispersion between stars and dense cores may be due to the influence of a 70 micro-Gauss magnetic field acting on the dense cores, or be the signature of a cluster with initial sub-structure undergoing global collapse. △ Less

Submitted 21 November, 2014; originally announced November 2014.

Comments: Accepted to ApJ. 23 pages, 9 figures and 3 tables

Time-monitoring Observations of Br$γ$ Emission from Young Stars

Authors: J. A. Eisner, G. H. Rieke, M. J. Rieke, K. M. Flaherty, Jordan M. Stone, T. J. Arnold, S. R. Cortes, E. Cox, C. Hawkins, A. Cole, S. Zajac, A. L. Rudolph

Abstract: We present multiple epochs of near-IR spectroscopy for a sample of 25 young stars, including T Tauri, Herbig Ae/Be, and FU Ori objects. Using the FSPEC instrument on the Bok 90-inch telescope, we obtained K-band spectra of the BrGamma transition of hydrogen, with a resolution of ~3500. Epochs were taken over a span of >1 year, sampling time-spacings of roughly one day, one month, and one year. The… ▽ More We present multiple epochs of near-IR spectroscopy for a sample of 25 young stars, including T Tauri, Herbig Ae/Be, and FU Ori objects. Using the FSPEC instrument on the Bok 90-inch telescope, we obtained K-band spectra of the BrGamma transition of hydrogen, with a resolution of ~3500. Epochs were taken over a span of >1 year, sampling time-spacings of roughly one day, one month, and one year. The majority of our targets show BrGamma emission, and in some cases these are the first published detections. Time-variability is seen in approximately half of the targets showing BrGamma emission. We compare the observed variability with expectations for rotationally-modulated accretion onto the central stars and time-variable continuum emission or extinction from matter in the inner disk. Our observations are not entirely consistent with models of rotationally-modulated magnetospheric accretion. Further monitoring, over a larger number of epochs, will facilitate more quantitative constraints on variability timescales and amplitudes, and a more conclusive comparison with theoretical models. △ Less

Submitted 19 November, 2014; originally announced November 2014.

Comments: 16 pages in emulateapj format. Accepted for publication in MNRAS

IN-SYNC I: Homogeneous Stellar Parameters from High Resolution APOGEE Spectra for Thousands of Pre-main Sequence Star

Authors: Michiel Cottaar, Kevin R. Covey, Michael R. Meyer, David L. Nidever, Keivan G. Stassun, Jonathan B. Foster, Jonathan C. Tan, S. Drew Chojnowski, Nicola da Rio, Kevin M. Flaherty, Peter M. Frinchaboy, Michael Skrutskie, Steven R. Majewski, John C. Wilson, Gail Zasowski

Abstract: Over two years 8,859 high-resolution H-band spectra of 3,493 young (1 - 10 Myr) stars were gathered by the multi-object spectrograph of the APOGEE project as part of the IN-SYNC ancillary program of that SDSS-III survey. Here we present the forward modeling approach used to derive effective temperatures, surface gravities, radial velocities, rotational velocities, and H-band veiling from these nea… ▽ More Over two years 8,859 high-resolution H-band spectra of 3,493 young (1 - 10 Myr) stars were gathered by the multi-object spectrograph of the APOGEE project as part of the IN-SYNC ancillary program of that SDSS-III survey. Here we present the forward modeling approach used to derive effective temperatures, surface gravities, radial velocities, rotational velocities, and H-band veiling from these near-infrared spectra. We discuss in detail the statistical and systematic uncertainties in these stellar parameters. In addition we present accurate extinctions by measuring the E(J-H) of these young stars with respect to the single-star photometric locus in the Pleiades. Finally we identify an intrinsic stellar radius spread of about 25% for late-type stars in IC 348 using three (nearly) independent measures of stellar radius, namely the extinction-corrected J-band magnitude, the surface gravity and the $R \sin i$ from the rotational velocities and literature rotation periods. We exclude that this spread is caused by uncertainties in the stellar parameters by showing that the three estimators of stellar radius are correlated, so that brighter stars tend to have lower surface gravities and larger $R \sin i$ than fainter stars at the same effective temperature. Tables providing the spectral and photometric parameters for the Pleiades and IC 348 have been provided online. △ Less

Submitted 29 August, 2014; originally announced August 2014.

Comments: online tables can be found at http://www.astro.ufl.edu/insync/

Young Stellar Object Variability (YSOVAR): Long Timescale Variations in the Mid-Infrared

Authors: L. M. Rebull, A. M. Cody, K. R. Covey, H. M. Guenther, L. A. Hillenbrand, P. Plavchan, K. Poppenhaeger, J. R. Stauffer, S. J. Wolk, R. Gutermuth, M. Morales-Calderon, I. Song, D. Barrado, A. Bayo, D. James, J. L. Hora, F. J. Vrba, C. Alves de Oliveira, J. Bouvier, S. J. Carey, J. M. Carpenter, F. Favata, K. Flaherty, J. Forbrich, J. Hernandez , et al. (10 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 28 August, 2014; originally announced August 2014.

Comments: Accepted to AJ; 38 figures, 93 pages

Connecting X-ray and Infrared Variability among Young Stellar Objects: Ruling out potential sources of disk fluctuations

Authors: Kevin Flaherty, James Muzerolle, Scott Wolk, George Rieke, Robert Gutermuth, Zoltan Balog, William Herbst, S. Thomas Megeath, Elise Furlan

Abstract: Variability in the infrared emission from disks around pre-main sequence stars over the course of days to weeks appears to be common, but the physical cause of the changes in disk structure are not constrained. Here we present coordinated monitoring of one young cluster with the Spitzer and Chandra space telescopes aimed at studying the physical source of the variability. In fall 2011 we obtained… ▽ More Variability in the infrared emission from disks around pre-main sequence stars over the course of days to weeks appears to be common, but the physical cause of the changes in disk structure are not constrained. Here we present coordinated monitoring of one young cluster with the Spitzer and Chandra space telescopes aimed at studying the physical source of the variability. In fall 2011 we obtained ten epochs of Chandra ACIS photometry over a period of 30 days with a roughly 3 day cadence contemporaneous with 20 epochs of Spitzer [3.6],[4.5] photometry over 40 days with a roughly 2 day cadence of the IC 348 cluster. This cadence allows us to search for week to month long responses of the infrared emission to changes in the high-energy flux. We find no strong evidence for a direct link between the X-ray and infrared variability on these timescales among 39 cluster members with circumstellar disks. There is no significant correlation between the shape of the infrared and X-ray light curves, or between the size of the X-ray and infrared variability. Among the stars with an X-ray flare none showed evidence of a correlated change in the infrared photometry on timescales of days to weeks following the flare. This lack of connection implies that X-ray heating of the planet forming region of the disk is not significant, although we cannot rule out rapid or instantaneous changes in infrared emission. △ Less

Submitted 18 July, 2014; originally announced July 2014.

Comments: 7 pages of text, 5 figures (+1 figure set)

The extraordinary far-infrared variation of a protostar: Herschel/PACS observations of LRLL54361

Authors: Zoltan Balog, James Muzerolle, Kevin Flaherty, Ors H. Detre, Jeroen Bouwmann, Elise Furlan, Rob Gutermuth, Attila Juhasz, John Bally, Markus Nielbock, Ulrich Klaas, Oliver Krause, Thomas Henning

Abstract: We report Herschel/PACS photometric observations at 70 μm and 160 μm of LRLL54361 - a suspected binary protostar that exhibits periodic (P=25.34 days) flux variations at shorter wavelengths (3.6 μm and 4.5 μm) thought to be due to pulsed accretion caused by binary motion. The PACS observations show unprecedented flux variation at these far-infrared wavelengths that are well cor- related with the v… ▽ More We report Herschel/PACS photometric observations at 70 μm and 160 μm of LRLL54361 - a suspected binary protostar that exhibits periodic (P=25.34 days) flux variations at shorter wavelengths (3.6 μm and 4.5 μm) thought to be due to pulsed accretion caused by binary motion. The PACS observations show unprecedented flux variation at these far-infrared wavelengths that are well cor- related with the variations at shorter wavelengths. At 70 μm the object increases its flux by a factor of six while at 160μm the change is about a factor of two, consistent with the wavelength dependence seen in the far-infrared spectra. The source is marginally resolved at 70 μm with varying FWHM. Deconvolved images of the sources show elongations exactly matching the outflow cavities traced by the scattered light observations. The spatial variations are anti-correlated with the flux variation indicating that a light echo is responsible for the changes in FWHM. The observed far-infrared flux variability indicates that the disk and en- velope of this source is periodically heated by the accretion pulses of the central source, and suggests that such long wavelength variability in general may provide a reasonable proxy for accretion variations in protostars. △ Less

Submitted 5 June, 2014; originally announced June 2014.

Comments: 20 pages, 4 figures accepted for publication in ApJ Letters

Time-monitoring Observations of the Ro-Vibrational Overtone CO bands in Young Stars

Authors: J. A. Eisner, G. H. Rieke, M. J. Rieke, K. M. Flaherty, T. J. Arnold, J. M. Stone, S. R. Cortes, E. Cox, C. Hawkins, A. Cole, S. Zajac, A. L. Rudolph

Abstract: We present near-IR spectra of a sample of T Tauri, Herbig Ae/Be, and FU Ori objects. Using the FSPEC instrument on the Bok 90-inch telescope, we obtained K-band spectra with a resolution of ~3500. Here we present spectra of the v=2->0 and v=3->1 bandheads of ro-vibrational transitions of carbon monoxide. We observed these spectra over multiple epochs spaced by a few days and approximately one mont… ▽ More We present near-IR spectra of a sample of T Tauri, Herbig Ae/Be, and FU Ori objects. Using the FSPEC instrument on the Bok 90-inch telescope, we obtained K-band spectra with a resolution of ~3500. Here we present spectra of the v=2->0 and v=3->1 bandheads of ro-vibrational transitions of carbon monoxide. We observed these spectra over multiple epochs spaced by a few days and approximately one month. Several of our targets show CO emission or absorption features. However we see little evidence of variability in these features across multiple epochs. We compare our results with previous observations, and discuss the physical implications of non-variable CO emission across the sampled timescales. △ Less

Submitted 11 June, 2013; originally announced June 2013.

Comments: 30 pages, 13 figures; Accepted for publication in MNRAS

Young Stellar Objects in Lynds 1641: Disks, Accretion, and Star Formation History

Authors: Min Fang, Jinyoung Serena Kim, Roy van Boekel, Aurora Sicilia-Aguilar, Thomas Henning, Kevin Flaherty

Abstract: We investigate the young stellar objects (YSOs) in the Lynds~1641 (L1641) cloud using multi-wavelength data including Spitzer, WISE, 2MASS, and XMM covering ~1390 YSOs across a range of evolutionary stages. In addition, we targeted a sub-sample of YSOs for optical spectroscopy. We use this data, along with archival photometric data, to derive spectral types, extinction values, masses, ages, as wel… ▽ More We investigate the young stellar objects (YSOs) in the Lynds~1641 (L1641) cloud using multi-wavelength data including Spitzer, WISE, 2MASS, and XMM covering ~1390 YSOs across a range of evolutionary stages. In addition, we targeted a sub-sample of YSOs for optical spectroscopy. We use this data, along with archival photometric data, to derive spectral types, extinction values, masses, ages, as well as accretion rates. We obtain a disk fraction of ~50% in L1641. The disk frequency is almost constant as a function of stellar mass with a slight peak at log(M_*/M_sun)\approx-0.25. The analysis of multi-epoch spectroscopic data indicates that the accretion variability of YSOs cannot explain the two orders of magnitude of scatter for YSOs with similar masses. Forty-six new transition disk (TD) objects are confirmed in this work, and we find that the fraction of accreting TDs is lower than for optically thick disks (40--45% vs. 77--79% respectively). We confirm our previous result that the accreting TDs have a similar median accretion rate to normal optically thick disks. We confirm that two star formation modes (isolated vs. clustered) exist in L1641. We find that the diskless YSOs are statistically older than the YSOs with optically-thick disks and the transition disk objects have a median age which is intermediate between the two populations. We tentatively study the star formation history in L1641 based on the age distribution and find that star formation started to be active 2--3 Myr ago. △ Less

Submitted 29 April, 2013; originally announced April 2013.

Comments: 112 pages, 46 figures, 6 tables, accepted by ApJs

Pulsed Accretion in a Variable Protostar

Authors: James Muzerolle, Elise Furlan, Kevin Flaherty, Zoltan Balog, Robert Gutermuth

Abstract: Periodic increases in luminosity arising from variable accretion rates have been predicted for some close pre-main sequence binary stars as they grow from circumbinary disks. The phenomenon is known as "pulsed accretion" and can affect the orbital evolution and mass distribution of young binaries, as well as the potential for planet formation in the circumbinary environment. Accretion variability… ▽ More Periodic increases in luminosity arising from variable accretion rates have been predicted for some close pre-main sequence binary stars as they grow from circumbinary disks. The phenomenon is known as "pulsed accretion" and can affect the orbital evolution and mass distribution of young binaries, as well as the potential for planet formation in the circumbinary environment. Accretion variability is a common feature of young stars, with a large range of amplitudes and timescales as measured from multi-epoch observations at optical and infrared wavelengths. Periodic variations consistent with pulsed accretion have been seen in only a few young binaries via optical accretion tracers, albeit intermittently with accretion luminosity variations ranging from 0-50 percent from orbit to orbit. Here we report on a young protostar (age ~10^5 yr) that exhibits periodic variability in which the infrared luminosity increases by a factor of 10 in roughly one week every 25.34 days. We attribute this to pulsed accretion associated with an unseen binary companion. The strength and regularity of this accretion signal is surprising; it may be related to the very young age of the system, which is a factor of ten younger than the other pulsed accretors previously studied. △ Less

Submitted 24 January, 2013; originally announced January 2013.

Comments: 25 pages, 9 figures

Journal ref: Nature, Volume 493, Issue 7432, pp. 378-380 (2013)

Kinks and Dents in Protoplanetary Disks: Rapid Infrared Variability as Evidence for Large Structural Perturbations

Authors: Kevin Flaherty, james Muzerolle, George Rieke, Robert Gutermuth, Zoltan Balog, William Herbst, S. T. Megeath

Abstract: We report on synoptic observations at 3.6 and 4.5micron of young stellar objects in IC 348 with 38 epochs covering 40 days. We find that among the detected cluster members, 338 at [3.6] and 269 at both [3.6] and [4.5], many are variable on daily to weekly timescales with typical fluctuations of ~0.1 mag. The fraction of variables ranges from 20% for the diskless pre-main sequence stars to 60% for… ▽ More We report on synoptic observations at 3.6 and 4.5micron of young stellar objects in IC 348 with 38 epochs covering 40 days. We find that among the detected cluster members, 338 at [3.6] and 269 at both [3.6] and [4.5], many are variable on daily to weekly timescales with typical fluctuations of ~0.1 mag. The fraction of variables ranges from 20% for the diskless pre-main sequence stars to 60% for the stars still surrounded by infalling envelopes. We also find that stars in the exposed cluster core are less variable than the stars in the dense, slightly younger, south-western ridge. This trend persists even after accounting for the underlying correlation with infrared SED type, suggesting that the change in variable fraction is not simply a reflection of the change in relative fraction of class I vs. class II sources across the cloud, but instead reflects a change in variability with age. We also see a strong correlation between infrared variability and X-ray luminosity among the class II sources. The observed variability most likely reflects large changes in the structure of the inner wall located at the dust sublimation radius. We explore the possibility that these structural perturbations could be caused by a hot spot on the star heating dust above the sublimation temperature, causing it to evaporate rapidly, and increasing the inner radius for a portion of the disk. Under a number of simplifying assumptions we show that this model can reproduce the size and timescale of the 3.6 and 4.5micron fluctuations. Regardless of its source, the infrared variability indicates that the inner disk is not a slowly evolving entity, but instead is a bubbling, warped, dented mass of gas and dust whose global size and shape fluctuate in a matter of days. △ Less

Submitted 7 December, 2012; originally announced December 2012.

Comments: Accepted to AJ. 25 pages, 24 figures (+75 additional online figures and 3 online tables)

The Spitzer Space Telescope Survey of the Orion A & B Molecular Clouds - Part I: A Census of Dusty Young Stellar Objects and a Study of their Mid-IR Variability

Authors: S. T. Megeath, R. Gutermuth, J. Muzerolle, E. Kryukova, K. Flaherty, J. Hora, L E. Allen, L. Hartmann, P. C. Myers, J. L. Pipher, J. Stauffer, E. T. Young, G. G. Fazio

Abstract: We present a survey of the Orion A and B molecular clouds undertaken with the IRAC and MIPS instruments onboard Spitzer. In total, five distinct fields were mapped covering 9 sq. degrees in five mid-IR bands spanning 3-24 microns. The survey includes the Orion Nebula Cluster, the Lynds 1641, 1630 and 1622 dark clouds, and the NGC 2023, 2024, 2068 and 2071 nebulae. These data are merged with the 2M… ▽ More We present a survey of the Orion A and B molecular clouds undertaken with the IRAC and MIPS instruments onboard Spitzer. In total, five distinct fields were mapped covering 9 sq. degrees in five mid-IR bands spanning 3-24 microns. The survey includes the Orion Nebula Cluster, the Lynds 1641, 1630 and 1622 dark clouds, and the NGC 2023, 2024, 2068 and 2071 nebulae. These data are merged with the 2MASS point source catalog to generate a catalog of eight band photometry. We identify 3479 dusty young stellar objects (YSOs) in the Orion molecular clouds by searching for point sources with mid-IR colors indicative of reprocessed light from dusty disks or infalling envelopes. The YSOs are subsequently classified on the basis of their mid-IR colors and their spatial distributions are presented. We classify 2991 of the YSOs as pre-main sequence stars with disks and 488 as likely protostars. Most of the sources were observed with IRAC in 2-3 epochs over 6 months; we search for variability between the epochs by looking for correlated variability in the 3.6 and 4.5 micron bands. We find that 50% of the dusty YSOs show variability. The variations are typically small (0.2 mag.) with the protostars showing a higher incidence of variability and larger variations. The observed correlations between the 3.6, 4.5, 5.8 and 8 micron variability suggests that we are observing variations in the heating of the inner disk due to changes in the accretion luminosity or rotating accretion hot spots. △ Less

Submitted 10 October, 2012; v1 submitted 17 September, 2012; originally announced September 2012.

Comments: Accepted to the Astronomical Journal. A full resolution version of the paper and the point source catalog can be found at http://astro1.physics.utoledo.edu/~megeath/Orion/The_Spitzer_Orion_Survey.html

Infrared Variability of Evolved Protoplanetary Disks: Evidence for Scale Height Variations in the Inner Disk

Authors: Kevin Flaherty, James Muzerolle, George Rieke, Robert Gutermuth, Zoltan Balog, William Herbst, S. T. Megeath, Maria Kun

Abstract: We present the results of a multi-wavelength multi-epoch survey of five evolved protoplanetary disks in the IC 348 cluster that show significant infrared variability. Using 3-8micron and 24micron photometry along with 5-40micron spectroscopy from the Spitzer Space Telescope, as well as ground-based 0.8-5micron spectroscopy, optical spectroscopy and near-infrared photometry, covering timescales of… ▽ More We present the results of a multi-wavelength multi-epoch survey of five evolved protoplanetary disks in the IC 348 cluster that show significant infrared variability. Using 3-8micron and 24micron photometry along with 5-40micron spectroscopy from the Spitzer Space Telescope, as well as ground-based 0.8-5micron spectroscopy, optical spectroscopy and near-infrared photometry, covering timescales of days to years, we examine the variability in the disk, stellar and accretion flux. We find substantial variations (10-60%) at all infrared wavelengths on timescales of weeks to months for all of these young stellar objects. This behavior is not unique when compared to other cluster members and is consistent with changes in the structure of the inner disk, most likely scale height fluctuations on a dynamical timescale. Previous observations, along with our near-infrared photometry, indicate that the stellar fluxes are relatively constant; stellar variability does not appear to drive the large changes in the infrared fluxes. Based on our near-infrared spectroscopy of the Pa-beta and Br-gamma lines we find that the accretion rates are variable in most of the evolved disks but the overall rates are probably too small to cause the infrared variability. We discuss other possible physical causes for the variability, including the influence of a companion, magnetic fields threading the disk, and X-ray flares. △ Less

Submitted 7 February, 2012; originally announced February 2012.

Comments: Accepted to ApJ. 33 pages, emulate apj format

The Highly Dynamic Behavior of the Innermost Dust and Gas in the Transition Disk Variable LRLL 31

Authors: Kevin Flaherty, James Muzerolle, George Rieke, Robert Gutermuth, Zoltan Balog, William Herbst, S. T. Megeath, Maria Kun

Abstract: We describe extensive synoptic multi-wavelength observations of the transition disk LRLL 31 in the young cluster IC 348. We combined four epochs of IRS spectra, nine epochs of MIPS photometry, seven epochs of cold-mission IRAC photometry and 36 epochs of warm mission IRAC photometry along with multi-epoch near-infrared spectra, optical spectra and polarimetry to explore the nature of the rapid var… ▽ More We describe extensive synoptic multi-wavelength observations of the transition disk LRLL 31 in the young cluster IC 348. We combined four epochs of IRS spectra, nine epochs of MIPS photometry, seven epochs of cold-mission IRAC photometry and 36 epochs of warm mission IRAC photometry along with multi-epoch near-infrared spectra, optical spectra and polarimetry to explore the nature of the rapid variability of this object. We find that the inner disk, as traced by the 2-5micron excess stays at the dust sublimation radius while the strength of the excess changes by a factor of 8 on weekly timescales, and the 3.6 and 4.5micron photometry shows a drop of 0.35 magnitudes in one week followed by a slow 0.5 magnitude increase over the next three weeks. The accretion rate, as measured by PaBeta and BrGamma emission lines, varies by a factor of five with evidence for a correlation between the accretion rate and the infrared excess. While the gas and dust in the inner disk are fluctuating the central star stays relatively static. Our observations allow us to put constraints on the physical mechanism responsible for the variability. The variabile accretion, and wind, are unlikely to be causes of the variability, but both are effects of the same physical process that disturbs the disk. The lack of periodicity in our infrared monitoring indicates that it is unlikely that there is a companion within ~0.4 AU that is perturbing the disk. The most likely explanation is either a companion beyond ~0.4 AU or a dynamic interface between the stellar magnetic field and the disk leading to a variable scale height and/or warping of the inner disk. △ Less

Submitted 3 March, 2011; originally announced March 2011.

Comments: Accepted to ApJ. 10 pages of text, plus 11 tables and 13 figures at the end

Young Stars in the Time Domain: A Cool Stars 16 Splinter Summary

Authors: Kevin R. Covey, Peter Plavchan, Fabienne Bastien, Ettore Flaccomio, Kevin Flaherty, Stephen Marsden, Maria Morales-Calderón, James Muzerolle, Neal J. Turner

Abstract: Variability is a defining characteristic of young stellar systems, and optical variability has been heavily studied to select and characterize the photospheric properties of young stars. In recent years, multi-epoch observations sampling a wider range of wavelengths and time-scales have revealed a wealth of time-variable phenomena at work during the star formation process. This splinter session wa… ▽ More Variability is a defining characteristic of young stellar systems, and optical variability has been heavily studied to select and characterize the photospheric properties of young stars. In recent years, multi-epoch observations sampling a wider range of wavelengths and time-scales have revealed a wealth of time-variable phenomena at work during the star formation process. This splinter session was convened to summarize recent progress in providing improved coverage and understanding of time-variable processes in young stars and circumstellar disks. We begin by summarizing results from several multi-epoch Spitzer campaigns, which have demonstrated that many young stellar objects evidence significant mid-IR variability. While some of these variations can be attributed to processes in the stellar photosphere, others appear to trace short time-scale changes in the circumstellar disk which can be successfully modeled with axisymmetric or non-axisymmetric structures. We also review recent studies probing variability at shorter wavelengths that provide evidence for high frequency pulsations associated with accretion outbursts, correlated optical/X-ray variability in Classical T Tauri stars, and magnetic reversals in young solar analogs. △ Less

Submitted 16 December, 2010; originally announced December 2010.

Comments: 4 figures, 13 pages in ASP format; to be published in the proceedings of the Cool Stars 16 Workshop

Modeling Mid-Infrared Variability of Circumstellar Disks with Non-Axisymmetric Structure

Authors: Kevin Flaherty, James Muzerolle

Abstract: Recent mid-infrared observations of young stellar objects have found significant variations possibly indicative of changes in the structure of the circumstellar disk. Previous models of this variability have been restricted to axisymmetric perturbations in the disk. We consider simple models of a non-axisymmetric variation in the inner disk, such as a warp or a spiral wave. We find that the preces… ▽ More Recent mid-infrared observations of young stellar objects have found significant variations possibly indicative of changes in the structure of the circumstellar disk. Previous models of this variability have been restricted to axisymmetric perturbations in the disk. We consider simple models of a non-axisymmetric variation in the inner disk, such as a warp or a spiral wave. We find that the precession of these non-axisymmetric structures produce negligible flux variations but a change in the height of these structures can lead to significant changes in the mid-infrared flux. Applying these models to observations of the young stellar object LRLL 31 suggests that the observed variability could be explained by a warped inner disk with variable scale height. This suggests that some of the variability observed in young stellar objects could be explained by non-axisymmetric disturbances in the inner disk and this variability would be easily observable in future studies. △ Less

Submitted 7 July, 2010; originally announced July 2010.

Comments: 9 pages plus 16 figures and 1 appendix. Accepted to ApJ

Evidence for Dynamical Changes in a Transitional Protoplanetary Disk with Mid-infrared Variability

Authors: James Muzerolle, Kevin Flaherty, Zoltan Balog, Elise Furlan, Paul S. Smith, Lori Allen, Nuria Calvet, Paola D'Alessio, S. Thomas Megeath, August Muench, George H. Rieke, William H. Sherry

Abstract: We present multi-epoch Spitzer Space Telescope observations of the transitional disk LRLL 31 in the 2-3 Myr-old star forming region IC 348. Our measurements show remarkable mid-infrared variability on timescales as short as one week. The infrared continuum emission exhibits systematic wavelength-dependent changes that suggest corresponding dynamical changes in the inner disk structure and variab… ▽ More We present multi-epoch Spitzer Space Telescope observations of the transitional disk LRLL 31 in the 2-3 Myr-old star forming region IC 348. Our measurements show remarkable mid-infrared variability on timescales as short as one week. The infrared continuum emission exhibits systematic wavelength-dependent changes that suggest corresponding dynamical changes in the inner disk structure and variable shadowing of outer disk material. We propose several possible sources for the structural changes, including a variable accretion rate or a stellar or planetary companion embedded in the disk. Our results indicate that variability studies in the infrared can provide important new constraints on protoplanetary disk behavior. △ Less

Submitted 29 September, 2009; originally announced September 2009.

Comments: 15 pages, 4 figures, accepted to ApJ Letters