Source author record

Lee Hartmann

Lee Hartmann appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.SR astro-ph.GA astro-ph.EP astro-ph astro-ph.IM

Catalog footprint

What is connected

66works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

The Structure and Kinematics of Three Class 0 Protostellar Jets from JWST

We present observations of jets within 2000 au of three deeply embedded protostars using 2.9-27 micron observations with JWST. These observations show the morphologies and kinematics of the collimated jets from three protostars, the low-mass Class 0 protostars B335 and HOPS 153, and the intermediate-mass protostar HOPS 370. These jets are traced by shock-ionized fine-structure line emission observed with the JWST NIRSpec and MIRI IFUs. We find that [Fe II] emission traces the full extent of the inner 1000 to 2000 au of the jets, depending on distance to the protostar, while other ions mostly trace isolated shocked knots. The jets show evidence of wiggling motion in the plane of the sky as well as asymmetries between blue and red-shifted lobes. The widths of the jets increase non-monotonically with distance from the central protostar, with opening angles ranging from 2.1 degrees to < 10.1 degrees for the three protostars in the sample. The jets have total velocities ranging from 147 to 184 km/s after correcting for disk inclination. For B335, an 8-month gap between NIRSpec and MIRI MRS observations enabled measurement of the tangential velocity of a shocked knot; in combination with the radial velocity, this shows that the jet has a different inclination than the outflow cavity. We find multiple knots before and during a recent outburst in B335, although the knots were more frequent during the burst. The asymmetries between blue- and red-shifted lobes strongly suggest complex interactions between the circumstellar disks and magnetic fields.

preprint2022arXiv

Anisotropic Infall and Substructure formation in Embedded Disks

The filamentary nature of accretion streams found around embedded sources suggest that protostellar disks experience heterogenous infall from the star-forming environment, consistent with the accretion behavior onto star-forming cores in top-down star-cluster formation simulations. This may produce disk substructures in the form of rings, gaps, and spirals continuing to be identified by high-resolution imaging surveys in both embedded Class 0/I and later Class II sources. We present a parameter study of anisotropic infall, informed by the properties of accretion flows onto protostellar cores in numerical simulations, and varying the relative specific angular momentum of incoming flows as well as their flow geometry. Our results show that anisotropic infall perturbs the disk and readily launches the Rossby wave instability (RWI). It forms vortices at the inner and outer edge of the infall zone where material is deposited. These vortices drive spiral waves and angular momentum transport, with some models able to drive stresses corresponding to a viscosity parameter on the order of $α\sim 10^{-2}$. The resulting azimuthal shear forms robust pressure bumps that act as barriers to radial drift of dust grains, as demonstrated by post-processing calculations of drift-dominated dust evolution. We discuss how a self-consistent model of anisotropic infall can account for the formation of millimeter rings in the outer disk as well as producing compact dust disks, consistent with observations of embedded sources.

preprint2022arXiv

Gemini-LIGHTS: Herbig Ae/Be and massive T-Tauri protoplanetary disks imaged with Gemini Planet Imager

We present the complete sample of protoplanetary disks from the Gemini- Large Imaging with GPI Herbig/T-tauri Survey (Gemini-LIGHTS) which observed bright Herbig Ae/Be stars and T-Tauri stars in near-infrared polarized light to search for signatures of disk evolution and ongoing planet formation. The 44 targets were chosen based on their near- and mid-infrared colors, with roughly equal numbers of transitional, pre-transitional, and full disks. Our approach explicitly did not favor well-known, "famous" disks or those observed by ALMA, resulting in a less-biased sample suitable to probe the major stages of disk evolution during planet formation. Our optimized data reduction allowed polarized flux as low as 0.002% of the stellar light to be detected, and we report polarized scattered light around 80% of our targets. We detected point-like companions for 47% of the targets, including 3 brown dwarfs (2 confirmed, 1 new), and a new super-Jupiter mass candidate around V1295 Aql. We searched for correlations between the polarized flux and system parameters, finding a few clear trends: presence of a companion drastically reduces the polarized flux levels, far-IR excess correlates with polarized flux for non-binary systems, and systems hosting disks with ring structures have stellar masses $<$ 3 Msun. Our sample also included four hot, dusty "FS CMa" systems and we detected large-scale ($>100$ au) scattered light around each, signs of extreme youth for these enigmatic systems. Science-ready images are publicly available through multiple distribution channels using a new FITS file standard jointly developed with members of the VLT/SPHERE team.

preprint2022arXiv

The ODYSSEUS Survey. Motivation and First Results: Accretion, Ejection, and Disk Irradiation of CVSO 109

The Hubble UV Legacy Library of Young Stars as Essential Standards (ULLYSES) Director's Discretionary Program of low-mass pre-main-sequence stars, coupled with forthcoming data from ALMA and JWST, will provide the foundation to revolutionize our understanding of the relationship between young stars and their protoplanetary disks. A comprehensive evaluation of the physics of disk evolution and planet formation requires understanding the intricate relationships between mass accretion, mass outflow, and disk structure. Here we describe the Outflows and Disks around Young Stars: Synergies for the Exploration of ULLYSES Spectra (ODYSSEUS) Survey and present initial results of the classical T Tauri Star CVSO 109 in Orion OB1b as a demonstration of the science that will result from the survey. ODYSSEUS will analyze the ULLYSES spectral database, ensuring a uniform and systematic approach in order to (1) measure how the accretion flow depends on the accretion rate and magnetic structures, (2) determine where winds and jets are launched and how mass-loss rates compare with accretion, and (3) establish the influence of FUV radiation on the chemistry of the warm inner regions of planet-forming disks. ODYSSEUS will also acquire and provide contemporaneous observations at X-ray, optical, NIR, and millimeter wavelengths to enhance the impact of the ULLYSES data. Our goal is to provide a consistent framework to accurately measure the level and evolution of mass accretion in protoplanetary disks, the properties and magnitudes of inner-disk mass loss, and the influence of UV radiation fields that determine ionization levels and drive disk chemistry.

preprint2022arXiv

The Rate, Amplitude and Duration of Outbursts from Class 0 Protostars in Orion

At least half of a protostar's mass is accreted in the Class 0 phase, when the central protostar is deeply embedded in a dense, infalling envelope. We present the first systematic search for outbursts from Class 0 protostars in the Orion clouds. Using photometry from Spitzer/IRAC spanning 2004 to 2017, we detect three outbursts from Class 0 protostars with $\ge 2$ mag changes at 3.6 or 4.5 $μ$m. This is comparable to the magnitude change of a known protostellar FU Ori outburst. Two are newly detected bursts from the protostars HOPS 12 and 124. The number of detections implies that Class 0 protostars burst every 438 yr, with a 95% confidence interval of 161 to 1884 yr. Combining Spitzer and WISE/NEOWISE data spanning 2004-2019, we show that the bursts persist for more than nine years with significant variability during each burst. Finally, we use $19-100$ $μ$m photometry from SOFIA, Spitzer and Herschel to measure the amplitudes of the bursts. Based on the burst interval, a duration of 15 yr, and the range of observed amplitudes, 3-100% of the mass accretion during the Class 0 phase occurs during bursts. In total, we show that bursts from Class 0 protostars are as frequent, or even more frequent, than those from more evolved protostars. This is consistent with bursts being driven by instabilities in disks triggered by rapid mass infall. Furthermore, we find that bursts may be a significant, if not dominant, mode of mass accretion during the Class 0 phase.

preprint2020arXiv

A Study of Millimeter Variability in FUor Objects

FU Orionis objects (FUors) are rapidly-accreting, pre-main sequence objects that are known to exhibit large outbursts at optical and near-infrared wavelengths, with post-eruption, small-scale photometric variability superimposed on longer-term trends. In contrast, little is known about the variability of FUors at longer wavelengths. To explore this further, we observed six FUor objects using the NOrthern Extended Millimeter Array (NOEMA) and for a subset of three objects we obtained coordinated observations with NOEMA and the Lowell Discovery Telescope (LDT). In combination with previously published NOEMA observations from 2014, our 2017 observations of V1735 Cyg provide the first detection of variability in an FUor object at 2.7 mm. In the absence of significant optical variability, we discount the possibility that the mm flux density changed as a result of irradiation from the central disk. In addition, a change in the dust mass due to infall is highly unlikely. A plausible explanation for the change in 2.7 mm flux density is variability in free-free emission due to changes in the object's jet/wind. Thus, it may be that free-free emission in some FUor objects is significant at $\sim$3 mm and must be considered when deriving disk masses in order to help constrain the mechanism responsible for triggering FUor outbursts.

preprint2020arXiv

A triple star system with a misaligned and warped circumstellar disk shaped by disk tearing

Young stars are surrounded by a circumstellar disk of gas and dust, within which planet formation can occur. Gravitational forces in multiple star systems can disrupt the disk. Theoretical models predict that if the disk is misaligned with the orbital plane of the stars, the disk should warp and break into precessing rings, a phenomenon known as disk tearing. We present observations of the triple star system GWOrionis, finding evidence for disk tearing. Our images show an eccentric ring that is misaligned with the orbital planes and the outer disk. The ring casts shadows on a strongly warped intermediate region of the disk. If planets can form within the warped disk, disk tearing could provide a mechanism for forming wide-separation planets on oblique orbits.

preprint2020arXiv

Angular Momenta, Magnetization, and Accretion of Protostellar Cores

Building on our previous hydrodynamic study of the angular momenta of cloud cores formed during gravitational collapse of star-forming molecular gas in our previous work, we now examine core properties assuming ideal magnetohydrodynamics (MHD). Using the same sink-patch implementation for the \emph{Athena} MHD code, we characterize the statistical properties of cores, including the mass accretion rates, specific angular momenta, and alignments between the magnetic field and the spin axis of the core on the $0.1 \ \mathrm{pc}$ scale. Our simulations, which reproduce the observed relation between magnetic field strength and gas density, show that magnetic fields can help collimate low density flows and help seed the locations of filamentary structures. Consistent with our previous purely hydrodynamic simulations, stars (sinks) form within the heterogeneous environments of filaments, such that accretion onto cores is highly episodic leading to short-term variability but no long-term monotonic growth of the specific angular momenta. With statistical characterization of protostellar cores properties and behaviors, we aim to provide a starting point for building more realistic and self-consistent disk formation models, helping to address whether magnetic fields can prevent the development of (large) circumstellar disks in the ideal MHD limit.

preprint2020arXiv

The interpretation of protoplanetary disc wind diagnostic lines from X-ray photoevaporation and analytical MHD models

High resolution spectra of typical wind diagnostics ([OI] 6300 Å and other forbidden emission lines) can often be decomposed into multiple components: high-velocity components with blueshifts up to several 100 km/s are usually attributed to fast jets, while narrow (NLVC) and broad (BLVC) low-velocity components are believed to trace slower disc winds. Under the assumption that the line-broadening is dominated by Keplerian rotation, several studies have found that the BLVCs should trace gas launched between 0.05 and 0.5 au and correlations between the properties of BLVCs and NLVCs have been interpreted as evidence for the emission tracing an extended MHD wind and not a photoevaporative wind. We calculated synthetic line profiles obtained from detailed photoionisation calculations of an X-ray photoevaporation model and a simple MHD wind model and analyzed the emission regions of different diagnostic lines and the resulting spectral profiles. The photoevaporation model reproduces most of the observed NLVCs but not the BLVCs or HVCs. The MHD model is able to reproduce all components but produces Keplerian double peaks at average inclinations that are rarely observed. The combination of MHD and photoevaporative winds could solve this problem. Our results suggest that the Gaussian decomposition does not allow for a clear distinction of flux from different wind regions and that the line broadening is often dominated by the velocity gradient in the outflow rather than by Keplerian rotation. We show that observed correlations between BLVC and NLVC do not necessarily imply a common origin in an extended MHD wind.

