Showing 1–50 of 170 results for author: Tobin, J J

Constraining the Stellar Masses and Origin of the Protostellar VLA 1623 System

Authors: Sarah I Sadavoy, Patrick Sheehan, John J. Tobin, Nadia M. Murillo, Richard Teague, Ian Stephens, Thomas Henning, Philip C. Myers, Edwin A. Bergin

Abstract: We present ALMA Band 7 molecular line observations of the protostars within the VLA 1623 system. We map C$^{17}$O (3 - 2) in the circumbinary disk around VLA 1623A and the outflow cavity walls of the collimated outflow. We further detect red-shifted and blue-shifted velocity gradients in the circumstellar disks around VLA 1623B and VLA 1623W that are consistent with Keplerian rotation. We use the… ▽ More We present ALMA Band 7 molecular line observations of the protostars within the VLA 1623 system. We map C$^{17}$O (3 - 2) in the circumbinary disk around VLA 1623A and the outflow cavity walls of the collimated outflow. We further detect red-shifted and blue-shifted velocity gradients in the circumstellar disks around VLA 1623B and VLA 1623W that are consistent with Keplerian rotation. We use the radiative transfer modeling code, pdspy, and simple flared disk models to measure stellar masses of $0.27 \pm 0.03$ M$_\odot$, $1.9^{+0.3}_{-0.2}$ M$_\odot$, and $0.64 \pm 0.06$ M$_\odot$ for the VLA 1623A binary, VLA 1623B, and VLA 1623W, respectively. These results represent the strongest constraints on stellar mass for both VLA 1623B and VLA 1623W, and the first measurement of mass for all stellar components using the same tracer and methodology. We use these masses to discuss the relationship between the young stellar objects (YSOs) in the VLA 1623 system. We find that VLA 1623W is unlikely to be an ejected YSO, as has been previously proposed. While we cannot rule out that VLA 1623W is a unrelated YSO, we propose that it is a true companion star to the VLA 1623A/B system and that the these stars formed in situ through turbulent fragmentation and have had only some dynamical interactions since their inception. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Comments: Accepted to A&A; 16 pages, 12 figures

Identification of hot gas around low-mass protostars

Authors: Merel L. R. van 't Hoff, Edwin A. Bergin, Penelope Riley, Sanil Mittal, Jes K. Jørgensen, John J. Tobin

Abstract: The low carbon content of Earth and primitive meteorites compared to the Sun and interstellar grains suggests that carbon-rich grains were destroyed in the inner few astronomical units of the young solar system. A promising mechanism to selectively destroy carbonaceous grains is thermal sublimation within the soot line at $\gtrsim$ 300 K. To address whether such hot conditions are common amongst l… ▽ More The low carbon content of Earth and primitive meteorites compared to the Sun and interstellar grains suggests that carbon-rich grains were destroyed in the inner few astronomical units of the young solar system. A promising mechanism to selectively destroy carbonaceous grains is thermal sublimation within the soot line at $\gtrsim$ 300 K. To address whether such hot conditions are common amongst low-mass protostars, we observe CH$_3$CN transitions at 1, 2 and 3 mm with the NOrthern Extended Millimeter Array (NOEMA) toward seven low-mass and one intermediate-mass protostar ($L_{\rm{bol}} \sim2-300 L_\odot$), as CH$_3$CN is an excellent temperature tracer. We find $>$ 300 K gas toward all sources, indicating that hot gas may be prevalent. Moreover, the excitation temperature for CH$_3$OH obtained with the same observations is always lower ($\sim$135-250 K), suggesting that CH$_3$CN and CH$_3$OH have a different spatial distribution. A comparison of the column densities at 1 and 3 mm shows a stronger increase at 3 mm for CH$_3$CN than for CH$_3$OH. Since the dust opacity is lower at longer wavelengths, this indicates that CH$_3$CN is enhanced in the hot gas compared to CH$_3$OH. If this CH$_3$CN enhancement is the result of carbon-grain sublimation, these results suggests that Earth's initial formation conditions may not be rare. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Comments: 19 pages, 9 figures, 4 tables (plus 13 pages appendix with 10 figures, 4 tables). Accepted for publication in ApJ

Early Planet Formation in Embedded Disks (eDisk) XV: Influence of Magnetic Field Morphology in Dense Cores on Sizes of Protostellar Disks

Authors: Hsi-Wei Yen, Jonathan P. Williams, Jinshi Sai, Patrick M. Koch, Ilseung Han, Jes K. Jørgensen, Woojin Kwon, Chang Won Lee, Zhi-Yun Li, Leslie W. Looney, Mayank Narang, Nagayoshi Ohashi, Shigehisa Takakuwa, John J. Tobin, Itziar de Gregorio-Monsalvo, Shih-Ping Lai, Jeong-Eun Lee, Kengo Tomida

Abstract: The magnetic field of a molecular cloud core may play a role in the formation of circumstellar disks in the core. We present magnetic field morphologies in protostellar cores of 16 targets in the Atacama Large Millimeter/submillimeter Array large program "Early Planet Formation in Embedded Disks (eDisk)", which resolved their disks with 7 au resolutions. The 0.1-pc scale magnetic field morphologie… ▽ More The magnetic field of a molecular cloud core may play a role in the formation of circumstellar disks in the core. We present magnetic field morphologies in protostellar cores of 16 targets in the Atacama Large Millimeter/submillimeter Array large program "Early Planet Formation in Embedded Disks (eDisk)", which resolved their disks with 7 au resolutions. The 0.1-pc scale magnetic field morphologies were inferred from the James Clerk Maxwell Telescope (JCMT) POL-2 observations. The mean orientations and angular dispersions of the magnetic fields in the dense cores are measured and compared with the radii of the 1.3 mm continuum disks and the dynamically determined protostellar masses from the eDisk program. We observe a significant correlation between the disk radii and the stellar masses. We do not find any statistically significant dependence of the disk radii on the projected misalignment angles between the rotational axes of the disks and the magnetic fields in the dense cores, nor on the angular dispersions of the magnetic fields within these cores. However, when considering the projection effect, we cannot rule out a positive correlation between disk radii and misalignment angles in three-dimensional space. Our results suggest that the morphologies of magnetic fields in dense cores do not play a dominant role in the disk formation process. Instead, the sizes of protostellar disks may be more strongly affected by the amount of mass that has been accreted onto star+disk systems, and possibly other parameters, for example, magnetic field strength, core rotation, and magnetic diffusivity. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: 14 pages, 4 figures, accepted by ApJ

JWST detections of amorphous and crystalline HDO ice toward massive protostars

Authors: Katerina Slavicinska, Ewine F. van Dishoeck, Łukasz Tychoniec, Pooneh Nazari, Adam E. Rubinstein, Robert Gutermuth, Himanshu Tyagi, Yuan Chen, Nashanty G. C. Brunken, Will R. M. Rocha, P. Manoj, Mayank Narang, S. Thomas Megeath, Yao-Lun Yang, Leslie W. Looney, John J. Tobin, Henrik Beuther, Tyler L. Bourke, Harold Linnartz, Samuel Federman, Dan M. Watson, Hendrik Linz

Abstract: This work aims to utilize the increased sensitivity and resolution of the JWST to quantify the HDO/H$_{2}$O ratio in ices toward young stellar objects (YSOs) and to determine if the HDO/H$_{2}$O ratios measured in the gas phase toward massive YSOs (MYSOs) are representative of the ratios in their ice envelopes. Two protostars observed in the Investigating Protostellar Accretion (IPA) program using… ▽ More This work aims to utilize the increased sensitivity and resolution of the JWST to quantify the HDO/H$_{2}$O ratio in ices toward young stellar objects (YSOs) and to determine if the HDO/H$_{2}$O ratios measured in the gas phase toward massive YSOs (MYSOs) are representative of the ratios in their ice envelopes. Two protostars observed in the Investigating Protostellar Accretion (IPA) program using JWST NIRSpec were analyzed: HOPS 370, an intermediate-mass YSO (IMYSO), and IRAS 20126+4104, a MYSO. The HDO ice toward these sources was detected above the 3$σ$ level and quantified via its 4.1 $μ$m band. The contributions from the CH$_{3}$OH combination modes to the observed optical depth in this spectral region were constrained via the CH$_{3}$OH 3.53 $μ$m band to ensure that the integrated optical depth of the HDO feature was not overestimated. H$_{2}$O ice was quantified via its 3 $μ$m band. From these fits, ice HDO/H$_{2}$O abundance ratios of 4.6$\pm$1.8$\times$10$^{-3}$ and 2.6$\pm$1.2$\times$10$^{-3}$ are obtained for HOPS 370 and IRAS 20126+4104, respectively. The simultaneous detections of both crystalline HDO and crystalline H$_{2}$O corroborate the assignment of the observed feature at 4.1 $μ$m to HDO ice. The ice HDO/H$_{2}$O ratios are similar to the highest reported gas HDO/H$_{2}$O ratios measured toward MYSOs as well as the hot inner regions of isolated low-mass protostars, suggesting that at least some of the gas HDO/H$_{2}$O ratios measured toward massive hot cores are representative of the HDO/H$_{2}$O ratios in ices. The need for an H$_{2}$O-rich CH$_{3}$OH component in the CH$_{3}$OH ice analysis supports recent experimental and observational results that indicate that some CH$_{3}$OH ice may form prior to the CO freeze-out stage in H$_{2}$O-rich ice layers. △ Less

Submitted 23 April, 2024; originally announced April 2024.

Comments: Accepted for publication in A&A. 23 pages, 17 figures, 10 tables

JWST/MIRI detection of suprathermal OH rotational emissions: probing the dissociation of the water by Lyman alpha photons near the protostar HOPS 370

Authors: David A. Neufeld, P. Manoj, Himanshu Tyagi, Mayank Narang, Dan M. Watson, S. Thomas Megeath, Ewine F. Van Dishoeck, Robert A. Gutermuth, Thomas Stanke, Yao-Lun Yang, Adam E. Rubinstein, Guillem Anglada, Henrik Beuther, Alessio Caratti o Garatti, Neal J. Evans II, Samuel Federman, William J. Fischer, Joel Green, Pamela Klaassen, Leslie W. Looney, Mayra Osorio, Pooneh Nazari, John J. Tobin, Lukasz Tychoniec, Scott Wolk

Abstract: Using the MIRI/MRS spectrometer on JWST, we have detected pure rotational, suprathermal OH emissions from the vicinity of the intermediate-mass protostar HOPS 370 (OMC2/FIR3). These emissions are observed from shocked knots in a jet/outflow, and originate in states of rotational quantum number as high as 46 that possess excitation energies as large as $E_U/k = 4.65 \times 10^4$ K. The relative str… ▽ More Using the MIRI/MRS spectrometer on JWST, we have detected pure rotational, suprathermal OH emissions from the vicinity of the intermediate-mass protostar HOPS 370 (OMC2/FIR3). These emissions are observed from shocked knots in a jet/outflow, and originate in states of rotational quantum number as high as 46 that possess excitation energies as large as $E_U/k = 4.65 \times 10^4$ K. The relative strengths of the observed OH lines provide a powerful diagnostic of the ultraviolet radiation field in a heavily-extinguished region ($A_V \sim 10 - 20$) where direct UV observations are impossible. To high precision, the OH line strengths are consistent with a picture in which the suprathermal OH states are populated following the photodissociation of water in its $\tilde B - X$ band by ultraviolet radiation produced by fast ($\sim 80\,\rm km\,s^{-1}$) shocks along the jet. The observed dominance of emission from symmetric ($A^\prime$) OH states over that from antisymmetric ($A^{\prime\prime}$) states provides a distinctive signature of this particular population mechanism. Moreover, the variation of intensity with rotational quantum number suggests specifically that Ly$α$ radiation is responsible for the photodissociation of water, an alternative model with photodissociation by a 10$^4$ K blackbody being disfavored at a high level of significance. Using measurements of the Br$α$ flux to estimate the Ly$α$ production rate, we find that $\sim 4\%$ of the Ly$α$ photons are absorbed by water. Combined with direct measurements of water emissions in the $ν_2 = 1 -0$ band, the OH observations promise to provide key constraints on future models for the diffusion of Ly$α$ photons in the vicinity of a shock front. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: 30 pages, 7 figures. Accepted for publication in ApJ Letters

An Observational View of Structure in Protostellar Systems

Authors: John J. Tobin, Patrick D. Sheehan

