Showing 51–100 of 282 results for author: Hartmann, L

VLBA Observations of Strong Anisotripic Radio Scattering towards the Orion Nebula

Authors: Marina Kounkel, Lee Hartmann, Laurent Loinard, Amy J. Mioduszewski, Luis F. Rodríguez, Gisela N. Ortiz-León, Michael D. Johnson, Rosa M. Torres, Cesar Briceño

Abstract: We present observations of VLBA 20, a radio source found towards the edge of the Orion Nebula Cluster (ONC). Nonthermal emission dominates the spectral energy distribution of this object from radio to mid-infrared regime, suggesting that VLBA 20 is extragalactic. This source is heavily scattered in the radio regime. Very Long Baseline Array observations resolve it to ~34x19 mas at 5 GHz, and the w… ▽ More We present observations of VLBA 20, a radio source found towards the edge of the Orion Nebula Cluster (ONC). Nonthermal emission dominates the spectral energy distribution of this object from radio to mid-infrared regime, suggesting that VLBA 20 is extragalactic. This source is heavily scattered in the radio regime. Very Long Baseline Array observations resolve it to ~34x19 mas at 5 GHz, and the wavelength dependence of the scattering disk is consistent with $ν^{-2}$ at other frequencies. The origin of the scattering is most likely the ionized X-ray emitting gas from the winds of the most massive stars of the ONC. The scattering is highly anisotropic, with the axis ratio of 2:1, higher than what is typically observed towards other sources. △ Less

Submitted 4 April, 2018; originally announced April 2018.

Comments: 6 pages, accepted to AJ

How do T Tauri stars accrete

Authors: Lee Hartmann, Jaehan Bae

Abstract: We conjecture that observed protoplanetary disc accretion rates may be explained with low viscosities which could be the result of hydrodynamic turbulence. We show that viscosities parameterized in the usual way with $α\gtrsim 10^{-4}$, comparable to values suggested for hydrodynamic turbulence, can explain the observed accretion rates and lifetimes with plausible inner disc surface densities. Our… ▽ More We conjecture that observed protoplanetary disc accretion rates may be explained with low viscosities which could be the result of hydrodynamic turbulence. We show that viscosities parameterized in the usual way with $α\gtrsim 10^{-4}$, comparable to values suggested for hydrodynamic turbulence, can explain the observed accretion rates and lifetimes with plausible inner disc surface densities. Our models are also in better agreement with surface density estimates of the minimum mass solar nebula than models with rapid transport for a given mass accretion rate, such as recent models of accretion driven by magnetic winds. The required surface densities are a natural result of the protostellar infall phase, as long as non-gravitational transport is limited. We argue that, in addition to possible non-ideal magnetic transport due to disc winds possibly modified by the Hall effect, the effects of low-viscosity hydrodynamic accretion deserve more consideration. △ Less

Submitted 24 October, 2017; originally announced October 2017.

Comments: 8 pages, 4 figures, to appear in MNRAS

Kinematics and Structure of Star-forming Regions: Insights from Cold Collapse Models

Authors: Aleksandra Kuznetsova, Lee Hartmann, Javier Ballesteros-Paredes

Abstract: The origin of the observed morphological and kinematic substructure of young star forming regions is a matter of debate. We offer a new analysis of data from simulations of globally gravitationally collapsing clouds of progenitor gas to answer questions about sub-structured star formation in the context of cold collapse. As a specific example, we compare our models to recent radial velocity survey… ▽ More The origin of the observed morphological and kinematic substructure of young star forming regions is a matter of debate. We offer a new analysis of data from simulations of globally gravitationally collapsing clouds of progenitor gas to answer questions about sub-structured star formation in the context of cold collapse. As a specific example, we compare our models to recent radial velocity survey data from the IN-SYNC survey of Orion and new observations of dense gas kinematics, and offer possible interpretations of kinematic and morphological signatures in the region. In the context of our model, we find the frequently-observed hub-filament morphology of the gas naturally arises during gravitational evolution, as well as the dynamically-distinct kinematic substructure of stars. We emphasize that the global and not just the local gravitational potential plays an important role in determining the dynamics of both clusters and filaments. △ Less

Submitted 22 September, 2017; originally announced September 2017.

Comments: 7 pages, 4 figures, accepted to MNRAS

On the domain of Dirac and Laplace type operators on stratified spaces

Authors: Luiz Hartmann, Matthias Lesch, Boris Vertman

Abstract: We consider a generalized Dirac operator on a compact stratified space with an iterated cone-edge metric. Assuming a spectral Witt condition, we prove its essential self-adjointness and identify its domain and the domain of its square with weighted edge Sobolev spaces. This sharpens previous results where the minimal domain is shown only to be a subset of an intersection of weighted edge Sobolev s… ▽ More We consider a generalized Dirac operator on a compact stratified space with an iterated cone-edge metric. Assuming a spectral Witt condition, we prove its essential self-adjointness and identify its domain and the domain of its square with weighted edge Sobolev spaces. This sharpens previous results where the minimal domain is shown only to be a subset of an intersection of weighted edge Sobolev spaces. Our argument does not rely on microlocal techniques and is very explicit. The novelty of our approach is the use of an abstract functional analytic notion of interpolation scales. Our results hold for the Gauss-Bonnet and spin Dirac operators satisfying a spectral Witt condition. △ Less

Submitted 31 May, 2018; v1 submitted 5 September, 2017; originally announced September 2017.

Comments: revised; to appear in Journal of Spectral Theory

MSC Class: 35J75; 58j52

Journal ref: J. Spectr. Theory 8 (2018), no. 4, 1295-1348

Brasselet number and Newton polygons

Authors: Thais M. Dalbelo, Luiz Hartmann

Abstract: We present a formula to compute the Brasselet number of $f:(Y,0)\to (\C, 0)$ where $Y\subset X$ is a non-degenerate complete intersection in a toric variety $X$. As applications we establish several results concerning invariance of the Brasselet number for families of non-degenerate complete intersections. Moreover, when $(X,0) = (\mathbb{C}^n,0)$ we derive sufficient conditions to obtain the inva… ▽ More We present a formula to compute the Brasselet number of $f:(Y,0)\to (\C, 0)$ where $Y\subset X$ is a non-degenerate complete intersection in a toric variety $X$. As applications we establish several results concerning invariance of the Brasselet number for families of non-degenerate complete intersections. Moreover, when $(X,0) = (\mathbb{C}^n,0)$ we derive sufficient conditions to obtain the invariance of the Euler obstruction for families of complete intersections with an isolated singularity which are contained in $X$. △ Less

Submitted 18 October, 2018; v1 submitted 24 August, 2017; originally announced August 2017.

Disk evolution and the fate of water

Authors: Lee Hartmann, Fred Ciesla, Oliver Gressel, Richard Alexander