preprint2019arXiv

Irregular dust features around intermediate-mass young stars with GPI: signs of youth or misaligned disks?

We are undertaking a large survey of over thirty disks using the Gemini Planet Imager (GPI) to see whether the observed dust structures match spectral energy distribution (SED) predictions and have any correlation with stellar properties. GPI can observe near-infrared light scattered from dust in circumstellar environments using high-resolution Polarimetric Differential Imaging (PDI) with coronagraphy and adaptive optics. The data have been taken in J and H bands over two years, with inner working angles of 0.08'' and 0.11'' respectively. Ahead of the release of the complete survey results, here we present five objects with extended and irregular dust structures within 2'' of the central star. These objects are: FU Ori; MWC 789; HD 45677; Hen 3-365; and HD 139614. The observed structures are consistent with each object being a pre-main-sequence star with protoplanetary dust. The five objects' circumstellar environments could result from extreme youth and complex initial conditions, from asymmetric scattering patterns due to shadows cast by misaligned disks, or in some cases from interactions with companions. We see complex U_phi structures in most objects that could indicate multiple scattering or result from the illumination of companions. Specific key findings include the first high-contrast observation of MWC 789 revealing a newly-discovered companion candidate and arc, and two faint companion candidates around Hen 3-365. These two objects should be observed further to confirm whether the companion candidates are co-moving. Further observations and modeling are required to determine the causes of the structures.

preprint2016arXiv

CSI 2264: Characterizing Young Stars in NGC 2264 with Stochastically Varying Light Curves

We provide CoRoT and Spitzer light curves, as well as broad-band multi-wavelength photometry and high resolution, multi- and single-epoch spectroscopy for 17 classical T Tauris in NGC 2264 whose CoRoT light curves (LCs) exemplify the "stochastic" LC class as defined in Cody et al. (2014). The most probable physical mechanism to explain the optical variability in this LC class is time-dependent mass accretion onto the stellar photosphere, producing transient hot spots. As evidence in favor of this hypothesis, multi-epoch high resolution spectra for a subset of these stars shows that their veiling levels also vary in time and that this veiling variability is consistent in both amplitude and timescale with the optical LC morphology. Furthermore, the veiling variability is well-correlated with the strength of the HeI 6678A emission line, a feature predicted by models to arise in accretion shocks on or near the stellar photosphere. Stars with accretion burst LC morphology (Stauffer et al. 2014) are also attributed to variable mass accretion. Both the stochastic and accretion burst LCs can be explained by a simple model of randomly occurring flux bursts, with the stochastic LC class having a higher frequency of lower amplitude events. Based on their UV excesses, veiling, and mean Ha equivalent widths, members of the stochastic LC class have only moderate time-averaged mass accretion rates. The most common feature of their Ha profiles is for them to exhibit blue-shifted absorption features, most likely originating in a disk wind. The lack of periodic signatures in the LCs suggests that little of the variability is due to long-lived hot spots rotating into or out of our line of sight; instead, the primary driver of the observed photometric variability is likely to be instabilities in the inner disk that lead to variable mass accretion.

preprint2016arXiv

Planetary Signatures in the SAO 206462 (HD 135344B) Disk: A Spiral Arm Passing Through Vortex?

The disk surrounding SAO 206462, an 8 Myr-old Herbig Ae star, has recently been reported to exhibit spiral arms, an asymmetric dust continuum, and a dust-depleted inner cavity. By carrying out two-dimensional, two-fluid hydrodynamic calculations, we find that a planetary-mass companion located at the outer disk could be responsible for these observed structures. In this model, the planet excites primary and secondary arms interior to its orbit. It also carves a gap and generates a local pressure bump at the inner gap edge where a vortex forms through Rossby wave instability. The vortex traps radially drifting dust particles, forming a dust-depleted cavity in the inner disk. We propose that the vortex is responsible for the brightest southwestern peak seen in infrared scattered light and sub-millimeter dust continuum emission. In particular, it is possible that the scattered light is boosted as one of the spiral arms passes through the high density vortex region, although the vortex alone may be able to explain the peak. We suggest that a planetary companion with a mass of 10-15 M_J is orbiting SAO 206462 at 100-120 au. Monitoring of the brightest peak over the next few years will help reveal its origin because the spiral arms and vortex will show distinguishable displacement.

preprint2016arXiv

Self-Destructing Spiral Waves: Global Simulations of a Spiral Wave Instability in Accretion Disks

We present results from a suite of three-dimensional global hydrodynamic simulations which show that spiral density waves propagating in circumstellar disks are unstable to the growth of a parametric instability that leads to break-down of the flow into turbulence. This spiral wave instability (SWI) arises from a resonant interaction between pairs of inertial waves, or inertial-gravity waves, and the background spiral wave. The development of the instability in the linear regime involves the growth of a broad spectrum of inertial modes, with growth rates on the order of the orbital time, and results in a nonlinear saturated state in which turbulent velocity perturbations are of a similar magnitude to those induced by the spiral wave. The turbulence induces angular momentum transport, and vertical mixing, at a rate that depends locally on the amplitude of the spiral wave (we obtain a stress parameter $α\sim 5 \times 10^{-4}$ in our reference model). The instability is found to operate in a wide-range of disk models, including those with isothermal or adiabatic equations of state, and in viscous disks where the dimensionless kinematic viscosity $ν\le 10^{-5}$. This robustness suggests that the instability will have applications to a broad range of astrophysical disk-related phenomena, including those in close binary systems, planets embedded in protoplanetary disks (including Jupiter in our own Solar System) and FU Orionis outburst models. Further work is required to determine the nature of the instability, and to evaluate its observational consequences, in physically more complete disk models than we have considered in this paper.

preprint2016arXiv

Spectroscopic Binaries in the Orion Nebula Cluster and NGC 2264

We examine the spectroscopic binary population for two massive nearby regions of clustered star formation, the Orion Nebula Cluster and NGC 2264, supplementing the data presented by Tobin et al. (2009, 2015) with more recent observations and more extensive analysis. The inferred multiplicity fraction up to 10 AU based on these observations is $5.3\pm 1.2$% for NGC 2264 and $5.8\pm 1.1$% for the ONC; they are consistent with the distribution of binaries in the field in the relevant parameter range. Eight of the multiple systems in the sample have enough epochs to make an initial fit for the orbital parameters. Two of these sources are double-lined spectroscopic binaries; for them we determine the mass ratio. Our reanalysis of the distribution of stellar radial velocities towards these clusters presents a significantly better agreement between stellar and gas kinematics than was previously thought.

preprint2016arXiv

The Spiral Wave Instability Induced by a Giant Planet: I. Particle Stirring in the Inner Regions of Protoplanetary Disks

We have recently shown that spiral density waves propagating in accretion disks can undergo a parametric instability by resonantly coupling with and transferring energy into pairs of inertial waves (or inertial-gravity waves when buoyancy is important). In this paper, we perform inviscid three-dimensional global hydrodynamic simulations to examine the growth and consequence of this instability operating on the spiral waves driven by a Jupiter-mass planet in a protoplanetary disk. We find that the spiral waves are destabilized via the spiral wave instability (SWI), generating hydrodynamic turbulence and sustained radially-alternating vertical flows that appear to be associated with long wavelength inertial modes. In the interval $0.3~R_p \leq R \leq 0.7~R_p$, where $R_p$ denotes the semi-major axis of the planetary orbit (assumed to be 5~au), the estimated vertical diffusion rate associated with the turbulence is characterized by $α_{\rm diff} \sim (0.2-1.2) \times 10^{-2}$. For the disk model considered here, the diffusion rate is such that particles with sizes up to several centimeters are vertically mixed within the first pressure scale height. This suggests that the instability of spiral waves launched by a giant planet can significantly disperse solid particles and trace chemical species from the midplane. In planet formation models where the continuous local production of chondrules/pebbles occurs over Myr time scales to provide a feedstock for pebble accretion onto these bodies, this stirring of solid particles may add a time constraint: planetary embryos and large asteroids have to form before a gas giant forms in the outer disk, otherwise the SWI will significantly decrease the chondrule/pebble accretion efficiency.

preprint2015arXiv

Are protoplanetary disks born with vortices? -- Rossby wave instability driven by protostellar infall

We carry out two-fluid, two-dimensional global hydrodynamic simulations to test whether protostellar infall can trigger Rossby wave instability (RWI) in protoplanetry disks. Our results show that infall can trigger the RWI and generate vortices near the outer edge of the mass landing on the disk (i.e. centrifugal radius). We find that the RWI is triggered under a variety of conditions, although the details depend on the disk parameters and the infall pattern. The common key feature of triggering the RWI is the steep radial gradient of the azimuthal velocity induced by the local increase in density at the outer edge of the infall region. Vortices form when the instability enters the nonlinear regime. In our standard model where self-gravity is neglected, vortices merge together to a single vortex within $\sim 20$ local orbital times, and the merged vortex survives for the remaining duration of the calculation ($> 170$ local orbital times). The vortex takes part in outward angular momentum transport, with a Reynolds stress of $\lesssim10^{-2}$. Our two-fluid calculations show that vortices efficiently trap dust particles with stopping times of the order of the orbital time, locally enhancing the dust to gas ratio for particles of the appropriate size by a factor of $\sim 40$ in our standard model. When self-gravity is considered, however, vortices tend to be impeded from merging and may eventually dissipate. We conclude it may well have that protoplanetary disks have favorable conditions for vortex formation during the protostellar infall phase, which might enhance early planetary core formation.

preprint2015arXiv

CSI 2264: Characterizing Young Stars in NGC 2264 with Short-Duration, Periodic Flux Dips in their Light Curves

We identify nine young stellar objects (YSOs) in the NGC 2264 star-forming region with optical {\em CoRoT} light curves exhibiting short-duration, shallow, periodic flux dips. All of these stars have infrared (IR) excesses that are consistent with their having inner disk walls near the Keplerian co-rotation radius. The repeating photometric dips have FWHM generally less than one day, depths almost always less than 15%, and periods (3<P<11 days) consistent with dust near the Keplerian co-rotation period. The flux dips vary considerably in their depth from epoch to epoch, but usually persist for several weeks and, in two cases, were present in data collected on successive years. For several of these stars, we also measure the photospheric rotation period and find that the rotation and dip periods are the same, as predicted by standard "disk-locking" models. We attribute these flux dips to clumps of material in or near the inner disk wall, passing through our line of sight to the stellar photosphere. In some cases, these dips are also present in simultaneous {\em Spitzer} IRAC light curves at 3.6 and 4.5 microns. We characterize the properties of these dips, and compare the stars with light curves exhibiting this behavior to other classes of YSO in NGC 2264. A number of physical mechanisms could locally increase the dust scale height near the inner disk wall, and we discuss several of those mechanisms; the most plausible mechanisms are either a disk warp due to interaction with the stellar magnetic field or dust entrained in funnel-flow accretion columns arising near the inner disk wall.