Abstract: The envelopes and disks that surround protostars reflect the initial conditions of star and planet formation and govern the assembly of stellar masses. Characterizing these structures requires observations that span the near-infrared to centimeter wavelengths. Consequently, the past two decades have seen progress driven by numerous advances in observational facilities across this spectrum, includi… ▽ More The envelopes and disks that surround protostars reflect the initial conditions of star and planet formation and govern the assembly of stellar masses. Characterizing these structures requires observations that span the near-infrared to centimeter wavelengths. Consequently, the past two decades have seen progress driven by numerous advances in observational facilities across this spectrum, including the \textit{Spitzer Space Telescope}, \textit{Herschel Space Observatory}, the Atacama Large Millimeter/submillimeter Array, and a host of other ground-based interferometers and single-dish radio telescopes. Nearly all protostars appear to have well-formed circumstellar disks that are likely to be rotationally-supported; the ability to detect a disk around a protostar is more a question of spatial resolution than whether or not a disk is present. The disks around protostars have inherently higher millimeter/submillimeter luminosities as compared to disks around more-evolved pre-main sequence stars, though there may be systematic variations between star forming regions. The envelopes around protostars are inherently asymmetric and streamers emphasize that mass flow through the envelopes to the disks may not be homogeneous. The current mass distribution of protostars may be impacted by selection bias given that it is skewed toward solar-mass protostars, inconsistent with the stellar initial mass function. △ Less

Submitted 22 March, 2024; originally announced March 2024.

Comments: To appear in Annual Reviews in Astronomy and Astrophysics; 42 pages, 12 Figures, 1 Table

Early Planet Formation in Embedded Disks (eDisk) XIII: Aligned Disks with Non-Settled Dust Around the Newly Resolved Class 0 Protobinary R CrA IRAS 32

Authors: Frankie J. Encalada, Leslie W. Looney, Shigehisa Takakuwa, John J. Tobin, Nagayoshi Ohashi, Jes K. Jørgensen, Zhi-Yun Li, Yuri Aikawa, Yusuke Aso, Patrick M. Koch, Woojin Kwon, Shih-Ping Lai, Chang Won Lee, Zhe-Yu Daniel Lin, Alejandro Santamarıa-Miranda, Itziar de Gregorio-Monsalvo, Nguyen Thi Phuong, Adele Plunkett, Jinshi Sai, Rajeeb Sharma, Hsi-Wei Yen, Ilseung Han

Abstract: Young protostellar binary systems, with expected ages less than $\sim$10$^5$ years, are little modified since birth, providing key clues to binary formation and evolution. We present a first look at the young, Class 0 binary protostellar system R CrA IRAS 32 from the Early Planet Formation in Embedded Disks (eDisk) ALMA large program, which observed the system in the 1.3 mm continuum emission,… ▽ More Young protostellar binary systems, with expected ages less than $\sim$10$^5$ years, are little modified since birth, providing key clues to binary formation and evolution. We present a first look at the young, Class 0 binary protostellar system R CrA IRAS 32 from the Early Planet Formation in Embedded Disks (eDisk) ALMA large program, which observed the system in the 1.3 mm continuum emission, $^{12}$CO (2-1), $^{13}$CO (2-1), C$^{18}$O (2-1), SO (6$_5$-5$_4$), and nine other molecular lines that trace disk, envelope, shocks, and outflows. With a continuum resolution of $\sim$0.03$^{\prime\prime}$ ($\sim$5 au, at a distance of 150 pc), we characterize the newly discovered binary system with a separation of 207 au, their circumstellar disks, and a circumbinary disk-like structure. The circumstellar disk radii are 26.9$\pm$0.3 and 22.8$\pm$0.3 au for sources A and B, respectively, and their circumstellar disk dust masses are estimated as 22.5$\pm$1.1 and 12.4$\pm$0.6 M$_{\Earth}$. The circumstellar disks and the circumbinary structure have well aligned position angles and inclinations, indicating formation in a smooth, ordered process such as disk fragmentation. In addition, the circumstellar disks have a near/far-side asymmetry in the continuum emission suggesting that the dust has yet to settle into a thin layer near the midplane. Spectral analysis of CO isotopologues reveals outflows that originate from both of the sources and possibly from the circumbinary disk-like structure. Furthermore, we detect Keplerian rotation in the $^{13}$CO isotopologues toward both circumstellar disks and likely Keplerian rotation in the circumbinary structure; the latter suggests that it is probably a circumbinary disk. △ Less

Submitted 21 March, 2024; originally announced March 2024.

Comments: 33 pages, 29 figures, 6 tables

On the magnetic field properties of protostellar envelopes in Orion

Authors: Bo Huang, Josep M. Girart, Ian W. Stephens, Manuel Fernandez-Lopez, Hector G. Arce, John M. Carpenter, Paulo Cortes, Erin G. Cox, Rachel Friesen, Valentin J. M. Le Gouellec, Charles L. H. Hull, Nicole Karnath, Woojin Kwon, Zhi-Yun Li, Leslie W. Looney, Tom Megeath, Philip C. Myers, Nadia M. Murillo, Jaime E. Pineda, Sarah Sadavoy, Alvaro Sanchez-Monge, Patricio Sanhueza, John J. Tobin, Qizhou Zhang, James M. Jackson , et al. (1 additional authors not shown)

Abstract: We present 870 um polarimetric observations toward 61 protostars in the Orion molecular clouds, with ~400 au (1") resolution using the Atacama Large Millimeter/submillimeter Array. We successfully detect dust polarization and outflow emission in 56 protostars, in 16 of them the polarization is likely produced by self-scattering. Self-scattering signatures are seen in several Class 0 sources, sugge… ▽ More We present 870 um polarimetric observations toward 61 protostars in the Orion molecular clouds, with ~400 au (1") resolution using the Atacama Large Millimeter/submillimeter Array. We successfully detect dust polarization and outflow emission in 56 protostars, in 16 of them the polarization is likely produced by self-scattering. Self-scattering signatures are seen in several Class 0 sources, suggesting that grain growth appears to be significant in disks at earlier protostellar phases. For the rest of the protostars, the dust polarization traces the magnetic field, whose morphology can be approximately classified into three categories: standard-hourglass, rotated-hourglass (with its axis perpendicular to outflow), and spiral-like morphology. 40.0% (+-3.0%) of the protostars exhibit a mean magnetic field direction approximately perpendicular to the outflow on several 100--1000 au scales. However, in the remaining sample, this relative orientation appears to be random, probably due to the complex set of morphologies observed. Furthermore, we classify the protostars into three types based on the C17O (3--2) velocity envelope's gradient: perpendicular to outflow, non-perpendicular to outflow, and unresolved gradient (<1.0~km/s/arcsec). In protostars with a velocity gradient perpendicular to outflow, the magnetic field lines are preferentially perpendicular to outflow, most of them exhibit a rotated hourglass morphology, suggesting that the magnetic field has been overwhelmed by gravity and angular momentum. Spiral-like magnetic fields are associated with envelopes having large velocity gradients, indicating that the rotation motions are strong enough to twist the field lines. All of the protostars with a standard-hourglass field morphology show no significant velocity gradient due to the strong magnetic braking. △ Less

Submitted 15 May, 2024; v1 submitted 11 February, 2024; originally announced February 2024.

Comments: 12 pages, 5 figures (16 pages and 7 figures for the supplementary materials) Accepted for publication on ApJ Letters

Journal ref: ApJL, 963:L31, 2024

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)

Early Planet Formation in Embedded Disks (eDisk) XIV: Flared Dust Distribution and Viscous Accretion Heating of the Disk around R CrA IRS 7B-a

Authors: Shigehisa Takakuwa, Kazuya Saigo, Miyu Kido, Nagayoshi Ohashi, John J. Tobin, Jes K. Jørgensen, Yuri Aikawa, Yusuke Aso, Sacha Gavino, Ilseung Han, Patrick M. Koch, Woojin Kwon, Chang Won Lee, Jeong-Eun Lee, Zhi-Yun Li, Zhe-Yu Daniel Lin, Leslie W. Looney, Shoji Mori, Jinshi Sai, Rajeeb Sharma, Patrick Sheehan, Kengo Tomida, Jonathan P. Williams, Yoshihide Yamato, Hsi-Wei Yen

Abstract: We performed radiative transfer calculations and observing simulations to reproduce the 1.3-mm dust-continuum and C$^{18}$O (2-1) images in the Class I protostar R CrA IRS7B-a, observed with the ALMA Large Program ``Early Planet Formation in Embedded Disks (eDisk)". We found that the dust disk model passively heated by the central protostar cannot reproduce the observed peak brightness temperature… ▽ More We performed radiative transfer calculations and observing simulations to reproduce the 1.3-mm dust-continuum and C$^{18}$O (2-1) images in the Class I protostar R CrA IRS7B-a, observed with the ALMA Large Program ``Early Planet Formation in Embedded Disks (eDisk)". We found that the dust disk model passively heated by the central protostar cannot reproduce the observed peak brightness temperature of the 1.3-mm continuum emission ($\sim$195 K), regardless of the assumptions about the dust opacity. Our calculation suggests that viscous accretion heating in the disk is required to reproduce the observed high brightness temperature. The observed intensity profile of the 1.3-mm dust-continuum emission along the disk minor axis is skewed toward the disk far side. Our modeling reveals that such an asymmetric intensity distribution requires flaring of the dust along the disk's vertical direction with the scale-height following $h/r \sim r^{0.3}$ as function of radius. These results are in sharp contrast to those of Class II disks, which show geometrically flat dust distributions and lower dust temperatures. From our modeling of the C$^{18}$O (2-1) emission, the outermost radius of the gas disk is estimated to be $\sim$80 au, larger than that of the dust disk ($\sim$62 au), to reproduce the observed distribution of the C$^{18}$O (2-1) emission in IRS 7B-a. Our modeling unveils a hot and thick dust disk plus a larger gas disk around one of the eDisk targets, which could be applicable to other protostellar sources in contrast to more evolved sources. △ Less

Submitted 16 January, 2024; originally announced January 2024.

Comments: 26 pages, 13 figures

The Disk Orientations of Perseus Protostellar Multiples at 8 au Resolution

Authors: Nickalas K. Reynolds, John J. Tobin, Patrick D. Sheehan, Sarah I. Sadavoy, Leslie W. Looney, Kaitlin M. Kratter, Zhi-Yun Li, Dominique M. Segura-Cox, Nathan A. Kaib

Abstract: We present a statistical characterization of circumstellar disk orientations toward 12 protostellar multiple systems in the Perseus molecular cloud using the Atacama Large Millimeter/submillimeter Array at Band 6 (1.3 mm) with a resolution of 25 mas (8 au). This exquisite resolution enabled us to resolve the compact inner disk structures surrounding the components of each multiple system and to de… ▽ More We present a statistical characterization of circumstellar disk orientations toward 12 protostellar multiple systems in the Perseus molecular cloud using the Atacama Large Millimeter/submillimeter Array at Band 6 (1.3 mm) with a resolution of 25 mas (8 au). This exquisite resolution enabled us to resolve the compact inner disk structures surrounding the components of each multiple system and to determine the projected 3-D orientation of the disks (position angle and inclination) to high precision. We performed a statistical analysis on the relative alignment of disk pairs to determine whether the disks are preferentially aligned or randomly distributed. We considered three subsamples of the observations selected by the companion separations, a <100 au, a >500 au, and a < 10,000 au. We found for the compact (< 100 au) subsample, the distribution of orientation angles is best described by an underlying distribution of preferentially aligned sources (within 30deg) but does not rule out distributions with 40% misaligned sources. The wide companion (>500 au) subsample appears to be consistent with a distribution of 40%-80% preferentially aligned sources. Similarly, the full sample of systems with companions (a< 10, 000 au) is most consistent with a fractional ratio of at most 80% preferentially aligned source and rules out purely randomly aligned distributions. Thus our results imply the compact sources (<100 au) and the wide companions (>500 au) are statistically different. △ Less

Submitted 20 December, 2023; originally announced December 2023.

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

The CO outflow components ejected by a recent accretion event in B335

Authors: Chul-Hwan Kim, Jeong-Eun Lee, Carlos Contreras Peña, Doug Johnstone, Gregory J. Herczeg, John J. Tobin, Neal J. Evans II

Abstract: Protostellar outflows often present a knotty appearance, providing evidence of sporadic accretion in stellar mass growth. To understand the direct relation between mass accretion and ejection, we analyze the contemporaneous accretion activity and associated ejection components in B335. B335 has brightened in the mid-IR by 2.5 mag since 2010, indicating increased luminosity, presumably due to incre… ▽ More Protostellar outflows often present a knotty appearance, providing evidence of sporadic accretion in stellar mass growth. To understand the direct relation between mass accretion and ejection, we analyze the contemporaneous accretion activity and associated ejection components in B335. B335 has brightened in the mid-IR by 2.5 mag since 2010, indicating increased luminosity, presumably due to increased mass accretion rate onto the protostar. ALMA observations of 12CO emission in the outflow reveal high-velocity emission, estimated to have been ejected 4.6 - 2 years before the ALMA observation and consistent with the jump in mid-IR brightness. The consistency in timing suggests that the detected high-velocity ejection components are directly linked to the most recent accretion activity. We calculated the kinetic energy, momentum, and force for the ejection component associated with the most recent accretion activity and found that at least, about 1.0% of accreted mass has been ejected. More accurate information on the jet inclination and the temperature of the ejected gas components will better constrain the ejected mass induced by the recently enhanced accretion event. △ Less

