Search for: All records

Abstract We use molecular line data from the Atacama Large Millimeter/submillimeter Array, Submillimeter Array, James Clerk Maxwell Telescope, and NANTEN2 to study the multiscale (∼15–0.005 pc) velocity statistics in the massive star formation region NGC 6334. We find that the nonthermal motions revealed by the velocity dispersion function (VDF) stay supersonic over scales of several orders of magnitude. The multiscale nonthermal motions revealed by different instruments do not follow the same continuous power law, which is because the massive star formation activities near central young stellar objects have increased the nonthermal motions in small-scale and high-density regions. The magnitudes of VDFs vary in different gas materials at the same scale, where the infrared dark clump N6334S in an early evolutionary stage shows a lower level of nonthermal motions than other more evolved clumps due to its more quiescent star formation activity. We find possible signs of small-scale-driven (e.g., by gravitational accretion or outflows) supersonic turbulence in clump N6334IV with a three-point VDF analysis. Our results clearly show that the scaling relation of velocity fields in NGC 6334 deviates from a continuous and universal turbulence cascade due to massive star formation activities. 

ABSTRACT We compare the structure of synthetic dust polarization with synthetic molecular line emission from radiative transfer calculations using a three-dimensional, turbulent collapsing-cloud magnetohydrodynamics simulation. The histogram of relative orientation (HRO) technique and the projected Rayleigh statistic (PRS) are considered. In our trans-Alfvénic (more strongly magnetized) simulation, there is a transition to perpendicular alignment at densities above ∼4 × 103 cm−3. This transition is recovered in most of our synthetic observations of optically thin molecular tracers; however, for 12CO it does not occur and the PRS remains in parallel alignment across the whole observer space. We calculate the physical depth of the optical depth τ = 1 surface and find that for 12CO it is largely located in front of the cloud mid-plane, suggesting that 12CO is too optically thick and instead mainly probes low-volume density gas. In our super-Alfvénic simulation, the magnetic field becomes significantly more tangled, and all observed tracers tend towards no preference for perpendicular or parallel alignment. An observable difference in alignment between optically thin and optically thick tracers may indicate the presence of a dynamically important magnetic field, though there is some degeneracy with viewing angle. We convolve our data with a Gaussian beam and compare it with HRO results of the Vela C molecular cloud. We find good agreement between these results and our sub-Alfvénic simulations when viewed with the magnetic field in the plane of the sky (especially when sensitivity limitations are considered), though the observations are also consistent with an intermediately inclined magnetic field. 

Abstract Magnetic fields play essential roles in protoplanetary disks. Magnetic fields in the disk atmosphere are of particular interest, as they are connected to the wind-launching mechanism. In this work, we study the polarization of the light scattered off of magnetically aligned grains in the disk atmosphere, focusing on the deviation of the polarization orientation from the canonical azimuthal direction, which may be detectable in near-IR polarimetry with instruments such as VLT/SPHERE. We show with a simple disk model that the polarization can even be oriented along the radial (rather than azimuthal) direction, especially in highly inclined disks with toroidally dominated magnetic fields. This polarization reversal is caused by the anisotropy in the polarizability of aligned grains and is thus a telltale sign of such grains. We show that the near-IR light is scattered mostly byμm-sized grains or smaller at theτ= 1 surface and such grains can be magnetically aligned if they contain superparamagnetic inclusions. For comparison with observations, we generate synthetic maps of the ratios ofU_ϕ/IandQ_ϕ/I, which can be used to infer the existence of (magnetically) aligned grains through a negativeQ_ϕ(polarization reversal) and/or a significant level ofU_ϕ/I. We show that two features observed in the existing data, an asymmetric distribution ofU_ϕwith respect to the disk minor axis and a spatial distribution ofU_ϕthat is predominantly positive or negative, are incompatible with scattering by spherical grains in an axisymmetric disk. They provide indirect evidence for scattering by aligned nonspherical grains. 