preprint2015arXiv

Kinematic and Spatial Substructure in NGC 2264

We present an expanded kinematic study of the young cluster NGC 2264 based upon optical radial velocities measured using multi-fiber echelle spectroscopy at the 6.5 meter MMT and Magellan telescopes. We report radial velocities for 695 stars, of which approximately 407 stars are confirmed or very likely members. Our results more than double the number of members with radial velocities from F{\H u}r{é}sz et al., resulting in a much better defined kinematic relationship between the stellar population and the associated molecular gas. In particular, we find that there is a significant subset of stars that are systematically blueshifted with respect to the molecular ($^{13}$CO) gas. The detection of Lithium absorption and/or infrared excesses in this blue-shifted population suggests that at least some of these stars are cluster members; we suggest some speculative scenarios to explain their kinematics. Our results also more clearly define the redshifted population of stars in the northern end of the cluster; we suggest that the stellar and gas kinematics of this region are the result of a bubble driven by the wind from O7 star S Mon. Our results emphasize the complexity of the spatial and kinematic structure of NGC 2264, important for eventually building up a comprehensive picture of cluster formation.

preprint2015arXiv

Signatures of star cluster formation by cold collapse

Sub-virial gravitational collapse is one mechanism by which star clusters may form. Here we investigate whether this mechanism can be inferred from observations of young clusters. To address this question, we have computed SPH simulations of the initial formation and evolution of a dynamically young star cluster through cold (sub-virial) collapse, starting with an ellipsoidal, turbulently seeded distribution of gas, and forming sink particles representing (proto)stars. While the initial density distributions of the clouds do not have large initial mass concentrations, gravitational focusing due to the global morphology leads to cluster formation. We use the resulting structures to extract observable morphological and kinematic signatures for the case of sub-virial collapse. We find that the signatures of the initial conditions can be erased rapidly as the gas and stars collapse, suggesting that kinematic observations need to be made either early in cluster formation and/or at larger scales, away from the growing cluster core. Our results emphasize that a dynamically young system is inherently evolving on short timescales, so that it can be highly misleading to use current-epoch conditions to study aspects such as star formation rates as a function of local density. Our simulations serve as a starting point for further studies of collapse including other factors such as magnetic fields and stellar feedback.

preprint2015arXiv

The number fraction of discs around brown dwarfs in Orion OB1a and the 25 Orionis group

We present a study of 15 new brown dwarfs belonging to the $\sim7$ Myr old 25 Orionis group and Orion OB1a sub-association with spectral types between M6 and M9 and estimated masses between $\sim0.07$M$_\odot$ and $\sim0.01$ M$_\odot$. By comparing them through a Bayesian method with low mass stars ($0.8\lesssim$ M/M$_\odot\lesssim0.1$) from previous works in the 25 Orionis group, we found statistically significant differences in the number fraction of classical T Tauri stars, weak T Tauri stars, class II, evolved discs and purely photospheric emitters at both sides of the sub-stellar mass limit. Particularly we found a fraction of $3.9^{+2.4}_{-1.6}~\%$ low mass stars classified as CTTS and class II or evolved discs, against a fraction of $33.3^{+10.8}_{-9.8}~\%$ in the sub-stellar mass domain. Our results support the suggested scenario in which the dissipation of discs is less efficient for decreasing mass of the central object.

preprint2014arXiv

A spectroscopic census in young stellar regions: the Sigma Orionis cluster

We present a spectroscopic survey of the stellar population of the Sigma Orionis cluster. We have obtained spectral types for 340 stars. Spectroscopic data for spectral typing come from several spectrographs with similar spectroscopic coverage and resolution. More than a half of stars of our sample are members confirmed by the presence of lithium in absorption, strong H$α$ in emission or weak gravity-sensitive features. In addition, we have obtained high resolution (R~34000) spectra in the H$α$ region for 169 stars in the region. Radial velocities were calculated from this data set. The radial velocity distribution for members of the cluster is in agreement with previous work. Analysis of the profile of the H$α$ line and infrared observations reveals two binary systems or fast rotators that mimic the H$α$ width expected in stars with accretion disks. On the other hand there are stars with optically thick disks and narrow H$α$ profile not expected in stars with accretion disks. This contribution constitutes the largest homogeneous spectroscopic data set of the Sigma Orionis cluster to date.

preprint2014arXiv

Accretion and Diffusion Timescales in Sheets and Filaments

A comparison of accretion and (turbulent) magnetic diffusion timescales for sheets and filaments demonstrates that dense star-forming clouds generally will -- under realistic conditions -- become supercritical due to mass accretion on timescales at least an order of magnitude shorter than ambipolar and/or turbulent diffusion timescales. Thus, ambipolar or turbulent diffusion -- while present -- is unlikely to control the formation of cores and stars.

preprint2014arXiv

Accretion Outbursts in Self-gravitating Protoplanetary Disks

We improve on our previous treatments of long-term evolution of protostellar disks by explicitly solving disk self-gravity in two dimensions. The current model is an extension of the one-dimensional layered accretion disk model of Bae et al. We find that gravitational instability (GI)-induced spiral density waves heat disks via compressional heating (i.e. $P\rm{d}V$ work), and can trigger accretion outbursts by activating the magnetorotational instability (MRI) in the magnetically inert disk dead-zone. The GI-induced spiral waves propagate well inside of gravitationally unstable region before they trigger outbursts at $R \lesssim 1$ AU where GI cannot be sustained. This long-range propagation of waves cannot be reproduced with the previously used local $α$ treatments for GI. In our standard model where zero dead-zone residual viscosity ($α_{\rm rd}$) is assumed, the GI-induced stress measured at the onset of outbursts is locally as large as $0.01$ in terms of the generic $α$ parameter. However, as suggested in our previous one-dimensional calculations, we confirm that the presence of a small but finite $α_{\rm rd}$ triggers thermally-driven bursts of accretion instead of the GI + MRI-driven outbursts that are observed when $α_{\rm rd}=0$. The inclusion of non-zero residual viscosity in the dead-zone decreases the importance of GI soon after mass feeding from the envelope cloud ceases. During the infall phase while the central protostar is still embedded, our models stay in a quiescent accretion phase with $\dot{M}_{acc}\sim10^{-8}-10^{-7}~M_{\odot}~{\rm yr^{-1}}$ over $60~\%$ of the time and spend less than $15~\%$ of the infall phase in accretion outbursts.

preprint2014arXiv

Characterizing the Youngest Herschel-detected Protostars I. Envelope Structure Revealed by CARMA Dust Continuum Observations

We present CARMA 2.9 mm dust continuum emission observations of a sample of 14 Herschel-detected Class 0 protostars in the Orion A and B molecular clouds, drawn from the PACS Bright Red Sources (PBRS) sample (Stutz et al.). These objects are characterized by very red 24 \micron\ to 70 \micron\ colors and prominent submillimeter emission, suggesting that they are very young Class 0 protostars embedded in dense envelopes. We detect all of the PBRS in 2.9 mm continuum emission and emission from 4 protostars and 1 starless core in the fields toward the PBRS; we also report 1 new PBRS source. The ratio of 2.9 mm luminosity to bolometric luminosity is higher by a factor of $\sim$5 on average, compared to other well-studied protostars in the Perseus and Ophiuchus clouds. The 2.9 mm visibility amplitudes for 6 of the 14 PBRS are very flat as a function of uv-distance, with more than 50\% of the source emission arising from radii $<$ 1500 AU. These flat visibility amplitudes are most consistent with spherically symmetric envelope density profiles with $ρ$~$\propto$~R$^{-2.5}$. Alternatively, there could be a massive unresolved structure like a disk or a high-density inner envelope departing from a smooth power-law. The large amount of mass on scales $<$ 1500 AU (implying high average central densities) leads us to suggest that that the PBRS with flat visibility amplitude profiles are the youngest PBRS and may be undergoing a brief phase of high mass infall/accretion and are possibly among the youngest Class 0 protostars. The PBRS with more rapidly declining visibility amplitudes still have large envelope masses, but could be slightly more evolved.

preprint2014arXiv

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

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

preprint2014arXiv

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

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

preprint2014arXiv

The Gould Belt Very Large Array Survey IV: The Taurus-Auriga complex

We present a multi-epoch radio study of the Taurus-Auriga star-forming complex made with the Karl G. Jansky Very Large Array at frequencies of 4.5 GHz and 7.5 GHz. We detect a total of 610 sources, 59 of which are related to young stellar objects and 18 to field stars. The properties of 56\% of the young stars are compatible with non-thermal radio emission. We also show that the radio emission of more evolved young stellar objects tends to be more non-thermal in origin and, in general, that their radio properties are compatible with those found in other star forming regions. By comparing our results with previously reported X-ray observations, we notice that young stellar objects in Taurus-Auriga follow a Güdel-Benz relation with $κ$=0.03, as we previously suggested for other regions of star formation. In general, young stellar objects in Taurus-Auriga and in all the previous studied regions seem to follow this relation with a dispersion of $\sim1$ dex. Finally, we propose that most of the remaining sources are related with extragalactic objects but provide a list of 46 unidentified radio sources whose radio properties are compatible with a YSO nature.

preprint2014arXiv

The Gould's Belt Very Large Array Survey III. The Orion region

We present results from a high-sensitivity (60 $μ$Jy), large-scale (2.26 square degree) survey obtained with the Karl G. Jansky Very Large Array as part of the Gould's Belt Survey program. We detected 374 and 354 sources at 4.5 and 7.5 GHz, respectively. Of these, 148 are associated with previously known Young Stellar Objects (YSOs). Another 86 sources previously unclassified at either optical or infrared wavelengths exhibit radio properties that are consistent with those of young stars. The overall properties of our sources at radio wavelengths such as their variability and radio to X-ray luminosity relation are consistent with previous results from the Gould's Belt Survey. Our detections provide target lists for followup VLBA radio observations to determine their distances as YSOs are located in regions of high nebulosity and extinction, making it difficult to measure optical parallaxes.

preprint2014arXiv

The low mass star and sub-stellar populations of the 25 Orionis group