Submitted 10 December, 2023; originally announced December 2023.

Comments: 11 pages, 5 figures, and 1 table

Early Planet Formation in Embedded Disks (eDisk) X: Compact Disks, Extended Infall, and a Fossil Outburst in the Class I Oph IRS43 Binary

Authors: Suchitra Narayanan, Jonathan P. Williams, John J. Tobin, Jes K. Jorgensen, Nagayoshi Ohashi, Zhe-Yu Daniel Lin, Merel L. R. van't Hoff, Zhi-Yun Li, Adele L. Plunkett, Leslie W. Looney, Shigehisa Takakuwa, Hsi-Wei Yen, Yusuke Aso, Christian Flores, Jeong-Eun Lee, Shih-Ping Lai, Woojin Kwon, Itziar de Gregorio-Monsalvo, Rajeeb Sharma, Chang Won Lee

Abstract: We present the first results from the Early Planet Formation in Embedded Disks (eDisk) ALMA Large Program toward Oph IRS43, a binary system of solar mass protostars. The 1.3 mm dust continuum observations resolve a compact disk, ~6au radius, around the northern component and show that the disk around the southern component is even smaller, <~3 au. CO, 13CO, and C18O maps reveal a large cavity in a… ▽ More We present the first results from the Early Planet Formation in Embedded Disks (eDisk) ALMA Large Program toward Oph IRS43, a binary system of solar mass protostars. The 1.3 mm dust continuum observations resolve a compact disk, ~6au radius, around the northern component and show that the disk around the southern component is even smaller, <~3 au. CO, 13CO, and C18O maps reveal a large cavity in a low mass envelope that shows kinematic signatures of rotation and infall extending out to ~ 2000au. An expanding CO bubble centered on the extrapolated location of the source ~130 years ago suggests a recent outburst. Despite the small size of the disks, the overall picture is of a remarkably large and dynamically active region. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: Paper 10 of the ALMA eDisk Large Program. Accepted for publication in ApJ

Early Planet Formation in Embedded Disks (eDisk) XII: Accretion streamers, protoplanetary disk, and outflow in the Class I source Oph IRS63

Authors: Christian Flores, Nagayoshi Ohashi, John J. Tobin, Jes K. Jørgensen, Shigehisa Takakuwa, Zhi-Yun Li, Zhe-Yu Daniel Lin, Merel L. R. van 't Hoff, Adele L. Plunkett, Yoshihide Yamato, Jinshi Sai, Patrick M. Koch, Hsi-Wei Yen, Yuri Aikawa, Yusuke Aso, Itziar de Gregorio-Monsalvo, Miyu Kido, Woojin Kwon, Jeong-Eun Lee, Chang Won Lee, Leslie W. Looney, Alejandro Santamaría-Miranda, Rajeeb Sharma, Travis J. Thieme, Jonathan P. Williams , et al. (3 additional authors not shown)

Abstract: We present ALMA observations of the Class I source Oph IRS63 in the context of the Early Planet Formation in Embedded Disks (eDisk) large program. Our ALMA observations of Oph IRS63 show a myriad of protostellar features, such as a shell-like bipolar outflow (in $^{12}$CO), an extended rotating envelope structure (in $^{13}$CO), a streamer connecting the envelope to the disk (in C$^{18}$O), and se… ▽ More We present ALMA observations of the Class I source Oph IRS63 in the context of the Early Planet Formation in Embedded Disks (eDisk) large program. Our ALMA observations of Oph IRS63 show a myriad of protostellar features, such as a shell-like bipolar outflow (in $^{12}$CO), an extended rotating envelope structure (in $^{13}$CO), a streamer connecting the envelope to the disk (in C$^{18}$O), and several small-scale spiral structures seen towards the edge of the dust continuum (in SO). By analyzing the velocity pattern of $^{13}$CO and C$^{18}$O, we measure a protostellar mass of $\rm M_\star = 0.5 \pm 0.2 $~$\rm M_\odot$ and confirm the presence of a disk rotating at almost Keplerian velocity that extends up to $\sim260$ au. These calculations also show that the gaseous disk is about four times larger than the dust disk, which could indicate dust evolution and radial drift. Furthermore, we model the C$^{18}$O streamer and SO spiral structures as features originating from an infalling rotating structure that continuously feeds the young protostellar disk. We compute an envelope-to-disk mass infall rate of $\sim 10^{-6}$~$\rm M_\odot \, yr^{-1}$ and compare it to the disk-to-star mass accretion rate of $\sim 10^{-8}$~$\rm M_\odot \, yr^{-1}$, from which we infer that the protostellar disk is in a mass build-up phase. At the current mass infall rate, we speculate that soon the disk will become too massive to be gravitationally stable. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: 26 pages and 17 figures

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

Early Planet Formation in Embedded Disks (eDisk). VIII. A Small Protostellar Disk around the Extremely Low-Mass and Young Class 0 Protostar, IRAS 15398-3359

Authors: Travis J. Thieme, Shih-Ping Lai, Nagayoshi Ohashi, John J. Tobin, Jes K. Jørgensen, Jinshi Sai, Yusuke Aso, Jonathan P. Williams, Yoshihide Yamato, Yuri Aikawa, Itziar de Gregorio-Monsalvo, Ilseung Han, Woojin Kwon, Chang Won Lee, Jeong-Eun Lee, Zhi-Yun Li, Zhe-Yu Daniel Lin, Leslie W. Looney, Suchitra Narayanan, Nguyen Thi Phuong, Adele L. Plunkett, Alejandro Santamaría-Miranda, Rajeeb Sharma, Shigehisa Takakuwa, Hsi-Wei Yen

Abstract: Protostellar disks are a ubiquitous part of the star formation process and the future sites of planet formation. As part of the Early Planet Formation in Embedded Disks (eDisk) large program, we present high-angular resolution dust continuum ($\sim40\,$mas) and molecular line ($\sim150\,$mas) observations of the Class 0 protostar, IRAS 15398-3359. The dust continuum is small, compact, and centrall… ▽ More Protostellar disks are a ubiquitous part of the star formation process and the future sites of planet formation. As part of the Early Planet Formation in Embedded Disks (eDisk) large program, we present high-angular resolution dust continuum ($\sim40\,$mas) and molecular line ($\sim150\,$mas) observations of the Class 0 protostar, IRAS 15398-3359. The dust continuum is small, compact, and centrally peaked, while more extended dust structures are found in the outflow directions. We perform a 2D Gaussian fitting to find the deconvolved size and $2σ$ radius of the dust disk to be $4.5\times2.8\,\mathrm{au}$ and $3.8\,\mathrm{au}$, respectively. We estimate the gas+dust disk mass assuming optically thin continuum emission to be $0.6-1.8\,M_\mathrm{jup}$, indicating a very low-mass disk. The CO isotopologues trace components of the outflows and inner envelope, while SO traces a compact, rotating disk-like component. Using several rotation curve fittings on the PV diagram of the SO emission, the lower limits of the protostellar mass and gas disk radius are $0.022\,M_\odot$ and $31.2\,\mathrm{au}$ from our Modified 2 single power-law fitting. A conservative upper limit of the protostellar mass is inferred to be $0.1\,M_\odot$. The protostellar mass-accretion rate and the specific angular momentum at the protostellar disk edge are found to be between $1.3-6.1\times10^{-6}\,M_\odot\,\mathrm{yr^{-1}}$ and $1.2-3.8\times10^{-4}\,\mathrm{km\,s^{-1}\,pc}$, respectively, with an age estimated between $0.4-7.5\times10^{4}\,$yr. At this young age with no clear substructures in the disk, planet formation would likely not yet have started. This study highlights the importance of high-resolution observations and systematic fitting procedures when deriving dynamical properties of deeply embedded Class 0 protostars. △ Less

Submitted 19 October, 2023; originally announced October 2023.

Comments: 28 pages, 16 figures. Accepted for publication in ApJ as one of the first-look papers of the eDisk ALMA Large Program

IRAS4A1: Multi-wavelength continuum analysis of a very flared Class 0 disk

Authors: O. M. Guerra-Alvarado, N. van der Marel, J. Di Francesco, L. W. Looney, J. J. Tobin, E. G. Cox, P. D. Sheehan, D. J. Wilner, E. Macías, C. Carrasco-González

Abstract: Understanding the formation of substructures in protoplanetary disks is vital for gaining insights into dust growth and the process of planet formation. Studying these substructures in highly embedded Class 0 objects using the Atacama Large Millimeter/submillimeter Array (ALMA), however, poses significant challenges. Nonetheless, it is imperative to do so to unravel the mechanisms and timing behin… ▽ More Understanding the formation of substructures in protoplanetary disks is vital for gaining insights into dust growth and the process of planet formation. Studying these substructures in highly embedded Class 0 objects using the Atacama Large Millimeter/submillimeter Array (ALMA), however, poses significant challenges. Nonetheless, it is imperative to do so to unravel the mechanisms and timing behind the formation of these substructures. In this study, we present high-resolution ALMA data at Bands 6 and 4 of the NGC1333 IRAS4A Class 0 protobinary system. This system consists of two components, A1 and A2, separated by 1.8" and located in the Perseus molecular cloud at $\sim$293 pc distance. To gain a comprehensive understanding of the dust properties and formation of substructures in the early stages, we conducted a multi-wavelength analysis of IRAS4A1. Additionally, we sought to address whether the lack of observed substructures in very young disks, could be attributed to factors such as high degrees of disk flaring and large scale heights. To explore this phenomenon, we employed radiative transfer models using RADMC-3D. Our multi-wavelength analysis of A1 discovered characteristics such as high dust surface density, substantial dust mass within the disk, and elevated dust temperatures. These findings suggest the presence of large dust grains compared to the ones in the interstellar medium (ISM), greater than 100 microns in size within the region. Furthermore, while there's no direct detection of any substructure, our models indicate that some, such as a small gap, must be present. In summary, this result implies that disk substructures may be masked or obscured by a large scale height in combination with a high degree of flaring in Class 0 disks. [Abridged] △ Less

Submitted 18 October, 2023; originally announced October 2023.

Comments: Accepted for publication in A&A, 13 pages, 7 figures

Journal ref: A&A 681, A82 (2024)

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)

Finding substructures in protostellar disks in Ophiuchus

Authors: Arnaud Michel, Sarah I. Sadavoy, Patrick D. Sheehan, Leslie W. Looney, Erin G. Cox, John J. Tobin, Nienke van der Marel, Dominique M. Segura-Cox

Abstract: High-resolution, millimeter observations of disks at the protoplanetary stage reveal substructures such as gaps, rings, arcs, spirals, and cavities. While many protoplanetary disks host such substructures, only a few at the younger protostellar stage have shown similar features. We present a detailed search for early disk substructures in ALMA 1.3 and 0.87~mm observations of ten protostellar disks… ▽ More High-resolution, millimeter observations of disks at the protoplanetary stage reveal substructures such as gaps, rings, arcs, spirals, and cavities. While many protoplanetary disks host such substructures, only a few at the younger protostellar stage have shown similar features. We present a detailed search for early disk substructures in ALMA 1.3 and 0.87~mm observations of ten protostellar disks in the Ophiuchus star-forming region. Of this sample, four disks have identified substructure, two appear to be smooth disks, and four are considered ambiguous. The structured disks have wide Gaussian-like rings ($σ_R/R_{\mathrm{disk}}\sim0.26$) with low contrasts ($C<0.2$) above a smooth disk profile, in comparison to protoplanetary disks where rings tend to be narrow and have a wide variety of contrasts ($σ_R/R_{\mathrm{disk}}\sim0.08$ and $C$ ranges from $0-1$). The four protostellar disks with the identified substructures are among the brightest sources in the Ophiuchus sample, in agreement with trends observed for protoplanetary disks. These observations indicate that substructures in protostellar disks may be common in brighter disks. The presence of substructures at the earliest stages suggests an early start for dust grain growth and, subsequently, planet formation. The evolution of these protostellar substructures is hypothesized in two potential pathways: (1) the rings are the sites of early planet formation, and the later observed protoplanetary disk ring-gap pairs are secondary features, or (2) the rings evolve over the disk lifetime to become those observed at the protoplanetary disk stage. △ Less

Submitted 14 September, 2023; originally announced September 2023.

Comments: Accepted by ApJ, 22 pages, 10 figures

Early Planet Formation in Embedded Disks (eDisk) VI: Kinematic Structures around the Very Low Mass Protostar IRAS 16253-2429

