Episodic accretion has been used to explain the wide range of protostellar luminosities, but its origin and influence on the star-forming process are not yet fully understood. We present an ALMA ...survey of N2H+ (1−0) and HCO+ (3−2) toward 39 Class 0 and Class I sources in the Perseus molecular cloud. N2H+ and HCO+ are destroyed via gas-phase reactions with CO and H2O, respectively, thus tracing the CO and H2O snowline locations. A snowline location at a much larger radius than that expected from the current luminosity suggests that an accretion burst has occurred in the past that has shifted the snowline outward. We identified 18/18 Class 0 and 9/10 Class I post-burst sources from N2H+ and 7/17 Class 0 and 1/8 Class I post-burst sources from HCO+. The accretion luminosities during the past bursts are found to be ∼10-100 L . This result can be interpreted as either evolution of burst frequency or disk evolution. In the former case, assuming that refreeze-out timescales are 1000 yr for H2O and 10,000 yr for CO, we found that the intervals between bursts increase from 2400 yr in the Class 0 stage to 8000 yr in the Class I stage. This decrease in the burst frequency may reflect that fragmentation is more likely to occur at an earlier evolutionary stage when the young stellar object is more prone to instability.
Abstract
We investigate ionization and heating of gas in the dense, shielded clumps/cores of molecular clouds bathed by an influx of energetic, charged cosmic rays (CRs). These molecular clouds have ...complex structures, with substantial variation in their physical properties over a wide range of length scales. The propagation and distribution of CRs is thus regulated accordingly, in particular, by the magnetic fields threaded through the clouds and into the dense regions within. We have found that a specific heating rate reaching 10
−26
erg cm
−3
s
−1
can be sustained in the dense clumps/cores for Galactic environments, and this rate increases with CR energy density. The propagation of CRs and heating rates in some star-forming filaments identified in IC 5146 are calculated, with the CR diffusion coefficients in these structures determined from magnetic field fluctuations inferred from optical and near-infrared polarizations of starlight, which is presumably a magnetic field tracer. Our calculations indicate that CR heating can vary by nearly three orders of magnitude between different filaments within a cloud due to different levels of CR penetration. The CR ionization rate among these filaments is similar. The equilibrium temperature that could be maintained by CR heating alone is of order 1 K in a Galactic environment, but this value would be higher in strongly star-forming environments, thus causing an increase in the Jeans mass of their molecular clouds.
Polarized emission is detected in two young nearly edge-on protostellar disks in 343 GHz continuum at ∼50 au (∼0 12) resolution with Atacama Large Millimeter/submillimeter Array. One disk is in HH ...212 (Class 0) and the other in the HH 111 (early Class I) protostellar system. The polarization fraction is ∼1%. The disk in HH 212 has a radius of ∼60 au. The emission is mainly detected from the nearside of the disk. The polarization orientations are almost perpendicular to the disk major axis, consistent with either self-scattering or emission by grains aligned with a poloidal field around the outer edge of the disk because of the optical depth effect and temperature gradient; the presence of a poloidal field would facilitate the launching of a disk wind, for which there is already tentative evidence in the same source. The disk of HH 111 VLA 1 has a larger radius of ∼220 au and is thus more resolved. The polarization orientations are almost perpendicular to the disk major axis in the nearside, but more along the major axis in the farside, forming roughly half of an elliptical pattern there. It appears that toroidal and poloidal magnetic field may explain the polarization on the near and far sides of the disk, respectively. However, it is also possible that the polarization is due to self-scattering. In addition, alignment of dust grains by radiation flux may play a role in the farside. Our observations reveal a diversity of disk polarization patterns that should be taken into account in future modeling efforts.
We present high angular resolution dust polarization and molecular line observations carried out with the Atacama Large Millimeter/submillimeter Array (ALMA) toward the Class 0 protostar Serpens ...SMM1. By complementing these observations with new polarization observations from the Submillimeter Array (SMA) and archival data from the Combined Array for Research in Millimeter-wave Astronomy (CARMA) and the James Clerk Maxwell Telescopes (JCMT), we can compare the magnetic field orientations at different spatial scales. We find major changes in the magnetic field orientation between large (∼0.1 pc) scales-where the magnetic field is oriented E-W, perpendicular to the major axis of the dusty filament where SMM1 is embedded-and the intermediate and small scales probed by CARMA (∼1000 au resolution), the SMA (∼350 au resolution), and ALMA (∼140 au resolution). The ALMA maps reveal that the redshifted lobe of the bipolar outflow is shaping the magnetic field in SMM1 on the southeast side of the source; however, on the northwestern side and elsewhere in the source, low-velocity shocks may be causing the observed chaotic magnetic field pattern. High-spatial-resolution continuum and spectral-line observations also reveal a tight (∼130 au) protobinary system in SMM1-b, the eastern component of which is launching an extremely high-velocity, one-sided jet visible in both ) and ); however, that jet does not appear to be shaping the magnetic field. These observations show that with the sensitivity and resolution of ALMA, we can now begin to understand the role that feedback (e.g., from protostellar outflows) plays in shaping the magnetic field in very young, star-forming sources like SMM1.
Aims. Although from a theoretical point of view magnetic fields are believed to play a significant role during the early stages of star formation, especially during the main accretion phase, the ...magnetic field strength and topology is poorly constrained in the youngest accreting Class 0 protostars that lead to the formation of solar-type stars. Methods. We carried out observations of the polarized dust continuum emission with the SMA interferometer at 0.87 mm to probe the structure of the magnetic field in a sample of 12 low-mass Class 0 envelopes in nearby clouds, including both single protostars and multiple systems. Our SMA observations probed the envelope emission at scales ~600 − 5000 au with a spatial resolution ranging from 600 to 1500 au depending on the source distance. Results. We report the detection of linearly polarized dust continuum emission in all of our targets with average polarization fractions ranging from 2% to 10% in these protostellar envelopes. The polarization fraction decreases with the continuum flux density, which translates into a decrease with the H2 column density within an individual envelope. Our analysis show that the envelope-scale magnetic field is preferentially observed either aligned or perpendicular to the outflow direction. Interestingly, our results suggest for the first time a relation between the orientation of the magnetic field and the rotational energy of envelopes, with a larger occurrence of misalignment in sources in which strong rotational motions are detected at hundreds to thousands of au scales. We also show that the best agreement between the magnetic field and outflow orientation is found in sources showing no small-scale multiplicity and no large disks at ~100 au scales.
We report Atacama Large Millimeter/submillimeter Array (ALMA) observations of polarized dust emission from the protostellar source Ser-emb 8 at a linear resolution of 140 au. Assuming models of ...dust-grain alignment hold, the observed polarization pattern gives a projected view of the magnetic field structure in this source. Contrary to expectations based on models of strongly magnetized star formation, the magnetic field in Ser-emb 8 does not exhibit an hourglass morphology. Combining the new ALMA data with previous observational studies, we can connect magnetic field structure from protostellar core (∼80,000 au) to disk (∼100 au) scales. We compare our observations with four magnetohydrodynamic gravo-turbulence simulations made with the AREPO code that have initial conditions ranging from super-Alfvénic (weakly magnetized) to sub-Alfvénic (strongly magnetized). These simulations achieve the spatial dynamic range necessary to resolve the collapse of protostars from the parsec scale of star-forming clouds down to the ∼100 au scale probed by ALMA. Only in the very strongly magnetized simulation do we see both the preservation of the field direction from cloud to disk scales and an hourglass-shaped field at <1000 au scales. We conduct an analysis of the relative orientation of the magnetic field and the density structure in both the Ser-emb 8 ALMA observations and the synthetic observations of the four AREPO simulations. We conclude that the Ser-emb 8 data are most similar to the weakly magnetized simulations, which exhibit random alignment, in contrast to the strongly magnetized simulation, where the magnetic field plays a role in shaping the density structure in the source. In the weak-field case, it is turbulence-not the magnetic field-that shapes the material that forms the protostar, highlighting the dominant role that turbulence can play across many orders of magnitude in spatial scale.
We determine the magnetic field strength in the OMC 1 region of the Orion A filament via a new implementation of the Chandrasekhar-Fermi method using observations performed as part of the James Clerk ...Maxwell Telescope (JCMT) B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. We combine BISTRO data with archival SCUBA-2 and HARP observations to find a plane-of-sky magnetic field strength in OMC 1 of mG, where mG represents a predominantly systematic uncertainty. We develop a new method for measuring angular dispersion, analogous to unsharp masking. We find a magnetic energy density of J m−3 in OMC 1, comparable both to the gravitational potential energy density of OMC 1 (∼10−7 J m−3) and to the energy density in the Orion BN/KL outflow (∼10−7 J m−3). We find that neither the Alfvén velocity in OMC 1 nor the velocity of the super-Alfvénic outflow ejecta is sufficiently large for the BN/KL outflow to have caused large-scale distortion of the local magnetic field in the ∼500 yr lifetime of the outflow. Hence, we propose that the hourglass field morphology in OMC 1 is caused by the distortion of a primordial cylindrically symmetric magnetic field by the gravitational fragmentation of the filament and/or the gravitational interaction of the BN/KL and S clumps. We find that OMC 1 is currently in or near magnetically supported equilibrium, and that the current large-scale morphology of the BN/KL outflow is regulated by the geometry of the magnetic field in OMC 1, and not vice versa.
The formation of filaments in molecular clouds is an important process in star formation. Hub-filament systems (HFSs) are a transition stage connecting parsec-scale filaments and protoclusters. ...Understanding the origin of HFSs is crucial to reveal how star formation proceeds from clouds to cores. Here we report James Clerk Maxwell telescope POL-2 850 m polarization and IRAM 30 m C18O (2-1) line observations toward the massive HFS G33.92+0.11. The 850 m continuum map reveals four major filaments converging to the center of G33.92+0.11 with numerous short filaments connecting to the major filaments at local intensity peaks. We estimate the local orientations of filaments, magnetic field, gravity, and velocity gradients from observations, and we examine their correlations based on their local properties. In the high-density areas, our analysis shows that the filaments tend to align with the magnetic field and local gravity. In the low-density areas, we find that the local velocity gradients tend to be perpendicular to both the magnetic field and local gravity, although the filaments still tend to align with local gravity. A global virial analysis suggests that the gravitational energy overall dominates the magnetic and kinematic energy. Combining local and global aspects, we conclude that the formation of G33.92+0.11 is predominantly driven by gravity, dragging and aligning the major filaments and magnetic field on the way to the inner dense center. Traced by local velocity gradients in the outer diffuse areas, ambient gas might be accreted onto the major filaments directly or via the short filaments.
We present optical and near-infrared stellar polarization observations toward the dark filamentary clouds associated with IC5146. The data allow us to investigate the dust properties (this paper) and ...the magnetic field structure (Paper II). A total of 2022 background stars were detected in the Rc, , H, and/or K bands to mag. The ratio of the polarization percentage at different wavelengths provides an estimate of , the wavelength of the peak polarization, which is an indicator of the small-size cutoff of the grain size distribution. The grain size distribution seems to significantly change at mag, where both the average and dispersion of decrease. In addition, we found m for mag, which is larger than the ∼0.55 m in the general interstellar medium (ISM), suggesting that grain growth has already started in low-AV regions. Our data also reveal that polarization efficiency ( ) decreases with AV as a power law in the Rc, , and K bands with indices of −0.71 0.10, −1.23 0.10, and −0.53 0.09. However, H-band data show a power index change; the PE varies with AV steeply (index of −0.95 0.30) when mag, but softly (index of −0.25 0.06) for greater AV values. The soft decay of PE in high-AV regions is consistent with the radiative alignment torque model, suggesting that our data trace the magnetic field to mag. Furthermore, the breakpoint found in the H band is similar to that for AV, where we found the dispersion significantly decreased. Therefore, the flat PE-AV in high-AV regions implies that the power-index changes result from additional grain growth.
The IC 5146 cloud is a nearby star-forming region in Cygnus, consisting of molecular gas filaments in a variety of evolutionary stages. We used optical and near-infrared polarization data toward the ...IC 5146 cloud, reported in the first paper of this series, to reveal the magnetic fields in this cloud. Using the newly released Gaia data, we found that the IC 5146 cloud may contain two separate clouds: a first cloud, including the densest main filament at a distance of ∼600 pc, and a second cloud, associated with the Cocoon Nebula at a distance of ∼800 pc. The spatially averaged H-band polarization map revealed a well-ordered magnetic field morphology, with the polarization segments perpendicular to the main filament but parallel to the nearby subfilaments, consistent with models assuming that the magnetic field is regulating cloud evolution. We estimated the magnetic field strength using the Davis-Chandrasekhar-Fermi method and found that the magnetic field strength scales with volume density with a power-law index of ∼0.5 in the density range from to 3000 cm−3, which indicates an anisotropic cloud contraction with a preferred direction along the magnetic field. In addition, the mass-to-flux ratio of the cloud gradually changes from subcritical to supercritical from the cloud envelope to the deep regions. These features are consistent with strong magnetic field star formation models and suggest that the magnetic field is important in regulating the evolution of the IC 5146 cloud.