Showing 1–50 of 243 results for author: Hartmann, L

JWST observations of $^{13}$CO$_{2}$ ice: Tracing the chemical environment and thermal history of ices in protostellar envelopes

Authors: Nashanty G. C. Brunken, Will R. M. Rocha, Ewine F. van Dishoeck, Robert Gutermuth, Himanshu Tyagi, Katerina Slavicinska, Pooneh Nazari, S. Thomas Megeath, Neal J. Evans II, Mayank Narang, P. Manoj, Adam E. Rubinstein, Dan M. Watson, Leslie W. Looney, Harold Linnartz, Alessio Caratti o Garatti, Henrik Beuther, Hendrik Linz, Pamela Klaassen, Charles A. Poteet, Samuel Federman, Guillem Anglada, Prabhani Atnagulov, Tyler L. Bourke, William J. Fischer , et al. (16 additional authors not shown)

Abstract: The structure and composition of simple ices can be modified during stellar evolution by protostellar heating. Key to understanding the involved processes are thermal and chemical tracers that can diagnose the history and environment of the ice. The 15.2 $μ$m bending mode of $^{12}$CO$_2$ has proven to be a valuable tracer of ice heating events but suffers from grain shape and size effects. A viab… ▽ More The structure and composition of simple ices can be modified during stellar evolution by protostellar heating. Key to understanding the involved processes are thermal and chemical tracers that can diagnose the history and environment of the ice. The 15.2 $μ$m bending mode of $^{12}$CO$_2$ has proven to be a valuable tracer of ice heating events but suffers from grain shape and size effects. A viable alternative tracer is the weaker $^{13}$CO$_2$ isotopologue band at 4.39 $μ$m which has now become accessible at high S/N with the $\textit{James Webb}$ Space Telescope (JWST). We present JWST NIRSpec observations of $^{13}$CO$_2$ ice in five deeply embedded Class 0 sources spanning a wide range in luminosities (0.2 - 10$^4$ L$_{\odot}$ ) taken as part of the Investigating Protostellar Accretion Across the Mass Spectrum (IPA) program. The band profiles vary significantly, with the most luminous sources showing a distinct narrow peak at 4.38 $μ$m. We first apply a phenomenological approach and show that a minimum of 3-4 Gaussian profiles are needed to fit the $^{13}$CO$_2$ absorption feature. We then combine these findings with laboratory data and show that a 15.2 $μ$m $^{12}$CO$_2$ band inspired five-component decomposition can be applied for the isotopologue band where each component is representative of CO$_2$ ice in a specific molecular environment. The final solution consists of cold mixtures of CO$_2$ with CH$_3$OH, H$_2$O and CO as well as segregated heated pure CO$_2$ ice. Our results are in agreement with previous studies of the $^{12}$CO$_2$ ice band, further confirming that $^{13}$CO$_{2}$ is a useful alternative tracer of protostellar heating events. We also propose an alternative solution consisting only of heated CO$_2$:CH$_3$OH and CO$_2$:H$_2$O ices and warm pure CO$_2$ ice for decomposing the ice profiles of the two most luminous sources in our sample. △ Less

Submitted 7 March, 2024; v1 submitted 6 February, 2024; originally announced February 2024.

Journal ref: A&A 685, A27 (2024)

IPA. Class 0 Protostars Viewed in CO Emission Using JWST/NIRSpec

Authors: Adam E. Rubinstein, Himanshu Tyagi, Pooneh Nazari, Robert Gutermuth, Samuel Federman, Mayank Narang, Will R. M. Rocha, Nashanty Brunken, Katie Slavicinska, Neal J. Evans II, Joel D. Green, Dan M. Watson, Henrik Beuther, Tyler Bourke, Alessio Caratti o Garatti, Lee Hartmann, Pamela Klaassen, Hendrik Linz, Leslie W. Looney, Puravankara Manoj, S. Thomas Megeath, James Muzerolle Page, Thomas Stanke, John J. Tobin, Ewine F. van Dishoeck , et al. (2 additional authors not shown)

Abstract: We investigate the bright CO fundamental emission in the central regions of five Class 0 protostars using the JWST's Near-Infrared Spectrograph (NIRSpec) and provide clues to what processes excite the gas. CO line emission images are extracted for a forest of $\sim$150 ro-vibrational transitions from two vibrational bands, $v=1-0$ and $v=2-1$. However, ${}^{13}$CO is not detected, and thus we can… ▽ More We investigate the bright CO fundamental emission in the central regions of five Class 0 protostars using the JWST's Near-Infrared Spectrograph (NIRSpec) and provide clues to what processes excite the gas. CO line emission images are extracted for a forest of $\sim$150 ro-vibrational transitions from two vibrational bands, $v=1-0$ and $v=2-1$. However, ${}^{13}$CO is not detected, and thus we can only statistically constrain the ${}^{12}$CO optical depth. Using noise measurements to determine upper limits to the ${}^{13}$CO emission, the flux ratio of ${}^{12}$CO/${}^{13}$CO indicates that the ${}^{12}$CO emission itself is not optically thick for ro-vibrational transitions with upper state rotational quantum number $J_u \geq 15$. We construct population diagrams to estimate the rotational temperature and number of molecules from extinction-corrected CO line fluxes assuming CO emission is optically thin. Two different temperature components are required for $v=1$ ($\sim600-1000$ K and $\sim1500-3500$ K), while one hotter component is required for $v=2$ ($\sim2000-6000$ K). The vibrational temperature is $\sim 900$ K among our sources and shows no trend with luminosity. Using vibrational temperatures and the inferred total amount of CO molecules for our sources, the total warm gas mass correlates strongly with luminosity ranging from $\sim$0.1 $\rm M_{Earth}$ for the low-mass protostars to $\sim$1 M$_{\rm sun}$ for the high-mass protostars. Interpreting the distribution of gas column densities and temperatures depends on radiative and chemical processes affecting CO. The presence of a $v=2$ population may indicate CO gas radiatively excited. Selective UV photodissociation of CO isotopologues around our high-mass sources may explain their depletion of ${}^{13}$CO. △ Less

Submitted 12 December, 2023; originally announced December 2023.

Comments: 30 pages, 6 figures, 5 tables, submitted to ApJ

Discovery of a collimated jet from the low luminosity protostar IRAS 16253$-$2429 in a quiescent accretion phase with the JWST

Authors: Mayank Narang, Manoj P., Himanshu Tyagi, Dan M. Watson, S. Thomas Megeath, Samuel Federman, Adam E. Rubinstein, Robert Gutermuth, Alessio Caratti o Garatti, Henrik Beuther, Tyler L. Bourke, Ewine F. Van Dishoeck, Neal J. Evans II, Guillem Anglada, Mayra Osorio, Thomas Stanke, James Muzerolle, Leslie W. Looney, Yao-Lun Yang, John J. Tobin, Pamela Klaassen, Nicole Karnath, Prabhani Atnagulov, Nashanty Brunken, William J. Fischer , et al. (14 additional authors not shown)