Authors: Yusuke Aso, Woojin Kwon, Nagayoshi Ohashi, Jes K. Jorgensen, John J. Tobin, Yuri Aikawa, Itziar de Gregorio-Monsalvo, Ilseung Han, Miyu Kido, Patrick M. Koch, Shih-Ping Lai, Chang Won Lee, Jeong-Eun Lee, Zhi-Yun Li, Zhe-Yu Daniel Lin, Leslie W. Looney, Suchitra Narayanan, Nguyen Thi Phuong, Jinshi Sai, Kazuya Saigo, Alejandro Santamaria-Miranda, Rajeeb Sharma, Shigehisa Takakuwa, Travis J. Thieme, Kengo Tomida , et al. (2 additional authors not shown)

Abstract: Precise estimates of protostellar masses are crucial to characterize the formation of stars of low masses down to brown-dwarfs (BDs; M* < 0.08 Msun). The most accurate estimation of protostellar mass uses the Keplerian rotation in the circumstellar disk around the protostar. To apply the Keplerian rotation method to a protostar at the low-mass end, we have observed the Class 0 protostar IRAS 16253… ▽ More Precise estimates of protostellar masses are crucial to characterize the formation of stars of low masses down to brown-dwarfs (BDs; M* < 0.08 Msun). The most accurate estimation of protostellar mass uses the Keplerian rotation in the circumstellar disk around the protostar. To apply the Keplerian rotation method to a protostar at the low-mass end, we have observed the Class 0 protostar IRAS 16253-2429 using the Atacama Large Millimeter/submillimeter Array (ALMA) in the 1.3 mm continuum at an angular resolution of 0.07" (10 au), and in the 12CO, C18O, 13CO (J=2-1), and SO (J_N = 6_5-5_4) molecular lines, as part of the ALMA Large Program Early Planet Formation in Embedded Disks (eDisk). The continuum emission traces a non-axisymmetric, disk-like structure perpendicular to the associated 12CO outflow. The position-velocity (PV) diagrams in the C18O and 13CO lines can be interpreted as infalling and rotating motions. In contrast, the PV diagram along the major axis of the disk-like structure in the 12CO line allows us to identify Keplerian rotation. The central stellar mass and the disk radius are estimated to be ~0.12-0.17 Msun and ~13-19 au, respectively. The SO line suggests the existence of an accretion shock at a ring (r~28 au) surrounding the disk and a streamer from the eastern side of the envelope. IRAS 16253-2429 is not a proto-BD but has a central stellar mass close to the BD mass regime, and our results provide a typical picture of such very low-mass protostars. △ Less

Submitted 4 September, 2023; originally announced September 2023.

Comments: 41 pages, 14 figures

Journal ref: ApJ 954 (2023) 101-117

Early Planet Formation in Embedded Disks (eDisk) IX: High-resolution ALMA Observations of the Class 0 Protostar R CrA IRS5N and its surrounding

Authors: Rajeeb Sharma, Jes K. Jørgensen, Sacha Gavino, Nagayoshi Ohashi, John J. Tobin, Zhe-Yu Daniel Lin, Zhi-Yun Li, Shigehisa Takakuwa, Chang Won Lee, Jinshi Sai, Woojin Kwon, Itziar de Gregorio-Monsalvo, Alejandro Santamaría-Miranda, Hsi-Wei Yen, Yuri Aikawa, Yusuke Aso, Shih-Ping Lai, Jeong-Eun Lee, Leslie W. Looney, Nguyen Thi Phuong, Travis J. Thieme, Jonathan P. Williams

Abstract: We present high-resolution, high-sensitivity observations of the Class 0 protostar RCrA IRS5N as part of the Atacama Large Milimeter/submilimeter Array (ALMA) large program Early Planet Formation in Embedded Disks (eDisk). The 1.3 mm continuum emission reveals a flattened continuum structure around IRS5N, consistent with a protostellar disk in the early phases of evolution. The continuum emission… ▽ More We present high-resolution, high-sensitivity observations of the Class 0 protostar RCrA IRS5N as part of the Atacama Large Milimeter/submilimeter Array (ALMA) large program Early Planet Formation in Embedded Disks (eDisk). The 1.3 mm continuum emission reveals a flattened continuum structure around IRS5N, consistent with a protostellar disk in the early phases of evolution. The continuum emission appears smooth and shows no substructures. However, a brightness asymmetry is observed along the minor axis of the disk, suggesting the disk is optically and geometrically thick. We estimate the disk mass to be between 0.007 and 0.02 M$_{\odot}$. Furthermore, molecular emission has been detected from various species, including C$^{18}$O (2$-$1), $^{12}$CO (2$-$1), $^{13}$CO (2$-$1), and H$_2$CO (3$_{0,3}-2_{0,2}$, 3$_{2,1}-2_{2,0}$, and 3$_{2,2}-2_{2,1}$). By conducting a position-velocity analysis of the C$^{18}$O (2$-$1) emission, we find that the disk of IRS5N exhibits characteristics consistent with Keplerian rotation around a central protostar with a mass of approximately 0.3 M$_{\odot}$. Additionally, we observe dust continuum emission from the nearby binary source, IRS5a/b. The emission in $^{12}$CO toward IRS5a/b seems to emanate from IRS5b and flow into IRS5a, suggesting material transport between their mutual orbits. The lack of a detected outflow and large-scale negatives in \tlvco~observed toward IRS5N suggests that much of the flux from IRS5N is being resolved out. Due to this substantial surrounding envelope, the central IRS5N protostar is expected to be significantly more massive in the future. △ Less

Submitted 1 September, 2023; originally announced September 2023.

Comments: 25 pages, 24 figures. Accepted for publication in ApJ

Early Planet Formation in Embedded Disks (eDisk) V: Possible Annular Substructure in a Circumstellar Disk in the Ced110 IRS4 System

Authors: Jinshi Sai, Hsi-Wei Yen, Nagayoshi Ohashi, John J. Tobin, Jes K. Jørgensen, Shigehisa Takakuwa, Kazuya Saigo, Yusuke Aso, Zhe-Yu Daniel Lin, Patrick M. Koch, Yuri Aikawa, Christian Flores, Itziar de Gregorio-Monsalvo, Ilseung Han, Miyu Kido, Woojin Kwon, Shih-Ping Lai, Chang Won Lee, Jeong-Eun Lee, Zhi-Yun Li, Leslie W. Looney, Shoji Mori, Nguyen Thi Phuong, Alejandro Santamaría-Miranda, Rajeeb Sharma , et al. (3 additional authors not shown)

Abstract: We have observed the Class 0/I protostellar system Ced110 IRS4 at an angular resolution of $0.05''$ ($\sim$10 au) as a part of the ALMA large program; Early Planet Formation in the Embedded Disks (eDisk). The 1.3 mm dust continuum emission reveals that Ced110 IRS4 is a binary system with a projected separation of $\sim$250 au. The continuum emissions associated with the main source and its compani… ▽ More We have observed the Class 0/I protostellar system Ced110 IRS4 at an angular resolution of $0.05''$ ($\sim$10 au) as a part of the ALMA large program; Early Planet Formation in the Embedded Disks (eDisk). The 1.3 mm dust continuum emission reveals that Ced110 IRS4 is a binary system with a projected separation of $\sim$250 au. The continuum emissions associated with the main source and its companion, named Ced110 IRS4A and IRS4B respectively, exhibit disk-like shapes and likely arise from dust disks around the protostars. The continuum emission of Ced110 IRS4A has a radius of $\sim$91.7 au ($\sim0.485''$), and shows bumps along its major axis with an asymmetry. The bumps can be interpreted as an shallow, ring-like structure at a radius of $\sim$40 au ($\sim0.2''$) in the continuum emission, as demonstrated from two-dimensional intensity distribution models. A rotation curve analysis on the C$^{18}$O and $^{13}$CO $J=2$-1 lines reveals the presence of a Keplerian disk within a radius of 120 au around Ced110 IRS4A, which supports the interpretation that the dust continuum emission arises from a disk. The ring-like structure in the dust continuum emission might indicate a possible, annular substructure in the surface density of the embedded disk, although the possibility that it is an apparent structure due to the optically thick continuum emission cannot be ruled out. △ Less

Submitted 31 August, 2023; v1 submitted 17 July, 2023; originally announced July 2023.

Comments: 32 pages, 23 figures. Accepted for publication in ApJ as one of the first-look papers of the eDisk ALMA Large Program

Complex Organic Molecules in a Very Young Hot Corino, HOPS 373SW

Authors: Jeong-Eun Lee, Giseon Baek, Seokho Lee, Jae-Hong Jeong, Chul-Hwan Kim, Yuri Aikawa, Gregory J. Herczeg, Doug Johnstone, John J. Tobin

Abstract: We present the spectra of Complex Organic Molecules (COMs) detected in HOPS 373SW with the Atacama Large Millimeter/submillimeter Array (ALMA). HOPS 373SW, which is a component of a protostellar binary with a separation of 1500 au, has been discovered as a variable protostar by the JCMT Transient monitoring survey with a modest ~30% brightness increase at submillimeter wavelengths. Our ALMA Target… ▽ More We present the spectra of Complex Organic Molecules (COMs) detected in HOPS 373SW with the Atacama Large Millimeter/submillimeter Array (ALMA). HOPS 373SW, which is a component of a protostellar binary with a separation of 1500 au, has been discovered as a variable protostar by the JCMT Transient monitoring survey with a modest ~30% brightness increase at submillimeter wavelengths. Our ALMA Target of Opportunity (ToO) observation at ~345 GHz for HOPS 373SW revealed extremely young chemical characteristics with strong deuteration of methanol. The dust continuum opacity is very high toward the source center, obscuring line emission from within 0.03 arcsec. The other binary component, HOPS 373NE, was detected only in C17O in our observation, implying a cold and quiescent environment. We compare the COMs abundances relative to CH3OH in HOPS 373SW with those of V883 Ori, which is an eruptive disk object, as well as other hot corinos, to demonstrate the chemical evolution from envelope to disk. High abundances of singly, doubly, and triply deuterated methanol (CH2DOH, CHD2OH, and CD3OH) and a low CH3CN abundance in HOPS 373SW compared to other hot corinos suggest a very early evolutionary stage of HOPS 373SW in the hot corino phase. Since the COMs detected in HOPS 373SW would have been sublimated very recently from grain surfaces, HOPS 373SW is a promising place to study the surface chemistry of COMs in the cold prestellar phase, before sublimation. △ Less

Submitted 29 June, 2023; originally announced June 2023.

Comments: Accepted for publication in ApJ

Early Planet Formation in Embedded Disks (eDisk). VII. Keplerian Disk, Disk Substructure, and Accretion Streamers in the Class 0 Protostar IRAS 16544-1604 in CB 68

Authors: Miyu Kido, Shigehisa Takakuwa, Kazuya Saigo, Nagayoshi Ohashi, John J. Tobin, Jes K, Jørgensen, Yuri Aikawa, Yusuke Aso, Frankie J. Encalada, Christian Flores, Sacha Gavino, Itziar de Gregorio-Monsalvo, Ilseung Han, Shingo Hirano, Patrick M. Koch, Woojin Kwon, Shih-Ping Lai, Chang Won Lee, Jeong-Eun Lee, Zhi-Yun Li, Zhe-Yu Daniel Lin, Leslie W. Looney, Shoji Mori, Suchitra Narayanan , et al. (12 additional authors not shown)

Abstract: We present observations of the Class 0 protostar IRAS 16544-1604 in CB 68 from the ''Early Planet Formation in Embedded Disks (eDisk)'' ALMA Large program. The ALMA observations target continuum and lines at 1.3-mm with an angular resolution of $\sim$5 au. The continuum image reveals a dusty protostellar disk with a radius of $\sim$30 au seen close to edge-on, and asymmetric structures both along… ▽ More We present observations of the Class 0 protostar IRAS 16544-1604 in CB 68 from the ''Early Planet Formation in Embedded Disks (eDisk)'' ALMA Large program. The ALMA observations target continuum and lines at 1.3-mm with an angular resolution of $\sim$5 au. The continuum image reveals a dusty protostellar disk with a radius of $\sim$30 au seen close to edge-on, and asymmetric structures both along the major and minor axes. While the asymmetry along the minor axis can be interpreted as the effect of the dust flaring, the asymmetry along the major axis comes from a real non-axisymmetric structure. The C$^{18}$O image cubes clearly show the gas in the disk that follows a Keplerian rotation pattern around a $\sim$0.14 $M_{\odot}$ central protostar. Furthermore, there are $\sim$1500 au-scale streamer-like features of gas connecting from North-East, North-North-West, and North-West to the disk, as well as the bending outflow as seen in the $^{12}$CO (2-1) emission. At the apparent landing point of NE streamer, there are SO (6$_5$-5$_4$) and SiO (5-4) emission detected. The spatial and velocity structure of NE streamer can be interpreted as a free-falling gas with a conserved specific angular momentum, and the detection of the SO and SiO emission at the tip of the streamer implies presence of accretion shocks. Our eDisk observations have unveiled that the Class 0 protostar in CB 68 has a Keplerian rotating disk with flaring and non-axisymmetric structure associated with accretion streamers and outflows. △ Less