Abstract: We review the general theoretical concepts and observational constraints on the distribution and evolution of water vapor and ice in protoplanetary disks, with a focus on the Solar System. Water is expected to freeze out at distances greater than 1-3 AU from solar-type central stars; more precise estimates are difficult to obtain due to uncertainties in the complex processes involved in disk evolu… ▽ More We review the general theoretical concepts and observational constraints on the distribution and evolution of water vapor and ice in protoplanetary disks, with a focus on the Solar System. Water is expected to freeze out at distances greater than 1-3 AU from solar-type central stars; more precise estimates are difficult to obtain due to uncertainties in the complex processes involved in disk evolution, including dust growth, settling, and radial drift, and the level of turbulence and viscous dissipation within disks. Interferometric observations are now providing constraints on the positions of CO snow lines, but extrapolation to the unresolved regions where water ice sublimates will require much better theoretical understanding of mass and angular momentum transport in disks as well as more refined comparison of observations with sophisticated disk models. △ Less

Submitted 10 August, 2017; originally announced August 2017.

Comments: 26 pages, 8 figures: in the proceedings of the ISSI workshop on "The delivery of water to proto-planets, planets and satellites", to appear in Space Science Reviews

Are fibres in molecular cloud filaments real objects?

Authors: Manuel Zamora-Avilés, Javier Ballesteros-Paredes, Lee W. Hartmann

Abstract: We analyse the morphology and kinematics of dense filamentary structures produced in a numerical simulation of a star--forming cloud of $1.4 \times 10^4 \, \Msun$ evolving under their own self--gravity in a magnetized medium. This study is motivated by recent observations of velocity--coherent substructures ("fibres") in star-forming filaments. We find such "fibres" ubiquitously in our simulated f… ▽ More We analyse the morphology and kinematics of dense filamentary structures produced in a numerical simulation of a star--forming cloud of $1.4 \times 10^4 \, \Msun$ evolving under their own self--gravity in a magnetized medium. This study is motivated by recent observations of velocity--coherent substructures ("fibres") in star-forming filaments. We find such "fibres" ubiquitously in our simulated filament. We found that a fibre in one projection is not necessarily a fibre in another projection, and thus, caution should be taken into account when considering them as real objects. We found that only the densest parts of the filament ($\sim$30\% of the densest volume, which contains $\sim$70\% of the mass) belong to fibres in 2 projections. Moreover, it is quite common that they are formed by separated density enhancements superposed along the line of sight. Observations of fibres can yield insight into the level of turbulent substructure driven by gravity, but care should be taken in interpreting the results given the problem of line of sight superposition. We also studied the morphology and kinematics of the 3D skeleton (spine), finding that subfilaments accrete structured material mainly along the magnetic field lines, which are preferentially perpendicular to the skeleton. The magnetic field is at the same time dragged by the velocity field due to the gravitational collapse. △ Less

Submitted 12 August, 2017; v1 submitted 4 August, 2017; originally announced August 2017.

Comments: Accepted for publication in MNRAS (a minor correction of the second author affiliations)

Kinematics of the Optically Visible YSOs Toward the Orion B Molecular Cloud

Authors: Marina Kounkel, Lee Hartmann, Mario Mateo, John I. Bailey, III

Abstract: We present results from high resolution optical spectra toward 66 young stars in the Orion B molecular cloud to study their kinematics and other properties. Observations of the H$α$ and Li I 6707 Å lines are used to check membership and accretion properties. While the stellar radial velocities of in NGC 2068 and L1622 show good agreement with that of the molecular gas, many of the stars in NGC 202… ▽ More We present results from high resolution optical spectra toward 66 young stars in the Orion B molecular cloud to study their kinematics and other properties. Observations of the H$α$ and Li I 6707 Å lines are used to check membership and accretion properties. While the stellar radial velocities of in NGC 2068 and L1622 show good agreement with that of the molecular gas, many of the stars in NGC 2024 show a considerable offset. This could be a signature of either expansion of the cluster, high degree of the ejection of the stars from the cluster through the dynamical interaction, or the acceleration of the gas due to stellar feedback. △ Less

Submitted 4 July, 2017; originally announced July 2017.

Comments: Accepted to ApJ; 8 pages, 9 figures

On the formation of multiple concentric rings and gaps in protoplanetary disks

Authors: Jaehan Bae, Zhaohuan Zhu, Lee Hartmann

Abstract: As spiral waves driven by a planet in a gaseous disk steepen into a shock, they deposit angular momentum, opening a gap in the disk. This has been well studied using both linear theory and numerical simulations, but so far, only for the primary spiral arm -- the one directly attached to the planet. Using two-dimensional hydrodynamic simulations, we show that the secondary and tertiary arms driven… ▽ More As spiral waves driven by a planet in a gaseous disk steepen into a shock, they deposit angular momentum, opening a gap in the disk. This has been well studied using both linear theory and numerical simulations, but so far, only for the primary spiral arm -- the one directly attached to the planet. Using two-dimensional hydrodynamic simulations, we show that the secondary and tertiary arms driven by a planet can also open gaps as they steepen into shocks. The depths of the secondary/tertiary gaps in surface density grow with time in a low viscosity disk ($α= 5\times10^{-5}$), so even low-mass planets (e.g., super-Earth or mini-Neptune-mass) embedded in the disk can open multiple observable gaps, provided that sufficient time has passed. Applying our results to the HL Tau disk, we show that a single 30 Earth-mass planet embedded in the ring at 68.8~au (B5) can reasonably well reproduce the positions of the two major gaps at 13.2 and 32.3 au (D1 and D2), and roughly reproduce two other major gaps at 64.2 and 74.7 au (D5 and D6) seen in the mm continuum. The positions of secondary/tertiary gaps are found to be sensitive to the planetary mass and the disk temperature profile, so with accurate observational measurements of the temperature structure the positions of multiple gaps can be used to constrain the mass of the planet. We also comment on the gaps seen in the TW Hya and HD~163296 disk. △ Less

Submitted 28 October, 2017; v1 submitted 9 June, 2017; originally announced June 2017.

Comments: 11 pages, 6 figures, Accepted for publication in the ApJ

Characterizing the stellar population of NGC 1980

Authors: Marina Kounkel, Lee Hartmann, Nuria Calvet, Tom Megeath

Abstract: NGC 1980 is a young cluster that is located about 0.5 degrees south of the Orion Nebula Cluster (ONC). Recent studies by Bouy et al. and Pillitteri et al. have suggested that NGC 1980 contains an older population of stars compared to a much younger ONC, and that it belongs to a foreground population that may be located in front of the Orion A molecular gas by as much as 40 pc. In this work we pres… ▽ More NGC 1980 is a young cluster that is located about 0.5 degrees south of the Orion Nebula Cluster (ONC). Recent studies by Bouy et al. and Pillitteri et al. have suggested that NGC 1980 contains an older population of stars compared to a much younger ONC, and that it belongs to a foreground population that may be located in front of the Orion A molecular gas by as much as 40 pc. In this work we present low-resolution spectra towards 148 young stars found towards the NGC 1980 region. We determine the spectral types of these stars, examine accretion signatures and measure the extinction towards them. We determine that based on these observations, the age of the population of NGC 1980 is indistinguishable from L1641, estimated to be ~3 Myr, comparable with the study by Fang et al. △ Less