We present the results of a survey of the low mass star and brown dwarf population of the 25 Orionis group. Using optical photometry from the CIDA Deep Survey of Orion, near IR photometry from the Visible and Infrared Survey Telescope for Astronomy and low resolution spectroscopy obtained with Hectospec at the MMT, we selected 1246 photometric candidates to low mass stars and brown dwarfs with estimated masses within $0.02 \lesssim M/M_\odot \lesssim 0.8$ and spectroscopically confirmed a sample of 77 low mass stars as new members of the cluster with a mean age of $\sim$7 Myr. We have obtained a system initial mass function of the group that can be well described by either a Kroupa power-law function with indices $α_3=-1.73\pm0.31$ and $α_2=0.68\pm0.41$ in the mass ranges $0.03\leq M/M_\odot\leq 0.08$ and $0.08\leq M/M_\odot\leq0.5$ respectively, or a Scalo log-normal function with coefficients $m_c=0.21^{+0.02}_{-0.02}$ and $σ=0.36\pm0.03$ in the mass range $0.03\leq M/M_\odot\leq0.8$. From the analysis of the spatial distribution of this numerous candidate sample, we have confirmed the East-West elongation of the 25 Orionis group observed in previous works, and rule out a possible southern extension of the group. We find that the spatial distributions of low mass stars and brown dwarfs in 25 Orionis are statistically indistinguishable. Finally, we found that the fraction of brown dwarfs showing IR excesses is higher than for low mass stars, supporting the scenario in which the evolution of circumstellar discs around the least massive objects could be more prolonged.

preprint2013arXiv

A Herschel and APEX Census of the Reddest Sources in Orion: Searching for the Youngest Protostars

We perform a census of the reddest, and potentially youngest, protostars in the Orion molecular clouds using data obtained with the PACS instrument onboard the Herschel Space Observatory and the LABOCA and SABOCA instruments on APEX as part of the Herschel Orion Protostar Survey (HOPS). A total of 55 new protostar candidates are detected at 70 um and 160 um that are either too faint (m24 > 7 mag) to be reliably classified as protostars or undetected in the Spitzer/MIPS 24 um band. We find that the 11 reddest protostar candidates with log (lambda F_lambda 70) / (lambda F_lambda 24) > 1.65 are free of contamination and can thus be reliably explained as protostars. The remaining 44 sources have less extreme 70/24 colors, fainter 70 um fluxes, and higher levels of contamination. Taking the previously known sample of Spitzer protostars and the new sample together, we find 18 sources that have log (lambda F_lambda 70) / (lambda F_lambda 24) > 1.65; we name these sources "PACS Bright Red sources", or PBRs. Our analysis reveals that the PBRs sample is composed of Class 0 like sources characterized by very red SEDs (T_bol < 45 K) and large values of sub-millimeter fluxes (L_smm/L_bol > 0.6%). Modified black-body fits to the SEDs provide lower limits to the envelope masses of 0.2 M_sun to 2 M_sun and luminosities of 0.7 L_sun to 10 L_sun. Based on these properties, and a comparison of the SEDs with radiative transfer models of protostars, we conclude that the PBRs are most likely extreme Class 0 objects distinguished by higher than typical envelope densities and hence, high mass infall rates.

preprint2013arXiv

HOPS 136: An Edge-On Orion Protostar Near the End of Envelope Infall

Edge-on protostars are valuable for understanding the disk and envelope properties of embedded young stellar objects, since the disk, envelope, and envelope cavities are all distinctly visible in resolved images and well constrained in modeling. Comparing 2MASS, WISE, Spitzer, Herschel, APEX, and IRAM photometry and limits from 1.25 to 1200 microns, Spitzer spectroscopy from 5 to 40 microns, and high-resolution Hubble imaging at 1.60 and 2.05 microns to radiative transfer modeling, we determine envelope and disk properties for the Class I protostar HOPS 136, an edge-on source in Orion's Lynds 1641 region. The source has a bolometric luminosity of 0.8 L_sun, a bolometric temperature of 170 K, and a ratio of submillimeter to bolometric luminosity of 0.8%. Via modeling, we find a total luminosity of 4.7 L_sun (larger than the observed luminosity due to extinction by the disk), an envelope mass of 0.06 M_sun, and a disk radius and mass of 450 AU and 0.002 M_sun. The stellar mass is highly uncertain but is estimated to fall between 0.4 and 0.5 M_sun. To reproduce the flux and wavelength of the near-infrared scattered-light peak in the spectral energy distribution, we require 5.4 x 10^-5 M_sun of gas and dust in each cavity. The disk has a large radius and a mass typical of more evolved T Tauri stars in spite of the significant remaining envelope. HOPS 136 appears to be a key link between the protostellar and optically revealed stages of star formation.

preprint2013arXiv

Long-term Evolution of Photoevaporating Protoplanetary Disks

We perform calculations of our one-dimensional, two-zone disk model to study the long-term evolution of the circumstellar disk. In particular, we adopt published photoevaporation prescriptions and examine whether the photoevaporative loss alone, coupled with a range of initial angular momenta of the protostellar cloud, can explain the observed decline of the frequency of optically-thick dusty disks with increasing age. In the parameter space we explore, disks have accreting and/or non-accreting transitional phases lasting of $\lesssim20 %$ of their lifetime, which is in reasonable agreement with observed statistics. Assuming that photoevaporation controls disk clearing, we find that initial angular momentum distribution of clouds needs to be weighted in favor of slowly rotating protostellar cloud cores. Again, assuming inner disk dispersal by photoevaporation, we conjecture that this skewed angular momentum distribution is a result of fragmentation into binary or multiple stellar systems in rapidly-rotating cores. Accreting and non-accreting transitional disks show different evolutionary paths on the $\dot{M}-R_{\rm wall}$ plane, which possibly explains the different observed properties between the two populations. However, we further find that scaling the photoevaporation rates downward by a factor of 10 makes it difficult to clear the disks on the observed timescales, showing that the precise value of the photoevaporative loss is crucial to setting the clearing times. While our results apply only to pure photoevaporative loss (plus disk accretion), there may be implications for models in which planets clear disks preferentially at radii of order 10 AU.

preprint2013arXiv

Modeling the Resolved Disk Around the Class 0 Protostar L1527

We present high-resolution sub/millimeter interferometric imaging of the Class 0 protostar L1527 IRS (IRAS 04368+2557) at 870 micron and 3.4 mm from the Submillimeter Array (SMA) and Combined Array for Research in Millimeter Astronomy (CARMA). We detect the signature of an edge-on disk surrounding the protostar with an observed diameter of 180 AU in the sub/millimeter images. The mass of the disk is estimated to be 0.007 M_sun, assuming optically thin, isothermal dust emission. The millimeter spectral index is observed to be quite shallow at all the spatial scales probed; alpha ~ 2, implying a dust opacity spectral index beta ~ 0. We model the emission from the disk and surrounding envelope using Monte Carlo radiative transfer codes, simultaneously fitting the sub/millimeter visibility amplitudes, sub/millimeter images, resolved L\arcmin\ image, spectral energy distribution, and mid-infrared spectrum. The best fitting model has a disk radius of R = 125 AU, is highly flared (H ~ R^1.3), has a radial density profile rho ~ R^-2.5, and has a mass of 0.0075 M_sun. The scale height at 100 AU is 48 AU, about a factor of two greater than vertical hydrostatic equilibrium. The resolved millimeter observations indicate that disks may grow rapidly throughout the Class 0 phase. The mass and radius of the young disk around L1527 is comparable to disks around pre-main sequence stars; however, the disk is considerably more vertically extended, possibly due to a combination of lower protostellar mass, infall onto the disk upper layers, and little settling of ~1 micron-sized dust grains.

preprint2013arXiv

Resolved Depletion Zones and Spatial Differentiation of N2H+ and N2D+

We present a study on the spatial distribution of N2D+ and N2H+ in thirteen protostellar systems. Eight of thirteen objects observed with the IRAM 30m telescope show relative offsets between the peak N2D+ (J=2-1) and N2H+ (J=1-0) emission. We highlight the case of L1157 using interferometric observations from the Submillimeter Array and Plateau de Bure Interferometer of the N2D+ (J=3-2) and N2H+ (J=1-0) transitions respectively. Depletion of N2D+ in L1157 is clearly observed inside a radius of ~2000 AU (7") and the N2H+ emission is resolved into two peaks at radii of ~1000 AU (3.5"), inside the depletion region of N2D+. Chemical models predict a depletion zone in N2H+ and N2D+ due to destruction of H2D+ at T ~ 20 K and the evaporation of CO off dust grains at the same temperature. However, the abundance offsets of 1000 AU between the two species are not reproduced by chemical models, including a model that follows the infall of the protostellar envelope. The average abundance ratios of N2D+ to N2H+ have been shown to decrease as protostars evolve by Emprechtinger et al., but this is the first time depletion zones of N2D+ have been spatially resolved. We suggest that the difference in depletion zone radii for N2H+ and N2D+ is caused by either the CO evaporation temperature being above 20 K or an H2 ortho-to-para ratio gradient in the inner envelope.

preprint2013arXiv

The Dependence of Star Formation Efficiency on Gas Surface Density

Studies by Lada (2010) and Heiderman (2010) have suggested that star formation mostly occurs above a threshold in gas surface density Sigma of Sigma_c = 120 Msun pc^{-2} (A_K = 0.8). Heiderman infer a threshold by combining low-mass star-forming regions, which show a steep increase in the star formation rate per unit area Sigma_SFR with increasing Sigma, and massive cores forming luminous stars which show a linear relation. We argue that these observations do not require a particular density threshold. The steep dependence of Sigma_SFR, approaching unity at protostellar core densities, is a natural result of the increasing importance of self-gravity at high densities along with the corresponding decrease in evolutionary timescales. The linear behavior of Sigma_SFR vs. Sigma in massive cores is consistent with probing dense gas in gravitational collapse, forming stars at a characteristic free-fall timescale given by the use of a particular molecular tracer. The low-mass and high-mass regions show different correlations between gas surface density and the area A spanned at that density, with A=Sigma^{-3} for low-mass regions and A=Sigma^{-1} for the massive cores; this difference, along with the use of differing techniques to measure gas surface density and star formation, suggests that connecting the low-mass regions with massive cores is problematic. We show that the approximately linear relationship between dense gas mass and stellar mass used by Lada similarly does not demand a particular threshold for star formation, and requires continuing formation of dense gas. Our results are consistent with molecular clouds forming by galactic hydrodynamic flows with subsequent gravitational collapse.

preprint2013arXiv

The Gould's Belt Very Large Array Survey I: The Ophiuchus complex

We present large-scale ($\sim$ 2000 square arcminutes), deep ($\sim$ 20 $μ$Jy), high-resolution ($\sim$ 1$''$) radio observations of the Ophiuchus star-forming complex obtained with the Karl G.\ Jansky Very Large Array at $λ$ = 4 and 6 cm. In total, 189 sources were detected, 56 of them associated with known young stellar sources, and 4 with known extragalactic objects; the other 129 remain unclassified, but most of them are most probably background quasars. The vast majority of the young stars detected at radio wavelengths have spectral types K or M, although we also detect 4 objects of A/F/B types and 2 brown dwarf candidates. At least half of these young stars are non-thermal (gyrosynchrotron) sources, with active coronas characterized by high levels of variability, negative spectral indices, and (in some cases) significant circular polarization. As expected, there is a clear tendency for the fraction of non-thermal sources to increase from the younger (Class 0/I or flat spectrum) to the more evolved (Class III or weak line T Tauri) stars. The young stars detected both in X-rays and at radio wavelengths broadly follow a Güdel-Benz relation, but with a different normalization than the most radio-active types of stars. Finally, we detect a $\sim$ 70 mJy compact extragalactic source near the center of the Ophiuchus core, which should be used as gain calibrator for any future radio observations of this region.

preprint2013arXiv