Abstract We present ALMA dust polarization and molecular line observations toward four clumps (I(N), I, IV, and V) in the massive star-forming region NGC 6334. In conjunction with large-scale dust polarization and molecular line data from JCMT, Planck, and NANTEN2, we make a synergistic analysis of relative orientations between magnetic fields (θ_B), column density gradients (θ_NG), local gravity (θ_LG), and velocity gradients (θ_VG) to investigate the multi-scale (from ∼30 to 0.003 pc) physical properties in NGC 6334. We find that the relative orientation betweenθ_Bandθ_NGchanges from statistically more perpendicular to parallel as column density ( $N_{{\mathrm{H}}_2}$) increases, which is a signature of trans-to-sub-Alfvénic turbulence at complex/cloud scales as revealed by previous numerical studies. Becauseθ_NGandθ_LGare preferentially aligned within the NGC 6334 cloud, we suggest that the more parallel alignment betweenθ_Bandθ_NGat higher $N_{{\mathrm{H}}_2}$is because the magnetic field line is dragged by gravity. At even higher $N_{{\mathrm{H}}_2}$, the angle betweenθ_Bandθ_NGorθ_LGtransits back to having no preferred orientation, or statistically slightly more perpendicular, suggesting that the magnetic field structure is impacted by star formation activities. A statistically more perpendicular alignment is found betweenθ_Bandθ_VGthroughout our studied $N_{{\mathrm{H}}_2}$range, which indicates a trans-to-sub-Alfvénic state at small scales as well, and this signifies that magnetic field has an important role in the star formation process in NGC 6334. The normalized mass-to-flux ratio derived from the polarization-intensity gradient (KTH) method increases with $N_{{\mathrm{H}}_2}$, but the KTH method may fail at high $N_{{\mathrm{H}}_2}$due to the impact of star formation feedback. 

ABSTRACT Recent high angular resolution ALMA observations have revealed numerous gaps in protoplanetary discs. A popular interpretation has been that planets open them. Most previous investigations of planet gap-opening have concentrated on viscous discs. Here, we carry out 2D (axisymmetric) global simulations of gap opening by a planet in a wind-launching non-ideal MHD disc with consistent thermochemistry. We find a strong concentration of poloidal magnetic flux in the planet-opened gap, where the gas dynamics are magnetically dominated. The magnetic field also drives a fast (nearly sonic) meridional gas circulation in the denser disc regions near the inner and outer edges of the gap, which may be observable through high-resolution molecular line observations. The gap is more ionized than its denser surrounding regions, with a better magnetic field–matter coupling. In particular, it has a much higher abundance of molecular ion HCO+, consistent with ALMA observations of the well-studied AS 209 protoplanetary disc that has prominent gaps and fast meridional motions reaching the local sound speed. Finally, we provide fitting formulae for the ambipolar and Ohmic diffusivities as a function of the disc local density, which can be used for future 3D simulations of planet gap-opening in non-ideal MHD discs where thermochemistry is too computationally expensive to evolve self-consistently with the magneto-hydrodynamics. 

ABSTRACT Recent observations indicate that mm/cm-sized grains may exist in the embedded protostellar discs. How such large grains grow from the micron size (or less) in the earliest phase of star formation remains relatively unexplored. In this study, we take a first step to model the grain growth in the protostellar environment, using 2D (axisymmetric) radiation hydrodynamic and grain growth simulations. We show that the grain growth calculations can be greatly simplified by the ‘terminal velocity approximation’, where the dust drift velocity relative to the gas is proportional to its stopping time, which is proportional to the grain size. We find that the grain–grain collision from size-dependent terminal velocity alone is too slow to convert a significant fraction of the initially micron-sized grains into mm/cm sizes during the deeply embedded Class 0 phase. Substantial grain growth is achieved when the grain–grain collision speed is enhanced by a factor of 4. The dust growth above and below the disc midplane enables the grains to settle faster towards the midplane, which increases the local dust-to-gas ratio, which, in turn, speeds up further growth there. How this needed enhancement can be achieved is unclear, although turbulence is a strong possibility that deserves further exploration. 