Submitted 22 May, 2017; originally announced May 2017.

Comments: Accepted for publication at AJ. 7 pages, 8 figures

The Herschel Orion Protostar Survey: Luminosity and Envelope Evolution

Authors: William J. Fischer, S. Thomas Megeath, Elise Furlan, Babar Ali, Amelia M. Stutz, John J. Tobin, Mayra Osorio, Thomas Stanke, P. Manoj, Charles A. Poteet, Joseph J. Booker, Lee Hartmann, Thomas L. Wilson, Philip C. Myers, Dan M. Watson

Abstract: The Herschel Orion Protostar Survey obtained well-sampled 1.2 - 870 micron spectral energy distributions (SEDs) of over 300 protostars in the Orion molecular clouds, home to most of the young stellar objects (YSOs) in the nearest 500 pc. We plot the bolometric luminosities and temperatures for 330 Orion YSOs, 315 of which have bolometric temperatures characteristic of protostars. The histogram of… ▽ More The Herschel Orion Protostar Survey obtained well-sampled 1.2 - 870 micron spectral energy distributions (SEDs) of over 300 protostars in the Orion molecular clouds, home to most of the young stellar objects (YSOs) in the nearest 500 pc. We plot the bolometric luminosities and temperatures for 330 Orion YSOs, 315 of which have bolometric temperatures characteristic of protostars. The histogram of bolometric temperature is roughly flat; 29% of the protostars are in Class 0. The median luminosity decreases by a factor of four with increasing bolometric temperature; consequently, the Class 0 protostars are systematically brighter than the Class I protostars, with a median luminosity of 2.3 L_sun as opposed to 0.87 L_sun. At a given bolometric temperature, the scatter in luminosities is three orders of magnitude. Using fits to the SEDs, we analyze how the luminosities corrected for inclination and foreground reddening relate to the mass in the inner 2500 AU of the best-fit model envelopes. The histogram of envelope mass is roughly flat, while the median corrected luminosity peaks at 15 L_sun for young envelopes and falls to 1.7 L_sun for late-stage protostars with remnant envelopes. The spread in luminosity at each envelope mass is three orders of magnitude. Envelope masses that decline exponentially with time explain the flat mass histogram and the decrease in luminosity, while the formation of a range of stellar masses explains the dispersion in luminosity. △ Less

Submitted 19 April, 2017; originally announced April 2017.

Comments: Accepted by ApJ, 16 pages, 14 figures

Polarized disk emission from Herbig Ae/Be stars observed using Gemini Planet Imager: HD 144432, HD 150193, HD 163296, and HD 169142

Authors: John D. Monnier, T. J. Harries, A. Aarnio, F. Adams, S. Andrews, N. Calvet, C. Espaillat, L. Hartmann, S. Hinkley, S. Kraus, M. McClure, R. Oppenheimer, M. Perrin, D. Wilner

Abstract: In order to look for signs of on-going planet formation in young disks, we carried out the first J-band polarized emission imaging of the Herbig Ae/Be stars HD 150193, HD 163296, and HD 169142 using the Gemini Planet Imager (GPI), along with new H band observations of HD 144432. We confirm the complex "double ring" structure for the nearly face-on system HD 169142 first seen in H-band, finding the… ▽ More In order to look for signs of on-going planet formation in young disks, we carried out the first J-band polarized emission imaging of the Herbig Ae/Be stars HD 150193, HD 163296, and HD 169142 using the Gemini Planet Imager (GPI), along with new H band observations of HD 144432. We confirm the complex "double ring" structure for the nearly face-on system HD 169142 first seen in H-band, finding the outer ring to be substantially redder than the inner one in polarized intensity. Using radiative transfer modeling, we developed a physical model that explains the full spectral energy distribution (SED) and J- and H-band surface brightness profiles, suggesting that the differential color of the two rings could come from reddened starlight traversing the inner wall and may not require differences in grain properties. In addition, we clearly detect an elongated, off-center ring in HD 163296 (MWC 275), locating the scattering surface to be 18 AU above the midplane at a radial distance of 77 AU, co-spatial with a ring seen at 1.3mm by ALMA linked to the CO snow line. Lastly, we report a weak tentative detection of scattered light for HD 150193 (MWC 863) and a non-detection for HD 144432; the stellar companion known for each of these targets has likely disrupted the material in the outer disk of the primary star. For HD 163296 and HD 169142, the prominent outer rings we detect could be evidence for giant planet formation in the outer disk or a manifestation of large-scale dust growth processes possibly related to snow-line chemistry. △ Less

Submitted 15 February, 2017; originally announced February 2017.

Comments: Accepted Astrophysical Journal

Gravitational Focusing and the Star Cluster Initial Mass Function

Authors: Aleksandra Kuznetsova, Lee Hartmann, Andreas Burkert

Abstract: We discuss the possibility that gravitational focusing, is responsible for the power-law mass function of star clusters $N(\log M) \propto M^{-1}$. This power law can be produced asymptotically when the mass accretion rate of an object depends upon the mass of the accreting body as $\dot{M} \propto M^2$. While Bondi-Hoyle-Littleton accretion formally produces this dependence on mass in a uniform m… ▽ More We discuss the possibility that gravitational focusing, is responsible for the power-law mass function of star clusters $N(\log M) \propto M^{-1}$. This power law can be produced asymptotically when the mass accretion rate of an object depends upon the mass of the accreting body as $\dot{M} \propto M^2$. While Bondi-Hoyle-Littleton accretion formally produces this dependence on mass in a uniform medium, realistic environments are much more complicated. However, numerical simulations in SPH allowing for sink formation yield such an asymptotic power-law mass function. We perform pure N-body simulations to isolate the effects of gravity from those of gas physics and to show that clusters naturally result with the power-law mass distribution. We also consider the physical conditions necessary to produce clusters on appropriate timescales. Our results help support the idea that gravitationally-dominated accretion is the most likely mechanism for producing the cluster mass function. △ Less

Submitted 30 January, 2017; originally announced February 2017.

Comments: 6 pages; accepted by ApJ

The Gould's Belt Distances Survey (GOBELINS) I. Trigonometric parallax distances and depth of the Ophiuchus complex