Abstract: Investigating Protostellar Accretion (IPA) is a JWST Cycle~1 GO program that uses NIRSpec IFU and MIRI MRS to obtain 2.9--28~$μ$m spectral cubes of young, deeply embedded protostars with luminosities of 0.2 to 10,000~L$_{\odot}$ and central masses of 0.15 to 12~M$_{\odot}$. In this Letter, we report the discovery of a highly collimated atomic jet from the Class~0 protostar IRAS~16253$-$2429, the l… ▽ More Investigating Protostellar Accretion (IPA) is a JWST Cycle~1 GO program that uses NIRSpec IFU and MIRI MRS to obtain 2.9--28~$μ$m spectral cubes of young, deeply embedded protostars with luminosities of 0.2 to 10,000~L$_{\odot}$ and central masses of 0.15 to 12~M$_{\odot}$. In this Letter, we report the discovery of a highly collimated atomic jet from the Class~0 protostar IRAS~16253$-$2429, the lowest luminosity source ($L_\mathrm{bol}$ = 0.2 $L_\odot$) in the IPA program. The collimated jet is detected in multiple [Fe~II] lines, [Ne~II], [Ni~II], and H~I lines, but not in molecular emission. The atomic jet has a velocity of about 169~$\pm$~15~km\,s$^{-1}$, after correcting for inclination. The width of the jet increases with distance from the central protostar from 23 to~60 au, corresponding to an opening angle of 2.6~$\pm$~0.5\arcdeg. By comparing the measured flux ratios of various fine structure lines to those predicted by simple shock models, we derive a shock {speed} of 54~km\,s$^{-1}$ and a preshock density of 2.0$\times10^{3}$~cm$^{-3}$ at the base of the jet. {From these quantities and using a suite of jet models and extinction laws we compute a mass loss rate between $0.4 -1.1\times10^{-10}~M_{\odot}$~yr~$^{-1}$.} The low mass loss rate is consistent with simultaneous measurements of low mass accretion rate ($2.4~\pm~0.8~\times~10^{-9}~M_{\odot}$~yr$^{-1}$) for IRAS~16253$-$2429 from JWST observations (Watson et al. in prep), indicating that the protostar is in a quiescent accretion phase. Our results demonstrate that very low-mass protostars can drive highly collimated, atomic jets, even during the quiescent phase. △ Less

Submitted 11 January, 2024; v1 submitted 21 October, 2023; originally announced October 2023.

Comments: Accepted to ApJL. Comments and feedback welcome

Investigating Protostellar Accretion-Driven Outflows Across the Mass Spectrum: JWST NIRSpec IFU 3-5~$μ$m Spectral Mapping of Five Young Protostars

Authors: Samuel Federman, S. Thomas Megeath, Adam E. Rubinstein, Robert Gutermuth, Mayank Narang, Himanshu Tyagi, P. Manoj, Guillem Anglada, Prabhani Atnagulov, Henrik Beuther, Tyler L. Bourke, Nashanty Brunken, Alessio Caratti o Garatti, Neal J. Evans II, William J. Fischer, Elise Furlan, Joel Green, Nolan Habel, Lee Hartmann, Nicole Karnath, Pamela Klaassen, Hendrik Linz, Leslie W. Looney, Mayra Osorio, James Muzerolle Page , et al. (13 additional authors not shown)

Abstract: Investigating Protostellar Accretion is a Cycle 1 JWST program using the NIRSpec+MIRI integral field units to obtain 2.9--28 $μ$m spectral cubes of five young protostars with luminosities of 0.2-10,000 L$_{\odot}$ in their primary accretion phase. This paper introduces the NIRSpec 2.9--5.3 $μ$m data of the inner 840-9000 au with spatial resolutions from 28-300 au. The spectra show rising continuum… ▽ More Investigating Protostellar Accretion is a Cycle 1 JWST program using the NIRSpec+MIRI integral field units to obtain 2.9--28 $μ$m spectral cubes of five young protostars with luminosities of 0.2-10,000 L$_{\odot}$ in their primary accretion phase. This paper introduces the NIRSpec 2.9--5.3 $μ$m data of the inner 840-9000 au with spatial resolutions from 28-300 au. The spectra show rising continuum emission; deep ice absorption; emission from H$_{2}$, H~I, and [Fe~II]; and the CO fundamental series in emission and absorption. Maps of the continuum emission show scattered light cavities for all five protostars. In the cavities, collimated jets are detected in [Fe~II] for the four $< 320$~L$_{\odot}$ protostars, two of which are additionally traced in Br-$α$. Knots of [Fe~II] emission are detected toward the most luminous protostar, and knots of [FeII] emission with dynamical times of $< 30$~yrs are found in the jets of the others. While only one jet is traced in H$_2$, knots of H$_2$ and CO are detected in the jets of four protostars. H$_2$ is seen extending through the cavities, showing that they are filled by warm molecular gas. Bright H$_2$ emission is seen along the walls of a single cavity, while in three cavities narrow shells of H$_2$ emission are found, one of which has an [Fe~II] knot at its apex. These data show cavities containing collimated jets traced in atomic/ionic gas surrounded by warm molecular gas in a wide-angle wind and/or gas accelerated by bow shocks in the jets. △ Less

Submitted 24 April, 2024; v1 submitted 5 October, 2023; originally announced October 2023.

Comments: 26 pages, 12 figures

Journal ref: ApJ 966 41 (2024)

Magnetically-activated accretion outbursts of pre-main sequence discs

Authors: Jacob Cleaver, Lee Hartmann, Jaehan Bae

Abstract: We investigate whether triggering of the magnetorotational instability (MRI) in protoplanetary discs can account for the wide diversity of observed accretion outbursts. We show that short-lived, relatively low accretion rate events probably result from triggering in the inner disc and can occur at low surface densities, comparable to or smaller than the minimum mass solar nebula, and thus are very… ▽ More We investigate whether triggering of the magnetorotational instability (MRI) in protoplanetary discs can account for the wide diversity of observed accretion outbursts. We show that short-lived, relatively low accretion rate events probably result from triggering in the inner disc and can occur at low surface densities, comparable to or smaller than the minimum mass solar nebula, and thus are very unlikely to result from MRI triggering by gravitational instability. We develop time-dependent accretion disc models using an $α$-viscosity approach and calculate light curves to compare with observations. Our modeling indicates that the lag time between infrared and optical bursts seen in Gaia 17bpi can be explained with an outside-in propagation with an $α\sim 0.1$ in the MRI-active region, consistent with other estimates. While outbursts in inner discs can show time delays of a few years between infrared and optical light curves, our models indicate that large, FU Ori-like bursts can exhibit infrared precursors decades before optical bursts. Detecting such precursors could enable analysis of the central star before it is overwhelmed by the rapid accreting material, as well as constraining outburst physics. Our results emphasize the importance of near-infrared monitoring of young stellar objects in addition to optical surveys. In addition, our findings emphasize the need for more sophisticated, three-dimensional, non-ideal magnetohydrodynamic simulations to fully exploit observational results. △ Less

Submitted 19 June, 2023; originally announced June 2023.

Comments: 13 pages, 20 figures, accepted by MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, stad1784, June 14 2023

Constraints on star formation in Orion A from Gaia

Authors: Jairo A. Alzate, Gustavo Bruzual, Marina Kounkel, Gladis Magris, Lee Hartmann, Nuria Calvet, Lyra Cao