Abstract Chemical models and experiments indicate that interstellar dust grains and their ice mantles play an important role in the production of complex organic molecules (COMs). To date, the most complex solid-phase molecule detected with certainty in the interstellar medium is methanol, but the James Webb Space Telescope (JWST) may be able to identify still larger organic species. In this study, we use a coupled chemodynamical model to predict new candidate species for JWST detection toward the young star-forming core Cha-MMS1, combining the gas–grain chemical kinetic code MAGICKAL with a 1D radiative hydrodynamics simulation using Athena++ . With this model, the relative abundances of the main ice constituents with respect to water toward the core center match well with typical observational values, providing a firm basis to explore the ice chemistry. Six oxygen-bearing COMs (ethanol, dimethyl ether, acetaldehyde, methyl formate, methoxy methanol, and acetic acid), as well as formic acid, show abundances as high as, or exceeding, 0.01% with respect to water ice. Based on the modeled ice composition, the infrared spectrum is synthesized to diagnose the detectability of the new ice species. The contribution of COMs to IR absorption bands is minor compared to the main ice constituents, and the identification of COM ice toward the core center of Cha-MMS1 with the JWST NIRCAM/Wide Field Slitless Spectroscopy (2.4–5.0 μ m) may be unlikely. However, MIRI observations (5–28 μ m) toward COM-rich environments where solid-phase COM abundances exceed 1% with respect to the column density of water ice might reveal the distinctive ice features of COMs. 

Abstract We present Atacama Large Millimeter Array band 6/7 (1.3 mm/0.87 mm) and Very Large Array Ka-band (9 mm) observations toward NGC 2071 IR, an intermediate-mass star-forming region. We characterize the continuum and associated molecular line emission toward the most luminous protostars, i.e., IRS1 and IRS3, on ∼100 au (0.″2) scales. IRS1 is partly resolved in the millimeter and centimeter continuum, which shows a potential disk. IRS3 has a well-resolved disk appearance in the 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 C 18 O, H 2 CO, SO, SO 2 , and complex organic molecules like CH 3 OH, 13 CH 3 OH, and CH 3 OCHO. We use an analytic method to fit the Keplerian rotation of the disks and give constraints on physical parameters with a Markov Chain Monte Carlo routine. The IRS3 binary system is estimated to have a total mass of 1.4–1.5 M ⊙ . IRS1 has a central mass of 3–5 M ⊙ 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 H 2 emission, are not always consistent, and for IRS1 these can be misaligned by ∼50°. 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. 

ABSTRACT Despite the rich observational results on interstellar magnetic fields in star-forming regions, it is still unclear how dynamically significant the magnetic fields are at varying physical scales, because direct measurement of the field strength is observationally difficult. The Davis–Chandrasekhar–Fermi (DCF) method has been the most commonly used method to estimate the magnetic field strength from polarization data. It is based on the assumption that gas turbulent motion is the driving source of field distortion via linear Alfvén waves. In this work, using MHD simulations of star-forming clouds, we test the validity of the assumption underlying the DCF method by examining its accuracy in the real 3D space. Our results suggest that the DCF relation between turbulent kinetic energy and magnetic energy fluctuation should be treated as a statistical result instead of a local property. We then develop and investigate several modifications to the original DCF method using synthetic observations, and propose new recipes to improve the accuracy of DCF-derived magnetic field strength. We further note that the biggest uncertainty in the DCF analysis may come from the linewidth measurement instead of the polarization observation, especially since the line-of-sight gas velocity can be used to estimate the gas volume density, another critical parameter in the DCF method. 

Abstract We present H -band (1.65 μ m) and SOFIA HAWC+ 154 μ m polarization observations of the low-mass core L483. Our H -band observations reveal a magnetic field that is overwhelmingly in the E–W direction, which is approximately parallel to the bipolar outflow that is observed in scattered IR light and in single-dish 12 CO observations. From our 154 μ m data, we infer a ∼45° twist in the magnetic field within the inner 5″ (1000 au) of L483. We compare these new observations with published single-dish 350 μ m polarimetry and find that the 10,000 au scale H -band data match the smaller-scale 350 μ m data, indicating that the collapse of L483 is magnetically regulated on these larger scales. We also present high-resolution 1.3 mm Atacama Large Millimeter/submillimeter Array data of L483 that reveals it is a close binary star with a separation of 34 au. The plane of the binary of L483 is observed to be approximately parallel to the twisted field in the inner 1000 au. Comparing this result to the ∼1000 au protostellar envelope, we find that the envelope is roughly perpendicular to the 1000 au HAWC+ field. Using the data presented, we speculate that L483 initially formed as a wide binary and the companion star migrated to its current position, causing an extreme shift in angular momentum thereby producing the twisted magnetic field morphology observed. More observations are needed to further test this scenario.