Authors: Gisela N. Ortiz-León, Laurent Loinard, Marina A. Kounkel, Sergio A. Dzib, Amy J. Mioduszewski, Luis F. Rodríguez, Rosa M. Torres, Rosa A. González-Lópezlira, Gerardo Pech, Juana L. Rivera, Lee Hartmann, Andrew F. Boden, Neal J. Evans II, Cesar Briceño, John J. Tobin, Phillip A. B. Galli, Donald Gudehus

Abstract: We present the first results of the Gould's Belt Distances Survey (GOBELINS), a project aimed at measuring the proper motion and trigonometric parallax of a large sample of young stars in nearby regions using multi-epoch Very Long Baseline Array (VLBA) radio observations. Enough VLBA detections have now been obtained for 16 stellar systems in Ophiuchus to derive their parallax and proper motion. T… ▽ More We present the first results of the Gould's Belt Distances Survey (GOBELINS), a project aimed at measuring the proper motion and trigonometric parallax of a large sample of young stars in nearby regions using multi-epoch Very Long Baseline Array (VLBA) radio observations. Enough VLBA detections have now been obtained for 16 stellar systems in Ophiuchus to derive their parallax and proper motion. This leads to distance determinations for individual stars with an accuracy of 0.3 to a few percent. In addition, the orbits of 6 multiple systems were modelled by combining absolute positions with VLBA (and in some cases, near infrared) angular separations. Twelve stellar systems are located in the dark cloud Lynds 1688, the individual distances for this sample are highly consistent with one another, and yield a mean parallax for Lynds 1688 of $\varpi=7.28\pm0.06$ mas, corresponding to a distance $d=137.3\pm1.2$ pc. This represents an accuracy better than 1%. Three systems for which astrometric elements could be measured are located in the eastern streamer (Lynds 1689) and yield an estimate of $\varpi=6.79\pm0.16$ mas, corresponding to a distance $d=147.3\pm3.4$ pc. This suggests that the eastern streamer is located about 10 pc farther than the core, but this conclusion needs to be confirmed by observations (currently being collected) of additional sources in the eastern streamer. From the measured proper motions, we estimate the one-dimensional velocity dispersion in Lynds 1688 to be 2.8$\pm$1.8 and 3.0$\pm$2.0 ${\rm km~s}^{-1}$, in R.A. and DEC., respectively, these are larger than, but still consistent within $1σ$, with those found in other studies. △ Less

Submitted 19 November, 2016; originally announced November 2016.

Comments: Accepted for publication in ApJ

Radio Measurements of the stellar proper motions in the core of the Orion Nebula Cluster

Authors: Sergio A. Dzib, Laurent Loinard, Luis F. Rodríguez, Laura Gómez, Jan Forbrich, Karl M. Menten, Marina A. Kounkel, Amy J. Mioduszewski, Lee Hartmann, John J. Tobin, Juana L. Rivera

Abstract: Using multi-epoch VLA observations, covering a time baseline of 29.1 years, we have measured the proper motions of 88 young stars with compact radio emission in the core of the Orion Nebula Cluster (ONC) and the neighboring BN/KL region. Our work increases the number of young stars with measured proper motion at radio frequencies by a factor of 2.5 and enables us to perform a better statistical an… ▽ More Using multi-epoch VLA observations, covering a time baseline of 29.1 years, we have measured the proper motions of 88 young stars with compact radio emission in the core of the Orion Nebula Cluster (ONC) and the neighboring BN/KL region. Our work increases the number of young stars with measured proper motion at radio frequencies by a factor of 2.5 and enables us to perform a better statistical analysis of the kinematics of the region than was previously possible. Most stars (79 out of 88) have proper motions consistent with a Gaussian distribution centered on $\overline{μ_α\cosδ}=1.07\pm0.09\quad{\rm mas\,yr^{-1}}$, and $\overline{μ_δ}=-0.84\pm0.16\quad{\rm mas\,yr^{-1}}$, with velocity dispersions of $σ_α=1.08\pm0.07\quad{\rm mas\,\,yr^{-1}},$ $σ_δ=1.27\pm0.15\quad{\rm mas\,\,yr^{-1}}$. We looked for organized movements of these stars but found no clear indication of radial expansion/contraction or rotation. The remaining nine stars in our sample show peculiar proper motions that differ from the mean proper motions of the ONC by more than 3-$σ$. One of these stars, V 1326 Ori, could have been expelled from the Orion Trapezium 7,000 years ago. Two could be related to the multi-stellar disintegration in the BN/KL region, in addition to the previously known sources BN, I and n. The others either have high uncertainties (so their anomalous proper motions are not firmly established) or could be foreground objects. △ Less

Submitted 11 November, 2016; originally announced November 2016.

Comments: 15 pages, 5 figures, 2 tables, 1 appendix, accepted in The Astrophysical Journal

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

Authors: Jaehan Bae, Richard P. Nelson, Lee Hartmann

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

Submitted 26 October, 2016; originally announced October 2016.

Comments: Accepted for publication in the The Astrophysical Journal, 19 pages, 12 figures, 1 table

The Gould's Belt Distances Survey (GOBELINS) III. The distance to the Serpens/Aquila Molecular Complex

Authors: Gisela N. Ortiz-León, Sergio A. Dzib, Marina A. Kounkel, Laurent Loinard, Amy J. Mioduszewski, Luis F. Rodríguez, Rosa M. Torres, Gerardo Pech, Juana L. Rivera, Lee Hartmann, Andrew F. Boden, Neal J. Evans II, Cesar Briceño, John J. Tobin, Phillip A. B. Galli

Abstract: We report on new distances and proper motions to seven stars across the Serpens/Aquila complex. The observations were obtained as part of the Gould's Belt Distances Survey (GOBELINS) project between September 2013 and April 2016 with the Very Long Baseline Array (VLBA). One of our targets is the proto-Herbig AeBe object EC 95, which is a binary system embedded in the Serpens Core. For this system,… ▽ More We report on new distances and proper motions to seven stars across the Serpens/Aquila complex. The observations were obtained as part of the Gould's Belt Distances Survey (GOBELINS) project between September 2013 and April 2016 with the Very Long Baseline Array (VLBA). One of our targets is the proto-Herbig AeBe object EC 95, which is a binary system embedded in the Serpens Core. For this system, we combined the GOBELINS observations with previous VLBA data to cover a total period of ~8 years, and derive the orbital elements and an updated source distance. The individual distances to sources in the complex are fully consistent with each other, and the mean value corresponds to a distance of $436.0\pm9.2$~pc for the Serpens/W40 complex. Given this new evidence, we argue that Serpens Main, W40 and Serpens South are physically associated and form a single cloud structure. △ Less

Submitted 10 October, 2016; originally announced October 2016.

Comments: Accepted for publication in ApJ

The Gould's Belt Distances Survey (GOBELINS) II. Distances and structure towards the Orion Molecular Clouds