Abstract: We develop statistical methods within a Bayesian framework to infer the star formation history from photometric surveys of pre-main sequence populations. Our procedures include correcting for biases due to extinction in magnitude-limited surveys, and using distributions from subsets of stars with individual extinction measurements. We also make modest corrections for unresolved binaries. We apply… ▽ More We develop statistical methods within a Bayesian framework to infer the star formation history from photometric surveys of pre-main sequence populations. Our procedures include correcting for biases due to extinction in magnitude-limited surveys, and using distributions from subsets of stars with individual extinction measurements. We also make modest corrections for unresolved binaries. We apply our methods to samples of populations with Gaia photometry in the Orion A molecular cloud. Using two well-established sets of evolutionary tracks, we find that, although our sample is incomplete at youngest ages due to extinction, star formation has proceeded in Orion A at a relatively constant rate between ages of about 0.3 and 5 Myr, in contrast to other studies suggesting multiple epochs of star formation. Similar results are obtained for a set of tracks that attempt to take the effects of strong magnetic fields into account. We also find no evidence for a well-constrained "birthline" that would result from low-mass stars appearing first along the deuterium-burning main sequence, especially using the magnetic evolutionary tracks. While our methods have been developed to deal with Gaia data, they may be useful for analyzing other photometric surveys of star-forming regions. △ Less

Submitted 19 May, 2023; originally announced May 2023.

Report number: 2305.1182

300: An ACA 870 $μ$m Continuum Survey of Orion Protostars and their Evolution

Authors: Samuel Federman, S. Thomas Megeath, John J. Tobin, Patrick D. Sheehan, Riwaj Pokhrel, Nolan Habel, Amelia M. Stutz, William J. Fischer, Lee Hartmann, Thomas Stanke, Mayank Narang, Mayra Osorio, Prabhani Atnagulov, Rohan Rahatgaonkar

Abstract: We present an 870 $μ$m continuum survey of 300 protostars from the Herschel Orion Protostar Survey using the Atacama Compact Array (ACA). These data measure protostellar flux densities on envelope scales $\leq$8000 au (20") and resolve the structure of envelopes with 1600 au (4") resolution, a factor of 3-5 improvement in angular resolution over existing single-dish 870 $μ$m observations. We compa… ▽ More We present an 870 $μ$m continuum survey of 300 protostars from the Herschel Orion Protostar Survey using the Atacama Compact Array (ACA). These data measure protostellar flux densities on envelope scales $\leq$8000 au (20") and resolve the structure of envelopes with 1600 au (4") resolution, a factor of 3-5 improvement in angular resolution over existing single-dish 870 $μ$m observations. We compare the ACA observations to Atacama Large Millimeter/submillimeter Array 12 m array observations at 870 $μ$m with $\sim$0.1 (40 au) resolution. Using the 12 m data to measure the fluxes from disks and the ACA data within 2500 au to measure the combined disk plus envelope fluxes, we calculate the 12 m/ACA 870 $μ$m flux ratios. Our sample shows a clear evolution in this ratio. Class 0 protostars are mostly envelope-dominated with ratios $<$0.5. In contrast, Flat Spectrum protostars are primarily disk-dominated with ratios near 1, although with a number of face-on protostars dominated by their envelopes. Class I protostars span the range from envelope to disk-dominated. The increase in ratio is accompanied by a decrease in the envelope fluxes and estimated mass infall rates. We estimate that 80$\%$ of the mass is accreted during the envelope-dominated phase. We find that the 12 m/ACA flux ratio is an evolutionary indicator that largely avoids the inclination and foreground extinction dependence of spectral energy distribution-based indicators. △ Less

Submitted 23 February, 2023; v1 submitted 14 October, 2022; originally announced October 2022.

Journal ref: ApJ 944 49 (2023)

Extension of HOPS Out to 500 ParSecs (eHOPS). I. Identification and Modeling of Protostars in the Aquila Molecular Clouds

Authors: Riwaj Pokhrel, S. Thomas Megeath, Robert A. Gutermuth, Elise Furlan, William J. Fischer, Samuel Federman, John J. Tobin, Amelia M. Stutz, Lee Hartmann, Mayra Osorio, Dan M. Watson, Thomas Stanke, P. Manoj, Mayank Narang, Prabhani Atnagulov, Nolan Habel, Wafa Zakri

Abstract: We present a Spitzer/Herschel focused survey of the Aquila molecular clouds ($d \sim 436$~pc) as part of the eHOPS (extension of HOPS Out to 500 ParSecs) census of nearby protostars. For every source detected in the Herschel/PACS bands, the eHOPS-Aquila catalog contains 1-850~$μ$m SEDs assembled from 2MASS, Spitzer, Herschel, WISE, and JCMT/SCUBA-2 data. Using a newly developed set of criteria, we… ▽ More We present a Spitzer/Herschel focused survey of the Aquila molecular clouds ($d \sim 436$~pc) as part of the eHOPS (extension of HOPS Out to 500 ParSecs) census of nearby protostars. For every source detected in the Herschel/PACS bands, the eHOPS-Aquila catalog contains 1-850~$μ$m SEDs assembled from 2MASS, Spitzer, Herschel, WISE, and JCMT/SCUBA-2 data. Using a newly developed set of criteria, we classify objects by their SEDs as protostars, pre-ms sequence stars with disks, and galaxies. A total of 172 protostars are found in Aquila, tightly concentrated in the molecular filaments that thread the clouds. Of these, 71 (42\%) are Class 0 protostars, 54 (31\%) are Class I protostars, 43 (25\%) are flat-spectrum protostars, and 4 (2\%) are Class II sources. Ten of the Class 0 protostars are young PACS Bright Red Sources similar to those discovered in Orion. We compare the SEDs to a grid of radiative transfer models to constrain the luminosities, envelope densities, and envelope masses of the protostars. A comparison of the eHOPS-Aquila to the HOPS protostars in Orion finds that the protostellar luminosity functions in the two star-forming regions are statistically indistinguishable, the bolometric temperatures/envelope masses of eHOPS-Aquila protostars are shifted to cooler temperatures/higher masses, and the eHOPS-Aquila protostars do not show the decline in luminosity with evolution found in Orion. We briefly discuss whether these differences are due to biases between the samples, diverging star formation histories, or the influence of environment on protostellar evolution. △ Less

Submitted 22 February, 2023; v1 submitted 24 September, 2022; originally announced September 2022.

Comments: Accepted in ApJS. NASA IPAC/IRSA Data Release DOI: 10.26131/IRSA553

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

Authors: Evan A. Rich, John D. Monnier, Alicia Aarnio, Anna S. E. Laws, Benjamin R. Setterholm, David J. Wilner, Nuria Calvet, Tim Harries, Chris Miller, Claire L. Davies, Fred C. Adams, Sean M. Andrews, Jaehan Bae, Catherine Espaillat, Alexandra Z. Greenbaum, Sasha Hinkley, Stefan Kraus, Lee Hartmann, Andrea Isella, Melissa McClure, Rebecca Oppenheimer, Laura M. Pérez, Zhaohuan Zhu

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

Submitted 12 June, 2022; originally announced June 2022.

Comments: 51 pages, 31 figures, 7 tables, accepted to AJ

Anisotropic Infall and Substructure formation in Embedded Disks

Authors: Aleksandra Kuznetsova, Jaehan Bae, Lee Hartmann, Mordecai-Mark Mac Low

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

Submitted 10 February, 2022; originally announced February 2022.

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

Gravity or turbulence V: Star forming regions undergoing violent relaxation

Authors: Andrea Bonilla-Barroso, J. Ballesteros-Paredes, Jesus Hernández, Luis Aguilar, Manuel Zamora, Lee W. Hartmann, Aleksandra Kuznetsova, Vianey Camacho, Verónica Lora