VLA and CARMA Observations of Protostars in the Cepheus Clouds: Sub-arcsecond Proto-Binaries Formed via Disk Fragmentation

We present observations of three Class 0/I protostars (L1157-mm, CB230 IRS1, and L1165-SMM1) using the Karl G. Jansky Very Large Array (VLA) and observations of two (L1165-SMM1 and CB230 IRS1) with the Combined Array for Research in Millimeter-wave Astronomy (CARMA). The VLA observations were taken at wavelengths of $λ= 7.3$ mm, 1.4 cm, 3.3 cm, 4.0 cm, and 6.5 cm with a best resolution of $\sim$0\farcs06 (18 AU) at 7.3 mm. The L1165-SMM1 CARMA observations were taken at $λ= 1.3$ mm with a best resolution of $\sim0\farcs3$ (100 AU), and the CB230 IRS1 observations were taken at $λ= 3.4$ mm with a best resolution of $\sim$3\arcsec\ (900 AU). We find that L1165-SMM1 and CB230 IRS1 have probable binary companions at separations of $\sim$0\farcs3 (100 AU) from detections of secondary peaks at multiple wavelengths. The position angles of these companions are nearly orthogonal to the direction of the observed bipolar outflows, consistent with the expected protostellar disk orientations. We suggest that these companions may have formed from disk fragmentation; turbulent fragmentation would not preferentially arrange the binary companions to be orthogonal to the outflow direction. For L1165-SMM1, both the 7.3 mm and 1.3 mm emission show evidence of a large (R $>$ 100 AU) disk. For the L1165-SMM1 primary protostar and the CB230 IRS1 secondary protostar, the 7.3 mm emission is resolved into structures consistent with $\sim20$ AU radius disks. For the other protostars, including L1157-mm, the emission is unresolved, suggesting disks with radii $< 20$ AU.

preprint2012arXiv

Complex Structure in Class 0 Protostellar Envelopes III: Velocity Gradients in Non-Axisymmetric Envelopes, Infall or Rotation?

We present an interferometric kinematic study of morphologically complex protostellar envelopes based on observations of the dense gas tracers N2H+ and NH3. The strong asymmetric nature of most envelopes in our sample leads us to question the common interpretation of velocity gradients as rotation, given the possibility of projection effects in the observed velocities. Several "idealized" sources with well-ordered velocity fields and envelope structures are now analyzed in more detail. We compare the interferometric data to position-velocity diagrams of kinematic models for spherical rotating collapse and filamentary rotating collapse. For this purpose, we developed a filamentary parametrization of the rotating collapse model to explore the effects of geometric projection on the observed velocity structures. We find that most envelopes in our sample have PV structures that can be reproduced by an infalling filamentary envelope projected at different angles within the plane of the sky. The infalling filament produces velocity shifts across the envelope that can mimic rotation, especially when viewed at single-dish resolutions and the axisymmetric rotating collapse model does not uniquely describe any dataset. Furthermore, if the velocities are assumed to reflect rotation, then the inferred centrifugal radii are quite large in most cases, indicating significant fragmentation potential or more likely another component to the line-center velocity. We conclude that ordered velocity gradients cannot be interpreted as rotation alone when envelopes are non-axisymmetric and that projected infall velocities likely dominate the velocity field on scales larger than 1000 AU.

preprint2012arXiv

Dust Filtration by Planet-Induced Gap Edges: Implications for Transitional Disks

By carrying out two-dimensional two-fluid global simulations, we have studied the response of dust to gap formation by a single planet in the gaseous component of a protoplanetary disk - the so-called "dust filtration" mechanism. We have found that a gap opened by a giant planet at 20 AU in a α=0.01, \dot{M}=10^{-8} Msun/yr disk can effectively stop dust particles larger than 0.1 mm drifting inwards, leaving a sub-millimeter dust cavity/hole. However, smaller particles are difficult to filter by a planet-induced gap due to 1) dust diffusion, and 2) a high gas accretion velocity at the gap edge. An analytic model is also derived to understand what size particles can be filtered by the gap edge. Finally, with our updated understanding of dust filtration, we have computed Monte-Carlo radiative transfer models with variable dust size distributions to generate the spectral energy distributions (SEDs) of disks with gaps. By comparing with transitional disk observations (e.g. GM Aur), we have found that dust filtration alone has difficulties to deplete small particles sufficiently to explain the near-IR deficit of transitional disks, except under some extreme circumstances. The scenario of gap opening by multiple planets studied previously suffers the same difficulty. One possible solution is by invoking both dust filtration and dust growth in the inner disk. In this scenario, a planet induced gap filters large dust particles in the disk, and the remaining small dust particles passing to the inner disk can grow efficiently without replenishment from fragmentation of large grains. Predictions for ALMA have also been made based on all these scenarios. We conclude that dust filtration with planet(s) in the disk is a promising mechanism to explain submm observations of transitional disks but it may need to be combined with other processes (e.g. dust growth) to explain the near-IR deficit.

preprint2012arXiv

On the nature of the Herbig B[e] star binary system V921 Scorpii: Discovery of a close companion and relation to the large-scale bipolar nebula