Authors: Marina Kounkel, Lee Hartmann, Laurent Loinard, Gisela N. Ortiz-León, Amy J. Mioduszewski, Luis F. Rodríguez, Sergio A. Dzib, Rosa M. Torres, Gerardo Pech, Phillip A. B. Galli, Juana L. Rivera, Andrew F. Boden, Neal J. Evans II, Cesar Briceño, John J. Tobin

Abstract: We present the results of the Gould's Belt Distances Survey (GOBELINS) of young star forming regions towards the Orion Molecular Cloud Complex. We detected 36 YSOs with the Very Large Baseline Array (VLBA), 27 of which have been observed in at least 3 epochs over the course of 2 years. At least half of these YSOs belong to multiple systems. We obtained parallax and proper motions towards these sta… ▽ More We present the results of the Gould's Belt Distances Survey (GOBELINS) of young star forming regions towards the Orion Molecular Cloud Complex. We detected 36 YSOs with the Very Large Baseline Array (VLBA), 27 of which have been observed in at least 3 epochs over the course of 2 years. At least half of these YSOs belong to multiple systems. We obtained parallax and proper motions towards these stars to study the structure and kinematics of the Complex. We measured a distance of 388$\pm$5 pc towards the Orion Nebula Cluster, 428$\pm$10 pc towards the southern portion L1641, 388$\pm$10 pc towards NGC 2068, and roughly $\sim$420 pc towards NGC 2024. Finally, we observed a strong degree of plasma radio scattering towards $λ$ Ori. △ Less

Submitted 4 October, 2016; v1 submitted 13 September, 2016; originally announced September 2016.

Comments: Accepted for publication to ApJ

Zeta-determinants of Sturm-Liouville operators with quadratic potentials at infinity

Authors: Luiz Hartmann, Matthias Lesch, Boris Vertman

Abstract: We consider Sturm-Liouville operators on a half line $[a,\infty), a>0$, with potentials that are growing at most quadratically at infinity. Such operators arise naturally in the analysis of hyperbolic manifolds, or more generally manifolds with cusps. We establish existence and a formula for the associated zeta-determinant in terms of the Wronski-determinant of a fundamental system of solutions ad… ▽ More We consider Sturm-Liouville operators on a half line $[a,\infty), a>0$, with potentials that are growing at most quadratically at infinity. Such operators arise naturally in the analysis of hyperbolic manifolds, or more generally manifolds with cusps. We establish existence and a formula for the associated zeta-determinant in terms of the Wronski-determinant of a fundamental system of solutions adapted to the boundary conditions. Despite being the natural objects in the context of hyperbolic geometry, spectral geometry of such operators has only recently been studied in the context of analytic torsion. △ Less

Submitted 5 December, 2016; v1 submitted 25 July, 2016; originally announced July 2016.

Comments: 31 pages

MSC Class: 58J52; 34B24

Journal ref: J. Differential Equations 262 (2017), no. 5, 3431-3465

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

Authors: Jaehan Bae, Richard P. Nelson, Lee Hartmann, Samuel Richard

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

Submitted 6 July, 2016; originally announced July 2016.

Comments: Accepted for publication in ApJ, 22 pages, 22 figures, 2 tables, some figures are reduced in size to meet the arxiv requirement

An HST Survey for 100-1000 AU Companions around Young Stellar Objects in the Orion Molecular Clouds: Evidence for Environmentally Dependent Multiplicity

Authors: M. Kounkel, S. T. Megeath, C. A. Poteet, W. J. Fischer, L. Hartmann

Abstract: We present a near-IR survey for the visual multiples in the Orion molecular clouds region at separations between 100 and 1000 AU. These data were acquired at 1.6~$μ$m with the NICMOS and WFC3 cameras on the Hubble Space Telescope. Additional photometry was obtained for some of the sources at 2.05~$μ$m with NICMOS and in the $L'$-band with NSFCAM2 on the IRTF. Towards 129 protostars and 197 pre-mai… ▽ More We present a near-IR survey for the visual multiples in the Orion molecular clouds region at separations between 100 and 1000 AU. These data were acquired at 1.6~$μ$m with the NICMOS and WFC3 cameras on the Hubble Space Telescope. Additional photometry was obtained for some of the sources at 2.05~$μ$m with NICMOS and in the $L'$-band with NSFCAM2 on the IRTF. Towards 129 protostars and 197 pre-main sequence stars with disks observed with WFC3, we detect 21 and 28 candidate companions between the projected separations of 100---1000 AU, of which less than 5 and 8, respectively, are chance line of sight coincidences. The resulting companion fraction ($CF$) after the correction for the line of sight contamination is 14.4$^{+1.1}_{-1.3}$% for protostars and 12.5$^{+1.2}_{-0.8}$% for the pre-main sequence stars. These values are similar to those found for main sequence stars, suggesting that there is little variation in the $CF$ with evolution, although several observational biases may mask a decrease in the $CF$ from protostars to the main sequence stars. After segregating the sample into two populations based on the surrounding surface density of YSOs, we find that the $CF$ in the high stellar density regions ($Σ_{YSO} > 45$~pc$^{-2}$) is approximately 50% higher than that found in the low stellar density regions ($Σ_{YSO} < 45$~pc$^{-2}$). We interpret this as evidence for the elevated formation of companions at 100 to 1000 AU in the denser environments of Orion. We discuss possible reasons for this elevated formation. △ Less

Submitted 24 February, 2016; originally announced February 2016.

Comments: 20 pages, 7 tables, 22 figures Accepted to ApJ

The Herschel Orion Protostar Survey: Spectral Energy Distributions and Fits Using a Grid of Protostellar Models

Authors: E. Furlan, W. J. Fischer, B. Ali, A. M. Stutz, T. Stanke, J. J. Tobin, S. T. Megeath, M. Osorio, L. Hartmann, N. Calvet, C. A. Poteet, J. Booker, P. Manoj, D. M. Watson, L. Allen

Abstract: We present key results from the Herschel Orion Protostar Survey (HOPS): spectral energy distributions (SEDs) and model fits of 330 young stellar objects, predominantly protostars, in the Orion molecular clouds. This is the largest sample of protostars studied in a single, nearby star-formation complex. With near-infrared photometry from 2MASS, mid- and far-infrared data from Spitzer and Herschel,… ▽ More We present key results from the Herschel Orion Protostar Survey (HOPS): spectral energy distributions (SEDs) and model fits of 330 young stellar objects, predominantly protostars, in the Orion molecular clouds. This is the largest sample of protostars studied in a single, nearby star-formation complex. With near-infrared photometry from 2MASS, mid- and far-infrared data from Spitzer and Herschel, and sub-millimeter photometry from APEX, our SEDs cover 1.2-870 $μ$m and sample the peak of the protostellar envelope emission at ~100 $μ$m. Using mid-IR spectral indices and bolometric temperatures, we classify our sample into 92 Class 0 protostars, 125 Class I protostars, 102 flat-spectrum sources, and 11 Class II pre-main-sequence stars. We implement a simple protostellar model (including a disk in an infalling envelope with outflow cavities) to generate a grid of 30400 model SEDs and use it to determine the best-fit model parameters for each protostar. We argue that far-IR data are essential for accurate constraints on protostellar envelope properties. We find that most protostars, and in particular the flat-spectrum sources, are well-fit. The median envelope density and median inclination angle decrease from Class 0 to Class I to flat-spectrum protostars, despite the broad range in best-fit parameters in each of the three categories. We also discuss degeneracies in our model parameters. Our results confirm that the different protostellar classes generally correspond to an evolutionary sequence with a decreasing envelope infall rate, but the inclination angle also plays a role in the appearance, and thus interpretation, of the SEDs. △ Less