Submitted 27 June, 2023; originally announced June 2023.

Comments: 30 pages, 24 figures, accepted for publication in The Astrophysical Journal as one of the first-look papers of the eDisk ALMA Large Program

Early Planet Formation in Embedded Disks (eDisk). II. Limited Dust Settling and Prominent Snow Surfaces in the Edge-on Class I Disk IRAS 04302+2247

Authors: Zhe-Yu Daniel Lin, Zhi-Yun Li, John J. Tobin, Nagayoshi Ohashi, Jes Kristian Jørgensen, Leslie W. Looney, Yusuke Aso, Shigehisa Takakuwa, Yuri Aikawa, Merel L. R. van 't Hoff, Itziar de Gregorio-Monsalvo, Frankie J. Encalada, Christian Flores, Sacha Gavino, Ilseung Han, Miyu Kido, Patrick M. Koch, Woojin Kwon, Shih-Ping Lai, Chang Won Lee, Jeong-Eun Lee, Nguyen Thi Phuong, Jinshi Sai, Rajeeb Sharma, Patrick Sheehan , et al. (4 additional authors not shown)

Abstract: While dust disks around optically visible, Class II protostars are found to be vertically thin, when and how dust settles to the midplane are unclear. As part of the Atacama Large Millimeter/submillimeter Array (ALMA) large program, Early Planet Formation in Embedded Disks, we analyze the edge-on, embedded, Class I protostar IRAS 04302+2247, also nicknamed the ``Butterfly Star." With a resolution… ▽ More While dust disks around optically visible, Class II protostars are found to be vertically thin, when and how dust settles to the midplane are unclear. As part of the Atacama Large Millimeter/submillimeter Array (ALMA) large program, Early Planet Formation in Embedded Disks, we analyze the edge-on, embedded, Class I protostar IRAS 04302+2247, also nicknamed the ``Butterfly Star." With a resolution of 0.05" (8~au), the 1.3 mm continuum shows an asymmetry along the minor axis which is evidence of an optically thick and geometrically thick disk viewed nearly edge-on. There is no evidence of rings and gaps, which could be due to the lack of radial substructure or the highly inclined and optically thick view. With 0.1" (16~au) resolution, we resolve the 2D snow surfaces, i.e., the boundary region between freeze-out and sublimation, for $^{12}$CO $J$=2--1, $^{13}$CO $J$=2--1, C$^{18}$O $J$=2--1, $H_{2}$CO $J$=$3_{0,3}$--$2_{0,2}$, and SO $J$=$6_{5}$--$5_{4}$, and constrain the CO midplane snow line to $\sim 130$ au. We find Keplerian rotation around a protostar of $1.6 \pm 0.4 M_{\odot}$ using C$^{18}$O. Through forward ray-tracing using RADMC-3D, we find that the dust scale height is $\sim 6$ au at a radius of 100~au from the central star and is comparable to the gas pressure scale height. The results suggest that the dust of this Class~I source has yet to vertically settle significantly. △ Less

Submitted 27 June, 2023; originally announced June 2023.

Comments: 33 pages, 21 figures. Accepted for publication in ApJ as one of the first-look papers of the eDisk ALMA Large Program

Early Planet Formation in Embedded Disks (eDisk). IV. The Ringed and Warped Structure of the Disk around the Class I Protostar L1489 IRS

Authors: Yoshihide Yamato, Yuri Aikawa, Nagayoshi Ohashi, John J. Tobin, Jes K. Jørgensen, Shigehisa Takakuwa, Yusuke Aso, Jinshi Sai, Christian Flores, Itziar de Gregorio-Monsalvo, Shingo Hirano, Ilseung Han, Miyu Kido, Patrick M. Koch, Woojin Kwon, Shih-Ping Lai, Chang Won Lee, Jeong-Eun Lee, Zhi-Yun Li, Zhe-Yu Daniel Lin, Leslie W. Looney, Shoji Mori, Suchitra Narayanan, Nguyen Thi Phuong, Kazuya Saigo , et al. (6 additional authors not shown)

Abstract: Constraining the physical and chemical structure of young embedded disks is crucial to understanding the earliest stages of planet formation. As part of the Early Planet Formation in Embedded Disks Atacama Large Millimeter/submillimeter Array Large Program, we present high spatial resolution ($\sim$0$.\!\!^{\prime\prime}$1 or $\sim$15 au) observations of the 1.3 mm continuum and $^{13}$CO $J=$ 2-1… ▽ More Constraining the physical and chemical structure of young embedded disks is crucial to understanding the earliest stages of planet formation. As part of the Early Planet Formation in Embedded Disks Atacama Large Millimeter/submillimeter Array Large Program, we present high spatial resolution ($\sim$0$.\!\!^{\prime\prime}$1 or $\sim$15 au) observations of the 1.3 mm continuum and $^{13}$CO $J=$ 2-1, C$^{18}$O $J=$ 2-1, and SO $J_N=$ $6_5$-$5_4$ molecular lines toward the disk around the Class I protostar L1489 IRS. The continuum emission shows a ring-like structure at 56 au from the central protostar and a tenuous, optically thin emission extending beyond $\sim$300 au. The $^{13}$CO emission traces the warm disk surface, while the C$^{18}$O emission originates from near the disk midplane. The coincidence of the radial emission peak of C$^{18}$O with the dust ring may indicate a gap-ring structure in the gaseous disk as well. The SO emission shows a highly complex distribution, including a compact, prominent component at $\lesssim$30 au, which is likely to originate from thermally sublimated SO molecules. The compact SO emission also shows a velocity gradient along a slightly ($\sim15^\circ$) tilted direction with respect to the major axis of the dust disk, which we interpret as an inner warped disk in addition to the warp around $\sim$200 au suggested by previous work. These warped structures may be formed by a planet or companion with an inclined orbit, or by a gradual change in the angular momentum axis during gas infall. △ Less

Submitted 27 June, 2023; originally announced June 2023.

Comments: 24 pages, 12 figures. Accepted for publication in The Astrophysical Journal as one of the first-look papers of the eDisk ALMA Large Program

Early Planet Formation in Embedded Disks (eDisk) III: A first high-resolution view of sub-mm continuum and molecular line emission toward the Class 0 protostar L1527 IRS

Authors: Merel L. R. van 't Hoff, John J. Tobin, Zhi-Yun Li, Nagayoshi Ohashi, Jes K. Jørgensen, Zhe-Yu Daniel Lin, Yuri Aikawa, Yusuke Aso, Itziar de Gregorio-Monsalvo, Sacha Gavino, Ilseung Han, Patrick M. Koch, Woojin Kwon, Chang Won Lee, Jeong-Eun Lee, Leslie W. Looney, Suchitra Narayanan, Adele Plunkett, Jinshi Sai, Alejandro Santamaría-Miranda, Rajeeb Sharma, Patrick D. Sheehan, Shigehisa Takakuwa, Travis J. Thieme, Jonathan P. Williams , et al. (3 additional authors not shown)