Belonging to the group of B[e] stars, V921 Scorpii is associated with a strong infrared excess and permitted and forbidden line emission, indicating the presence of low- and high-density circumstellar gas and dust. Many aspects of V921 Sco and other B[e] stars still remain mysterious, including their evolutionary state and the physical conditions resulting in the class-defining characteristics. In this paper, we employ VLTI/AMBER spectro-interferometry in order to reconstruct high-resolution (lambda/2B=0.0013") model-independent interferometric images for three wavelength bands around 1.65, 2.0, and 2.3 micrometer. In our images, we discover a close (25.0+/-0.8 milliarcsecond, corresponding to 29+/-0.9 AU at 1.15 kpc) companion around V921 Sco. Between two epochs in 2008 and 2009, we measure orbital motion of 7 degrees, implying an orbital period of about 35 years (for a circular orbit). Around the primary star, we detect a disk-like structure with indications for a radial temperature gradient. The polar axis of this AU-scale disk is aligned with the arcminute-scale bipolar nebula in which V921 Sco is embedded. Using Magellan/IMACS imaging, we detect multi-layered arc-shaped sub-structure in the nebula, suggesting episodic outflow activity from the system with a period of about 25 years, roughly matching the estimated orbital period of the companion. Our study supports the hypothesis that the B[e] phenomenon is related to dynamical interaction in a close binary system.

preprint2012arXiv

On the nature of the Herbig B[e] star binary system V921 Scorpii: Geometry and kinematics of the circumprimary disk on sub-AU scales

V921 Scorpii is a close binary system (separation 0.025") showing the B[e]-phenomenon. The system is surrounded by an enigmatic bipolar nebula, which might have been shaped by episodic mass-loss events, possibly triggered by dynamical interactions between the companion and the circumprimary disk (Kraus et al. 2012a). In this paper, we investigate the spatial structure and kinematics of the circumprimary disk, with the aim to obtain new insights into the still strongly debated evolutionary stage. For this purpose, we combine, for the first time, infrared spectro-interferometry (VLTI/AMBER, R=12,000) and spectro-astrometry (VLT/CRIRES, R=100,000), which allows us to study the AU-scale distribution of circumstellar gas and dust with an unprecedented velocity resolution of 3 km*s^-1. Using a model-independent photocenter analysis technique, we find that the Br-gamma-line emission rotates in the same plane as the dust disk. We can reproduce the wavelength-differential visibilities and phases and the double-peaked line profile using a Keplerian-rotating disk model. The derived mass of the central star is 5.4+/-0.4 M_sun*(d/1150 pc), which is considerably lower than expected from the spectral classification, suggesting that V921 Sco might be more distant (d approx 2kpc) than commonly assumed. Using the geometric information provided by our Br-gamma spectro-interferometric data and Paschen, Brackett, and Pfund line decrement measurements in 61 hydrogen recombination line transitions, we derive the density of the line-emitting gas (N_e=2...6*10^19 m^-3). Given that our measurements can be reproduced with a Keplerian velocity field without outflowing velocity component and the non-detection of age-indicating spectroscopic diagnostics, our study provides new evidence for the pre-main-sequence nature of V921 Sco.

preprint2012arXiv

The Low-mass Stellar Population in L1641: Evidence for Environmental Dependence of the Stellar Initial Mass Function

We present results from an optical photometric and spectroscopic survey of the young stellar population in L1641, the low-density star-forming region of the Orion A cloud south of the Orion Nebula Cluster (ONC). Our goal is to determine whether L1641 has a large enough low-mass population to make the known lack of high-mass stars a statistically-significant demonstration of environmental dependence of the upper mass stellar initial mass function (IMF). Our spectroscopic sample consists of IR-excess objects selected from the Spitzer/IRAC survey and non-excess objects selected from optical photometry. We have spectral confirmation of 864 members, with another 98 probable members; of the confirmed members, 406 have infrared excesses and 458 do not. Assuming the same ratio of stars with and without IR excesses in the highly-extincted regions, L1641 may contain as many as ~1600 stars down to ~0.1 solar mass, comparable within a factor of two to the the ONC. Compared to the standard models of the IMF, L1641 is deficient in O and early B stars to a 3-4 sigma significance level, assuming that we know of all the massive stars in L1641. With a forthcoming survey of the intermediate-mass stars, we will be in a better position to make a direct comparison with the neighboring, dense ONC, which should yield a stronger test of the dependence of the high-mass end of the stellar initial mass function upon environment.

preprint2012arXiv

Variable Accretion Outbursts in Protostellar Evolution

We extend the one-dimensional, two-zone models of long-term protostellar disk evolution with infall of Zhu et al. to consider the potential effects of a finite viscosity in regions where the ionization is too low for the magnetorotational instability (MRI) to operate (the "dead zone"). We find that the presence of a small but finite dead zone viscosity, as suggested by simulations of stratified disks with MRI-active outer layers, can trigger inside-out bursts of accretion, starting at or near the inner edge of the disk, instead of the previously-found outside-in bursts with zero dead zone viscosity, which originate at a few AU in radius. These inside-out bursts of accretion bear a qualitative resemblance to the outburst behavior of one FU Ori object, V1515 Cyg, in contrast to the outside-in burst models which more closely resemble the accretion events in FU Ori and V1057 Cyg. Our results suggest that the type and frequency of outbursts are potentially a probe of transport efficiency in the dead zone. Simulations must treat the inner disk regions, $R \lesssim 0.5$ AU, to show the detailed time evolution of accretion outbursts in general and to observe the inside-out bursts in particular.

preprint2011arXiv

A Spitzer-IRS Detection of Crystalline Silicates in a Protostellar Envelope

We present the Spitzer Space Telescope Infrared Spectrograph spectrum of the Orion A protostar HOPS-68. The mid-infrared spectrum reveals crystalline substructure at 11.1, 16.1, 18.8, 23.6, 27.9, and 33.6 microns superimposed on the broad 9.7 and 18 micron amorphous silicate features; the substructure is well matched by the presence of the olivine end-member forsterite. Crystalline silicates are often observed as infrared emission features around the circumstellar disks of Herbig Ae/Be stars and T Tauri stars. However, this is the first unambiguous detection of crystalline silicate absorption in a cold, infalling, protostellar envelope. We estimate the crystalline mass fraction along the line-of-sight by first assuming that the crystalline silicates are located in a cold absorbing screen and secondly by utilizing radiative transfer models. The resulting crystalline mass fractions of 0.14 and 0.17, respectively, are significantly greater than the upper limit found in the interstellar medium (< 0.02-0.05). We propose that the amorphous silicates were annealed within the hot inner disk and/or envelope regions and subsequently transported outward into the envelope by entrainment in a protostellar outflow

preprint2011arXiv

Challenges in Forming Planets by Gravitational Instability: Disk Irradiation and Clump Migration, Accretion & Tidal Destructio

We present two-dimensional hydrodynamic simulations of self-gravitating protostellar disks subject to axisymmetric infall from envelopes and irradiation from the central star, to explore disk fragmentation due to gravitational instability (GI), and the fragmented clump evolution. We assume that the disk is built gradually and smoothly by the infall, resulting in good numerical convergence. We confirm that for disks around solar-mass stars, infall at high rates at radii beyond ~50 AU leads to disk fragmentation. At lower infall rates <1e-5 Msun/yr, however, irradiation suppresses fragmentation. We find that, once formed, the fragments or clumps migrate inward on typical type-I time scales of ~2e3 yr initially, but later migration deviates from the type-I time scale when the clump becomes more massive than the local disk mass, and/or when they starts to open gaps. As they migrate, the clumps accrete from the disk at a rate 1e-3 to 1e-1 MJupiter/yr, consistent with analytic estimates that assume a 1-2 Hill radii cross section. Most clumps can grow to their isolation masses >0.1 Msun quickly. The eventual fates of these clumps, however, diverges depending on the migration speed: 3 out of 13 clumps become massive enough (brown dwarf mass) to open gaps in the disk and essentially stop migrating; 4 are tidally destroyed during inward migration; 6 migrate across the inner simulated disk boundary. A simple analytic model for clump evolution is derived to explain these different fates. Overall, our results indicate that fast migration, accretion, and tidal destruction of the clumps pose challenges to the scenario of giant planet formation by GI in situ, but may provide a formation mechanism for close binary systems.

preprint2011arXiv

Complex Structure in Class 0 Protostellar Envelopes II: Kinematic Structure from Single-Dish and Interferometric Molecular Line Mapping

We present a study of dense molecular gas kinematics in seventeen nearby protostellar systems using single-dish and interferometric molecular line observations. The non-axisymmetric envelopes around a sample of Class 0/I protostars were mapped in the N2H+ (J=1-0) tracer with the IRAM 30m, CARMA and PdBI as well as NH3 (1,1) with the VLA. The molecular line emission is used to construct line-center velocity and linewidth maps for all sources to examine the kinematic structure in the envelopes on spatial scales from 0.1 pc to ~1000 AU. The direction of the large-scale velocity gradients from single-dish mapping is within 45 degrees of normal to the outflow axis in more than half the sample. Furthermore, the velocity gradients are often quite substantial, the average being ~2.3 km\s\pc. The interferometric data often reveal small-scale velocity structure, departing from the more gradual large-scale velocity gradients. In some cases, this likely indicates accelerating infall and/or rotational spin-up in the inner envelope; the median velocity gradient from the interferometric data is ~10.7 km/s/pc. In two systems, we detect high-velocity HCO+ (J=1-0) emission inside the highest-velocity \nthp\ emission. This enables us to study the infall and rotation close to the disk and estimate the central object masses. The velocity fields observed on large and small-scales are more complex than would be expected from rotation alone, suggesting that complex envelope structure enables other dynamical processes (i.e. infall) to affect the velocity field.

preprint2011arXiv

Evolution of X-ray and FUV Disk-Dispersing Radiation Fields

We present new X-ray and Far Ultraviolet (FUV) observations of T Tauri stars covering the age range 1 to 10 Myr. Our goals are to observationally constrain the intensity of radiation fields responsible for evaporating gas from the circumstellar disk and to assess the feasibility of current photoevaporation models, focusing on X-ray and UV radiation. We greatly increase the number of 7-10 Myr old T Tauri stars observed in the X-rays by including observations of the well populated 25 Ori aggregate in the Orion OB1a subassociation. With these new 7-10 Myr objects, we confirm that X-ray emission remains constant from 1-10 Myr. We also show, for the first time, observational evidence for the evolution of FUV radiation fields with a sample of 56 accreting and non-accreting young stars spanning 1 Myr to 1 Gyr. We find that the FUV emission decreases on timescales consistent with the decline of accretion in classical T Tauri stars until reaching the chromospheric level in weak T Tauri stars and debris disks. Overall, we find that the observed strength of high energy radiation is consistent with that required by photoevaporation models to dissipate the disks in timescales of approximately 10 Myr. Finally, we find that the high energy fields that affect gas evolution are not similarly affecting dust evolution; in particular, we find that disks with inner clearings, the transitional disks, have similar levels of FUV emission as full disks.

preprint2011arXiv

On Rapid Disk Accretion and Initial Conditions in Protostellar Evolution

Low-mass protostars may accrete most of their material through short-lived episodes of rapid disk accretion; yet until recently evolutionary tracks for these protostars assumed only constant or slowly-varying accretion. Important initial steps toward examining the potential effects of rapid accretion were recently made by Baraffe, Chabrier, & Gallardo, who showed that in the limit of low-temperature ("cold") accretion, protostars may have much smaller radii than found in previous treatments. Such small radii at the end of protostellar accretion would have the effect of making some young stars appear much older - perhaps as much as 10 Myr - than they really are. However, we argue that very rapid disk accretion is unlikely to be cold, because observations of the best-studied pre-main sequence disks with rapid disk accretion outbursts - the FU Ori objects - have spectral energy distributions which imply large, not small, protostellar radii. In addition, theory indicates models that at high accretion rates, protostellar disks become internally hot and geometrically thick, making it much more likely that hot material is added to the star. In addition, the very large luminosity of the accretion disk is likely to irradiate the central star strongly, heating up the outer layers and potentially expanding them. Nevertheless, the Baraffe et al. calculations emphasize the importance of initial protostellar radii for subsequent evolution.

preprint2011arXiv

Rapid Star Formation and Global Gravitational Collapse

Most young stars in nearby molecular clouds have estimated ages of 1-2 Myr, suggesting that star formation is rapid. However, small numbers of stars in these regions with inferred ages of >= 5-10 Myr have been cited to argue that star formation is instead a slow, quasi-static process. When considering these alternative pictures it is important to recognize that the age spread in a given star-forming cloud is necessarily an upper limit to the timescales of local collapse, as not all spatially-distinct regions will start contracting at precisely the same instant. Moreover, star-forming clouds may dynamically evolve on timescales of a few Myr; in particular, global gravitational contraction will tend to yield increasing star formation rates with time due to generally increasing local gas densities. We show that two different numerical simulations of dynamic, flow-driven molecular cloud formation and evolution 1) predict age spreads for the main stellar population roughly consistent with observations, and 2) raise the possibility of forming small numbers of stars early in cloud evolution, before global contraction concentrates the gas and the bulk of the stellar population is produced. In general, the existence of a small number of older stars among a generally much-younger population is consistent with the picture of dynamic star formation, and may even provide clues to the time evolution of star-forming clouds.

preprint2010arXiv

Accretion in Evolved and Transitional Disks in Cep OB2: Looking for the Origin of the Inner Holes

We present accretion rates for a large number of solar-type stars in the Cep OB2 region, based on U band observations. Our study comprises 95 members of the ~4 Myr-old cluster Tr 37 (including 20 "transition" objects; TO), as well as the only CTTS in the ~12 Myr-old cluster NGC 7160. The stars show different disk morphologies, with the majority of them having evolved and flattened disks. The typical accretion rates are about one order of magnitude lower than in regions aged 1-2 Myr, and we find no strong correlation between disk morphology and accretion rates. Although half of the TO are not accreting, the median accretion rates of normal CTTS and accreting "transition" disks are similar (~3 10^{-9} and 2 10^{-9} Msun/yr, respectively). Comparison with other regions suggests that the TO observed at different ages do not necessarily represent the same type of objects, which is consistent with the fact that the different processes that can lead to reduced IR excess/inner disk clearing (e.g., binarity, dust coagulation/settling, photoevaporation, giant planet formation) do not operate on the same timescales. Accreting TO in Tr 37 are probably suffering strong dust coagulation/settling. Regarding the equally large number of non-accreting TO in the region, other processes, like photoevaporation, the presence of stellar/substellar companions, and/or giant planet formation may account for their "transitional" SEDs and negligible accretion rates.

preprint2010arXiv

Competitive Accretion in Sheet Geometry and the Stellar IMF

We report a set of numerical experiments aimed at addressing the applicability of competitive accretion to explain the high-mass end of the stellar initial mass function in a sheet geometry with shallow gravitational potential, in contrast to most previous simulations which have assumed formation in a cluster gravitational potential. Our flat cloud geometry is motivated by models of molecular cloud formation due to large-scale flows in the interstellar medium. The experiments consisted of SPH simulations of gas accretion onto sink particles formed rapidly from Jeans-unstable dense clumps placed randomly in the finite sheet. These simplifications allow us to study accretion with a minimum of free parameters, and to develop better statistics on the resulting mass spectra. We considered both clumps of equal mass and gaussian distributions of masses, and either uniform or spatially-varying gas densities. In all cases, the sink mass function develops a power law tail at high masses, with $dN/dlog M \propto M^{-Γ}$. The accretion rates of individual sinks follow $\dot{M} \propto M^2$ at high masses; this results in a continual flattening of the slope of the mass function towards an asymptotic form $Γ\sim 1$ (where the Salpeter slope is $Γ= 1.35$). The asymptotic limit is most rapidly reached when starting from a relatively broad distribution of initial sink masses. In general the resulting upper mass slope is correlated with the maximum sink mass; higher sink masses are found in simulations with flatter upper mass slopes. Although these simulations are of a highly idealized situation, the results suggest that competitive accretion may be relevant in a wider variety of environments than previously considered, and in particular that the upper mass distribution may generally evolve towards a limiting value of $Γ\sim 1$.

preprint2010arXiv

Complex Structure in Class 0 Protostellar Envelopes

We use archived IRAC images from the Spitzer Space Telescope to show that many Class 0 protostars exhibit complex, irregular, and non-axisymmetric structure within their dusty envelopes. Our 8 $μ$m extinction maps probe some of the densest regions in these protostellar envelopes. Many of the systems are observed to have highly irregular and non-axisymmetric morphologies on scales >= 1000 AU, with a quarter of the sample exhibiting filamentary or flattened dense structures. Complex envelope structure is observed in regions spatially distinct from outflow cavities, and the densest structures often show no systematic alignment perpendicular to the cavities. These results indicate that mass ejection is not responsible for much of the irregular morphologies we detect; rather, we suggest that the observed envelope complexity is mostly the result of collapse from protostellar cores with initially non-equilibrium structures. The striking non-axisymmetry in many envelopes could provide favorable conditions for the formation of binary systems. We also note that protostars in the sample appear to be formed preferentially near the edges of clouds or bends in filaments, suggesting formation by gravitational focusing.