Submitted 30 March, 2016; v1 submitted 23 February, 2016; originally announced February 2016.

Comments: Accepted by ApJS; 87 pages; 2 online-only figure sets; corrected typos and other minor fixes

Spectroscopic Binaries in the Orion Nebula Cluster and NGC 2264

Authors: Marina Kounkel, Lee Hartmann, John J. Tobin, Mario Mateo, John I. Bailey III., Meghin Spencer

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

Submitted 17 February, 2016; originally announced February 2016.

Comments: Accepted to ApJ 7 pages, 15 figures, 6 tables

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

Authors: Jaehan Bae, Zhaohuan Zhu, Lee Hartmann

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

Submitted 22 February, 2016; v1 submitted 19 January, 2016; originally announced January 2016.

Comments: proof version

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

Authors: John Stauffer, Ann Marie Cody, Luisa Rebull, Lynne A. Hillenbrand, Neal J. Turner, John Carpenter, Sean Carey, Susan Terebey, Maria Morales-Calderon, Silvia H. P. Alencar, Pauline McGinnis, Alana Sousa, Jerome Bouvier, Laura Venuti, Lee Hartmann, Nuria Calvet, Giusi Micela, Ettore Flaccomio, Inseok Song, Rob Gutermuth, David Barrado, Frederick J. Vrba, Kevin Covey, William Herbst, Edward Gillen , et al. (3 additional authors not shown)

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

Submitted 13 January, 2016; originally announced January 2016.

Comments: 36 pages, 19 figures. Accepted to AJ

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

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

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

Submitted 3 November, 2015; originally announced November 2015.

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

Signatures of star cluster formation by cold collapse

Authors: Aleksandra Kuznetsova, Lee Hartmann, Javier Ballesteros-Paredes

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

Submitted 13 November, 2015; v1 submitted 28 October, 2015; originally announced October 2015.

Comments: 12 pages, accepted to ApJ

Bondi-Hoyle-Littleton accretion and the upper mass stellar IMF

Authors: Javier Ballesteros-Paredes, Lee W. Hartmann, Nadia Perez-Goytia, Aleksandra Kuznetsova

Abstract: We report on a series of numerical simulations of gas clouds with self-gravity forming sink particles, adopting an isothermal equation of state to isolate the effects of gravity from thermal physics on the resulting sink mass distributions. Simulations starting with supersonic velocity fluctuations develop sink mass functions with a high-mass power-law tail $dN/d\log M \propto M^Γ$,… ▽ More We report on a series of numerical simulations of gas clouds with self-gravity forming sink particles, adopting an isothermal equation of state to isolate the effects of gravity from thermal physics on the resulting sink mass distributions. Simulations starting with supersonic velocity fluctuations develop sink mass functions with a high-mass power-law tail $dN/d\log M \propto M^Γ$, $Γ= -1 \pm 0.1$, independent of the initial Mach number of the velocity field. Similar results but with weaker statistical significance hold for a simulation starting with initial density fluctuations. This mass function power-law dependence agrees with the asymptotic limit found by Zinnecker assuming Bondi-Hoyle-Littleton (BHL) accretion, even though the mass accretion rates of individual sinks show significant departures from the predicted $\mdot \propto M^2$ behavior. While BHL accretion is not strictly applicable due to the complexity of the environment, we argue that the final mass functions are the result of a {\em relative} $M^2$ dependence resulting from gravitationally-focused accretion. Our simulations may show the power-law mass function particularly clearly compared with others because our adoption of an isothermal equation of state limits the effects of thermal physics in producing a broad initial fragmentation spectrum; $Γ\rightarrow -1$ is an asymptotic limit found only when sink masses grow well beyond their initial values. While we have purposely eliminated many additional physical processes (radiative transfer, feedback) which can affect the stellar mass function, our results emphasize the importance of gravitational focusing for massive star formation. △ Less

Submitted 8 June, 2015; originally announced June 2015.

Comments: Accepted by MNRAS. 10 pages, 14 figures

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

Authors: Juan José Downes, Carlos Román-Zúñiga, Javier Ballesteros-Paredes, Cecilia Mateu, César Briceño, Jesús Hernández, Monika G. Petr-Gotzens, Nuria Calvet, Lee Hartmann, Karina Mauco

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

Submitted 20 April, 2015; originally announced April 2015.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society. 15 pages, 9 figures and 5 tables

The Gould's Belt Very Large Array Survey II: The Serpens region

Authors: Gisela N. Ortiz-León, Laurent Loinard, Amy J. Mioduszewski, Sergio A. Dzib, Luis F. Rodríguez, Gerardo Pech, Juana Rivera, Rosa M. Torres, Andrew Boden, Lee W. Hartmann, Neal J. Evans, Cesar Briceño II, John J. Tobin, Marina A. Kounkel, Rosa A. González-Lópezlira

Abstract: We present deep ($\sim 17~μ$Jy) radio continuum observations of the Serpens molecular cloud, the Serpens south cluster, and the W40 region obtained using the Very Large Array in its A configuration. We detect a total of 146 sources, 29 of which are young stellar objects (YSOs), 2 are BV stars and 5 more are associated with phenomena related to YSOs. Based on their radio variability and spectral in… ▽ More We present deep ($\sim 17~μ$Jy) radio continuum observations of the Serpens molecular cloud, the Serpens south cluster, and the W40 region obtained using the Very Large Array in its A configuration. We detect a total of 146 sources, 29 of which are young stellar objects (YSOs), 2 are BV stars and 5 more are associated with phenomena related to YSOs. Based on their radio variability and spectral index, we propose that about 16 of the remaining 110 unclassified sources are also YSOs. For approximately 65% of the known YSOs detected here as radio sources, the emission is most likely non-thermal, and related to stellar coronal activity. As also recently observed in Ophiuchus, our sample of YSOs with X-ray counterparts lies below the fiducial Güdel & Benz relation. Finally, we analyze the proper motions of 9 sources in the W40 region. This allows us to better constrain the membership of the radio sources in the region. △ Less

Submitted 7 April, 2015; v1 submitted 12 March, 2015; originally announced March 2015.