Abstract: Using numerical simulations of the formation and evolution of stellar clusters within molecular clouds, we show that the stars in clusters formed within collapsing molecular cloud clumps exhibit a constant velocity dispersion regardless of their mass, as expected in a violent relaxation processes. In contrast, clusters formed in turbulence-dominated environments exhibit an {\it inverse} mass segre… ▽ More Using numerical simulations of the formation and evolution of stellar clusters within molecular clouds, we show that the stars in clusters formed within collapsing molecular cloud clumps exhibit a constant velocity dispersion regardless of their mass, as expected in a violent relaxation processes. In contrast, clusters formed in turbulence-dominated environments exhibit an {\it inverse} mass segregated velocity dispersion, where massive stars exhibit larger velocity dispersions than low-mass cores, consistent with massive stars formed in massive clumps, which in turn, are formed through strong shocks. We furthermore use Gaia EDR3 to show that the stars in the Orion Nebula Cluster exhibit a constant velocity dispersion as a function of mass, suggesting that it has been formed by collapse within one free-fall time of its parental cloud, rather than in a turbulence-dominated environment during many free-fall times of a supported cloud. Additionally, we have addressed several of the criticisms of models of collapsing star forming regions: namely, the age spread of the ONC, the comparison of the ages of the stars to the free-fall time of the gas that formed it, the star formation efficiency, and the mass densities of clouds vs the mass densities of stellar clusters, showing that observational and numerical data are consistent with clusters forming in clouds undergoing a process of global, hierarchical and chaotic collapse, rather than been supported by turbulence. △ Less

Submitted 26 January, 2022; originally announced January 2022.

Comments: Accepted for publication in MNRAS

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

Authors: C. C. Espaillat, G. J. Herczeg, T. Thanathibodee, C. Pittman, N. Calvet, N. Arulanantham, K. France, Javier Serna, J. Hernandez, A. Kospal, F. M. Walter, A. Frasca, W. J. Fischer, C. M. Johns-Krull, P. C. Schneider, C. Robinson, Suzan Edwards, P. Abraham, Min Fang, J. Erkal, C. F. Manara, J. M. Alcala, E. Alecian, R. D. Alexander, J. Alonso-Santiago , et al. (37 additional authors not shown)

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

Submitted 17 January, 2022; originally announced January 2022.

Comments: accepted to ApJ

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

Authors: W. Zakri, S. T. Megeath, W. J. Fischer, Robert Gutermuth, Elise Furlan, Lee Hartmann, Nicole Karnath, Mayra Osorio, Emily Safron, Thomas Stanke, Amelia M. Stutz, John J. Tobin, Thomas S. Allen, Sam Federman, Nolan Habel, P. Manoj, Mayank Narang, Riwaj Pokhrel, Luisa Rebull, Patrick D. Sheehan, Dan M. Watson

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

Submitted 12 January, 2022; originally announced January 2022.

Comments: Accepted to ApJL

A gravitational and dynamical framework of star formation: The Orion Nebula

Authors: Marina Kounkel, Keivan G. Stassun, Kevin Covey, Lee Hartmann

Abstract: The Orion Nebula Cluster (ONC) is the most massive region of active star formation within a kpc of the Sun. Using Gaia DR3 parallaxes and proper motions, we examine the bulk motions of stars radially and tangentially relative to the cluster center. We find an age gradient with distance to the stars in the ONC, from 385 pc for the oldest stars, to 395 pc for the younger stars, indicating that the s… ▽ More The Orion Nebula Cluster (ONC) is the most massive region of active star formation within a kpc of the Sun. Using Gaia DR3 parallaxes and proper motions, we examine the bulk motions of stars radially and tangentially relative to the cluster center. We find an age gradient with distance to the stars in the ONC, from 385 pc for the oldest stars, to 395 pc for the younger stars, indicating that the star forming front is propagating into the cloud. We find an organized signature of rotation of the central cluster, but it is present only in stars younger than 2 Myr. We also observe a net infall of young stars into the center of the ONC's deep gravitational potential well. The infalling sources lie preferentially along the filament, on the other hand, outflowing sources are distributed spherically around the cluster, and they have larger velocity dispersion. We further propose a solution to a long-standing question of why the ONC shows a weak signature of expansion even though the cluster is likely bound: much of this expansion may be driven by unstable N-body interactions among stars, resulting in low-velocity ejections. Though analyzing signatures imprinted on stellar dynamics across different spatial scales, these observation shed new light on the signatures of formation and evolution of young clusters. △ Less

Submitted 15 September, 2022; v1 submitted 1 November, 2021; originally announced November 2021.

Comments: 14 pages, accepted to MNRAS

PENELLOPE: the ESO data legacy program to complement the Hubble UV Legacy Library of Young Stars (ULLYSES) I. Survey presentation and accretion properties of Orion OB1 and $σ$-Orionis

Authors: C. F. Manara, A. Frasca, L. Venuti, M. Siwak, G. J. Herczeg, N. Calvet, J. Hernandez, Ł. Tychoniec, M. Gangi, J. M. Alcalá, H. M. J. Boffin, B. Nisini, M. Robberto, C. Briceno, J. Campbell-White, A. Sicilia-Aguilar, P. McGinnis, D. Fedele, Á. Kóspál, P. Ábrahám, J. Alonso-Santiago, S. Antoniucci, N. Arulanantham, F. Bacciotti, A. Banzatti , et al. (47 additional authors not shown)