preprint2010arXiv

Long-term Evolution of Protostellar and Protoplanetary Disks. I. Outbursts

As an initial investigation into the long-term evolution of protostellar disks, we explore the conditions required to explain the large outbursts of disk accretion seen in some young stellar objects. We use one-dimensional time-dependent disk models with a phenomenological treatment of the magnetorotational instability (MRI) and gravitational torques to follow disk evolution over long timescales. Comparison with our previous two-dimensional disk model calculations (Zhu et al. 2009b, Z2009b) indicates that the neglect of radial effects and two-dimensional disk structure in the one-dimensional case makes only modest differences in the results; this allows us to use the simpler models to explore parameter space efficiently. We find that the mass infall rates typically estimated for low-mass protostars generally result in AU-scale disk accretion outbursts, as predicted by our previous analysis (Zhu et al. 2009a,Z2009a). We also confirm quasi-steady accretion behavior for high mass infall rates if the values of $α$-parameter for the magnetorotational instability is small, while at this high accretion rate convection from the thermal instability may lead to some variations. We further constrain the combinations of the $α$-parameter and the MRI critical temperature, which can reproduce observed outburst behavior. Our results suggest that dust sublimation may be connected with full activation of the MRI. This is consistent with the idea that small dust captures ions and electrons to suppress the MRI. In a later paper we will explore both long-term outburst and disk evolution with this model, allowing for infall from protostellar envelopes with differing angular momenta.

preprint2010arXiv

Long-term Evolution of Protostellar and Protoplanetary Disks. II. Layered Accretion with Infall

We use one-dimensional two-zone time-dependent accretion disk models to study the long-term evolution of protostellar disks subject to mass addition from the collapse of a rotating cloud core. Our model consists of a constant surface density magnetically coupled active layer, with transport and dissipation in inactive regions only via gravitational instability. We start our simulations after a central protostar has formed, containing ~ 10% of the mass of the protostellar cloud. Subsequent evolution depends on the angular momentum of the accreting envelope. We find that disk accretion matches the infall rate early in the disk evolution because much of the inner disk is hot enough to couple to the magnetic field. Later infall reaches the disk beyond ~10 AU, and the disk undergoes outbursts of accretion in FU Ori-like events as described in Zhu et al. 2009c. If the initial cloud core is moderately rotating most of the central star's mass is built up by these outburst events. Our results suggest that the protostellar "luminosity problem" is eased by accretion during these FU Ori-like outbursts. After infall stops the disk enters the T Tauri phase. An outer, viscously evolving disk has structure that is in reasonable agreement with recent submillimeter studies and its surface density evolves from $Σ\propto R^{-1}$ to $R^{-1.5}$. An inner, massive belt of material-- the "dead zone" -- would not have been observed yet but should be seen in future high angular resolution observations by EVLA and ALMA. This high surface density belt is a generic consequence of low angular momentum transport efficiency at radii where the disk is magnetically decoupled, and would strongly affect planet formation and migration.

preprint2010arXiv

Morphological Complexity of Protostellar Envelopes

Extinction maps at 8 micron from the Spitzer Space Telescope show that many Class 0 protostars exhibit complex, irregular, and on-axisymmetric structure within the densest regions of their dusty envelopes. Many of the systems have highly irregular and on-axisymmetric morphologies on scales $\sim$1000 AU, with a quarter of the sample exhibiting filamentary or flattened dense structures. Complex envelope structure is observed in regions spatially distinct from outflow cavities, and the densest structures often show no systematic alignment perpendicular to the cavities. We suggest that the observed envelope complexity is the result of collapse from protostellar cores with initially non-equilibrium structures. The striking non-axisymmetry in many envelopes could provide favorable conditions for the formation of binary systems. We then show that the kinematics around L1165 as probed with N2H+ are indicative of asymmetric infall; the velocity gradient is not perpendicular to the outflow.

preprint2010arXiv

Spitzer Spectroscopy of the Circumprimary Disk in the Binary Brown Dwarf 2MASS J04414489+2301513

Using the Spitzer Infrared Spectrograph, we have performed mid-infrared spectroscopy on the young binary brown dwarf 2MASS J04414489+2301513 (15 AU) in the Taurus star-forming region. The spectrum exhibits excess continuum emission that likely arises from a circumstellar disk around the primary. Silicate emission is not detected in these data, indicating the presence of significant grain growth. This is one of the few brown dwarf disks at such a young age (~1 Myr) that has been found to lack silicate emission. To quantitatively constrain the properties of the disk, we have compared the spectral energy distribution of 2MASS J04414489+2301513 to the predictions of our vertical structure codes for irradiated accretion disks. Our models suggest that the remaining atmospheric grains of moderately depleted layers may have grown to a size of $\gtrsim5$ micron. In addition, our model fits indicate an outer radius of 0.2-0.3 AU for the disk. The small size of this circumprimary disk could be due to truncation by the secondary. The absence of an outer disk containing a reservoir of small, primordial grains, combined with a weak turbulent mechanism, may be responsible for the advanced grain growth in this disk.

preprint2010arXiv

The Inner Envelope and Disk of L1527 Revealed: Gemini L'-band Scattered Light Imaging

We present high-resolution L'-band imaging of the inner scattered light structure of Class 0 protostar L1527 IRS (IRAS 04368+2557) taken with the Gemini North telescope. The central point-source like feature seen in Spitzer Space Telescope IRAC images is resolved in the Gemini image into a compact bipolar structure with a narrow dark lane in the center. Two scattered light lobes are extended ~1.8" (200 AU) perpendicular to the direction of the outflow and ~2.5" (350 AU) along the outflow axis; the narrow dark lane between the scattered light lobes is ~0.45" (60 AU) thick. The observations are consistent with our initial modeling of a bright inner cavity separated by a dark lane due to extinction along the line of sight of the central protostar by the disk (Tobin et al. 2008). The bright, compact scattered light might be due to complex inner structure generated by the outflow, as suggested in our first paper, or it may more likely be the upper layers of the disk forming from infalling matter.

preprint2010arXiv

Transitional and Pre-Transitional disks: Gap Opening by Multiple Planets?

We use 2D hydrodynamic simulations of viscous disks to examine whether dynamically-interacting multiple giant planets can explain the large gaps (spanning over one order of magnitude in radius) inferred for the transitional and pre-transitional disks around T Tauri stars. In the absence of inner disk dust depletion, we find that it requires three to four giant planets to open up large enough gaps to be consistent with inferences from spectral energy distributions, because the gap width is limited by the tendency of the planets to be driven together into 2:1 resonances. With very strong tidal torques and/or rapid planetary accretion, fewer planets can also generate a large cavity interior to the locally formed gap(s) by preventing outer disk material from moving in. In these cases, however, the reduction of surface density produces a corresponding reduction in the inner disk accretion rate onto the star; this makes it difficult to explain the observed accretion rates of the pre/transitional disks. We find that even with four planets in disks, additional substantial dust depletion is required to explain observed disk gaps/holes. Substantial dust settling and growth, with consequent significant reductions in optical depths, is inferred for typical T Tauri disks in any case, and an earlier history of dust growth is consistent with the hypothesis that pre/transitional disks are explained by the presence of giant planets. We conclude that the depths and widths of gaps, and disk accretion rates in pre/transitional disks cannot be reproduced by a planet-induced gap opening scenario alone. Significant dust depletion is also required within the gaps/holes. Order of magnitude estimates suggest the mass of small dust particles (<1 micron) relative to the gas must be depleted to 1e-5 -- 1e-2 of the interstellar medium value, implying a very efficient mechanism of small dust removal or dust growth.

preprint2009arXiv

2-D simulations of FU Orionis disk outbursts

We have developed time-dependent models of FU Ori accretion outbursts to explore the physical properties of protostellar disks. Our two-dimensional, axisymmetric models incorporate full vertical structure with a new treatment of the radiative boundary condition for the disk photosphere. We find that FU Ori-type outbursts can be explained by a slow accumulation of matter due to gravitational instability. Eventually this triggers the magnetorotational instability, which leads to rapid accretion. The thermal instability is triggered in the inner disk but this instability is not necessary for the outburst. An accurate disk vertical structure, including convection, is important for understanding the outburst behavior. Large convective eddies develop during the high state in the inner disk. The models are in agreement with Spitzer IRS spectra and also with peak accretion rates and decay timescales of observed outbursts, though some objects show faster rise timescale. We also propose that convection may account for the observed mild-supersonic turbulence and the short-timescale variations of FU Orionis objects.

preprint2009arXiv

Far-Ultraviolet H2 Emission from Circumstellar Disks

We analyze the far-ultraviolet (FUV) spectra of 33 classical T Tauri stars (CTTS), including 20 new spectra obtained with the Advanced Camera for Surveys Solar Blind Channel (ACS/SBC) on the Hubble Space Telescope. Of the sources, 28 are in the ~1 Myr old Taurus-Auriga complex or Orion Molecular Cloud, 4 in the 8-10 Myr old Orion OB1a complex and one, TW Hya, in the 10 Myr old TW Hydrae Association. We also obtained FUV ACS/SBC spectra of 10 non-accreting sources surrounded by debris disks with ages between 10 and 125 Myr. We use a feature in the FUV spectra due mostly to electron impact excitation of \h2 to study the evolution of the gas in the inner disk. We find that the \h2 feature is absent in non-accreting sources, but is detected in the spectra of CTTS and correlates with accretion luminosity. Since all young stars have active chromospheres which produce strong X-ray and UV emission capable of exciting \h2 in the disk, the fact that the non-accreting sources show no \h2 emission implies that the \h2 gas in the inner disk has dissipated in the non-accreting sources, although dust (and possibly gas) remains at larger radii. Using the flux at 1600 Å, we estimate that the column density of \h2 left in the inner regions of the debris disks in our sample is less than ~ 3x10^-6 g cm^-2, nine orders of magnitude below the surface density of the minimum mass solar nebula at 1 AU.

preprint2009arXiv

Gravitational Collapse and Filament Formation: Comparison with the Pipe Nebula

Recent models of molecular cloud formation and evolution suggest that such clouds are dynamic and generally exhibit gravitational collapse. We present a simple analytic model of global collapse onto a filament and compare this with our numerical simulations of the flow-driven formation of an isolated molecular cloud to illustrate the supersonic motions and infall ram pressures expected in models of gravity-driven cloud evolution. We apply our results to observations of the Pipe Nebula, an especially suitable object for our purposes as its low star formation activity implies insignifcant perturbations from stellar feedback. We show that our collapsing cloud model can explain the magnitude of the velocity dispersions seen in the $^{13}$CO filamentary structure by Onishi et al. and the ram pressures required by Lada et al. to confine the lower-mass cores in the Pipe nebula. We further conjecture that higher-resolution simulations will show small velocity dispersions in the densest core gas, as observed, but which are infall motions and not supporting turbulence. Our results point out the inevitability of ram pressures as boundary conditions for molecular cloud filaments, and the possibility that especially lower-mass cores still can be accreting mass at significant rates, as suggested by observations.

preprint2009arXiv

Probing the Protostellar Envelope around L1157: the Dust and Gas Connection