Comments: Accepted in The Astrophysical Journal

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

Authors: Jaehan Bae, Lee Hartmann, Zhaohuan Zhu

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

Submitted 9 March, 2015; originally announced March 2015.

Comments: 14 pages, 14 figures (figures 3, 4, 9, 11, 13, 14 with reduced quality to meet arXiv requirements), 2 tables, accepted for publication in ApJ

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

Authors: John Stauffer, Ann Marie Cody, Pauline McGinnis, Luisa Rebull, Lynne A. Hillenbrand, Neal J. Turner, John Carpenter, Peter Plavchan, Sean Carey, Susan Terebey, María Morales-Calderón, Silvia H. P. Alencar, Jerome Bouvier, Laura Venuti, Lee Hartmann, Nuria Calvet, Giusi Micela, Ettore Flaccomio, Inseok Song, Rob Gutermuth, David Barrado, Frederick J. Vrba, Kevin Covey, Debbie Padgett, William Herbst , et al. (5 additional authors not shown)

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

Submitted 26 January, 2015; originally announced January 2015.

Comments: 38 pages emulateapj; 30 figures; 5 tables. Accepted to AJ

Kinematic and Spatial Substructure in NGC 2264

Authors: John J. Tobin, Lee Hartmann, Gabor Furesz, Wen-Hsin Hsu, Mario Mateo

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

Submitted 16 January, 2015; v1 submitted 13 January, 2015; originally announced January 2015.

Comments: Accepted to AJ. 38 pages, 5 Figures 3 Tables

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

Authors: Sergio A. Dzib, Laurent Loinard, Luis F. Rodríguez, Amy J. Mioduszewski, Gisela N. Ortiz-León, Marina A. Kounkel, Gerardo Pech, Juana L. Rivera, Rosa M. Torres, Andrew F. Boden, Lee Hartmann, Neal J. Evans II, Cesar Briceño, John Tobin

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

Submitted 19 December, 2014; originally announced December 2014.

Comments: 43 pages, 6 figures, accepted in The Astrophysical Journal

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

Authors: John J. Tobin, Amelia M. Stutz, S. Thomas Megeath, William J. Fischer, Thomas Henning, Sarah E. Ragan, Babar Ali, Thomas Stanke, P. Manoj, Nuria Calvet, Lee Hartmann

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

Submitted 4 November, 2014; originally announced November 2014.

Comments: Accepted to ApJ, 40 pages, 9 Figures, 4 Tables

Accretion Outbursts in Self-gravitating Protoplanetary Disks

Authors: Jaehan Bae, Lee Hartmann, Zhaohuan Zhu, Richard P. Nelson

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

Submitted 12 September, 2014; originally announced September 2014.

Comments: 16 pages, 17 figures, 1 table, Accepted for publication in ApJ

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

Authors: Jesús Hernández, Nuria Calvet, Alice Perez, Cesar Briceño, Lorenzo Olguin, Maria E. Contreras, Lee Hartmann, Lori Allen, Catherine Espaillat, Ramírez Hernan

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

Submitted 6 August, 2014; v1 submitted 1 August, 2014; originally announced August 2014.

Comments: 56 pages, 14 figures. Accepted in The Astrophysical Journal

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

Authors: Juan José Downes, César Briceño, Cecilia Mateu, Jesús Hernández, Anna Katherina Vivas, Nuria Calvet, Lee Hartmann, Monika G. Petr-Gotzens, Lori Allen

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

Submitted 31 July, 2014; originally announced August 2014.

Comments: Accepted for publication in MNRAS

Accretion and Diffusion Timescales in Sheets and Filaments

Authors: Fabian Heitsch, Lee Hartmann

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

Submitted 9 June, 2014; originally announced June 2014.

Comments: 12 pages, 6 figures, accepted by MNRAS

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

Authors: Marina Kounkel, Lee Hartmann, Laurent Loinard, Amy J. Mioduszewski, Sergio A. Dzib, Gisela N. Ortiz-León, Luis F. Rodríguez, Gerardo Pech, Juana L. Rivera, Rosa M. Torres, Andrew F. Boden, Neal J. II Evans, Cesar Briceño, John Tobin

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

Submitted 6 June, 2014; originally announced June 2014.

Comments: Accepted for publication in ApJ; 51 pages, 15 figures, 5 tables

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

Authors: John Stauffer, Ann Marie Cody, Annie Baglin, Silvia H. P. Alencar, Luisa Rebull, Lynne A. Hillenbrand, Laura Venuti, Neal J. Turner, John Carpenter, Peter Plavchan, Krzysztof Findeisen, Sean Carey, Susan Terebey, María Morales-Calderón, Jerome Bouvier, Giusi Micela, Ettore Flaccomio, Inseok Song, Rob Gutermuth, Lee Hartmann, Nuria Calvet, Barbara Whitney, David Barrado, Frederick J. Vrba, Kevin Covey , et al. (3 additional authors not shown)

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

Submitted 25 January, 2014; originally announced January 2014.

Comments: Accepted for publication in AJ. 39 pages; 6 tables; 25 figures, many of which are highly degraded to meet size limits. Please download the regular resolution version at http://web.ipac.caltech.edu/staff/amc/staufferetal2014.pdf

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

Authors: Ann Marie Cody, John Stauffer, Annie Baglin, Giuseppina Micela, Luisa M. Rebull, Ettore Flaccomio, María Morales-Calderón, Suzanne Aigrain, Jèrôme Bouvier, Lynne A. Hillenbrand, Robert Gutermuth, Inseok Song, Neal Turner, Silvia H. P. Alencar, Konstanze Zwintz, Peter Plavchan, John Carpenter, Krzysztof Findeisen, Sean Carey, Susan Terebey, Lee Hartmann, Nuria Calvet, Paula Teixeira, Frederick J. Vrba, Scott Wolk , et al. (20 additional authors not shown)

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

Submitted 14 March, 2014; v1 submitted 25 January, 2014; originally announced January 2014.

Comments: Published in AJ. 59 pages; 4 tables; 49 figures, most of which are highly degraded to fit size limits. Author name typo corrected. For a better resolution version, please visit http://web.ipac.caltech.edu/staff/amc/codyetal2014.pdf

Journal ref: 2014 AJ, 147, 82

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

Authors: William J. Fischer, S. Thomas Megeath, John J. Tobin, Lee Hartmann, Amelia M. Stutz, Marina Kounkel, Charles A. Poteet, Babar Ali, Mayra Osorio, P. Manoj, Ian Remming, Thomas Stanke, Dan M. Watson

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

Submitted 20 December, 2013; originally announced December 2013.

Comments: Accepted by ApJ, 12 pages, 6 figures

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