Abstract: The evolution of young stars and disks is driven by the interplay of several processes, notably accretion and ejection of material. Critical to correctly describe the conditions of planet formation, these processes are best probed spectroscopically. About five-hundred orbits of the Hubble Space Telescope (HST) are being devoted in 2020-2022 to the ULLYSES public survey of about 70 low-mass (M<2Msu… ▽ More The evolution of young stars and disks is driven by the interplay of several processes, notably accretion and ejection of material. Critical to correctly describe the conditions of planet formation, these processes are best probed spectroscopically. About five-hundred orbits of the Hubble Space Telescope (HST) are being devoted in 2020-2022 to the ULLYSES public survey of about 70 low-mass (M<2Msun) young (age<10 Myr) stars at UV wavelengths. Here we present the PENELLOPE Large Program that is being carried out at the ESO Very Large Telescope (VLT) to acquire, contemporaneous to HST, optical ESPRESSO/UVES high-resolution spectra to investigate the kinematics of the emitting gas, and UV-to-NIR X-Shooter medium-resolution flux-calibrated spectra to provide the fundamental parameters that HST data alone cannot provide, such as extinction and stellar properties. The data obtained by PENELLOPE have no proprietary time, and the fully reduced spectra are made available to the whole community. Here, we describe the data and the first scientific analysis of the accretion properties for the sample of thirteen targets located in the Orion OB1 association and in the sigma-Orionis cluster, observed in Nov-Dec 2020. We find that the accretion rates are in line with those observed previously in similarly young star-forming regions, with a variability on a timescale of days of <3. The comparison of the fits to the continuum excess emission obtained with a slab model on the X-Shooter spectra and the HST/STIS spectra shows a shortcoming in the X-Shooter estimates of <10%, well within the assumed uncertainty. Its origin can be either a wrong UV extinction curve or due to the simplicity of this modelling, and will be investigated in the course of the PENELLOPE program. The combined ULLYSES and PENELLOPE data will be key for a better understanding of the accretion/ejection mechanisms in young stars. △ Less

Submitted 6 April, 2021; v1 submitted 23 March, 2021; originally announced March 2021.

Comments: Accepted for publication on Astronomy & Astrophysics. 15 pages + appendix, language edited version

Journal ref: A&A 650, A196 (2021)

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

Authors: Michael L. Weber, Barbara Ercolano, Giovanni Picogna, Lee Hartmann, Peter J. Rodenkirch

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

Submitted 28 May, 2020; originally announced May 2020.

Comments: Accepted for publication in MNRAS. 20+3 pages, 14+2 figures

A Study of Millimeter Variability in FUor Objects

Authors: John Wendeborn, Catherine C. Espaillat, Enrique Macias, Orsolya Feher, A. Kospal, Lee Hartmann, Zhaohuan Zhu, Michael M. Dunham, Marina Kounkel

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

Submitted 20 May, 2020; originally announced May 2020.

Comments: 13 pages, 8 figures, to be published in The Astrophysical Journal

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

Authors: Stefan Kraus, Alexander Kreplin, Alison K. Young, Matthew R. Bate, John D. Monnier, Tim J. Harries, Henning Avenhaus, Jacques Kluska, Anna S. E. Laws, Evan A. Rich, Matthew Willson, Alicia N. Aarnio, Fred C. Adams, Sean M. Andrews, Narsireddy Anugu, Jaehan Bae, Theo ten Brummelaar, Nuria Calvet, Michel Curé, Claire L. Davies, Jacob Ennis, Catherine Espaillat, Tyler Gardner, Lee Hartmann, Sasha Hinkley , et al. (7 additional authors not shown)

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

Submitted 4 September, 2020; v1 submitted 2 April, 2020; originally announced April 2020.

Comments: 63 pages, 4+13 Figures, 6 Tables, published at Science

Angular Momenta, Magnetization, and Accretion of Protostellar Cores

Authors: Aleksandra Kuznetsova, Lee Hartmann, Fabian Heitsch

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

Submitted 30 March, 2020; originally announced March 2020.

Comments: 13 pages, 9 figures, accepted by ApJ

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

Authors: Anna S. E. Laws, Tim J. Harries, Benjamin R. Setterholm, John D. Monnier, Evan A. Rich, Alicia N. Aarnio, Fred C. Adams, Sean Andrews, Jaehan Bae, Nuria Calvet, Catherine Espaillat, Lee Hartmann, Sasha Hinkley, Andrea Isella, Stefan Kraus, David Wilner, Zhaohuan Zhu

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

Submitted 11 November, 2019; originally announced November 2019.

Principles Of Heliophysics: a textbook on the universal processes behind planetary habitability

Authors: Karel Schrijver, Fran Bagenal, Tim Bastian, Juerg Beer, Mario Bisi, Tom Bogdan, Steve Bougher, David Boteler, Dave Brain, Guy Brasseur, Don Brownlee, Paul Charbonneau, Ofer Cohen, Uli Christensen, Tom Crowley, Debrah Fischer, Terry Forbes, Tim Fuller-Rowell, Marina Galand, Joe Giacalone, George Gloeckler, Jack Gosling, Janet Green, Nick Gross, Steve Guetersloh , et al. (37 additional authors not shown)

Abstract: Heliophysics is the system science of the physical connections between the Sun and the solar system. As the physics of the local cosmos, it embraces space weather and planetary habitability. The wider view of comparative heliophysics forms a template for conditions in exoplanetary systems and provides a view over time of the aging Sun and its magnetic activity, of the heliosphere in different sett… ▽ More Heliophysics is the system science of the physical connections between the Sun and the solar system. As the physics of the local cosmos, it embraces space weather and planetary habitability. The wider view of comparative heliophysics forms a template for conditions in exoplanetary systems and provides a view over time of the aging Sun and its magnetic activity, of the heliosphere in different settings of the interstellar medium and subject to stellar impacts, of the space physics over evolving planetary dynamos, and of the long-term influence on planetary atmospheres by stellar radiation and wind. Based on a series of NASA-funded summer schools for early-career researchers, this textbook is intended for students in physical sciences in later years of their university training and for beginning graduate students in fields of solar, stellar, (exo-)planetary, and planetary-system sciences. The book emphasizes universal processes from a perspective that draws attention to what provides Earth (and similar (exo-)planets) with a relatively stable setting in which life as we know it could thrive. The text includes 200 "Activities" in the form of exercises, explorations, literature readings, "what if" challenges, and group discussion topics; many of the Activities provide additional information complementing the main text. Solutions and discussions are included in an Appendix for a selection of the exercises. △ Less

Submitted 27 March, 2024; v1 submitted 30 October, 2019; originally announced October 2019.

Comments: 429 pages, 121 figures, and 200 "activities" in the form of problems, exercises, literature readings, and "what if" challenges; V2.0: activities modified; added 'solutions' to selected activities; clarified text and corrected typos; added subject index; V2.1 added Activity hyperlinks and corrected typos (see 'Version history' in the manuscript)

Deuterated methanol toward NGC 7538-IRS1

Authors: Juan Ospina-Zamudio, Cécile Favre, Marina Kounkel, Li-Hong Xu, Justin Neill, Bertrand Lefloch, Alexandre Faure, Edwin Bergin, Davide Fedele, Lee Hartmann

Abstract: We investigate the deuteration of methanol towards the high-mass star forming region NGC 7538-IRS1. We have carried out a multi-transition study of CH$_3$OH, $^{13}$CH$_3$OH and of the deuterated fllavors, CH$_2$DOH and CH$_3$OD, between 1.0--1.4 mm with the IRAM-30~m antenna. In total, 34 $^{13}$CH$_3$OH, 13 CH$_2$DOH lines and 20 CH$_3$OD lines spanning a wide range of upper-state energies (E… ▽ More We investigate the deuteration of methanol towards the high-mass star forming region NGC 7538-IRS1. We have carried out a multi-transition study of CH$_3$OH, $^{13}$CH$_3$OH and of the deuterated fllavors, CH$_2$DOH and CH$_3$OD, between 1.0--1.4 mm with the IRAM-30~m antenna. In total, 34 $^{13}$CH$_3$OH, 13 CH$_2$DOH lines and 20 CH$_3$OD lines spanning a wide range of upper-state energies (E$_{up}$) were detected. From the detected transitions, we estimate that the measured D/H does not exceed 1$\%$, with a measured CH$_2$DOH/CH$_3$OH and CH$_3$OD/CH$_3$OH of about (32$\pm$8)$\times$10$^{-4}$ and (10$\pm$4)$\times$10$^{-4}$, respectively. This finding is consistent with the hypothesis of a short-time scale formation during the pre-stellar phase. We find a relative abundance ratio CH$_2$DOH/CH$_3$OD of 3.2 $\pm$ 1.5. This result is consistent with a statistical deuteration. We cannot exclude H/D exchanges between water and methanol if water deuteration is of the order 0.1$\%$, as suggested by recent Herschel observations. △ Less

Submitted 23 May, 2019; originally announced May 2019.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 627, A80 (2019)

The Origins of Protostellar Core Angular Momenta

Authors: Aleksandra Kuznetsova, Lee Hartmann, Fabian Heitsch

Abstract: We present the results of a suite of numerical simulations designed to explore the origin of the angular momenta of protostellar cores. Using the hydrodynamic grid code \emph{Athena} with a sink implementation, we follow the formation of protostellar cores and protostars (sinks) from the subvirial collapse of molecular clouds on larger scales to investigate the range and relative distribution of c… ▽ More We present the results of a suite of numerical simulations designed to explore the origin of the angular momenta of protostellar cores. Using the hydrodynamic grid code \emph{Athena} with a sink implementation, we follow the formation of protostellar cores and protostars (sinks) from the subvirial collapse of molecular clouds on larger scales to investigate the range and relative distribution of core properties. We find that the core angular momenta are relatively unaffected by large-scale rotation of the parent cloud; instead, we infer that angular momenta are mainly imparted by torques between neighboring mass concentrations and exhibit a log-normal distribution. Our current simulation results are limited to size scales $\sim 0.05$~pc ($\sim 10^4 \rm AU$), but serve as first steps toward the ultimate goal of providing initial conditions for higher-resolution studies of core collapse to form protoplanetary disks. △ Less

Submitted 22 March, 2019; originally announced March 2019.

Comments: 12 pages, 9 figures, accepted to ApJ

On the nature of the compact sources in IRAS 16293-2422 seen in at centimeter to sub-millimeter wavelengths

Authors: Antonio Hernández-Gómez, Laurent Loinard, Claire J. Chandler, Luis F. Rodríguez, Luis A. Zapata, David J. Wilner, Paul T. P. Ho, Emmanuel Caux, David Quénard, Sandrine Bottinelli, Crystal L. Brogan, Lee Hartmann, Karl M. Menten

Abstract: We present multi-epoch continuum observations of the Class 0 protostellar system IRAS 16293-2422 taken with the Very Large Array (VLA) at multiple wavelengths between 7 mm and 15 cm (41 GHz down to 2 GHz), as well as single-epoch Atacama Large Millimeter/submillimeter Array (ALMA) continuum observations covering the range from 0.4 to 1.3 mm (700 GHz down to 230 GHz). The new VLA observations confi… ▽ More We present multi-epoch continuum observations of the Class 0 protostellar system IRAS 16293-2422 taken with the Very Large Array (VLA) at multiple wavelengths between 7 mm and 15 cm (41 GHz down to 2 GHz), as well as single-epoch Atacama Large Millimeter/submillimeter Array (ALMA) continuum observations covering the range from 0.4 to 1.3 mm (700 GHz down to 230 GHz). The new VLA observations confirm that source A2 is a protostar driving episodic mass ejections, and reveal the complex relative motion between A2 and A1. The spectrum of component B can be described by a single power law ($S_ν\propto ν^{2.28}$) over the entire range from 3 to 700 GHz (10 cm down to 0.4 mm), suggesting that the emission is entirely dominated by dust even at $λ$ = 10 cm. Finally, the size of source B appears to increase with frequency up to 41 GHz, remaining roughly constant (at $0''.39$ $\equiv$ 55 AU) at higher frequencies. We interpret this as evidence that source B is a dusty structure of finite size that becomes increasingly optically thick at higher frequencies until, in the millimeter regime, the source becomes entirely optically thick. The lack of excess free-free emission at long wavelengths, combined with the absence of high-velocity molecular emission indicates that source B does not drive a powerful outflow, and might indicate that source B is at a particularly early stage of its evolution. △ Less

Submitted 6 March, 2019; originally announced March 2019.

Structure and Expansion Law of HII Regions in structured Molecular Clouds

Authors: Manuel Zamora-Avilés, Enrique Vázquez-Semadeni, Ricardo F. González, José Franco, Steven N. Shore, Lee W. Hartmann, Javier Ballesteros-Paredes, Robi Banerjee, Bastian Körtgen

Abstract: We present radiation-magnetohydrodynamic simulations aimed at studying evolutionary properties of H\,{\normalsize II} regions in turbulent, magnetised, and collapsing molecular clouds formed by converging flows in the warm neutral medium. We focus on the structure, dynamics and expansion laws of these regions. Once a massive star forms in our highly structured clouds, its ionising radiation eventu… ▽ More We present radiation-magnetohydrodynamic simulations aimed at studying evolutionary properties of H\,{\normalsize II} regions in turbulent, magnetised, and collapsing molecular clouds formed by converging flows in the warm neutral medium. We focus on the structure, dynamics and expansion laws of these regions. Once a massive star forms in our highly structured clouds, its ionising radiation eventually stops the accretion (through filaments) toward the massive star-forming regions. The new over-pressured H\,{\normalsize II} regions push away the dense gas, thus disrupting the more massive collapse centres. Also, because of the complex density structure in the cloud, the H\,{\normalsize II} regions expand in a hybrid manner: they virtually do not expand toward the densest regions (cores), while they expand according to the classical analytical result towards the rest of the cloud, and in an accelerated way, as a blister region, towards the diffuse medium. Thus, the ionised regions grow anisotropically, and the ionising stars generally appear off-centre of the regions. Finally, we find that the hypotheses assumed in standard H\,{\normalsize II}-region expansion models (fully embedded region, blister-type, or expansion in a density gradient) apply simultaneously in different parts of our simulated H\,{\normalsize II} regions, producing a net expansion law ($R \propto t^α$, with $α$ in the range of 0.93-1.47 and a mean value of $1.2 \pm 0.17$) that differs from any of those of the standard models. △ Less

Submitted 17 May, 2019; v1 submitted 4 March, 2019; originally announced March 2019.

Comments: Accepted for publication in MNRAS

Multiple spiral arms in the disk around intermediate-mass binary HD 34700A

Authors: John D. Monnier, Tim J. Harries, Jaehan Bae, Benjamin R. Setterholm, Anna Laws, Alicia Aarnio, Fred C. Adams, Sean Andrews, Nuria Calvet, Catherine Espaillat, Lee Hartmann, Stefan Kraus, Melissa McClure, Chris Miller, Rebecca Oppenheimer, David Wilner, Zhaohuan Zhu

Abstract: We present the first images of the transition disk around the close binary system HD 34700A in polarized scattered light using the Gemini Planet Imager instrument on Gemini South. The J and H band images reveal multiple spiral-arm structures outside a large (R = 0.49" = 175 au) cavity along with a bluish spiral structure inside the cavity. The cavity wall shows a strong discontinuity and we clearl… ▽ More We present the first images of the transition disk around the close binary system HD 34700A in polarized scattered light using the Gemini Planet Imager instrument on Gemini South. The J and H band images reveal multiple spiral-arm structures outside a large (R = 0.49" = 175 au) cavity along with a bluish spiral structure inside the cavity. The cavity wall shows a strong discontinuity and we clearly see significant non-azimuthal polarization Uphi consistent with multiple scattering within a disk at an inferred inclination ~42deg. Radiative transfer modeling along with a new Gaia distance suggest HD 37400A is a young (~5 Myr) system consisting of two intermediate-mass (~2Msun) stars surrounded by a transitional disk and not a solar-mass binary with a debris disk as previously classified. Conventional assumptions of the dust-to-gas ratio would rule out a gravitational instability origin to the spirals while hydrodynamical models using the known external companion or a hypothetical massive protoplanet in the cavity both have trouble reproducing the relatively large spiral arm pitch angles (~30deg) without fine tuning of gas temperature. We explore the possibility that material surrounding a massive protoplanet could explain the rim discontinuity after also considering effects of shadowing by an inner disk. Analysis of archival Hubble Space Telescope data suggests the disk is rotating counterclockwise as expected from the spiral arm structure and revealed a new low-mass companion at 6.45" separation. We include an appendix which sets out clear definitions of Q, U, Qphi, Uphi, correcting some confusion and errors in the literature. △ Less

Submitted 8 January, 2019; originally announced January 2019.

Comments: Accepted by the Astrophysical Journal on Dec 15, 2018

Gaia-DR2 confirms VLBA parallaxes in Ophiuchus, Serpens and Aquila

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

Abstract: We present Gaia-DR2 astrometry of a sample of YSO candidates in Ophiuchus, Serpens Main and Serpens South/W40 in the Aquila Rift, which had been mainly identified by their infrared excess with Spitzer. We compare the Gaia-DR2 parallaxes against published and new parallaxes obtained from our Very Long Baseline Array (VLBA) program GOBELINS. We obtain consistent results between Gaia and the VLBA for… ▽ More We present Gaia-DR2 astrometry of a sample of YSO candidates in Ophiuchus, Serpens Main and Serpens South/W40 in the Aquila Rift, which had been mainly identified by their infrared excess with Spitzer. We compare the Gaia-DR2 parallaxes against published and new parallaxes obtained from our Very Long Baseline Array (VLBA) program GOBELINS. We obtain consistent results between Gaia and the VLBA for the mean parallaxes in each of the regions analyzed here. We see small offsets, when comparing mean values, of a few tens of micro-arcseconds in the parallaxes, which are either introduced by the Gaia zero-point error or due to a selection effect by Gaia toward the brightest, less obscured stars. Gaia-DR2 data alone conclusively places Serpens Main and Serpens South at the same distance, as we first inferred from VLBA data alone in a previous publication. Thus, Serpens Main, Serpens South and W40 are all part of the same complex of molecular clouds, located at a mean distance of 436+/-9 pc. In Ophiuchus, both Gaia and VLBA suggest a small parallax gradient across the cloud, and the distance changes from 144.2+/-1.3 pc to 138.4+/-2.6 pc when going from L1689 to L1688. △ Less

Submitted 6 December, 2018; originally announced December 2018.

Comments: Accepted for publication in ApJL

Disc wind models for FU Ori objects

Authors: Kelly Milliner, James H. Matthews, Knox S. Long, Lee Hartmann

Abstract: We present disc wind models aimed at reproducing the main features of the strong Na I resonance line P-Cygni profiles in the rapidly-accreting pre-main sequence FU Ori objects. We conducted Monte Carlo radiative transfer simulations for a standard magnetocentrifugally driven wind (MHD) model and our own "Genwind" models, which allows for a more flexible wind parameterisation. We find that the fidu… ▽ More We present disc wind models aimed at reproducing the main features of the strong Na I resonance line P-Cygni profiles in the rapidly-accreting pre-main sequence FU Ori objects. We conducted Monte Carlo radiative transfer simulations for a standard magnetocentrifugally driven wind (MHD) model and our own "Genwind" models, which allows for a more flexible wind parameterisation. We find that the fiducial MHD wind and similar Genwind models, which have flows emerging outward from the inner disc edge, and thus have polar cavities with no absorbing gas, cannot reproduce the deep, wide Na I absorption lines in FU Ori objects viewed at low inclination. We find that it is necessary to include an "inner wind" to fill this polar cavity to reproduce observations. In addition, our models assuming pure scattering source functions in the Sobolev approximation at intermediate viewing angles ($30^{\circ} \lesssim i \lesssim 60^{\circ}$) do not yield sufficiently deep line profiles. Assuming complete absorption yields better agreement with observations, but simple estimates strongly suggest that pure scattering should be a much better approximation. The discrepancy may indicate that the Sobolev approximation is not applicable, possibly due to turbulence or non-monotonic velocity fields; there is some observational evidence for the latter. Our results provide guidance for future attempts to constrain FU Ori wind properties using full MHD wind simulations, by pointing to the importance of the boundary conditions necessary to give rise to an inner wind, and by suggesting that the winds must be turbulent to produce sufficiently deep line profiles. △ Less

Submitted 29 November, 2018; originally announced November 2018.

Comments: 12 pages, 17 figures, accepted for publication in MNRAS

Science with an ngVLA: Resolving the Radio Complexity of EXor and FUor-type Systems with the ngVLA

Authors: Jacob Aaron White, Marc Audard, Péter Ábrahám, Lucas Cieza, Fernando Cruz-Sáenz de Miera, Michael M. Dunham, Joel D. Green, Manuel Güdel, Nicolas Grosso, Antonio Hales, Lee Hartmann, Kundan Kadam, Joel H. Kastner, Ágnes Kóspál, Sebastian Perez, Andreas Postel, Dary Ruiz-Rodriguez, Christian Rab, Eduard I. Vorobyov, Zhaohuan Zhu

Abstract: Episodic accretion may be a common occurrence in the evolution of young pre-main sequence stars and has important implications for our understanding of star and planet formation. Many fundamental aspects of what drives the accretion physics, however, are still unknown. The ngVLA will be a key tool in understanding the nature of these events. The high spatial resolution, broad spectral coverage, an… ▽ More Episodic accretion may be a common occurrence in the evolution of young pre-main sequence stars and has important implications for our understanding of star and planet formation. Many fundamental aspects of what drives the accretion physics, however, are still unknown. The ngVLA will be a key tool in understanding the nature of these events. The high spatial resolution, broad spectral coverage, and unprecedented sensitivity will allow for the detailed analysis of outburst systems. The proposed frequency range of the ngVLA allows for observations of the gas, dust, and non-thermal emission from the star and disk. △ Less

Submitted 15 October, 2018; originally announced October 2018.

Comments: 8 pages, 1 figure, To be published in the ASP Monograph Series, "Science with a Next-Generation VLA", ed. E. J. Murphy (ASP, San Francisco, CA)

The Role of Gravity in Producing Power-Law Mass Functions

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

Abstract: Numerical simulations of star formation have found that a power-law mass function can develop at high masses. In a previous paper, we employed isothermal simulations which created large numbers of sinks over a large range in masses to show that the power law exponent of the mass function, $dN/d\log M \propto M^Γ$, asymptotically and accurately approaches $Γ= -1.$ Simple analytic models show that s… ▽ More Numerical simulations of star formation have found that a power-law mass function can develop at high masses. In a previous paper, we employed isothermal simulations which created large numbers of sinks over a large range in masses to show that the power law exponent of the mass function, $dN/d\log M \propto M^Γ$, asymptotically and accurately approaches $Γ= -1.$ Simple analytic models show that such a power law can develop if the mass accretion rate $\dot{M} \propto M^2$, as in Bondi-Hoyle accretion; however, the sink mass accretion rates in the simulations show significant departures from this relation. In this paper we show that the expected accretion rate dependence is more closely realized provided the gravitating mass is taken to be the sum of the sink mass and the mass in the near environment. This reconciles the observed mass functions with the accretion rate dependencies, and demonstrates that power-law upper mass functions are essentially the result of gravitational focusing, a mechanism present in, for example, the competitive accretion model. △ Less

Submitted 8 October, 2018; originally announced October 2018.

Comments: 11 pages, 10 figures, accepted by ApJ

On estimating angular momenta of infalling protostellar cores from observations

Authors: Shangjia Zhang, Lee Hartmann, Manuel Zamora, Aleksandra Kuznetsova

Abstract: We use numerical simulations of molecular cloud formation in the colliding flow scenario to investigate the reliability of observational estimates of the angular momenta of early-state, low-mass protostellar cores. We show that, with suitable corrections for projection factors, molecular line observations of velocity gradients in NH$_3$ can be used to provide reasonable estimates of core angular m… ▽ More We use numerical simulations of molecular cloud formation in the colliding flow scenario to investigate the reliability of observational estimates of the angular momenta of early-state, low-mass protostellar cores. We show that, with suitable corrections for projection factors, molecular line observations of velocity gradients in NH$_3$ can be used to provide reasonable estimates of core angular momenta within a factor of two to three, with a few large underestimates due to unfavorable viewing angles. Our results differ from previous investigations which suggested that observations might overestimate true angular momenta by as much as an order of magnitude; the difference is probably due to the much smoother velocity field on small scales in our simulations, which result from allowing turbulence to decay and gravitational infall to dominate. The results emphasize the importance of understanding the nature of 'turbulent' velocities, with implications for the formation of protostellar disks during core collapse. △ Less

Submitted 14 August, 2018; originally announced August 2018.

Comments: 10 pages, 9 figures, accepted to MNRAS

The Gould's Belt Distances Survey (GOBELINS). V. Distances and Kinematics of the Perseus molecular cloud

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

Abstract: We derive the distance and structure of the Perseus molecular cloud by combining trigonometric parallaxes from Very Long Baseline Array (VLBA) observations, taken as part of the GOBELINS survey, and Gaia Data Release 2. Based on our VLBA astrometry, we obtain a distance of 321+/-10 pc for IC 348. This is fully consistent with the mean distance of 320+/-26 measured by Gaia. The VLBA observations to… ▽ More We derive the distance and structure of the Perseus molecular cloud by combining trigonometric parallaxes from Very Long Baseline Array (VLBA) observations, taken as part of the GOBELINS survey, and Gaia Data Release 2. Based on our VLBA astrometry, we obtain a distance of 321+/-10 pc for IC 348. This is fully consistent with the mean distance of 320+/-26 measured by Gaia. The VLBA observations toward NGC 1333 are insufficient to claim a successful distance measurement to this cluster. Gaia parallaxes, on the other hand, yield a mean distance of 293+/-22 pc. Hence, the distance along the line of sight between the eastern and western edges of the cloud is ~30 pc, which is significantly smaller than previously inferred. We use Gaia proper motions and published radial velocities to derive the spatial velocities of a selected sample of stars. The average velocity vectors with respect to the LSR are (u,v,w) = (-6.1+/-1.6, 6.8+/-1.1, -0.9+/-1.2) and (-6.4+/-1.0, 2.1+/-1.4, -2.4+/-1.0) km/s for IC 348 and NGC 1333, respectively. Finally, our analysis of the kinematics of the stars has shown that there is no clear evidence of expansion, contraction, or rotational motions within the clusters. △ Less

Submitted 10 August, 2018; originally announced August 2018.

Comments: Accepted for publication in ApJ

The Gould's Belt Distances Survey (GOBELINS). IV. Distance, Depth and Kinematics of the Taurus Star-Forming Region

Authors: Phillip A. B. Galli, Laurent Loinard, Gisela N. Ortiz-Leon, Marina Kounkel, Sergio A. Dzib, Amy J. Mioduszewski, Luis F. Rodriguez, Lee Hartmann, Ramachrisna Teixeira, Rosa M. Torres, Juana L. Rivera, Andrew F. Boden, Neal J. Evans II, Cesar Briceno, John J. Tobin, Mark Heyer

Abstract: We present new trigonometric parallaxes and proper motions of young stellar objects in the Taurus molecular cloud complex from observations collected with the Very Long Baseline Array as part of the Gould's Belt Distances Survey (GOBELINS). We detected 26 young stellar objects and derived trigonometric parallaxes for 18 stars with an accuracy of 0.3$\%$ to a few percent. We modeled the orbits of s… ▽ More We present new trigonometric parallaxes and proper motions of young stellar objects in the Taurus molecular cloud complex from observations collected with the Very Long Baseline Array as part of the Gould's Belt Distances Survey (GOBELINS). We detected 26 young stellar objects and derived trigonometric parallaxes for 18 stars with an accuracy of 0.3$\%$ to a few percent. We modeled the orbits of six binaries and determined the dynamical masses of the individual components in four of these systems (V1023 Tau, T Tau S, V807 Tau and V1000 Tau). Our results are consistent with the first trigonometric parallaxes delivered by the Gaia satellite and reveal the existence of significant depth effects. We find that the central portion of the dark cloud Lynds 1495 is located at $d=129.5\pm 0.3$ pc while the B 216 clump in the filamentary structure connected to it is at $d=158.1\pm 1.2$ pc. The closest and remotest stars in our sample are located at $d=126.6\pm 1.7$ pc and $d=162.7\pm 0.8$ pc yielding a distance difference of about 36 pc. We also provide a new distance estimate for HL Tau that was recently imaged. Finally, we compute the spatial velocity of the stars with published radial velocity and investigate the kinematic properties of the various clouds and gas structures in this region. △ Less

Submitted 23 May, 2018; originally announced May 2018.

Comments: ApJ (published online)

The CIDA Variability Survey of Orion OB1 II: demographics of the young, low-mass stellar populations

Authors: Cesar Briceno, Nuria Calvet, Jesus Hernandez, Anna K. Vivas, Cecilia Mateu, Juan Jose Downes, Jaqueline Loerincs, Alice Perez-Blanco, Perry Berlind, Catherine Espaillat, Lori Allen, Lee Hartmann, Mario Mateo, John Bailey III

Abstract: We present results of our large scale, optical, multi-epoch photometric survey of ~180 square degrees across the Orion OB1 association, complemented with extensive follow up spectroscopy. We map and characterize in an uniform way the off-cloud, low-mass, pre-main sequence populations. We report 2064, mostly K and M-type, confirmed T Tauri members. Most (59%) are located in the OB1a subassociation,… ▽ More We present results of our large scale, optical, multi-epoch photometric survey of ~180 square degrees across the Orion OB1 association, complemented with extensive follow up spectroscopy. We map and characterize in an uniform way the off-cloud, low-mass, pre-main sequence populations. We report 2064, mostly K and M-type, confirmed T Tauri members. Most (59%) are located in the OB1a subassociation, 27% in OB1b, and 14% within the confines of the A and B molecular clouds. There is significant structure in the spatial distribution of the young stars. We characterize two new clusterings of T Tauri stars, HD 35762 and HR 1833 in the OB1a subassociation, and two stellar overdensities in OB1b. There is indication of two populations of young stars in the OB1b region, located at two different distances, which may be due to the OB1a subassociation overlapping on front of OB1b. The various groups and regions can be ordered in an age sequence that agrees with the long standing picture of star formation starting in Orion OB1a some 10-15 Myr ago. We define a new type of T Tauri star, the C/W class, objects we propose may be nearing the end of their accretion phase. We detect the observational signature of Li depletion in young K and M stars with a timescale of 8.5 Myr. The decline of the accretion fraction from ~2 - 10 Myr, implies an accretion e-folding timescale of 2.1 Myr. Finally, the median amplitude of the V-band variability shows the decline of stellar activity, from accreting Classical T Tauri stars to the least active field dwarfs. △ Less

Submitted 2 May, 2018; originally announced May 2018.

Comments: 41 pages, 30 figures, submitted to the Astronomical Journal on May 02, 2018. Abridged abstract

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

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