We present observations of the Class 0 protostar L1157-mm using the Combined Array for Research in Millimeter-wave Astronomy (CARMA) in 3 mm dust continuum and N2H+ line emission. In the N2H+ line, we detect a large-scale envelope extended over a linear size of ~20,000AU flattened in the direction perpendicular to the outflow. This N2H+ feature coincides with the outer envelope seen in the 8 micron extinction by Looney et al. Meanwhile, the dust continuum traces the compact, nearly spherical structure of the inner envelope, where N2H+ becomes depleted. This highly flattened N2H+ envelope also shows dynamical signatures consistent with gravitational infall in the inner region, but a slow, solid-body rotation at large scales. This flattened structure is not a rotationally supported circumstellar disk; instead, it resembles a prestellar core both morphologically and kinematically, representing the early phase of a Class 0 system. In this paper, we construct a simple model to interpret both the dust continuum and N2H+ emission and suggest a possible dynamical scenario for the overall properties of the envelope.

preprint2008arXiv

Non-steady Accretion in Protostars

Observations indicate that mass accretion rates onto low-mass protostars are generally lower than the rates of infall to their disks; this suggests that much of the protostellar mass must be accreted during rare, short outbursts of rapid accretion. We explore when protostellar disk accretion is likely to be highly variable. While constant $α$ disks can in principle adjust their accretion rates to match infall rates, protostellar disks are unlikely to have constant $α$. In particular we show that neither models with angular momentum ransport due solely to the magnetorotational instability (MRI) nor ravitational instability (GI) are likely to transport disk mass at rotostellar infall rates over the large range of radii needed to move infalling envelope material down to the central protostar. We show that the MRI and GI are likely to combine to produce outbursts of rapid accretion starting at a few AU. Our analysis is consistent with the time-dependent models of Armitage, Livio, & Pringle (2001) and agrees with our observational study of the outbursting object FU Ori.

preprint2007arXiv

25 Orionis: A Kinematically Distinct 10 Myr Old Group in Orion OB1a

We report here on the photometric and kinematic properties of a well defined group of nearly 200 low-mass pre-main sequence stars, concentrated within ~ 1 deg of the early-B star 25 Ori, in the Orion OB1a sub-association. We refer to this stellar aggregate as the 25 Orionis group. The group also harbors the Herbig Ae/Be star V346 Ori and a dozen other early type stars with photometry, parallaxes, and some with IR excess emission, consistent with group membership. The number of high and low-mass stars is in agreement with expectations from a standard Initial Mass Function. The velocity distribution for the young stars in 25 Ori shows a narrow peak centered at 19.7 km/s, very close to the velocity of the star 25 Ori. Our results provide new and compelling evidence that the 25 Ori group is a distinct kinematic entity, and that considerable space and velocity structure is present in the Ori OB1a sub-association. The low-mass members follow a well defined band in the color-magnitude diagram, consistent with an isochronal age of ~ 7-10 Myr, depending on the assumed evolutionary model. The highest density of members is located near the star 25 Ori, but the actual extent of the cluster cannot be well constrained with our present data. In a simple-minded kinematic evolution scenario, the 25 Ori group may represent the evolved counterpart of a younger aggregate like the sigma Ori cluster. The 25 Ori stellar aggregate is the most populous ~ 10 Myr sample yet known within 500 pc, setting it as an excellent laboratory to study the evolution of solar-like stars and protoplanetary disks.

preprint2004arXiv

The Truncated Disk of CoKu Tau/4

We present a model of a dusty disk with an inner hole which accounts for the Spitzer Space Telescope Infrared Spectrograph observations of the low-mass pre-main sequence star CoKu Tau/4. We have modeled the mid-IR spectrum (between 8 and 25 mic) as arising from the inner wall of a disk. Our model disk has an evacuated inner zone of radius ~ 10 AU, with a dusty inner ``wall'', of half-height ~ 2 AU, that is illuminated at normal incidence by the central star. The radiative equilibrium temperature decreases from the inner disk edge outward through the optically-thick disk; this temperature gradient is responsible for the emission of the silicate bands at 10 and 20 mic. The observed spectrum is consistent with being produced by Fe-Mg amorphous glassy olivine and/or pyroxene, with no evidence of a crystalline component. The mid-infrared spectrum of CoKu Tau/4 is reminiscent of that of the much older star TW Hya, where it has been suggested that the significant clearing of its inner disk is due to planet formation. However, no inner disk remains in CoKu Tau/4, consistent with the star being a weak-emission (non-accreting) T Tauri star. The relative youth of CoKu Tau/4 (~ 1 Myr) may indicate much more rapid planet formation than typically assumed.

preprint1999arXiv

Turbulent Flow-Driven Molecular Cloud Formation: A Solution to the Post-T Tauri Problem?

We suggest that molecular clouds can be formed on short time scales by compressions from large scale streams in the interstellar medium (ISM). In particular, we argue that the Taurus-Auriga complex, with filaments of 10-20 pc $\times$ 2-5 pc, most have been formed by H I flows in $\lesssim 3$Myr, explaining the absence of post-T Tauri stars in the region with ages $\gtrsim 3$ Myr. Observations in the 21 cm line of the H I `halos' around the Taurus molecular gas show many features (broad asymmetric profiles, velocity shifts of H I relative to $^{12}$CO) predicted by our MHD numerical simulations, in which large-scale H I streams collide to produce dense filamentary structures. This rapid evolution is possible because the H I flows producing and disrupting the cloud have much higher velocities (5-10 kms) than present in the molecular gas resulting from the colliding flows. The simulations suggest that such flows can occur from the global ISM turbulence without requiring a single triggering event such as a SN explosion.

Lee Hartmann

What is connected

Connect this record

See the researcher in context

Building this map preview

66 published item(s)

The Structure and Kinematics of Three Class 0 Protostellar Jets from JWST

Anisotropic Infall and Substructure formation in Embedded Disks

Gemini-LIGHTS: Herbig Ae/Be and massive T-Tauri protoplanetary disks imaged with Gemini Planet Imager

The ODYSSEUS Survey. Motivation and First Results: Accretion, Ejection, and Disk Irradiation of CVSO 109

The Rate, Amplitude and Duration of Outbursts from Class 0 Protostars in Orion

A Study of Millimeter Variability in FUor Objects

A triple star system with a misaligned and warped circumstellar disk shaped by disk tearing

Angular Momenta, Magnetization, and Accretion of Protostellar Cores

The interpretation of protoplanetary disc wind diagnostic lines from X-ray photoevaporation and analytical MHD models

Irregular dust features around intermediate-mass young stars with GPI: signs of youth or misaligned disks?

CSI 2264: Characterizing Young Stars in NGC 2264 with Stochastically Varying Light Curves

Planetary Signatures in the SAO 206462 (HD 135344B) Disk: A Spiral Arm Passing Through Vortex?

Self-Destructing Spiral Waves: Global Simulations of a Spiral Wave Instability in Accretion Disks

Spectroscopic Binaries in the Orion Nebula Cluster and NGC 2264

The Spiral Wave Instability Induced by a Giant Planet: I. Particle Stirring in the Inner Regions of Protoplanetary Disks

Are protoplanetary disks born with vortices? -- Rossby wave instability driven by protostellar infall

CSI 2264: Characterizing Young Stars in NGC 2264 with Short-Duration, Periodic Flux Dips in their Light Curves

Kinematic and Spatial Substructure in NGC 2264

Signatures of star cluster formation by cold collapse

The number fraction of discs around brown dwarfs in Orion OB1a and the 25 Orionis group

A spectroscopic census in young stellar regions: the Sigma Orionis cluster

Accretion and Diffusion Timescales in Sheets and Filaments

Accretion Outbursts in Self-gravitating Protoplanetary Disks

Characterizing the Youngest Herschel-detected Protostars I. Envelope Structure Revealed by CARMA Dust Continuum Observations

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

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

The Gould Belt Very Large Array Survey IV: The Taurus-Auriga complex

The Gould's Belt Very Large Array Survey III. The Orion region

The low mass star and sub-stellar populations of the 25 Orionis group

A Herschel and APEX Census of the Reddest Sources in Orion: Searching for the Youngest Protostars

HOPS 136: An Edge-On Orion Protostar Near the End of Envelope Infall

Long-term Evolution of Photoevaporating Protoplanetary Disks

Modeling the Resolved Disk Around the Class 0 Protostar L1527

Resolved Depletion Zones and Spatial Differentiation of N2H+ and N2D+

The Dependence of Star Formation Efficiency on Gas Surface Density

The Gould's Belt Very Large Array Survey I: The Ophiuchus complex

VLA and CARMA Observations of Protostars in the Cepheus Clouds: Sub-arcsecond Proto-Binaries Formed via Disk Fragmentation

Complex Structure in Class 0 Protostellar Envelopes III: Velocity Gradients in Non-Axisymmetric Envelopes, Infall or Rotation?

Dust Filtration by Planet-Induced Gap Edges: Implications for Transitional Disks

On the nature of the Herbig B[e] star binary system V921 Scorpii: Discovery of a close companion and relation to the large-scale bipolar nebula

On the nature of the Herbig B[e] star binary system V921 Scorpii: Geometry and kinematics of the circumprimary disk on sub-AU scales

The Low-mass Stellar Population in L1641: Evidence for Environmental Dependence of the Stellar Initial Mass Function

Variable Accretion Outbursts in Protostellar Evolution

A Spitzer-IRS Detection of Crystalline Silicates in a Protostellar Envelope

Challenges in Forming Planets by Gravitational Instability: Disk Irradiation and Clump Migration, Accretion & Tidal Destructio

Complex Structure in Class 0 Protostellar Envelopes II: Kinematic Structure from Single-Dish and Interferometric Molecular Line Mapping

Evolution of X-ray and FUV Disk-Dispersing Radiation Fields

On Rapid Disk Accretion and Initial Conditions in Protostellar Evolution

Rapid Star Formation and Global Gravitational Collapse

Accretion in Evolved and Transitional Disks in Cep OB2: Looking for the Origin of the Inner Holes

Competitive Accretion in Sheet Geometry and the Stellar IMF

Complex Structure in Class 0 Protostellar Envelopes

Long-term Evolution of Protostellar and Protoplanetary Disks. I. Outbursts

Long-term Evolution of Protostellar and Protoplanetary Disks. II. Layered Accretion with Infall

Morphological Complexity of Protostellar Envelopes

Spitzer Spectroscopy of the Circumprimary Disk in the Binary Brown Dwarf 2MASS J04414489+2301513

The Inner Envelope and Disk of L1527 Revealed: Gemini L'-band Scattered Light Imaging

Transitional and Pre-Transitional disks: Gap Opening by Multiple Planets?

2-D simulations of FU Orionis disk outbursts

Far-Ultraviolet H2 Emission from Circumstellar Disks

Gravitational Collapse and Filament Formation: Comparison with the Pipe Nebula

Probing the Protostellar Envelope around L1157: the Dust and Gas Connection

Non-steady Accretion in Protostars

25 Orionis: A Kinematically Distinct 10 Myr Old Group in Orion OB1a

The Truncated Disk of CoKu Tau/4

Turbulent Flow-Driven Molecular Cloud Formation: A Solution to the Post-T Tauri Problem?