Abstract: Studying the physical and chemical conditions of young embedded disks is crucial to constrain the initial conditions for planet formation. Here, we present Atacama Large Millimeter/submillimeter Array (ALMA) observations of dust continuum at $\sim$0.06" (8 au) resolution and molecular line emission at $\sim$0.17" (24 au) resolution toward the Class 0 protostar L1527 IRS from the Large Program eDis… ▽ More Studying the physical and chemical conditions of young embedded disks is crucial to constrain the initial conditions for planet formation. Here, we present Atacama Large Millimeter/submillimeter Array (ALMA) observations of dust continuum at $\sim$0.06" (8 au) resolution and molecular line emission at $\sim$0.17" (24 au) resolution toward the Class 0 protostar L1527 IRS from the Large Program eDisk (Early Planet Formation in Embedded Disks). The continuum emission is smooth without substructures, but asymmetric along both the major and minor axes of the disk as previously observed. The detected lines of $^{12}$CO, $^{13}$CO, C$^{18}$O, H$_2$CO, c-C$_3$H$_2$, SO, SiO, and DCN trace different components of the protostellar system, with a disk wind potentially visible in $^{12}$CO. The $^{13}$CO brightness temperature and the H$_2$CO line ratio confirm that the disk is too warm for CO freeze out, with the snowline located at $\sim$350 au in the envelope. Both molecules show potential evidence of a temperature increase around the disk-envelope interface. SO seems to originate predominantly in UV-irradiated regions such as the disk surface and the outflow cavity walls rather than at the disk-envelope interface as previously suggested. Finally, the continuum asymmetry along the minor axis is consistent with the inclination derived from the large-scale (100" or 14,000 au) outflow, but opposite to that based on the molecular jet and envelope emission, suggesting a misalignment in the system. Overall, these results highlight the importance of observing multiple molecular species in multiple transitions to characterize the physical and chemical environment of young disks. △ Less

Submitted 27 June, 2023; originally announced June 2023.

Comments: 27 pages, 16 figures, 2 tables, 10 pages appendix with 12 figures. Accepted for publication in ApJ as one of the first-look papers of the eDisk ALMA Large Program

Early Planet Formation in Embedded Disks (eDisk). I. Overview of the Program and First Results

Authors: Nagayoshi Ohashi, John J. Tobin, Jes K. Jørgensen, Shigehisa Takakuwa, Patrick Sheehan, Yuri Aikawa, Zhi-Yun Li, Leslie W. Looney, Jonathan P. Willians, Yusuke Aso, Rajeeb Sharma, Jinshi Sai, Yoshihide Yamato, Jeong-Eun Lee, Kengo Tomida, Hsi-Wei Yen, Frankie J Encalada, Christian Flores, Sacha Gavino, Miyu Kido, Ilseung Han, Zhe-Yu Daniel Lin, Suchitra Narayanan, Nguyen Thi Phuong, Alejandro Santamaría-Miranda , et al. (12 additional authors not shown)

Abstract: We present an overview of the Large Program, ``Early Planet Formation in Embedded Disks (eDisk)'', conducted with the Atacama Large Millimeter/submillimeter Array (ALMA). The ubiquitous detections of substructures, particularly rings and gaps, in protoplanetary disks around T Tauri stars raise the possibility that at least some planet formation may have already started during the embedded stages o… ▽ More We present an overview of the Large Program, ``Early Planet Formation in Embedded Disks (eDisk)'', conducted with the Atacama Large Millimeter/submillimeter Array (ALMA). The ubiquitous detections of substructures, particularly rings and gaps, in protoplanetary disks around T Tauri stars raise the possibility that at least some planet formation may have already started during the embedded stages of star formation. In order to address exactly how and when planet formation is initiated, the program focuses on searching for substructures in disks around 12 Class 0 and 7 Class I protostars in nearby ($< $200 pc) star-forming regions through 1.3 mm continuum observations at a resolution of $\sim7$ au (0.04"). The initial results show that the continuum emission, mostly arising from dust disks around the sample protostars, has relatively few distinctive substructures, such as rings and spirals, in marked contrast to Class II disks. The dramatic difference may suggest that substructures quickly develop in disks when the systems evolve from protostars to Class II sources or alternatively that high optical depth of the continuum emission could obscure internal structures. Kinematic information obtained through CO isotopologue lines and other lines reveals the presence of Keplerian disks around protostars, providing us with crucial physical parameters, in particular, the dynamical mass of the central protostars. We describe the background of the eDisk program, the sample selection and their ALMA observations, the data reduction, and also highlight representative first-look results. △ Less

Submitted 27 June, 2023; originally announced June 2023.

Comments: This is a publication of a series of eDisk ALMA large program first-look papers

Binary Formation in a 100 $μ$m-dark Massive Core

Authors: Shuo Kong, Héctor G. Arce, John J. Tobin, Yichen Zhang, María José Maureira, Kaitlin M. Kratter, Thushara G. S. Pillai

Abstract: We report high-resolution ALMA observations toward a massive protostellar core C1-Sa ($\sim$30 M$_\odot$) in the Dragon Infrared Dark Cloud. At the resolution of 140 AU, the core fragments into two kernels (C1-Sa1 and C1-Sa2) with a projected separation of $\sim$1400 AU along the elongation of C1-Sa, consistent with a Jeans length scale of $\sim$1100 AU. Radiative transfer modeling using RADEX ind… ▽ More We report high-resolution ALMA observations toward a massive protostellar core C1-Sa ($\sim$30 M$_\odot$) in the Dragon Infrared Dark Cloud. At the resolution of 140 AU, the core fragments into two kernels (C1-Sa1 and C1-Sa2) with a projected separation of $\sim$1400 AU along the elongation of C1-Sa, consistent with a Jeans length scale of $\sim$1100 AU. Radiative transfer modeling using RADEX indicates that the protostellar kernel C1-Sa1 has a temperature of $\sim$75 K and a mass of 0.55 M$_\odot$. C1-Sa1 also likely drives two bipolar outflows, one being parallel to the plane-of-the-sky. C1-Sa2 is not detected in line emission and does not show any outflow activity but exhibits ortho-H$_2$D$^+$ and N$_2$D$^+$ emission in its vicinity, thus it is likely still starless. Assuming a 20 K temperature, C1-Sa2 has a mass of 1.6 M$_\odot$. At a higher resolution of 96 AU, C1-Sa1 begins to show an irregular shape at the periphery, but no clear sign of multiple objects or disks. We suspect that C1-Sa1 hosts a tight binary with inclined disks and outflows. Currently, one member of the binary is actively accreting while the accretion in the other is significantly reduced. C1-Sa2 shows hints of fragmentation into two sub-kernels with similar masses, which requires further confirmation with higher sensitivity. △ Less

Submitted 3 May, 2023; originally announced May 2023.

Comments: 7 figures, 2 tables, accepted by ApJ

The Evolution of Protostellar Outflow Cavities, Kinematics, and Angular Distribution of Momentum and Energy in Orion A: Evidence for Dynamical Cores

Authors: Cheng-Han Hsieh, Héctor G. Arce, Zhi-Yun Li, Michael Dunham, Stella Offner, Ian W. Stephens, Amelia Stutz, Tom Megeath, Shuo Kong, Adele Plunkett, John J. Tobin, Yichen Zhang, Diego Mardones, Jaime E. Pineda, Thomas Stanke, John Carpenter

Abstract: We present Atacama Large Millimeter/submillimeter Array observations of the $\sim$10 kAU environment surrounding 21 protostars in the Orion A molecular cloud tracing outflows. Our sample is composed of Class 0 to flat-spectrum protostars, spanning the full $\sim$1 Myr lifetime. We derive the angular distribution of outflow momentum and energy profiles and obtain the first two-dimensional instantan… ▽ More We present Atacama Large Millimeter/submillimeter Array observations of the $\sim$10 kAU environment surrounding 21 protostars in the Orion A molecular cloud tracing outflows. Our sample is composed of Class 0 to flat-spectrum protostars, spanning the full $\sim$1 Myr lifetime. We derive the angular distribution of outflow momentum and energy profiles and obtain the first two-dimensional instantaneous mass, momentum, and energy ejection rate maps using our new approach: the Pixel Flux-tracing Technique (PFT). Our results indicate that by the end of the protostellar phase, outflows will remove $\sim$2 to 4 M$_\odot$ from the surrounding $\sim$1 M$_\odot$ low-mass core. These high values indicate that outflows remove a significant amount of gas from their parent cores and continuous core accretion from larger scales is needed to replenish core material for star formation. This poses serious challenges to the concept of ``cores as well-defined mass reservoirs", and hence to the simplified core-to-star conversion prescriptions. Furthermore, we show that cavity opening angles, and momentum and energy distributions all increase with the protostar evolutionary stage. This is clear evidence that even garden-variety protostellar outflows: (a) effectively inject energy and momentum into their environments on $10$ kAU scales, and (b) significantly disrupt their natal cores, ejecting a large fraction of the mass that would have otherwise fed the nascent star. Our results support the conclusion that protostellar outflows have a direct impact on how stars get their mass, and that the natal sites of individual low-mass star formation are far more dynamic than commonly accepted theoretical paradigms. △ Less

Submitted 10 February, 2023; v1 submitted 6 February, 2023; originally announced February 2023.

Comments: 76 pages, 43 figures. Accepted by ApJ 2023.2.6

The first interferometric measurements of NH$_2$D/NH$_3$ ratio in hot corinos

Authors: Yoshihide Yamato, Kenji Furuya, Yuri Aikawa, Magnus V. Persson, John J. Tobin, Jes K. Jørgensen, Mihkel Kama

Abstract: The nitrogen chemical evolution during star and planet formation is still not fully understood. Ammonia (NH$_3$) is a key specie in the understanding of the molecular evolution in star-forming clouds and nitrogen isotope fractionation. In this paper, we present high spatial resolution observations of multiple emission lines of NH$_3$ toward the protobinary system NGC1333 IRAS4A with Karl G. Jansky… ▽ More The nitrogen chemical evolution during star and planet formation is still not fully understood. Ammonia (NH$_3$) is a key specie in the understanding of the molecular evolution in star-forming clouds and nitrogen isotope fractionation. In this paper, we present high spatial resolution observations of multiple emission lines of NH$_3$ toward the protobinary system NGC1333 IRAS4A with Karl G. Jansky Very Large Array (VLA). We spatially resolved the binary (hereafter 4A1 and 4A2) and detected compact emission of NH$_3$ transitions with high excitation energies ($\gtrsim$100 K) from the vicinity of the protostars, indicating the NH$_3$ ice has sublimated at the inner hot region. The NH$_3$ column density is estimated to be $\sim 10^{17}-10^{18}$ cm$^{-2}$. We also detected two NH$_2$D transitions, allowing us to constrain the deuterium fractionation of ammonia. The NH$_2$D/NH$_3$ ratios are as high as $\sim 0.3-1$ in both 4A1 and 4A2. From the comparisons with the astrochemical models in the literature, the high NH$_2$D/NH$_3$ ratios suggest that the formation of NH$_3$ ices mainly started in the prestellar phase after the formation of bulk water ice finished, and that the primary nitrogen reservoir in the star-forming cloud could be atomic nitrogen (or N atoms) rather than nitrogen-bearing species such as N$_2$ and NH$_3$. The implications on the physical properties of IRAS4A cores are discussed as well. △ Less

Submitted 28 October, 2022; originally announced October 2022.

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

The Disk Population in a Distant Massive Protocluster

Authors: Yu Cheng, Jonathan C. Tan, John J. Tobin, Ruben Fedriani, Morten Andersen, Junfeng Wang

Abstract: The unprecedented angular resolution and sensitivity of ALMA makes it possible to unveil disk populations in distant ($>$2 kpc), embedded young cluster environments. We have conducted an observation towards the central region of the massive protocluster G286.21+0.16 at 1.3 mm. With a spatial resolution of 23 mas and a sensitivity of 15 $\rm μJy~beam^{-1}$, we detect a total of 38 protostellar disk… ▽ More The unprecedented angular resolution and sensitivity of ALMA makes it possible to unveil disk populations in distant ($>$2 kpc), embedded young cluster environments. We have conducted an observation towards the central region of the massive protocluster G286.21+0.16 at 1.3 mm. With a spatial resolution of 23 mas and a sensitivity of 15 $\rm μJy~beam^{-1}$, we detect a total of 38 protostellar disks. These disks have dust masses ranging from about 53 to 1825 $M_\oplus$, assuming a dust temperature of 20 K. This sample is not closely associated with previously identified dense cores, as would be expected for disks around Class 0 protostars. Thus, we expect our sample, being flux limited, to be mainly composed of Class I/flat-spectrum source disks, since these are typically more massive than Class II disks. Furthermore, we find that the distributions of disk masses and radii are statistically indistinguishable with those of the Class I/flat-spectrum objects in the Orion molecular cloud, indicating similar processes are operating in G286.21+0.16 to regulate disk formation and evolution. The cluster center appears to host a massive protostellar system composed of three sources within 1200 au, including a potential binary with 600 au projected separation. Relative to this center, there is no evidence for widespread mass segregation in the disk population. We do find a tentative trend of increasing disk radius versus distance from the cluster center, which may point to the influence of dynamical interactions being stronger in the central regions. △ Less

Submitted 18 October, 2022; originally announced October 2022.

Comments: 23 pages, 7 Figures, accepted by ApJ

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

CORINOS I: JWST/MIRI Spectroscopy and Imaging of a Class 0 protostar IRAS 15398-3359

Authors: Yao-Lun Yang, Joel D. Green, Klaus M. Pontoppidan, Jennifer B. Bergner, L. Ilsedore Cleeves, Neal J. Evans II, Robin T. Garrod, Mihwa Jin, Chul Hwan Kim, Jaeyeong Kim, Jeong-Eun Lee, Nami Sakai, Christopher N. Shingledecker, Brielle Shope, John J. Tobin, Ewine van Dishoeck

Abstract: The origin of complex organic molecules (COMs) in young Class 0 protostars has been one of the major questions in astrochemistry and star formation. While COMs are thought to form on icy dust grains via gas-grain chemistry, observational constraints on their formation pathways have been limited to gas-phase detection. Sensitive mid-infrared spectroscopy with JWST enables unprecedented investigatio… ▽ More The origin of complex organic molecules (COMs) in young Class 0 protostars has been one of the major questions in astrochemistry and star formation. While COMs are thought to form on icy dust grains via gas-grain chemistry, observational constraints on their formation pathways have been limited to gas-phase detection. Sensitive mid-infrared spectroscopy with JWST enables unprecedented investigation of COM formation by measuring their ice absorption features. We present an overview of JWST/MIRI MRS spectroscopy and imaging of a young Class 0 protostar, IRAS 15398-3359, and identify several major solid-state absorption features in the 4.9-28 $μ$m wavelength range. These can be attributed to common ice species, such as H$_2$O, CH$_3$OH, NH$_3$, and CH$_4$, and may have contributions from more complex organic species, such as C$_2$H$_5$OH and CH$_3$CHO. The MRS spectra show many weaker emission lines at 6-8 $μ$m, which are due to warm CO gas and water vapor, possibly from a young embedded disk previously unseen. Finally, we detect emission lines from [Fe II], [Ne II], [S I], and H$_2$, tracing a bipolar jet and outflow cavities. MIRI imaging serendipitously covers the south-western (blue-shifted) outflow lobe of IRAS 15398-3359, showing four shell-like structures similar to the outflows traced by molecular emission at sub-mm wavelengths. This overview analysis highlights the vast potential of JWST/MIRI observations and previews scientific discoveries in the coming years. △ Less

Submitted 14 November, 2022; v1 submitted 22 August, 2022; originally announced August 2022.

Comments: Accepted to ApJL. The methods of data reduction and spectral extraction are updated

A VLA View of the Flared, Asymmetric Disk Around the Class 0 Protostar L1527 IRS

Authors: Patrick D. Sheehan, John J. Tobin, Zhi-Yun Li, Merel L. R. van 't Hoff, Jes K. Jørgensen, Woojin Kwon, Leslie W. Looney, Nagayoshi Ohashi, Shigehisa Takakuwa, Jonathan P. Williams, Yusuke Aso, Sacha Gavino, Itziar de Gregorio-Monsalvo, Ilseung Han, Chang Won Lee, Adele Plunkett, Rajeeb Sharma, Yuri Aikawa, Shih-Ping Lai, Jeong-Eun Lee, Zhe-Yu Daniel Lin, Kazuya Saigo, Kengo Tomida, Hsi-Wei Yen

Abstract: We present high resolution Karl G. Jansky Very Large Array (VLA) observations of the protostar L1527 IRS at 7 mm, 1.3 cm, and 2 cm wavelengths. We detect the edge-on dust disk at all three wavelengths and find that it is asymmetric, with the southern side of the disk brighter than the northern side. We confirm this asymmetry through analytic modeling and also find that the disk is flared at 7 mm.… ▽ More We present high resolution Karl G. Jansky Very Large Array (VLA) observations of the protostar L1527 IRS at 7 mm, 1.3 cm, and 2 cm wavelengths. We detect the edge-on dust disk at all three wavelengths and find that it is asymmetric, with the southern side of the disk brighter than the northern side. We confirm this asymmetry through analytic modeling and also find that the disk is flared at 7 mm. We test the data against models including gap features in the intensity profile, and though we cannot rule such models out, they do not provide a statistically significant improvement in the quality of fit to the data. From these fits, we can however place constraints on allowed properties of any gaps that could be present in the true, underlying intensity profile. The physical nature of the asymmetry is difficult to associate with physical features due to the edge-on nature of the disk, but could be related to spiral arms or asymmetries seen in other imaging of more face-on disks. △ Less

Submitted 27 June, 2022; originally announced June 2022.

Comments: 20 pages, 5 figures, 3 tables; accepted for publication in the Astrophysical Journal

Disks and Outflows in the Intermediate-mass Star Forming Region NGC 2071 IR

Authors: Yu Cheng, John J. Tobin, Yao-Lun Yang, Merel L. R. van't Hoff, Sarah I. Sadavoy, Mayra Osorio, Ana Karla Díaz-Rodríguez, Guillem Anglada, Nicole Karnath, Patrick D. Sheehan, Zhi-Yun Li, Nickalas Reynolds, Nadia M. Murillo, Yichen Zhang, S. Thomas Megeath, Łukasz Tychoniec

Abstract: We present ALMA band 6/7 (1.3 mm/0.87 mm) and VLA Ka band (9 mm) observations toward NGC 2071 IR, an intermediate-mass star forming region. We characterize the continuum and associated molecular line emission towards the most luminous protostars, i.e., IRS1 and IRS3, on ~100 au (0. 2") scales. IRS1 is partly resolved in millimeter and centimeter continuum, which shows a potential disk. IRS3 has a… ▽ More We present ALMA band 6/7 (1.3 mm/0.87 mm) and VLA Ka band (9 mm) observations toward NGC 2071 IR, an intermediate-mass star forming region. We characterize the continuum and associated molecular line emission towards the most luminous protostars, i.e., IRS1 and IRS3, on ~100 au (0. 2") scales. IRS1 is partly resolved in millimeter and centimeter continuum, which shows a potential disk. IRS3 has a well resolved disk appearance in millimeter continuum and is further resolved into a close binary system separated by ~40 au at 9 mm. Both sources exhibit clear velocity gradients across their disk major axes in multiple spectral lines including C18O, H2CO, SO, SO2, and complex organic molecules like CH3OH, 13CH3OH and CH3OCHO. We use an analytic method to fit the Keplerian rotation of the disks, and give constraints on physical parameters with a MCMC routine. The IRS3 binary system is estimated to have a total mass of 1.4-1.5$M_\odot$. IRS1 has a central mass of 3-5$M_\odot$ based on both kinematic modeling and its spectral energy distribution, assuming that it is dominated by a single protostar. For both IRS1 and IRS3, the inferred ejection directions from different tracers, including radio jet, water maser, molecular outflow, and H2 emission, are not always consistent, and for IRS1, these can be misaligned by ~50$^{\circ}$. IRS3 is better explained by a single precessing jet. A similar mechanism may be present in IRS1 as well but an unresolved multiple system in IRS1 is also possible. △ Less

Submitted 30 May, 2022; originally announced May 2022.

Comments: 36 pages, 21 figures, accepted by ApJ

Protostellar and Protoplanetary Disk Masses in the Serpens Region

Authors: Alexa R. Anderson, Jonathan P. Williams, Nienke van der Marel, Charles J. Law, Luca Ricci, John J. Tobin, Simin Tong

Abstract: We present the results from an Atacama Large Millimeter/Submillimeter Array (ALMA) 1.3 mm continuum and $^{12}$CO ($J=2-1$) line survey spread over 10 square degrees in the Serpens star-forming region of 320 young stellar objects, 302 of which are likely members of Serpens (16 Class I, 35 Flat spectrum, 235 Class II, and 16 Class III). From the continuum data, we derive disk dust masses and show t… ▽ More We present the results from an Atacama Large Millimeter/Submillimeter Array (ALMA) 1.3 mm continuum and $^{12}$CO ($J=2-1$) line survey spread over 10 square degrees in the Serpens star-forming region of 320 young stellar objects, 302 of which are likely members of Serpens (16 Class I, 35 Flat spectrum, 235 Class II, and 16 Class III). From the continuum data, we derive disk dust masses and show that they systematically decline from Class I to Flat spectrum to Class II sources. Grouped by stellar evolutionary state, the disk mass distributions are similar to other young ($<3$ Myr) regions, indicating that the large scale environment of a star-forming region does not strongly affect its overall disk dust mass properties. These comparisons between populations reinforce previous conclusions that disks in the Ophiuchus star-forming region have anomalously low masses at all evolutionary stages. Additionally, we find a single deeply embedded protostar that has not been documented elsewhere in the literature and, from the CO line data, 15 protostellar outflows which we catalog here. △ Less

Submitted 19 April, 2022; originally announced April 2022.

Comments: 14 pages, 9 figures, 3 tables, accepted to ApJ

The young embedded disk L1527 IRS: constraints on the water snowline and cosmic ray ionization rate from HCO+ observations

Authors: Merel L. R. van 't Hoff, Margot Leemker, John J. Tobin, Daniel Harsono, Jes K. Jørgensen, Edwin A. Bergin

Abstract: The water snowline in circumstellar disks is a crucial component in planet formation, but direct observational constraints on its location remain sparse due to the difficulty of observing water in both young embedded and mature protoplanetary disks. Chemical imaging provides an alternative route to locate the snowline, and HCO$^+$ isotopologues have been shown to be good tracers in protostellar en… ▽ More The water snowline in circumstellar disks is a crucial component in planet formation, but direct observational constraints on its location remain sparse due to the difficulty of observing water in both young embedded and mature protoplanetary disks. Chemical imaging provides an alternative route to locate the snowline, and HCO$^+$ isotopologues have been shown to be good tracers in protostellar envelopes and Herbig disks. Here we present $\sim$0.5$^{\prime\prime}$ resolution ($\sim$35 au radius) Atacama Large Millimeter/submillimeter Array (ALMA) observations of HCO$^+$ $J=4-3$ and H$^{13}$CO$^+$ $J=3-2$ toward the young (Class 0/I) disk L1527 IRS. Using a source-specific physical model with the midplane snowline at 3.4 au and a small chemical network, we are able to reproduce the HCO$^+$ and H$^{13}$CO$^+$ emission, but for HCO$^+$ only when the cosmic ray ionization rate is lowered to $10^{-18}$ s$^{-1}$. Even though the observations are not sensitive to the expected HCO$^+$ abundance drop across the snowline, the reduction in HCO$^+$ above the snow surface and the global temperature structure allow us to constrain a snowline location between 1.8 and 4.1 au. Deep observations are required to eliminate the envelope contribution to the emission and to derive more stringent constraints on the snowline location. Locating the snowline in young disks directly with observations of H$_2$O isotopologues may therefore still be an alternative option. With a direct snowline measurement, HCO$^+$ will be able to provide constraints on the ionization rate. △ Less

Submitted 18 April, 2022; originally announced April 2022.

Comments: Accepted for publication in ApJ, 15 pages, 6 figures and appendix

The Origin and Evolution of Multiple Star Systems

Authors: Stella S. R. Offner, Maxwell Moe, Kaitlin M. Kratter, Sarah I. Sadavoy, Eric L. N. Jensen, John J. Tobin

Abstract: Observational advances over the last decade have enabled high-resolution, interferometric studies of forming multiple systems, statistical surveys of multiplicity in star-forming regions, and new insights into disk evolution and planetary architectures in these systems. In this review, we compile the results of observational and theoretical studies of stellar multiplicity. We summarize the populat… ▽ More Observational advances over the last decade have enabled high-resolution, interferometric studies of forming multiple systems, statistical surveys of multiplicity in star-forming regions, and new insights into disk evolution and planetary architectures in these systems. In this review, we compile the results of observational and theoretical studies of stellar multiplicity. We summarize the population statistics spanning system evolution from the protostellar phase through the main-sequence phase and evaluate the influence of the local environment. In short, most stars are born in multiple stellar systems, most main sequence stars are members of multiple systems, but most star systems are single. We describe current models for the origin of stellar multiplicity and review the landscape of numerical simulations and assess their consistency with observations. We review the properties of disks and discuss the impact of multiplicity on planet formation and system architectures. Finally, we summarize open questions and discuss the technical requirements for future observational and theoretical progress. △ Less

Submitted 7 December, 2022; v1 submitted 18 March, 2022; originally announced March 2022.

Comments: 36 pages, 9 figures, Protostars and Planets VII, Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura

The VLA/ALMA Nascent Disk and Multiplicity (VANDAM) Survey of Orion Protostars VI. Insights from Radiative Transfer Modeling

Authors: Patrick D. Sheehan, John J. Tobin, Leslie L. Looney, S. Thomas Megeath

Abstract: We present Markov Chain Monte Carlo radiative transfer modeling of a joint ALMA 345 GHz and spectral energy distribution dataset for a sample of 97 protostellar disks from the VLA and ALMA Nascent Disk and Multiplicity Survey of Orion Protostars. From this modeling, we derive disk and envelope properties for each protostar, allowing us to examine the bulk properties of a population of young protos… ▽ More We present Markov Chain Monte Carlo radiative transfer modeling of a joint ALMA 345 GHz and spectral energy distribution dataset for a sample of 97 protostellar disks from the VLA and ALMA Nascent Disk and Multiplicity Survey of Orion Protostars. From this modeling, we derive disk and envelope properties for each protostar, allowing us to examine the bulk properties of a population of young protostars. We find that disks are small, with a median dust radius of $ 29.4^{+ 4.1}_{- 2.7}$ au and a median dust mass of $ 5.8^{+ 4.6}_{- 2.7}$ M$_{\oplus}$. We find no statistically significant difference between most properties of Class 0, I, and Flat Spectrum sources with the exception of envelope dust mass and inclination. The distinction between inclination is an indication that the Class 0/I/Flat Spectrum system may be difficult to tie uniquely to the evolutionary state of protostars. When comparing with Class II disk dust masses in Taurus from similar radiative transfer modeling, we further find that the trend of disk dust mass decreasing from Class 0 to Class II disks is no longer present, though it remains unclear whether such a comparison is fair due to differences in star forming region and modeling techniques. Moreover, the disks we model are broadly gravitationally stable. Finally, we compare disk masses and radii with simulations of disk formation and find that magnetohydrodynamical effects may be important for reproducing the observed properties of disks. △ Less

Submitted 28 February, 2022; originally announced March 2022.

Comments: 31 pages, 15 figures, 2 tables, accepted for publication in ApJ

A 16 au Binary in the Class 0 Protostar L1157 MMS

Authors: John J. Tobin, Erin G. Cox, Leslie W. Looney

Abstract: We present VLA observations toward the Class 0 protostar L1157 MMS at 6.8 mm and 9 mm with a resolution of ~0.04" (14 au). We detect two sources within L1157 MMS and interpret these sources as a binary protostar with a separation of ~16 au. The material directly surrounding the binary system within the inner 50 au radius of the system has an estimated mass of 0.11 M_sun, calculated from the observ… ▽ More We present VLA observations toward the Class 0 protostar L1157 MMS at 6.8 mm and 9 mm with a resolution of ~0.04" (14 au). We detect two sources within L1157 MMS and interpret these sources as a binary protostar with a separation of ~16 au. The material directly surrounding the binary system within the inner 50 au radius of the system has an estimated mass of 0.11 M_sun, calculated from the observed dust emission. We interpret the observed binary system in the context of previous observations of its flattened envelope structure, low rates of envelope rotation from 5000 to 200 au scales, and an ordered, poloidal magnetic field aligned with the outflow. Thus, L1157 MMS is a prototype system for magnetically-regulated collapse and the presence of a compact binary within L1157 MMS demonstrates that multiple star formation can still occur within envelopes that likely have dynamically important magnetic fields. △ Less

Submitted 15 February, 2022; originally announced February 2022.

Comments: 21 pages, 6 figures, 4 tables, accepted to ApJ

Dissecting the different components of the modest accretion bursts of the very young protostar HOPS 373

Authors: Sung-Yong Yoon, Gregory J. Herczeg, Jeong-Eun Lee, Ho-Gyu Lee, Doug Johnstone, Watson Varricatt, John J. Tobin, Carlos Contreras Peña, Steve Mairs, Klaus Hodapp, P. Manoj, Mayra Osorio, S. Thomas Megeath, the JCMT Transient Team

Abstract: Observed changes in protostellar brightness can be complicated to interpret. In our JCMT~Transient monitoring survey, we discovered that a young binary protostar, HOPS 373, is undergoing a modest $30\%$ brightness increase at 850 $μ$m, caused by a factor of 1.8$-$3.3 enhancement in the accretion rate. The initial burst occurred over a few months, with a sharp rise and then shallower decay. A secon… ▽ More Observed changes in protostellar brightness can be complicated to interpret. In our JCMT~Transient monitoring survey, we discovered that a young binary protostar, HOPS 373, is undergoing a modest $30\%$ brightness increase at 850 $μ$m, caused by a factor of 1.8$-$3.3 enhancement in the accretion rate. The initial burst occurred over a few months, with a sharp rise and then shallower decay. A second rise occurred soon after the decay, and the source is still bright one year later. The mid-IR emission, the small-scale CO outflow mapped with ALMA, and the location of variable maser emission indicate that the variability is associated with the SW component. The near-infrared and NEOWISE $W1$ and $W2$ emission is located along the blueshifted CO outflow, spatially offset by $\sim3$ to $4^{\prime\prime}$ from the SW component. The $K$-band emission imaged by UKIRT shows a compact H$_2$ emission source at the edge of the outflow, with a tail tracing the outflow back to the source. The $W1$ emission, likely dominated by scattered light, brightens by 0.7 mag, consistent with expectations based on the sub-mm lightcurve. The signal of continuum variability in $K$-band and $W2$ is masked by stable H$_2$ emission, as seen in our Gemini/GNIRS spectrum, and perhaps by CO emission. These differences in emission sources complicate infrared searches for variability of the youngest protostars. △ Less

Submitted 11 February, 2022; originally announced February 2022.

Comments: 29 pages, 14 figures, accepted for publication in 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

The VLA/ALMA Nascent Disk and Multiplicity (VANDAM) Survey of Orion Protostars V. A Characterization of Protostellar Multiplicity

Authors: John J. Tobin, Stella S. R. Offner, Kaitlin M. Kratter, S. Thomas Megeath, Patrick D. Sheehan, Leslie W. Looney, Ana Karla Diaz-Rodriguez, Mayra Osorio, Guillem Anglada, Sarah I. Sadavoy, Elise Furlan, Dominique Segura-Cox, Nicole Karnath, Merel L. R. van 't Hoff, Ewine F. van Dishoeck, Zhi-Yun Li, Rajeeb Sharma, Amelia M. Stutz, Lukasz Tychoniec

Abstract: We characterize protostellar multiplicity in the Orion molecular clouds using ALMA 0.87~mm and VLA 9~mm continuum surveys toward 328 protostars. These observations are sensitive to projected spatial separations as small as $\sim$20~au, and we consider source separations up to 10$^4$~au as potential companions. The overall multiplicity fraction (MF) and companion fraction (CF) for the Orion protost… ▽ More We characterize protostellar multiplicity in the Orion molecular clouds using ALMA 0.87~mm and VLA 9~mm continuum surveys toward 328 protostars. These observations are sensitive to projected spatial separations as small as $\sim$20~au, and we consider source separations up to 10$^4$~au as potential companions. The overall multiplicity fraction (MF) and companion fraction (CF) for the Orion protostars are 0.30$\pm$0.03 and 0.44$\pm$0.03, respectively, considering separations from 20 to 10$^4$~au. The MFs and CFs are corrected for potential contamination by unassociated young stars using a probabilistic scheme based on the surface density of young stars around each protostar. The companion separation distribution as a whole is double peaked and inconsistent with the separation distribution of solar-type field stars, while the separation distribution of Flat Spectrum protostars is consistent solar-type field stars. The multiplicity statistics and companion separation distributions of the Perseus star-forming region are consistent with those of Orion. Based on the observed peaks in the Class 0 separations at $\sim$100~au and $\sim$10$^3$~au, we argue that multiples with separations $<$500~au are likely produced by both disk fragmentation and turbulent fragmentation with migration, and those at $\ga$10$^3$~au result primarily from turbulent fragmentation. We also find that MFs/CFs may rise from Class 0 to Flat Spectrum protostars between 100 and 10$^3$~au in regions of high YSO density. This finding may be evidence for migration of companions from $>$10$^3$~au to $<$10$^3$~au, and that some companions between 10$^3$ and 10$^4$~au must be (or become) unbound. △ Less

Submitted 8 November, 2021; originally announced November 2021.

Comments: 76 pages, 20 Figures, 10 Tables, accepted by the Astrophysical Journal

Imaging the water snowline around protostars with water and HCO$^+$ isotopologues

Authors: Merel L. R. van 't Hoff, Daniel Harsono, Martijn L. van Gelder, Tien-Hao Hsieh, John J. Tobin, Sigurd S. Jensen, Naomi Hirano, Jes K. Jørgensen, Edwin A. Bergin, Ewine F. van Dishoeck

Abstract: The water snowline location in protostellar envelopes provides crucial information about the thermal structure and the mass accretion process as it can inform about the occurrence of recent ($\lesssim$1,000 yr) accretion bursts. In addition, the ability to image water emission makes these sources excellent laboratories to test indirect snowline tracers such as H$^{13}$CO$^+$. We study the water sn… ▽ More The water snowline location in protostellar envelopes provides crucial information about the thermal structure and the mass accretion process as it can inform about the occurrence of recent ($\lesssim$1,000 yr) accretion bursts. In addition, the ability to image water emission makes these sources excellent laboratories to test indirect snowline tracers such as H$^{13}$CO$^+$. We study the water snowline in five protostellar envelopes in Perseus using a suite of molecular line observations taken with the Atacama Large Millimeter/submillimeter Array (ALMA) at $\sim$0.2$^{\prime\prime}-$0.7$^{\prime\prime}$ (60--210 au) resolution. B1-c provides a textbook example of compact H$_2^{18}$O ($3_{1,3}-2_{2,0}$) and HDO ($3_{1,2}-2_{2,1}$) emission surrounded by a ring of H$^{13}$CO$^+$ ($J=2-1$) and HC$^{18}$O$^+$ ($J=3-2$). Compact HDO surrounded by H$^{13}$CO$^+$ is also detected toward B1-bS. The optically thick main isotopologue HCO$^+$ is not suited to trace the snowline and HC$^{18}$O$^+$ is a better tracer than H$^{13}$CO$^+$ due to a lower contribution from the outer envelope. However, since a detailed analysis is needed to derive a snowline location from H$^{13}$CO$^+$ or HC$^{18}$O$^+$ emission, their true value as snowline tracer will lie in the application in sources where water cannot be readily detected. For protostellar envelopes, the most straightforward way to locate the water snowline is through observations of H$_2^{18}$O or HDO. Including all sub-arcsecond resolution water observations from the literature, we derive an average burst interval of $\sim$10,000 yr, but high-resolution water observations of a larger number of protostars is required to better constrain the burst frequency. △ Less

Submitted 15 October, 2021; originally announced October 2021.

Comments: Accepted for publication in ApJ

Which molecule traces what: chemical diagnostics of protostellar sources

Authors: Łukasz Tychoniec, Ewine F. van Dishoeck, Merel L. R. van 't Hoff, Martijn L. van Gelder, Benoît Tabone, Yuan Chen, Daniel Harsono, Charles L. H. Hull, Michiel R. Hogerheijde, Nadia M. Murillo, John J. Tobin

Abstract: The physical and chemical conditions in Class 0/I protostars are fundamental in unlocking the protostellar accretion process and its impact on planet formation. The aim is to determine which physical components are traced by different molecules at sub-arcsecond scales (100 - 400 au). We use a suite of Atacama Large Millimeter/submillimeter Array (ALMA) datasets in Band 6 (1 mm), Band 5 (1.8 mm) an… ▽ More The physical and chemical conditions in Class 0/I protostars are fundamental in unlocking the protostellar accretion process and its impact on planet formation. The aim is to determine which physical components are traced by different molecules at sub-arcsecond scales (100 - 400 au). We use a suite of Atacama Large Millimeter/submillimeter Array (ALMA) datasets in Band 6 (1 mm), Band 5 (1.8 mm) and Band 3 (3 mm) at spatial resolutions 0.5 - 3" for 16 protostellar sources. The protostellar envelope is well traced by C$^{18}$O, DCO$^+$ and N$_2$D$^+$, with the freeze-out of CO governing the chemistry at envelope scales. Molecular outflows are seen in classical shock tracers like SiO and SO, but ice-mantle products such as CH$_3$OH and HNCO released with the shock are also observed. The molecular jet is prominent not only in SiO and SO but also occasionally in H$_2$CO. The cavity walls show tracers of UV-irradiation such as C$_2$H c-C$_3$H$_2$ and CN. The hot inner envelope, apart from showing emission from complex organic molecules (COMs), also presents compact emission from small molecules like H$_2$S, SO, OCS and H$^{13}$CN, most likely related to ice sublimation and high-temperature chemistry. Sub-arcsecond millimeter-wave observations allow to identify those (simple) molecules that best trace each of the physical components of a protostellar system. COMs are found both in the hot inner envelope (high excitation lines) and in the outflows (lower-excitation lines) with comparable abundances. COMs can coexist with hydrocarbons in the same protostellar sources, but they trace different components. In the near future, mid-IR observations with JWST-MIRI will provide complementary information about the hottest gas and the ice mantle content, at unprecedented sensitivity and at resolutions comparable to ALMA for the same sources. △ Less

Submitted 13 July, 2021; v1 submitted 8 July, 2021; originally announced July 2021.

Comments: 38 pages, 33 figures, accepted for publication to A\&A

Journal ref: A&A 655, A65 (2021)

870 Micron Dust Continuum of the Youngest Protostars in Ophiuchus

Authors: Frankie J. Encalada, Leslie W. Looney, John J. Tobin, Sarah I. Sadavoy, Dominique Segura-Cox, Erin Cox, Zhi-Yun Li, Giles Novak

Abstract: We present a 0.15$^{\prime\prime}$ resolution (21 au) ALMA 870 $μ$m continuum survey of 25 pointings containing 31 young stellar objects in the Ophiuchus molecular clouds. Using the dust continuum as a proxy for dust mass and circumstellar disk radius in our sample, we report a mean mass of 2.8$^{+2.1}_{-1.3}$ and 2.5$^{+9.2}_{-1.1}$ M$_{\oplus}$ and a mean radii of 23.5$^{+1.8}_{-1.2}$ and 16.5… ▽ More We present a 0.15$^{\prime\prime}$ resolution (21 au) ALMA 870 $μ$m continuum survey of 25 pointings containing 31 young stellar objects in the Ophiuchus molecular clouds. Using the dust continuum as a proxy for dust mass and circumstellar disk radius in our sample, we report a mean mass of 2.8$^{+2.1}_{-1.3}$ and 2.5$^{+9.2}_{-1.1}$ M$_{\oplus}$ and a mean radii of 23.5$^{+1.8}_{-1.2}$ and 16.5$^{+2.8}_{-0.9}$ au, for Class I and Flat spectrum protostars, respectively. In addition, we calculate the multiplicity statistics of the dust surrounding young stellar objects in Ophiuchus. The multiplicity fraction (MF) and companion star fraction (CSF) of the combined Class I and Flats based solely on this work is 0.25 $\pm$ 0.09 and 0.33 $\pm$ 0.10, respectively, which are consistent with the values for Perseus and Orion. While we see clear differences in mass and radius between the Ophiuchus and Perseus/Orion protostellar surveys, we do not see any significant differences in the multiplicities of the various regions. We posit there are some differences in the conditions for star formation in Ophiuchus that strongly affects disk size (and consequently disk mass), but does not affect system multiplicity, which could imply important variation in planet formation processes. △ Less

Submitted 21 March, 2024; v1 submitted 2 May, 2021; originally announced May 2021.

Comments: 13 pages, 7 figures, 3 tables. Accepted to ApJ. Corrected typos

An ALMA Survey of Protoplanetary Disks in Lynds 1641

Authors: Sierra L. Grant, Catherine C. Espaillat, John Wendeborn, John J. Tobin, Enrique Macías, Anneliese Rilinger, Álvaro Ribas, S. Thomas Megeath, William J. Fischer, Nuria Calvet, Kyoung Hee Kim

Abstract: We present ALMA observations of 101 protoplanetary disks within the star-forming region Lynds 1641 in the Orion Molecular Cloud A. Our observations include 1.33 mm continuum emission and spectral windows covering the J=2-1 transition of $^{12}$CO, $^{13}$CO, and C$^{18}$O. We detect 89 protoplanetary disks in the dust continuum at the 4$σ$ level ($\sim$88% detection rate) and 31 in $^{12}$CO, 13 i… ▽ More We present ALMA observations of 101 protoplanetary disks within the star-forming region Lynds 1641 in the Orion Molecular Cloud A. Our observations include 1.33 mm continuum emission and spectral windows covering the J=2-1 transition of $^{12}$CO, $^{13}$CO, and C$^{18}$O. We detect 89 protoplanetary disks in the dust continuum at the 4$σ$ level ($\sim$88% detection rate) and 31 in $^{12}$CO, 13 in $^{13}$CO, and 4 in C$^{18}$O. Our sample contains 23 transitional disks, 20 of which are detected in the continuum. We target infrared-bright Class II objects, which biases our sample towards massive disks. We determine dust masses or upper limits for all sources in our sample and compare our sample to protostars in this region. We find a decrease in dust mass with evolutionary state. We also compare this sample to other regions surveyed in the (sub-)millimeter and find that Lynds 1641 has a relatively massive dust disk population compared to regions of similar and older ages, with a median dust mass of 11.1$^{+32.9}_{-4.6}$ $M_\oplus$ and 27% with dust masses equal to or greater than the minimum solar nebula dust mass value of $\sim$30 $M_\oplus$. We analyze the disk mass-accretion rate relationship in this sample and find that the viscous disk lifetimes are similar to the age of the region, however with a large spread. One object, [MGM2012] 512, shows large-scale ($>$5000 AU) structure in both the dust continuum and the three gas lines. We discuss potential origins for this emission, including an accretion streamer with large dust grains. △ Less

Submitted 30 March, 2021; originally announced March 2021.

Comments: 21 pages, 15 figures, 3 tables