Authors: John J. Tobin, Claire J. Chandler, David J. Wilner, Leslie W. Looney, Laurent Loinard, Hsin-Fang Chiang, Lee Hartmann, Nuria Calvet, Paola D'Alessio, Tyler L. Bourke, Woojin Kwon

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

Submitted 18 October, 2013; originally announced October 2013.

Comments: 44 pages, 13 Figures, 8 Tables, Accepted to ApJ

The analytic torsion of the finite metric cone over a compact manifold

Authors: L. Hartmann, M. Spreafico

Abstract: We give an explicit formula for the $L^2$ analytic torsion of the finite metric cone over an oriented compact connected Riemannian manifold. We provide an interpretation of the different factors appearing in this formula. We prove that the analytic torsion of the cone is the finite part of the limit obtained collapsing one of the boundaries, of the ratio of the analytic torsion of the frustum to a… ▽ More We give an explicit formula for the $L^2$ analytic torsion of the finite metric cone over an oriented compact connected Riemannian manifold. We provide an interpretation of the different factors appearing in this formula. We prove that the analytic torsion of the cone is the finite part of the limit obtained collapsing one of the boundaries, of the ratio of the analytic torsion of the frustum to a regularising factor. We show that the regularising factor comes from the set of the non square integrable eigenfunctions of the Laplace Beltrami operator on the cone. △ Less

Submitted 12 July, 2015; v1 submitted 25 August, 2013; originally announced August 2013.

Comments: To appear in Journal of the Mathematical Society of Japan

Curved walls: Grain growth, settling, and composition patterns in T Tauri disk dust sublimation fronts

Authors: M. K. McClure, P. D'Alessio, N. Calvet, C. Espaillat, L. Hartmann, B. Sargent, D. M. Watson, L. Ingleby, J. Hernandez

Abstract: The dust sublimation walls of disks around T Tauri stars represent a directly observable cross-section through the disk atmosphere and midplane. Their emission properties can probe the grain size distribution and composition of the innermost regions of the disk, where terrestrial planets form. Here we calculate the inner dust sublimation wall properties for four classical T Tauri stars with a narr… ▽ More The dust sublimation walls of disks around T Tauri stars represent a directly observable cross-section through the disk atmosphere and midplane. Their emission properties can probe the grain size distribution and composition of the innermost regions of the disk, where terrestrial planets form. Here we calculate the inner dust sublimation wall properties for four classical T Tauri stars with a narrow range of spectral types and inclination angles and a wide range of mass accretion rates to determine the extent to which the walls are radially curved. Best-fits to the near- and mid-IR excesses are found for curved, 2-layer walls in which the lower layer contains larger, hotter, amorphous pyroxene grains with Mg/(Mg+Fe)=0.6 and the upper layer contains submicron, cooler, mixed amorphous olivine and forsterite grains. As the mass accretion rates decrease from 10^(-8) to 10^(-10) Msol/yr, the maximum grain size in the lower layer decreases from 3 to 0.5 microns. We attribute this to a decrease in fragmentation and turbulent support for micron-sized grains with decreasing viscous heating. The atmosphere of these disks is depleted of dust with dust-gas mass ratios 1x10^(-4) of the ISM value, while the midplane is enhanced to 8 times the ISM value. For all accretion rates, the wall contributes at least half of the flux in the optically thin 10 micron silicate feature. Finally, we find evidence for an iron gradient in the disk, suggestive of that found in our solar system. △ Less

Submitted 19 August, 2013; v1 submitted 14 August, 2013; originally announced August 2013.

Comments: 18 pages, 12 figures, 1 table, accepted to ApJ

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

Authors: Sergio A. Dzib, Laurent Loinard, Amy J. Mioduszewski, Luis F. Rodríguez, Gisela N. Ortiz-León, Gerardo Pech, Juana L. Rivera, Rosa M. Torres, Andrew F. Boden, Lee Hartmann, Neal J. Evans II, Cesar Briceño, John Tobin

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

Submitted 18 July, 2013; originally announced July 2013.

Comments: Accepted for publication in The Astrophysical Journal. 31 pages, 8 figure, 5 tables

Long-term Evolution of Photoevaporating Protoplanetary Disks

Authors: Jaehan Bae, Lee Hartmann, Zhaohuan Zhu, Charles Gammie

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

Submitted 10 July, 2013; v1 submitted 9 July, 2013; originally announced July 2013.

Comments: 38 pages, 16 figures, 2 tables, Accepted for publication in ApJ

Modeling the Resolved Disk Around the Class 0 Protostar L1527

Authors: John J. Tobin, Lee Hartmann, Hsin-Fang Chiang, David J. Wilner, Leslie W. Looney, Laurent Loinard, Nuria Calvet, Paola D'Alessio

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

Submitted 15 May, 2013; originally announced May 2013.

Comments: 45 pages, 11 Figures, 7 Tables, Accepted to ApJ

Characterizing the stellar photospheres and near-infrared excesses in accreting T Tauri systems

Authors: M. K. McClure, N. Calvet, C. Espaillat, L. Hartmann, J. Hernandez, L. Ingleby, K. L. Luhman, P. D'Alessio, B. Sargent

Abstract: Using NASA IRTF SpeX data from 0.8 to 4.5 $μ$m, we determine self-consistently the stellar properties and excess emission above the photosphere for a sample of classical T Tauri stars (CTTS) in the Taurus molecular cloud with varying degrees of accretion. This process uses a combination of techniques from the recent literature as well as observations of weak-line T Tauri stars (WTTS) to account fo… ▽ More Using NASA IRTF SpeX data from 0.8 to 4.5 $μ$m, we determine self-consistently the stellar properties and excess emission above the photosphere for a sample of classical T Tauri stars (CTTS) in the Taurus molecular cloud with varying degrees of accretion. This process uses a combination of techniques from the recent literature as well as observations of weak-line T Tauri stars (WTTS) to account for the differences in surface gravity and chromospheric activity between the TTS and dwarfs, which are typically used as photospheric templates for CTTS. Our improved veiling and extinction estimates for our targets allow us to extract flux-calibrated spectra of the excess in the near-infrared. We find that we are able to produce an acceptable parametric fit to the near-infrared excesses using a combination of up to three blackbodies. In half of our sample, two blackbodies at temperatures of 8000 K and 1600 K suffice. These temperatures and the corresponding solid angles are consistent with emission from the accretion shock on the stellar surface and the inner dust sublimation rim of the disk, respectively. In contrast, the other half requires three blackbodies at 8000, 1800, and 800 K, to describe the excess. We interpret the combined two cooler blackbodies as the dust sublimation wall with either a contribution from the disk surface beyond the wall or curvature of the wall itself, neither of which should have single-temperature blackbody emission. In these fits, we find no evidence of a contribution from optically thick gas inside the inner dust rim. △ Less

Submitted 30 March, 2013; originally announced April 2013.

Comments: Accepted to ApJ, 22 pages, 2 appendices, 12 figures, 7 tables